Genomic mechanisms of climate adaptation in polyploid bioenergy switchgrass.

Long-term climate change and periodic environmental extremes threaten food and fuel security1 and global crop productivity2-4. Although molecular and adaptive breeding strategies can buffer the effects of climatic stress and improve crop resilience5, these approaches require sufficient knowledge of the genes that underlie productivity and adaptation6-knowledge that has been limited to a small number of well-studied model systems. Here we present the assembly and annotation of the large and complex genome of the polyploid bioenergy crop switchgrass (Panicum virgatum). Analysis of biomass and survival among 732 resequenced genotypes, which were grown across 10 common gardens that span 1,800 km of latitude, jointly revealed extensive genomic evidence of climate adaptation. Climate-gene-biomass associations were abundant but varied considerably among deeply diverged gene pools. Furthermore, we found that gene flow accelerated climate adaptation during the postglacial colonization of northern habitats through introgression of alleles from a pre-adapted northern gene pool. The polyploid nature of switchgrass also enhanced adaptive potential through the fractionation of gene function, as there was an increased level of heritable genetic diversity on the nondominant subgenome. In addition to investigating patterns of climate adaptation, the genome resources and gene-trait associations developed here provide breeders with the necessary tools to increase switchgrass yield for the sustainable production of bioenergy.

The impact of multifactorial stress combination on reproductive tissues and grain yield of a crop plant.

Global warming, climate change, and industrial pollution are altering our environment subjecting plants, microbiomes, and ecosystems to an increasing number and complexity of abiotic stress conditions, concurrently or sequentially. These conditions, termed, "multifactorial stress combination" (MFSC), can cause a significant decline in plant growth and survival. However, the impacts of MFSC on reproductive tissues and yield of major crop plants are largely unknown. We subjected soybean (Glycine max) plants to a MFSC of up to five different stresses (water deficit, salinity, low phosphate, acidity, and cadmium), in an increasing level of complexity, and conducted integrative transcriptomic-phenotypic analysis of their reproductive and vegetative tissues. We reveal that MFSC has a negative cumulative effect on soybean yield, that each set of MFSC condition elicits a unique transcriptomic response (that is different between flowers and leaves), and that selected genes expressed in leaves or flowers of soybean are linked to the effects of MFSC on different vegetative, physiological, and/or reproductive parameters. Our study identified networks and pathways associated with reactive oxygen species, ascorbic acid and aldarate, and iron/copper signaling/metabolism as promising targets for future biotechnological efforts to augment the resilience of reproductive tissues of major crop plants to MFSC. In addition, we provide unique phenotypic and transcriptomic datasets for dissecting the mechanistic effects of MFSC on the vegetative, physiological, and reproductive processes of a crop plant.

The effects of multifactorial stress combination on rice and maize.

The complexity of environmental factors affecting crops in the field is gradually increasing due to climate change-associated weather events, such as droughts or floods combined with heat waves, coupled with the accumulation of different environmental and agricultural pollutants. The impact of multiple stress conditions on plants was recently termed "multifactorial stress combination" (MFSC) and defined as the occurrence of 3 or more stressors that impact plants simultaneously or sequentially. We recently reported that with the increased number and complexity of different MFSC stressors, the growth and survival of Arabidopsis (Arabidopsis thaliana) seedlings declines, even if the level of each individual stress is low enough to have no significant effect on plants. However, whether MFSC would impact commercial crop cultivars is largely unknown. Here, we reveal that a MFSC of 5 different low-level abiotic stresses (salinity, heat, the herbicide paraquat, phosphorus deficiency, and the heavy metal cadmium), applied in an increasing level of complexity, has a significant negative impact on the growth and biomass of a commercial rice (Oryza sativa) cultivar and a maize (Zea mays) hybrid. Proteomics, element content, and mixOmics analyses of MFSC in rice identified proteins that correlate with the impact of MFSC on rice seedlings, and analysis of 42 different rice genotypes subjected to MFSC revealed substantial genetic variability in responses to this unique state of stress combination. Taken together, our findings reveal that the impacts of MFSC on 2 different crop species are severe and that MFSC may substantially affect agricultural productivity.

A generalist-specialist trade-off between switchgrass cytotypes impacts climate adaptation and geographic range.

