Transcripts of wheat at a target locus on chromosome 6B associated with increased yield, leaf mass and chlorophyll index under combined drought and heat stress

Drought and heat stress constrain wheat (Triticum aestivum L.) yields globally. To identify putative mechanisms and candidate genes associated with combined drought and heat stress tolerance, we developed bread wheat near-isogenic lines (NILs) targeting a quantitative trait locus (QTL) on chromosome 6B which was previously associated with combined drought and heat stress tolerance in a diverse panel of wheats. Genotyping-by-sequencing was used to identify additional regions that segregated in allelic pairs between the recurrent and the introduced exotic parent, genome-wide. NILs were phenotyped in a gravimetric platform with precision irrigation and exposed to either drought or to combined drought and heat stress from three days after anthesis. An increase in grain weight in NILs carrying the exotic allele at 6B locus was associated with thicker, greener leaves, higher photosynthetic capacity and increased water use index after re-watering. RNA sequencing of developing grains at early and later stages of treatment revealed 75 genes that were differentially expressed between NILs across both treatments and timepoints. Differentially expressed genes coincided with the targeted QTL on chromosome 6B and regions of genetic segregation on chromosomes 1B and 7A. Pathway enrichment analysis showed the involvement of these genes in cell and gene regulation, metabolism of amino acids and transport of carbohydrates. The majority of these genes have not been characterized previously under drought or heat stress and they might serve as candidate genes for improved abiotic stress tolerance.

Tolerance of Combined Drought and Heat Stress Is Associated With Transpiration Maintenance and Water Soluble Carbohydrates in Wheat Grains

Wheat (Triticum aestivum L.) production is increasingly challenged by simultaneous drought and heatwaves. We assessed the effect of both stresses combined on whole plant water use and carbohydrate partitioning in eight bread wheat genotypes that showed contrasting tolerance. Plant water use was monitored throughout growth, and water-soluble carbohydrates (WSC) and starch were measured following a 3-day heat treatment during drought. Final grain yield was increasingly associated with aboveground biomass and total water use with increasing stress intensity. Combined drought and heat stress immediately reduced daily water use in some genotypes and altered transpiration response to vapor pressure deficit during grain filling, compared to drought only. In grains, glucose and fructose concentrations measured 12 days after anthesis explained 43 and 40% of variation in final grain weight in the main spike, respectively. Starch concentrations in grains offset the reduction in WSC following drought or combined drought and heat stress in some genotypes, while in other genotypes both stresses altered the balance between WSC and starch concentrations. WSC were predominantly allocated to the spike in modern Australian varieties (28-50% of total WSC in the main stem), whereas the stem contained most WSC in older genotypes (67-87%). Drought and combined drought and heat stress increased WSC partitioning to the spike in older genotypes but not in the modern varieties. Ability to maintain transpiration, especially following combined drought and heat stress, appears essential for maintaining wheat productivity.

The Global Durum Wheat Panel (GDP): An International Platform to Identify and Exchange Beneficial Alleles

Representative, broad and diverse collections are a primary resource to dissect genetic diversity and meet pre-breeding and breeding goals through the identification of beneficial alleles for target traits. From 2,500 tetraploid wheat accessions obtained through an international collaborative effort, a Global Durum wheat Panel (GDP) of 1,011 genotypes was assembled that captured 94-97% of the original diversity. The GDP consists of a wide representation of Triticum turgidum ssp. durum modern germplasm and landraces, along with a selection of emmer and primitive tetraploid wheats to maximize diversity. GDP accessions were genotyped using the wheat iSelect 90K SNP array. Among modern durum accessions, breeding programs from Italy, France and Central Asia provided the highest level of genetic diversity, with only a moderate decrease in genetic diversity observed across nearly 50 years of breeding (1970-2018). Further, the breeding programs from Europe had the largest sets of unique alleles. LD was lower in the landraces (0.4 Mbp) than in modern germplasm (1.8 Mbp) at r ² = 0.5. ADMIXTURE analysis of modern germplasm defined a minimum of 13 distinct genetic clusters (k), which could be traced to the breeding program of origin. Chromosome regions putatively subjected to strong selection pressure were identified from fixation index (F _st ) and diversity reduction index (DRI) metrics in pairwise comparisons among decades of release and breeding programs. Clusters of putative selection sweeps (PSW) were identified as co-localized with major loci controlling phenology (Ppd and Vrn), plant height (Rht) and quality (gliadins and glutenins), underlining the role of the corresponding genes as driving elements in modern breeding. Public seed availability and deep genetic characterization of the GDP make this collection a unique and ideal resource to identify and map useful genetic diversity at loci of interest to any breeding program.

