Specificities and Dynamics of Transposable Elements in Land Plants

A strategy for studying epigenetic diversity in natural populations: proof of concept in poplar and oak

In the last 20 years, several techniques have been developed for quantifying DNA methylation, the most studied epigenetic marks in eukaryotes, including the gold standard method, whole-genome bisulfite sequencing (WGBS). WGBS quantifies genome-wide DNA methylation but has several inconveniences rendering it less suitable for population-scale epigenetic studies. The high cost of deep sequencing and the large amounts of data generated prompted us to seek an alternative approach. Restricting studies to parts of the genome would be a satisfactory alternative had there not been a major limitation: the need to select upstream targets corresponding to differentially methylated regions as targets. Given the need to study large numbers of samples, we propose a strategy for investigating DNA methylation variation in natural populations, taking into account the structural complexity of genomes, their size, and their content in unique coding regions versus repeated regions as transposable elements. We first identified regions of highly variable DNA methylation in a subset of genotypes representative of the biological diversity in the population by WGBS. We then analysed the variations of DNA methylation in these targeted regions at the population level by sequencing capture bisulfite (SeqCapBis). The entire strategy was then validated by applying it to another species. Our strategy was developed as a proof of concept on natural populations of two forest species: Populus nigra and Quercus petraea.

DNAJ protein gene expansion mechanism in Panicoideae and PgDNAJ functional identification in pearl millet

KEY MESSAGE

Analyze the evolutionary pattern of DNAJ protein genes in the Panicoideae, including pearl millet, to identify and characterize the biological function of PgDNAJ genes in pearl millet. Global warming has become a major factor threatening food security and human development. It is urgent to analyze the heat-tolerant mechanism of plants and cultivate crops that are adapted to high temperature conditions. The Panicoideae are the second largest subfamily of the Poaceae, widely distributed in warm temperate and tropical regions. Many of these species have been reported to have strong adaptability to high temperature stress, such as pearl millet, foxtail millet and sorghum. The evolutionary differences in DNAJ protein genes among 12 Panicoideae species and 10 other species were identified and analyzed. Among them, 79% of Panicoideae DNAJ protein genes were associated with retrotransposon insertion. Analysis of the DNAJ protein pan-gene family in six pearl millet accessions revealed that the non-core genes contained significantly more TEs than the core genes. By identifying and analyzing the distribution and types of TEs near the DNAJ protein genes, it was found that the insertion of Copia and Gypsy retrotransposons provided the source of expansion for the DNAJ protein genes in the Panicoideae. Based on the analysis of the evolutionary pattern of DNAJ protein genes in Panicoideae, the PgDNAJ was obtained from pearl millet through identification. PgDNAJ reduces the accumulation of reactive oxygen species caused by high temperature by activating ascorbate peroxidase (APX), thereby improving the heat resistance of plants. In summary, these data provide new ideas for mining potential heat-tolerant genes in Panicoideae, and help to improve the heat tolerance of other crops.

Evolutionary insights into the organization of chromatin structure and landscape of transcriptional regulation in plants

Development of complex traits necessitates the functioning and coordination of intricate regulatory networks involving multiple genes. Understanding 3D chromatin structure can facilitate insight into the regulation of gene expression by regulatory elements. This potential, of visualizing the role of chromatin organization in the evolution and function of regulatory elements, remains largely unexplored. Here, we describe new perspectives that arise from the dual considerations of sequence variation of regulatory elements and chromatin structure, with a special focus on whole-genome doubling or polyploidy. We underscore the significance of hierarchical chromatin organization in gene regulation during evolution. In addition, we describe strategies for exploring chromatin organization in future investigations of regulatory evolution in plants, enabling insights into the evolutionary influence of regulatory elements on gene expression and, hence, phenotypes.

Complexity of responses to ionizing radiation in plants, and the impact on interacting biotic factors

In nature, plants are simultaneously exposed to different abiotic (e.g., heat, drought, and salinity) and biotic (e.g., bacteria, fungi, and insects) stresses. Climate change and anthropogenic pressure are expected to intensify the frequency of stress factors. Although plants are well equipped with unique and common defense systems protecting against stressors, they may compromise their growth and development for survival in such challenging environments. Ionizing radiation is a peculiar stress factor capable of causing clustered damage. Radionuclides are both naturally present on the planet and produced by human activities. Natural and artificial radioactivity affects plants on molecular, biochemical, cellular, physiological, populational, and transgenerational levels. Moreover, the fitness of pests, pathogens, and symbionts is concomitantly challenged in radiologically contaminated areas. Plant responses to artificial acute ionizing radiation exposure and laboratory-simulated or field chronic exposure are often discordant. Acute or chronic ionizing radiation exposure may occasionally prime the defense system of plants to better tolerate the biotic stress or could often exhaust their metabolic reserves, making plants more susceptible to pests and pathogens. Currently, these alternatives are only marginally explored. Our review summarizes the available literature on the responses of host plants, biotic factors, and their interaction to ionizing radiation exposure. Such systematic analysis contributes to improved risk assessment in radiologically contaminated areas.

