Evolutionary Epigenomics of Retrotransposon-Mediated Methylation Spreading in Rice

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

Identification and evolution of a diterpenoid phytoalexin oryzalactone biosynthetic gene in the genus Oryza

The natural variation of plant-specialized metabolites represents the evolutionary adaptation of plants to their environments. However, the molecular mechanisms that account for the diversification of the metabolic pathways have not been fully clarified. Rice plants resist attacks from pathogens by accumulating diterpenoid phytoalexins. It has been confirmed that the composition of rice phytoalexins exhibits numerous natural variations. Major rice phytoalexins (momilactones and phytocassanes) are accumulated in most cultivars, although oryzalactone is a cultivar-specific compound. Here, we attempted to reveal the evolutionary trajectory of the diversification of phytoalexins by analyzing the oryzalactone biosynthetic gene in Oryza species. The candidate gene, KSLX-OL, which accounts for oryzalactone biosynthesis, was found around the single-nucleotide polymorphisms specific to the oryzalactone-accumulating cultivars in the long arm of chromosome 11. The metabolite analyses in Nicotiana benthamiana and rice plants overexpressing KSLX-OL indicated that KSLX-OL is responsible for the oryzalactone biosynthesis. KSLX-OL is an allele of KSL8 that is involved in the biosynthesis of another diterpenoid phytoalexin, oryzalexin S and is specifically distributed in the AA genome species. KSLX-NOL and KSLX-bar, which encode similar enzymes but are not involved in oryzalactone biosynthesis, were also found in AA genome species. The phylogenetic analyses of KSLXs, KSL8s, and related pseudogenes (KSL9s) indicated that KSLX-OL was generated from a common ancestor with KSL8 and KSL9 via gene duplication, functional differentiation, and gene fusion. The wide distributions of KSLX-OL and KSL8 in AA genome species demonstrate their long-term coexistence beyond species differentiation, suggesting a balancing selection between the genes.

Frequency of cytosine methylation in the adjacent regions of soybean retrotransposon SORE-1 depends on chromosomal location

Mobilization of transposable elements (TEs) is suppressed by epigenetic mechanisms involving cytosine methylation. However, few studies have focused on clarifying relationships between epigenetic influences of TEs on the adjacent DNA regions and time after insertion of TEs into the genome and/or their chromosomal location. Here we addressed these issues using soybean retrotransposon SORE-1. We analyzed SORE-1, inserted in exon 1 of the GmphyA2 gene, one of the newest insertions in this family so far identified. Cytosine methylation was detected in this element but was barely present in the adjacent regions. These results were correlated, respectively, with the presence and absence of the production of short interfering RNAs. Cytosine methylation profiles of 74 SORE-1 elements in the Williams 82 reference genome indicated that methylation frequency in the adjacent regions of SORE-1 was profoundly higher in pericentromeric regions than in euchromatic chromosome arms and was only weakly correlated with the length of time after insertion into the genome. Notably, the higher level of methylation in the 5' adjacent regions of SORE-1 coincided with the presence of repetitive elements in pericentromeric regions. Together, these results suggest that epigenetic influence of SORE-1 on the adjacent regions is influenced by its location on the chromosome.

Evolutionary insights into 3D genome organization and epigenetic landscape of Vigna mungo

Eukaryotic genomes show an intricate three-dimensional (3D) organization within the nucleus that regulates multiple biological processes including gene expression. Contrary to animals, understanding of 3D genome organization in plants remains at a nascent stage. Here, we investigate the evolution of 3D chromatin architecture in legumes. By using cutting-edge PacBio, Illumina, and Hi-C contact reads, we report a gap-free, chromosome-scale reference genome assembly of Vigna mungo, an important minor legume cultivated in Southeast Asia. We spatially resolved V. mungo chromosomes into euchromatic, transcriptionally active A compartment and heterochromatic, transcriptionally-dormant B compartment. We report the presence of TAD-like-regions throughout the diagonal of the HiC matrix that resembled transcriptional quiescent centers based on their genomic and epigenomic features. We observed high syntenic breakpoints but also high coverage of syntenic sequences and conserved blocks in boundary regions than in the TAD-like region domains. Our findings present unprecedented evolutionary insights into spatial 3D genome organization and epigenetic patterns and their interaction within the V. mungo genome. This will aid future genomics and epigenomics research and breeding programs of V. mungo.

Genomic variations combined with epigenetic modifications rewire open chromatin in rice

Cis-regulatory elements (CREs) fine-tune gene transcription in eukaryotes. CREs with sequence variations play vital roles in driving plant or crop domestication. However, how global sequence and structural variations (SVs) are responsible for multilevel changes between indica and japonica rice (Oryza sativa) is still not fully elucidated. To address this, we conducted multiomic studies using MNase hypersensitivity sequencing (MH-seq) in combination with RNA sequencing (RNA-seq), chromatin immunoprecipitation sequencing (ChIP-seq), and bisulfite sequencing (BS-seq) between the japonica rice variety Nipponbare (NIP) and indica rice variety 93-11. We found that differential MNase hypersensitive sites (MHSs) exhibited some distinct intrinsic genomic sequence features between NIP and 93-11. Notably, through MHS-genome-wide association studies (GWAS) integration, we found that key sequence variations may be associated with differences of agronomic traits between NIP and 93-11, which is partly achieved by MHSs harboring CREs. In addition, SV-derived differential MHSs caused by transposable element (TE) insertion, especially by noncommon TEs among rice varieties, were associated with genes with distinct functions, indicating that TE-driven gene neo- or subfunctionalization is mediated by changes of chromatin openness. This study thus provides insights into how sequence and genomic SVs control agronomic traits of NIP and 93-11; it also provides genome-editing targets for molecular breeding aiming at improving favorable agronomic properties.