Novel Informatic Tools to Support Functional Annotation of the Durum Wheat Genome

Seed mutagenesis is one strategy to create a population with thousands of useful mutations for the direct selection of desirable traits, to introduce diversity into varietal improvement programs, or to generate a mutant collection to support gene functional analysis. However, phenotyping such large collections, where each individual may carry many mutations, is a bottleneck for downstream analysis. Targeting Induced Local Lesions in Genomes (TILLinG), when coupled with next-generation sequencing allows high-throughput mutation discovery and selection by genotyping. We mutagenized an advanced durum breeding line, UAD0951096_F2:5 and performed short-read (2x125 bp) Illumina sequencing of the exome of 100 lines using an available exome capture platform. To improve variant calling, we generated a consolidated exome reference using the recently available genome sequences of the cultivars Svevo and Kronos to facilitate the alignment of reads from the UAD0951096_F2:5 derived mutants. The resulting exome reference was 484.4 Mbp. We also developed a user-friendly, searchable database and bioinformatic analysis pipeline that allowed us to predict zygosity of the mutations discovered and extracts flanking sequences for rapid marker development. Here, we present these tools with the aim of allowing researchers fast and accurate downstream selection of mutations discovered by TILLinG by sequencing to support functional annotation of the durum wheat genome.

Novel Alleles for Combined Drought and Heat Stress Tolerance in Wheat

Drought and heat waves commonly co-occur in many wheat-growing regions causing significant crop losses. The identification of stress associated quantitative trait loci, particularly those for yield, is problematic due to their association with plant phenology and the high genetic × environment interaction. Here we studied a panel of 315 diverse, spring type accessions of bread wheat (Triticum aestivum) in pots in a semi-controlled environment under combined drought and heat stress over 2 years. Importantly, we treated individual plants according to their flowering time. We found 134 out of the 145 identified loci for grain weight that were not associated with either plant phenology or plant height. The majority of loci uncovered here were novel, with favorable alleles widespread in Asian and African landraces providing opportunities for their incorporation into modern varieties through breeding. Using residual heterozygosity in lines from a nested association mapping population, we were able to rapidly develop near-isogenic lines for important target loci. One target locus on chromosome 6A contributed to higher grain weight, harvest index, thousand kernel weight, and grain number under drought and heat stress in field conditions consistent with allelic effects demonstrated in the genome-wide association study.

The physiological and genetic basis of combined drought and heat tolerance in wheat

Drought and heat stress cause losses in wheat productivity in major growing regions worldwide, and both the occurrence and the severity of these events are likely to increase with global climate change. Water deficits and high temperatures frequently occur simultaneously at sensitive growth stages, reducing wheat yields by reducing grain number or weight. Although genetic variation and underlying quantitative trait loci for either individual stress are known, the combination of the two stresses has rarely been studied. Complex and often antagonistic physiology means that genetic loci underlying tolerance to the combined stress are likely to differ from those for drought or heat stress tolerance alone. Here, we review what is known of the physiological traits and genetic control of drought and heat tolerance in wheat and discuss potential physiological traits to study for combined tolerance. We further place this knowledge in the context of breeding for new, more tolerant varieties and discuss opportunities and constraints. We conclude that a fine control of water relations across the growing cycle will be beneficial for combined tolerance and might be achieved through fine management of spatial and temporal gas exchange.

The impact of drought on wheat leaf cuticle properties

BACKGROUND

