Alien chromosome segment from Aegilops speltoides and Dasypyrum villosum increases drought tolerance in wheat via profuse and deep root system
BMC Plant Biology volume 19, Article number: 242 (2019)
Recurrent drought associated with climate change is a major constraint to wheat (Triticum aestivum L.) productivity. This study aimed to (i) quantify the effects of addition/substitution/translocation of chromosome segments from wild relatives of wheat on the root, physiological and yield traits of hexaploid wheat under drought, and (ii) understand the mechanism(s) associated with drought tolerance or susceptibility in wheat-alien chromosome lines.
A set of 48 wheat-alien chromosome lines (addition/substitution/translocation lines) with Chinese Spring background were used. Seedling root traits were studied on solid agar medium. To understand the influence of drought on the root system of adult plants, these 48 lines were grown in 150-cm columns for 65 d under full irrigation or withholding water for 58 d. To quantify the effect of drought on physiological and yield traits, the 48 lines were grown in pots under full irrigation until anthesis; after that, half of the plants were drought stressed by withholding water for 16 d before recording physiological and yield-associated traits.
The alien chromosome lines exhibited altered root architecture and decreased photochemical efficiency and seed yield and its components under drought. The wheat-alien chromosome lines T5DS·5S#3L (TA5088) with a chromosome segment from Aegilops speltoides (5S) and T5DL.5 V#3S (TA5638) with a chromosome segment from Dasypyrum villosum (5 V) were identified as drought tolerant, and the drought tolerance mechanism was associated with a deep, thin and profuse root system.
The two germplasm lines (TA5088 and TA5638) could be used in wheat breeding programs to improve drought tolerance in wheat and understand the underlying molecular genetic mechanisms of root architecture and drought tolerance.
Wheat (Triticum aestivum L.) is a major staple crop, and its production needs to increase by ~ 38% by 2050 to feed the growing population . Among the environmental changes expected in the future, climate models predict an increase in the variability of precipitation leading to an increased frequency and intensity of drought across the globe . Globally, ~ 79% of wheat harvesting regions show grain yield variability due to changes in precipitation and temperature . Therefore, the enhancement of wheat drought tolerance (i.e., maintenance of high yields) is a key challenge in wheat improvement programs .
Bread wheat is an allohexaploid species (AABBDD genomes) that arose ∼8000 years ago , from spontaneous hybridization of the tetraploid wheat T. turgidum L. (AABB genomes) with the diploid goatgrass Aegilops tauschii Coss. (DD genomes) [6, 7]. Wheat has several primary, secondary and tertiary relatives spread across several genera including Triticum, Aegilops (Ae.), Dasypyrum, and Thinopyrum/Agropyron. Wild relatives of wheat are native to semi-arid zones of West and Central Asia and are therefore well adapted to various abiotic stresses [8, 9]. Studies have shown that the introduction of alien chromosome segments from wild relatives into wheat have increased tolerance or resistance to drought , high temperatures , salinity , and water-logging . Similarly, introduction of alien chromosome segments from wild relatives into wheat have improved pest resistance and yield [14, 15]. Waines and Ehdaie  and Yediay et al.  have sucessufully introgressed the genes associated with stem rust (Puccinia graminis) and powdery mildew (Blumeria graminis f. sp. tritici) from rye (Secale cereale L.) to wheat germplasm, which resulted in increased stem rust and powdery mildew tolerance. The wheat lines showing rye-wheat translocation (1RS) had a positive performance on yield, root morphology, and water and nitrogen use efficiency . Sequencing, expression studies, functional annotations, and high-throughput genomics analysis can accelerate allele mining for several traits in sets of chromosome segment substitution lines. Though genetic diversity in wild wheat is a useful resource for trait discovery, only limited numbers of wild relatives have been exploited due to the weedy morphology and low fertility of interspecific hybrids . To utilize the genes and alleles conferring abiotic and biotic stress tolerance from wild relatives of wheat, different genetic materials in the form of addition, substitution, and translocation lines have been developed after laborious efforts in the last few decades . Some translocation lines of wheat-Agropyron elongatum and wheat-rye have been studied for their response to drought stress [4, 10, 19], and lines with rye translocation 1BL-1RS have been used in cultivars across the world to enhance drought tolerance.
Roots appear to be the most relevant organ for breeding drought tolerance, yet limited research is available on root traits due to the difficulty of phenotyping and measuring under both field and controlled environments. Root system architecture refers to the spatial and temporal configuration of roots in the soil. Understanding the variability and contribution of specific root traits can help in the development of drought–tolerant genotypes. In most crops, genetic variation for root traits has been reported using mapping populations or wild relatives [20, 21]. To our knowledge, genetic variation for root architecture in the alien chromosome addition/substitution or translocation lines involving Aegilops species and/or D. villosum under control and/or drought stress is not well understood. Placido et al.  reported that Ag. elongatum 7DL.7EL translocations in wheat improved seed yield under water limiting conditions by increasing root biomass. Similarly, a rye–wheat centric chromosome translocation 1RS.1BL increased yield under drought by enhancing root biomass [10, 19]. The yield advantage of the 1RS translocation line under water-limited conditions was partly associated with deceased root diameter, increased root length density and biomass . Lukaszewski  observed three centric translocations, namely 1RS.1AL, 1RS.1BL, and 1RS.1DL, in Pavon 76 wheat background with greater root biomass and higher grain yield under irrigated and drought conditions. The genetic analyses of 1RS.1BS recombinant breakpoints in Pavon 76 indicated that the distal 15% of the physical length of chromosome 1RS might carry the gene(s) for better rooting ability and root morphological traits .
Relatively large root systems under drought can increase water uptake to alleviate the drought stress effect . In contrast, since roots are a major sink for assimilates, reducing root biomass can increase the availability of assimilates for aboveground parts including grain . Wheat production in India, Australia, and the United States represents a cross-section of global spring wheat production. In India, wheat is grown during winter (in the post-rainy season) and hence, is dependent on water stored in deep soil layers. Similarly, in north eastern Australia and the United States, spring wheat relies largely on stored soil moisture . Hence, targeting access to deeper soil moisture with selected root traits in wheat is critical. Several root morphological traits like root angle, root diameter, and root length density have been associated with increased root system depth and water uptake [25,26,27]. Research into the physiological basis of drought tolerance in wheat is well established; however, the contribution of alien chromosome addition/substitution or translocation lines to drought tolerance is not fully understood. Osipova et al.  mapped quantitative trait locus (QTL) underlying chlorophyll fluorescence parameters and antioxidant enzyme on chromosome 7D of wheat under drought. Bobo et al.  reported that a spring wheat substitution line with the 3D chromosome from winter wheat had reduced quantum yield of photosystem II (ΦPSII) under low light intensities. The substitution lines of durum–Chinese Spring (1B with 1D and 3B with 3D chromosomes) had increased photosynthetic rates compared to the check, indicating that the D genome had a positive interaction with photosynthetic rate . However, Haour-Lurton and Planchon  identified inhibitory effects of specific chromosomes of the D genome on photosynthesis. Wheat–barley addition lines (7H and 7HL) had higher ΦPSII, stomatal conductance and photosynthetic rates under salinity stress . In another study, under terminal high temperature stress, the substitution lines of Chinese Spring (CS)–D. villosum (4 V.3 L, 6 V.3 L, and 5S.3L) had increased chlorophyll index, ΦPSII, individual seed weight and seed yield per plant than the check namely Chinese Spring . Drought during gametogenesis causes a maximum reduction in grain number by inducing floret sterility in wheat and other crops . However, drought during grain filling decreases individual grain weight [33, 34].
A combination of measuring and quantifying root phenotype and whole-plant physiological traits was used in this study to understand the drought-adaptive advantage introduced by the alien chromosome segment into wheat that can benefit wheat breeding programs focused on enhancing drought tolerance. The objectives of this research were to (i) quantify the effects of addition/substitution/translocation of chromosome segments from wheat wild relatives on the root, physiological and yield traits of hexaploid wheat under drought; and (ii) understand the mechanism(s) associated with drought tolerance or susceptibility in wheat-alien chromosome lines.
Principal component analyses
The first two principal components (PC) from the principal component analysis (PCA) represented 54% of the variability, with PC1 explaining 40% of the variance. For PC1, the major contributors to the variance were seed yield per plant (13%) and total root length (12%). For PC2, the major contributors were maximum root length:shoot length ratio (28%) and maximum root length (26%) (Fig. 1a). The 48 wheat-alien chromosome lines and the CS check were grouped into four (I to IV) groups based on the PCA (Fig. 1b). The lines in group I (TA5088 and TA5638) and II (e.g., TA7619 and TA7598) were classified as tolerant and moderately tolerant, respectively, to drought (Fig. 1b). The lines in group III (e.g., TA3583 and TA5584) and IV (CS and TA7659) were classified as susceptible and moderately susceptible, respectively to drought. Among the lines, the highest variation along PC1 was caused by TA5088 (7%) followed by TA5638 (6.5%), and in PC2 the highest variation was caused by TA5088 (3.5%) followed by TA7659 (2.5%) (Fig. 1b). The representative tolerant (TA5088 and TA5638) and susceptible genotypes (TA3583 and TA5584) for different traits were compared with a background check (CS) for root, physiological and yield traits to understand the mechanism of tolerance or susceptibility.
Experiment 1. Genetic variability in seedling root traits
The size of a plant’s root system is a key trait that can affect the uptake of water from the soil. The 48 lines had significant (P > 0.001) genetic variability for the number of seminal roots (2.6–6.3), seedling root length (6.8–23.7 cm) and root angle (38.6–59.3°); their corresponding mean values were 4.2, 18.1 cm and 48.7° (Table 1). Among the accessions, TA5088 and TA5638 had more seminal roots (6.3 vs. 2.6) and a lower root angle (av. 39 vs. 46°) than the CS check (Table 2). Accessions TA3583 and TA5584 had longer seedling roots (av. 22 vs. 10.9 cm) and a greater root angle (59 vs. 46°) than the CS check but no significant difference in the number of seminal roots (Table 2).
Experiment 2. Genetic variability in shoot and root traits under drought
There were significant (P > 0.001) effects of genotype, drought, and their interaction for all growth and root-related traits; except the genotype × drought interactions for plant height and shoot dry weight per plant. Under drought, maximum root length ranged from 58.5–164.5 cm, total root length ranged from 1826 to 10,873 cm, maximum root length:shoot length ratio ranged from 0.95–2.90 and root volume ranged from 0.55–4.60 cm3 in the 48 substitution lines (Table 3). Averaged across the lines, drought significantly (P > 0.001) reduced plant height (7%), maximum root length (26%), maximum root length:shoot length ratio (21%), total root length (7%), root diameter (12%), root volume (23%), number of tillers per plant (33%), and stem dry weight per plant (42%), relative to the control indicating that both shoot and root growth are affected by drought stress (Table 3).
Accessions TA5088 and TA5638 had significantly (P > 0.001) higher maximum root length (av. 161 vs. 124 cm), maximum root length:shoot length ratio (av. 2.8 vs. 2.2) and total root length (av. 10,517 vs. 3768 cm) than the CS check under drought (Table 4). However, the lines TA3583 and TA5584 had lower maximum root length (av. 73 vs. 124 cm), maximum root length:root length ratio (av. 1.2 vs. 2.2) and total root length (av. 2459 vs. 3768 cm) than the CS check under drought (Table 4). Compared to CS, higher maximum root length and total root length was observed in the accessions TA5088 and TA5638 due to its narrow root angle.
Experiment 3. Genetic variability in physiological and yield traits under drought
There were significant (P > 0.001) effects of genotype, drought, and their interaction for chlorophyll index (SPAD units), thylakoid membrane damage (Fo/Fm ratio; unitless), ΦPSII (unitless), electron transport rate (μmol electrons m− 2 s− 1), seed set percentage, seed yield (g spike− 1), seed number plant− 1, individual seed weight (mg seed− 1), seed yield (g plant− 1) and harvest index (%) (Table 5). There was wide genetic variability for various physiological and yield traits under control and drought conditions. Among the 48 lines, thylakoid membrane damage and the ΦPSII ranged from 0.255–0.425, and 0.105–0.465, respectively, under drought. Similarly, seed set percentage and number of seeds per plant ranged between 6.7–59.8% and 2.0–89.0 plant− 1, respectively. Individual seed weights ranged from 3.9–59.9 mg with a mean of 22.8 mg under drought. Seed yield and harvest index ranged between 0.02–2.0 g plant− 1 and from 0.18–29.5%, respectively, under drought. Irrespective of the line, drought significantly (P > 0.001) reduced chlorophyll index (26%), ΦPSII (43%), electron transport rate (44%), seed set percentage (52%), seed yield per spike (63%), number of seeds (60%), individual seed weight (23%), seed yield per plant (66%), and harvest index (54%) but increased thylakoid membrane damage (67%) compared to the control. Overall, the result indicates drought stress had negative effect on both photosynthetic efficiency and yield associated traits. Comparing both photosynthetic efficiency and yield traits the later was found to be more sensitive than former.
The ranking of wheat-alien chromosome lines based on the numerical values of different physiological traits are in Table 6. The data indicate that drought stress increases damage to the thylakoid membrane, and reduces chlorophyll content and ΦPSII in both the deep and shallow rooting accessions, however, the damage was lower in the deep rooting accessions than in the shallow rooting accessions. The lines TA5088 and TA5638 had a higher chlorophyll index and ΦPSII and lower thylakoid membrane damage than the CS check under drought (Table 6), and higher seed set percentages (55.8 and 59.8%, respectively), individual seed weights (44.5 and 45.5 mg seed− 1, respectively) and seed yields (2.27 and 1.61 g plant− 1, respectively) than the CS check (34.0%, 13.2 mg seed− 1 and 0.51 g plant− 1, respectively); both lines were grouped with the six best-performing genotypes under drought (Table 7). The lines TA3583 and TA5584 had lower seed set percentages (15.9 and 15.7%, respectively), individual seed weights (9.4 and 7.8 mg seed− 1, respectively), and seed yields (0.07 and 0.07 g plant− 1, respectively) than the CS check under drought. The higher seed set percentage and individual seed weight was observed in the accessions TA5088 and TA5638 compared to CS, resulting in higher seed yields.
Relationship among root, physiological and yield traits
Total root length was positively associated with chlorophyll index (r2 = 0.21; Fig. 2a), ΦPSII (r2 = 0.39; Fig. 2c), electron transport rate (r2 = 0.60; Fig. 2d), seed set percentage (r2 = 0.52; Fig. 2e), number of seeds per plant (r2 = 0.16; Fig. 2f), individual seed weight (r2 = 0.35; Fig. 2g) and seed yield per plant (r2 = 0.54; Fig. 2h) in the wheat–alien chromosome lines under drought. However, thylakoid membrane damage had a negative association with total root length (r2 = 0.39; Fig. 2b). The lines with more total root length had higher seed set percentage, individual seed weights and seed yields per plant under drought (Fig. 2e, g, h). Although all the above-mentioned traits were associated with total root length, the association between total root length and electron transport rate and seed yield per plant was higher than other traits indicating that these traits might be physiologically related.
The ΦPSII had a positive relationship [co-efficient of determination (r2) ≥ 0.39] with seed set percentage, individual seed weight, seed yield per plant and harvest index (Fig. 3c-f) but a negative relationship with thylakoid membrane damage (r2 = 0.33; Fig. 3b) under drought in the wheat-alien chromosome lines. The lines with a higher ΦPSII had higher seed set percentages, individual seed weights and seed yields per plant (Fig. 3c-e).
This study demonstrated that (i) alien segments in the wheat background have altered root architecture and largely reduced photochemical efficiency and seed yield components under drought, (ii) accessions TA5088 and TA5638 with alien chromosome segments from Ae. speltoides (5S) and D. villosum (5 V) in CS wheat, respectively were drought tolerant, and (iii) the drought tolerance mechanism is associated with a deep, thin and profuse root system.
Root traits play a significant role in plants exposed to drought . Plant capacity to extract soil moisture is a key factor determining drought adaptation, which is likely to result from the improved exploration of soil water . PCA identified accessions TA5088 and TA5638 as drought tolerant (Fig. 1b) and accessions TA3583 and TA5584 as drought susceptible (Fig. 1b) based on 16 root, shoot and yield traits. Accession TA5088 has the long arm of chromosome 5S from Ae. speltoides and TA5638 have the short arm of chromosome 5 V from D. villosum. Ae. speltoides and D. villosum are wild relatives of wheat and rich sources of genetic variation for resistance to drought [37, 38]. The effect of 5S and 5 V translocation in bread wheat was evident from the rooting depth, total root length, and root diameter. In general, cereal roots do not have a vascular cambium and must develop more root length to generate new xylem tissues for water transport . Hence, in water-limiting situations, the genotypes with increased root length had more opportunity for water uptake. This has been confirmed in rice, wheat and maize as deep rooting with profuse branching and small root diameters are associated with drought tolerance [25,26,27, 40]. In contrast, accession TA3583 with a monosomic addition for chromosome 4S from Ae. searsii and TA5584 with the chromosome 7 J arm translocated from Th. Intermedium had wider root angles with sparse root branching at depth and increased root diameter, which is associated with topsoil foraging that eventually results in drought susceptibility . Hence, it is postulated that the gene(s) associated with deep rooting and profuse branching are present in 5S of Ae. speltoides and 5 V of D. villosum and the gene(s) associated with shallow rooting and sparse branching are localized at 4S and 7 J of Ae. searsii and Th. intermedium, respectively.
Studying root architecture and associating it with crop performance under drought can help to identify proxy traits for improving drought tolerance in wheat [36, 40]. The present study clearly showed that total root length and ΦPSII could be used as a proxy trait to evaluate drought tolerance in wheat because of its strong, positive relationship between seed set percentage and seed yield per plant under drought (Fig. 3a). Under drought, the chlorophyll index, ΦPSII and electron transport rate decreased in all 48 lines but the thylakoid membrane damage increased (Table 5). The structural and functional alterations in chloroplast under drought namely damages to thylakoid membranes may lead to chlorophyll loss because chlorophyll is primarily located in the thylakoid membranes . The decreased ΦPSII and electron transport rate under drought shows the occurrence of photoinhibition [42, 43]; and it could be associated with the generation of excess excitation energy . The observed genetic variability for chlorophyll index and ΦPSII in the wheat-alien chromosome lines falls within the range reported by Pour-Aboughadareh et al. . Under drought stress, the ΦPSII in TA5088 and TA5638 increased more (av. 25%) than CS check. Chromosome substitution lines developed from synthetic wheat in the CS background showed that genes regulating Fo, Fm and ΦPSII values might be located on the 3A, 4D and 7A chromosomes . Higher chlorophyll content and lower membrane damage under drought stress in the tolerant genotype of wheat have been reported . Overall, the study suggests that photoinhibition of PSII under drought lead to reduced photochemical efficiency in the wheat-alien chromosome lines and the CS check.
Drought stress during flowering significantly reduced seed set percentage and individual seed weight in all wheat-alien chromosome lines (Table 7), which is attributed to the loss of gametic function and decreased seed filling rate and duration, respectively . Significant genotypic differences were observed for seed set percentage and individual seed weight, which reflects the ability of the alien chromosome segment to endure drought. The wheat-alien chromosome lines TA5088 and TA5638 had higher seed set percentage and individual seed weight than the CS check. This is in accordance with the finding of Fang et al.  that drought-tolerant wheat genotypes had more seeds per plant and higher individual seed weights than drought-susceptible genotypes.
In wheat, the root biomass in subsoil correlates positively with individual seed weight and grain yield under drought stress . Similarly, the ΦPSII correlated positively with grain yield under drought stress . The robust association of these variables with grain yield ratifies that yield is a function of water content and the photosynthetic process. The alien chromosome segments (5S of Ae. speltoides and 5 V of D. villosum) altered the root system and increased PSII photochemistry, which increased grain yield under drought stress. Two wheat-alien chromosome lines (TA5088 and TA5638) will be valuable germplasm for breeding for drought tolerance because of its deep root system and high reproductive success. Breeders can mobilize these translocated segments into adapted local germplasm and estimate the genetic value of these traits in their environment. It is possible that the same alien fragment introgressed into different wheat genotypes could exert different effects on resistance to drought. Therefore, the genetic background and interaction of these genetic factors need to be further investigated. Further, genetic and molecular studies need to be conducted to unravel the genetic factors controlling root- and shoot- related traits in the long arm of chromosome 5S of Ae. speltoides and the short arm of chromosome 5 V of D. villosum, and their contribution to drought tolerance in wheat.
Alien chromosome segments altered root architecture and decreased photochemical efficiency, seed set percentage, individual seed weight, and seed yield per plant in 48 lines of wheat under drought. The wheat-alien chromosome lines TA5088 and TA5638, having chromosome or chromosome segment from Ae. speltoides (5S) and D. villosum (5 V), respectively, were identified as drought tolerant. These two lines had deep, thin and profuse root system under drought stress, which can help to alleviate drought stress by enhancing access to water. Total root length and the ΦPSII were associated with higher seed set percentage and seed yield per plant under drought stress; hence, these traits can be used as proxy traits for improving drought tolerance in wheat. TA5088 and TA5638 lines will be valuable germplasm for understanding the molecular mechanism(s) and breeding of wheat for improved drought tolerance.
Seeds of 48 CS wheat-alien chromosome lines (chromosome addition/substitution /translocation) belonging to Ae. speltoides (4), Ae. searsii (3), Ae. longissima (5), Ae. peregrina (8), Ae. geniculata (14), Th. intermedium (3), L. racemosus (1) and D. villosum (10) and the background CS (check) were obtained from the Wheat Genetic Resources Center, Kansas State University, Manhattan, Kansas. The wheat–alien chromosome lines are derivatives of CS (T. aestivum)–Ae. speltoides, Ae. searsii, Ae. longissima, Ae. peregrina, Ae. geniculata, Th. Intermedium, L. racemosus and D. villosum. Details of the lines used in this study are in Additional file 1: Table S1.
Experiment 1. Genetic variability in seedling root traits
Twenty seeds of each of the 48 lines and the check were surface-sterilized using 10% (v/v) sodium hypochlorite for 5 min and then washed with deionized water for three times. The seeds were germinated in Petri plates using filter paper (Whatman no 42) moistened with five mL of deionized water for two days. Square Petri plates (12 × 12 × 1.7 cm, L × W × H) was used in this study. A slit was made using a scissor on the sides of Petri plates top and bottom and covered with cellophane tape (Staples® Invisible Tape, 2 × 3200 cm, Staples, Manhattan, KS). Sterilized agar (Sigma Type A; 2% w/v) was poured into Petri plates. The Petri plates were sealed with cellophane tape. On the third day, uniformly sized seedlings (radicle emerged) were selected and placed one per Petri plates containing agar in the slit with the radicle facing downward. The Petri plates were incubated at 25 ± 1 °C for 5 d . After the stipulated time, the root angle of individual root axes of seminal roots, counting upwards from the primary seminal root (or radicle), was measured at 3 cm from the seed relative to a vertical line passing through the stem base . Seedling root length was estimated using the cigar roll method .
Experiment 2. Genetic variability in the shoot and root traits under drought
The experiment was conducted in the greenhouse facilities at the Department of Agronomy, Kansas State University, Manhattan, KS to evaluate the variability of root system characteristics among 48 CS wheat–alien chromosome lines and a CS check. Before starting the experiment, the greenhouse was fumigated for one hour using an automated sulfur vaporizer (Rosemania, Franklin, TN) to avoid powdery mildew attack. The rooting medium was Turface MVP® (PROFILE Products LLC, Buffalo Grove, IL), which had a bulk density of 577 ± 32 kg m− 3. Turface is a non-swelling illite and silica clay that allows easy separation of roots. Plants were grown in polyvinyl chloride (PVC) columns with an inside diameter of 7.5 cm and height of 150 cm. The bottom of the PVC columns had plastic caps with a central hole of 0.5 cm diameter for drainage. Before sowing, each PVC column was filled with Turface and fertilized with 4 g of Osmocote (a slow-release fertilizer with 19:6:12 gravimetric percentages of N:P2O5:K2O; Scotts, Marysville, OH, USA) and 1 g of Marathon 1% G (granular; a.i.: Imidacloprid:1-[(6-chloro-3-pyridinyl)methyl]-N-nitro-2-imidazolidinimine; OHP, Inc., Mainland, PA, USA) which were evenly mixed with the Turface in the top 2 cm. Three seeds of a single genotype were sown at 4 cm depth in each PVC column. After emergence, the columns were thinned to one plant per column. Plants were maintained at 24/14 °C (daytime maximum/nighttime minimum temperature) from sowing to harvest (65 d after sowing) at a photoperiod of 16 h (natural light and supplemental fluorescent lights). The fungicide, Bumper 41.8 EC (emulsifiable concentrate; a.i.: Propiconazole:1-[[2-(2,4 dichlorophenyl)-4-propyl-1,3-dioxolan-2-yl]methyl]-1H-1,2,4–triazole; 1.2 mL L− 1; Makhteshim Agan of North America, Inc., Raleigh, NC, USA) was applied on 20 d after sowing to prevent powdery mildew attack. The control plants were maintained at 100% field capacity from sowing to final harvest with drip irrigation. For the drought treatment, plants were stressed by withholding water from day 7 to final harvest (65 days). The duration of drought stress was 58 days at which the genotypes were at the booting stage (Feekes’ stage 10). Our earlier experiments on wheat have indicated that plants grown in Turface have shown the drought stress effects (leaf rolling symptoms) after 45 days of stress imposition . In another experiment withholding water for 41 days in common bean grown in Turface have decreased chlorophyll index by 5% and increased thylakoid membrane damage by 13% . Therefore, in the present study water was withheld for 58 days (first drought symptoms of leaf rolling were observed about 45 d after withholding water) to cause a significant effect on biomass and root growth.
Plant height and number of tillers per plant were measured one day before harvest. Plant height was determined as the distance from the Turface level to the ligule of the youngest leaf and expressed in cm. At harvest, the PVC columns were gently inverted at about 140° to allow the contents (Turface and plants with entire root system) of the column to slip out. The shoots were cut at the base; and the aboveground biomass was oven dried at 65 °C to constant weight, weighed, and expressed as g plant− 1. The roots were carefully separated from the Turface without breaking the root system. The roots were laid on a flat surface and straightened to measure the maximum root length (rooting depth; from the base of the stem to the tip of the root system) and expressed in cm. The root system was carefully washed in water to remove any adhering Turface, placed between the moist paper towels, sealed in Ziploc bags (S.C. Johnson & Sons, Inc. Racine, WI, USA), transported to the laboratory, and stored at 4 °C.
The root system of each plant was sliced into 30-cm-long portions; each portion was submerged in water in a tray (20 × 15 × 2 cm; L × W × H), carefully spread to minimize root overlap, and scanned using an Epson photo scanner (Epson Perfection V700 with 600 dpi resolution, Epson, Long Beach, CA, USA). Images of scanned roots were analyzed using the WinRHIZO Pro image system (Regent Instruments, Inc., Quebec City, QC, Canada) to estimate total root length, root diameter and root volume as explained by McPhee  and Singh et al. . Total root length, root diameter, and root volume were expressed in cm, mm and cm3, respectively. The shoots were oven-dried at 60 °C for 7 d to determine shoot dry weight. Maximum root length:shoot length ratio for each genotype was calculated as the ratio of maximum root length to plant height .
Experiment 3. Genetic variability in physiological and yield traits under drought
Experiment 3 was conducted to evaluate the variability in seed yield and its associated component traits among 48 CS wheat-alien chromosome lines and a CS check under drought conditions. Three seeds of a single genotype were sown in 1.6-L plastic pots [14 cm (height) × 50 cm (top perimeter) × 36 cm (bottom perimeter)] containing a 4:1 mixture of soil [well-drained Kennebec silt loam (a fine-silty, mixed, superactive, mesic Cumulic Hapludoll)] and sand, and 4 g of Osmocote Plus (N:P2O5:K2O = 15:9:12; Scotts, Marysville, OH, USA). Two large indoor growth chambers (Conviron Model CMP 3244, Winnipeg, Manitoba, Canada) were used for this experiment, each being 136 cm wide, 246 cm long, and 180 cm high. Twenty-five lines were randomly placed within each growth chamber (24 CS wheat–alien chromosome lines and a CS check), with four pots per line. The growth chambers were maintained at 24/14 °C (daytime maximum/nighttime minimum temperature) from sowing to physiological maturity at a photoperiod of 16 h and 70% relative humidity (RH), conditions that were optimal for the growth and development of the chromosome substitution lines . The chambers were set at 70% RH to avoid drought due to rapid evapotranspiration. Air temperature and RH were monitored at 20-min intervals in the growth chambers throughout the experiment. Photosynthetically active radiation (PAR) of 680 μmol m− 2 s− 1 was provided by cool white fluorescent lamps (Philips Lighting Co., Somerset, NJ, USA). The PAR was monitored once a month with a Field Scout Light Sensor (Spectrum Technologies, Inc., Plainfield, IL, USA). Fourteen days after seedling emergence, plants were thinned and staked, leaving two plants per pot. Granular Marathon 1% (a.i.: Imidacloprid, 1-((6Chloro-3-pyridinyl) methyl)-N-nitro-2-imidazolidinimine) pesticide was applied to avoid infestation by sucking insect pests. Within a chamber, pots were randomly moved every 7 d to avoid any positional effects. The pots were kept in trays containing ~ 2 cm of water to avoid drought; i.e., the plants were maintained at 100% pot capacity. At the booting stage, the main stem of each plant in the pots was tagged for recording yield and associated traits. Miracle-Gro, a water-soluble fertilizer (N:P2O5:K2O = 24:8:16; Scotts Miracle-Gro Products, Inc., Marysville, OH, USA) was added to the irrigation water (according to manufacturer’s instructions) once in every 7 d until anthesis (Feekes’ 10.5.1 stage). At anthesis, two pots of each line within each growth chamber were randomly assigned to the drought treatment. The drought stress was imposed by withholding water for 16 d. The other two pots were continuously irrigated and served as the control. At the end of the stress, the plants were re-watered and kept well-watered until physiological maturity.
Various physiological traits viz., chlorophyll index, thylakoid membrane damage, ΦPSII, and electron transport rate were measured on tagged plants on 12 d after the drought treatment. Chlorophyll index was measured with a self-calibrating chlorophyll meter (SPAD-502, Spectrum Technologies, Plainfield, IL, USA) on the fully expanded flag leaf of the tagged main stem and expressed in SPAD units. Each time, data were taken in triplicate from the middle portion of the leaf, and the readings were averaged. Chlorophyll a fluorescence parameters were measured using a modulated fluorometer (OS-30p, Opti-Science Inc., Hudson, NH, USA). The minimal fluorescence (Fo) and maximum fluorescence (Fm) was measured in 30-min dark-adapted tagged flag leaves. Thylakoid membrane damage was determined as the ratio of Fo/Fm (unitless). For other fluorescence measurements, the flag leaves were dark adapted for 2 h; the leaves were then continuously irradiated with white actinic light to measure the initial fluorescence in leaves acclimated to irradiation (Fo’), steady-state fluorescence yield (Fs), and maximum fluorescence yield (Fms) of irradiated leaves. Using the above parameters, the ΦPSII (ΦPSII = [Fms – Fs] / Fms; unitless) and apparent rate of photochemical transport of electrons through PSII (ETR = ΦPSII × PAR × 0.5 × 0.84; μmol electrons m− 2 s− 1) were calculated using the instrument software [57, 58], where ETR is electron transport rate, PAR is incident photosynthetically active radiation on a leaf, 0.5 corresponds to the proportion of absorbed quanta used by PSII reaction centers, and 0.84 represents the proportion of incident irradiance absorbed by the leaf .
The spike from the main tiller was tagged at the heading stage. At physiological maturity, the tagged and remaining spikes were harvested separately and dried in an incubator at 40 °C until constant weight. The vegetative biomass per plant was the weight of the oven dried (65 °C for 10 d) plant material without spikes and roots and expressed in gram. Individual spikelets were checked for grain by pressing the floret between the thumb and the index finger. Seed set percentage was estimated as the ratio of spikelets with grain to the total number of spikelets. The tagged spikes were hand threshed after drying, and the number of filled and unfilled grains were estimated for each spike. The remaining spikes were hand threshed, and the grains from these spikes counted and weighed to determine the number of grains per plant, grain weight per spike (g) and grain weight per plant (g). Individual grain weight was calculated by dividing grain weight per plant by number of grains per plant and expressed as mg seed− 1. Harvest index (%) was estimated as the ratio of grain yield to aboveground biomass.
Statistical analyses were performed with SAS 9.4 . Experiment 1 had a completely randomized design with four replications, and Experiments 2 and 3 had a split-plot design in randomized complete block design with two replications. The main plots were water regimes and sub-plots were accessions. Experiments 1 and 3 were repeated. Data from experiment 1 and 3 and their corresponding repeats were statistically analysed independently and found that there were no significant differences. Therefore, the data from both the experiments were pooled together for combined statistical analyses and the mean responses are presented. The PROC GLM procedure of SAS was used for data analysis. Standard errors are shown as an estimate of variability, and the means of various variables are separated for significance by Fisher’s least significant difference (LSD) at 5% significance level. The REG procedure in SAS was used to regress total root length and ΦPSII against other traits. Principal component analysis (PCA) based on the correlation matrix, was performed using XLSTAT-Pro software (AddinSoft, Inc., NY, USA) to identify influential traits under drought stress . PCA biplots were plotted for the drought stress conditions using XLSTAT-Pro software to show the relationships among studied genotypes based on recorded traits.
Electron transport rate
- Fm :
- Fms :
Maximum fluorescence yield
- Fo :
- Fs :
Steady-state fluorescence yield
Least significant difference
Photosynthetically active radiation
Principal component analysis
Quantitative trait locus
Soil Plant Analysis Development
Quantum yield of photosystem II
Ray DK, Mueller ND, West PC, Foley JA. Yield trends are insufficient to double global crop production by 2050. PLoS One. 2013;8(6):66428.
Dai A. Increasing drought under global warming in observations and models. Nat Clim Chang. 2013;3:52–8.
Ray DK, Gerber JS, MacDonald GK, West PC. Climate variation explains a third of global crop yield variability. Nature Comm. 2015;6:5989.
Placido DF, Campbell MT, Folsom JJ, Cui X, Kruger GR, Baezinger PS, et al. Introgression of novel traits from a wild wheat relative improves drought adaptation in wheat. Plant Physiol. 2013;161:1806–19.
Dvorak J, Terlizzi P, Zhang HB, Resta P. The evolution of polyploid wheats: identification of the a genome donor species. Genome. 1993;36:21–31.
Kihara H. Discovery of the DD-analyser, one of the ancestors of Triticum vulgare. Agric Hort. 1944;19:13–4.
Dvorak J, Zhang HB. Variation in repeated nucleotide sequences sheds light on the phylogeny of the wheat B and G genomes. PNAS. 1990;87:9640–4.
Valkoun JJ. Wheat pre-breeding using wild progenitors. Euphytica. 2001;119:17–23.
Liu W, Koo DH, Friebe B, Gill BS. A set of Triticum aestivum-Aegilops speltoides Robertsonian translocation lines. Theor Appl Genet. 2016;129:2359–68.
Ehdaie B, Whitkus RW, Waines JG. Root biomass, water-use efficiency, and performance of wheat-rye translocations of chromosomes 1 and 2 in spring bread wheat ‘Pavon’. Crop Sci. 2003;43:710–7.
Pradhan GP, Prasad PVV. Evaluation of wheat chromosome translocation lines for high temperature stress tolerance at grain filling stage. PLoS One. 2015;10(2):e0116620.
Dvorak J, Edge M, Ross K. On the evolution of the adaptation of Lophopyrum elongatum to growth in saline environments. PNAS. 1988;85:3805–9.
McDonald MP, Galwey NW, Ellneskog-Staam P, Colmer TD. Evaluation of Lophopyrum elongatum as a source of genetic diversity to increase the waterlogging tolerance of hexaploid wheat (Triticum aestivum). New Phytol. 2001;151:369–80.
Monneveux P, Reynolds MP, Gonzalez-Aguilar J, Singh RP, Weber WE. Effects of the 7DL.7Ag translocation from Lophopyrum elongatum on wheat yield and related morphophysiological traits under different environments. Plant Breed. 2003;122:379–84.
Waines JG, Ehdaie B. Domestication and crop physiology: roots of green-revolution wheat. Ann Bot. 2007;100:991–8.
Yediay F, Baloch F, Kilian B, Ozkan H. Testing of rye-specific markers located on 1RS chromosome and distribution of 1AL.RS and 1BL.RS translocations in Turkish wheat (Triticum aestivum L. durum Desf.) varieties and landraces. Genet Resour Crop Evol. 2010;57:119–29.
Yuan B, Cao X, Lv A. Gene introgression from common wheat into Aegilops L. Saudi J Biol Sci. 2017;24:813–6.
Friebe BR, Tuleen NA, Gill BS. Development and identification of a complete set of Triticum aestivum-Aegilops geniculata chromosome addition lines. Genome. 1999;42:374–80.
Sharma S, Bhat P, Ehdaie B, Close T, Lukaszewski A, Waines J. Integrated genetic map and genetic analysis of a region associated with root traits on the short arm of rye chromosome 1 in bread wheat. Theor Appl Genet. 2009;119:783–93.
Liang H, Yongliang Y, Hongqi Y, Lanjie X, Wei D, Hua D, et al. Inheritance and QTL mapping of related root traits in soybean at the seedling stage. Theor Appl Genet. 2014;127:2127–37.
Zhang H, Mittal N, Leamy LJ, Barazani O, Song BH. Back into the wild-apply untapped genetic diversity of wild relatives for crop improvement. Evol Appl. 2017;10:5–24.
Lukaszewski AJ. Reconstruction of complete chromosomes 1B and 1R from the 1RS.1BL translocation of ‘Kavkaz’ origin. Genome. 1993;36:821–4.
Ehdaie B, Layne AP, Waines JG. Root system plasticity to drought influences grain yield in bread wheat. Euphytica. 2012;186:219–32.
Siddique KHM, Belford RK, Tennant D. Root: shoot ratios of old and modern, tall and semi-dwarf wheats in a Mediterranean environment. Plant Soil. 1990;121:89–98.
Wasson AP, Richards RA, Chatrath R, Misra SC, Prasad SV, Rebetzke GK, et al. Traits and selection strategies to improve root systems and water uptake in water-limited wheat crops. J Exp Bot. 2012;63:3485–98.
Comas LH, Becker SR, Curz VMV, Byrne PF, Dierig DA. Root traits contributing to plant productivity under drought. Frontiers in Plant Sci. 2013;4:442.
Uga Y, Sugimoto K, Ogawa S, Rane J, Ishitani M, Hara N, et al. Control of root system architecture by DEEPER ROOTING 1 increases rice yield under drought conditions. Nat Genet. 2013;45:1097–102.
Osipova S, Permyakov A, Permyakova M, Pshenichnikova T, Verkhoturov V, Rudikovsky A, et al. Regions of the bread wheat D genome associated with variation in key photosynthesis traits and shoot biomass under both well-watered and water deficient conditions. J Appl Genet. 2016;57:151–63.
Bobo MS, Planchon C, Morris R. Chromosome 3D influences photosystem II quantum efficiency in winter wheat. Crop Sci. 1992;32:958–61.
Watanabe N, Ogawa A, Kitaya T, Furuta Y. Effect of substituted D genome chromosomes on photosynthesis rate of durum wheat (Triticum aestivum L.). Euphytica. 1994;72:127–31.
Haour-Lurton B, Planchon C. Role of D-genome chromosome in photosynthesis expression in wheats. Theor Appl Genet. 1985;69:443–6.
Darko E, Janda T, Majlath I, Szopko D, Dulai S, Molnar I, et al. Salt stress response of wheat-barley addition lines carrying chromosomes from the winter barley “Manas”. Euphytica. 2015;203:491–504.
P.V.V. Prasad, S.A. Staggenborg, Z. Ristic, Impacts of drought and/or heat stress on physiological, developmental, growth and yield processes of crop plants. In: “Response of Crops to Limited Water: Understanding and Modeling Water Stress Effects on Plant Growth Processes”. (Eds L.H. Ahuja, S.A. Saseendran) (Advances in Agricultural Systems Modeling Series 1, ASA-CSSA: Madison, WI, USA) (2008) pp 301–355.
Saini HS, Aspinall D. Effect of water deficit on sporogenesis in wheat (Triticum aestivum L). Ann Bot. 1981;48:623–33.
Sharma S, Demason DA, Ehdaie B, Lukaszewski AJ, Waines JG. Dosage effect of the short arm of chromosome 1 of rye on root morphology and anatomy in bread wheat. J Exp Bot. 2010;61:2623–33.
Reynolds MP, Quilligan E, Aggarwal PK, Bansal KC, Cavalieri AJ, Chapman SC, et al. An integrated approach to maintaining cereal productivity under climate change. Global Food Security. 2016;8:9–18.
He XL, Wang JW, Li WX, Chen ZZ, Wu J, Zhao JX, et al. An intronless sucrose:fructan-6-fructosyltransferase (6-SFT) gene from Dasypyrum villosum enhances abiotic tolerance in tobacco. Biol Plant. 2017;61:235–45.
Yudina RS, Leonova IN, Salina EA, Khlestkina EK. Effect of alien genomic introgressions on the osmotic stress resistance of wheat. Russian J Genet Appl Res. 2015;5:168–73.
Rich SM, Watt M. Soil conditions and cereal root system architecture: review and considerations for linking Darwin and weaver. J Exp Bot. 2013;64:1193–208.
Manschadi AM, Christopher J, deVoil P, Hammer GL. The role of root architectural traits in adaptation of wheat to water-limited environments. Funct Plant Biol. 2006;33:823–37.
Chen YE, Su YQ, Zhang CM, Ma J, Mao HT, Yang ZH, Yuan M, Zhang ZW, Yuan S, Zhang HY. Comparison of photosynthetic characteristics and antioxidant systems in different wheat strains. J Plant Growth Reg. 2018;37:347–59.
Havaux M, Greppin H, Strasser RJ. Functioning of photosystems I and II in pea leaves exposed to heat stress in the presence or absence of light. Planta. 1991;186:88–98.
Tezara W, Marin O, Rengifo E, Martinez D, Herrera A. Photosynthesis and photoinhibition in two xerophytic shrubs during drought. Photosynthetica. 2005;43:37–45.
Pour-Aboughadareh A, Ahmadi J, Mehrabi AA, Etminan A, Moghaddam M, Siddique KHM. Physiological responses to drought stress in wild relatives of wheat: implications for wheat improvement. Acta Physiol Plant. 2017;39:106. https://doi.org/10.1007/s11738-017-2403-z.
Bai ZY, Li CD, Zhao JF, Wu TY, Zheng JF, Bi CR, et al. Effect and preliminary analysis of chromosomal control on the chlorophyll fluorescence parameters of wheat substitution lines between synthetic hexaploid wheat and Chinese spring under drought stress. Sci Ag Sin. 2011;44:47–57 in Chinese.
Sairam RK, Shukla DS, Saxena DC. Stress induced injury and antioxidant enzymes in relation to drought tolerance in wheat genotypes. Biol Plant. 1997;40:357–64.
Alqudah AM, Samarah NH, Mullen RE. Drought stress effect on crop pollination, seed set, yield and quality. In: Lichtfouse E, editor. Alternative farming systems, biotechnology, drought stress and ecological fertilization. Dordrecht, the Netherlands: Springer; 2011. p. 193–213.
Fang Y, Du Y, Wang J, Wu A, Qiao S, Xu B. Moderate drought stress affected root growth and grain yield in old, modern and newly released cultivars of winter wheat. Front Plant Sci. 2017;8:672.
Paul K, Pauk J, Deak Z, Sass L, Vass I. Contrasting response of biomass and grain yield to severe drought in Cappelle Desprez and plainsman V wheat cultivars. Peer J. 2016;4:e1708.
Manschadi AM, Hammer GL, Christopher JT, de Voil P. Genotypic variation in seedling root architectural traits and implications for drought adaptation in wheat (Triticum aestivum L.). Plant Soil. 2008;303:115–29.
Zhu J, Kaeppler SM, Lynch JP. Mapping of QTLs for lateral root branching and length in maize (Zea mays L.) under differential phosphorus supply. Theor App Genet. 2005;111:688–95.
Djanaguiraman M, Prasad PVV, Kumari J, Rengel Z. Root length and root lipid composition contribute to drought tolerance of winter and spring wheat. Plant Soil. 2018. https://doi.org/10.1007/s11104-018-3794-3.
Sofi PA, Djanaguiraman M, Siddique KHM, Prasad PVV. Reproductive fitness in common bean (Phaseolus vulgaris L.) under drought stress is associated with root length and volume. Ind J Plant Physiol. 2018;23:796–809.
McPhee K. Variation for seedling root architecture in the core collection of pea germplasm. Crop Sci. 2005;45:1758–63.
Singh V, van Oosterom EJ, Jordan DR, Hunt CH, Hammer GL. Genetic variability and control of nodal root angle in sorghum. Crop Sci. 2011;51:2011–20.
Tomar RSS, Tiwari S, Vinod-Naik BK, Chand S, Deshmukh R, et al. Molecular and morpho-agronomical characterization of root architecture at seedling and reproductive stages for drought tolerance in wheat. PLoS One. 2016;11(6):e0156528.
Van Kooten O, Snel JFH. The use of chlorophyll fluorescence nomenclature in plant stress physiology. Photosynth Res. 1990;25:147–50.
Maxwell K, Johnson GN. Chlorophyll fluorescence- a practical guide. J Exp Bot. 2000;51:659–68.
Baker NR. Chlorophyll fluorescence: a probe of photosynthesis in vivo. Ann Rev Plant Biol. 2008;59:89–113.
SAS Institute. The SAS users guide. In: Version 9, vol. 3. USA: SAS Inst, Cary, NC; 2008.
Benson NU, Asuquo FE, Williams AB, Essien JP, Ekong CI, et al. Source evaluation and trace metal contamination in benthic sediments from equatorial ecosystems using multivariate statistical techniques. PLoS One. 2016;11:e0156485.
Friebe B, Qi LL, Liu C, Gill BS. Genetic compensation abilities of Aegilops speltoides chromosomes for homoeologous B-genome chromosomes of polyploid wheat in disomic S(B) chromosome substitution lines. Cytogenet Genome Res. 2011;134:144–50.
Friebe B, Tuleen NA, Gill BS. Standard karyotype of Triticum searsii and its relationship with other S-genome species and common wheat. Theor Appl Genet. 1995;91:248–54.
Friebe B, Tuleen N, Jiang J, Gill BS. Standard karyotype of Triticum longissimum and its cytogenetic relationship with T. aestivum. Genome. 1993;36:731–42.
Friebe B, Tuleen NA, Badaeva ED, Gill BS. Cytogenetic identification of Aegilops peregrinum chromosomes added to common wheat. Genome. 1996b;39:272–6.
Koo H, Liu W, Friebe B, Gill BS. Homoeologous recombination in the presence of Ph1 gene in wheat. Chromosoma. 2017;126:531–40.
Liu C, Qi L, Liu W, Zhao W, Wilson J, et al. Development of a set of compensating Triticum aestivum-Dasypyrum villosum Robertsonian translocation lines. Genome. 2011;54:836–44.
Friebe B, Jiang J, Raupp WJ, McIntosh RA, Gill BS. Characterization of wheat-alien translocations conferring resistance to diseases and pests: current status. Euphytica. 1996a;91:59–87.
Qi LL, Pumphrey MO, Friebe B, Zhang P, Qian C, Bowden RL, Rouse MN, Jin Y, Gill BS. A novel robertsonian translocation event leads to transfer of a stem rust resistance gene (Sr52) effective against race Ug99 from Dasypyrum villosum into bread wheat. Theor Appl Genet. 2011;123:159–67.
Qi LL, Pumphrey MO, Friebe B, Chen PD, Gill BS. Molecular cytogenetic characterization of alien introgressions with gene Fhb3 for resistance to Fusarium head blight disease of wheat. Theor Appl Genet. 2008;117:1155–16.
The senior author (MD) thanks Tamil Nadu Agricultural University, India, for permitting him to perform post-doctoral research at Kansas State University, USA. Mention of trademark or proprietary product does not constitute a guarantee or warranty of the product by Kansas State University and does not imply its approval to the exclusion of other products, which may also be suitable. This publication is a Contribution Number 18-498-J from the Kansas Agricultural Experiment Station.
We thank the Triticeae Coordinated Agricultural Project (Triticeae-CAP), grant no. 2011–68002-30029 and 2017–67007-25939, from the USDA NIFA, United States Agency for International Development (USAID) Feed the Future Innovation Lab for Climate Resilient Wheat (Grant no. AID-0AA-A-13-00008); Feed the Future Innovation Lab for Collaborative Research on Sustainable Intensification (Grant no. AID-0AA-L-14-00006), and Kansas Wheat Commission. The funding from sponsors was used for experimental cost and publication fees. The experimental design, collection and analysis of data, interpretation of data and writing of this manuscript was carried out by the contributing authors.
Availability of data and materials
The data and materials are available from the corresponding authors upon request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Djanaguiraman, M., Prasad, P.V.V., Kumari, J. et al. Alien chromosome segment from Aegilops speltoides and Dasypyrum villosum increases drought tolerance in wheat via profuse and deep root system. BMC Plant Biol 19, 242 (2019). https://doi.org/10.1186/s12870-019-1833-8