Chromosomal-level genome assembly of the high-quality Xian/Indica rice (Oryza sativa L.) Xiangyaxiangzhan

The indica rice variety XYXZ carries elite traits including appearance and eating quality. Here, we report the de novo assembly of XYXZ using Illumine paired-end whole-genome shotgun sequencing and Nanopore sequencing. We annotated 39,722 protein-coding genes in the 395.04 Mb assembly. In comparison to other cultivars, XYXZ showed a larger gene size including the transcripts and introns, and more exons per gene. And hundreds of ultra-long genes were also detected. A total of 4362 complete LTRs were annotated, and among them, many were located next to or in protein-coding genes including several genes related to rice quality. We observed the different distributions of LTRs in these genes among XYXZ, Nipponbare, and R498, implying these LTRs might potentially affect expressions of the proximal genes and rice quality. Overall, This chromosome-length genome assembly of XYXZ provides a valuable resource for gene discovery, genetic variation and evolution, and the breeding of high-quality rice.

Off-target piRNA gene silencing in Drosophila melanogaster rescued by a transposable element insertion

Transposable elements (TE) are selfish genetic elements that can cause harmful mutations. In Drosophila, it has been estimated that half of all spontaneous visible marker phenotypes are mutations caused by TE insertions. Several factors likely limit the accumulation of exponentially amplifying TEs within genomes. First, synergistic interactions between TEs that amplify their harm with increasing copy number are proposed to limit TE copy number. However, the nature of this synergy is poorly understood. Second, because of the harm posed by TEs, eukaryotes have evolved systems of small RNA-based genome defense to limit transposition. However, as in all immune systems, there is a cost of autoimmunity and small RNA-based systems that silence TEs can inadvertently silence genes flanking TE insertions. In a screen for essential meiotic genes in Drosophila melanogaster, a truncated Doc retrotransposon within a neighboring gene was found to trigger the germline silencing of ald, the Drosophila Mps1 homolog, a gene essential for proper chromosome segregation in meiosis. A subsequent screen for suppressors of this silencing identified a new insertion of a Hobo DNA transposon in the same neighboring gene. Here we describe how the original Doc insertion triggers flanking piRNA biogenesis and local gene silencing. We show that this local gene silencing occurs in cis and is dependent on deadlock, a component of the Rhino-Deadlock-Cutoff (RDC) complex, to trigger dual-strand piRNA biogenesis at TE insertions. We further show how the additional Hobo insertion leads to de-silencing by reducing flanking piRNA biogenesis triggered by the original Doc insertion. These results support a model of TE-mediated gene silencing by piRNA biogenesis in cis that depends on local determinants of transcription. This may explain complex patterns of off-target gene silencing triggered by TEs within populations and in the laboratory. It also provides a mechanism of sign epistasis among TE insertions, illuminates the complex nature of their interactions and supports a model in which off-target gene silencing shapes the evolution of the RDC complex.

Methylation in the CHH Context Allows to Predict Recombination in Rice

DNA methylation is the most studied epigenetic trait. It is considered a key factor in regulating plant development and physiology, and has been associated with the regulation of several genomic features, including transposon silencing, regulation of gene expression, and recombination rates. Nonetheless, understanding the relation between DNA methylation and recombination rates remains a challenge. This work explores the association between recombination rates and DNA methylation for two commercial rice varieties. The results show negative correlations between recombination rates and methylated cytosine counts for all contexts tested at the same time, and for CG and CHG contexts independently. In contrast, a positive correlation between recombination rates and methylated cytosine count is reported in CHH contexts. Similar behavior is observed when considering only methylated cytosines within genes, transposons, and retrotransposons. Moreover, it is shown that the centromere region strongly affects the relationship between recombination rates and methylation. Finally, machine learning regression models are applied to predict recombination using the count of methylated cytosines in the CHH context as the entrance feature. These findings shed light on the understanding of the recombination landscape of rice and represent a reference framework for future studies in rice breeding, genetics, and epigenetics.

Species-specific chromatin landscape determines how transposable elements shape genome evolution

Transposable elements (TEs) are selfish genetic parasites that increase their copy number at the expense of host fitness. The ‘success’, or genome-wide abundance, of TEs differs widely between species. Deciphering the causes for this large variety in TE abundance has remained a central question in evolutionary genomics. We previously proposed that species-specific TE abundance could be driven by the inadvertent consequences of host-direct epigenetic silencing of TEs—the spreading of repressive epigenetic marks from silenced TEs into adjacent sequences. Here, we compared this TE-mediated local enrichment of repressive marks, or ‘the epigenetic effect of TEs’, in six species in the Drosophila melanogaster subgroup to dissect step-by-step the role of such effect in determining genomic TE abundance. We found that TE-mediated local enrichment of repressive marks is prevalent and substantially varies across and even within species. While this TE-mediated effect alters the epigenetic states of adjacent genes, we surprisingly discovered that the transcription of neighboring genes could reciprocally impact this spreading. Importantly, our multi-species analysis provides the power and appropriate phylogenetic resolution to connect species-specific host chromatin regulation, TE-mediated epigenetic effects, the strength of natural selection against TEs, and genomic TE abundance unique to individual species. Our findings point toward the importance of host chromatin landscapes in shaping genome evolution through the epigenetic effects of a selfish genetic parasite.