Study of Dispersed Repeats in the Cyanidioschyzon merolae Genome

In this study, we applied the iterative procedure (IP) method to search for families of highly diverged dispersed repeats in the genome of Cyanidioschyzon merolae, which contains over 16 million bases. The algorithm included the construction of position weight matrices (PWMs) for repeat families and the identification of more dispersed repeats based on the PWMs using dynamic programming. The results showed that the C. merolae genome contained 20 repeat families comprising a total of 33,938 dispersed repeats, which is significantly more than has been previously found using other methods. The repeats varied in length from 108 to 600 bp (522.54 bp in average) and occupied more than 72% of the C. merolae genome, whereas previously identified repeats, including tandem repeats, have been shown to constitute only about 28%. The high genomic content of dispersed repeats and their location in the coding regions suggest a significant role in the regulation of the functional activity of the genome.

Profiling genome-wide methylation in two maples: Fine-scale approaches to detection with nanopore technology

DNA methylation is critical to the regulation of transposable elements and gene expression and can play an important role in the adaptation of stress response mechanisms in plants. Traditional methods of methylation quantification rely on bisulfite conversion that can compromise accuracy. Recent advances in long-read sequencing technologies allow for methylation detection in real time. The associated algorithms that interpret these modifications have evolved from strictly statistical approaches to Hidden Markov Models and, recently, deep learning approaches. Much of the existing software focuses on methylation in the CG context, but methylation in other contexts is important to quantify, as it is extensively leveraged in plants. Here, we present methylation profiles for two maple species across the full range of 5mC sequence contexts using Oxford Nanopore Technologies (ONT) long-reads. Hybrid and reference-guided assemblies were generated for two new Acer accessions: Acer negundo (box elder; 65x ONT and 111X Illumina) and Acer saccharum (sugar maple; 93x ONT and 148X Illumina). The ONT reads generated for these assemblies were re-basecalled, and methylation detection was conducted in a custom pipeline with the published Acer references (PacBio assemblies) and hybrid assemblies reported herein to generate four epigenomes. Examination of the transposable element landscape revealed the dominance of LTR Copia elements and patterns of methylation associated with different classes of TEs. Methylation distributions were examined at high resolution across gene and repeat density and described within the broader angiosperm context, and more narrowly in the context of gene family dynamics and candidate nutrient stress genes.

Regulatory mechanism of heat-active retrotransposons by the SET domain protein SUVH2

New transposon insertions are deleterious to genome stability. The RNA-directed DNA methylation (RdDM) pathway evolved to regulate transposon activity via DNA methylation. However, current studies have not yet clearly described the transposition regulation. ONSEN is a heat-activated retrotransposon that is activated at 37°C. The plant-specific SUPPRESSOR OF VARIEGATION 3-9 HOMOLOG (SUVH) family proteins function downstream of the RdDM pathway. The SUVH protein families are linked to TE silencing by two pathways, one through DNA methylation and the other through chromatin remodeling. In this study, we analyzed the regulation of ONSEN activity by SUVH2. We observed that ONSEN transcripts were increased; however, there was no transpositional activity in Arabidopsis suvh2 mutant. The suvh2 mutant produced siRNAs from the ONSEN locus under heat stress, suggesting that siRNAs are involved in suppressing transposition. These results provide new insights into the regulatory mechanisms of retrotransposons that involve siRNA in the RdDM pathway.

High-fidelity (repeat) consensus sequences from short reads using combined read clustering and assembly

BACKGROUND

Despite the many cheap and fast ways to generate genomic data, good and exact genome assembly is still a problem, with especially the repeats being vastly underrepresented and often misassembled. As short reads in low coverage are already sufficient to represent the repeat landscape of any given genome, many read cluster algorithms were brought forward that provide repeat identification and classification. But how can trustworthy, reliable and representative repeat consensuses be derived from unassembled genomes?

RESULTS

