Long-Read cDNA Sequencing Enables a "Gene-Like" Transcript Annotation of Transposable Elements

Long-read sequencing of extrachromosomal circular DNA and genome assembly of a Solanum lycopersicum breeding line revealed active LTR retrotransposons originating from S. Peruvianum L. introgressions

Transposable elements (TEs) are a major force in the evolution of plant genomes. Differences in the transposition activities and landscapes of TEs can vary substantially, even in closely related species. Interspecific hybridization, a widely employed technique in tomato breeding, results in the creation of novel combinations of TEs from distinct species. The implications of this process for TE transposition activity have not been studied in modern cultivars. In this study, we used nanopore sequencing of extrachromosomal circular DNA (eccDNA) and identified two highly active Ty1/Copia LTR retrotransposon families of tomato (Solanum lycopersicum), called Salsa and Ketchup. Elements of these families produce thousands of eccDNAs under controlled conditions and epigenetic stress. EccDNA sequence analysis revealed that the major parts of eccDNA produced by Ketchup and Salsa exhibited low similarity to the S. lycopersicum genomic sequence. To trace the origin of these TEs, whole-genome nanopore sequencing and de novo genome assembly were performed. We found that these TEs occurred in a tomato breeding line via interspecific introgression from S. peruvianum. Our findings collectively show that interspecific introgressions can contribute to both genetic and phenotypic diversity not only by introducing novel genetic variants, but also by importing active transposable elements from other species.

Forces driving transposable element load variation during Arabidopsis range expansion

Genetic load refers to the accumulated and potentially life-threatening deleterious mutations in populations. Understanding the mechanisms underlying genetic load variation of transposable element (TE) insertion, a major large-effect mutation, during range expansion is an intriguing question in biology. Here, we used 1,115 global natural accessions of Arabidopsis (Arabidopsis thaliana) to study the driving forces of TE load variation during its range expansion. TE load increased with range expansion, especially in the recently established Yangtze River basin population. Effective population size, which explains 62.0% of the variance in TE load, high transposition rate, and selective sweeps contributed to TE accumulation in the expanded populations. We genetically mapped and identified multiple candidate causal genes and TEs, and revealed the genetic architecture of TE load variation. Overall, this study reveals the variation in TE genetic load during Arabidopsis expansion and highlights the causes of TE load variation from the perspectives of both population genetics and quantitative genetics.

Identification of epigenetically regulated genes involved in plant-virus interaction and their role in virus-triggered induced resistance

BACKGROUND

Plant responses to a wide range of stresses are known to be regulated by epigenetic mechanisms. Pathogen-related investigations, particularly against RNA viruses, are however scarce. It has been demonstrated that Arabidopsis thaliana plants defective in some members of the RNA-directed DNA methylation (RdDM) or histone modification pathways presented differential susceptibility to the turnip mosaic virus. In order to identify genes directly targeted by the RdDM-related RNA Polymerase V (POLV) complex and the histone demethylase protein JUMONJI14 (JMJ14) during infection, the transcriptomes of infected mutant and control plants were obtained and integrated with available chromatin occupancy data for various epigenetic proteins and marks.

RESULTS

A comprehensive list of virus-responsive gene candidates to be regulated by the two proteins was obtained. Twelve genes were selected for further characterization, confirming their dynamic regulation during the course of infection. Several epigenetic marks on their promoter sequences were found using in silico data, raising confidence that the identified genes are actually regulated by epigenetic mechanisms. The altered expression of six of these genes in mutants of the methyltransferase gene CURLY LEAF and the histone deacetylase gene HISTONE DEACETYLASE 19 suggests that some virus-responsive genes may be regulated by multiple coordinated epigenetic complexes. A temporally separated multiple plant virus infection experiment in which plants were transiently infected with one virus and then infected by a second one was designed to investigate the possible roles of the identified POLV- and JMJ14-regulated genes in wild-type (WT) plants. Plants that had previously been stimulated with viruses were found to be more resistant to subsequent virus challenge than control plants. Several POLV- and JMJ14-regulated genes were found to be regulated in virus induced resistance in WT plants, with some of them poisoned to be expressed in early infection stages.

CONCLUSIONS

A set of confident candidate genes directly regulated by the POLV and JMJ14 proteins during virus infection was identified, with indications that some of them may be regulated by multiple epigenetic modules. A subset of these genes may also play a role in the tolerance of WT plants to repeated, intermittent virus infections.

Cas9-targeted Nanopore sequencing rapidly elucidates the transposition preferences and DNA methylation profiles of mobile elements in plants

Transposable element insertions (TEIs) are an important source of genomic innovation by contributing to plant adaptation, speciation, and the production of new varieties. The often large, complex plant genomes make identifying TEIs from short reads difficult and expensive. Moreover, rare somatic insertions that reflect mobilome dynamics are difficult to track using short reads. To address these challenges, we combined Cas9-targeted Nanopore sequencing (CANS) with the novel pipeline NanoCasTE to trace both genetically inherited and somatic TEIs in plants. We performed CANS of the EVADÉ (EVD) retrotransposon in wild-type Arabidopsis thaliana and rapidly obtained up to 40× sequence coverage. Analysis of hemizygous T-DNA insertion sites and genetically inherited insertions of the EVD transposon in the ddm1 (decrease in DNA methylation 1) genome uncovered the crucial role of DNA methylation in shaping EVD insertion preference. We also investigated somatic transposition events of the ONSEN transposon family, finding that genes that are downregulated during heat stress are preferentially targeted by ONSENs. Finally, we detected hypomethylation of novel somatic insertions for two ONSENs. CANS and NanoCasTE are effective tools for detecting TEIs and exploring mobilome organization in plants in response to stress and in different genetic backgrounds, as well as screening T-DNA insertion mutants and transgenic plants.

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.