Polyploidy results from whole-genome duplication and is a unique form of heritable variation with pronounced evolutionary implications. Different ploidy levels, or cytotypes, can exist within a single species, and such systems provide an opportunity to assess how ploidy variation alters phenotypic novelty, adaptability, and fitness, which can, in turn, drive the development of unique ecological niches that promote the coexistence of multiple cytotypes. Switchgrass, Panicum virgatum, is a widespread, perennial C4 grass in North America with multiple naturally occurring cytotypes, primarily tetraploids (4×) and octoploids (8×). Using a combination of genomic, quantitative genetic, landscape, and niche modeling approaches, we detect divergent levels of genetic admixture, evidence of niche differentiation, and differential environmental sensitivity between switchgrass cytotypes. Taken together, these findings support a generalist (8×)­specialist (4×) trade-off. Our results indicate that the 8× represent a unique combination of genetic variation that has allowed the expansion of switchgrass' ecological niche and thus putatively represents a valuable breeding resource.

The transcriptome of soybean reproductive tissues subjected to water deficit, heat stress, and a combination of water deficit and heat stress.

Global warming and climate change are driving an alarming increase in the frequency and intensity of extreme climate events, such as droughts, heat waves, and their combination, inflicting heavy losses to agricultural production. Recent studies revealed that the transcriptomic responses of different crops to water deficit (WD) or heat stress (HS) are very different from that to a combination of WD + HS. In addition, it was found that the effects of WD, HS, and WD + HS are significantly more devastating when these stresses occur during the reproductive growth phase of crops, compared to vegetative growth. As the molecular responses of different reproductive and vegetative tissues of plants to WD, HS, or WD + HS could be different from each other and these differences could impact many current and future attempts to enhance the resilience of crops to climate change through breeding and/or engineering, we conducted a transcriptomic analysis of different soybean (Glycine max) tissues to WD, HS, and WD + HS. Here we present a reference transcriptomic dataset that includes the response of soybean leaf, pod, anther, stigma, ovary, and sepal to WD, HS, and WD + HS conditions. Mining this dataset for the expression pattern of different stress response transcripts revealed that each tissue had a unique transcriptomic response to each of the different stress conditions. This finding is important as it suggests that enhancing the overall resilience of crops to climate change could require a coordinated approach that simultaneously alters the expression of different groups of transcripts in different tissues in a stress-specific manner.

Differential transpiration between pods and leaves during stress combination in soybean.

Climate change is causing an increase in the frequency and intensity of droughts, heat waves, and their combinations, diminishing agricultural productivity and destabilizing societies worldwide. We recently reported that during a combination of water deficit (WD) and heat stress (HS), stomata on leaves of soybean (Glycine max) plants are closed, while stomata on flowers are open. This unique stomatal response was accompanied by differential transpiration (higher in flowers, while lower in leaves) that cooled flowers during a combination of WD + HS. Here, we reveal that developing pods of soybean plants subjected to a combination of WD + HS use a similar acclimation strategy of differential transpiration to reduce internal pod temperature by approximately 4 °C. We further show that enhanced expression of transcripts involved in abscisic acid degradation accompanies this response and that preventing pod transpiration by sealing stomata causes a significant increase in internal pod temperature. Using an RNA-Seq analysis of pods developing on plants subjected to WD + HS, we also show that the response of pods to WD, HS, or WD + HS is distinct from that of leaves or flowers. Interestingly, we report that although the number of flowers, pods, and seeds per plant decreases under conditions of WD + HS, the seed mass of plants subjected to WD + HS increases compared to plants subjected to HS, and the number of seeds with suppressed/aborted development is lower in WD + HS compared to HS. Taken together, our findings reveal that differential transpiration occurs in pods of soybean plants subjected to WD + HS and that this process limits heat-induced damage to seed production.

Identification of QTLs for symbiotic nitrogen fixation and related traits in a soybean recombinant inbred line population.