Here, we combine methods from repeat identification and genome assembly to derive these robust consensuses. We test several use cases, such as (1) consensus building from clustered short reads of non-model genomes, (2) from genome-wide amplification setups, and (3) specific repeat-centred questions, such as the linked vs. unlinked arrangement of ribosomal genes. In all our use cases, the derived consensuses are robust and representative. To evaluate overall performance, we compare our high-fidelity repeat consensuses to RepeatExplorer2-derived contigs and check, if they represent real transposable elements as found in long reads. Our results demonstrate that it is possible to generate useful, reliable and trustworthy consensuses from short reads by a combination from read cluster and genome assembly methods in an automatable way.

CONCLUSION

We anticipate that our workflow opens the way towards more efficient and less manual repeat characterization and annotation, benefitting all genome studies, but especially those of non-model organisms.

The Genomic Shock Hypothesis: Genetic and Epigenetic Alterations of Transposable Elements after Interspecific Hybridization in Plants

Transposable elements (TEs) are major components of plant genomes with the ability to change their position in the genome or to create new copies of themselves in other positions in the genome. These can cause gene disruption and large-scale genomic alterations, including inversions, deletions, and duplications. Host organisms have evolved a set of mechanisms to suppress TE activity and counter the threat that they pose to genome integrity. These includes the epigenetic silencing of TEs mediated by a process of RNA-directed DNA methylation (RdDM). In most cases, the silencing machinery is very efficient for the vast majority of TEs. However, there are specific circumstances in which TEs can evade such silencing mechanisms, for example, a variety of biotic and abiotic stresses or in vitro culture. Hybridization is also proposed as an inductor of TE proliferation. In fact, the discoverer of the transposons, Barbara McClintock, first hypothesized that interspecific hybridization provides a "genomic shock" that inhibits the TE control mechanisms leading to the mobilization of TEs. However, the studies carried out on this topic have yielded diverse results, showing in some cases a total absence of mobilization or being limited to only some TE families. Here, we review the current knowledge about the impact of interspecific hybridization on TEs in plants and the possible implications of changes in the epigenetic mechanisms.

Epigenetic control of transposons during plant reproduction: From meiosis to hybrid seeds

The regulation of transposable elements (TEs) requires overlapping epigenetic modifications that must be reinforced every cell division and generation. In plants, this is achieved by multiple pathways including small RNAs, DNA methylation, and repressive histone marks that act together to control TE expression and activity throughout the entire life cycle. However, transient TE activation is observed during reproductive transitions as a result of epigenome reprogramming, thus providing windows of opportunity for TE proliferation and epigenetic novelty. Ultimately, these events may originate complex TE-driven transcriptional networks or cell-to-cell communication strategies via mobile small RNAs. In this review, we discuss recent findings and current understanding of TE regulation during sexual plant reproduction, and its implications for fertility, early seed development, and epigenetic inheritance.

DNA hypomethylation of the host tree impairs interaction with mutualistic ectomycorrhizal fungus

Ectomycorrhizas are an intrinsic component of tree nutrition and responses to environmental variations. How epigenetic mechanisms might regulate these mutualistic interactions is unknown. By manipulating the level of expression of the chromatin remodeler DECREASE IN DNA METHYLATION 1 (DDM1) and two demethylases DEMETER-LIKE (DML) in Populus tremula × Populus alba lines, we examined how host DNA methylation modulates multiple parameters of the responses to root colonization with the mutualistic fungus Laccaria bicolor. We compared the ectomycorrhizas formed between transgenic and wild-type (WT) trees and analyzed their methylomes and transcriptomes. The poplar lines displaying lower mycorrhiza formation rate corresponded to hypomethylated overexpressing DML or RNAi-ddm1 lines. We found 86 genes and 288 transposable elements (TEs) differentially methylated between WT and hypomethylated lines (common to both OX-dml and RNAi-ddm1) and 120 genes/1441 TEs in the fungal genome suggesting a host-induced remodeling of the fungal methylome. Hypomethylated poplar lines displayed 205 differentially expressed genes (cis and trans effects) in common with 17 being differentially methylated (cis). Our findings suggest a central role of host and fungal DNA methylation in the ability to form ectomycorrhizas including not only poplar genes involved in root initiation, ethylene and jasmonate-mediated pathways, and immune response but also terpenoid metabolism.

Epigenetic Landscape Is Largely Shaped by Diversiform Transposons in Aegilops tauschii

Transposons (TEs) account for more than 80% of the wheat genome, the highest among all known crop species. They play an important role in shaping the elaborate genomic landscape, which is the key to the speciation of wheat. In this study, we analyzed the association between TEs, chromatin states, and chromatin accessibility in Aegilops tauschii, the D genome donor of bread wheat. We found that TEs contributed to the complex but orderly epigenetic landscape as chromatin states showed diverse distributions on TEs of different orders or superfamilies. TEs also contributed to the chromatin state and openness of potential regulatory elements, affecting the expression of TE-related genes. Some TE superfamilies, such as hAT-Ac, carry active/open chromatin regions. In addition, the histone mark H3K9ac was found to be associated with the accessibility shaped by TEs. These results suggest the role of diversiform TEs in shaping the epigenetic landscape and in gene expression regulation in Aegilops tauschii. This has positive implications for understanding the transposon roles in Aegilops tauschii or the wheat D genome.