Taming transposable elements in livestock and poultry: a review of their roles and applications

Livestock and poultry play a significant role in human nutrition by converting agricultural by-products into high-quality proteins. To meet the growing demand for safe animal protein, genetic improvement of livestock must be done sustainably while minimizing negative environmental impacts. Transposable elements (TE) are important components of livestock and poultry genomes, contributing to their genetic diversity, chromatin states, gene regulatory networks, and complex traits of economic value. However, compared to other species, research on TE in livestock and poultry is still in its early stages. In this review, we analyze 72 studies published in the past 20 years, summarize the TE composition in livestock and poultry genomes, and focus on their potential roles in functional genomics. We also discuss bioinformatic tools and strategies for integrating multi-omics data with TE, and explore future directions, feasibility, and challenges of TE research in livestock and poultry. In addition, we suggest strategies to apply TE in basic biological research and animal breeding. Our goal is to provide a new perspective on the importance of TE in livestock and poultry genomes.

Long-read direct RNA sequencing reveals epigenetic regulation of chimeric gene-transposon transcripts in Arabidopsis thaliana

Transposable elements (TEs) are accumulated in both intergenic and intragenic regions in plant genomes. Intragenic TEs often act as regulatory elements of associated genes and are also co-transcribed with genes, generating chimeric TE-gene transcripts. Despite the potential impact on mRNA regulation and gene function, the prevalence and transcriptional regulation of TE-gene transcripts are poorly understood. By long-read direct RNA sequencing and a dedicated bioinformatics pipeline, ParasiTE, we investigated the transcription and RNA processing of TE-gene transcripts in Arabidopsis thaliana. We identified a global production of TE-gene transcripts in thousands of A. thaliana gene loci, with TE sequences often being associated with alternative transcription start sites or transcription termination sites. The epigenetic state of intragenic TEs affects RNAPII elongation and usage of alternative poly(A) signals within TE sequences, regulating alternative TE-gene isoform production. Co-transcription and inclusion of TE-derived sequences into gene transcripts impact regulation of RNA stability and environmental responses of some loci. Our study provides insights into TE-gene interactions that contributes to mRNA regulation, transcriptome diversity, and environmental responses in plants.

Effects of Wolbachia on Transposable Element Expression Vary Between Drosophila melanogaster Host Genotypes

Transposable elements (TEs) are repetitive DNA sequences capable of changing position in host genomes, thereby causing mutations. TE insertions typically have deleterious effects but they can also be beneficial. Increasing evidence of the contribution of TEs to adaptive evolution further raises interest in understanding what factors impact TE activity. Based on previous studies associating the bacterial endosymbiont Wolbachia with changes in the abundance of piRNAs, a mechanism for TE repression, and to transposition of specific TEs, we hypothesized that Wolbachia infection would interfere with TE activity. We tested this hypothesis by studying the expression of 14 TEs in a panel of 25 Drosophila melanogaster host genotypes, naturally infected with Wolbachia and annotated for TE insertions. The host genotypes differed significantly in Wolbachia titers inside individual flies, with broad-sense heritability around 20%, and in the number of TE insertions, which depended greatly on TE identity. By removing Wolbachia from the target host genotypes, we generated a panel of 25 pairs of Wolbachia-positive and Wolbachia-negative lines in which we quantified transcription levels for our target TEs. We found variation in TE expression that was dependent on Wolbachia status, TE identity, and host genotype. Comparing between pairs of Wolbachia-positive and Wolbachia-negative flies, we found that Wolbachia removal affected TE expression in 21.1% of the TE-genotype combinations tested, with up to 2.3 times differences in the median level of transcript. Our data show that Wolbachia can impact TE activity in host genomes, underscoring the importance this endosymbiont can have in the generation of genetic novelty in hosts.

Extensive de novo activity stabilizes epigenetic inheritance of CG methylation in Arabidopsis transposons

Cytosine methylation within CG dinucleotides (mCG) can be epigenetically inherited over many generations. Such inheritance is thought to be mediated by a semiconservative mechanism that produces binary present/absent methylation patterns. However, we show here that, in Arabidopsis thaliana h1ddm1 mutants, intermediate heterochromatic mCG is stably inherited across many generations and is quantitatively associated with transposon expression. We develop a mathematical model that estimates the rates of semiconservative maintenance failure and de novo methylation at each transposon, demonstrating that mCG can be stably inherited at any level via a dynamic balance of these activities. We find that DRM2-the core methyltransferase of the RNA-directed DNA methylation pathway-catalyzes most of the heterochromatic de novo mCG, with de novo rates orders of magnitude higher than previously thought, whereas chromomethylases make smaller contributions. Our results demonstrate that stable epigenetic inheritance of mCG in plant heterochromatin is enabled by extensive de novo methylation.

Long-lasting memory of jasmonic acid-dependent immunity requires DNA demethylation and ARGONAUTE1