KEY MESSAGE: The genetic architecture of symbiotic N fixation and related traits was investigated in the field. QTLs were identified for percent N derived from the atmosphere, shoot [N] and C to N ratio. Soybean [Glycine max (L.) Merr.] is cultivated worldwide and is the most abundant source of plant-based protein. Symbiotic N2 fixation (SNF) in legumes such as soybean is of great importance; however, yields may still be limited by N in both high yielding and stressful environments. To better understand the genetic architecture of SNF and facilitate the development of high yielding cultivars and sustainable soybean production in stressful environments, a recombinant inbred line population consisting of 190 lines, developed from a cross between PI 442012A and PI 404199, was evaluated for N derived from the atmosphere (Ndfa), N concentration ([N]), and C to N ratio (C/N) in three environments. Significant genotype, environment and genotype × environment effects were observed for all three traits. A linkage map was constructed containing 3309 single nucleotide polymorphism (SNP) markers. QTL analysis was performed for additive effects of QTLs, QTL × environment interactions, and QTL × QTL interactions. Ten unique additive QTLs were identified across all traits and environments. Of these, two QTLs were detected for Ndfa and eight for C/N. Of the eight QTLs for C/N, four were also detected for [N]. Using QTL × environment analysis, six QTLs were detected, of which five were also identified in the additive QTL analysis. The QTL × QTL analysis identified four unique epistatic interactions. The results of this study may be used for genomic selection and introgression of favorable alleles for increased SNF, [N], and C/N via marker-assisted selection.

The genetic basis of the root economics spectrum in a perennial grass.

Construction economics of plant roots exhibit predictable relationships with root growth, death, and nutrient uptake strategies. Plant taxa with inexpensively constructed roots tend to more precisely explore nutrient hotspots than do those with costly constructed roots but at the price of more frequent tissue turnover. This trade-off underlies an acquisitive to conservative continuum in resource investment, described as the "root economics spectrum (RES)." Yet the adaptive role and genetic basis of RES remain largely unclear. Different ecotypes of switchgrass (Panicum virgatum) display root features exemplifying the RES, with costly constructed roots in southern lowland and inexpensively constructed roots in northern upland ecotypes. We used an outbred genetic mapping population derived from lowland and upland switchgrass ecotypes to examine the genetic architecture of the RES. We found that absorptive roots (distal first and second orders) were often "deciduous" in winter. The percentage of overwintering absorptive roots was decreased by northern upland alleles compared with southern lowland alleles, suggesting a locally-adapted conservative strategy in warmer and acquisitive strategy in colder regions. Relative turnover of absorptive roots was genetically negatively correlated with their biomass investment per unit root length, suggesting that the key trade-off in framing RES is genetically facilitated. We also detected strong genetic correlations among root morphology, root productivity, and shoot size. Overall, our results reveal the genetic architecture of multiple traits that likely impacts the evolution of RES and plant aboveground-belowground organization. In practice, we provide genetic evidence that increasing switchgrass yield for bioenergy does not directly conflict with enhancing its root-derived carbon sequestration.

Overexpression of LaGRAS enhances phosphorus acquisition via increased root growth of phosphorus-deficient white lupin.

The GRAS transcription factors play an indispensable role in plant growth and responses to environmental stresses. The GRAS gene family has extensively been explored in various plant species; however, the comprehensive investigation of GRAS genes in white lupin remains insufficient. In this study, bioinformatics analysis of white lupin genome revealed 51 LaGRAS genes distributed into 10 distinct phylogenetic clades. Gene structure analyses revealed that LaGRAS proteins were considerably conserved among the same subfamilies. Notably, 25 segmental duplications and a single tandem duplication showed that segmental duplication was the major driving force for the expansion of GRAS genes in white lupin. Moreover, LaGRAS genes exhibited preferential expression in young cluster root and mature cluster roots and may play key roles in nutrient acquisition, particularly phosphorus (P). To validate this, RT-qPCR analysis of white lupin plants grown under +P (normal P) and -P (P deficiency) conditions elucidated significant differences in the transcript level of GRAS genes. Among them, LaGRAS38 and LaGRAS39 were identified as potential candidates with induced expression in MCR under -P. Additionally, white lupin transgenic hairy root overexpressing OE-LaGRAS38 and OE-LaGRAS39 showed increased root growth, and P concentration in root and leaf compared to those with empty vector control, suggesting their role in P acquisition. We believe this comprehensive analysis of GRAS members in white lupin is a first step in exploring their role in the regulation of root growth, tissue development, and ultimately improving P use efficiency in legume crops under natural environments.

Differential regulation of flower transpiration during abiotic stress in annual plants.