All the instructions required for life are encoded in the set of DNA present in a cell. It therefore seems natural to think that every bit of this genetic information should serve the organism. And yet most species carry parasitic ‘transposable’ sequences, or transposons, whose only purpose is to multiply and insert themselves at other positions in the genome.

It is possible for cells to suppress these selfish elements. Chemical marks can be deposited onto the DNA to temporarily ‘silence’ transposons and prevent them from being able to move and replicate. However, this sometimes comes at a cost: the repressive chemical modifications can spread to nearby genes that are essential for the organism and perturb their function.

Strangely, the prevalence of transposons varies widely across the tree of life. These sequences form the majority of the genome of certain species – in fact, they represent about half of the human genetic information. But their abundance is much lower in other organisms, forming a measly 6% of the genome of puffer fish for instance. Even amongst fruit fly species, the prevalence of transposable elements can range between 2% and 25%. What explains such differences?

Huang et al. set out to examine this question through the lens of transposon silencing, systematically comparing how this process impacts nearby regions in six species of fruit flies. This revealed variations in the strength of the side effects associated with transposon silencing, resulting in different levels of perturbation on neighbouring genes. A stronger impact was associated with the species having fewer transposons in its genome, suggesting that an evolutionary pressure is at work to keep the abundance of transposons at a low level in these species. Further analyses showed that the genes which determine how silencing marks are distributed may also be responsible for the variations in the impact of transposon silencing. They could therefore be the ones driving differences in the abundance of transposons between species.

Overall, this work sheds light on the complex mechanisms shaping the evolution of genomes, and it may help to better understand how transposons are linked to processes such as aging and cancer.

Characterization of Transposon-Derived Accessible Chromatin Regions in Rice (Oryza Sativa)

Growing evidence indicates that transposons or transposable elements (TEs)-derived accessible chromatin regions (ACRs) play essential roles in multiple biological processes by interacting with trans-acting factors. However, the function of TE-derived ACRs in the regulation of gene expression in the rice genome has not been well characterized. In this study, we examined the chromatin dynamics in six types of rice tissues and found that ~8% of ACRs were derived from TEs and exhibited distinct levels of accessibility and conservation as compared to those without TEs. TEs exhibited a TE subtype-dependent impact on ACR formation, which can be mediated by changes in the underlying DNA methylation levels. Moreover, we found that tissue-specific TE-derived ACRs might function in the tissue development through the modulation of nearby gene expression. Interestingly, many genes in domestication sweeps were found to overlap with TE-derived ACRs, suggesting their potential functions in the rice domestication. In addition, we found that the expression divergence of 1070 duplicate gene pairs were associated with TE-derived ACRs and had distinct distributions of TEs and ACRs around the transcription start sites (TSSs), which may experience different selection pressures. Thus, our study provides some insights into the biological implications of TE-derived ACRs in the rice genome. Our results imply that these ACRs are likely involved in the regulation of tissue development, rice domestication and functional divergence of duplicated genes.

Whole-genome DNA methylation patterns of Oryza sativa (L.) and Oryza glumaepatula (Steud) genotypes associated with aluminum response

Epigenetic mechanisms in crops have emerged as a fundamental factor in plant adaptation and acclimation to biotic and abiotic stresses. Among described epigenetic mechanisms, DNA methylation has been defined as the most studied epigenetic modification involved in several developmental processes. It has been shown that contrasting methylation marks are associated with gene expression variations between cultivated and wild crop species. In this study, we analyzed single-base resolution methylome maps for Oryza sativa (a cultivated species) and Oryza glumaepatula (a wild species) genotypes grown under control conditions. Our results showed that overall, genome-wide methylation profiles are mainly conserved between both species, nevertheless, there are several differentially methylated regions with species-specific methylation patterns. In addition, we analyzed the association of identified DNA methylation marks in relation with Aluminum-tolerance levels of studied genotypes. We found several differentially methylated regions (DMRs) and DMR-associated genes (DAGs) that are linked with Al tolerance. Some of these DAGs have been previously reported as differentially expressed under Al exposure in O. sativa. Complementarily a Transposable Elements (TE) analysis revealed that specific aluminum related genes have associated-TEs potentially regulated by DNA methylation. Interestingly, the DMRs and DAGs between Al-tolerant and susceptible genotypes were different between O. sativa and O. glumaepatula, suggesting that methylation patterns related to Al responses are unique for each rice species. Our findings provide novel insights into DNA methylation patterns in wild and cultivated rice genotypes and their possible role in the regulation of plant stress responses.

TE Density: a tool to investigate the biology of transposable elements

BACKGROUND