Stress can have long-lasting impacts on plants. Here we report the long-term effects of the stress hormone jasmonic acid (JA) on the defence phenotype, transcriptome and DNA methylome of Arabidopsis. Three weeks after transient JA signalling, 5-week-old plants retained induced resistance (IR) against herbivory but showed increased susceptibility to pathogens. Transcriptome analysis revealed long-term priming and/or upregulation of JA-dependent defence genes but repression of ethylene- and salicylic acid-dependent genes. Long-term JA-IR was associated with shifts in glucosinolate composition and required MYC2/3/4 transcription factors, RNA-directed DNA methylation, the DNA demethylase ROS1 and the small RNA (sRNA)-binding protein AGO1. Although methylome analysis did not reveal consistent changes in DNA methylation near MYC2/3/4-controlled genes, JA-treated plants were specifically enriched with hypomethylated ATREP2 transposable elements (TEs). Epigenomic characterization of mutants and transgenic lines revealed that ATREP2 TEs are regulated by RdDM and ROS1 and produce 21 nt sRNAs that bind to nuclear AGO1. Since ATREP2 TEs are enriched with sequences from IR-related defence genes, our results suggest that AGO1-associated sRNAs from hypomethylated ATREP2 TEs trans-regulate long-lasting memory of JA-dependent immunity.

Epigenetic Stress and Long-Read cDNA Sequencing of Sunflower (Helianthus annuus L.) Revealed the Origin of the Plant Retrotranscriptome

Transposable elements (TEs) contribute not only to genome diversity but also to transcriptome diversity in plants. To unravel the sources of LTR retrotransposon (RTE) transcripts in sunflower, we exploited a recently developed transposon activation method ('TEgenesis') along with long-read cDNA Nanopore sequencing. This approach allows for the identification of 56 RTE transcripts from different genomic loci including full-length and non-autonomous RTEs. Using the mobilome analysis, we provided a new set of expressed and transpositional active sunflower RTEs for future studies. Among them, a Ty3/Gypsy RTE called SUNTY3 exhibited ongoing transposition activity, as detected by eccDNA analysis. We showed that the sunflower genome contains a diverse set of non-autonomous RTEs encoding a single RTE protein, including the previously described TR-GAG (terminal repeat with the GAG domain) as well as new categories, TR-RT-RH, TR-RH, and TR-INT-RT. Our results demonstrate that 40% of the loci for RTE-related transcripts (nonLTR-RTEs) lack their LTR sequences and resemble conventional eucaryotic genes encoding RTE-related proteins with unknown functions. It was evident based on phylogenetic analysis that three nonLTR-RTEs encode GAG (HadGAG1-3) fused to a host protein. These HadGAG proteins have homologs found in other plant species, potentially indicating GAG domestication. Ultimately, we found that the sunflower retrotranscriptome originated from the transcription of active RTEs, non-autonomous RTEs, and gene-like RTE transcripts, including those encoding domesticated proteins.

Ongoing endeavors to detect mobilization of transposable elements

Transposable elements (TEs) are DNA sequences capable of mobilization from one location to another in the genome. Since the discovery of ‘Dissociation (Dc) locus’ by Barbara McClintock in maize (1), mounting evidence in the era of genomics indicates that a significant fraction of most eukaryotic genomes is composed of TE sequences, involving in various aspects of biological processes such as development, physiology, diseases and evolution. Although technical advances in genomics have discovered numerous functional impacts of TE across species, our understanding of TEs is still ongoing process due to challenges resulted from complexity and abundance of TEs in the genome. In this mini-review, we briefly summarize biology of TEs and their impacts on the host genome, emphasizing importance of understanding TE landscape in the genome. Then, we introduce recent endeavors especially in vivo retrotransposition assays and long read sequencing technology for identifying de novo insertions/TE polymorphism, which will broaden our knowledge of extraordinary relationship between genomic cohabitants and their host.

Comprehensive analysis of both long and short read transcriptomes of a clonal and a seed-propagated model species reveal the prerequisites for transcriptional activation of autonomous and non-autonomous transposons in plants

Background

Transposable element (TE) transcription is a precursor to its mobilisation in host genomes. However, the characteristics of expressed TE loci, the identification of self-competent transposon loci contributing to new insertions, and the genomic conditions permitting their mobilisation remain largely unknown.

Results

Using Vitis vinifera embryogenic callus, we explored the impact of biotic stressors on transposon transcription through the exposure of the callus to live cultures of an endemic grapevine yeast, Hanseniaspora uvarum. We found that only 1.7–2.5% of total annotated TE loci were transcribed, of which 5–10% of these were full-length, and the expressed TE loci exhibited a strong location bias towards expressed genes. These trends in transposon transcription were also observed in RNA-seq data from Arabidopsis thaliana wild-type plants but not in epigenetically compromised Arabidopsis ddm1 mutants. Moreover, differentially expressed TE loci in the grapevine tended to share expression patterns with co-localised differentially expressed genes. Utilising nanopore cDNA sequencing, we found a strong correlation between the inclusion of intronic TEs in gene transcripts and the presence of premature termination codons in these transcripts. Finally, we identified low levels of full-length transcripts deriving from structurally intact TE loci in the grapevine model.

Conclusion

Our observations in two disparate plant models representing clonally and seed propagated plant species reveal a closely connected transcriptional relationship between TEs and co-localised genes, particularly when epigenetic silencing is not compromised. We found that the stress treatment alone was insufficient to induce large-scale full-length transcription from structurally intact TE loci, a necessity for non-autonomous and autonomous mobilisation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00271-5.

Patterns of selection in the evolution of a transposable element

