Additional ORFs in Plant LTR-Retrotransposons

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.

Genomic object detection: An improved approach for transposable elements detection and classification using convolutional neural networks

Analysis of eukaryotic genomes requires the detection and classification of transposable elements (TEs), a crucial but complex and time-consuming task. To improve the performance of tools that accomplish these tasks, Machine Learning approaches (ML) that leverage computer resources, such as GPUs (Graphical Processing Unit) and multiple CPU (Central Processing Unit) cores, have been adopted. However, until now, the use of ML techniques has mostly been limited to classification of TEs. Herein, a detection-classification strategy (named YORO) based on convolutional neural networks is adapted from computer vision (YOLO) to genomics. This approach enables the detection of genomic objects through the prediction of the position, length, and classification in large DNA sequences such as fully sequenced genomes. As a proof of concept, the internal protein-coding domains of LTR-retrotransposons are used to train the proposed neural network. Precision, recall, accuracy, F1-score, execution times and time ratios, as well as several graphical representations were used as metrics to measure performance. These promising results open the door for a new generation of Deep Learning tools for genomics. YORO architecture is available at https://github.com/simonorozcoarias/YORO.

The CERV protein of Cer1, a C. elegans LTR retrotransposon, is required for nuclear export of viral genomic RNA and can form giant nuclear rods

Retroviruses and closely related LTR retrotransposons export full-length, unspliced genomic RNA (gRNA) for packaging into virions and to serve as the mRNA encoding GAG and POL polyproteins. Because gRNA often includes splice acceptor and donor sequences used to splice viral mRNAs, retroelements must overcome host mechanisms that retain intron-containing RNAs in the nucleus. Here we examine gRNA expression in Cer1, an LTR retrotransposon in C. elegans which somehow avoids silencing and is highly expressed in germ cells. Newly exported Cer1 gRNA associates rapidly with the Cer1 GAG protein, which has structural similarity with retroviral GAG proteins. gRNA export requires CERV (C. elegans regulator of viral expression), a novel protein encoded by a spliced Cer1 mRNA. CERV phosphorylation at S214 is essential for gRNA export, and phosphorylated CERV colocalizes with nuclear gRNA at presumptive sites of transcription. By electron microscopy, tagged CERV proteins surround clusters of distinct, linear fibrils that likely represent gRNA molecules. Single fibrils, or groups of aligned fibrils, also localize near nuclear pores. During the C. elegans self-fertile period, when hermaphrodites fertilize oocytes with their own sperm, CERV concentrates in two nuclear foci that are coincident with gRNA. However, as hermaphrodites cease self-fertilization, and can only produce cross-progeny, CERV undergoes a remarkable transition to form giant nuclear rods or cylinders that can be up to 5 microns in length. We propose a novel mechanism of rod formation, in which stage-specific changes in the nucleolus induce CERV to localize to the nucleolar periphery in flattened streaks of protein and gRNA; these streaks then roll up into cylinders. The rods are a widespread feature of Cer1 in wild strains of C. elegans, but their function is not known and might be limited to cross-progeny. We speculate that the adaptive strategy Cer1 uses for the identical self-progeny of a host hermaphrodite might differ for heterozygous cross-progeny sired by males. For example, mating introduces male chromosomes which can have different, or no, Cer1 elements.

The discovery of multiple active mys-related LTR-retroelements within the Neotominae subfamily of cricetid rodents

Retrotransposon families in the rodent family Cricetidae have been understudied in contrast to Muridae, both taxa classified within the superfamily Muroidea. Therefore, we carried out a study to advance our knowledge of the unique mys LTR-retroelement identified in Peromyscus leucopus, by incorporating intra-ORF PCR, quantitative dot blots, DNA and protein library screens, the generation of molecular phylogenies, and analyses of orthologous LTR-retroelement loci. These analyses led to the discovery of three additional related families of LTR-retroelements, which include a 2900 bp full-length element of mys-related sequences (mysRS), an 8000 bp element containing the mys ORF1 sequence (mORF1) with ERV-related sequences downstream in the reverse orientation, as well as an 1800 bp element primarily consisting of mys ORF2 (mORF2) related sequences flanked by LTRs. Our data revealed only a few full-length mys elements among genera of the Neotominae subfamily of cricetid rodents, most existing as partial copies. The mysRS and mORF1 elements are also limited to the genomes of the Neotominae subfamily, whereas mORF2 appears to be restricted to the Peromyscus genus. Molecular phylogenies demonstrating concerted evolution along with an assessment of orthologous loci in Peromyscus for the presence or absence of elements are consistent with activity of these novel LTR-retroelement families within this genus. Together with known activity of various families of non-LTR retroelements in Peromyscus species, we propose that retrotransposons have been continually contributing to the dynamics of Peromyscus genomes promoting genomic diversity and may be correlated with the evolution of more than 50 identified Peromyscus species.

