Acclimatization of Photosynthetic Apparatus of Tor Grass (Brachypodium pinnatum) during Expansion

Abstract

The aim of this study was to understand the acclimatization mechanisms of photosynthetic apparatus in Brachypodium pinnatum (L.) P. Beauv grass during its expansion. Twelve populations differentiated by age: young (30–50 years old), intermediate age (ca. 100 y) and old (>300 y) were studied. It was confirmed that the decrease of the number of genotypes as a result of environmental stress and competition were reflected in changes in chlorophyll fluorescence (ChlF) parameters. The old stands were dominated by a few genotypes which seem to be the best acclimatized to the self-shading/competition by lowering their photosynthetic performance during light-phase of photosynthesis. On the other hand, the 'high-speed' photosynthetic rate observed in the young populations can be seen as acclimatization to very adverse conditions. Our results clearly confirm that ChlF is a powerful method of inferring physiological mechanisms of the expansion of tor grass. The Principal Component and Redundancy Analyses, followed with k-means classification, allowed to find the differentiation of groups of distinct ChlF parameters and enabled us to relate them to changes in genotypic diversity of populations. We conclude that the plastic morphological and physiological response to changeable habitat light conditions with its optimum in half-shade refers to its forest-steppe origin.

About the Authors:

Wojciech Bąba

* E-mail: [email protected]

Affiliation: Department of Plant Ecology, Institute of Botany, Jagiellonian University, Lubicz 46, 31–512, Kraków, Poland

Hazem M. Kalaji

Affiliation: Department of Plant Physiology, Warsaw University of Life Sciences SGGW, Nowoursynowska 159, 02–776, Warsaw, Poland

Agnieszka Kompała-Bąba

Affiliation: Department of Botany and Nature Protection, University of Silesia, Jagiellonska 28, 40–032, Katowice, Poland

Vasilij Goltsev

Affiliation: Department of Biophysics and Radiobiology, Faculty of Biology, St. Kliment Ohridski University of Sofia, 8 Dr. Tzankov Blvd., 1164, Sofia, Bulgaria

Introduction

Native species, similar to invasive alien (IAS) ones, can also have negative ecological and economic impact. They spread within their natural range attaining in some cases extreme abundances and exert effect on native vegetation [1]. Brachypodium pinnatum (L.) P. Beauv. belongs to the 'native invaders', which naturally occur in low abundance within calcareous grasslands. Its expansion, however, heavily reduces the biodiversity of calcareous grasslands [2,3] and often ends with formation of nearly monodominant stands [4]. The spread of false brome in Europe was triggered at the end of the 1960s as a result of the abandonment of traditional land use such as grazing and mowing [5,6]. Deposition of airborne nitrogen and phosphorus is considered to be a second possible reason for its expansion in Western Europe [7,8]. Numerous aspects of biology have been reported to date including the plastic response to light conditions [8], biomass and nutrient allocation, nutrient cycling [9–11], mycorrhizal colonization [12] karyotype evolution [13] and the population genetic diversity of this species along colonization gradients or under different management [6,14,15]. Moreover, several attempts have been made to control the Brachypodium expansion [16,17]. However, there is still a need for a quick, easy-to-perform and low cost method of monitoring the viability of large numbers of individuals, aimed at assessing its physiological acclimatization to changes in environmental conditions during expansion It is crucial in elaborating the proper management schemes.

Chlorophyll fluorescence (ChlF) measurements, which have been frequently performed [18–20], could be a method that meets these requirements and enables to assess different stress effects on photosynthetic apparatus using this technique, particularly on Photosystem II (PSII) status and linear electron transport rate. ChlF is a naturally occurring phenomenon, characteristic to all photosynthetic organisms. The ca 1–8% of the sun’s energy that is not used to drive photosynthesis is dissipated as heat radiation or re-emitted as light photons [21]. The analysis based on high time-resolution measurements of the ChlF transient represents a method for gaining detailed information about PSII photochemical activity, electron transport events and the different regulatory processes. Fast ChlF kinetics data are derived from the time dependent increase in fluorescence intensity achieved after application of bright light to a dark adapted sample. The resulting curve is called the Kautsky curve or ChlF transient [19].

The fluorescence parameters obtained in this way, called OJIP-test, allow the quantification of the stepwise flow of energy through PSII, using input data from the fluorescence transient and are formulated with a simplified model of the energy fluxes incorporating the parameters that define each type of flux [22,23]. The energy fluxes consist of an absorbed flux (ABS), trapping flux (TR), electron transport flux (ET) and the flux defining the dissipation of non-trapped energy as heat (DI)–a flux quantifying the reduction of Photosystem I (PSI) end acceptor (RE) was later introduced [19,24–26].

The aim of this research was to understand the mechanisms of acclimatization of photosynthetic apparatus during expansion of tor grass populations. Particularly, the following questions were asked: (i) How the ChlF parameters, related to different aspects of PSII functioning, change with habitat age and (ii) is there any relationship between genotypic/genetic diversity and ChlF parameters. It was expected that ChlF parameters can be used in a quick assessment of genotype acclimatization and can be helpful in management planning intended to prevent the future expansion of this species.

Materials and Methods

Species

Brachypodium pinnatum (L.) P. Beauv (subsequently referred to as Brachypodium) is a rhizomatous, perennial grass with widespread distribution in temperate regions of Northern Hemisphere. It forms clones, usually up to 1.5 (max 4) m in diameter with ca. 80% short rhizomes (1–10 mm) and dense tillers in form of 'clumps of shoots' (3–5 cm diameter) connected with longer rhizomes (100–200 mm) [15]. This enables it the long-term occupation of a given area and quick lateral spread. Brachypodium is both a diploid (2n = 18) and allotetraploid species (2n = 28) [27,28]. The stem height varies from 35 to 120 cm. It starts flowering at the end of June and set seeds from late July to September. The mean seed production on the study area varied from 2 to 45 per inflorescence and from 600 to 5000 per 1 m² and was lowest on the old grasslands.

Study area

The study area lies in Cracow-Częstochowa Upland, southern Poland. It consists of the system of the Upper Jurassic limestone valleys and hills on the plateau between them. The annual average temperature is 7.5°C in vegetation period and mean annual precipitation reaches 773 mm. Brachypodium occurs on the S, SW and SE slopes, and sometimes forms nearly monodominant stands [15,29]. 12 Brachypodium-dominated calcareous grasslands were selected for establishment of the study plots (250 m² 50 x 50 m). The sites were classified according to grassland age and genotypic diversity of Brachypodium stands. Based on cadastral maps (1791–1970) and aerial photographs (1957–2009), the sites were classified as young– 30–50 years old: Racławice [19°40′ 37.06′′E, 50°12′ 02.10′′N], Powroźnikowa 1 [19°40′08.48′′E, 50°11′44.06′′N], Powroźnikowa 2 [19°40′04.53′′E, 50° 12′02.39′′N], Willisowe Skały [19°42′30.76′′E, 50°11′23.49′′N], intermediate age–ca. 100 years old: Dolina Kobylańska [19°45′42.06′′E, 50°9′27.33′′N], Skala Żytnia [19°47′54.96′′E, 50°11′14.72′′N], WielkieSkały [19°48′20.17′E 50°11′23.49′′N], Bolechowice [19°46′58.08′′E, 50°9′11.28′′N], and old>300 years: Grodzisko [19°49′46.55′′E, 50°13′36.39′′N], Grodzisko-Onobrychis [19°49′36.86′′E, 50°13′40.54′′N], Dolina Będkowska [19°44′21.28′′E, 50°8′49.44′′N] and Dolina Kluczwody [19°49′7.04′′E, 50°9′49.94′′N]. Moreover, according to the results of the previous population genetic analyses based on AFLP markers, it was possible to relate age to genotypic richness (G) and the percentage of polymorphic loci (PPL), which decreased with the habitat age [6,15]. The other criteria for classification of Brachypodium populations were the differences in genotypic diversity and morphological parameters of its clones. According to the results of previous population genetic analyses based on AFLP markers, it was possible to relate age to genotypic richness (G) and the percentage of polymorphic loci (PPL), which both decreased with the habitat age [6,15]. Moreover, the average shoots and 'clumps of shoots' density/1m²[9],the number of leaves in clumps of shoots/1m^2, (given in S1–S3 Figs), roots and rhizome dry biomass increased with habitat age S4 and S5 Figs). In contrast, the mean density of generative shoots/1 m² (S6 Fig) and the average seed production/1m² [15] were lowest on the old grasslands. We found this as a sign of a competitive exclusion among genotypes during the Brachypodium expansion [15]. The Authors confirmed that the field studies did not involved endangered or protected species. The sampling sites were located on areas with no specific permissions required for locations and activities.