Transposable elements (TEs) are powerful creators of genotypic and phenotypic diversity due to their inherent mutagenic capabilities and in this way they serve as a deep reservoir of sequences for genomic variation. As agents of genetic disruption, a TE's potential to impact phenotype is partially a factor of its location in the genome. Previous research has shown TEs' ability to impact the expression of neighboring genes, however our understanding of this trend is hampered by the exceptional amount of diversity in the TE world, and a lack of publicly available computational methods that quantify the presence of TEs relative to genes.

RESULTS

Here, we have developed a tool to more easily quantify TE presence relative to genes through the use of only a gene and TE annotation, yielding a new metric we call TE Density. Briefly defined as the proportion of TE-occupied base-pairs relative to a window-size of the genome. This new pipeline reports TE density for each gene in the genome, for each type descriptor of TE (order and superfamily), and for multiple positions and distances relative to the gene (upstream, intragenic, and downstream) over sliding, user-defined windows. In this way, we overcome previous limitations to the study of TE-gene relationships by focusing on all TE types present in the genome, utilizing flexible genomic distances for measurement, and reporting a TE presence metric for every gene in the genome.

CONCLUSIONS

Together, this new tool opens up new avenues for studying TE-gene relationships, genome architecture, comparative genomics, and the tremendous diversity present of the TE world. TE Density is open-source and freely available at: https://github.com/sjteresi/TE_Density .

The important contribution of transposable elements to phenotypic variation and evolution

Transposable elements (TEs) are responsible for significant genomic variation in plants. Our understanding of the evolutionary forces shaping TE polymorphism has lagged behind other mutations because of the difficulty of accurately identifying TE polymorphism in short-read population genomic data. However, new approaches allow us to quantify TE polymorphisms in population datasets and address fundamental questions about the evolution of these polymorphisms. Here, we discuss how insertional biases shape where, when, and how often TEs insert throughout the genome. Next, we examine mechanisms by which TEs can affect phenotype. Finally, we evaluate current evidence for selection on TE polymorphisms. All together, it is clear that TEs are important, but underappreciated, contributors to intraspecific phenotypic variation, and that understanding the dynamics governing TE polymorphism is crucial for evolutionary biologists interested in the maintenance of variation.

Epigenome and Epitranscriptome: Potential Resources for Crop Improvement

Crop breeding faces the challenge of increasing food demand, especially under climatic changes. Conventional breeding has relied on genetic diversity by combining alleles to obtain desired traits. In recent years, research on epigenetics and epitranscriptomics has shown that epigenetic and epitranscriptomic diversity provides additional sources for crop breeding and harnessing epigenetic and epitranscriptomic regulation through biotechnologies has great potential for crop improvement. Here, we review epigenome and epitranscriptome variations during plant development and in response to environmental stress as well as the available sources for epiallele formation. We also discuss the possible strategies for applying epialleles and epitranscriptome engineering in crop breeding.

Impact of transposable elements on the evolution of complex living systems and their epigenetic control

Transposable elements (TEs) contribute to genomic innovations, as well as genome instability, across a wide variety of species. Popular designations such as 'selfish DNA' and 'junk DNA,' common in the 1980s, may be either inaccurate or misleading, while a more enlightened view of the TE-host relationship covers a range from parasitism to mutualism. Both plant and animal hosts have evolved epigenetic mechanisms to reduce the impact of TEs, both by directly silencing them and by reducing their ability to transpose in the genome. However, TEs have also been co-opted by both plant and animal genomes to perform a variety of physiological functions, ranging from TE-derived proteins acting directly in normal biological functions to innovations in transcription factor activity and also influencing gene expression. Their presence, in fact, can affect a range of features at genome, phenotype, and population levels. The impact TEs have had on evolution is multifaceted, and many aspects still remain unexplored. In this review, the epigenetic control of TEs is contextualized according to the evolution of complex living systems.

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.

Epigenome plasticity in plants

Plant intra-individual and inter-individual variation can be determined by the epigenome, a set of covalent modifications of DNA and chromatin that can alter genome structure and activity without changes to the genome sequence. The epigenome of plant cells is plastic, that is, it can change in response to internal or external cues, such as during development or due to environmental changes, to create a memory of such events. Ongoing advances in technologies to read and write epigenomic patterns with increasing resolution, scale and precision are enabling the extent of plant epigenome variation to be more extensively characterized and functionally interrogated. In this Review, we discuss epigenome dynamics and variation within plants during development and in response to environmental changes, including stress, as well as between plants. We review known or potential functions of such plasticity and emphasize the importance of investigating the causality of epigenomic changes. Finally, we discuss emerging technologies that may underpin future research into plant epigenome plasticity.

Assessing the regulatory potential of transposable elements using chromatin accessibility profiles of maize transposons

Transposable elements (TEs) have the potential to create regulatory variation both through the disruption of existing DNA regulatory elements and through the creation of novel DNA regulatory elements. In a species with a large genome, such as maize, many TEs interspersed with genes create opportunities for significant allelic variation due to TE presence/absence polymorphisms among individuals. We used information on putative regulatory elements in combination with knowledge about TE polymorphisms in maize to identify TE insertions that interrupt existing accessible chromatin regions (ACRs) in B73 as well as examples of polymorphic TEs that contain ACRs among four inbred lines of maize including B73, Mo17, W22, and PH207. The TE insertions in three other assembled maize genomes (Mo17, W22, or PH207) that interrupt ACRs that are present in the B73 genome can trigger changes to the chromatin, suggesting the potential for both genetic and epigenetic influences of these insertions. Nearly 20% of the ACRs located over 2 kb from the nearest gene are located within an annotated TE. These are regions of unmethylated DNA that show evidence for functional importance similar to ACRs that are not present within TEs. Using a large panel of maize genotypes, we tested if there is an association between the presence of TE insertions that interrupt, or carry, an ACR and the expression of nearby genes. While most TE polymorphisms are not associated with expression for nearby genes, the TEs that carry ACRs exhibit enrichment for being associated with higher expression of nearby genes, suggesting that these TEs may contribute novel regulatory elements. These analyses highlight the potential for a subset of TEs to rewire transcriptional responses in eukaryotic genomes.