Transgressive and parental dominant gene expression and cytosine methylation during seed development in Brassica napus hybrids

KEY MESSAGE

Transcriptomic and epigenomic profiling of gene expression and small RNAs during seed and seedling development reveals expression and methylation dominance levels with implications on early stage heterosis in oilseed rape. The enhanced performance of hybrids through heterosis remains a key aspect in plant breeding; however, the underlying mechanisms are still not fully elucidated. To investigate the potential role of transcriptomic and epigenomic patterns in early expression of hybrid vigor, we investigated gene expression, small RNA abundance and genome-wide methylation in hybrids from two distant Brassica napus ecotypes during seed and seedling developmental stages using next-generation sequencing. A total of 31117, 344, 36229 and 7399 differentially expressed genes, microRNAs, small interfering RNAs and differentially methylated regions were identified, respectively. Approximately 70% of the differentially expressed or methylated features displayed parental dominance levels where the hybrid followed the same patterns as the parents. Via gene ontology enrichment and microRNA-target association analyses during seed development, we found copies of reproductive, developmental and meiotic genes with transgressive and paternal dominance patterns. Interestingly, maternal dominance was more prominent in hypermethylated and downregulated features during seed formation, contrasting to the general maternal gamete demethylation reported during gametogenesis in angiosperms. Associations between methylation and gene expression allowed identification of putative epialleles with diverse pivotal biological functions during seed formation. Furthermore, most differentially methylated regions, differentially expressed siRNAs and transposable elements were in regions that flanked genes without differential expression. This suggests that differential expression and methylation of epigenomic features may help maintain expression of pivotal genes in a hybrid context. Differential expression and methylation patterns during seed formation in an F₁ hybrid provide novel insights into genes and mechanisms with potential roles in early heterosis.

Epigenetic management of self and non-self: lessons from 40 years of transgenic plants

Plant varieties exhibiting unstable or variegated phenotypes, or showing virus recovery have long remained a mystery. It is only with the development of transgenic plants 40 years ago that the epigenetic features underlying these phenomena were elucidated. Indeed, the study of transgenic plants that did not express the introduced sequences revealed that transgene loci sometimes undergo transcriptional gene silencing (TGS) or post-transcriptional gene silencing (PTGS) by activating epigenetic defenses that naturally control transposable elements, duplicated genes or viruses. Even when they do not trigger TGS or PTGS spontaneously, stably expressed transgenes driven by viral promoters set apart from endogenous genes in their epigenetic regulation. As a result, transgenes driven by viral promoters are capable of undergoing systemic PTGS throughout the plant, whereas endogenous genes can only undergo local PTGS in cells where RNA quality control is impaired. Together, these results indicate that the host genome distinguishes self from non-self at the epigenetic level, allowing PTGS to eliminate non-self, and preventing PTGS to become systemic and kill the plant when it is locally activated against deregulated self.

Transposable elements are associated with genome-specific gene expression in bread wheat

INTRODUCTION

Recent studies in wheat emphasized the importance of TEs, which occupy ~85% of the wheat genome, as a major source of intraspecific genetic variation due to their recent activity and involvement in genomic rearrangements. The contribution of TEs to structural and functional variations in bread wheat genes is not fully understood.

METHODS

Here, publicly available RNA-Seq databases of bread wheat were integrated to identify TE insertions within gene bodies (exons\ introns) and assess the impact of TE insertions on gene expression variations of homoeologs gene groups. Overall, 70,818 homoeologs genes were analyzed: 55,170 genes appeared in each one of the three subgenomes (termed ABD), named triads; 12,640 genes appeared in two of the three subgenomes (in A and B only, termed AB; or in A and D only, termed AD; or in B and D only, termed BD);, named dyads; and 3,008 genes underwent duplication in one of the three subgenomes (two copies in: subgenome A, termed AABD; subgenome B, termed ABBD; or subgenome D, termed ABDD), named tetrads.

RESULTS