Miniature Inverted-Repeat Transposable Elements: Small DNA Transposons That Have Contributed to Plant MICRORNA Gene Evolution

Angiosperms form the largest phylum within the Plantae kingdom and show remarkable genetic variation due to the considerable difference in the nuclear genome size of each species. Transposable elements (TEs), mobile DNA sequences that can amplify and change their chromosome position, account for much of the difference in nuclear genome size between individual angiosperm species. Considering the dramatic consequences of TE movement, including the complete loss of gene function, it is unsurprising that the angiosperms have developed elegant molecular strategies to control TE amplification and movement. Specifically, the RNA-directed DNA methylation (RdDM) pathway, directed by the repeat-associated small-interfering RNA (rasiRNA) class of small regulatory RNA, forms the primary line of defense to control TE activity in the angiosperms. However, the miniature inverted-repeat transposable element (MITE) species of TE has at times avoided the repressive effects imposed by the rasiRNA-directed RdDM pathway. MITE proliferation in angiosperm nuclear genomes is due to their preference to transpose within gene-rich regions, a pattern of transposition that has enabled MITEs to gain further transcriptional activity. The sequence-based properties of a MITE results in the synthesis of a noncoding RNA (ncRNA), which, after transcription, folds to form a structure that closely resembles those of the precursor transcripts of the microRNA (miRNA) class of small regulatory RNA. This shared folding structure results in a MITE-derived miRNA being processed from the MITE-transcribed ncRNA, and post-maturation, the MITE-derived miRNA can be used by the core protein machinery of the miRNA pathway to regulate the expression of protein-coding genes that harbor homologous MITE insertions. Here, we outline the considerable contribution that the MITE species of TE have made to expanding the miRNA repertoire of the angiosperms.

The role of LTR retrotransposons in plant genetic engineering: how to control their transposition in the genome

We briefly discuss that the similarity of LTR retrotransposons to retroviruses is a great opportunity for the development of a genetic engineering tool that exploits intragenic elements in the plant genome for plant genetic improvement. Long terminal repeat (LTR) retrotransposons are very similar to retroviruses but do not have the property of being infectious. While spreading between its host cells, a retrovirus inserts a DNA copy of its genome into the cells. The ability of retroviruses to cause infection with genome integration allows genes to be delivered to cells and tissues. Retrovirus vectors are, however, only specific to animals and insects, and, thus, are not relevant to plant genetic engineering. However, the similarity of LTR retrotransposons to retroviruses is an opportunity to explore the former as a tool for genetic engineering. Although recent long-read sequencing technologies have advanced the knowledge about transposable elements (TEs), the integration of TEs is still unable either to control them or to direct them to specific genomic locations. The use of existing intragenic elements to achieve the desired genome composition is better than using artificial constructs like vectors, but it is not yet clear how to control the process. Moreover, most LTR retrotransposons are inactive and unable to produce complete proteins. They are also highly mutable. In addition, it is impossible to find a full active copy of a LTR retrotransposon out of thousands of its own copies. Theoretically, if these elements were directly controlled and turned on or off using certain epigenetic mechanisms (inducing by stress or infection), LTR retrotransposons could be a great opportunity to develop a genetic engineering tool using intragenic elements in the plant genome. In this review, the recent developments in uncovering the nature of LTR retrotransposons and the possibility of using these intragenic elements as a tool for plant genetic engineering are briefly discussed.

Identification, characterization, and expression profiling of the putative U-box E3 ubiquitin ligase gene family in Sorghum bicolor

The U-box family is one of the main E3 ubiquitin ligase families in plants. The U-box family has been characterized in several species. However, genome-wide gene identification and expression profiling of the U-box family in response to abiotic stress in Sorghum bicolor remain unclear. In this study, we broadly identified 68 U-box genes in the sorghum genome, including 2 CHIP genes, and 1 typical UFD2 (Ub fusion degradation 2) gene. The U-box gene family was divided into eight subclasses based on homology and conserved domain characteristics. Evolutionary analysis identified 14, 66, and 82 U-box collinear gene pairs in sorghum compared with arabidopsis, rice, and maize, respectively, and a unique tandem repeat pair (SbPUB26/SbPUB27) is present in the sorghum genome. Gene Ontology (GO) enrichment analysis showed that U-box proteins were mainly related to ubiquitination and modification, and various stress responses. Comprehensive analysis of promoters, expression profiling, and gene co-regulation networks also revealed that many sorghum U-box genes may be correlated with multiple stress responses. In summary, our results showed that sorghum contains 68 U-box genes, which may be involved in multiple abiotic stress responses. The findings will support future gene functional studies related to ubiquitination in sorghum.

