Genome-wide analysis of horizontal transfer in non-model wild species from a natural ecosystem reveals new insights into genetic exchange in plants.

Horizontal transfer (HT) refers to the exchange of genetic material between divergent species by mechanisms other than reproduction. In recent years, several studies have demonstrated HTs in eukaryotes, particularly in the context of parasitic relationships and in model species. However, very little is known about HT in natural ecosystems, especially those involving non-parasitic wild species, and the nature of the ecological relationships that promote these HTs. In this work, we conducted a pilot study investigating HTs by sequencing the genomes of 17 wild non-model species from a natural ecosystem, the Massane forest, located in southern France. To this end, we developed a new computational pipeline called INTERCHANGE that is able to characterize HTs at the whole genome level without prior annotation and directly in the raw sequencing reads. Using this pipeline, we identified 12 HT events, half of which occurred between lianas and trees. We found that mainly low copy number LTR-retrotransposons from the Copia superfamily were transferred between these wild plant species, especially those of the Ivana and Ale lineages. This study revealed a possible new route for HTs between non-parasitic plants and provides new insights into the genomic characteristics of horizontally transferred DNA in plant genomes.

Diverse and mobile: eccDNA-based identification of carrot low-copy-number LTR retrotransposons active in callus cultures.

Long terminal repeat retrotransposons (LTR-RTs) are mobilized via an RNA intermediate using a 'copy and paste' mechanism, and account for the majority of repetitive DNA in plant genomes. As a side effect of mobilization, the formation of LTR-RT-derived extrachromosomal circular DNAs (eccDNAs) occurs. Thus, high-throughput sequencing of eccDNA can be used to identify active LTR-RTs in plant genomes. Despite the release of a reference genome assembly, carrot LTR-RTs have not yet been thoroughly characterized. LTR-RTs are abundant and diverse in the carrot genome. We identified 5976 carrot LTR-RTs, 2053 and 1660 of which were attributed to Copia and Gypsy superfamilies, respectively. They were further classified into lineages, families and subfamilies. More diverse LTR-RT lineages, i.e. lineages comprising many low-copy-number subfamilies, were more frequently associated with genic regions. Certain LTR-RT lineages have been recently active in Daucus carota. In particular, low-copy-number LTR-RT subfamilies, e.g. those belonging to the DcAle lineage, have significantly contributed to carrot genome diversity as a result of continuing activity. We utilized eccDNA sequencing to identify and characterize two DcAle subfamilies, Alex1 and Alex3, active in carrot callus. We documented 14 and 32 de novo insertions of Alex1 and Alex3, respectively, which were positioned in non-repetitive regions.

Large tandem duplications affect gene expression, 3D organization, and plant-pathogen response.

Rapid plant genome evolution is crucial to adapt to environmental changes. Chromosomal rearrangements and gene copy number variation (CNV) are two important tools for genome evolution and sources for the creation of new genes. However, their emergence takes many generations. In this study, we show that in Arabidopsis thaliana, a significant loss of ribosomal RNA (rRNA) genes with a past history of a mutation for the chromatin assembly factor 1 (CAF1) complex causes rapid changes in the genome structure. Using long-read sequencing and microscopic approaches, we have identified up to 15 independent large tandem duplications in direct orientation (TDDOs) ranging from 60 kb to 1.44 Mb. Our data suggest that these TDDOs appeared within a few generations, leading to the duplication of hundreds of genes. By subsequently focusing on a line only containing 20% of rRNA gene copies (20rDNA line), we investigated the impact of TDDOs on 3D genome organization, gene expression, and cytosine methylation. We found that duplicated genes often accumulate more transcripts. Among them, several are involved in plant-pathogen response, which could explain why the 20rDNA line is hyper-resistant to both bacterial and nematode infections. Finally, we show that the TDDOs create gene fusions and/or truncations and discuss their potential implications for the evolution of plant genomes.

Chromosome-level Thlaspi arvense genome provides new tools for translational research and for a newly domesticated cash cover crop of the cooler climates.