The plant cuticle is the outermost layer covering aerial tissues and is composed of cutin and waxes. The cuticle plays an important role in protection from environmental stresses and glaucousness, the bluish-white colouration of plant surfaces associated with cuticular waxes, has been suggested as a contributing factor in crop drought tolerance. However, the cuticle structure and composition is complex and it is not clear which aspects are important in determining a role in drought tolerance. Therefore, we analysed residual transpiration rates, cuticle structure and epicuticular wax composition under well-watered conditions and drought in five Australian bread wheat genotypes, Kukri, Excalibur, Drysdale, RAC875 and Gladius, with contrasting glaucousness and drought tolerance.

RESULTS

Significant differences were detected in residual transpiration rates between non-glaucous and drought-sensitive Kukri and four glaucous and drought-tolerant lines. No simple correlation was found between residual transpiration rates and the level of glaucousness among glaucous lines. Modest differences in the thickness of cuticle existed between the examined genotypes, while drought significantly increased thickness in Drysdale and RAC875. Wax composition analyses showed various amounts of C31 β-diketone among genotypes and increases in the content of alkanes under drought in all examined wheat lines.

CONCLUSIONS

The results provide new insights into the relationship between drought stress and the properties and structure of the wheat leaf cuticle. In particular, the data highlight the importance of the cuticle's biochemical makeup, rather than a simple correlation with glaucousness or stomatal density, for water loss under limited water conditions.

Differential expression of microRNAs and potential targets under drought stress in barley

Drought is a crucial environmental constraint limiting crop production in many parts of the world. microRNA (miRNA) based gene regulation has been shown to act in several pathways, including crop response to drought stress. Sequence based profiling and computational analysis have revealed hundreds of miRNAs and their potential targets in different plant species under various stress conditions, but few have been biologically verified. In this study, 11 candidate miRNAs were tested for their expression profiles in barley. Differences in accumulation of only four miRNAs (Ath-miR169b, Osa-miR1432, Hv-miRx5 and Hv-miR166b/c) were observed between drought-treated and well-watered barley in four genotypes. miRNA targets were predicted using degradome analysis of two, different genotypes, and genotype-specific target cleavage was observed. Inverse correlation of mature miRNA accumulation with miRNA target transcripts was also genotype dependent under drought treatment. Drought-responsive miRNAs accumulated predominantly in mesophyll tissues. Our results demonstrate genotype-specific miRNA regulation under drought stress and evidence for their role in mediating expression of target genes for abiotic stress response in barley.

Drought-inducible expression of Hv-miR827 enhances drought tolerance in transgenic barley

Drought is one of the major abiotic stresses reducing crop yield. Since the discovery of plant microRNAs (miRNAs), considerable progress has been made in clarifying their role in plant responses to abiotic stresses, including drought. miR827 was previously reported to confer drought tolerance in transgenic Arabidopsis. We examined barley (Hordeum vulgare L. 'Golden Promise') plants over-expressing miR827 for plant performance under drought. Transgenic plants constitutively expressing CaMV-35S::Ath-miR827 and drought-inducible Zm-Rab17::Hv-miR827 were phenotyped by non-destructive imaging for growth and whole plant water use efficiency (WUE_wp). We observed that the growth, WUE_wp, time to anthesis and grain weight of transgenic barley plants expressing CaMV-35S::Ath-miR827 were negatively affected in both well-watered and drought-treated growing conditions compared with the wild-type plants. In contrast, transgenic plants over-expressing Zm-Rab17::Hv-miR827 showed improved WUE_wp with no growth or reproductive timing change compared with the wild-type plants. The recovery of Zm-Rab17::Hv-miR827 over-expressing plants also improved following severe drought stress. Our results suggest that Hv-miR827 has the potential to improve the performance of barley under drought and that the choice of promoter to control the timing and specificity of miRNA expression is critical.

Identification of reference genes for quantitative expression analysis of microRNAs and mRNAs in barley under various stress conditions

