Genetic and epigenetic variation in transposable element expression responses to abiotic stress in maize

Understanding the Regulation Activities of Transposons in Driving the Variation and Evolution of Polyploid Plant Genome

Transposon is the main component of the eukaryotic genome, and more and more plant genome data show that transposons are diverse in regulating genome structure, variation, function and evolution, with different transposition mechanisms in the genome. Hybridization and polyploidy play an important role in promoting plant speciation and evolution, and recent studies have shown that polyploidy is usually accompanied by the expansion of transposons, which affect the genome size and structure of polyploid plants. Transposons can insert into genes and intergenic regions, resulting in great differences in the overall genome structure of closely related plant species, and it can also capture gene segments in the genome to increase the copy number of genes. In addition, transposons influence the epigenetic modification state of the genome and regulate the expression of the gene, while plant phenotype, biological and abiotic stress response are also regulated by transposons. Overall, transposons play an important role in the plant genome, especially polyploid plant genome, adaptation and evolution.

Epigenetic variation in maize agronomical traits for breeding and trait improvement

Epigenetics-mediated breeding (epibreeding) involves engineering crop traits and stress responses through the targeted manipulation of key epigenetic features to enhance agricultural productivity. While conventional breeding methods raise concerns about reduced genetic diversity, epibreeding propels crop improvement through epigenetic variations that regulate gene expression, ultimately impacting crop yield. Epigenetic regulation in crops encompasses various modes, including histone modification, DNA modification, RNA modification, non-coding RNA, and chromatin remodeling. This review summarizes the epigenetic mechanisms underlying major agronomic traits in maize and identifies candidate epigenetic landmarks in the maize breeding process. We propose a valuable strategy for improving maize yield through epibreeding, combining CRISPR/Cas-based epigenome editing technology and Synthetic Epigenetics (SynEpi). Finally, we discuss the challenges and opportunities associated with maize trait improvement through epibreeding.

Structural Variants Contribute to Phenotypic Variation in Maize

Comprehensively identifying the loci shaping trait variation has been challenging, in part because standard approaches often miss many types of genetic variants. Structural variants (SVs), especially transposable elements (TEs), are likely to affect phenotypic variation but we lack methods that can detect polymorphic SVs and TEs using short-read sequencing data. Here, we used a whole genome alignment between two maize genotypes to identify polymorphic SVs and then genotyped a large maize diversity panel for these variants using short-read sequencing data. After characterising SV variation in the panel, we identified SV polymorphisms that are associated with life history traits and genotype-by-environment (GxE) interactions. While most of the SVs associated with traits contained TEs, only two of the SVs had boundaries that clearly matched TE breakpoints indicative of a TE insertion, while the other polymorphisms were likely caused by deletions. One of the SVs that appeared to be caused by a TE insertion had the most associations with gene expression compared to other trait-associated SVs. All of the SVs associated with traits were in linkage disequilibrium with nearby single nucleotide polymorphisms (SNPs), suggesting that the approach used here did not identify unique associations that would have been missed in a SNP association study. Overall, we have (1) created a technique to genotype SV polymorphisms across a large diversity panel using support from genomic short-read sequencing alignments and (2) connected this presence/absence SV variation to diverse traits and GxE interactions.

Structural variants contribute to phenotypic variation in maize

Comprehensively identifying the loci shaping trait variation has been challenging, in part because standard approaches often miss many types of genetic variants. Structural variants (SVs), especially transposable elements (TEs), are likely to affect phenotypic variation but we lack methods that can detect polymorphic structural variants and TEs using short-read sequencing data. Here, we used a whole genome alignment between two maize genotypes to identify polymorphic structural variants and then genotyped a large maize diversity panel for these variants using short-read sequencing data. After characterizing SV variation in the panel, we identified SV polymorphisms that are associated with life history traits and genotype-by-environment (GxE) interactions. While most of the SVs associated with traits contained TEs, only two of the SVs had boundaries that clearly matched TE breakpoints indicative of a TE insertion, while the other polymorphisms were likely caused by deletions. One of the SVs that appeared to be caused by a TE insertion had the most associations with gene expression compared to other trait-associated SVs. All of the SVs associated with traits were in linkage disequilibrium with nearby single nucleotide polymorphisms (SNPs), suggesting that the approach used here did not identify unique associations that would have been missed in a SNP association study. Overall, we have created a technique to genotype SV polymorphisms across a large diversity panel using support from genomic short-read sequencing alignments and connecting this presence/absence SV variation to diverse traits and GxE interactions.