Thlaspi arvense (field pennycress) is being domesticated as a winter annual oilseed crop capable of improving ecosystems and intensifying agricultural productivity without increasing land use. It is a selfing diploid with a short life cycle and is amenable to genetic manipulations, making it an accessible field-based model species for genetics and epigenetics. The availability of a high-quality reference genome is vital for understanding pennycress physiology and for clarifying its evolutionary history within the Brassicaceae. Here, we present a chromosome-level genome assembly of var. MN106-Ref with improved gene annotation and use it to investigate gene structure differences between two accessions (MN108 and Spring32-10) that are highly amenable to genetic transformation. We describe non-coding RNAs, pseudogenes and transposable elements, and highlight tissue-specific expression and methylation patterns. Resequencing of forty wild accessions provided insights into genome-wide genetic variation, and QTL regions were identified for a seedling colour phenotype. Altogether, these data will serve as a tool for pennycress improvement in general and for translational research across the Brassicaceae.

RNAi suppression of DNA methylation affects the drought stress response and genome integrity in transgenic poplar.

Trees are long-lived organisms that continuously adapt to their environments, a process in which epigenetic mechanisms are likely to play a key role. Via downregulation of the chromatin remodeler DECREASED IN DNA METHYLATION 1 (DDM1) in poplar (Populus tremula × Populus alba) RNAi lines, we examined how DNA methylation coordinates genomic and physiological responses to moderate water deficit. We compared the growth and drought response of two RNAi-ddm1 lines to wild-type (WT) trees under well-watered and water deficit/rewatering conditions, and analyzed their methylomes, transcriptomes, mobilomes and phytohormone contents in the shoot apical meristem. The RNAi-ddm1 lines were more tolerant to drought-induced cavitation but did not differ in height or stem diameter growth. About 5000 differentially methylated regions were consistently detected in both RNAi-ddm1 lines, colocalizing with 910 genes and 89 active transposable elements. Under water deficit conditions, 136 differentially expressed genes were found, including many involved in phytohormone pathways; changes in phytohormone concentrations were also detected. Finally, the combination of hypomethylation and drought led to the mobility of two transposable elements. Our findings suggest major roles for DNA methylation in regulation of genes involved in hormone-related stress responses, and the maintenance of genome integrity through repression of transposable elements.

The somatic piRNA pathway controls germline transposition over generations.

Transposable elements (TEs) are parasitic DNA sequences that threaten genome integrity by replicative transposition in host gonads. The Piwi-interacting RNAs (piRNAs) pathway is assumed to maintain Drosophila genome homeostasis by downregulating transcriptional and post-transcriptional TE expression in the ovary. However, the bursts of transposition that are expected to follow transposome derepression after piRNA pathway impairment have not yet been reported. Here, we show, at a genome-wide level, that piRNA loss in the ovarian somatic cells boosts several families of the endogenous retroviral subclass of TEs, at various steps of their replication cycle, from somatic transcription to germinal genome invasion. For some of these TEs, the derepression caused by the loss of piRNAs is backed up by another small RNA pathway (siRNAs) operating in somatic tissues at the post transcriptional level. Derepressed transposition during 70 successive generations of piRNA loss exponentially increases the genomic copy number by up to 10-fold.

Sequencing the extrachromosomal circular mobilome reveals retrotransposon activity in plants.

Retrotransposons are mobile genetic elements abundant in plant and animal genomes. While efficiently silenced by the epigenetic machinery, they can be reactivated upon stress or during development. Their level of transcription not reflecting their transposition ability, it is thus difficult to evaluate their contribution to the active mobilome. Here we applied a simple methodology based on the high throughput sequencing of extrachromosomal circular DNA (eccDNA) forms of active retrotransposons to characterize the repertoire of mobile retrotransposons in plants. This method successfully identified known active retrotransposons in both Arabidopsis and rice material where the epigenome is destabilized. When applying mobilome-seq to developmental stages in wild type rice, we identified PopRice as a highly active retrotransposon producing eccDNA forms in the wild type endosperm. The mobilome-seq strategy opens new routes for the characterization of a yet unexplored fraction of plant genomes.

LTR-TEs abundance, timing and mobility in Solanum commersonii and S. tuberosum genomes following cold-stress conditions.