The Underlying Nature of Epigenetic Variation: Origin, Establishment, and Regulatory Function of Plant Epialleles

In plants, the gene expression and associated phenotypes can be modulated by dynamic changes in DNA methylation, occasionally being fixed in certain genomic loci and inherited stably as epialleles. Epiallelic variations in a population can occur as methylation changes at an individual cytosine position, methylation changes within a stretch of genomic regions, and chromatin changes in certain loci. Here, we focus on methylated regions, since it is unclear whether variations at individual methylated cytosines can serve any regulatory function, and the evidence for heritable chromatin changes independent of genetic changes is limited. While DNA methylation is known to affect and regulate wide arrays of plant phenotypes, most epialleles in the form of methylated regions have not been assigned any biological function. Here, we review how epialleles can be established in plants, serve a regulatory function, and are involved in adaptive processes. Recent studies suggest that most epialleles occur as byproducts of genetic variations, mainly from structural variants and Transposable Element (TE) activation. Nevertheless, epialleles that occur spontaneously independent of any genetic variations have also been described across different plant species. Here, we discuss how epialleles that are dependent and independent of genetic architecture are stabilized in the plant genome and how methylation can regulate a transcription relative to its genomic location.

The epiallelic potential of transposable elements and its evolutionary significance in plants

DNA provides the fundamental framework for heritability, yet heritable trait variation need not be completely ‘hard-wired’ into the DNA sequence. In plants, the epigenetic machinery that controls transposable element (TE) activity, and which includes DNA methylation, underpins most known cases of inherited trait variants that are independent of DNA sequence changes. Here, we review our current knowledge of the extent, mechanisms and potential adaptive contribution of epiallelic variation at TE-containing alleles in this group of species. For the purpose of this review, we focus mainly on DNA methylation, as it provides an easily quantifiable readout of such variation. The picture that emerges is complex. On the one hand, pronounced differences in DNA methylation at TE sequences can either occur spontaneously or be induced experimentally en masse across the genome through genetic means. Many of these epivariants are stably inherited over multiple sexual generations, thus leading to transgenerational epigenetic inheritance. Functional consequences can be significant, yet they are typically of limited magnitude and although the same epivariants can be found in nature, the factors involved in their generation in this setting remain to be determined. On the other hand, moderate DNA methylation variation at TE-containing alleles can be reproducibly induced by the environment, again usually with mild effects, and most of this variation tends to be lost across generations. Based on these considerations, we argue that TE-containing alleles, rather than their inherited epiallelic variants, are the main targets of natural selection. Thus, we propose that the adaptive contribution of TE-associated epivariation, whether stable or not, lies predominantly in its capacity to modulate TE mobilization in response to the environment, hence providing hard-wired opportunities for the flexible exploration of the phenotypic space.

This article is part of the theme issue ‘How does epigenetics influence the course of evolution?’

Transcriptional regulation of dosage compensation in Carica papaya

Sex chromosome evolution results in the disparity in gene content between heterogametic sex chromosomes and creates the need for dosage compensation to counteract the effects of gene dose imbalance of sex chromosomes in males and females. It is not known at which stage of sex chromosome evolution dosage compensation would evolve. We used global gene expression profiling in male and female papayas to assess gene expression patterns of sex-linked genes on the papaya sex chromosomes. By analyzing expression ratios of sex-linked genes to autosomal genes and sex-linked genes in males relative to females, our results showed that dosage compensation was regulated on a gene-by-gene level rather than whole sex-linked region in papaya. Seven genes on the papaya X chromosome exhibited dosage compensation. We further compared gene expression ratios in the two evolutionary strata. Y alleles in the older evolutionary stratum showed reduced expression compared to X alleles, while Y alleles in the younger evolutionary stratum showed elevated expression compared to X alleles. Reduced expression of Y alleles in the older evolutionary stratum might be caused by accumulation of deleterious mutations in regulatory regions or transposable element-mediated methylation spreading. Most X-hemizygous genes exhibited either no or very low expression, suggesting that gene silencing might play a role in maintaining transcriptional balance between females and males.

DNA Modification Patterns within the Transposable Elements of the Fig (Ficus carica L.) Genome