Unveiling host-parasite relationships through conserved MITEs in prokaryote and viral genomes

Transposable elements (TEs) play a pivotal role in the evolution of genomes across all life domains. 'Miniature Inverted-repeat Transposable-Elements' (MITEs) are non-autonomous TEs mainly located in intergenic regions, relying on external transposases for mobilization. The extent of MITEs' mobilome was explored across nearly 1700 prokaryotic genera, 183 232 genomes, revealing a broad distribution. MITEs were identified in 56.5% of genomes, totaling over 1.4 million cMITEs (cellular MITEs). Cluster analysis revealed that 97.4% of cMITEs were specific within genera boundaries, with up to 23% being species-specific. Subsequently, this genus-specificity was evaluated as a method to link microbial host to their viruses. A total of 51 655 cMITEs had counterparts in viral sequences, termed vMITEs (viral MITEs), resulting in the identification of 2500 viral sequences with them. Among these, 1501 sequences were positively assigned to a previously known host (41.8% were isolated viruses and 12.3% were assigned through CRISPR data), while 379 new host-virus associations were predicted. Deeper analysis in Neisseria and Bacteroidota groups allowed the association of 242 and 530 new viral sequences, respectively. MITEs are proposed as a novel approach to establishing valid virus-host relationships.

Genetic variation in the aquaporin TONOPLAST INTRINSIC PROTEIN 4;3 modulates maize cold tolerance

Cold stress is a major abiotic stress that threatens maize (Zea mays L.) production worldwide. Understanding the molecular mechanisms underlying cold tolerance is crucial for breeding resilient maize varieties. Tonoplast intrinsic proteins (TIPs) are a subfamily of aquaporins in plants. Here, we report that TIP family proteins are involved in maize cold tolerance. The expression of most TIP genes was responsive to cold stress. Overexpressing TIP2;1, TIP3;2 or TIP4;3 reduced the cold tolerance of maize seedlings, while loss-of-function mutants of TIP4;3 exhibited enhanced cold tolerance. Candidate gene-based association analysis revealed that a 328-bp transposon insertion in the promoter region of TIP4;3 was strongly associated with maize cold tolerance. This transposon insertion conferred cold tolerance by repressing TIP4;3 expression through increased methylation of its promoter region. Moreover, TIP4;3 was found to suppress stomatal closure and facilitate reactive oxygen species (ROS) accumulation under cold stress, thereby inhibiting the expression of cold-responsive genes, including DEHYDRATION-RESPONSIVE ELEMENT BINDING FACTOR 1 (DREB1) genes and a subset of peroxidase genes, ultimately attenuating maize cold tolerance. This study thus elucidates the mechanism underlying TIP-mediated cold tolerance and identifies a favourable TIP4;3 allele as a potential genetic resource for breeding cold-tolerant maize varieties.

A Transferable Machine Learning Framework for Predicting Transcriptional Responses of Genes Across Species

Leveraging existing resources in studied species to predict gene functions has the potential to rapidly expand understanding of annotated genes in other, less well-studied, species with assembled genomes. However, orthology is not a reliable predictor for the transcriptional responses of genes to stress. Machine learning methods can quantitatively estimate expression patterns and gene functions using known annotations and collections of features describing each gene. In this chapter, we describe a supervised machine learning framework to predict stress-responsive genes across species using only features derived from nucleotide sequences, using the example of cold stress-responsive genes in different Panicoid grass species.

Mapping responsive genomic elements to heat stress in a maize diversity panel

BACKGROUND

Many plant species exhibit genetic variation for coping with environmental stress. However, there are still limited approaches to effectively uncover the genomic region that regulates distinct responsive patterns of the gene across multiple varieties within the same species under abiotic stress.

RESULTS

By analyzing the transcriptomes of more than 100 maize inbreds, we reveal many cis- and trans-acting eQTLs that influence the expression response to heat stress. The cis-acting eQTLs in response to heat stress are identified in genes with differential responses to heat stress between genotypes as well as genes that are only expressed under heat stress. The cis-acting variants for heat stress-responsive expression likely result from distinct promoter activities, and the differential heat responses of the alleles are confirmed for selected genes using transient expression assays. Global footprinting of transcription factor binding is performed in control and heat stress conditions to document regions with heat-enriched transcription factor binding occupancies.