Transposable elements are a major component of most eukaryotic genomes. Here, we present a new approach which allows us to study patterns of natural selection in the evolution of transposable elements over short time scales. The method uses the alignment of all elements with intact gag/pol genes of a transposable element family from a single genome. We predict that the ratio of nonsynonymous to synonymous variants in the alignment should decrease as a function of the frequency of the variants, because elements with nonsynonymous variants that reduce transposition will have fewer progeny. We apply our method to Sirevirus long-terminal repeat retrotransposons that are abundant in maize and other plant species and show that nonsynonymous to synonymous variants declines as variant frequency increases, indicating that negative selection is acting strongly on the Sirevirus genome. The asymptotic value of nonsynonymous to synonymous variants suggests that at least 85% of all nonsynonymous mutations in the transposable element reduce transposition. Crucially, these patterns in nonsynonymous to synonymous variants are only predicted to occur if the gene products from a particular transposable element insertion preferentially promote the transposition of the same insertion. Overall, by using large numbers of intact elements, this study sheds new light on the selective processes that act on transposable elements.

Fast co-evolution of anti-silencing systems shapes the invasiveness of Mu-like DNA transposons in eudicots

Transposable elements (TEs) constitute a major threat to genome stability and are therefore typically silenced by epigenetic mechanisms. In response, some TEs have evolved counteracting systems to suppress epigenetic silencing. In the model plant Arabidopsis thaliana, two such anti-silencing systems have been identified and found to be mediated by the VANC DNA-binding proteins encoded by VANDAL transposons. Here, we show that anti-silencing systems have rapidly diversified since their origin in eudicots by gaining and losing VANC-containing domains, such as DUF1985, DUF287, and Ulp1, as well as target sequence motifs. We further demonstrate that these motifs determine anti-silencing specificity by sequence, density, and helical periodicity. Moreover, such rapid diversification yielded at least 10 distinct VANC-induced anti-silencing systems in Arabidopsis. Strikingly, anti-silencing of non-autonomous VANDALs, which can act as reservoirs of 24-nt small RNAs, is critical to prevent the demise of cognate autonomous TEs and to ensure their propagation. Our findings illustrate how complex co-evolutionary dynamics between TEs and host suppression pathways have shaped the emergence of new epigenetic control mechanisms.

RNA sequencing and its applications in cancer and rare diseases

With the invention of RNA sequencing over a decade ago, diagnosis and identification of the gene-related diseases entered a new phase that enabled more accurate analysis of the diseases that are difficult to approach and analyze. RNA sequencing has availed in-depth study of transcriptomes in different species and provided better understanding of rare diseases and taxonomical classifications of various eukaryotic organisms. Development of single-cell, short-read, long-read and direct RNA sequencing using both blood and biopsy specimens of the organism together with recent advancement in computational analysis programs has made the medical professional’s ability in identifying the origin and cause of genetic disorders indispensable. Altogether, such advantages have evolved the treatment design since RNA sequencing can detect the resistant genes against the existing therapies and help medical professions to take a further step in improving methods of treatments towards higher effectiveness and less side effects. Therefore, it is of essence to all researchers and scientists to have deeper insight in all available methods of RNA sequencing while taking a step-in therapy design.

Cis-regulatory sequences in plants: Their importance, discovery, and future challenges

The identification and characterization of cis-regulatory DNA sequences and how they function to coordinate responses to developmental and environmental cues is of paramount importance to plant biology. Key to these regulatory processes are cis-regulatory modules (CRMs), which include enhancers and silencers. Despite the extraordinary advances in high-quality sequence assemblies and genome annotations, the identification and understanding of CRMs, and how they regulate gene expression, lag significantly behind. This is especially true for their distinguishing characteristics and activity states. Here, we review the current knowledge on CRMs and breakthrough technologies enabling identification, characterization, and validation of CRMs; we compare the genomic distributions of CRMs with respect to their target genes between different plant species, and discuss the role of transposable elements harboring CRMs in the evolution of gene expression. This is an exciting time to study cis-regulomes in plants; however, significant existing challenges need to be overcome to fully understand and appreciate the role of CRMs in plant biology and in crop improvement.

Tracing Mobile DNAs: From Molecular to Population Scales

Transposable elements (TEs, transposons) are mobile DNAs that are prevalent in most eukaryotic genomes. In plants, their mobility has vastly contributed to genetic diversity which is essential for adaptive changes and evolution of a species. Such mobile nature of transposon has been also actively exploited in plant science research by generating genetic mutants in non-model plant systems. On the other hand, transposon mobilization can bring about detrimental effects to host genomes and they are therefore mostly silenced by the epigenetic mechanisms. TEs have been studied as major silencing targets and acted a main feature in the remarkable growth of the plant epigenetics field. Despite the importance of transposon in plant biology and biotechnology, their mobilization and the underlying mechanisms are largely left unanswered. This is mainly because of the sequence repetitiveness of transposons, which makes their detection and analyses difficult and complicated. Recently, some attempts have been made to develop new experimental methods detecting active transposons and their mobilization behavior. These techniques reveal TE mobility in various levels, including the molecular, cellular, organismal and population scales. In this review, we will highlight the novel technical approaches in the study of mobile genetic elements and discuss how these techniques impacted on the advancement of transposon research and broadened our understanding of plant genome plasticity.

Assembling Plant Genomes with Long-Read Sequencing

Continuous improvements in long-read sequencing allow us to tackle increasingly big and complex genomes. Here we present the principles of long-read genome assembly, taking Solanum pennellii nanopore sequencing as an example.

Transposons Hidden in Arabidopsis thaliana Genome Assembly Gaps and Mobilization of Non-Autonomous LTR Retrotransposons Unravelled by Nanotei Pipeline

