Variability in drought gene expression datasets highlight the need for community standardization.

Physiologically relevant drought stress is difficult to apply consistently, and the heterogeneity in experimental design, growth conditions, and sampling schemes make it challenging to compare water deficit studies in plants. Here, we re-analyzed hundreds of drought gene expression experiments across diverse model and crop species and quantified the variability across studies. We found that drought studies are surprisingly uncomparable, even when accounting for differences in genotype, environment, drought severity, and method of drying. Many studies, including most Arabidopsis work, lack high-quality phenotypic and physiological datasets to accompany gene expression, making it impossible to assess the severity or in some cases the occurrence of water deficit stress events. From these datasets, we developed supervised learning classifiers that can accurately predict if RNA-seq samples have experienced a physiologically relevant drought stress, and suggest this can be used as a quality control for future studies. Together, our analyses highlight the need for more community standardization, and the importance of paired physiology data to quantify stress severity for reproducibility and future data analyses.

Topological data analysis reveals a core gene expression backbone that defines form and function across flowering plants.

Since they emerged approximately 125 million years ago, flowering plants have evolved to dominate the terrestrial landscape and survive in the most inhospitable environments on earth. At their core, these adaptations have been shaped by changes in numerous, interconnected pathways and genes that collectively give rise to emergent biological phenomena. Linking gene expression to morphological outcomes remains a grand challenge in biology, and new approaches are needed to begin to address this gap. Here, we implemented topological data analysis (TDA) to summarize the high dimensionality and noisiness of gene expression data using lens functions that delineate plant tissue and stress responses. Using this framework, we created a topological representation of the shape of gene expression across plant evolution, development, and environment for the phylogenetically diverse flowering plants. The TDA-based Mapper graphs form a well-defined gradient of tissues from leaves to seeds, or from healthy to stressed samples, depending on the lens function. This suggests that there are distinct and conserved expression patterns across angiosperms that delineate different tissue types or responses to biotic and abiotic stresses. Genes that correlate with the tissue lens function are enriched in central processes such as photosynthetic, growth and development, housekeeping, or stress responses. Together, our results highlight the power of TDA for analyzing complex biological data and reveal a core expression backbone that defines plant form and function.

Core cellular and tissue-specific mechanisms enable desiccation tolerance in Craterostigma.

Resurrection plants can survive prolonged life without water (anhydrobiosis) in regions with seasonal drying. This desiccation tolerance requires the coordination of numerous cellular processes across space and time, and individual plant tissues face unique constraints related to their function. Here, we analyzed the complex, octoploid genome of the model resurrection plant Craterostigma (C. plantagineum), and surveyed spatial and temporal expression dynamics to identify genetic elements underlying desiccation tolerance. Homeologous genes within the Craterostigma genome have divergent expression profiles, suggesting the subgenomes contribute differently to desiccation tolerance traits. The Craterostigma genome contains almost 200 tandemly duplicated early light-induced proteins, a hallmark trait of desiccation tolerance, with massive upregulation under water deficit. We identified a core network of desiccation-responsive genes across all tissues, but observed almost entirely unique expression dynamics in each tissue during recovery. Roots and leaves have differential responses related to light and photoprotection, autophagy and nutrient transport, reflecting their divergent functions. Our findings highlight a universal set of likely ancestral desiccation tolerance mechanisms to protect cellular macromolecules under anhydrobiosis, with secondary adaptations related to tissue function.

Cross-species predictive modeling reveals conserved drought responses between maize and sorghum.