Transposable element activity can be harmful to the host’s genome integrity, but it can also provide selective advantages. One strategy to cope with transposons is epigenetic control through DNA base modifications. We report the non-canonic DNA modification dynamics of fig (Ficus carica L.) by exploiting high-quality genome reference and related N4-methylcytosine (4mC) and N6-methyladenine (6mA) data. Overall, 1.49% of transposon nucleotides showed either 4mC or 6mA modifications: the 4mC/6mA ratio was similar in Class I and Class II transposons, with a prevalence of 4mC, which is comparable to coding genes. Different percentages of 4mC or 6mA were observed among LTR-retrotransposon lineages and sub-lineages. Furthermore, both the Copia and Gypsy retroelements showed higher modification rates in the LTR and coding regions compared with their neighbour regions. Finally, the unconventional methylation of retrotransposons is unrelated to the number of close genes, suggesting that the 4mC and 6mA frequency in LTR-retrotransposons should not be related to transcriptional repression in the adjacency of the element. In conclusion, this study highlighted unconventional DNA modification patterns in fig transposable elements. Further investigations will focus on functional implications, in regards to how modified retroelements affect the expression of neighbouring genes, and whether these epigenetic markers can spread from repeats to genes, shaping the plant phenotype.

Evolutionary Genomics of Structural Variation in Asian Rice (Oryza sativa) Domestication

Structural variants (SVs) are a largely unstudied feature of plant genome evolution, despite the fact that SVs contribute substantially to phenotypes. In this study, we discovered SVs across a population sample of 347 high-coverage, resequenced genomes of Asian rice (Oryza sativa) and its wild ancestor (O. rufipogon). In addition to this short-read data set, we also inferred SVs from whole-genome assemblies and long-read data. Comparisons among data sets revealed different features of genome variability. For example, genome alignment identified a large (∼4.3 Mb) inversion in indica rice varieties relative to japonica varieties, and long-read analyses suggest that ∼9% of genes from the outgroup (O. longistaminata) are hemizygous. We focused, however, on the resequencing sample to investigate the population genomics of SVs. Clustering analyses with SVs recapitulated the rice cultivar groups that were also inferred from SNPs. However, the site-frequency spectrum of each SV type-which included inversions, duplications, deletions, translocations, and mobile element insertions-was skewed toward lower frequency variants than synonymous SNPs, suggesting that SVs may be predominantly deleterious. Among transposable elements, SINE and mariner insertions were found at especially low frequency. We also used SVs to study domestication by contrasting between rice and O. rufipogon. Cultivated genomes contained ∼25% more derived SVs and mobile element insertions than O. rufipogon, indicating that SVs contribute to the cost of domestication in rice. Peaks of SV divergence were enriched for known domestication genes, but we also detected hundreds of genes gained and lost during domestication, some of which were enriched for traits of agronomic interest.

Impact of Transposable Elements on Methylation and Gene Expression across Natural Accessions of Brachypodium distachyon

Transposable elements (TEs) constitute a large fraction of plant genomes and are mostly present in a transcriptionally silent state through repressive epigenetic modifications, such as DNA methylation. TE silencing is believed to influence the regulation of adjacent genes, possibly as DNA methylation spreads away from the TE. Whether this is a general principle or a context-dependent phenomenon is still under debate, pressing for studying the relationship between TEs, DNA methylation, and nearby gene expression in additional plant species. Here, we used the grass Brachypodium distachyon as a model and produced DNA methylation and transcriptome profiles for 11 natural accessions. In contrast to what is observed in Arabidopsis thaliana, we found that TEs have a limited impact on methylation spreading and that only few TE families are associated with a low expression of their adjacent genes. Interestingly, we found that a subset of TE insertion polymorphisms is associated with differential gene expression across accessions. Thus, although not having a global impact on gene expression, distinct TE insertions may contribute to specific gene expression patterns in B. distachyon.

Comparative analysis of genome-wide DNA methylation in Neurospora

DNA methylation is an epigenetic mark that plays an important role in genetic regulation in eukaryotes. Major progress has been made in dissecting the molecular pathways that regulate DNA methylation. Yet, little is known about DNA methylation variation over evolutionary time. Here we present an investigation of the variation of DNA methylation and transposable element (TE) content in species of the filamentous ascomycetes Neurospora. We generated genome-wide DNA methylation data at single-base resolution, together with genomic TE content and gene expression data, of 10 individuals representing five closely related Neurospora species. We found that the methylation levels were low (ranging from 1.3% to 2.5%) and varied among the genomes in a species-specific way. Furthermore, we found that the TEs over 400 bp long were targeted by DNA methylation, and in all genomes, high methylation correlated with low GC, confirming a conserved link between DNA methylation and Repeat Induced Point (RIP) mutations in this group of fungi. Both TE content and DNA methylation pattern showed phylogenetic signal, and the species with the highest TE load (N. crassa) also exhibited the highest methylation level per TE. Our results suggest that DNA methylation is an evolvable trait and indicate that the genomes of Neurospora are shaped by an evolutionary arms race between TEs and host defence.

Double-edged sword: The evolutionary consequences of the epigenetic silencing of transposable elements

Transposable elements (TEs) are genomic parasites that selfishly replicate at the expense of host fitness. Fifty years of evolutionary studies of TEs have concentrated on the deleterious genetic effects of TEs, such as their effects on disrupting genes and regulatory sequences. However, a flurry of recent work suggests that there is another important source of TEs’ harmful effects—epigenetic silencing. Host genomes typically silence TEs by the deposition of repressive epigenetic marks. While this silencing reduces the selfish replication of TEs and should benefit hosts, a picture is emerging that the epigenetic silencing of TEs triggers inadvertent spreading of repressive marks to otherwise expressed neighboring genes, ultimately jeopardizing host fitness. In this Review, we provide a long-overdue overview of the recent genome-wide evidence for the presence and prevalence of TEs’ epigenetic effects, highlighting both the similarities and differences across mammals, insects, and plants. We lay out the current understanding of the functional and fitness consequences of TEs’ epigenetic effects, and propose possible influences of such effects on the evolution of both hosts and TEs themselves. These unique evolutionary consequences indicate that TEs’ epigenetic effect is not only a crucial component of TE biology but could also be a significant contributor to genome function and evolution.