Recent Acquisition of Functional m6A RNA Demethylase Domain in Orchid Ty3/Gypsy Elements

Long terminal repeats (LTR) retrotransposons are transposable elements (TEs) representing major components of most plant genomes. The fixation of additional conserved protein domains in their genomes is considered a rare event in the course of their evolution. Such changes can bring novel functions and increase their fitness by playing a role in the regulation of their replicative cycle or by affecting their integration landscape so that the detection of new domains can in turn reveal important aspects of host-TE interactions. We have mined angiosperm genomes for the presence of additional domains in LTR retrotransposons. We report a lineage of large (25 kbp) Gypsy-type elements in the genomes of Phalaenopsis orchids that contain an additional open reading frame containing a 2-ODD domain with close similarity to those responsible for m6A RNA demethylase activity in AlkB proteins. By performing in vitro assays, we demonstrate the RNA binding capability and the demethylase activity of the Gypsy-encoded AlkB protein, suggesting it could be functional against cognate TE mRNA or any cellular RNA in planta. In line with recent literature, we propose that the fixation of an RNA demethylase in this lineage of LTR retrotransposons may reflect an important role for epitranscriptomic control in host surveillance against TEs.

The structural diversity of CACTA transposons in genomes of Chenopodium (Amaranthaceae, Caryophyllales) species: specific traits and comparison with the similar elements of angiosperms

Background

CACTA transposable elements (TEs) comprise one of the most abundant superfamilies of Class 2 (cut-and-paste) transposons. Over recent decades, CACTA elements were widely identified in species from the plant, fungi, and animal kingdoms, but sufficiently studied in the genomes of only a few model species although non-model genomes can bring additional and valuable information. It primarily concerned the genomes of species belonging to clades in the base of large taxonomic groups whose genomes, to a certain extent, can preserve relict and/or possesses specific traits. Thus, we sought to investigate the genomes of Chenopodium (Amaranthaceae, Caryophyllales) species to unravel the structural variability of CACTA elements. Caryophyllales is a separate branch of Angiosperms and until recently the diversity of CACTA elements in this clade was unknown.

Results

Application of the short-read genome assembly algorithm followed by analysis of detected complete CACTA elements allowed for the determination of their structural diversity in the genomes of 22 Chenopodium album aggregate species. This approach yielded knowledge regarding: (i) the coexistence of two CACTA transposons subtypes in single genome; (ii) gaining of additional protein conserved domains within the coding sequence; (iii) the presence of captured gene fragments, including key genes for flower development; and (iv)) identification of captured satDNA arrays. Wide comparative database analysis revealed that identified events are scattered through Angiosperms in different proportions.

Conclusions

Our study demonstrated that while preserving the basic element structure a wide range of coding and non-coding additions to CACTA transposons occur in the genomes of C. album aggregate species. Ability to relocate additions inside genome in combination with the proposed novel functional features of structural-different CACTA elements can impact evolutionary trajectory of the host genome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00265-3.

Genomic structural variation in tomato and its role in plant immunity

It is well known that large genomic variations can greatly impact the phenotype of an organism. Structural Variants (SVs) encompass any genomic variation larger than 30 base pairs, and include changes caused by deletions, inversions, duplications, transversions, and other genome modifications. Due to their size and complex nature, until recently, it has been difficult to truly capture these variations. Recent advances in sequencing technology and computational analyses now permit more extensive studies of SVs in plant genomes. In tomato, advances in sequencing technology have allowed researchers to sequence hundreds of genomes from tomatoes, and tomato relatives. These studies have identified SVs related to fruit size and flavor, as well as plant disease response, resistance/susceptibility, and the ability of plants to detect pathogens (immunity). In this review, we discuss the implications for genomic structural variation in plants with a focus on its role in tomato immunity. We also discuss how advances in sequencing technology have led to new discoveries of SVs in more complex genomes, the current evidence for the role of SVs in biotic and abiotic stress responses, and the outlook for genetic modification of SVs to advance plant breeding objectives.

Comparative Transcriptome Sequencing of Taro Corm Development With a Focus on the Starch and Sucrose Metabolism Pathway