Drought tolerance is a highly complex trait controlled by numerous interconnected pathways with substantial variation within and across plant species. This complexity makes it difficult to distill individual genetic loci underlying tolerance, and to identify core or conserved drought-responsive pathways. Here, we collected drought physiology and gene expression datasets across diverse genotypes of the C4 cereals sorghum and maize and searched for signatures defining water-deficit responses. Differential gene expression identified few overlapping drought-associated genes across sorghum genotypes, but using a predictive modeling approach, we found a shared core drought response across development, genotype, and stress severity. Our model had similar robustness when applied to datasets in maize, reflecting a conserved drought response between sorghum and maize. The top predictors are enriched in functions associated with various abiotic stress-responsive pathways as well as core cellular functions. These conserved drought response genes were less likely to contain deleterious mutations than other gene sets, suggesting that core drought-responsive genes are under evolutionary and functional constraints. Our findings support a broad evolutionary conservation of drought responses in C4 grasses regardless of innate stress tolerance, which could have important implications for developing climate resilient cereals.

A comparative genomics examination of desiccation tolerance and sensitivity in two sister grass species.

Desiccation tolerance is an ancient and complex trait that spans all major lineages of life on earth. Although important in the evolution of land plants, the mechanisms that underlay this complex trait are poorly understood, especially for vegetative desiccation tolerance (VDT). The lack of suitable closely related plant models that offer a direct contrast between desiccation tolerance and sensitivity has hampered progress. We have assembled high-quality genomes for two closely related grasses, the desiccation-tolerant Sporobolus stapfianus and the desiccation-sensitive Sporobolus pyramidalis Both species are complex polyploids; S. stapfianus is primarily tetraploid, and S. pyramidalis is primarily hexaploid. S. pyramidalis undergoes a major transcriptome remodeling event during initial exposure to dehydration, while S. stapfianus has a muted early response, with peak remodeling during the transition between 1.5 and 1.0 grams of water (gH2O) g-1 dry weight (dw). Functionally, the dehydration transcriptome of S. stapfianus is unrelated to that for S. pyramidalis A comparative analysis of the transcriptomes of the hydrated controls for each species indicated that S. stapfianus is transcriptionally primed for desiccation. Cross-species comparative analyses indicated that VDT likely evolved from reprogramming of desiccation tolerance mechanisms that evolved in seeds and that the tolerance mechanism of S. stapfianus represents a recent evolution for VDT within the Chloridoideae. Orthogroup analyses of the significantly differentially abundant transcripts reconfirmed our present understanding of the response to dehydration, including the lack of an induction of senescence in resurrection angiosperms. The data also suggest that failure to maintain protein structure during dehydration is likely critical in rendering a plant desiccation sensitive.

Evolutionary innovations driving abiotic stress tolerance in C4 grasses and cereals.

Grasslands dominate the terrestrial landscape, and grasses have evolved complex and elegant strategies to overcome abiotic stresses. The C4 grasses are particularly stress tolerant and thrive in tropical and dry temperate ecosystems. Growing evidence suggests that the presence of C4 photosynthesis alone is insufficient to account for drought resilience in grasses, pointing to other adaptations as contributing to tolerance traits. The majority of grasses from the Chloridoideae subfamily are tolerant to drought, salt, and desiccation, making this subfamily a hub of resilience. Here, we discuss the evolutionary innovations that make C4 grasses so resilient, with a particular emphasis on grasses from the Chloridoideae (chloridoid) and Panicoideae (panicoid) subfamilies. We propose that a baseline level of resilience in chloridoid ancestors allowed them to colonize harsh habitats, and these environments drove selective pressure that enabled the repeated evolution of abiotic stress tolerance traits. Furthermore, we suggest that a lack of evolutionary access to stressful environments is partially responsible for the relatively poor stress resilience of major C4 crops compared to their wild relatives. We propose that chloridoid crops and the subfamily more broadly represent an untapped reservoir for improving resilience to drought and other abiotic stresses in cereals.

Intertwined signatures of desiccation and drought tolerance in grasses.