Long-read data is a great tool to discover new active transposable elements (TEs). However, no ready-to-use tools were available to gather this information from low coverage ONT datasets. Here, we developed a novel pipeline, nanotei, that allows detection of TE-contained structural variants, including individual TE transpositions. We exploited this pipeline to identify TE insertion in the Arabidopsis thaliana genome. Using nanotei, we identified tens of TE copies, including ones for the well-characterized ONSEN retrotransposon family that were hidden in genome assembly gaps. The results demonstrate that some TEs are inaccessible for analysis with the current A. thaliana (TAIR10.1) genome assembly. We further explored the mobilome of the ddm1 mutant with elevated TE activity. Nanotei captured all TEs previously known to be active in ddm1 and also identified transposition of non-autonomous TEs. Of them, one non-autonomous TE derived from (AT5TE33540) belongs to TR-GAG retrotransposons with a single open reading frame (ORF) encoding the GAG protein. These results provide the first direct evidence that TR-GAGs and other non-autonomous LTR retrotransposons can transpose in the plant genome, albeit in the absence of most of the encoded proteins. In summary, nanotei is a useful tool to detect active TEs and their insertions in plant genomes using low-coverage data from Nanopore genome sequencing.

Histone H1 prevents non-CG methylation-mediated small RNA biogenesis in Arabidopsis heterochromatin

Flowering plants utilize small RNA (sRNA) molecules to guide DNA methyltransferases to genomic sequences. This RNA-directed DNA methylation (RdDM) pathway preferentially targets euchromatic transposable elements. However, RdDM is thought to be recruited by methylation of histone H3 at lysine 9 (H3K9me), a hallmark of heterochromatin. How RdDM is targeted to euchromatin despite an affinity for H3K9me is unclear. Here, we show that loss of histone H1 enhances heterochromatic RdDM, preferentially at nucleosome linker DNA. Surprisingly, this does not require SHH1, the RdDM component that binds H3K9me. Furthermore, H3K9me is dispensable for RdDM, as is CG DNA methylation. Instead, we find that non-CG methylation is specifically associated with sRNA biogenesis, and without H1 sRNA production quantitatively expands to non-CG-methylated loci. Our results demonstrate that H1 enforces the separation of euchromatic and heterochromatic DNA methylation pathways by excluding the sRNA-generating branch of RdDM from non-CG-methylated heterochromatin.

Cells adapt to different roles by turning different groups of genes on and off. One way cells control which genes are on or off is by creating regions of active and inactive DNA, which are created and maintained by different groups of proteins. Genes in active DNA regions can be turned on, while genes in inactive regions are switched off or silenced. Silenced DNA regions also turn off ‘transposable elements’: pieces of DNA that can copy themselves and move to other regions of the genome if they become active. Transposons can be dangerous if they are activated, because they can disrupt genes or regulatory sequences when they move.

There are different types of active and inactive DNA, but it is not always clear why these differences exist, or how they are maintained over time. In plants, such as the commonly-studied weed Arabidopsis thaliana, there are two types of inactive DNA, called E and H, that can silence transposons. In both types, DNA has small chemicals called methyl groups attached to it, which help inactivate the DNA. Type E DNA is methylated by a process called RNA-directed DNA methylation (RdDM), but RdDM is rarely seen in type H DNA.

Choi, Lyons and Zilberman showed that RdDM is attracted to E and H regions by previously existing methylated DNA. However, in the H regions, a protein called histone H1 blocks RdDM from attaching methyl groups. This helps focus RdDM onto E regions where it is most needed, because E regions contain the types of transposons RdDM is best suited to silence.

When Choi, Lyons and Zilberman examined genetically modified A. thaliana plants that do not produce histone H1, they found that RdDM happened in both E and H regions. There are many more H regions than E regions, so stretching RdDM across both made it less effective at silencing DNA.

This work shows how different DNA silencing processes are focused onto specific genetic regions, helping explain why there are different types of active and inactive DNA within cells. RdDM has been studied as a way to affect crop growth and yield by altering DNA methylation. These results may help such studies by explaining how RdDM is naturally targeted.

Regulation of retrotransposition in Arabidopsis

Plant genomes are largely comprised of retrotransposons which can replicate through 'copy and paste' mechanisms. Long terminal repeat (LTR) retrotransposons are the major class of retrotransposons in plant species, and importantly they broadly affect the expression of nearby genes. Although most LTR retrotransposons are non-functional, active retrotranspositions have been reported in plant species or mutants under normal growth condition and environmental stresses. With the well-defined reference genome and numerous mutant alleles, Arabidopsis studies have significantly expanded our understanding of retrotransposon regulation. Active LTR retrotransposon loci produce virus-like particles to perform reverse transcription, and their complementary DNA can be inserted into new genomic loci. Due to the detrimental consequences of retrotransposition, plants like animals, have developed transcriptional and post-transcriptional silencing mechanisms. Recently several different genome-wide techniques have been developed to understand LTR retrotransposition in Arabidopsis and different plant species. Transposome, methylome, transcriptome, translatome and small RNA sequencing data have revealed how host silencing mechanisms can affect multiple steps of retrotransposition. These recent advances shed light on future mechanistic studies of retrotransposition as well as retrotransposon diversity.

The Evolutionary Volte-Face of Transposable Elements: From Harmful Jumping Genes to Major Drivers of Genetic Innovation