Rice 3D chromatin structure correlates with sequence variation and meiotic recombination rate

Genomes of many eukaryotic species have a defined three-dimensional architecture critical for cellular processes. They are partitioned into topologically associated domains (TADs), defined as regions of high chromatin inter-connectivity. While TADs are not a prominent feature of A. thaliana genome organization, they have been reported for other plants including rice, maize, tomato and cotton and for which TAD formation appears to be linked to transcription and chromatin epigenetic status. Here we show that in the rice genome, sequence variation and meiotic recombination rate correlate with the 3D genome structure. TADs display increased SNP and SV density and higher recombination rate compared to inter-TAD regions. We associate the observed differences with the TAD epigenetic landscape, TE composition and an increased incidence of meiotic crossovers.

Golicz et al. report an increase in single nucleotide polymorphisms and structural variations across and within Topologically Associated Domains (TADs) in the rice genome, which is different to the pattern observed in the human genome. They show that this may be due to epigenetic modifications, transposable elements composition, and meiotic crossovers in the TAD regions.

Genomic history and ecology of the geographic spread of rice

Rice (Oryza sativa) is one of the world's most important food crops, and is comprised largely of japonica and indica subspecies. Here, we reconstruct the history of rice dispersal in Asia using whole-genome sequences of more than 1,400 landraces, coupled with geographic, environmental, archaeobotanical and paleoclimate data. Originating around 9,000 yr ago in the Yangtze Valley, rice diversified into temperate and tropical japonica rice during a global cooling event about 4,200 yr ago. Soon after, tropical japonica rice reached Southeast Asia, where it rapidly diversified, starting about 2,500 yr BP. The history of indica rice dispersal appears more complicated, moving into China around 2,000 yr BP. We also identify extrinsic factors that influence genome diversity, with temperature being a leading abiotic factor. Reconstructing the dispersal history of rice and its climatic correlates may help identify genetic adaptations associated with the spread of a key domesticated species.

Current status of genomic resources on wild relatives of rice

Rice is a food crop of global importance, cultivated in diverse agro-climatic zones of the world. However, in the process of domestication many beneficial alleles have been eroded from the gene pool of the rice cultivated globally and eventually has made it vulnerable to a plethora of stresses. In contrast, the wild relatives of rice, despite being agronomically inferior, have inherited a potential of surviving in a range of geographical habitats. These adaptations enrich them with novel traits that upon introgression to modern cultivated varieties offer tremendous potential of increasing yield and adaptability. But, due to the unavailability of their genetic as well as genomic resources, identification and characterisation of these novel beneficial alleles has been a challenging task. Nevertheless, with the unprecedented surge in the area of conservation genomics, researchers have now shifted their focus towards these natural repositories of beneficial traits. Presently, there are several generic and specialized databases harboring genome-wide information on wild species of rice, and are acting as a useful resource for identification of novel genes and alleles, designing of molecular markers, comparative analysis and evolutionary biology studies. In this review, we introduce the key features of these databases focusing on their utility in rice breeding programs.

Transcriptional regulation of genes bearing intronic heterochromatin in the rice genome

Intronic regions of eukaryotic genomes accumulate many Transposable Elements (TEs). Intronic TEs often trigger the formation of transcriptionally repressive heterochromatin, even within transcription-permissive chromatin environments. Although TE-bearing introns are widely observed in eukaryotic genomes, their epigenetic states, impacts on gene regulation and function, and their contributions to genetic diversity and evolution, remain poorly understood. In this study, we investigated the genome-wide distribution of intronic TEs and their epigenetic states in the Oryza sativa genome, where TEs comprise 35% of the genome. We found that over 10% of rice genes contain intronic heterochromatin, most of which are associated with TEs and repetitive sequences. These heterochromatic introns are longer and highly enriched in promoter-proximal positions. On the other hand, introns also accumulate hypomethylated short TEs. Genes with heterochromatic introns are implicated in various biological functions. Transcription of genes bearing intronic heterochromatin is regulated by an epigenetic mechanism involving the conserved factor OsIBM2, mutation of which results in severe developmental and reproductive defects. Furthermore, we found that heterochromatic introns evolve rapidly compared to non-heterochromatic introns. Our study demonstrates that heterochromatin is a common epigenetic feature associated with actively transcribed genes in the rice genome.

Nanopore sequencing-based genome assembly and evolutionary genomics of circum-basmati rice

Background

The circum-basmati group of cultivated Asian rice (Oryza sativa) contains many iconic varieties and is widespread in the Indian subcontinent. Despite its economic and cultural importance, a high-quality reference genome is currently lacking, and the group’s evolutionary history is not fully resolved. To address these gaps, we use long-read nanopore sequencing and assemble the genomes of two circum-basmati rice varieties.

Results