Heat waves occurring during droughts can have a devastating impact on yield, especially if they happen during the flowering and seed set stages of the crop cycle. Global warming and climate change are driving an alarming increase in the frequency and intensity of combined drought and heat stress episodes, critically threatening global food security. Because high temperature is detrimental to reproductive processes, essential for plant yield, we measured the inner temperature, transpiration, sepal stomatal aperture, hormone concentrations and transcriptomic response of closed soybean flowers developing on plants subjected to a combination of drought and heat stress. Here, we report that, during a combination of drought and heat stress, soybean plants prioritize transpiration through flowers over transpiration through leaves by opening their flower stomata, while keeping their leaf stomata closed. This acclimation strategy, termed 'differential transpiration', lowers flower inner temperature by about 2-3°C, protecting reproductive processes at the expense of vegetative tissues. Manipulating stomatal regulation, stomatal size and/or stomatal density of flowers could serve as a viable strategy to enhance the yield of different crops and mitigate some of the current and future impacts of global warming and climate change on agriculture.

A multiple ion-uptake phenotyping platform reveals shared mechanisms affecting nutrient uptake by roots.

Nutrient uptake is critical for crop growth and is determined by root foraging in soil. Growth and branching of roots lead to effective root placement to acquire nutrients, but relatively little is known about absorption of nutrients at the root surface from the soil solution. This knowledge gap could be alleviated by understanding sources of genetic variation for short-term nutrient uptake on a root length basis. A modular platform called RhizoFlux was developed for high-throughput phenotyping of multiple ion-uptake rates in maize (Zea mays L.). Using this system, uptake rates were characterized for the crop macronutrients nitrate, ammonium, potassium, phosphate, and sulfate among the Nested Association Mapping (NAM) population founder lines. The data revealed substantial genetic variation for multiple ion-uptake rates in maize. Interestingly, specific nutrient uptake rates (nutrient uptake rate per length of root) were found to be both heritable and distinct from total uptake and plant size. The specific uptake rates of each nutrient were positively correlated with one another and with specific root respiration (root respiration rate per length of root), indicating that uptake is governed by shared mechanisms. We selected maize lines with high and low specific uptake rates and performed an RNA-seq analysis, which identified key regulatory components involved in nutrient uptake. The high-throughput multiple ion-uptake kinetics pipeline will help further our understanding of nutrient uptake, parameterize holistic plant models, and identify breeding targets for crops with more efficient nutrient acquisition.

Physiological trait networks enhance understanding of crop growth and water use in contrasting environments.

Plant function arises from a complex network of structural and physiological traits. Explicit representation of these traits, as well as their connections with other biophysical processes, is required to advance our understanding of plant-soil-climate interactions. We used the Terrestrial Regional Ecosystem Exchange Simulator (TREES) to evaluate physiological trait networks in maize. Net primary productivity (NPP) and grain yield were simulated across five contrasting climate scenarios. Simulations achieving high NPP and grain yield in high precipitation environments featured trait networks conferring high water use strategies: deep roots, high stomatal conductance at low water potential ("risky" stomatal regulation), high xylem hydraulic conductivity and high maximal leaf area index. In contrast, high NPP and grain yield was achieved in dry environments with low late-season precipitation via water conserving trait networks: deep roots, high embolism resistance and low stomatal conductance at low leaf water potential ("conservative" stomatal regulation). We suggest that our approach, which allows for the simultaneous evaluation of physiological traits, soil characteristics and their interactions (i.e., networks), has potential to improve our understanding of crop performance in different environments. In contrast, evaluating single traits in isolation of other coordinated traits does not appear to be an effective strategy for predicting plant performance.

The genetic basis for panicle trait variation in switchgrass (Panicum virgatum).

KEY MESSAGE: We investigate the genetic basis of panicle architecture in switchgrass in two mapping populations across a latitudinal gradient, and find many stable, repeatable genetic effects and limited genetic interactions with the environment. Grass species exhibit large diversity in panicle architecture influenced by genes, the environment, and their interaction. The genetic study of panicle architecture in perennial grasses is limited. In this study, we evaluate the genetic basis of panicle architecture including panicle length, primary branching number, and secondary branching number in an outcrossed switchgrass QTL population grown across ten field sites in the central USA through multi-environment mixed QTL analysis. We also evaluate genetic effects in a diversity panel of switchgrass grown at three of the ten field sites using genome-wide association (GWAS) and multivariate adaptive shrinkage. Furthermore, we search for candidate genes underlying panicle traits in both of these independent mapping populations. Overall, 18 QTL were detected in the QTL mapping population for the three panicle traits, and 146 unlinked genomic regions in the diversity panel affected one or more panicle trait. Twelve of the QTL exhibited consistent effects (i.e., no QTL by environment interactions or no QTL × E), and most (four of six) of the effects with QTL × E exhibited site-specific effects. Most (59.3%) significant partially linked diversity panel SNPs had significant effects in all panicle traits and all field sites and showed pervasive pleiotropy and limited environment interactions. Panicle QTL co-localized with significant SNPs found using GWAS, providing additional power to distinguish between true and false associations in the diversity panel.

