Original Research ARTICLE
Photosystem II Repair and Plant Immunity: Lessons Learned from Arabidopsis Mutant Lacking the THYLAKOID LUMEN PROTEIN 18.3
- 1Molecular Plant Biology, Department of Biochemistry, University of Turku, Turku, Finland
- 2Plant Biology, Department of Biosciences, University of Helsinki, Helsinki, Finland
- 3Molecular Biomimetics, Department of Chemistry—Ångström Laboratory, Uppsala University, Uppsala, Sweden
Chloroplasts play an important role in the cellular sensing of abiotic and biotic stress. Signals originating from photosynthetic light reactions, in the form of redox and pH changes, accumulation of reactive oxygen and electrophile species or stromal metabolites are of key importance in chloroplast retrograde signaling. These signals initiate plant acclimation responses to both abiotic and biotic stresses. To reveal the molecular responses activated by rapid fluctuations in growth light intensity, gene expression analysis was performed with Arabidopsis thaliana wild type and the tlp18.3 mutant plants, the latter showing a stunted growth phenotype under fluctuating light conditions (Biochem. J, 406, 415–425). Expression pattern of genes encoding components of the photosynthetic electron transfer chain did not differ between fluctuating and constant light conditions, neither in wild type nor in tlp18.3 plants, and the composition of the thylakoid membrane protein complexes likewise remained unchanged. Nevertheless, the fluctuating light conditions repressed in wild-type plants a broad spectrum of genes involved in immune responses, which likely resulted from shade-avoidance responses and their intermixing with hormonal signaling. On the contrary, in the tlp18.3 mutant plants there was an imperfect repression of defense-related transcripts upon growth under fluctuating light, possibly by signals originating from minor malfunction of the photosystem II (PSII) repair cycle, which directly or indirectly modulated the transcript abundances of genes related to light perception via phytochromes. Consequently, a strong allocation of resources to defense reactions in the tlp18.3 mutant plants presumably results in the stunted growth phenotype under fluctuating light.
Photosystem II (PSII), embedded in the thylakoid membranes, catalyzes light-dependent water splitting with concomitant oxygen evolution and electron transfer to the plastoquinone pool. PSII consists of the chloroplast-encoded core subunits D1, D2, CP43, and CP47, as well as numerous other subunits, encoded by both the chloroplast and nuclear genomes. Of these proteins, the nuclear-encoded proteins PsbO, PsbP, and PsbQ together with the manganese-calcium cluster form the so called oxygen-evolving complex (OEC), located at the lumenal surface of the PSII complex. In higher plants, the functional PSII complex is formed as a PSII dimer, to which nuclear-encoded light-harvesting complex (LHC) II proteins, Lhcb1-6, are tightly connected forming PSII-LHCII supercomplexes.
Photosynthetic water splitting and evolution of one oxygen molecule require four sequential excitations and subsequent charge separations in the reaction center chlorophyll (Chl) P680, thus producing extremely oxidizing, and potentially hazardous reactive oxygen species (ROS), which enhance oxidative damage to PSII as well as to other thylakoid proteins (Krieger-Liszkay et al., 2008; Pospísil, 2009). Despite the existence of detoxification systems for scavenging of ROS, damage to PSII is unavoidable (Aro et al., 1993; Tyystjärvi and Aro, 1996; Takahashi and Badger, 2011). In particular, the PSII core protein D1 is prone to light-induced damage, and thus an efficient repair cycle has evolved for PSII, which includes proteolytic degradation of damaged D1 protein and its replacement with a newly-synthetized D1 copy (reviewed in Baena-Gonzalez and Aro, 2002; Edelman and Mattoo, 2008; Nixon et al., 2010). These processes involve reversible monomerization of the PSII-LHCII supercomplexes (Danielsson et al., 2006), as well as dynamic changes in grana diameter and in lumen volume (Kirchhoff et al., 2011; Herbstova et al., 2012). A vast number of auxiliary proteins, such as kinases, phosphatases, proteases, transporters, and chaperones have been shown to assist the PSII repair cycle (reviewed in Mulo et al., 2008; Chi et al., 2012; Nickelsen and Rengstl, 2013; Järvi et al., 2015). One of these, the THYLAKOID LUMEN PROTEIN OF 18.3 kDa (TLP18.3) has been shown to be required for efficient degradation of the damaged D1 protein and dimerization of the PSII complex (Sirpiö et al., 2007). Notably, high light treatment challenging the PSII repair cycle triggered only a moderate damage of PSII in tlp18.3 plants (Sirpiö et al., 2007), which suggest that TLP18.3 is not a crucial component of the repair cycle but instead plays a role in fine tuning the repair cycle. Based on structural data, TLP18.3 has been suggested to be an acidic phosphatase, but only low phosphatase activity was measured for TLP18.3 (Wu et al., 2011). Recently, the regulatory role of the PSII repair cycle has been extended to include the maintenance of photosystem I (PSI) and indeed, insufficient regulation of the PSII repair cycle seems to exert an effect also on the function of PSI (Tikkanen et al., 2014). Moreover, PSII is crucial for plant immunity through production of ROS, which are not only damaging the components of the photosynthetic electron transfer chain, but also act as important retrograde signaling molecules (Rodríguez-Herva et al., 2012; de Torres Zabala et al., 2015). In line with this, a functional connection between PSII repair and regulation of cell death in tobacco leaves infected by tobacco mosaic virus has been established (Seo et al., 2000).
While the exact role of photosynthetic components in sensing and signaling the pathogen infection is only emerging, a wealth of information has accumulated during the past few years on the consequences of fluctuating light on the activity of the photosynthetic machinery (Grieco et al., 2012; Suorsa et al., 2012; Allahverdiyeva et al., 2013; Kono and Terashima, 2014). Nevertheless, we still lack knowledge on how the rapid fluctuations in growth light intensity affect the acclimation processes at the level of nuclear gene expression, and even less is known about potential cross-talk between light acclimation, the PSII repair cycle and disease resistance under fluctuating light. Here, we investigated how the constantly fluctuating growth light intensity modulates the transcript profile of wild-type Arabidopsis thaliana (hereafter Arabidopsis) plants, and how such an acclimation response is further affected by the deficiency of the thylakoid lumen protein TLP18.3. Five-week old plants grown either under constant or fluctuating light conditions for their entire life span were used as material to study the late stage of the acclimation process.
Materials and Methods
Plant Material and Growth Conditions
Arabidopsis, ecotype Columbia 0, wild-type and tlp18.3 (GABI-Kat 459D12) plants (Sirpiö et al., 2007) were used in all experiments. Plants were grown in 8 h light regime at 23°C either under a photon flux density of 120 μmol photons m−2 s−1 or under fluctuating light intensities, in which plants were exposed to 50 μmol photons m−2 s−1 for 5 min and subsequently to high-light of 500 μmol photons m−2 s−1 for 1 min (Tikkanen et al., 2010), the cycles being repeated during the entire photoperiod. Osram HQI-BT 400 W/D Metal Halide lamps with spectral power distribution from 350 to 800 nm were used as a light source. Five-week-old plants were used for all experiments.
Gene Expression Analyses
Microarray analyses of wild-type and tlp18.3 plants were performed essentially as in Konert et al. (2015). In short, leaf material was harvested 4 h after the onset of the light period in order to be sure that the plants were in a photosynthetically active state and that the PSII repair cycle was properly ongoing and immediately frozen in liquid nitrogen. RNA was isolated using an Agilent Plant RNA isolation mini kit according to manufacturer's instructions. Cy-3 labeled RNA samples were hybridized to Agilent Arabidopsis Gene Expression Microarrays, 4 × 44 K (Design ID 021169) and scanned with Agilent Technologies Scanner G2565CA with a profile AgilentHD _GX_1Color. Numeric data were produced with Agilent Feature Extraction program, version 10.7.3.
Pre-processing of microarrays was performed using Limma's normexp background correction method to avoid negative or zero corrected intensities, followed by between-array normalization using the quantile method to make all array distributions to have the same empirical distribution. Control probes were filtered and then within-array replicate spots were replaced with their average. Pair-wise comparisons between groups were conducted using the Linear Models for Microarray Data (Limma) package Version 3.26.1 from Bioconductor (http://www.bioconductor.org/). The false discovery rate of differentially expressed genes for treatment/control and between-treatment comparisons was based on the Benjamini and Hochberg (BH) procedure. Genes with a score below an adjusted p-value threshold of 0.01 and which also showed a minimum of twofold change in expression between conditions or genotype were selected as significantly differentially expressed genes. Gene annotations were obtained from the Arabidopsis Information Resource (TAIR; http://www.arabidopsis.org/). Functional clustering and analysis was performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) (http://david.abcc.ncifcrf.gov/home.jsp) version 6.7. Differentially expressed genes were compared against gene sets collected from various sources such as publications using the Plant GeneSet Enrichment Analysis Toolkit (PlantGSEA) (http://structuralbiology.cau.edu.cn/PlantGSEA/).
To detect co-regulated gene sets, a cluster analysis of the differentially expressed genes was carried out using data from (Georgii et al., 2012), consisting of microarray data downloaded from NASCArrays (ftp://uiftparabid.nottingham.ac.uk/NASCarrays/By_Experiment_ID/), ArrayExpress (http://www.ebi.ac.uk/microarrayas/ae/), Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/), and The Integrated Microarray Database System (http://ausubellab.mgh.harvard.edu/). Arrays were normalized with Robust Multi-array Average (RMA), and log2 ratio of the mean of treatment and control expressions across biological replicates was computed. Bayesian Hierarchical Clustering was carried out using R package BHC (Cooke et al., 2011) using log2 fold change ±1 as discretization threshold. Gene set enrichment analysis of the co-regulated gene clusters was carried out using StringDB (http://string-db.org/; Szklarczyk et al., 2015).
Isolation of the Thylakoid Membrane and Separation of Protein Complexes
Thylakoid isolation and blue native (BN)-PAGE were performed essentially as described in Järvi et al. (2011). Sodium fluoride was included in thylakoid isolation buffers for samples intended for BN-PAGE, whilst excluded from thylakoids used for spectroscopy analyses (see below). For BN-PAGE, the thylakoid membrane (4 μg Chl) was resuspended into ice-cold 25BTH20G buffer [25 mM BisTris/HCl (pH 7.0), 20% (w/v) glycerol and 0.25 mg ml−1 Pefabloc] to a Chl concentration of 1.0 mg ml−1. An equal volume of 2.0% (w/v) detergent (n-dodecyl β-D-maltoside, Sigma) solution (diluted in 25BTH20G) was added to the sample and thylakoid membrane was solubilized in darkness for 5 min on ice. Traces of insoluble material were removed by centrifugation at 18,000 g at 4°C for 20 min. Prior to loading, the samples were supplemented with a one-tenth volume of Serva Blue G buffer [100 mM BisTris/HCl (pH 7.0), 0.5 M ACA, 30% (w/v) sucrose, and 50 mg ml−1 Serva Blue G].
Spectroscopic Quantitation of PSI and PSII
Room temperature continuous wave electron paramagnetic resonance (EPR) spectroscopy was performed essentially as described in Danielsson et al. (2004) and Suorsa et al. (2015). Measurements were performed at the Chl concentration of 2 mg ml−1.
Photosynthetic Activity Measurements
The Dual-PAM-100 (Walz, http://www.walz.com/) was used for the measurement of PSII quantum yields. Quantum yields of PSII (FV/FM, ΦII, ΦNPQ, and ΦNO) were determined from leaves dark adapted for 30 min before the measurements. Saturating pulse (800 ms, 6000 μmol photons m−2s−1) was applied to determine the maximal fluorescence. Measurements were done in actinic red light of 50, 120, or 500 μmol photons m−2s−1.
The numerical data were subjected to statistical analysis by Student's t-test with statistical significance at the p < 0.05.
Fluctuating Growth Light Only Slightly Modified the Photosynthetic Light Reactions
Accumulating evidence during recent years has demonstrated that sudden, abrupt changes in light intensity threaten particularly PSI, not PSII (Grieco et al., 2012; Suorsa et al., 2012; Allahverdiyeva et al., 2013; Kono and Terashima, 2014). Indeed, quantitation of the functional PSI/PSII ratios from wild-type plants with EPR revealed a PSI/PSII ratio of 1.12 for plants grown under constant light conditions (Suorsa et al., 2015), whereas plants grown under fluctuating light conditions exhibited a clearly lower value, 1.02.
The tlp18.3 plants showed a distinct stunted phenotype upon growth under fluctuating white light and the dry weight of the tlp18.3 plants (12.2 ± 5.7 mg) was markedly decreased as compared to wild type (29.9 ± 4.7 mg; n = 6). This observation prompted us to monitor whether the oligomeric structure of the thylakoid membrane protein complexes of wild-type and tlp18.3 plants grown either under constant or fluctuating light conditions is altered. Malfunction of the PSII repair cycle is often evidenced by a low amount of the most active PSII complexes, the PSII-LHCII complexes, accompanied by a high amount of PSII monomers, which are under the repair cycle (Danielsson et al., 2006). To that end, the BN-PAGE separation of thylakoid protein complexes according to their molecular mass was applied. In line with earlier results (Sirpiö et al., 2007), the tlp18.3 thylakoids accumulated slightly less of the PSII-LHCII complexes under constant light (Figure 1). Similar result was also evident under fluctuating light intensities, the amount of PSII-LHCII being somewhat lower in tlp18.3 plants as compared to wild type. However, no significant differences were observed in heterogeneity of the photosynthetic protein complexes, when wild-type and mutant plants grown either under constant or fluctuating light were compared (Figure 1). A previous report has shown that the maximal PSII quantum yield is not changed in tlp18.3 plants grown under constant growth light conditions as compared to wild type (Sirpiö et al., 2007). In line with this, the maximum quantum yield and effective quantum yields of PSII remained rather similar, when the tlp18.3 and wild-type plants grown their entire life span under fluctuating light were compared (Table 1). Indeed, the PSII activity was only slightly down-regulated in tlp18.3 plants as compared to wild type. Thus, the growth defect shown by the tlp18.3 plants under fluctuating light intensities does not originate from the diminished pool of active PSII complexes.
Figure 1. Accumulation of thylakoid protein complexes in wild-type and tlp18.3 plants. Plants were grown in 8 h light regime either in a photon flux density of 120 μmol photons m−2 s−1 (constant growth light; CL) or 50 μmol photons m−2s−1 for 5 min and 500 μmol μmol photons m−2s−1 for 1 min (FL, fluctuating light). sc. supercomplex. A representative example from three independent biological replications is shown.
Consequences of Fluctuating Growth Light Intensity on Gene Expression
To further characterize plant acclimation to fluctuating light, we performed transcript profiling of the wild-type and tlp18.3 plants grown under constant and fluctuating light intensities and compared the four datasets: (i) wild-type plants grown under fluctuating vs. constant growth light, (ii) tlp18.3 plants grown under fluctuating vs. constant growth light, (iii) tlp18.3 vs. wild-type plants grown under fluctuating light, and (iv) tlp18.3 vs. wild-type plants grown under constant light. Gene enrichment analysis and functional annotation clustering of differentially expressed genes were performed using the DAVID bioinformatic resource (the cutoff was set to logFC > 1 and the adjusted p-value threshold to a minimum of 0.01).
Wild-type plants grown under fluctuating light showed significantly different transcript abundance for 406 genes as compared to wild type grown under constant light, whereas in the tlp18.3 mutant, 321 genes responded differentially to fluctuating light as compared to growth light (Figure 2). When the transcript abundances between the genotypes was compared, 237 genes showed significantly different transcript abundance in tlp18.3 plants compared to wild type when grown under fluctuating light conditions, whereas under constant growth light the number of differentially expressed genes between wild type and the tlp18.3 mutant was 102 (Figure 2). Thus, it can be concluded that the growth light condition altered the number of differentially regulated genes more pronouncedly than the genotype. Moreover, the wild-type plants showed more profound changes at their gene expression level as a response to fluctuating growth light than the tlp18.3 plants.
Figure 2. Venn diagram showing the overlap of significantly differentially regulated genes (logFC > 1) in response to either fluctuating light (FL) as compared to constant growth light (CL) or deficient function of the TLP18.3 protein.
Plants Grown under Fluctuating Light did not Show Differential Abundance of Photosynthesis Related Transcripts
Examination of differentially expressed genes revealed no photosynthesis-related gene ontologies in any of the four datasets analyzed (Tables 2, 3). Indeed, no gene ontologies related to photosynthetic light reactions, Calvin-Benson-Bassham cycle, or biosynthesis of photosynthetic pigments was observed in the gene enrichment analysis. Presumably, regulation of the photosynthetic machinery at transcriptional level does not play an important role during acclimation to relatively mild light intensity fluctuations, being designed such that the total amount of photons hitting the leaf remained nearly unchanged during the 8 h light period, when constant and fluctuating light conditions were compared. Likewise, deficient function of the TLP18.3 protein had only minor effects on transcript abundance of various photosynthesis genes.
Table 2. Classification of significantly differently expressed genes base on gene enrichment analysis of plants grown either under fluctuating light (FL) or constant growth light (CL): (A) Gene enrichment analysis of wild-type plants grown either under fluctuating or constant light; (B) Gene enrichment analysis of tlp18.3 plants grown either under fluctuating or constant light.
Table 3. Classification of significantly differentially expressed genes base on gene enrichment analysis in wild-type and tlp18.3 plants: (A) Gene enrichment analysis of in tlp18.3 plants as compared to wild-type plants grown under fluctuating light (FL); (B) Gene enrichment analysis of in tlp18.3 plants as compared to wild-type plants grown under constant light (CL).
Fluctuating Light Conditions Induced Transcriptional Adjustments in Immunity Related Genes Both in Wild-Type and tlp18.3 Plants
Bioinformatic analysis revealed that the majority of differentially expressed gene ontologies between plants grown under fluctuating and constant light conditions were linked to biotic or abiotic stress responses (Tables 2A,B). In wild type, growth under fluctuating light resulted in decreased transcript abundance within numerous gene ontologies related to plant immunity, as compared to wild type grown under constant light (Table 2A). These genes included mitogen-activated protein kinases (MAPKs) involved in early defense signaling, Toll/Interleukin-1 receptor-nucleotide binding site (TIR-NBS) class resistance (R) proteins mediating effector-triggered immunity (ETI) as well as pathogen related defense proteins, such as plant defensins (Supplementary Table 1). In contrast, the tlp18.3 mutant showed both decreased and increased transcript abundance within gene ontologies related to plant immunity, when fluctuating and constant light grown plants were compared to each other (Table 2B). For example, ankyrin BDA1 (AT5G54610), which is induced by salicylic acid (SA) and is involved in innate immunity (Blanco et al., 2005; Yang et al., 2012) showed cumulative repression in the transcript abundance in response to fluctuating light and deficient function of the TLP18.3 protein. In contrast, plant defensin PDF2.1 (AT2G02120) and defensin-like (AT2G43535) genes, which are activated in response to fungal infection, were induced in tlp18.3 plants under fluctuating light.
With respect to abiotic stress, gene ontologies “response to UV” and “response to light stimulus” were enriched in the transcriptome of tlp18.3 leaves, when plants grown under fluctuating and constant light were compared (Table 2B). For example, increased abundance of transcripts for EARLY LIGHT-INDUCED PROTEIN2 (ELIP2; AT4G14690), which modulates Chl biosynthesis to prevent photo-oxidative stress (Tzvetkova-Chevolleau et al., 2007; Hayami et al., 2015), was observed in the fluctuating-light-grown tlp18.3 plants (Supplementary Table 1). In contrast, no gene ontologies related to light perception showed differential expression in the wild-type plants as a response to fluctuating light (Table 2A). Decreased transcript abundance of gene ontologies associated with lipid localization and lipid transport were also observed as response to fluctuating light specifically in tlp18.3 leaves. Several genes encoding lipid-transfer proteins such as LIPID TRANSFER PROTEIN 3 (LTP3; AT5G59320), which mediates freezing and drought stress in Arabidopsis (Guo et al., 2013), were down-regulated in the tlp18.3 mutant, when plants were grown under fluctuating light as compared to constant growth light (Supplementary Table 1).
When fluctuating-light-grown tlp18.3 and wild-type plants were compared to each other, increased transcript abundance of genes related to the defense mechanisms in the tlp18.3 mutant was again the most prominent result (Table 3A). Enrichment analysis and functional annotation clustering of the differentially expressed gene ontologies in tlp18.3 and wild-type plants also revealed that several gene clusters related to abiotic stresses were differentially expressed in tlp18.3 plants as compared to wild type under fluctuating light. Decreased transcript abundance of gene ontologies “response to light stimulus” and “response to oxidative stress” was observed in the tlp18.3 mutant as compared to wild type. Closer look at the genes among these categories pinpointed that the transcript abundance for cytosolic and chloroplastic COPPER/ZINC SUPEROXIDE DISMUTASES 1 (AT1G08830) and 2 (AT2G28190), respectively, was repressed in tlp18.3 plants as compared to wild type under fluctuating light conditions (Supplementary Table 1).
Finally, when constant-light-grown tlp18.3 and wild-type plants were compared, only a few gene ontologies related to biotic or abiotic stresses were identified (Table 3B). This result is consistent with the postulated role of TLP18.3 specifically during the dynamic light acclimation process, as evidenced by the distinct growth phenotype of the mutant plants under fluctuating light.
Adjustments in Immunity-Related Genes under Fluctuating Light are Linked to Plant Hormones
Plant acclimation to various stresses, including light stress, is regulated by signaling cascades, which include plant hormones as central components (Karpinski et al., 2013; Müller and Munné-Bosch, 2015). In wild-type plants, growth under fluctuating light as compared to constant light resulted in reduced transcript abundance of several genes related to SA signaling cascades (Table 2A). For example, expression of a gene encoding SYSTEMIC ACQUIRED RESISTANCE DEFICIENT 1 (SARD1; AT1G73805), a key regulator of ISOCHORISMATE SYNTHASE 1, a rate-limiting enzyme in pathogen-induced SA biosynthesis (Zhang et al., 2010), was shown to be down-regulated in wild-type plants grown under fluctuating light. Also expression of a gene encoding BENZOIC ACID/SA CARBOXYL METHYLTRANSFERASE 1 (BSMT1; AT3G11480), which synthetizes methyl salicylate (a mobile signal molecule for plant systemic acquired resistance) from SA (Park et al., 2007), was down-regulated in fluctuating light. In line with these results, WALL-ASSOCIATED KINASE 2 (WAK2; AT1G21270) and L-TYPE LECTIN RECEPTOR KINASE IV.1 (LecRK-IV.1; AT2G37710), which are both induced by SA, showed reduced transcript abundance in wild-type plants as response to fluctuating light (He et al., 1999; Blanco et al., 2005) (Supplementary Table 1). Also the tlp18.3 plants grown under fluctuating light showed decreased abundance of gene transcripts related to SA signaling as compared to plants grown under constant light (Table 2B). However, the number of repressed genes was lower in the tlp18.3 mutant as compared to wild type and no differential expression of SARD1 or BSMT1 were observed in tlp18.3 plants as response to fluctuating light (Table 2, Supplementary Table 1). Decreased amount of transcripts related to SA signaling was also evident when tlp18.3 plants grown under constant light were compared to wild type (Table 3B), while no difference in SA signaling was observed between tlp18.3 and wild-type plants grown under fluctuating light (Table 3A). To that end, the fluctuating light condition and to a lesser extent deficient function of the TLP18.3 protein repressed the SA responsive genes.
Similarly, ethylene (ET)- and jasmonate (JA)-related defense pathways showed reduced transcript abundance in wild-type plants grown under fluctuating light as compared to constant light (Table 2A), while in the tlp18.3 mutant no difference was observed in ET/JA defense reactions between the light conditions (Table 2B). It seems that the repression of ET/JA responsive gene expression under fluctuating light is blocked in the tlp18.3 mutants, which became apparent when ET/JA responses between fluctuating light grown tlp18.3 and wild-type plants were compared (Table 3A).
The most prominent alteration in the gene ontology level, when the transcript abundances of constant light grown tlp18.3 and wild-type plants were compared, was an increase in transcripts of six genes encoding CC-type glutaredoxins (ROXY 5, ROXY 11-15) and two of those, ROXY 5 and ROXY 13, were up-regulated in tlp18.3 as compared to wild type also under fluctuating light (Tables 3, 4, Supplementary Table 1). As CC-type glutaredoxins have been suggested to be capable of suppressing the JA and ET-induced defense genes (Zander et al., 2012), a causal connection might exist between expression of JA and ET-responsive genes and differential expression of ROXY genes. It can be concluded that alteration in the gene expression patterns of SA, ET, and JA signaling are taking place during plant acclimation to fluctuating light and that these alterations are strongly affected by the deficient function of the TLP18.3 protein.
Table 4. List of genes which are significantly differentially expressed in tlp18.3 plants as compared to wild type both under fluctuating (FL) and constant light (CL) conditions (logFC > 1).
Phytochrome-Mediated Light Signaling is Likely to be Altered in tlp18.3 Plants
Next, we wanted to further explore which Arabidopsis genes showed a differential expression pattern in the tlp18.3 plants both under constant and fluctuating light conditions. In addition to ROXY5 and ROXY13 located in the endomembrane system, genes encoding cold (DELTA-9 DESATURASE 1) and drought-repressed (DROUGHT-REPRESSED 4) proteins, acid phosphatase (AT4G29270), and two putative membrane transporters (AT5G62730, AT2G16660) showed differential expression in the tlp18.3 mutant. Interestingly, two genes encoding bHLH class phytochrome A-signaling components, LONG HYPOCOTYL IN FAR-RED 1 (HFR1; AT1G02340) and PHYTOCHROME INTERACTING FACTOR 3-LIKE 1 (PIL1; AT2G46970; Fairchild et al., 2000; Salter et al., 2003), showed decreased transcript abundance in tlp18.3 plants as compared to wild type (Table 4). Instead, expression of the gene encoding EARLY FLOWERING 4 (ELF4; AT2G40080), a phytochrome-controlled regulator of circadian clock was induced in the tlp18.3 mutant as compared to wild type. Taken together, the deficient function of TLP18.3 is likely to change the phytochrome-mediated light signaling both under constant and fluctuating light intensities.
Decreased Transcript Abundance of Dark-Induced Genes Suggest that Nitrogen to Carbon and/or Phosphorus to Carbon Ratios Might be Altered in tlp18.3 Plants under Fluctuating Light
Nutrient availability plays an important regulatory role in growth and development of plants, but also cross-talk between nutrient availability and disease resistance exist (Huber, 1980; Hermans et al., 2006). Interestingly, GLUTAMINE-DEPENDENT ASPARAGINE SYNTHASE 1/DARK-INDUCED 6 (ASN1/DIN6; AT3G47340) and DARK-INDUCED 1/SENESCENCE 1 (DIN1/SEN1; AT4G35770) genes showed strong down-regulation in fluctuating light grown tlp18.3 plants as compared to either fluctuating light grown wild type or constant light grown tlp18.3 plants (Supplementary Table 1). ASN1/DIN6 regulates the flow of nitrogen into asparagine, which acts as a nitrogen storage and transport compound in darkness and its gene expression is regulated by the nitrogen to carbon ratio (Lam et al., 1994). DIN1/SEN1, which has been suggested to contribute to enhanced susceptibility to plant viruses, is induced by phosphate starvation and repressed by sugars (Fernández-Calvino et al., 2015). The differential expression of ASN1/DIN6 and DIN1/SEN1 is linked to deficient function of TLP18.3 under fluctuating light but the exact mechanism behind transcriptional repression of these two genes remains to be verified.
Cluster Analysis of Genes whose Expression in Fluctuating Light Requires Functionality of TLP18.3
Finally, to shed light on gene expression changes that depend on the functionality of TLP18.3 under fluctuating light, the expression profiles of genes differentially expressed in wild type but not in tlp18.3 upon growth under fluctuating light were clustered using publicly available datasets (Figure 3). These wild-type specific genes grouped into 13 co-expression clusters, which were further analyzed for enrichment of gene ontology categories (Supplementary Table 2). Clusters 3-13 contained genes with increased transcript abundance in different abiotic stress conditions including salinity and drought as well as methyl viologen (Paraquat; PQ) and the SA analog BTH (Figure 3). Under UV-B stress, in contrast, the expression of these genes was generally down-regulated (Figure 3). This pattern of gene expression was particularly evident within the gene clusters 5, 6, and 9, which showed significant enrichment of gene ontology categories related to plant immunity, such as “response to chitin,” “ethylene-activated signaling pathway,” or “systemic acquired resistance” (Supplementary Table 2). In wild type the genes belonging to clusters 5, 6, and 9 were generally down-regulated, showing a similar pattern to UV-B stress.
Figure 3. Cluster analysis of genes differentially expressed in the wild-type but not in tlp18.3 plants in response to fluctuating light as compared to constant growth light. Bayesian hierarchical clustering of genes, which are significantly differentially regulated (logFC > 1) in wild type under fluctuating light as compared to constant light, is presented. Data sets used include abiotic and biotic stress experiments. Blue and red indicate decreased and increased expression as compared to untreated plants, respectively.
During the past few years evidence has been accumulated concerning the role of photosynthesis in plant immunity. Here, we have provided new insights into the linkage between light acclimation and plant immunity at the level of gene expression as well as addressed the role of the TLP18.3 protein within these processes. Chloroplasts, in addition to their main task in conversion of solar energy into chemical energy, participate in a number of other reactions like biosynthesis of amino acids, hormones, and secondary metabolites as well as cellular sensing of abiotic and biotic stress signals. Indeed, signals originating from the photosynthetic light reactions such as redox state of the electron transfer chain, accumulation of stromal metabolites as well as ROS and reactive electrophilic species are key components of chloroplast retrograde signaling (Fey et al., 2005; Piippo et al., 2006; Queval and Foyer, 2012; Szechyńska-Hebda and Karpiński, 2013; Bobik and Burch-Smith, 2015; Gollan et al., 2015). These signals respond rapidly to changes in perception of light by the two photosystems.
Here, we focused on plants grown under either constant or fluctuating light conditions for their entire life span in order to unravel how the rapid fluctuations in the growth light intensity affect the acclimation processes at the level of nuclear gene expression. In short, neither photosynthesis-related genes nor the photosynthetic protein complexes showed significant alterations as a response to fluctuating light (Figure 1, Tables 1–3). Instead, EPR spectroscopy revealed that the relative amount of functional PSI complexes was lowered in fluctuating light as compared to plants grown under constant light. Most prominently, in wild-type plants fluctuations in growth light suppressed the expression of genes related to defense reactions (Table 2A). Despite the high-light peaks of 1 min, the low-light phase is dominant in our fluctuating light setup. Hence, it is highly likely that decreased transcript abundance of the defense genes in wild-type Arabidopsis under fluctuating light is linked to shade-avoidance and is mediated by plant hormones (Vandenbussche et al., 2005; Wit et al., 2013). The experimental setup, in which the gene expression was studied from plants grown their entire life span either under constant or fluctuating light did not allow us to identify specific immune responses activated by the fluctuations in the growth light intensity. Instead, this experimental setup shed light into late stages of the plant acclimation process, in which a vast number of defense pathways were affected.
Contrary to wild type, in the tlp18.3 mutant the alterations in the overall gene expression pattern, as a response to fluctuating light, were less evident and indeed, the tlp18.3 plants were less capable of turning off the gene expression related to plant immunity under fluctuating light conditions (Table 2B, Figures 2, 3). It is known that the photoreceptor-derived signals activate the shade-avoidance responses and reduce the defense reactions against pathogens and pests to save resources for the growth of the plant (Ballare, 2014). Interestingly, the gene expression of two components of phytochrome-mediated light signaling, HFR1 and PIL1, was shown to be altered in tlp18.3 leaves (Table 4). HFR1 and PIL1 genes are involved in transcriptional regulation pathways downstream of phytochromes, which integrate light and hormonal signals and play a role in shade-avoidance responses (Jiao et al., 2007). Of these, HFR1 also contributes to the crosstalk between light signaling and plant innate immunity (Tan et al., 2015). Based on these results, it is evident that the functionality of TLP18.3 protein modifies the light perception and/or signaling network, and possibly also the signaling related to nutrient availability (Supplementary Table 1). Allocation of resources to defense reactions in the tlp18.3 mutant is likely associated with the lower biomass of mutant plants as compared to wild-type plants under low-light dominant fluctuating light. It should be noted that the tlp18.3 plants also had lower biomass as compared to wild type when grown under high-light dominant fluctuating light with longer, 1 h light pulses (Sirpiö et al., 2007). It remains to be studied whether the growth phenotype of tlp18.3 plants under high-light dominant fluctuating light originates directly from the diminished pool of active PSII complexes. Indeed, duration, frequency, and intensity of fluctuating light regimes have been shown to affect the acclimation responses in Arabidopsis (Alter et al., 2012). To that end, it would be interesting to compare how the gene expression patterns of low-light and high-light dominant fluctuating light conditions differ from each other.
Defective degradation of the D1 core protein of PSII in tlp18.3 plants is a promising system for the search of chloroplast-derived retrograde signals which affect gene expression related to plant immunity. In line with this, low amount of the D1 degrading protease FtsH has been earlier observed to accelerate the hypersensitive reaction in tobacco (Seo et al., 2000). Recently, a link between PsbS-mediated photoprotection and pathogen resistance has also been shown to exist (Göhre et al., 2012; Johansson Jänkänp et al., 2013). Further, as the PSII repair cycle and maintenance of PSI are interconnected (Tikkanen et al., 2014), also PSI and/or PSI electron acceptors might act as a source of retrograde signaling components under fluctuating light. It should be noted that the pool of active PSII was not changed in tlp18.3 plants as compared to wild type under low-light dominant fluctuating light (Table 1) and thus the effect might be indirect. We postulate that the compensation mechanisms activated in the tlp18.3 mutant are likely to alter the chloroplast-derived retrograde signals. Taken together, our results demonstrate that light acclimation and plant immunity are interconnected and the proper repair cycle of PSII plays a key role in the process.
SJ, JI, SK, JS, and FM contributed to acquisition, analysis, and drafting the work, while MS and EA designed the work and contributed to acquisition, analysis, and drafting the work.
Our research was financially supported by the Academy of Finland (project numbers 272424, 271832, and 273870), TEKES LIF 40128/14, the Swedish Research Council, the Swedish Energy Agency, the Knut and Alice Wallenberg Foundation and the Initial Training Networks (ITN) CALIPSO (607607), and PHOTOCOMM (317184).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Microarray and sequencing unit of the Turku Centre for Biotechnology is thanked for assistance with microarray hybridizations. Kurt Ståle, Mika Keränen, Virpi Paakkarinen, Marjaana Rantala, Sanna Rantala, Ville Käpylä, and Saara Mikola are acknowledged for their excellent technical assistance.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.00405
Allahverdiyeva, Y., Mustila, H., Ermakova, M., Bersanini, L., Richaud, P., Ajlani, G., et al. (2013). Flavodiiron proteins Flv1 and Flv3 enable cyanobacterial growth and photosynthesis under fluctuating light. Proc. Natl. Acad. Sci. U.S.A. 110, 4111–4116. doi: 10.1073/pnas.1221194110
Alter, P., Dreissen, A., Luo, F. L., and Matsubara, S. (2012). Acclimatory responses of Arabidopsis to fluctuating light environment: comparison of different sunfleck regimes and accessions. Photosynth. Res. 113, 221–237. doi: 10.1007/s11120-012-9757-2
Baena-Gonzalez, E., and Aro, E. M. (2002). Biogenesis, assembly and turnover of photosystem II units. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 357, 1451–1459. discussion: 1459–1460. doi: 10.1098/rstb.2002.1141
Blanco, F., Garreton, V., Frey, N., Dominguez, C., Perez-Acle, T., Van der Straeten, D., et al. (2005). Identification of NPR1-dependent and independent genes early induced by salicylic acid treatment in Arabidopsis. Plant Mol. Biol. 59, 927–944. doi: 10.1007/s11103-005-2227-x
Cooke, E. J., Savage, R. S., Kirk, P. D., Darkins, R., and Wild, D. L. (2011). Bayesian hierarchical clustering for microarray time series data with replicates and outlier measurements. BMC Bioinformatics 12:399. doi: 10.1186/1471-2105-12-399
Danielsson, R., Albertsson, P. A., Mamedov, F., and Styring, S. (2004). Quantification of photosystem I and II in different parts of the thylakoid membrane from spinach. Biochim. Biophys. Acta 1608, 53–61. doi: 10.1016/j.bbabio.2003.10.005
Danielsson, R., Suorsa, M., Paakkarinen, V., Albertsson, P. A., Styring, S., Aro, E. M., et al. (2006). Dimeric and monomeric organization of photosystem II. Distribution of five distinct complexes in the different domains of the thylakoid membrane. J. Biol. Chem. 281, 14241–14249. doi: 10.1074/jbc.M600634200
de Torres Zabala, M., Littlejohn, G., Jayaraman, S., Studholme, D., Bailey, T., Lawson, T., et al. (2015). Chloroplasts play a central role in plant defence and are targeted by pathogen effectors. Nat. Plants 1:15074. doi: 10.1038/nplants.2015.74
Fernández-Calvino, L., Guzmán-Benito, I., Toro, F. J., Donaire, L., Castro-Sanz, A. B., Ruíz-Ferrer, V., et al. (2015). Activation of senescence-associated Dark-inducible (DIN) genes during infection contributes to enhanced susceptibility to plant viruses. Mol. Plant. Pathol. 17, 3–15. doi: 10.1111/mpp.12257
Fey, V., Wagner, R., Brautigam, K., Wirtz, M., Hell, R., Dietzmann, A., et al. (2005). Retrograde plastid redox signals in the expression of nuclear genes for chloroplast proteins of Arabidopsis thaliana. J. Biol. Chem. 280, 5318–5328. doi: 10.1074/jbc.M406358200
Georgii, E., Salojärvi, J., Brosche, M., Kangasjärvi, J., and Kaski, S. (2012). Targeted retrieval of gene expression measurements using regulatory models. Bioinformatics 28, 2349–2356. doi: 10.1093/bioinformatics/bts361
Göhre, V., Jones, A. M., Sklenár, J., Robatzek, S., and Weber, A. P. (2012). Molecular crosstalk between PAMP-triggered immunity and photosynthesis. Mol. Plant Microbe Interact. 25, 1083–1092. doi: 10.1094/MPMI-11-11-0301
Grieco, M., Tikkanen, M., Paakkarinen, V., Kangasjärvi, S., and Aro, E. M. (2012). Steady-state phosphorylation of light-harvesting complex II proteins preserves photosystem I under fluctuating white light. Plant Physiol. 160, 1896–1910. doi: 10.1104/pp.112.206466
Guo, L., Yang, H., Zhang, X., and Yang, S. (2013). Lipid transfer protein 3 as a target of MYB96 mediates freezing and drought stress in Arabidopsis. J. Exp. Bot. 64, 1755–1767. doi: 10.1093/jxb/ert040
Hayami, N., Sakai, Y., Kimura, M., Saito, T., Tokizawa, M., Iuchi, S., et al. (2015). The responses of Arabidopsis early light-induced protein2 to ultraviolet B, high light, and cold stress are regulated by a transcriptional regulatory unit composed of two elements. Plant Physiol. 169, 840–855. doi: 10.1104/pp.15.00398
He, Z., Cheeseman, I., He, D., and Kohorn, B. D. (1999). A cluster of five cell wall-associated receptor kinase genes, Wak1–5, are expressed in specific organs of Arabidopsis. Plant Mol. Biol. 39, 1189–1196. doi: 10.1023/A:1006197318246
Herbstova, M., Tietz, S., Kinzel, C., Turkina, M. V., and Kirchhoff, H. (2012). Architectural switch in plant photosynthetic membranes induced by light stress. Proc. Natl. Acad. Sci. U.S.A. 109, 20130–20135. doi: 10.1073/pnas.1214265109
Hermans, C., Hammond, J. P., White, P. J., and Verbruggen, N. (2006). How do plants respond to nutrient shortage by biomass allocation? Trends Plant Sci. 11, 610–617. doi: 10.1016/j.tplants.2006.10.007
Huber, D. (1980). “The role of mineral nutrition in defense,” in Plant Disease: An Advanced Treatise: How Plants Defend Themselves, ed J. G. Horsfall (New York, NY: Academic Press), 381–406. doi: 10.1016/B978-0-12-356405-4.50028-9
Järvi, S., Suorsa, M., and Aro, E. M. (2015). Photosystem II repair in plant chloroplasts—regulation, assisting proteins and shared components with photosystem II biogenesis. Biochim. Biophys. Acta 1847, 900–909. doi: 10.1016/j.bbabio.2015.01.006
Järvi, S., Suorsa, M., Paakkarinen, V., and Aro, E. M. (2011). Optimized native gel systems for separation of thylakoid protein complexes: novel super- and mega-complexes. Biochem. J. 439, 207–214. doi: 10.1042/BJ20102155
Johansson Jänkänpää, H., Frenkel, M., Zulfugarov, I., Reichelt, M., Krieger-Liszkay, A., Mishra, Y., et al. (2013). Non-photochemical quenching capacity in Arabidopsis thaliana affects herbivore behaviour. PLoS ONE 8:e53232. doi: 10.1371/journal.pone.0053232
Karpinski, S., Szechynska-Hebda, M., Wituszynska, W., and Burdiak, P. (2013). Light acclimation, retrograde signalling, cell death and immune defences in plants. Plant Cell Environ. 36, 736–744. doi: 10.1111/pce.12018
Kirchhoff, H., Hall, C., Wood, M., Herbstova, M., Tsabari, O., Nevo, R., et al. (2011). Dynamic control of protein diffusion within the granal thylakoid lumen. Proc. Natl. Acad. Sci. U.S.A. 108, 20248–20253. doi: 10.1073/pnas.1104141109
Konert, G., Rahikainen, M., Trotta, A., Durian, G., Salojärvi, J., Khorobrykh, S., et al. (2015). Subunits B'gamma and B'zeta of protein phosphatase 2A regulate photo-oxidative stress responses and growth in Arabidopsis thaliana. Plant Cell Environ. 38, 2641–2651. doi: 10.1111/pce.12575
Kono, M., and Terashima, I. (2014). Long-term and short-term responses of the photosynthetic electron transport to fluctuating light. J. Photochem. Photobiol. B. 137, 89–99. doi: 10.1016/j.jphotobiol.2014.02.016
Lam, H. M., Peng, S. S., and Coruzzi, G. M. (1994). Metabolic regulation of the gene encoding glutamine-dependent asparagine synthetase in Arabidopsis thaliana. Plant Physiol. 106, 1347–1357. doi: 10.1104/pp.106.4.1347
Park, S. W., Kaimoyo, E., Kumar, D., Mosher, S., and Klessig, D. F. (2007). Methyl salicylate is a critical mobile signal for plant systemic acquired resistance. Science 318, 113–116. doi: 10.1126/science.1147113
Piippo, M., Allahverdiyeva, Y., Paakkarinen, V., Suoranta, U. M., Battchikova, N., and Aro, E. M. (2006). Chloroplast-mediated regulation of nuclear genes in Arabidopsis thaliana in the absence of light stress. Physiol. Genomics 25, 142–152. doi: 10.1152/physiolgenomics.00256.2005
Rodríguez-Herva, J. J., González-Melendi, P., Cuartas-Lanza, R., Antúnez-Lamas, M., Río-Alvarez, I., Li, Z., et al. (2012). A bacterial cysteine protease effector protein interferes with photosynthesis to suppress plant innate immune responses. Cell. Microbiol. 14, 669–681. doi: 10.1111/j.1462-5822.2012.01749.x
Seo, S., Okamoto, M., Iwai, T., Iwano, M., Fukui, K., Isogai, A., et al. (2000). Reduced levels of chloroplast FtsH protein in tobacco mosaic virus-infected tobacco leaves accelerate the hypersensitive reaction. Plant Cell 12, 917–932. doi: 10.1105/tpc.12.6.917
Sirpiö, S., Allahverdiyeva, Y., Suorsa, M., Paakkarinen, V., Vainonen, J., Battchikova, N., et al. (2007). TLP18.3, a novel thylakoid lumen protein regulating photosystem II repair cycle. Biochem. J. 406, 415–425. doi: 10.1042/BJ20070460
Suorsa, M., Järvi, S., Grieco, M., Nurmi, M., Pietrzykowska, M., Rantala, M., et al. (2012). PROTON GRADIENT REGULATION5 is essential for proper acclimation of Arabidopsis photosystem I to naturally and artificially fluctuating light conditions. Plant Cell 24, 2934–2948. doi: 10.1105/tpc.112.097162
Suorsa, M., Rantala, M., Mamedov, F., Lespinasse, M., Trotta, A., Grieco, M., et al. (2015). Light acclimation involves dynamic re-organization of the pigment-protein megacomplexes in non-appressed thylakoid domains. Plant J. 84, 360–373. doi: 10.1111/tpj.13004
Szklarczyk, D., Franceschini, A., Wyder, S., Forslund, K., Heller, D., Huerta-Cepas, J., et al. (2015). STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 43, D447–D452. doi: 10.1093/nar/gku1003
Tan, C. M., Li, M. Y., Yang, P. Y., Chang, S. H., Ho, Y. P., Lin, H., et al. (2015). Arabidopsis HFR1 is a potential nuclear substrate regulated by the Xanthomonas type III effector XopD(Xcc8004). PLoS ONE 10:e0117067. doi: 10.1371/journal.pone.0117067
Tikkanen, M., Grieco, M., Kangasjärvi, S., and Aro, E. M. (2010). Thylakoid protein phosphorylation in higher plant chloroplasts optimizes electron transfer under fluctuating light. Plant Physiol. 152, 723–735. doi: 10.1104/pp.109.150250
Tikkanen, M., Mekala, N. R., and Aro, E. M. (2014). Photosystem II photoinhibition-repair cycle protects Photosystem I from irreversible damage. Biochim. Biophys. Acta 1837, 210–215. doi: 10.1016/j.bbabio.2013.10.001
Tyystjärvi, E., and Aro, E. M. (1996). The rate constant of photoinhibition, measured in lincomycin-treated leaves, is directly proportional to light intensity. Proc. Natl. Acad. Sci. U.S.A. 93, 2213–2218. doi: 10.1073/pnas.93.5.2213
Tzvetkova-Chevolleau, T., Franck, F., Alawady, A. E., Dall'Osto, L., Carrière, F., Bassi, R., et al. (2007). The light stress-induced protein ELIP2 is a regulator of chlorophyll synthesis in Arabidopsis thaliana. Plant J. 50, 795–809. doi: 10.1111/j.1365-313X.2007.03090.x
Wit, M., Spoel, S. H., Sanchez-Perez, G. F., Gommers, C. M., Pieterse, C. M., Voesenek, L. A., et al. (2013). Perception of low red: far-red ratio compromises both salicylic acid-and jasmonic acid-dependent pathogen defences in Arabidopsis. Plant J. 75, 90–103. doi: 10.1111/tpj.12203
Wu, H. Y., Liu, M. S., Lin, T. P., and Cheng, Y. S. (2011). Structural and functional assays of AtTLP18.3 identify its novel acid phosphatase activity in thylakoid lumen. Plant Physiol. 157, 1015–1025. doi: 10.1104/pp.111.184739
Yang, Y., Zhang, Y., Ding, P., Johnson, K., Li, X., and Zhang, Y. (2012). The ankyrin-repeat transmembrane protein BDA1 functions downstream of the receptor-like protein SNC2 to regulate plant immunity. Plant Physiol. 159, 1857–1865. doi: 10.1104/pp.112.197152
Zander, M., Chen, S., Imkampe, J., Thurow, C., and Gatz, C. (2012). Repression of the Arabidopsis thaliana jasmonic acid/ethylene-induced defense pathway by TGA-interacting glutaredoxins depends on their C-terminal ALWL motif. Mol. Plant. 5, 831–840. doi: 10.1093/mp/ssr113
Zhang, Y., Xu, S., Ding, P., Wang, D., Cheng, Y. T., He, J., et al. (2010). Control of salicylic acid synthesis and systemic acquired resistance by two members of a plant-specific family of transcription factors. Proc. Natl. Acad. Sci. U.S.A. 107, 18220–18225. doi: 10.1073/pnas.1005225107
Keywords: Arabidopsis thaliana, defense, photosynthesis, photosystem II repair cycle, thylakoid lumen, transcriptomics
Citation: Järvi S, Isojärvi J, Kangasjärvi S, Salojärvi J, Mamedov F, Suorsa M and Aro E-M (2016) Photosystem II Repair and Plant Immunity: Lessons Learned from Arabidopsis Mutant Lacking the THYLAKOID LUMEN PROTEIN 18.3. Front. Plant Sci. 7:405. doi: 10.3389/fpls.2016.00405
Received: 15 December 2015; Accepted: 16 March 2016;
Published: 31 March 2016.
Edited by:Julian Eaton-Rye, University of Otago, New Zealand
Reviewed by:Masaru Kono, The University of Tokyo, Japan
Tiago Toscano Selão, Nanyang Technological University, Singapore
Copyright © 2016 Järvi, Isojärvi, Kangasjärvi, Salojärvi, Mamedov, Suorsa and Aro. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Eva-Mari Aro, firstname.lastname@example.org