CONCLUSIONS

Footprints enriched near proximal regions of characterized heat-responsive genes in a large association panel can be utilized for prioritizing functional genomic regions that regulate genotype-specific responses under heat stress.

The evolution and function of transposons in epigenetic regulation in response to the environment

Transposable elements (TEs) make up a major proportion of plant genomes. Despite their prevalence genome-wide, TEs are often tossed aside as "junk DNA" since they rarely cause phenotypes, and epigenetic mechanisms silence TEs to prevent them from causing deleterious mutations through movement. While this bleak picture of TEs in genomes is true on average, a growing number of examples across many plant species point to TEs as drivers of phenotypic diversity and novel stress responses. Examples of TE-influenced phenotypes illustrate the many ways that novel transposition events can alter local gene expression and how this relates to potential variation in plant responses to environmental stress. Since TE families and insertions at the locus level lack evolutionary conservation, advancements in the field will require TE experts across diverse species to identify and utilize TE variation in their own systems as a means of crop improvement.

DNA Methylation in Plant Responses and Adaption to Abiotic Stresses

Due to their sessile state, plants are inevitably affected by and respond to the external environment. So far, plants have developed multiple adaptation and regulation strategies to abiotic stresses. One such system is epigenetic regulation, among which DNA methylation is one of the earliest and most studied regulatory mechanisms, which can regulate genome functioning and induce plant resistance and adaption to abiotic stresses. In this review, we outline the most recent findings on plant DNA methylation responses to drought, high temperature, cold, salt, and heavy metal stresses. In addition, we discuss stress memory regulated by DNA methylation, both in a transient way and the long-term memory that could pass to next generations. To sum up, the present review furnishes an updated account of DNA methylation in plant responses and adaptations to abiotic stresses.

Multi-omics data integration reveals link between epigenetic modifications and gene expression in sugar beet (Beta vulgaris subsp. vulgaris) in response to cold

Background

DNA methylation is thought to influence the expression of genes, especially in response to changing environmental conditions and developmental changes. Sugar beet (Beta vulgaris ssp. vulgaris), and other biennial or perennial plants are inevitably exposed to fluctuating temperatures throughout their lifecycle and might even require such stimulus to acquire floral competence. Therefore, plants such as beets, need to fine-tune their epigenetic makeup to ensure phenotypic plasticity towards changing environmental conditions while at the same time steering essential developmental processes. Different crop species may show opposing reactions towards the same abiotic stress, or, vice versa, identical species may respond differently depending on the specific kind of stress.

Results

In this study, we investigated common effects of cold treatment on genome-wide DNA methylation and gene expression of two Beta vulgaris accessions via multi-omics data analysis. Cold exposure resulted in a pronounced reduction of DNA methylation levels, which particularly affected methylation in CHH context (and to a lesser extent CHG) and was accompanied by transcriptional downregulation of the chromomethyltransferase CMT2 and strong upregulation of several genes mediating active DNA demethylation.

Conclusion

Integration of methylomic and transcriptomic data revealed that, rather than methylation having directly influenced expression, epigenetic modifications correlated with changes in expression of known players involved in DNA (de)methylation. In particular, cold triggered upregulation of genes putatively contributing to DNA demethylation via the ROS1 pathway. Our observations suggest that these transcriptional responses precede the cold-induced global DNA-hypomethylation in non-CpG, preparing beets for additional transcriptional alterations necessary for adapting to upcoming environmental changes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08312-2.

Mechanisms of temperature-regulated growth and thermotolerance in crop species

Temperature is a major environmental factor affecting the development and productivity of crop species. The ability to cope with periods of high temperatures, also known as thermotolerance, is becoming an increasingly indispensable trait for the future of agriculture owing to the current trajectory of average global temperatures. From temperature sensing to downstream transcriptional changes, here, we review recent findings involving the thermal regulation of plant growth and the effects of heat on hormonal pathways, reactive oxygen species, and epigenetic regulation. We also highlight recent approaches and strategies that could be integrated to confront the challenges in sustaining crop productivity in future decades.

Whole-genome variation of transposable element insertions in a maize diversity panel