Climatic data

The climatic data: monthly precipitation and mean air temperatures for young, intermediate and old populations were assessed on the basis of climatic service, http://klimat.icm.edu.pl/serv_climate.php (S1 Table).

Soil and plant analyses

Five soil samples were collected from 0–30 horizons at each site for analysis of the physicochemical parameters. The soils samples were pooled, mixed, air-dried, and then grounded and sieved through 2 mm mesh. The particle size fractions (sand, silt and clay) were determined by sieving and sedimentation method (Prószyński method). Soil pH was measured electrometrically in water suspension following extraction with 1 M KCl and H₂O (1:2.5). Organic matter content was estimated as loss on ignition (LOI) at 550°C of soil samples dried at 105°C for 12 h and expressed as a percentage of dry weight. Total N content was assessed applying the Kjeldahl method using Automatic Kjeldahl Digestion Units and UDK 129 Kjeldahl DKL Distillation Unit (VELP Scientifica, Italy). The total content of Ca, Mg, Mn, Cu and Fe in the soil and plant material was determined by flame atomic absorption spectrometry (Varian Spectra AA 330) after hot digestion of 3 mg of soil in mixture of 65% HNO₃ and 35%H₂O₂ (8:2) with ETHOS ONE microwave system (Milestone, S.r.L., Italy) and exchangeable forms of metals in soil were estimated in 0.5 M HCl soil extracts [30].

Leaf morpho-anatomical traits

Fully developed leaf blades from each of the populations were collected randomly in order to compare the morpho-anatomical traits among individuals of Brachypodium from the old and young populations. The leaf width (LW), height (LH), leaf area (LA), specific leaf area (SLA) were measured according to standard methods [31]. Anatomical measurements of leaves were done on slides, prepared following the standard method for fluorescence microscopy. Cross-sections of leaves were analyzed with Nikon Eclipse with a DS-Fi1-U2 camera (Ni-U, Nikon Co., Japan) to acquire microscopic images at 10× 20× magnification. The thickness of the leaf blades, width, height and area of the central vascular bundle of leaf, number of the sclerenchyma strands on the adaxial and abaxial sides of the tiller leaf were compared among the leaves from old and young populations.

Measurement of Chl a fluorescence

The first fully developed leaves were collected at each location on July 2014, the optimum period for vegetative development of the Brachypodium [8]. The sampling was performed between 6–11 a.m. on cloudy days. 100 fully developed leaves per site were collected from randomly chosen plants. The leaves were put into the paper envelopes, sealed, placed in cooler bags and transported immediately to the laboratory. ChlF measurements were performed on the middle part of abaxial leaf blades away from the main leaf vein after additional dark adaptation (30 min) in a dark room using leaf clips. Fluorescence measurements were performed with the PocketPEA fluorimeter (Hansatech Instruments, King's Lynn, Norfolk, UK). For induction of fluorescence red actinic light was used (wavelength at peak 650 nm; spectral line half-width 22 nm) with the intensity of 3500 μmol m –² s –¹, and 1 second of transient fluorescence was recorded [26,32–34]. The fluorescence signal was collected with a maximum frequency of 10⁵ points s –¹ (each 10 μs) within 0–0.3 ms, after which the frequency of recording gradually decreased, collecting a total of 118 points within 1 s. ChlF transient data were used to calculate basic parameters and the parameters needed for the OJIP-test (Table 1). The F_O level was measured as the fluorescence at 50 μs. The collected data were used for the calculation of basic parameters, while the fluorescence intensities determined at O-50 μs, J-2ms, I-30 ms and maximum fluorescence, P ~ 300 ms (F_M)were used for the calculation of the OJIP test parameters (Table 1)[19,35,36]. To visualize the K and L bands, the collected data points were double normalized as relative variable fluorescence between points O-I, O-J and O-K [37]. Then, the kinetic differences between the old grasslands (treated as control) vs. interm and young were calculated (Table 1). This procedure helped to reveal bands that are normally hidden between the O and P steps on relative variable fluorescence (Table 1) [38].

[Figure omitted. See PDF.]

Table 1. Summary of measured and calculated Chl a fluorescence parameters.

https://doi.org/10.1371/journal.pone.0156201.t001

Determination of Chl a and b content

The leaf greenness index as the average of five readings for each of 30 leaves per site was obtained using a portable chlorophyll meter (SPAD-502 Konica-Minolta, Japan). The SPAD recordings were obtained from a leaf disc of area 169.72 mm^2. After recording, the plant tissue was stored in a glass tube containing 5 mL DMSO (96% Dimethyl sulfoxide) [39]. The test tubes were incubated at 70°C for 48 hours. After cooling the extract in the dark, 3 mL aliquot was analyzed spectrophotometrically at 470, 647 and 663 nm wavelength with the DR 5000 spectrophotometer (Hach Lange). The chlorophyll a (Chl a) and b (Chl b) content was determined according to the following formulas [40]:(1)(2)where A is the absorbance of wavelength, after the correction for scattering at 750 nm and D is the optical thickness of the cuvette [41]. Then, the chlorophyll content per unit leaf area (mg per m²) was calculated.

Statistical analyses

The values shown in the tables represent the mean of (i) all fluorescence parameters, (ii) the values obtained in OJIP test, (iii) the content of selected mineral compounds in the soil and in the plants and (iv) leaf traits, were calculated for all habitat age. The differences in all these variables among the habitats of different age were tested with one-way ANOVA. Prior to the analyses, the assumption of ANOVA: (i) homogeneity and (ii) normality, symmetry of distribution and outliers were visually assessed and tested with (i) Levenne and (ii) QQ plots, boxplots and Shapiro-Wilk W tests [42]. The Welch ANOVA was used in the case of violation of homogeneity assumption. If the F-test was significant, a pairwise comparison of means was calculated using Tukey's test [42]. To visualize the divergence in overall changes of Chl a fluorescence traits among the grasslands of different age, a Principal Component Analysis (PCA) was performed. To relate the fluorescence parameters to genetic diversity and leaf Chl content, redundancy analysis (RDA) was performed. Genetic diversity: percentage of polymorphic loci (PPL), number of genotypes (G) and distribution of frequency of genotypes (Pareto, beta), was assessed using AFLP markers (details in Bąba et al. [6,15]. The Pareto index was high when many genotypes of comparable size occurred in population, and the lowest, when one of the few genotypes dominated the population. It could be used as a proxy value of competition among genotypes [43]. The significance of these variables was tested with the Monte Carlo permutation test (N = 999 permutations). In order to reveal the pattern of changes of ChlF parameters during expansion of Brachypodium, the classification of Chl a fluorescence parameters was performed with k-means clustering. The optimal number of groups was estimated based on the Calinski-Harabasz criterion [44]. All the statistical calculations were performed with R 3.2.0 packages MASS, stats, agricolae and vegan [45].

Results

Soil conditions

In field experiments, as opposed to laboratory ones, the crucial thing is the careful selection of study locations to make them as comparable as possible in terms of all factors, beside the ones of interest (i.e. chlorophyll fluorescence). The sites under study did not differ significantly in soil physico-chemical parameters: skeleton, soil particle fractions percentage content, soil reaction and most of micro- and macroelements content. It was only a slight (but not significant) increase in total or exchangeable forms of K, P, Fe, Cu and decrease in total Mg with grassland age (Table 2). Additionally, the mean content of selected macro—and microelements in the Brachypodium tissue did not show any pattern (Table 2).

[Figure omitted. See PDF.]

Table 2. The physico-chemical parameters of soil and plant tissue concentration of selected micro- and macroelements.

https://doi.org/10.1371/journal.pone.0156201.t002

Variation of morpho-anatomical and chemical characteristics of leaves of Brachypodium

The significant differences were found in morphometric traits of leaf blades during the expansion of Brachypodium. The plants from old populations had higher values of leaf dry mass (LDM), leaf area (LA), leaf length (LL) and leaf width (LW) than those which came from interm and young populations of Brachypodium (Table 3). Individuals from old populations had significantly higher number of leaves per ramet (= rooted shoot sensu Falińska et al [46]), than in the other populations (11.9±2.5 vs. 10.2±4.1 and 10.1±4.1, Welch ANOVA, F = 22.9, p<0.001). However, there were no significant differences in specific leaf area (SLA) among the studied populations.

[Figure omitted. See PDF.]

Table 3. The leaf traits, leaf chlorophyll content and parameters of the Chl a fluorescence tor grass Brachypodium pinnatum from population of different age: young (30–50 years old), intermediate age (ca. 100 years) and old(> 300 years).

The means of for each populations are presented. and means±SE for the particular age classes were compared with ANOVA, followed with Tukey post-hoc test. Values with the same letters are not significantly different at p < 0.05 level. Abbreviations: LA—leaf area, LDM—Leaf dry matter content, SLA—specific leaf area, LL—leaf length, LW—leaf width. Other abbreviations in Table 1. Locality: Rac–Racławice, Pow1 –Powroznikowa 1, Pow2 –Powroznikowa 2, WilSk–Wilisowe Skały, Kob–Dolina Kobylanska, Żyt–Zytnia, WielSk–Wielskie Skały, Bolech–Bolechowice, Gro_Onob–Grodzisko Onobrychis, Bedk–Dolina Bedkowska, DKluc–Dolina Kluczwody, Gro–Grodzisko.

https://doi.org/10.1371/journal.pone.0156201.t003

The mean thickness of the tiller leaf blades of plants from young populations and width of central rib were lower than those from old ones (224 vs. 269 μm and 170 vs. 238μm respectively; Fig 1). Moreover, the differences were found in width (84 vs. 114μm) height (123 vs. 129μm) and area (8574 vs 10050μm²) of the central vascular bundle of leaf, the number and shape of bulliform cells (6–8 vs. 5–6), the number of the sclerenchyma strands on adaxial (5–7 vs. 4–5) and abaxial (4–5 vs. 6–7) sides of the tiller leaf when compared among the leaves of individuals from old and young populations (Fig 1). The chloroplasts in the inner part of mesophyll aligned in vertical columns along the plant cell walls, which is a well-known mechanism of avoidance of photodamage in plants [47]; Fig 1). However, the distribution of chloroplasts in leaves from young populations was more uniform across mesophyll, while those from old ones had bigger chloroplasts located close to the upper and lower epidermis (Fig 1).