MAIN CONCLUSION: Copia/Ale is the youngest lineage in both Solanum tuberosum and S. commersonii. Within it, we identified nightshade, a new LTR element active in the cultivated potato. From an evolutionary perspective, long-terminal repeat retrotransposons (LTR-RT) activity during stress may be viewed as a mean by which organisms can keep up rates of genetic adaptation to changing conditions. Potato is one of the most important crop consumed worldwide, but studies on LTR-RT characterization are still lacking. Here, we assessed the abundance, insertion time and activity of LTR-RTs in both cultivated Solanum tuberosum and its cold-tolerant wild relative S. commersonii genomes. Gypsy elements were more abundant than Copia ones, suggesting that the former was somehow more successful in colonizing potato genomes. However, Copia elements, and in particular, the Ale lineage, are younger than Gypsy ones, since their insertion time was in average ~ 2 Mya. Due to the ability of LTR-RTs to be circularized by the host DNA repair mechanisms, we identified via mobilome-seq a Copia/Ale element (called nightshade, informal name used for potato family) active in S. tuberosum genome. Our analyses represent a valuable resource for comparative genomics within the Solanaceae, transposon-tagging and for the design of cultivar-specific molecular markers in potato.

Dicer-like and RNA-dependent RNA polymerase gene family identification and annotation in the cultivated Solanum tuberosum and its wild relative S. commersonii.

MAIN CONCLUSION: We provide advances in DCL and RDR gene diversity in Solanaceae. We also shed light on DCL and RDR gene expression in response to cold stress. DICER-like (DCL) and RNA-dependent RNA polymerase (RDR) genes form the core components to trigger small non-coding RNA (ncRNA) production. In spite of this, little is known about the two gene families in non-model plant species. As their genome sequences are now available, the cultivated potato (Solanum tuberosum) and its cold-tolerant wild relative Solanum commersonii offer a valuable opportunity to advance our understanding of the above genes. To determine the extent of diversification and evolution of DCLs and RDRs in these species, we performed a comparative analysis. Seven DCLs were identified in the two species, whereas seven and six RDR genes were found in S. tuberosum and S. commersonii, respectively. Based on phylogenetic analysis with DCLs and RDRs from several species, we provide evidence for an increase in their number in both potato species. We also disclosed that tandem duplications played a major role in the evolution of these gene families in Solanaceae. DCL and RDR expression was investigated in different tissues and under cold and virus stresses, with divergent profiles of the tandem duplicated genes being found in different tissues. DCL paralogs showed a contrasting expression in S. tuberosum and S. commersonii following cold stress and virus infection. By contrast, no change in RDR transcript activity was detected following both stresses. Overall, this study provides the first comparative genomic analysis of the core components of the RNAi machinery in Solanaceae and offers a scaffold for future functional analysis of these gene families.

Detection of active transposable elements in Arabidopsis thaliana using Oxford Nanopore Sequencing technology.

BACKGROUND: Transposables elements (TEs) contribute to both structural and functional dynamics of most eukaryotic genomes. Because of their propensity to densely populate plant and animal genomes, the precise estimation of the impact of transposition on genomic diversity has been considered as one of the main challenges of today's genomics. The recent development of NGS (next generation sequencing) technologies has open new perspectives in population genomics by providing new methods for high throughput detection of Transposable Elements-associated Structural Variants (TEASV). However, these have relied on Illumina platform that generates short reads (up to 350 nucleotides). This limitation in size of sequence reads can cause high false discovery rate (FDR) and therefore limit the power of detection of TEASVs, especially in the case of large, complex genomes. The newest sequencing technologies, such as Oxford Nanopore Technologies (ONT) can generate kilobases-long reads thus representing a promising tool for TEASV detection in plant and animals. RESULTS: We present the results of a pilot experiment for TEASV detection on the model plant species Arabidopsis thaliana using ONT sequencing and show that it can be used efficiently to detect TE movements. We generated a ~0.8X genome coverage of a met1-derived epigenetic recombinant inbred line (epiRIL) using a MinIon device with R7 chemistry. We were able to detect nine new copies of the LTR-retrotransposon Evadé (EVD). We also evidenced the activity of the DNA transposon CACTA, CAC1. CONCLUSIONS: Even at a low sequence coverage (0.8X), ONT sequencing allowed us to reliably detect several TE insertions in Arabidopsis thaliana genome. The long read length allowed a precise and un-ambiguous mapping of the structural variations caused by the activity of TEs. This suggests that the trade-off between read length and genome coverage for TEASV detection may be in favor of the former. Should the technology be further improved both in terms of lower error rate and operation costs, it could be efficiently used in diversity studies at population level.