Transposable elements (TEs) are self-replicating DNA elements that constitute major fractions of eukaryote genomes. Their ability to transpose can modify the genome structure with potentially deleterious effects. To repress TE activity, host cells have developed numerous strategies, including epigenetic pathways, such as DNA methylation or histone modifications. Although TE neo-insertions are mostly deleterious or neutral, they can become advantageous for the host under specific circumstances. The phenomenon leading to the appropriation of TE-derived sequences by the host is known as TE exaptation or co-option. TE exaptation can be of different natures, through the production of coding or non-coding DNA sequences with ultimately an adaptive benefit for the host. In this review, we first give new insights into the silencing pathways controlling TE activity. We then discuss a model to explain how, under specific environmental conditions, TEs are unleashed, leading to a TE burst and neo-insertions, with potential benefits for the host. Finally, we review our current knowledge of coding and non-coding TE exaptation by providing several examples in various organisms and describing a method to identify TE co-option events.

The Dynamism of Transposon Methylation for Plant Development and Stress Adaptation

Plant development processes are regulated by epigenetic alterations that shape nuclear structure, gene expression, and phenotypic plasticity; these alterations can provide the plant with protection from environmental stresses. During plant growth and development, these processes play a significant role in regulating gene expression to remodel chromatin structure. These epigenetic alterations are mainly regulated by transposable elements (TEs) whose abundance in plant genomes results in their interaction with genomes. Thus, TEs are the main source of epigenetic changes and form a substantial part of the plant genome. Furthermore, TEs can be activated under stress conditions, and activated elements cause mutagenic effects and substantial genetic variability. This introduces novel gene functions and structural variation in the insertion sites and primarily contributes to epigenetic modifications. Altogether, these modifications indirectly or directly provide the ability to withstand environmental stresses. In recent years, many studies have shown that TE methylation plays a major role in the evolution of the plant genome through epigenetic process that regulate gene imprinting, thereby upholding genome stability. The induced genetic rearrangements and insertions of mobile genetic elements in regions of active euchromatin contribute to genome alteration, leading to genomic stress. These TE-mediated epigenetic modifications lead to phenotypic diversity, genetic variation, and environmental stress tolerance. Thus, TE methylation is essential for plant evolution and stress adaptation, and TEs hold a relevant military position in the plant genome. High-throughput techniques have greatly advanced the understanding of TE-mediated gene expression and its associations with genome methylation and suggest that controlled mobilization of TEs could be used for crop breeding. However, development application in this area has been limited, and an integrated view of TE function and subsequent processes is lacking. In this review, we explore the enormous diversity and likely functions of the TE repertoire in adaptive evolution and discuss some recent examples of how TEs impact gene expression in plant development and stress adaptation.

Genomic approaches to trace the history of human brain evolution with an emerging opportunity for transposon profiling of ancient humans

Transposable elements (TEs) significantly contribute to shaping the diversity of the human genome, and lines of evidence suggest TEs as one of driving forces of human brain evolution. Existing computational approaches, including cross-species comparative genomics and population genetic modeling, can be adapted for the study of the role of TEs in evolution. In particular, diverse ancient and archaic human genome sequences are increasingly available, allowing reconstruction of past human migration events and holding the promise of identifying and tracking TEs among other evolutionarily important genetic variants at an unprecedented spatiotemporal resolution. However, highly degraded short DNA templates and other unique challenges presented by ancient human DNA call for major changes in current experimental and computational procedures to enable the identification of evolutionarily important TEs. Ancient human genomes are valuable resources for investigating TEs in the evolutionary context, and efforts to explore ancient human genomes will potentially provide a novel perspective on the genetic mechanism of human brain evolution and inspire a variety of technological and methodological advances. In this review, we summarize computational and experimental approaches that can be adapted to identify and validate evolutionarily important TEs, especially for human brain evolution. We also highlight strategies that leverage ancient genomic data and discuss unique challenges in ancient transposon genomics.

Reinforcement of transcriptional silencing by a positive feedback between DNA methylation and non-coding transcription

Non-coding transcription is an important determinant of heterochromatin formation. In Arabidopsis thaliana a specialized RNA polymerase V (Pol V) transcribes pervasively and produces long non-coding RNAs. These transcripts work with small interfering RNA to facilitate locus-specific establishment of RNA-directed DNA methylation (RdDM). Subsequent maintenance of RdDM is associated with elevated levels of Pol V transcription. However, the impact of DNA methylation on Pol V transcription remained unresolved. We found that DNA methylation strongly enhances Pol V transcription. The level of Pol V transcription is reduced in mutants defective in RdDM components working downstream of Pol V, indicating that RdDM is maintained by a mutual reinforcement of DNA methylation and Pol V transcription. Pol V transcription is affected only on loci that lose DNA methylation in all sequence contexts in a particular mutant, including mutants lacking maintenance DNA methyltransferases, which suggests that RdDM works in a complex crosstalk with other silencing pathways.

Expression Quantitative Trait Loci (eQTLs) Associated with Retrotransposons Demonstrate their Modulatory Effect on the Transcriptome

Transposable elements (TEs) are repetitive elements that belong to a variety of functional classes and have an important role in shaping genome evolution. Around 50% of the human genome contains TEs, and they have been termed the "dark matter" of the genome because relatively little is known about their function. While TEs have been shown to participate in aberrant gene regulation and the pathogenesis of diseases, only a few studies have explored the systemic effect of TEs on gene expression. In the present study, we analysed whole genome sequences and blood whole transcriptome data from 570 individuals within the Parkinson's Progressive Markers Initiative (PPMI) cohort to identify expression quantitative trait loci (eQTL) regulating genome-wide gene expression associated with TEs. We identified 2132 reference TEs that were polymorphic for their presence or absence in our study cohort. The presence or absence of the TE element could change the expression of the gene or gene clusters from zero to tens of thousands of copies of RNA. The main finding is that many TEs possess very strong regulatory effects, and they have the potential to modulate large genetic networks with hundreds of target genes over the genome. We illustrate the plethora of regulatory mechanisms using examples of their action at the HLA gene cluster and data showing different TEs' convergence to modulate WFS1 gene expression. In conclusion, the presence or absence of polymorphisms of TEs has an eminent genome-wide regulatory function with large effect size at the level of the whole transcriptome. The role of TEs in explaining, in part, the missing heritability for complex traits is convincing and should be considered.