QTL × environment interactions underlie adaptive divergence in switchgrass across a large latitudinal gradient.

Local adaptation is the process by which natural selection drives adaptive phenotypic divergence across environmental gradients. Theory suggests that local adaptation results from genetic trade-offs at individual genetic loci, where adaptation to one set of environmental conditions results in a cost to fitness in alternative environments. However, the degree to which there are costs associated with local adaptation is poorly understood because most of these experiments rely on two-site reciprocal transplant experiments. Here, we quantify the benefits and costs of locally adaptive loci across 17° of latitude in a four-grandparent outbred mapping population in outcrossing switchgrass (Panicum virgatum L.), an emerging biofuel crop and dominant tallgrass species. We conducted quantitative trait locus (QTL) mapping across 10 sites, ranging from Texas to South Dakota. This analysis revealed that beneficial biomass (fitness) QTL generally incur minimal costs when transplanted to other field sites distributed over a large climatic gradient over the 2 y of our study. Therefore, locally advantageous alleles could potentially be combined across multiple loci through breeding to create high-yielding regionally adapted cultivars.

The impact of multifactorial stress combination on plant growth and survival.

Climate change-driven extreme weather events, combined with increasing temperatures, harsh soil conditions, low water availability and quality, and the introduction of many man-made pollutants, pose a unique challenge to plants. Although our knowledge of the response of plants to each of these individual conditions is vast, we know very little about how a combination of many of these factors, occurring simultaneously, that is multifactorial stress combination, impacts plants. Seedlings of wild-type and different mutants of Arabidopsis thaliana plants were subjected to a multifactorial stress combination of six different stresses, each applied at a low level, and their survival, physiological and molecular responses determined. Our findings reveal that, while each of the different stresses, applied individually, had a negligible effect on plant growth and survival, the accumulated impact of multifactorial stress combination on plants was detrimental. We further show that the response of plants to multifactorial stress combination is unique and that specific pathways and processes play a critical role in the acclimation of plants to multifactorial stress combination. Taken together our findings reveal that further polluting our environment could result in higher complexities of multifactorial stress combinations that in turn could drive a critical decline in plant growth and survival.

Meta-analysis of drought and heat stress combination impact on crop yield and yield components.

Episodes of prolonged drought coupled with heat waves (i.e. drought and heat combination) can have a devastating impact on agricultural production and crop yield. It is therefore not surprising that improving tolerance to drought and heat combination has been a major goal for breeders and biotech companies. Although much is known about the physiological and molecular responses of vegetative tissues to a combination of drought and heat stress, less is known about the impact of this stress combination on yield and different yield components. Here, we used a meta-analysis approach to synthesize results from over 120 published case studies of crop responses to combined drought and heat stress. Our findings reveal that drought and heat stress combination significantly impacts yield by decreasing harvest index, shortening the life cycle of crops, and altering seed number, size and composition. Furthermore, these impacts are more severe when the stress combination is applied during the reproductive stage of plants. We further identify differences in how legumes and cereals respond to the stress combination and reveal that utilizing C3 or C4 metabolism may not provide an advantage to plants during stress combinations. Taken together our study highlights a need to focus future studies, as well as breeding efforts, on crop responses to drought and heat combination at the reproductive stage of different crop species.

The impact of water deficit and heat stress combination on the molecular response, physiology, and seed production of soybean.