To this end, we found that ~36% of the 70,818 genes contained at least one TE insertion within the gene body, mostly in triads. Analysis of 14,258 triads revealed that the presence of TE insertion in at least one of the triad genes (7,439 triads) was associated with balanced expression (similar expression levels) between the homoeolog genes. TE insertions within the exon or in the untranslated regions (UTRs) of one or more of the homoeologs in a triad were significantly associated with homoeolog expression bias. Furthermore, we found a statistically significant correlation between the presence\absence of TEs insertions belonging to six TE superfamilies and 17 TE subfamilies and the suppression of a single homoeolog gene. A significant association was observed between the presence of TE insertions from specific superfamilies and the expression of genes that are associated with biotic and abiotic stress responses.

CONCLUSION

Our data strongly indicate that TEs might play a prominent role in controlling gene expression in a genome-specific manner in bread wheat.

Inpactor2: a software based on deep learning to identify and classify LTR-retrotransposons in plant genomes

LTR-retrotransposons are the most abundant repeat sequences in plant genomes and play an important role in evolution and biodiversity. Their characterization is of great importance to understand their dynamics. However, the identification and classification of these elements remains a challenge today. Moreover, current software can be relatively slow (from hours to days), sometimes involve a lot of manual work and do not reach satisfactory levels in terms of precision and sensitivity. Here we present Inpactor2, an accurate and fast application that creates LTR-retrotransposon reference libraries in a very short time. Inpactor2 takes an assembled genome as input and follows a hybrid approach (deep learning and structure-based) to detect elements, filter partial sequences and finally classify intact sequences into superfamilies and, as very few tools do, into lineages. This tool takes advantage of multi-core and GPU architectures to decrease execution times. Using the rice genome, Inpactor2 showed a run time of 5 minutes (faster than other tools) and has the best accuracy and F1-Score of the tools tested here, also having the second best accuracy and specificity only surpassed by EDTA, but achieving 28% higher sensitivity. For large genomes, Inpactor2 is up to seven times faster than other available bioinformatics tools.

Integration of Repeatomic and Cytogenetic Data on Satellite DNA for the Genome Analysis in the Genus Salvia (Lamiaceae)

Within the complicated and controversial taxonomy of cosmopolitan genus Salvia L. (Lamiaceae) are valuable species Salvia officinalis L. and Salvia sclarea L., which are important for the pharmaceutical, ornamental horticulture, food, and perfume industries. Genome organization and chromosome structure of these essential oil species remain insufficiently studied. For the first time, the comparative repeatome analysis of S. officinalis and S. sclarea was performed using the obtained NGS data, RepeatExplorer/TAREAN pipelines and FISH-based chromosome mapping of the revealed satellite DNA families (satDNAs). In repeatomes of these species, LTR retrotransposons made up the majority of their repetitive DNA. Interspecific variations in genome abundance of Class I and Class II transposable elements, ribosomal DNA, and satellite DNA were revealed. Four (S. sclarea) and twelve (S. officinalis) putative satDNAs were identified. Based on patterns of chromosomal distribution of 45S rDNA; 5S rDNA and the revealed satDNAs, karyograms of S. officinalis and S. sclarea were constructed. Promising satDNAs which can be further used as chromosome markers to assess inter- and intraspecific chromosome variability in Salvia karyotypes were determined. The specific localization of homologous satDNA and 45S rDNA on chromosomes of the studied Salvia species confirmed their common origin, which is consistent with previously reported molecular phylogenetic data.

Genome-Wide Profiling of DNA Methylome and Transcriptome Reveals Epigenetic Regulation of Potato Response to DON Stress

Potato is an important food crop that occupies lesser area but has greater production than rice and wheat. However, potato production is affected by numerous biotic and abiotic stresses, among which Fusarium dry rot is a disease that has significant effect on potato production, storage, and processing. However, the role of DNA methylation in regulating potato response to Fusarium toxin deoxynivalenol (DON) stress is still not fully understood. In this study, we performed DNA methylome and transcriptome analyses of potato tubers treated with five concentrations of DON. The global DNA methylation levels in potato tubers treated with different concentrations of DON showed significant changes relative to those in the control. In particular, the 20 ng/ml treatment showed the largest decrease in all three contexts of methylation levels, especially CHH contexts in transposon regions. The differentially methylated region (DMR)-associated differentially expressed genes (DEGs) were significantly enriched in resistance-related metabolic pathways, indicating that DNA methylation plays an essential role in potato response to DON stress. Furthermore, we examined lesions on potato tubers infested with Fusarium after treatment. Furthermore, the potato tubers treated with 5 and 35 ng/ml DON had lesions of significantly smaller diameters than those of the control, indicating that DON stress may induce resistance. We speculate that this may be related to epigenetic memory created after DNA methylation changes. The detailed DNA methylome and transcriptome profiles suggest that DNA methylation plays a vital role in potato disease resistance and has great potential for enhancing potato dry rot resistance.