Taro (Colocasia esculenta) is an important tuber crop and staple food. Taro corms have higher nutritional value and starch contents as compared to most of the other root/tuber crops. However, the growth and development of the taro rhizome have not been critically examined in terms of transcriptomic signatures in general or specific to carbohydrates (starch and sucrose) accumulation. In current study, we have conducted a comprehensive survey of transcripts in taro corms aged 1, 2, 3, 4, 5, and 8 months. In this context, we have employed a whole transcriptome sequencing approach for identification of mRNAs, CircRNAs, and miRNAs in corms and performed functional enrichment analysis of the screened differentially expressed RNAs. A total of 11,203 mRNAs, 245 CircRNAs, and 299 miRNAs were obtained from six developmental stages. The mRNAs included 139 DEGs associated with 24 important enzymes of starch and sucrose metabolism. The expression of genes encoding key enzymes of starch and sucrose metabolism pathway (GBSS, AGPase, UGPase, SP, SSS, βFRUCT and SuSy) demonstrated significant variations at the stage of 4 months (S4). A total of 191 CircRNAs were differentially expressed between the studied comparisons of growth stages and 99 of these were associated with those miRNA (or target genes) that were enriched in starch and sucrose metabolism pathway. We also identified 205 miRNAs including 46 miRNAs targeting DEGs enriched in starch and sucrose biosynthesis pathway. The results of current study provide valuable resources for future exploration of the molecular mechanisms involved in the starch properties of Taro.

An 82 bp tandem repeat family typical of 3' non-coding end of Gypsy/TAT LTR retrotransposons is conserved in Coffea spp. pericentromeres

Coffea spp. chromosomes are very small and accumulate a variety of repetitive DNA families around the centromeres. However, the proximal regions of Coffea chromosomes remain poorly understood, especially regarding the nature and organisation of the sequences. Taking advantage of the genome sequences of C. arabica (2n = 44), C. canephora, and C. eugenioides (C. arabica progenitors with 2n = 22) and good coverage genome sequencing of dozens of other wild Coffea spp., repetitive DNA sequences were identified, and the genomes were compared to decipher particularities of pericentromeric structures. The searches revealed a short tandem repeat (82 bp length) typical of Gypsy/TAT LTR retrotransposons, named Coffea_sat11. This repeat organises clusters with fragments of other transposable elements, comprising regions of non-coding RNA production. Cytogenomic analyses showed that Coffea_sat11 extends from the pericentromeres towards the middle of the chromosomal arms. This arrangement was observed in the allotetraploid C. arabica chromosomes, as well as in its progenitors. This study improves our understanding of the role of the Gypsy/TAT LTR retrotransposon lineage in the organisation of Coffea pericentromeres, as well as the conservation of Coffea_sat11 within the genus. The relationships between fragments of other transposable elements and the functional aspects of these sequences on the pericentromere chromatin were also evaluated. Highlights: A scattered short tandem repeat, typical of Gypsy/TAT LTR retrotransposons, associated with several fragments of other transposable elements, accumulates in the pericentromeres of Coffea chromosomes. This arrangement is preserved in all clades of the genus and appears to have a strong regulatory role in the organisation of chromatin around centromeres.

Genome-wide characterization of LTR retrotransposons in the non-model deep-sea annelid Lamellibrachia luymesi

Background

Long Terminal Repeat retrotransposons (LTR retrotransposons) are mobile genetic elements composed of a few genes between terminal repeats and, in some cases, can comprise over half of a genome’s content. Available data on LTR retrotransposons have facilitated comparative studies and provided insight on genome evolution. However, data are biased to model systems and marine organisms, including annelids, have been underrepresented in transposable elements studies. Here, we focus on genome of Lamellibrachia luymesi, a vestimentiferan tubeworm from deep-sea hydrocarbon seeps, to gain knowledge of LTR retrotransposons in a deep-sea annelid.

Results

We characterized LTR retrotransposons present in the genome of L. luymesi using bioinformatic approaches and found that intact LTR retrotransposons makes up about 0.1% of L. luymesi genome. Previous characterization of the genome has shown that this tubeworm hosts several known LTR-retrotransposons. Here we describe and classify LTR retrotransposons in L. luymesi as within the Gypsy, Copia and Bel-pao superfamilies. Although, many elements fell within already recognized families (e.g., Mag, CSRN1), others formed clades distinct from previously recognized families within these superfamilies. However, approximately 19% (41) of recovered elements could not be classified. Gypsy elements were the most abundant while only 2 Copia and 2 Bel-pao elements were present. In addition, analysis of insertion times indicated that several LTR-retrotransposons were recently transposed into the genome of L. luymesi, these elements had identical LTR’s raising possibility of recent or ongoing retrotransposon activity.

Conclusions

Our analysis contributes to knowledge on diversity of LTR-retrotransposons in marine settings and also serves as an important step to assist our understanding of the potential role of retroelements in marine organisms. We find that many LTR retrotransposons, which have been inserted in the last few million years, are similar to those found in terrestrial model species. However, several new groups of LTR retrotransposons were discovered suggesting that the representation of LTR retrotransposons may be different in marine settings. Further study would improve understanding of the diversity of retrotransposons across animal groups and environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07749-1.