Grasses are among the most resilient plants, and some can survive prolonged desiccation in semiarid regions with seasonal rainfall. However, the genetic elements that distinguish grasses that are sensitive versus tolerant to extreme drying are largely unknown. Here, we leveraged comparative genomic approaches with the desiccation-tolerant grass Eragrostis nindensis and the related desiccation-sensitive cereal Eragrostis tef to identify changes underlying desiccation tolerance. These analyses were extended across C4 grasses and cereals to identify broader evolutionary conservation and divergence. Across diverse genomic datasets, we identified changes in chromatin architecture, methylation, gene duplications, and expression dynamics related to desiccation in E. nindensis It was previously hypothesized that transcriptional rewiring of seed desiccation pathways confers vegetative desiccation tolerance. Here, we demonstrate that the majority of seed-dehydration-related genes showed similar expression patterns in leaves of both desiccation-tolerant and -sensitive species. However, we identified a small set of seed-related orthologs with expression specific to desiccation-tolerant species. This supports a broad role for seed-related genes, where many are involved in typical drought responses, with only a small subset of crucial genes specifically induced in desiccation-tolerant plants.

Exceptional subgenome stability and functional divergence in the allotetraploid Ethiopian cereal teff.

Teff (Eragrostis tef) is a cornerstone of food security in the Horn of Africa, where it is prized for stress resilience, grain nutrition, and market value. Here, we report a chromosome-scale assembly of allotetraploid teff (variety Dabbi) and patterns of subgenome dynamics. The teff genome contains two complete sets of homoeologous chromosomes, with most genes maintaining as syntenic gene pairs. TE analysis allows us to estimate that the teff polyploidy event occurred ~1.1 million years ago (mya) and that the two subgenomes diverged ~5.0 mya. Despite this divergence, we detect no large-scale structural rearrangements, homoeologous exchanges, or biased gene loss, in contrast to many other allopolyploids. The two teff subgenomes have partitioned their ancestral functions based on divergent expression across a diverse expression atlas. Together, these genomic resources will be useful for accelerating breeding of this underutilized grain crop and for fundamental insights into polyploid genome evolution.

Transcriptome-Based Prediction of Complex Traits in Maize.

The ability to predict traits from genome-wide sequence information (i.e., genomic prediction) has improved our understanding of the genetic basis of complex traits and transformed breeding practices. Transcriptome data may also be useful for genomic prediction. However, it remains unclear how well transcript levels can predict traits, particularly when traits are scored at different development stages. Using maize (Zea mays) genetic markers and transcript levels from seedlings to predict mature plant traits, we found that transcript and genetic marker models have similar performance. When the transcripts and genetic markers with the greatest weights (i.e., the most important) in those models were used in one joint model, performance increased. Furthermore, genetic markers important for predictions were not close to or identified as regulatory variants for important transcripts. These findings demonstrate that transcript levels are useful for predicting traits and that their predictive power is not simply due to genetic variation in the transcribed genomic regions. Finally, genetic marker models identified only 1 of 14 benchmark flowering-time genes, while transcript models identified 5. These data highlight that, in addition to being useful for genomic prediction, transcriptome data can provide a link between traits and variation that cannot be readily captured at the sequence level.

Massive Tandem Proliferation of ELIPs Supports Convergent Evolution of Desiccation Tolerance across Land Plants.

Desiccation tolerance was a critical adaptation for the colonization of land by early nonvascular plants. Resurrection plants have maintained or rewired these ancestral protective mechanisms, and desiccation-tolerant species are dispersed across the land plant phylogeny. Although common physiological, biochemical, and molecular signatures are observed across resurrection plant lineages, features underlying the recurrent evolution of desiccation tolerance are unknown. Here we used a comparative approach to identify patterns of genome evolution and gene duplication associated with desiccation tolerance. We identified a single gene family with dramatic expansion in all sequenced resurrection plant genomes and no expansion in desiccation-sensitive species. This gene family of early light-induced proteins (ELIPs) expanded in resurrection plants convergent through repeated tandem gene duplication. ELIPs are universally highly expressed during desiccation in all surveyed resurrection plants and may play a role in protecting against photooxidative damage of the photosynthetic apparatus during prolonged dehydration. Photosynthesis is particularly sensitive to dehydration, and the increased abundance of ELIPs may help facilitate the rapid recovery observed for most resurrection plants. Together, these observations support convergent evolution of desiccation tolerance in land plants through tandem gene duplication.