TrancriptomeReconstructoR: data-driven annotation of complex transcriptomes

BACKGROUND

The quality of gene annotation determines the interpretation of results obtained in transcriptomic studies. The growing number of genome sequence information calls for experimental and computational pipelines for de novo transcriptome annotation. Ideally, gene and transcript models should be called from a limited set of key experimental data.

RESULTS

We developed TranscriptomeReconstructoR, an R package which implements a pipeline for automated transcriptome annotation. It relies on integrating features from independent and complementary datasets: (i) full-length RNA-seq for detection of splicing patterns and (ii) high-throughput 5' and 3' tag sequencing data for accurate definition of gene borders. The pipeline can also take a nascent RNA-seq dataset to supplement the called gene model with transient transcripts. We reconstructed de novo the transcriptional landscape of wild type Arabidopsis thaliana seedlings and Saccharomyces cerevisiae cells as a proof-of-principle. A comparison to the existing transcriptome annotations revealed that our gene model is more accurate and comprehensive than the most commonly used community gene models, TAIR10 and Araport11 for A.thaliana and SacCer3 for S.cerevisiae. In particular, we identify multiple transient transcripts missing from the existing annotations. Our new annotations promise to improve the quality of A.thaliana and S.cerevisiae genome research.

CONCLUSIONS

Our proof-of-concept data suggest a cost-efficient strategy for rapid and accurate annotation of complex eukaryotic transcriptomes. We combine the choice of library preparation methods and sequencing platforms with the dedicated computational pipeline implemented in the TranscriptomeReconstructoR package. The pipeline only requires prior knowledge on the reference genomic DNA sequence, but not the transcriptome. The package seamlessly integrates with Bioconductor packages for downstream analysis.

Male gametophyte development in flowering plants: A story of quarantine and sacrifice

Over 160 years ago, scientists made the first microscopic observations of angiosperm pollen. Unlike in animals, male meiosis in angiosperms produces a haploid microspore that undergoes one asymmetric division to form a vegetative cell and a generative cell. These two cells have distinct fates: the vegetative cell exits the cell cycle and elongates to form a tip-growing pollen tube; the generative cell divides once more in the pollen grain or within the growing pollen tube to form a pair of sperm cells. The concept that male germ cells are less active than the vegetative cell came from biochemical analyses and pollen structure anatomy early in the last century and is supported by the pollen transcriptome data of the last decade. However, the mechanism of how and when the transcriptional repression in male germ cells occurs is still not fully understood. In this review, we provide a brief account of the cytological and metabolic differentiation between the vegetative cell and male germ cells, with emphasis on the role of temporary callose walls, dynamic nuclear pore density, transcription repression, and histone variants. We further discuss the intercellular movement of small interfering RNA (siRNA) derived from transposable elements (TEs) and reexamine the function of TE expression in male germ cells.

Recent advancement of NGS technologies to detect active transposable elements in plants

BACKGROUND

Unlike peoples' belief that transposable elements (TEs) are "junk DNAs" or "genomic parasites", TEs are essential genomic elements that bring about genetic diversity and enable evolution of a species. In fact, transposons are major constituent of chromosome in crop genomes, particularly in major cereal crops, the primary type of which is long terminal repeat (LTR) retrotransposon. Since TE mobilization can be controlled by specific environmental stimulation and as the result can generate novel genetic variations, it has been suggested that controlled mobilization of TEs can be a plausible method for crop breeding. To achieve this goal, series of sequencing techniques have been recently established to identify TEs that are active in mobility. These methods target and detect extrachromosomal DNAs (ecDNAs), which are final products of integration. The newly identified TEs by these methods exhibit strong transpositional activity which can generate novel genetic diversity and provide useful breeding resources.

CONCLUSIONS

In this mini review, we summarize and introduce ALE-seq, mobilome-seq, and VLP DNA-seq techniques employed to detect active TEs in plants.

Gene capture by transposable elements leads to epigenetic conflict in maize

Transposable elements (TEs) regularly capture fragments of genes. When the host silences these TEs, siRNAs homologous to the captured regions may also target the genes. This epigenetic crosstalk establishes an intragenomic conflict: silencing the TEs has the cost of silencing the genes. If genes are important, however, natural selection may maintain function by moderating the silencing response, which may also advantage the TEs. In this study, we examined this model by focusing on Helitrons, Pack-MULEs, and Sirevirus LTR retrotransposons in the maize genome. We documented 1263 TEs containing exon fragments from 1629 donor genes. Consistent with epigenetic conflict, donor genes mapped more siRNAs and were more methylated than genes with no evidence of capture. However, these patterns differed between syntelog versus translocated donor genes. Syntelogs appeared to maintain function, as measured by gene expression, consistent with moderation of silencing for functionally important genes. Epigenetic marks did not spread beyond their captured regions and 24nt crosstalk siRNAs were linked with CHH methylation. Translocated genes, in contrast, bore the signature of silencing. They were highly methylated and less expressed, but also overrepresented among donor genes and located away from chromosomal arms, which suggests a link between capture and gene movement. Splitting genes into potential functional categories based on evolutionary constraint supported the synteny-based findings. TE families captured genes in different ways, but the evidence for their advantage was generally less obvious; nevertheless, TEs with captured fragments were older, mapped fewer siRNAs, and were slightly less methylated than TEs without captured fragments. Collectively, our results argue that TE capture triggers an intragenomic conflict that may not affect the function of important genes but may lead to the pseudogenization of less-constrained genes.