[Figure omitted. See PDF.]

Fig 1. Brachypodium pinnatum leaf anatomy from fluorescence microscopy.

Transverse cross section of leaf blade of individuals from old (A, C) and young (B, D) populations. Notice the differences in number and size of bulliform cells occurred on the adaxial side of leaf blade. Transverse section in the midrib at median level (A, B). The differences in (i) thickness and shape of the leaf blades in zone of the central rib, width of the central rib of tiller leaf, (ii) surface, height and width of the central vascular bundle and(iii) number of the sclerenchyma strands on abaxial side of leaf and (iv) distribution of are visible. Detailed leaf measurements were presented in the Results. A, B—20× magnification, C, D—10× magnification. Bars on each of the pictures indicate 100 μm.

https://doi.org/10.1371/journal.pone.0156201.g001

Leaves of individuals from old and intermediate age populations had significantly higher Chl a, and Chl b content, than young ones. However, it was only slightly higher Chl a/b ratio in young populations (3.04 vs. 2.84 in interm and 2.88 in old populations respectively; Table 3). The increase in total leaf Chl content with population age was also reflected in pattern of leaf greenness index measured by SPAD (Table 3).

Chl a fluorescence changes along habitat age gradient

The apparent differences in Brachypodium leaf traits were reflected in the state of Photosystem II, measured by ChlF. The mean prompt ChlF (OJIP) curves both for particular populations and for populations of different age differed both in the shape and amplitude, especially among the young vs. interm and old populations (Fig 2). This suggests possible differences in energy fluxes at the donor as well as at the acceptor side of PSII [23,41]. Moreover, it was also confirmed by mean values of measured and calculated fluorescence parameters (Table 3). Significant differences in formal characteristics of ChlF rise were found: F_O, was significantly lower in interm populations, F_K and F_Jand its values increased from young to old populations.

[Figure omitted. See PDF.]

Fig 2. The chlorophyll fluorescence (ChlF) transients from dark-adapted leaves of expansive grass Brachypodium pinnatum from young (30–50 years), intermediate age (ca. 100 years) and old (>300 years) calcareous grasslands.

The results were plotted on logarithmic scale from F_O(50 μs) to 1 s. A. The fluorescence curves for all 12 populations under study divided into age classes. The time points important for the calculation of JIP test were marked: O–fluorescence intensity recorded at F_O (50 μs), L–at 150 μs, K–at 300 μs J–2 ms, I–at 30 ms, P–maximum fluorescence intensity (F_M) at ca. 1 s. B. The curves of average ChlF values for each age classes.

https://doi.org/10.1371/journal.pone.0156201.g002

The fluorescence rise at the OJ phase was light dependent and provided the information about antenna size and connectivity of PSII reaction centers [48]. The significant peak on the double-normalized fluorescence curve was observed at J (V_J) in young populations(Table 3). It was more pronounced when the course of the Fl rise was expressed as a difference between the analyzed curves. As a referent curve, the normalized fluorescence transient for plants from old populations was accepted as it was assumed that it as an end-point on the expansion gradient (Fig 3).

[Figure omitted. See PDF.]

Fig 3. Variation in parameters reflecting morphological, chemical and photosynthetic conversions in PSII in leaves of Brachypodium plants from population with different age.

The values of each parameter are normalized to values corresponding parameter from old-age population.

https://doi.org/10.1371/journal.pone.0156201.g003

This indicated a slight limitation of electron transport from Q_A to Q_B as indicated by the values 1-V_J, a probability of trapped PSII electron transfer from reduced Q_A to Q_B [41]. The double normalizations of fluorescence curves enabled us also to reveal the less pronounced positive L (V_L) and K (V_K) bands of fluorescence rise in young and interm populations as compared to old ones (Table 3, Fig 4). If the curves were normalized at O and I points (Fig 4) the main two bands appeared in “young–old” curve reflecting the well pronounced J band and well expressed K peak. For the “interm–old” differential curve the K peak at 0.3 ms became dominant. Separately, the effect of the population age on L—and K-bands could be better expressed in the differential curves at double normalization at time intervals “0.02–0.3 ms” and “0.02–2.0 ms”, respectively (see Fig 4).

[Figure omitted. See PDF.]

Fig 4. The chlorophyll fluorescence transients from dark-adapted leaves of expansive grass Brachypodium pinnatum.