Ecological plant epigenetics: Evidence from model and non-model species, and the way forward.

Growing evidence shows that epigenetic mechanisms contribute to complex traits, with implications across many fields of biology. In plant ecology, recent studies have attempted to merge ecological experiments with epigenetic analyses to elucidate the contribution of epigenetics to plant phenotypes, stress responses, adaptation to habitat, and range distributions. While there has been some progress in revealing the role of epigenetics in ecological processes, studies with non-model species have so far been limited to describing broad patterns based on anonymous markers of DNA methylation. In contrast, studies with model species have benefited from powerful genomic resources, which contribute to a more mechanistic understanding but have limited ecological realism. Understanding the significance of epigenetics for plant ecology requires increased transfer of knowledge and methods from model species research to genomes of evolutionarily divergent species, and examination of responses to complex natural environments at a more mechanistic level. This requires transforming genomics tools specifically for studying non-model species, which is challenging given the large and often polyploid genomes of plants. Collaboration among molecular geneticists, ecologists and bioinformaticians promises to enhance our understanding of the mutual links between genome function and ecological processes.

Widespread and frequent horizontal transfers of transposable elements in plants.

Vertical, transgenerational transmission of genetic material occurs through reproduction of living organisms. In addition to vertical inheritance, horizontal gene transfer between reproductively isolated species has recently been shown to be an important, if not dominant, mechanism in the evolution of prokaryotic genomes. In contrast, only a few horizontal transfer (HT) events have been characterized so far in eukaryotes and mainly concern transposable elements (TEs). Whether these are frequent and have a significant impact on genome evolution remains largely unknown. We performed a computational search for highly conserved LTR retrotransposons among 40 sequenced eukaryotic genomes representing the major plant families. We found that 26 genomes (65%) harbor at least one case of horizontal TE transfer (HTT). These transfers concern species as distantly related as palm and grapevine, tomato and bean, or poplar and peach. In total, we identified 32 cases of HTTs, which could translate into more than 2 million among the 13,551 monocot and dicot genera. Moreover, we show that these TEs have remained functional after their transfer, occasionally causing a transpositional burst. This suggests that plants can frequently exchange genetic material through horizontal transfers and that this mechanism may be important in TE-driven genome evolution.

An siRNA pathway prevents transgenerational retrotransposition in plants subjected to stress.

Eukaryotic genomes consist to a significant extent of retrotransposons that are suppressed by host epigenetic mechanisms, preventing their uncontrolled propagation. However, it is not clear how this is achieved. Here we show that in Arabidopsis seedlings subjected to heat stress, a copia-type retrotransposon named ONSEN (Japanese 'hot spring') not only became transcriptionally active but also synthesized extrachromosomal DNA copies. Heat-induced ONSEN accumulation was stimulated in mutants impaired in the biogenesis of small interfering RNAs (siRNAs); however, there was no evidence of transposition occurring in vegetative tissues. After stress, both ONSEN transcripts and extrachromosomal DNA gradually decayed and were no longer detected after 20-30 days. Surprisingly, a high frequency of new ONSEN insertions was observed in the progeny of stressed plants deficient in siRNAs. Insertion patterns revealed that this transgenerational retrotransposition occurred during flower development and before gametogenesis. Therefore in plants with compromised siRNA biogenesis, memory of stress was maintained throughout development, priming ONSEN to transpose during differentiation of generative organs. Retrotransposition was not observed in the progeny of wild-type plants subjected to stress or in non-stressed mutant controls, pointing to a crucial role of the siRNA pathway in restricting retrotransposition triggered by environmental stress. Finally, we found that natural and experimentally induced variants in ONSEN insertions confer heat responsiveness to nearby genes, and therefore mobility bursts may generate novel, stress-responsive regulatory gene networks.

Compromised stability of DNA methylation and transposon immobilization in mosaic Arabidopsis epigenomes.