Intact transposable elements (TEs) account for 65% of the maize genome and can impact gene function and regulation. Although TEs comprise the majority of the maize genome and affect important phenotypes, genome-wide patterns of TE polymorphisms in maize have only been studied in a handful of maize genotypes, due to the challenging nature of assessing highly repetitive sequences. We implemented a method to use short-read sequencing data from 509 diverse inbred lines to classify the presence/absence of 445,418 nonredundant TEs that were previously annotated in four genome assemblies including B73, Mo17, PH207, and W22. Different orders of TEs (i.e., LTRs, Helitrons, and TIRs) had different frequency distributions within the population. LTRs with lower LTR similarity were generally more frequent in the population than LTRs with higher LTR similarity, though high-frequency insertions with very high LTR similarity were observed. LTR similarity and frequency estimates of nested elements and the outer elements in which they insert revealed that most nesting events occurred very near the timing of the outer element insertion. TEs within genes were at higher frequency than those that were outside of genes and this is particularly true for those not inserted into introns. Many TE insertional polymorphisms observed in this population were tagged by SNP markers. However, there were also 19.9% of the TE polymorphisms that were not well tagged by SNPs (R2 < 0.5) that potentially represent information that has not been well captured in previous SNP-based marker-trait association studies. This study provides a population scale genome-wide assessment of TE variation in maize and provides valuable insight on variation in TEs in maize and factors that contribute to this variation.

Predicting transcriptional responses to cold stress across plant species

The same gene is often regulated differently in response to stress in even closely related plant species. Directly measuring stress-responsive gene expression can be financially and logistically challenging in nonmodel species. Here, we show that models trained using data on which genes respond to cold in one species can predict which genes will respond to cold in related species, even when the training and target species vary in their degree of tolerance to cold. The prediction models we used require only genomic sequence and gene models. As a result, data from well-studied model species may be used to predict which genes will respond to stress in less-studied species with sequenced genomes.

Although genome-sequence assemblies are available for a growing number of plant species, gene-expression responses to stimuli have been cataloged for only a subset of these species. Many genes show altered transcription patterns in response to abiotic stresses. However, orthologous genes in related species often exhibit different responses to a given stress. Accordingly, data on the regulation of gene expression in one species are not reliable predictors of orthologous gene responses in a related species. Here, we trained a supervised classification model to identify genes that transcriptionally respond to cold stress. A model trained with only features calculated directly from genome assemblies exhibited only modest decreases in performance relative to models trained by using genomic, chromatin, and evolution/diversity features. Models trained with data from one species successfully predicted which genes would respond to cold stress in other related species. Cross-species predictions remained accurate when training was performed in cold-sensitive species and predictions were performed in cold-tolerant species and vice versa. Models trained with data on gene expression in multiple species provided at least equivalent performance to models trained and tested in a single species and outperformed single-species models in cross-species prediction. These results suggest that classifiers trained on stress data from well-studied species may suffice for predicting gene-expression patterns in related, less-studied species with sequenced genomes.

Stability of DNA methylation and chromatin accessibility in structurally diverse maize genomes

Accessible chromatin and unmethylated DNA are associated with many genes and cis-regulatory elements. Attempts to understand natural variation for accessible chromatin regions (ACRs) and unmethylated regions (UMRs) often rely upon alignments to a single reference genome. This limits the ability to assess regions that are absent in the reference genome assembly and monitor how nearby structural variants influence variation in chromatin state. In this study, de novo genome assemblies for four maize inbreds (B73, Mo17, Oh43, and W22) are utilized to assess chromatin accessibility and DNA methylation patterns in a pan-genome context. A more complete set of UMRs and ACRs can be identified when chromatin data are aligned to the matched genome rather than a single reference genome. While there are UMRs and ACRs present within genomic regions that are not shared between genotypes, these features are 6- to 12-fold enriched within regions between genomes. Characterization of UMRs present within shared genomic regions reveals that most UMRs maintain the unmethylated state in other genotypes with only ∼5% being polymorphic between genotypes. However, the majority (71%) of UMRs that are shared between genotypes only exhibit partial overlaps suggesting that the boundaries between methylated and unmethylated DNA are dynamic. This instability is not solely due to sequence variation as these partially overlapping UMRs are frequently found within genomic regions that lack sequence variation. The ability to compare chromatin properties among individuals with structural variation enables pan-epigenome analyses to study the sources of variation for accessible chromatin and unmethylated DNA.