Time to Wake Up: Epigenetic and Small-RNA-Mediated Regulation during Seed Germination

Plants make decisions throughout their lifetime based on complex networks. Phase transitions during seed growth are not an exception. From embryo development through seedling growth, several molecular pathways control genome stability, environmental signal transduction and the transcriptional landscape. Particularly, epigenetic modifications and small non-coding RNAs (sRNAs) have been extensively studied as significant handlers of these processes in plants. Here, we review key epigenetic (histone modifications and methylation patterns) and sRNA-mediated regulatory networks involved in the progression from seed maturation to germination, their relationship with seed traits and crosstalk with environmental inputs.

RJunBase: a database of RNA splice junctions in human normal and cancerous tissues

Splicing is an essential step of RNA processing for multi-exon genes, in which introns are removed from a precursor RNA, thereby producing mature RNAs containing splice junctions. Here, we develope the RJunBase (www.RJunBase.org), a web-accessible database of three types of RNA splice junctions (linear, back-splice, and fusion junctions) that are derived from RNA-seq data of non-cancerous and cancerous tissues. The RJunBase aims to integrate and characterize all RNA splice junctions of both healthy or pathological human cells and tissues. This new database facilitates the visualization of the gene-level splicing pattern and the junction-level expression profile, as well as the demonstration of unannotated and tumor-specific junctions. The first release of RJunBase contains 682 017 linear junctions, 225 949 back-splice junctions and 34 733 fusion junctions across 18 084 non-cancerous and 11 540 cancerous samples. RJunBase can aid researchers in discovering new splicing-associated targets and provide insights into the identification and assessment of potential neoepitopes for cancer treatment.

Nanopore RNA Sequencing Revealed Long Non-Coding and LTR Retrotransposon-Related RNAs Expressed at Early Stages of Triticale SEED Development

The intergenic space of plant genomes encodes many functionally important yet unexplored RNAs. The genomic loci encoding these RNAs are often considered “junk”, DNA as they are frequently associated with repeat-rich regions of the genome. The latter makes the annotations of these loci and the assembly of the corresponding transcripts using short RNAseq reads particularly challenging. Here, using long-read Nanopore direct RNA sequencing, we aimed to identify these “junk” RNA molecules, including long non-coding RNAs (lncRNAs) and transposon-derived transcripts expressed during early stages (10 days post anthesis) of seed development of triticale (AABBRR, 2n = 6x = 42), an interspecific hybrid between wheat and rye. Altogether, we found 796 lncRNAs and 20 LTR retrotransposon-related transcripts (RTE-RNAs) expressed at this stage, with most of them being previously unannotated and located in the intergenic as well as intronic regions. Sequence analysis of the lncRNAs provide evidence for the frequent exonization of Class I (retrotransposons) and class II (DNA transposons) transposon sequences and suggest direct influence of “junk” DNA on the structure and origin of lncRNAs. We show that the expression patterns of lncRNAs and RTE-related transcripts have high stage specificity. In turn, almost half of the lncRNAs located in Genomes A and B have the highest expression levels at 10–30 days post anthesis in wheat. Detailed analysis of the protein-coding potential of the RTE-RNAs showed that 75% of them carry open reading frames (ORFs) for a diverse set of GAG proteins, the main component of virus-like particles of LTR retrotransposons. We further experimentally demonstrated that some RTE-RNAs originate from autonomous LTR retrotransposons with ongoing transposition activity during early stages of triticale seed development. Overall, our results provide a framework for further exploration of the newly discovered lncRNAs and RTE-RNAs in functional and genome-wide association studies in triticale and wheat. Our study also demonstrates that Nanopore direct RNA sequencing is an indispensable tool for the elucidation of lncRNA and retrotransposon transcripts.

Genomic and Transcriptomic Survey Provides New Insight into the Organization and Transposition Activity of Highly Expressed LTR Retrotransposons of Sunflower (Helianthus annuus L.)

LTR retrotransposons (RTEs) play a crucial role in plant genome evolution and adaptation. Although RTEs are generally silenced in somatic plant tissues under non-stressed conditions, some expressed RTEs (exRTEs) escape genome defense mechanisms. As our understanding of exRTE organization in plants is rudimentary, we systematically surveyed the genomic and transcriptomic organization and mobilome (transposition) activity of sunflower (Helianthus annuus L.) exRTEs. We identified 44 transcribed RTEs in the sunflower genome and demonstrated their distinct genomic features: more recent insertion time, longer open reading frame (ORF) length, and smaller distance to neighboring genes. We showed that GAG-encoding ORFs are present at significantly higher frequencies in exRTEs, compared with non-expressed RTEs. Most exRTEs exhibit variation in copy number among sunflower cultivars and one exRTE Gagarin produces extrachromosomal circular DNA in seedling, demonstrating recent and ongoing transposition activity. Nanopore direct RNA sequencing of full-length RTE RNA revealed complex patterns of alternative splicing in RTE RNAs, resulting in isoforms that carry ORFs for distinct RTE proteins. Together, our study demonstrates that tens of expressed sunflower RTEs with specific genomic organization shape the hidden layer of the transcriptome, pointing to the evolution of specific strategies that circumvent existing genome defense mechanisms.