For accurate and reliable gene expression analysis using quantitative real-time reverse transcription PCR (qPCR), the selection of appropriate reference genes as an internal control for normalization is crucial. We hypothesized that non-coding, small nucleolar RNAs (snoRNAs) would be stably expressed in different barley varieties and under different experimental treatments, in different tissues and at different developmental stages of plant growth and therefore might prove to be suitable reference genes for expression analysis of both microRNAs (miRNAs) and mRNAs. In this study, we examined the expression stability of ten candidate reference genes in six barley genotypes under five experimental stresses, drought, fungal infection, boron toxicity, nutrient deficiency and salinity. We compared four commonly used housekeeping genes; Actin (ACT), alpha-Tubulin (α-TUB), Glycolytic glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ADP-ribosylation factor 1-like protein (ADP), four snoRNAs; (U18, U61, snoR14 and snoR23) and two microRNAs (miR168, miR159) as candidate reference genes. We found that ADP, snoR14 and snoR23 were ranked as the best of these candidates across diverse samples. For accurate and reliable gene expression analysis using quantitative real-time reverse transcription PCR (qPCR), the selection of appropriate reference genes as an internal control for normalization is crucial. We hypothesized that non-coding, small nucleolar RNAs (snoRNAs) would be stably expressed in different barley varieties and under different experimental treatments, in different tissues and at different developmental stages of plant growth and therefore might prove to be suitable reference genes for expression analysis of both microRNAs (miRNAs) and mRNAs. In this study, we examined the expression stability of ten candidate reference genes in six barley genotypes under five experimental stresses, drought, fungal infection, boron toxicity, nutrient deficiency and salinity. We compared four commonly used housekeeping genes; Actin (ACT), alpha-Tubulin (α-TUB), Glycolytic glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ADP-ribosylation factor 1-like protein (ADP), four snoRNAs; (U18, U61, snoR14 and snoR23) and two microRNAs (miR168, miR159) as candidate reference genes. We found that ADP, snoR14 and snoR23 were ranked as the best of these candidates across diverse samples. Additionally, we found that miR168 was a suitable reference gene for expression analysis in barley. Finally, we validated the performance of our stable and unstable candidate reference genes for both mRNA and miRNA qPCR data normalization under different stress conditions and demonstrated the superiority of the stable candidates. Our data demonstrate the suitability of barley snoRNAs and miRNAs as potential reference genes for miRNA and mRNA qPCR data normalization under different stress treatments.

Transgenerational inheritance or resetting of stress-induced epigenetic modifications: two sides of the same coin

The transgenerational inheritance of stress-induced epigenetic modifications is still controversial. Despite several examples of defense "priming" and induced genetic rearrangements, the involvement and persistence of transgenerational epigenetic modifications is not known to be general. Here I argue that non-transmission of epigenetic marks through meiosis may be regarded as an epigenetic modification in itself, and that we should understand the implications for plant evolution in the context of both selection for and selection against transgenerational epigenetic memory. Recent data suggest that both epigenetic inheritance and resetting are mechanistically directed and targeted. Stress-induced epigenetic modifications may buffer against DNA sequence-based evolution to maintain plasticity, or may form part of plasticity's adaptive potential. To date we have tended to concentrate on the question of whether and for how long epigenetic memory persists. I argue that we should now re-direct our question to investigate the differences between where it persists and where it does not, to understand the higher order evolutionary methods in play and their contribution.

Pre-conditioning the epigenetic response to high vapor pressure deficit increases the drought tolerance of Arabidopsis thaliana

Epigenetic modification of the genome via cytosine methylation is a dynamic process that responds to changes in the growing environment. This modification can also be heritable. The combination of both properties means that there is the potential for the life experiences of the parental generation to modify the methylation profiles of their offspring and so potentially to “pre-condition” them to better accommodate abiotic conditions encountered by their parents. We recently identified high vapor pressure deficit (vpd)-induced DNA methylation at 2 gene loci in the stomatal development pathway and an associated reduction in leaf stomatal frequency.¹ Here, we test whether this epigenetic modification pre-conditioned parents and their offspring to the more severe water stress of periodic drought. We found that 3 generations of high vpd-grown plants were better able to withstand periodic drought stress over 2 generations. This resistance was not directly associated with de novo methylation of the target stomata genes, but was associated with the cmt3 mutant’s inability to maintain asymmetric sequence context methylation. If our finding applies widely, it could have significant implications for evolutionary biology and breeding for stressful environments.

Transgenerational, dynamic methylation of stomata genes in response to low relative humidity