We generate two high-quality, chromosome-level reference genomes that represent the 12 chromosomes of Oryza. The assemblies show a contig N50 of 6.32 Mb and 10.53 Mb for Basmati 334 and Dom Sufid, respectively. Using our highly contiguous assemblies, we characterize structural variations segregating across circum-basmati genomes. We discover repeat expansions not observed in japonica—the rice group most closely related to circum-basmati—as well as the presence and absence variants of over 20 Mb, one of which is a circum-basmati-specific deletion of a gene regulating awn length. We further detect strong evidence of admixture between the circum-basmati and circum-aus groups. This gene flow has its greatest effect on chromosome 10, causing both structural variation and single-nucleotide polymorphism to deviate from genome-wide history. Lastly, population genomic analysis of 78 circum-basmati varieties shows three major geographically structured genetic groups: Bhutan/Nepal, India/Bangladesh/Myanmar, and Iran/Pakistan.

Conclusion

The availability of high-quality reference genomes allows functional and evolutionary genomic analyses providing genome-wide evidence for gene flow between circum-aus and circum-basmati, describes the nature of circum-basmati structural variation, and reveals the presence/absence variation in this important and iconic rice variety group.

Cytosine methylation of rice mitochondrial DNA from grain and leaf tissues

The rice leaf mitochondrial DNA is  more methylated compared to the rice grain mitochondrial DNA. The old rice leaf mitochondrial DNA has also a higher methylation level than the young rice leaf mitochondrial DNA. The presence of DNA methylation in rice organelles has not been well characterized. We have previously shown that cytosine methylation of chloroplast DNA is different between leaf and grain, and varies between young and old leaves in rice. However, the variation in cytosine methylation of mitochondrial DNA is still poorly characterized. In this study, we have investigated cytosine methylation of mitochondrial DNA in the rice grain and leaf. Based on CpG, CHG, and CHH methylation analyses, the leaf mitochondrial DNA was found to be  more methylated compared to the grain mitochondrial DNA. The methylation of the leaf mitochondrial DNA was also higher in old compared to young leaves. Differences in methylation were observed at different cytosine positions of the mitochondrial DNA between grain and leaf, although there were also positions with a similar level of high methylation in all the tissues examined. The differentially methylated cytosine positions in rice mitochondrial DNA were observed mostly in the intergenic region and in some mitochondrial-specific genes involved in ATP production, transcription, and translation. The functional importance of cytosine methylation in the life cycle of rice mitochondria is still to be determined.

How Stress Facilitates Phenotypic Innovation Through Epigenetic Diversity

Climate adaptation through phenotypic innovation will become the main challenge for plants during global warming. Plants exhibit a plethora of mechanisms to achieve environmental and developmental plasticity by inducing dynamic alterations of gene regulation and by maximizing natural variation through large population sizes. While successful over long evolutionary time scales, most of these mechanisms lack the short-term adaptive responsiveness that global warming will require. Here, we review our current understanding of the epigenetic regulation of plant genomes, with a focus on stress-response mechanisms and transgenerational inheritance. Field and laboratory-scale experiments on plants exposed to stress have revealed a multitude of temporally controlled, mechanistic strategies integrating both genetic and epigenetic changes on the genome level. We analyze inter- and intra-species population diversity to discuss how methylome differences and transposon activation can be harnessed for short-term adaptive efforts to shape co-evolving traits in response to qualitatively new climate conditions and environmental stress.

Monitoring the interplay between transposable element families and DNA methylation in maize

DNA methylation and epigenetic silencing play important roles in the regulation of transposable elements (TEs) in many eukaryotic genomes. A majority of the maize genome is derived from TEs that can be classified into different orders and families based on their mechanism of transposition and sequence similarity, respectively. TEs themselves are highly methylated and it can be tempting to view them as a single uniform group. However, the analysis of DNA methylation profiles in flanking regions provides evidence for distinct groups of chromatin properties at different TE families. These differences among TE families are reproducible in different tissues and different inbred lines. TE families with varying levels of DNA methylation in flanking regions also show distinct patterns of chromatin accessibility and modifications within the TEs. The differences in the patterns of DNA methylation flanking TE families arise from a combination of non-random insertion preferences of TE families, changes in DNA methylation triggered by the insertion of the TE and subsequent selection pressure. A set of nearly 70,000 TE polymorphisms among four assembled maize genomes were used to monitor the level of DNA methylation at haplotypes with and without the TE insertions. In many cases, TE families with high levels of DNA methylation in flanking sequence are enriched for insertions into highly methylated regions. The majority of the >2,500 TE insertions into unmethylated regions result in changes in DNA methylation in haplotypes with the TE, suggesting the widespread potential for TE insertions to condition altered methylation in conserved regions of the genome. This study highlights the interplay between TEs and the methylome of a major crop species.

Evolutionary and ecological functional genomics, from lab to the wild

Plant phenotypes are the result of both genetic and environmental forces that act to modulate trait expression. Over the last few years, numerous approaches in functional genomics and systems biology have led to a greater understanding of plant phenotypic variation and plant responses to the environment. These approaches, and the questions that they can address, have been loosely termed evolutionary and ecological functional genomics (EEFG), and have been providing key insights on how plants adapt and evolve. In particular, by bringing these studies from the laboratory to the field, EEFG studies allow us to gain greater knowledge of how plants function in their natural contexts.