Transgenerational epigenetic inheritance has been defined by the study of relatively few loci. We examined a population of recombinant inbred lines with epigenetically mosaic chromosomes consisting of wild-type and CG methylation-depleted segments (epiRILs). Surprisingly, transposons that were immobile in the parental lines displayed stochastic movement in 28% of the epiRILs. Although analysis after eight generations of inbreeding, supported by genome-wide DNA methylation profiling, identified recombined parental chromosomal segments, these were interspersed with unexpectedly high frequencies of nonparental methylation polymorphism. Hence, epigenetic inheritance in hybrids derived from parents with divergent epigenomes permits long-lasting epi-allelic interactions that violate Mendelian expectations. Such persistently "unstable" epigenetic states complicate linkage-based epigenomic mapping. Thus, future epigenomic analyses should consider possible genetic instabilities and alternative mapping strategies.

Parent-of-origin control of transgenerational retrotransposon proliferation in Arabidopsis.

Retrotransposons are ubiquitous mobile genetic elements constituting a major part of eukaryotic genomes. Yet, monitoring retrotransposition and subsequent copy number increases in multicellular eukaryotes is intrinsically difficult. By following the transgenerational accumulation of a newly activated retrotransposon EVADE (EVD) in Arabidopsis, we noticed fast expansion of activated elements transmitted through the paternal germ line but suppression when EVD-active copies are maternally inherited. This parent-of-origin effect on EVD proliferation was still observed when gametophytes carried mutations for key epigenetic regulators previously shown to restrict EVD mobility. Therefore, the main mechanism preventing active EVD proliferation seems to act through epigenetic control in sporophytic tissues in the mother plant. In consequence, once activated, this retrotransposon proliferates in plant populations owing to suppressed epigenetic control during paternal transmission. This parental gateway might contribute to the occasional bursts of retrotransposon mobilization deduced from the genome sequences of many plant species.

Selective epigenetic control of retrotransposition in Arabidopsis.

Retrotransposons are mobile genetic elements that populate chromosomes, where the host largely controls their activities. In plants and mammals, retrotransposons are transcriptionally silenced by DNA methylation, which in Arabidopsis is propagated at CG dinucleotides by METHYLTRANSFERASE 1 (MET1). In met1 mutants, however, mobilization of retrotransposons is not observed, despite their transcriptional activation. A post-transcriptional mechanism therefore seems to be preventing retrotransposition. Here we show that a copia-type retrotransposon (Evadé, French for 'fugitive') evaded suppression of its movement during inbreeding of hybrid epigenomes consisting of met1- and wild-type-derived chromosomes. Evadé (EVD) reinsertions caused a series of developmental mutations that allowed its identification. Genetic testing of host control of the EVD life cycle showed that transcriptional suppression occurred by CG methylation supported by RNA-directed DNA methylation. On transcriptional reactivation, subsequent steps of the EVD cycle were inhibited by plant-specific RNA polymerase IV/V and the histone methyltransferase KRYPTONITE (KYP). Moreover, genome resequencing demonstrated retrotransposition of EVD but no other potentially active retroelements when this combination of epigenetic mechanisms was compromised. Our results demonstrate that epigenetic control of retrotransposons extends beyond transcriptional suppression and can be individualized for particular elements.

Loss of DNA methylation affects the recombination landscape in Arabidopsis.

During sexual reproduction, one-half of the genetic material is deposited in gametes, and a complete set of chromosomes is restored upon fertilization. Reduction of the genetic information before gametogenesis occurs in meiosis, when cross-overs (COs) between homologous chromosomes secure an exchange of their genetic information. COs are not evenly distributed along chromosomes and are suppressed in chromosomal regions encompassing compact, hypermethylated centromeric and pericentromeric DNA. Therefore, it was postulated that DNA hypermethylation is inhibitory to COs. Here, when analyzing meiotic recombination in mutant plants with hypomethylated DNA, we observed unexpected and counterintuitive effects of DNA methylation losses on CO distribution. Recombination was further promoted in the hypomethylated chromosome arms while it was inhibited in heterochromatic regions encompassing pericentromeric DNA. Importantly, the total number of COs was not affected, implying that loss of DNA methylation led to a global redistribution of COs along chromosomes. To determine by which mechanisms altered levels of DNA methylation influence recombination--whether directly in cis or indirectly in trans by changing expression of genes encoding recombination components--we analyzed CO distribution in wild-type lines with randomly scattered and well-mapped hypomethylated chromosomal segments. The results of these experiments, supported by expression profiling data, suggest that DNA methylation affects meiotic recombination in cis. Because DNA methylation exhibits significant variation even within a single species, our results imply that it may influence the evolution of plant genomes through the control of meiotic recombination.