Transgenerational inheritance of abiotic stress-induced epigenetic modifications in plants has potential adaptive significance and might condition the offspring to improve the response to the same stress, but this is at least partly dependent on the potency, penetrance and persistence of the transmitted epigenetic marks. We examined transgenerational inheritance of low Relative Humidity-induced DNA methylation for two gene loci in the stomatal developmental pathway in Arabidopsis thaliana and the abundance of associated short-interfering RNAs (siRNAs). Heritability of low humidity-induced methylation was more predictable and penetrative at one locus (SPEECHLESS, entropy ≤ 0.02; χ² < 0.001) than the other (FAMA, entropy ≤ 0.17; χ² ns). Methylation at SPEECHLESS correlated positively with the continued presence of local siRNAs (r² = 0.87; p = 0.013) which, however, could be disrupted globally in the progeny under repeated stress. Transgenerational methylation and a parental low humidity-induced stomatal phenotype were heritable, but this was reversed in the progeny under repeated treatment in a previously unsuspected manner.

Low relative humidity triggers RNA-directed de novo DNA methylation and suppression of genes controlling stomatal development

Environmental cues influence the development of stomata on the leaf epidermis, and allow plants to exert plasticity in leaf stomatal abundance in response to the prevailing growing conditions. It is reported that Arabidopsis thaliana 'Landsberg erecta' plants grown under low relative humidity have a reduced stomatal index and that two genes in the stomatal development pathway, SPEECHLESS and FAMA, become de novo cytosine methylated and transcriptionally repressed. These environmentally-induced epigenetic responses were abolished in mutants lacking the capacity for de novo DNA methylation, for the maintenance of CG methylation, and in mutants for the production of short-interfering non-coding RNAs (siRNAs) in the RNA-directed DNA methylation pathway. Induction of methylation was quantitatively related to the induction of local siRNAs under low relative humidity. Our results indicate the involvement of both transcriptional and post-transcriptional gene suppression at these loci in response to environmental stress. Thus, in a physiologically important pathway, a targeted epigenetic response to a specific environmental stress is reported and several of its molecular, mechanistic components are described, providing a tractable platform for future epigenetics experiments. Our findings suggest epigenetic regulation of stomatal development that allows for anatomical and phenotypic plasticity, and may help to explain at least some of the plant's resilience to fluctuating relative humidity.

Epigenomic plasticity within populations: its evolutionary significance and potential

Epigenetics has progressed rapidly from an obscure quirk of heredity into a data-heavy 'omic' science. Our understanding of the molecular mechanisms of epigenomic regulation, and the extent of its importance in nature, are far from complete, but in spite of such drawbacks, population-level studies are extremely valuable: epigenomic regulation is involved in several processes central to evolutionary biology including phenotypic plasticity, evolvability and the mediation of intragenomic conflicts. The first studies of epigenomic variation within populations suggest high levels of phenotypically relevant variation, with the patterns of epigenetic regulation varying between individuals and genome regions as well as with environment. Epigenetic mechanisms appear to function primarily as genome defences, but result in the maintenance of plasticity together with a degree of buffering of developmental programmes; periodic breakdown of epigenetic buffering could potentially cause variation in rates of phenotypic evolution.

Adaptation of tree growth to elevated CO2: quantitative trait loci for biomass in Populus

* Information on the genetic variation of plant response to elevated CO(2) (e[CO(2)]) is needed to understand plant adaptation and to pinpoint likely evolutionary response to future high atmospheric CO(2) concentrations. * Here, quantitative trait loci (QTL) for above- and below-ground tree growth were determined in a pedigree - an F(2) hybrid of poplar (Populus trichocarpa and Populus deltoides), following season-long exposure to either current day ambient CO(2) (a[CO(2)]) or e[CO(2)] at 600 microl l(-1), and genotype by environment interactions investigated. * In the F(2) generation, both above- and below-ground growth showed a significant increase in e[CO(2)]. Three areas of the genome on linkage groups I, IX and XII were identified as important in determining above-ground growth response to e[CO(2)], while an additional three areas of the genome on linkage groups IV, XVI and XIX appeared important in determining root growth response to e[CO(2)]. * These results quantify and identify genetic variation in response to e[CO(2)] and provide an insight into genomic response to the changing environment.