A combination of drought and heat stress, occurring at the vegetative or reproductive growth phase of many different crops can have a devastating impact on yield. In soybean (Glycine max), a considerable effort has been made to develop genotypes with enhanced yield production under conditions of drought or heat stress. However, how these genotypes perform in terms of growth, physiological responses, and most importantly seed production, under conditions of drought and heat combination is mostly unknown. Here, we studied the impact of water deficit and heat stress combination on the physiology, seed production, and yield per plant of two soybean genotypes, Magellan and Plant Introduction (PI) 548313, that differ in their reproductive responses to heat stress. Our findings reveal that although PI 548313 produced more seeds than Magellan under conditions of heat stress, under conditions of water deficit, and heat stress combination its seed production decreased. Because the number of flowers and pollen germination of PI 548313 remained high under heat or water deficit and heat combination, the reduced seed production exhibited by PI 548313 under the stress combination could be a result of processes that occur at the stigma, ovaries and/or other parts of the flower following pollen germination.

Geographic variation in the genetic basis of resistance to leaf rust between locally adapted ecotypes of the biofuel crop switchgrass (Panicum virgatum).

Local adaptation is an important process in plant evolution, which can be impacted by differential pathogen pressures along environmental gradients. However, the degree to which pathogen resistance loci vary in effect across space and time is incompletely described. To understand how the genetic architecture of resistance varies across time and geographic space, we quantified rust (Puccinia spp.) severity in switchgrass (Panicum virgatum) plantings at eight locations across the central USA for 3 yr and conducted quantitative trait locus (QTL) mapping for rust progression. We mapped several variable QTLs, but two large-effect QTLs which we have named Prr1 and Prr2 were consistently associated with rust severity in multiple sites and years, particularly in northern sites. By contrast, there were numerous small-effect QTLs at southern sites, indicating a genotype-by-environment interaction in rust resistance loci. Interestingly, Prr1 and Prr2 had a strong epistatic interaction, which also varied in the strength and direction of effect across space. Our results suggest that abiotic factors covarying with latitude interact with the genetic loci underlying plant resistance to control rust infection severity. Furthermore, our results indicate that segregating genetic variation in epistatically interacting loci may play a key role in determining response to infection across geographic space.

Jasmonic Acid Is Required for Plant Acclimation to a Combination of High Light and Heat Stress.

In the field, plants experience high light (HL) intensities that are often accompanied by elevated temperatures. Such conditions are a serious threat to agriculture production, because photosynthesis is highly sensitive to both HL intensities and high-temperature stress. One of the potential cellular targets of HL and heat stress (HS) combination is PSII because its degree of photoinhibition depends on the balance between the rate of PSII damage (induced by light stress), and the rate of PSII repair (impaired under HS). Here, we studied the responses of Arabidopsis (Arabidopsis thaliana) plants to a combination of HL and HS (HL+HS) conditions. Combined HL+HS was accompanied by irreversible damage to PSII, decreased D1 (PsbA) protein levels, and an enhanced transcriptional response indicative of PSII repair activation. We further identified several unique aspects of this stress combination that included enhanced accumulation of jasmonic acid (JA) and JA-Ile, elevated expression of over 2,200 different transcripts that are unique to the stress combination (including many that are JA-associated), and distinctive structural changes to chloroplasts. A mutant deficient in JA biosynthesis (allene oxide synthase) displayed enhanced sensitivity to combined HL+HS and further analysis revealed that JA is required for regulating several transcriptional responses unique to the stress combination. Our study reveals that JA plays an important role in the acclimation of plants to a combination of HL+HS.

Signal transduction networks during stress combination.

Episodes of heat waves combined with drought can have a devastating impact on agricultural production worldwide. These conditions, as well as many other types of stress combinations, impose unique physiological and developmental demands on plants and require the activation of dedicated pathways. Here, we review recent RNA sequencing studies of stress combination in plants, and conduct a meta-analysis of the transcriptome response of plants to different types of stress combination. Our analysis reveals that each different stress combination is accompanied by its own set of stress combination-specific transcripts, and that the response of different transcription factor families is unique to each stress combination. The alarming rate of increase in global temperatures, coupled with the predicted increase in future episodes of extreme weather, highlight an urgent need to develop crop plants with enhanced tolerance to stress combination. The uniqueness and complexity of the physiological and molecular response of plants to each different stress combination, highlighted here, demonstrate the daunting challenge we face in accomplishing this goal. Dedicated efforts combining field experimentation, omics, and network analyses, coupled with advanced phenotyping and breeding methods, will be needed to address specific crops and particular stress combinations relevant to maintaining our future food chain secured.