Transposable elements, a treasure trove to decipher epigenetic variation: insights from Arabidopsis and crop epigenomes.

In the past decade, plant biologists and breeders have developed a growing interest in the field of epigenetics, which is defined as the study of heritable changes in gene expression that cannot be explained by changes in the DNA sequence. Epigenetic marks can be responsive to the environment, and evolve faster than genetic changes. Therefore, epigenetic diversity may represent an unexplored resource of natural variation that could be used in plant breeding programmes. On the other hand, crop genomes are largely populated with transposable elements (TEs) that are efficiently targeted by epigenetic marks, and part of the epigenetic diversity observed might be explained by TE polymorphisms. Characterizing the degree to which TEs influence epigenetic variation in crops is therefore a major goal to better use epigenetic variation. To date, epigenetic analyses have been mainly focused on the model plant Arabidopsis thaliana, and have provided clues on epigenome features, components that silence pathways, and effects of silencing impairment. But to what extent can Arabidopsis be used as a model for the epigenomics of crops? In this review, we discuss the similarities and differences between the epigenomes of Arabidopsis and crops. We explore the relationship between TEs and epigenomes, focusing on TE silencing control and escape, and the impact of TE mobility on epigenomic variation. Finally, we provide insights into challenges to tackle, and future directions to take in the route towards using epigenetic diversity in plant breeding programmes.

Extrachromosomal circular DNA and structural variants highlight genome instability in Arabidopsis epigenetic mutants.

Abundant extrachromosomal circular DNA (eccDNA) is associated with transposable element (TE) activity. However, how the eccDNA compartment is controlled by epigenetic regulations and what is its impact on the genome is understudied. Here, using long reads, we sequence both the eccDNA compartment and the genome of Arabidopsis thaliana mutant plants affected in DNA methylation and post-transcriptional gene silencing. We detect a high load of TE-derived eccDNA with truncated and chimeric forms. On the genomic side, on top of truncated and full length TE neo-insertions, we detect complex structural variations (SVs) notably at a disease resistance cluster being a natural hotspot of SV. Finally, we serendipitously identify large tandem duplications in hypomethylated plants, suggesting that SVs could have been overlooked in epigenetic mutants. We propose that a high eccDNA load may alter DNA repair pathways leading to genome instability and the accumulation of SVs, at least in plants.

45S rDNA Diversity In Natura as One Step towards Ribosomal Heterogeneity in Arabidopsis thaliana.

The keystone of ribosome biogenesis is the transcription of 45S rDNA. The Arabidopsis thaliana genome contains hundreds of 45S rDNA units; however, they are not all transcribed. Notably, 45S rDNA units contain insertions/deletions revealing the existence of heterogeneous rRNA genes and, likely, heterogeneous ribosomes for rRNAs. In order to obtain an overall picture of 45S rDNA diversity sustaining the synthesis of rRNAs and, subsequently, of ribosomes in natura, we took advantage of 320 new occurrences of Arabidopsis thaliana as a metapopulation named At66, sampled from 0 to 1900 m of altitude in the eastern Pyrenees in France. We found that the 45S rDNA copy number is very dynamic in natura and identified new genotypes for both 5' and 3' External Transcribed Spacers (ETS). Interestingly, the highest 5'ETS genotype diversity is found in altitude while the highest 3'ETS genotype diversity is found at sea level. Structural analysis of 45S rDNA also shows conservation in natura of specific 5'ETS and 3'ETS sequences/features required to control rDNA expression and the processing of rRNAs. In conclusion, At66 is a worthwhile natural laboratory, and unraveled 45S rDNA diversity represents an interesting starting material to select subsets for rDNA transcription and alter the rRNA composition of ribosomes both intra- and inter-site.