The transcriptome of Populus in elevated CO2

The consequences of increasing atmospheric carbon dioxide for long-term adaptation of forest ecosystems remain uncertain, with virtually no studies undertaken at the genetic level. A global analysis using cDNA microarrays was conducted following 6 yr exposure of Populus x euramericana (clone I-214) to elevated [CO(2)] in a FACE (free-air CO(2) enrichment) experiment. Gene expression was sensitive to elevated [CO(2)] but the response depended on the developmental age of the leaves, and < 50 transcripts differed significantly between different CO(2) environments. For young leaves most differentially expressed genes were upregulated in elevated [CO(2)], while in semimature leaves most were downregulated in elevated [CO(2)]. For transcripts related only to the small subunit of Rubisco, upregulation in LPI 3 and downregulation in LPI 6 leaves in elevated CO(2) was confirmed by anova. Similar patterns of gene expression for young leaves were also confirmed independently across year 3 and year 6 microarray data, and using real-time RT-PCR. This study provides the first clues to the long-term genetic expression changes that may occur during long-term plant response to elevated CO(2).

Stomatal conductance and not stomatal density determines the long-term reduction in leaf transpiration of poplar in elevated CO2

Using a free-air CO2 enrichment (FACE) experiment, poplar trees (Populus x euramericana clone I214) were exposed to either ambient or elevated [CO2] from planting, for a 5-year period during canopy development, closure, coppice and re-growth. In each year, measurements were taken of stomatal density (SD, number mm(-2)) and stomatal index (SI, the proportion of epidermal cells forming stomata). In year 5, measurements were also taken of leaf stomatal conductance (gs, micromol m(-2) s(-1)), photosynthetic CO2 fixation (A, mmol m(-2) s(-1)), instantaneous water-use efficiency (A/E) and the ratio of intercellular to atmospheric CO2 (Ci:Ca). Elevated [CO2] caused reductions in SI in the first year, and in SD in the first 2 years, when the canopy was largely open. In following years, when the canopy had closed, elevated [CO2] had no detectable effects on stomatal numbers or index. In contrast, even after 5 years of exposure to elevated [CO2], gs was reduced, A/E was stimulated, and Ci:Ca was reduced relative to ambient [CO2]. These outcomes from the long-term realistic field conditions of this forest FACE experiment suggest that stomatal numbers (SD and SI) had no role in determining the improved instantaneous leaf-level efficiency of water use under elevated [CO2]. We propose that altered cuticular development during canopy closure may partially explain the changing response of stomata to elevated [CO2], although the mechanism for this remains obscure.

Spatial and temporal effects of free-air CO2 enrichment (POPFACE) on leaf growth, cell expansion, and cell production in a closed canopy of poplar

Leaf expansion in the fast-growing tree, Populus x euramericana was stimulated by elevated [CO(2)] in a closed-canopy forest plantation, exposed using a free air CO(2) enrichment technique enabling long-term experimentation in field conditions. The effects of elevated [CO(2)] over time were characterized and related to the leaf plastochron index (LPI), and showed that leaf expansion was stimulated at very early (LPI, 0-3) and late (LPI, 6-8) stages in development. Early and late effects of elevated [CO(2)] were largely the result of increased cell expansion and increased cell production, respectively. Spatial effects of elevated [CO(2)] were also marked and increased final leaf size resulted from an effect on leaf area, but not leaf length, demonstrating changed leaf shape in response to [CO(2)]. Leaves exhibited a basipetal gradient of leaf development, investigated by defining seven interveinal areas, with growth ceasing first at the leaf tip. Interestingly, and in contrast to other reports, no spatial differences in epidermal cell size were apparent across the lamina, whereas a clear basipetal gradient in cell production rate was found. These data suggest that the rate and timing of cell production was more important in determining leaf shape, given the constant cell size across the leaf lamina. The effect of elevated [CO(2)] imposed on this developmental gradient suggested that leaf cell production continued longer in elevated [CO(2)] and that basal increases in cell production rate were also more important than altered cell expansion for increased final leaf size and altered leaf shape in elevated [CO(2)].