Desiccation Tolerance Evolved through Gene Duplication and Network Rewiring in Lindernia.

Although several resurrection plant genomes have been sequenced, the lack of suitable dehydration-sensitive outgroups has limited genomic insights into the origin of desiccation tolerance. Here, we utilized a comparative system of closely related desiccation-tolerant (Lindernia brevidens) and -sensitive (Lindernia subracemosa) species to identify gene- and pathway-level changes associated with the evolution of desiccation tolerance. The two high-quality Lindernia genomes we assembled are largely collinear, and over 90% of genes are conserved. L. brevidens and L. subracemosa have evidence of an ancient, shared whole-genome duplication event, and retained genes have neofunctionalized, with desiccation-specific expression in L. brevidens Tandem gene duplicates also are enriched in desiccation-associated functions, including a dramatic expansion of early light-induced proteins from 4 to 26 copies in L. brevidens A comparative differential gene coexpression analysis between L. brevidens and L. subracemosa supports extensive network rewiring across early dehydration, desiccation, and rehydration time courses. Many LATE EMBRYOGENESIS ABUNDANT genes show significantly higher expression in L. brevidens compared with their orthologs in L. subracemosa Coexpression modules uniquely upregulated during desiccation in L. brevidens are enriched with seed-specific and abscisic acid-associated cis-regulatory elements. These modules contain a wide array of seed-associated genes that have no expression in the desiccation-sensitive L. subracemosa Together, these findings suggest that desiccation tolerance evolved through a combination of gene duplications and network-level rewiring of existing seed desiccation pathways.

Extreme haplotype variation in the desiccation-tolerant clubmoss Selaginella lepidophylla.

Plant genome size varies by four orders of magnitude, and most of this variation stems from dynamic changes in repetitive DNA content. Here we report the small 109 Mb genome of Selaginella lepidophylla, a clubmoss with extreme desiccation tolerance. Single-molecule sequencing enables accurate haplotype assembly of a single heterozygous S. lepidophylla plant, revealing extensive structural variation. We observe numerous haplotype-specific deletions consisting of largely repetitive and heavily methylated sequences, with enrichment in young Gypsy LTR retrotransposons. Such elements are active but rapidly deleted, suggesting "bloat and purge" to maintain a small genome size. Unlike all other land plant lineages, Selaginella has no evidence of a whole-genome duplication event in its evolutionary history, but instead shows unique tandem gene duplication patterns reflecting adaptation to extreme drying. Gene expression changes during desiccation in S. lepidophylla mirror patterns observed across angiosperm resurrection plants.

A chromosome-scale assembly of the model desiccation tolerant grass Oropetium thomaeum.

Oropetium thomaeum is an emerging model for desiccation tolerance and genome size evolution in grasses. A draft genome of Oropetium was recently sequenced, but the lack of a chromosome-scale assembly has hindered comparative analyses and downstream functional genomics. Here, we reassembled Oropetium, and anchored the genome into 10 chromosomes using high-throughput chromatin conformation capture (Hi-C) based chromatin interactions. A combination of high-resolution RNAseq data and homology-based gene prediction identified thousands of new, conserved gene models that were absent from the V1 assembly. This includes thousands of new genes with high expression across a desiccation timecourse. Comparison between the Sorghum and Oropetium genomes revealed a surprising degree of chromosome-level collinearity, and several chromosome pairs have near perfect synteny. Other chromosomes are collinear in the gene rich chromosome arms but have experienced pericentric translocations. Together, these resources will be useful for the grass-comparative genomic community and further establish Oropetium as a model resurrection plant.