A. ChlF curves represent the increase of the relative variable fluorescence in young and intermediate age population relative to old ones (treated as control) between points F_O and F_K; ΔW_OK = (F_t-F_O)/(F_K-F_O); L–band. B. ChlF curves represent the increase of the relative variable fluorescence in young and intermediate age population relative to old (control) ones between points F_O and F_J (ΔW_OJ = (F_t-F_O)/(F_J-F_O); K-band, C. The curves represent the increase of the relative variable fluorescence in young and intermediate age population relative to old (control) ones between points F_O and F_I (ΔW_OI = (F_t-F_O)/(F_I-F_O).

https://doi.org/10.1371/journal.pone.0156201.g004

Increasing positive peaks in these time regions in plants from interm and from young populations as compare to old population reflected more ungrouped PSII units in photosynthetic membranes and less activity of PSII oxygen evolving complex (OEC). The J–P rise of the prompt ChlF curve was attributed to the thermal phase of the fluorescence transient. It reflected a reduction in the electron transport chain [49,50] and, represented the electron transport from Q_A beyond PSI. The V_I values increased significantly with age of populations (Table 3).

The average values of maximum quantum efficiency of PSII φ_Po for all the populations under study was lower than typical for healthy plants in all populations (ca 0.83)[51] and ranged from 0.73 to 0.77 and were slightly, but significantly reduced in the young populations (Table 3; Fig 2; Fig 3). Moreover, the maximum fluorescence value F_M, and the variable part of the ChlF, F_v were highest in old populations (Table 3; Fig 2; Fig 3). The parameters, reflecting the size of electron acceptor pools available on the reducing side of Photosystem II (S_M)and at PSI acceptors (N) were lowest in the old populations. Moreover, the significant decrease in the standardized area (S_M) was observed with age of population (Table 3; Fig 2; Fig 3).

The value of RC/CS₀ (reflecting the relative density of active PSII reaction centers) was significantly lower (of about 40%) in plants from young as compare to old populations. In contrast of this, all parameters presenting the energy fluxes in one active reaction center (ABS/RC, TR₀/RC, ET₀/RC, RE₀/RC and DI/RC) were visibly increased in the young populations. This increase was especially expressed for parameters correlation with the electron transfer site within PSI–from PQH₂ to PSI end acceptors: RE₀/RC, φ_Ro and δ_Ro (Fig 3; Table 3). The most sensitive to the population age parameter was the total Performance Index (PI_total), which combines the efficiency of energy conversion both in PSII and PSI, significantly decreased with age of populations.

The Principal Component Analysis (PCA, Fig 5; Table 4, S2 Table) was used in order to reduce the multivariate fluorescence data into the few principal components and to find the pattern of changes in fluorescence parameters during the expansion of Brachypodium. The classification of the ChlF parameters with the k-means method, based on the Calinski-Harabasz criterion, revealed n = 3 optimal number of groups, which were separated by the first and second PCA axes (Fig 5). The three first PCA axes explained 91.35% of variation in the dataset. The first axis explained 54.45% variation and clearly separated three of four young populations from the interm and old population with higher value of relative fluorescence at point I (V_I). The PCA confirmed that the young populations were characterized by higher values of ChlF related to the plastoquinone (PQ) size pool (S_M), higher electron transport rate and overall PSII performance (TR₀/RC, ET₀/RC,PI_ABS, PI_total), a higher number of Q_A turnover (N) and higher probability with which an electron from the intersystem carriers moves to reduce end electron acceptors at the PSI acceptor side (RE₀/RC).

[Figure omitted. See PDF.]

Fig 5. Principal Component (PCA) ordination diagram displaying the relationships between the Chl a fluorescence measurements and parameters.

The classification of Chl a fluorescence parameters were performed with k-means clustering. The optimal number of groups (3) was estimated based on Caliński-Harabasz criterion. The result of classification was superimposed onto the ordination diagram. The ellipses represent 95% confidence intervals for the groups

https://doi.org/10.1371/journal.pone.0156201.g005

[Figure omitted. See PDF.]

Table 4. The results of Principal Component Analysis (PCA) and Redundancy Analysis (RDA).

The correlations of variables with the first three principal components (PC1-PC3, RDA1-RDA3) and cannonical correlations of leaf chlorophyll content and genetic diversity with Chl a fluorescence parameters were shown. Values ca. 0.7 are marked in bold.

https://doi.org/10.1371/journal.pone.0156201.t004

They also characterized with increase values of parameter related to the probability of grouping of PSII reaction centers (P_G, Table 3; Fig 5). The second axis, which explained 22.73% variation in the data, separated the group of parameters related to the ChlF rise (V_K) and quantum yield (at t = 0) of energy dissipation (φ_Do), from the group of parameters related to fluorescence intensity F_M, and maximum quantum efficiency of PSII (φ_Po). The third PCA axis, which explained 14.16% of variation, separated the old populations with higher F_O, F_M and RC/CS₀values, from interm age populations (Fig 5).

Redundancy Analysis (RDA, Figs 6–7, and Tables 4 and 5) is a two-table method in which the gradient found in the fluorescence data, observed in PCA, could be directly related to the external variables. The first three RDA axes explained 66.5% variance in the fluorescence data (Table 4). The main gradient of variability along the first RDA axis, which explained 28.9% variation, could be related, as in PCA, to differences among young and other populations (Fig 6). The ChlF parameters recorded for young populations: lower activity of PSII oxygen evolving complex (OEC), lower values of maximum quantum efficiency of PSII (φ_Po) and significantly reduced value RC/CS₀, reflecting the relative density of active reaction centers, higher size of electron acceptor pools available on the reducing side of PSII (S_M), between the both photosystems and at PSI acceptors (Area and N) and higher energy fluxes per one active reaction center are highly positively correlated to the value of the Pareto (beta) index high population genotypic diversity (G) and, weekly—percentage polymorphic (PPL) loci and Chl a/b ratio, (Table 5). This could be explained by populations of Brachypodium colonizing the new habitats, consisting of high numbers of equal size genotypes, characterized, on average, by high efficiency of PSII and high connectivity among PSII units. In the course of expansion some clones started to dominate others. They differed in structural traits expressed by higher Chl_tot content in leaves (Fig 6, Tables 4 and 5).

[Figure omitted. See PDF.]

Fig 6.

Redundancy Analysis (RDA) ordination diagram displaying the relationships between the (i) Chl a fluorescence measurements and parameters and (ii) genetic, genotypic diversity and chlorophyll content in leaves. The first two RDA axes explained of 28.9 and 24.8% of variation in the data respectively. The classification of Chl a fluorescence parameters were performed with k-means clustering. The optimal number of groups (3) was estimated based on Caliński-Harabasz criterion. The result of classification was superimposed onto the ordination diagram. The ellipses represent 95% confidence intervals for the groups.

https://doi.org/10.1371/journal.pone.0156201.g006

[Figure omitted. See PDF.]

Fig 7. The optimal number of groups based on Calinski-Harabasz index.

A. differences in Calinski-Harabasz index. B. The peak in 3b panel point 3 as a optimal number of groups.

https://doi.org/10.1371/journal.pone.0156201.g007

[Figure omitted. See PDF.]

Table 5. Results of Redundancy Analysis (RDA) that relates selected genetic diversity (G—number of AFLP genotypes, Pareto beta index) and leaf chlorophyll content with chlorophyll fluorescence parameters.

Pareto beta is index measuring the spatial distribution of genotypes and could be related to inter- and intraclonal competition. The highest value of beta were recorded in population with many equall-sized genotypes while the lowest with population dominated with few large clones. The results of Monte Carlo permutation test (with 999 permutations) of the variables is presented: Var—variance explained, F—pseudo-F, and significance level

https://doi.org/10.1371/journal.pone.0156201.t005

The second RDA axes, which explained 24.8% variation in the data, differentiated groups of parameters related to the fluorescence rise from those related to maximal efficiency of PSII (φ_Po). On the other hand, the third RDA axis, which explained 12.8%, separated the interm and old populations with similar efficiency of PSII but different F_O, F_M and RC/CS₀ values

Discussion

Understanding of the processes related to establishment and expansion of species in the novel habitats is crucial to nature conservation, for example for development of management plans in order to protect the biodiversity of species-rich grasslands or preventing the colonization and spreading of invasive species. There is ample evidence that the probability of successful establishment in novel environment increases with the number of individuals in founder groups which come from different sources and with the number of repeated introductions [52–55]. However, many theoretical and empirical studies seem to indicate that founder effects and bottlenecks as a scenario of many colonizations, due to geographical separation of populations It have resulted in genetically uniform populations of many expansive species [54,56–59]. As confirmed previously by genetic (AFLP) analyses [15], the colonization of Brachypodium pinnatum follows the first model. It started from seeds from multiple locations and during the first 30–50 years it often produced monodominant stands especially on former fallow lands or on abandoned (i.e. not grazed or mowed) calcareous grasslands. At this stage, the grass formed stands consisted of many genotypes. On the intermediate and old grasslands (from 100 up to >300 years old), a strong decrease in the number of genotypes (G) and slight decrease in the percentage of polymorphic loci was observed (PPL)[6,15]. Although the data on (mostly neutral) genetic variability could be important for distinguishing 'individuals' of Brachypodium in the field conditions, there is the need for additional ecological or physiological data for the assessment of probability maintenance its viability or further expansion [55, 60]After establishment, individuals face a series of environmental stresses (high light, drought, high or low temperature and nutrient deficiency), inter- and intraspecific competition, introduced predators and novel or changed diseases [57]. The performed researches have revealed that the inter–and intraclonal competition and environmental conditions result in strong changes in the ecological properties of Brachypodium individuals and populations: increase in stem density per 1m², clone diameter, lateral vegetative spread, mean number of leaves per ramet and strong decrease in seed production [15]. Moreover, during expansion, an increase in leaf length, leaf width, LA and overall leaf area are also observed. These changes are reflected in the leaf anatomy and showed a plastic acclimatization of Brachypodiumto possibly high-light in young populations and high-light competition in old ones. The leaves of Brachypodium from young populations possesses traits typical for scleromorphic, sun leaves with big chloroplasts, distributed uniformly across mesophyll, while the other ones from old grasslands, additionally to strongly marked scleromorphism, display the optimal distribution of chloroplasts typical for leaves grown in dense stands [61]. The big chloroplasts are distributed close to the lower and upper epidermis.

Apart from changes in population, genotypic diversity and ecological traits during Brachypodium expansions, leaf structural changes and in ChlF parameters related to various aspects of PSII functioning were also visible. Our results confirmed the increase in Chla, Chlb, total Chl content, and leaf greenness index (SPAD)with population age (Table 3). The values of Chl a/b ratio recorded in young populations were close to 3:1, often assumed as optimal to C3 plant [62], while slightly lowered in interm and old populations (2.84–2.88). Our results confirm those obtained by Fabell et all [63], Stroch et all [64] and Živčák et al [41], who observed relatively low differences in Chl a/b ratio in plants growing in high light and shade conditions. The changes in fast Chl a fluorescence kinetics could be related to photosynthesis performance since a light saturation point was found at ~ 1000μmol photons m^–2 s^–2 in young populations and ~ 1200 μmol photons m^–2 s^–2 in old ones (unpublished data—not shown). On the other hand, the increase in total Chl content in old populations is a typical reaction of shade leaves. Pattanayak et al [65] related this phenomenon to an increase in stand density, overall leaf area and possibly due to self-shadowing of leaves within Brachypodium stands. However, Murchie and Horton [66] and Murchie and Lawson [67] had found that in shade-grown plants there was no change in the Chl a/b ratio, while the Chl content decreased. Thus, the decrease in Chl a/b ratio in low light conditions does not seem to be a universal phenomenon, and the level of its dependence on light intensity strongly depends on plant species [41].

The size of the PSII antenna per active reaction center (ABS/RC) was also larger in young populations and these differences are preserved after correction for connectivity [41].The relative density of active PSII reaction centers (RC/CS₀) is significantly lower (of about 40%) in plants from young populations compared to old ones. Hence, both pigment composition and prompt ChlF induction analysis indicate that Brachypodium pinnatum belongs to a group of plants with changeable antenna size [68]. Moreover, in this study, signs of reorganization of PSII units with population age were observed. The O-J part of prompt fluorescence kinetic curve is used to estimate the connectivity parameter among PSII units [23,26,69–71]. The calculated parameters associated with connectivity: probability of connectivity among PSII units–p [36], were 150% and 40% higher in leaves from young populations compared to the interm and old ones. The other connectivity parameters (P_G) were also slightly higher in young populations (Table 3). The value of p parameter ranged between 0.53 to 0.67, which confirmed that the antenna organization of all populations studied followed the “connected unit” model, since the connectivity parameter p obtained ranged between 0 and 1 [69]. This means that the excitation energy of closed RCs can be transferred to a number of nearby open RCs [41].

Biotic stress affects plant growth through reduction of photosynthesis [72]. The initial slope of variable fluorescence M₀, within rapid ChlF kinetics, indicates more rapid initial accumulation of closed RCs in leaves from young populations as compared to the interm and old ones. Moreover, the higher values of ChlF at the J and the I steps, and hence higher V_J and V_I values in the plants from old populations point to limited number of electron carriers on the PSII acceptor side [26,73]. Detailed analysis, based on the selected parameters (S_m, N), suggests a decreased size of the pool of PSII and PSI electron carriers (from Q_A to ferredoxin), as well as, a decrease in the number of Q_A turnovers between F_O and F_M. PSII is one of the most susceptible components of the photosynthetic machinery. Abiotic stresses, such as drought or high light, results in an over-reduction of the electron transport chain (ETC) [74–77]. The specific energy fluxes in one active reaction centers (ABS/RC, TR₀/RC, ET₀/RC, RE₀/RC and DI₀/RC) are visibly decreased with population age. This decrease is especially expressed for parameters correlation with the electron transfer site within PSI–from PQH₂ to PSI end acceptors: RE₀/RC, φ_Ro and δ_Ro. These could be mechanisms by which the PSII of plants from old populations reduce the rate of electron transport chain by converting the excess of absorbed light into thermal energy. [26,78,79]. The diminished efficiency of each RC (PI_ABS and PI_total) is compensated by an increased RC/CS₀ (reflecting active reaction center density) per leaf area as well as by higher total Chl concentration in photosynthesizing tissue of plants from old populations.

The above mentioned results demonstrate the very plastic response of Brachypodium to changeable environmental conditions, reflected in reduction of the maximum quantum yield of PSII (φ_Po) in the early stages of expansion while increasing the quantum yield and probability for electron transport from Q_A- to PQ (φ_Eo; ψ_O). In addition, the γ_Rc decrease in old and interm populations, caused a reduction in the amount of light harvesting complexes in PSII.

Among the ChlF parameters used in this study, the performance index (PI) provides the information on the general state of plants and their vitality [51]. It combines the information about the concentration of the fully active reaction center per chlorophyll, primary photochemistry and electron transport (Strasser et al. 2004). P_ABS is related to energy conservation of the photons absorbed by PSII in the form of reduced intersystem electron acceptors, and PI_total in the form of reduced acceptors of PS I. Changes in PI susceptible to changes in antenna properties, electron trapping efficiency and transport beyond Q_A [51]. Živčák et al [80] pointed out PI as very sensitive index to prolonged drought stress in winter wheat. In this study, significantly higher values of performance indices were found: higher by 50% for PI_ABS, and 300% for PI_total, in young populations.

Conclusions

The results presented above confirm that the strong decrease of genotypes as a result of environmental stress occurs at the early and intermediate stages of Brachypodiumexpansion and is clearly reflected in changes in ChlF parameters. The old stands are dominated by a few genotypes which seem to be the best acclimatized to the self-shading/competition by lowering their photosynthetic performance during light-phase of photosynthesis. On the other hand the 'high-speed' photosynthetic observed in the young populations can be seen as acclimatization to very adverse conditions (probably combination of high light, high temperature, low nutrient water) in which the plants try to catch the sunflecks in the 'windows of opportunity' i.e. early morning or on cloudy days. The high capacity of Brachypodium for plastic morphological and physiological adjustment to changeable habitat light environment refers to its original forest-steppe coenological affinity [8]. Thus, the species has its physiological optimum in half-shade conditions of open forest from which it spreads on fallows or clearings.

The population genetic and ecological analyses are often costly, destructive to vegetation and time-consuming. The presented results clearly confirm that ChlF is a powerful method for inferring physiological mechanisms of expansion of tor grass. The PCA/RDA analyses followed with k-means classification based on the Calinski-Harabasz criterion allowed to distinguish groups of distinct ChlF parameters and allowed them to be related to the changes along stress-competition gradient, occurred during different stages of Brachypodium expansion.

Supporting Information

[Figure omitted. See PDF.]

S1 Fig. Box and whiskers plots of average shoot density of Brachypodium pinnatum/1m².

The shoots were recorded on 20 1m² plots on each Brachypodium populations. The values are averaged over three years: 2013–2015. On the charts the median (line inside the box), box (i.e. inter-quartile range, IQR) and whiskers, defined as 1.5*IQR, are presented. The points are the values beyond the norm (outliers).

https://doi.org/10.1371/journal.pone.0156201.s001

(PDF)

S2 Fig. Box and whiskers plots of average 'clumps of shoots' (see the Materials and Methods) density of Brachypodium pinnatum/1m².

The recordings were performed on 20 1m² plots within each Brachypodium populations. The values are averaged over three years: 2013–2015. On the charts, the median (line inside the box), box (inter-quartile range, IQR) and whiskers (defined as 1.5*IQR) are presented. The points are the values beyond the norm (outliers).

https://doi.org/10.1371/journal.pone.0156201.s002

(PDF)

S3 Fig. Box and whiskers plots of average number of leaves in 'clumps of shoots' of Brachypodium pinnatum/1m².

The shoots were recorded on 20 1m² plots within each Brachypodium populations. The values are averaged over three years: 2013–2015. On the charts, the median (line inside the box), box (inter-quartile range, IQR) and whiskers (defined as 1.5*IQR) are presented. The points are the values beyond the norm (outliers).

https://doi.org/10.1371/journal.pone.0156201.s003

(PDF)

S4 Fig. Box and whiskers plots of average rhizome dry biomass of Brachypodium pinnatum/1m². in 2014, on each Brachypodium population 10 soil samples 20 x 20 x 20 cm were collected randomly within stands, and rhizomes were washed out.

On the charts the median (line inside the box), box (inter-quartile range, IQR) and whiskers (defined as 1.5*IQR) are presented. The points are the values beyond the norm (outliers).

https://doi.org/10.1371/journal.pone.0156201.s004

(PDF)

S5 Fig. Box and whiskers plots of average root dry biomass of Brachypodium pinnatum/1m².

In 2014 on each Brachypodium population 10 soil samples 20 x 20 x 20 cm were collected randomly within stands, and roots were washed out. On the charts, the median (line inside the box), box (inter-quartile range, IQR) and whiskers (defined as 1.5*IQR) are presented. The points are the values beyond the norm (outliers).

https://doi.org/10.1371/journal.pone.0156201.s005

(PDF)

S6 Fig. Box and whiskers plots of average generative shoot density of Brachypodium pinnatum/1m².

The shoots were recorded on each Brachypodium populations on 20 1m² plots. The values are averaged over three years: 2013–2015. On the charts the median (line inside the box), box (inter-quartile range, IQR) and whiskers (defined as 1.5*IQR) are presented. The points are the values beyond the norm (outliers).

https://doi.org/10.1371/journal.pone.0156201.s006

(PDF)

S1 Table. Monthly precipitation [mm] and mean air temperature[°C] in the year of the study.

The values are averaged within age classes.

https://doi.org/10.1371/journal.pone.0156201.s007

(DOCX)

S2 Table. The results of Principal Component Analysis (PCA) and Redundancy Analysis (RDA).

The percentage contributions of Chl of fluorescence parameters to the first three principal components (PC1-PC3) are shown.

https://doi.org/10.1371/journal.pone.0156201.s008

(DOCX)

Acknowledgments

We are very grateful to J. Mitka and I. Franiel for fruitful discussion and two anonymous Reviewers for their valuable comments and suggestions to the manuscript.

Author Contributions

Conceived and designed the experiments: WB HMK. Performed the experiments: WB AKB. Analyzed the data: WB VG. Contributed reagents/materials/analysis tools: WB HMK AKB VG. Wrote the paper: WB HMK AKB VG. Designed the study: WB HMK. Conducted the field study—collecting the materials, performed leaf morphometric, chlorophyll content, SPAD and soil analyses: WB AKB. Performed the statistical analysis; WB VG. Wrote the initial draft of manuscript: WB. Did the editing of the manuscript: WB, AKB, HMK VG. Author contributions in percentage: WB (70%) HMK (10%) AKB (10%) VG (10%).

Citation: Bąba W, Kalaji HM, Kompała-Bąba A, Goltsev V (2016) Acclimatization of Photosynthetic Apparatus of Tor Grass (Brachypodium pinnatum) during Expansion. PLoS ONE 11(6): e0156201. https://doi.org/10.1371/journal.pone.0156201

References

1. Carey MP, Sanderson BL, Barnas KA, Olden JD. Native invaders–challenges for science, management, policy, and society. Front Ecol Environ. 2012;10: 373–381.

2. Bobbink R, Willems JH. Increasing dominance of Brachypodium pinnatum (L.) Beauv. in chalk grasslands: a threat to a species-rich ecosystem. Biol Conserv. 1987;40: 301–314.

3. Bąba W. The species composition and dynamics in well-preserved and restored calcareous xerothermic grasslands (South Poland). Biologia (Bratisl). 2004;59: 447–456.

4. Dujardin G, Bureau F, Decaens T, Langlois E. Morphological and reproductive responses of dominant plant species to local conditions in herbaceous successional stages of a calcareous hillside. Flora. 2011;206: 1030–1039.

5. Bobbink R, Willems JH. Impact of different cutting regimes on the performance of Brachypodium pinnatum in Dutch chalk grassland. Biol Conserv. 1991;56: 1–21.

6. Bąba W, Kurowska M, Kompala-Bąba A, Wilczek A, Długosz J, Szarejko I. Genetic diversity of populations of Brachypodium pinnatum (L.) P. Beauv.: expansive grass in a fragmented landscape. Pol J Ecol. 2012;60: 31–40. Available at: https://www.researchgate.net/publication/236684543_Genetic_diversity_of_populations_of_Brachypodium_pinnatum_L_P_Beauv_Expansive_grass_in_a_fragmented_landscape

7. Bobbink R. Impacts of tropospheric ozone and airborne nitrogenous pollutants on natural and seminatural ecosystems: a commentary. New Phytol. 1998;139: 161–168.

8. Mojzes A, Kalapos T, Viragh K. Plasticity of leaf and shoot morphology and leaf photochemistry for Brachypodium pinnatum (L.) Beauv. growing in contrasting microenvironments in a semiarid loess forest-steppe vegetation mosaic. Flora. 2003;198: 304–320.

9. De Kroon H, Knops J. Habitat exploration through morphological plasticity in two chalk grassland perennials. Oikos. 1990;59: 39–49.

10. Ryser P, Lamberts H. Root and leaf attributes accounting for the performance of fast-growing and slow-growing grasses at different grasses supply. Plant Soil. 1995;170: 251–265.

11. Kahlert B, Ryser P, Edwards P. Leaf phenology of three dominant limestone grassland plants matching the disturbance regime. J Veg Sci. 2005;16: 433–442.

12. Füzy A, Bothe H, Molnár E, Biró B. Mycorrhizal symbiosis effects on growth of chalk false-brome (Brachypodium pinnatum) are dependent on the environmental light regime. J Plant Physiol. 2014;171: 1–6.

13. Idziak D, Hazuka I, Poliwczak B, Wiszyńska A, Wolny E, Hasterok R. Insight into the karyotype evolution of Brachypodium species using comparative chromosome barcoding. PLoS ONE 2014;9: e93503. pmid:24675822

14. Schlapfer F, Fischer M. An isozyme study of clone diversity and relative importance of sexual and vegetative recruitment in the grass Brachypodium pinnatum. Ecography. 1998;21: 351–360.

15. Bąba W, Kurowska M, Kompala-Bąba A, Wilczek A, Dlugosz J, Szarejko I. Genetic diversity of the expansive grass Brachypodium pinnatum in a changing landscape: effect of habitat age. Flora 2012;207: 346–353.

16. Farmer AM, Baxter D. A review of management options for the control of Brachypodium pinnatum in calcareous grassland in England. J Pract Ecol Conserv. 1998;2: 9–18.

17. Hurst A, John E. The effectiveness of glyphosate for controlling Brachypodium pinnatum in chalk grassland. Biol Conserv. 1999;89: 261–265.

18. Krause G, Weis E. Chlorophyll fluorescence and photosynthesis—the basics. Annu Rev Plant Physiol Plant Mol Biol. 1991;42: 313–349.

19. Brestič M, Živčák M, Kalaji HM, Carpentier R, Allakhverdiev SI. Photosystem II thermostability in situ: environmentally induced acclimation and genotype-specific reactions in Triticum aestivum L. Plant Physiol Biochem. 2012;57: 93–105. pmid:22698752

20. Rapacz M, Sasal M, Kalaji HM, Kościelniak J. Is the OJIP test a reliable indicator of winter hardiness and freezing tolerance of common wheat and Triticale under variable winter environments? PLoS ONE. 2015;10: e0134820. pmid:26230839

21. Żurek G, Rybka K, Pogrzeba M, Krzyżak J, Prokopiuk K. Chlorophyll a fluorescence in evaluation of the effect of heavy metal soil contamination on perennial grasses. PLoS ONE. 2014;9.

22. Force L, Critchley C, van Rensen JJS. New fluorescence parameters for monitoring photosynthesis in plants. Photosynth Res. 2003;78: 17–33. pmid:16245061

23. Strasser RJ, Tsimilli-Michael M, Srivastava A. Analysis of the Chlorophyll a Fluorescence Transient. In: Papageorgiou GC, Govindjee , editors. Chlorophyll a Fluorescence. Dordrecht: Springer Netherlands; 2004. pp. 321–362. Available: http://link.springer.com/10.1007/978-1-4020-3218-9_12

24. Schansker G, Toth SZ, Kovacs L, Holzwarth AR, Garab G. Evidence for a fluorescence yield change driven by a light-induced conformational change within photosystem II during the fast chlorophyll a fluorescence rise. Biochim Biophys Acta-Bioenerg. 2011;1807: 1032–1043.

25. Kalaji HM, Carpentier R, Allakhverdiev SI, Bosa K. Fluorescence parameters as early indicators of light stress in barley. J Photochem Photobiol B-Biol. 2012;112: 1–6.

26. Kalaji HM, Schansker G, Ladle RJ, Goltsev V, Bosa K, Allakhverdiev SI, et al. Frequently asked questions about in vivo chlorophyll fluorescence: practical issues. Photosynth Res. 2014;122: 121–158. pmid:25119687

27. Pogan M, Wcisło H. Chromosome numbers of Polish Angiosperms. Acta Biol Cracoviensia Ser Bot. 1991;32: 1–188.

28. Khan MA, Stace CA. Breeding relationship in the genus Brachypodium (Poaceae: Pooidae). Nord J Bot. 1999;19: 257–269.

29. Bąba W. Changes in the structure and floristic composition of the limestone grasslands after cutting trees and shrubs and mowing. Acta Soc Bot Pol. 2003;72: 61–69.

30. Sapek A, Sapek B. Metody analizy chemicznej gleb organicznych. Falenty: Wydawnictwo IMUZ; 1997.

31. Pérez-Harguindeguy N, Díaz S, Garnier E, Lavorel S, Poorter H, Jaureguiberry P, et al. New handbook for standardised measurement of plant functional traits worldwide. Aust J Bot. 2013;61: 167.

32. Strasser RJ, Tsimilli-Michael M. Experimental approach to estimate in vivo the balance of the mechanisms for performance of light reactions, dark reactions, energetic cooperativity, non-Q(A) and non-Q(B)-reducing and non-oxygen-evolving PSII centers. Acta Physiol Plant. 2004;26: 94.

33. Tsimilli-Michael M, Strasser RJ. Photosystem II adaptation to changing light conditions as a process leading to optimality/stability: a model using different types of conformational/state changes. Acta Physiol Plant. 2004;26: 176–177.

34. Stirbet A, Govindjee . On the relation between the Kautsky effect (chlorophyll a fluorescence induction) and Photosystem II: Basics and applications of the OJIP fluorescence transient. J Photochem Photobiol B-Biol. 2011;104: 236–257.

35. Mathis P, editor. Photosynthesis: from light to biosphere: proceedings of the Xth International Photosynthesis Congress, Montpellier, France, 20–25 August 1995. Dordrecht; Boston: Kluwer Academic Publishers; 1995. https://doi.org/10.1007/978-94-009-0173-5

36. Strasser RJ, Stirbet AD. Estimation of the energetic connectivity of PSII centres in plants using the fluorescence rise O-J-I-P—Fitting of experimental data to three different PSII models. Math Comput Simul. 2001;56: 451–461.

37. Oukarroum A, El Madidi S, Schansker G, Strasser RJ. Probing the responses of barley cultivars (Hordeum vulgare L.) by chlorophyll a fluorescence OLKJIP under drought stress and re-watering. Environ Exp Bot. 2007;60: 438–446.

38. Papageorgiou GC, Tsimilli-Michael M, Stamatakis K. The fast and slow kinetics of chlorophyll a fluorescence induction in plants, algae and cyanobacteria: a viewpoint. Photosynth Res. 2007;94: 275–290. pmid:17665151

39. Hiscox JD, Israelstam GF. A method for the extraction of chlorophyll from leaf tissue without maceration. Can J Bot. 1979;57: 1332–1334.

40. Wellburn A. The spectral determination of chlorophyll a and b using various solvents with spectrophotometers with differenr resolutions. J Plant Physiol. 1994;144: 307–313.

41. Živčák M, Brestič M, Kalaji HM, . Photosynthetic responses of sun- and shade-grown barley leaves to high light: is the lower PSII connectivity in shade leaves associated with protection against excess of light? Photosynth Res. 2014;119: 339–354. pmid:24445618

42. Sokal R, Rohlf F. Biometry. The principles and practice of statistics in biological research. New York: W.H.Freeman and Company; 1995.

43. Arnaud-Haond S. Standardizing methods to address clonality in population studies. Mol. Ecol. 17. 2007;5115–5139. doi: https://doi.org/10.1111/j.1365-294X.2007.03535.x

44. Caliński T, Harabasz J. A dendrite method for cluster analysis. Commun Stat—Theory Methods. 1974;3: 1–27.

45. R Core team: A language and environment for statistical computing.R Foundation for Statistical Computing Vienna, Austria [Internet]. 2008. Available: http://www.R-project.org

46. Falińska K, Lembicz M, Jarmołowski A, Borkowska L. Patterns of genetic diversity in populations of Filipendula ulmaria (L.) at different stages of succession on a meadow abandoned for 30 years. Pol J Ecol. 2010;58: 27–40.

47. Kasahara M, Kagawa T, Oikawa K, Suetsugu N, Miyao M, Wada M. Chloroplast avoidance movement reduces photodamage in plants. Nature. 2002;420: 829–832. pmid:12490952

48. Kalaji HM, Oukarroum A, Alexandrov V, Kouzmanova M, Brestič M, Živčák M, et al. Identification of nutrient deficiency in maize and tomato plants by in vivo chlorophyll a fluorescence measurements. Plant Physiol Biochem. 2014;81: 16–25. pmid:24811616

49. Schansker G, Toth SZ, Strasser RJ. Methylviologen and dibromothymoquinone treatments of pea leaves reveal the role of photosystem I in the Chl a fluorescence rise OJIP. Biochim Biophys Acta-Bioenerg. 2005;1706: 250–261.

50. Schansker G, Toth SZ, Strasser RJ. Dark recovery of the Chl a fluorescence transient (OJIP) after light adaptation: The qT-component of non-photochemical quenching is related to an activated photosystem I acceptor side. Biochim Biophys Acta-Bioenerg. 2006;1757: 787–797.

51. Kalaji HM, Jajoo A, Oukarroum A, Brestič M, Živčák M, Samborska IA, et al. The use of chlorophyll fluorescence kinetics analysis to study the performance of photosynthetic machinery in plants. Emerging technologies and management of crop stress tolerance. Elsevier; 2014. pp. 347–384. Available: http://linkinghub.elsevier.com/retrieve/pii/B9780128008751000156

52. Rosenthal DM, Ramakrishnan AP, Cruzan MB. Evidence for multiple sources of invasion and intraspecific hybridization in Brachypodium sylvaticum (Hudson) Beauv. in North America. Mol Ecol. 2008;17: 4657–4669. pmid:18627455

53. Doorduin LJ, Van Den Hof K, Vrieling K, Joshi J. The lack of genetic bottleneck in invasive Tansy Ragwort populations suggests multiple source populations. Basic Appl Ecol. 2010;11: 244–250.

54. Jacquemyn H, Roldan-Ruiz I, Honnay O. Evidence for demographic bottlenecks and limited gene flow leading to low genetic diversity in a rare thistle. Conserv Genet. 2010;11: 1979–1987.

55. Forsman A. Effects of genotypic and phenotypic variation on establishment are important for conservation, invasion, and infection biology. Proc Natl Acad Sci. 2014;111: 302–307. pmid:24367109

56. Charlesworth D, Charlesworth B. Inbreeding Depression and its Evolutionary Consequences. Annu Rev Ecol Syst. 1987;18: 237–268.

57. Frankham R. Stress and adaptation in conservation genetics. J Evol Biol. 2005;18: 750–755. pmid:16033545

58. Crawford KM, Whitney KD. Population genetic diversity influences colonization success. Mol Ecol. 2010;19: 1253–1263. pmid:20456227

59. Sutkowska A, Anamthawat-Jónsson K, Magnússon B, Bąba W, Mitka J. ISSR analysis of two founding plant species on the volcanic island Surtsey, Iceland: grass versus shrub. Surtsey Res. 2015;13: 17–30. Available: http://www.surtsey.is/SRS_publ/2015-XIII/high_res/2%20Surtsey%202015_13%20Sutkowska%2017_30_HRes.pdf

60. Holderegger R, Kamm U, Gugerli F. Adaptive vs. neutral genetic diversity: implications for landscape genetics. Landsc Ecol. 2006;21: 797–807.

61. Way DA, Pearcy RW. Sunflecks in trees and forests: from photosynthetic physiology to global change biology. Tree Physiol. 2012;32: 1066–1081. pmid:22887371

62. Larcher W. Physiological plant ecology: ecophysiology and stress physiology of functional groups. 4th ed. Berlin; New York: Springer; 2003. Available: http://www.springer.com/us/book/9783540435167

63. Falbel TG, Meehl JB, Staehelin LA. Severity of mutant phenotype in a series of chlorophyll-deficient wheat mutants depends on light intensity and the severity of the block in chlorophyll synthesis. Plant Physiol. 1996;112: 821–832. Available: http://www.jstor.org/stable/4277389 pmid:8883392

64. Stroch M, Spunda V, Kurasova I. Non-radiative dissipation of absorbed excitation energy within photosynthetic apparatus of higher plants. Photosynthetica. 2004;42: 323–337.

65. Pattanayak GK, Biswal AK, Reddy VS, Tripathy BC. Light-dependent regulation of chlorophyll b biosynthesis in chlorophyllide a oxygenase overexpressing tobacco plants. Biochem Biophys Res Commun. 2005;326: 466–471. pmid:15582600

66. Murchie E, Horton P. Acclimation of photosynthesis to irradiance and spectral quality in British plant species: Chlorophyll content, photosynthetic capacity and habitat preference. PLANT CELL Environ. 1997;20: 438–448.

67. Murchie EH, Lawson T. Chlorophyll fluorescence analysis: a guide to good practice and understanding some new applications. J Exp Bot. 2013;64: 3983–3998. pmid:23913954

68. Tanaka R, Tanaka A. Chlorophyll b is not just an accessory pigment but a regulator of the photosynthetic antenna. Porphyrins. 2000;9: 240–245.

69. Joliot A, Joliot P. [Kinetic study of the photochemival reaction liberating oxygen during photosynthesis]. Comptes Rendus Hebd Séances Académie Sci. 1964;258: 4622–4625.

70. Joliot P, Joliot A. Excitation transfer between photosynthetic units: the 1964 experiment. In: Govindjee , Beatty JT, Gest H, Allen JF, editors. Discoveries in Photosynthesis. Berlin/Heidelberg: Springer-Verlag; 2005. pp. 187–191. Available: http://link.springer.com/10.1007/1-4020-3324-9_18

71. Stirbet A. Excitonic connectivity between photosystem II units: what is it, and how to measure it? Photosynth Res. 2013;116: 189–214. pmid:23794168

72. Bilgin DD, Zavala JA, Zhu J, Clough SJ, Ort DR, DeLucia EH. Biotic stress globally downregulates photosynthesis genes: Biotic stress downregulates photosynthesis. Plant Cell Environ. 2010;33: 1597–1613. pmid:20444224

73. Lazar D. The polyphasic chlorophyll a fluorescence rise measured under high intensity of exciting light. Funct Plant Biol. 2006;33: 9–30.

74. Foyer CH, Noctor G. Oxidant and antioxidant signalling in plants: a re-evaluation of the concept of oxidative stress in a physiological context. Plant Cell Environ. 2005;28: 1056–1071.

75. Foyer CH, Neukermans J, Queval G, Noctor G, Harbinson J. Photosynthetic control of electron transport and the regulation of gene expression. J Exp Bot. 2012;63: 1637–1661. pmid:22371324

76. Kangasjarvi S, Neukermans J, Li S, Aro E-M, Noctor G. Photosynthesis, photorespiration, and light signalling in defence responses. J Exp Bot. 2012;63: 1619–1636. pmid:22282535

77. Nishiyama Y, Murata N. Revised scheme for the mechanism of photoinhibition and its application to enhance the abiotic stress tolerance of the photosynthetic machinery. Appl Microbiol Biotechnol. 2014;98: 8777–8796. pmid:25139449

78. Tikkanen M, Grieco M, Aro E-M. Novel insights into plant light-harvesting complex II phosphorylation and “state transitions.” Trends Plant Sci. 2011;16: 126–131. pmid:21183394

79. Rochaix J-D. Reprint of: Regulation of photosynthetic electron transport. Biochim Biophys Acta BBA—Bioenerg. 2011;1807: 878–886.

80. Živčák M, Brestič M, Olšovská K, Slamka P. Performance index as a sensitive indicator of water stress in Triticum aestivum L. Plant Soil Environ. 2008;54: 133–139. Available: http://www.agriculturejournals.cz/publicFiles/01163.pdf