The diversity of retrotransposons and the properties of their reverse transcriptases

Diversity and evolution of transposable elements in the plant-parasitic nematodes

BACKGROUND

Transposable elements (TEs) are mobile DNA sequences that propagate within genomes, occupying a significant portion of eukaryotic genomes and serving as a source of genetic variation and innovation. TEs can impact genome dynamics through their repetitive nature and mobility. Nematodes are incredibly versatile organisms, capable of thriving in a wide range of environments. The plant-parasitic nematodes are able to infect nearly all vascular plants, leading to significant crop losses and management expenses worldwide. It is worth noting that plant parasitism has evolved independently at least three times within this nematode group. Furthermore, the genome size of plant-parasitic nematodes can vary substantially, spanning from 41.5 Mbp to 235 Mbp. To investigate genome size variation and evolution in plant-parasitic nematodes, TE composition, diversity, and evolution were analysed in 26 plant-parasitic nematodes from 9 distinct genera in Clade IV.

RESULTS

Interestingly, despite certain species lacking specific types of DNA transposons or retrotransposon superfamilies, they still exhibit a diverse range of TE content. Identification of species-specific TE repertoire in nematode genomes provides a deeper understanding of genome evolution in plant-parasitic nematodes. An intriguing observation is that plant-parasitic nematodes possess extensive DNA transposons and retrotransposon insertions, including recent sightings of LTR/Gypsy and LTR/Pao superfamilies. Among them, the Gypsy superfamilies were found to encode Aspartic proteases in the plant-parasitic nematodes.

CONCLUSIONS

The study of the transposable element (TE) composition in plant-parasitic nematodes has yielded insightful discoveries. The findings revealed that certain species exhibit lineage-specific variations in their TE makeup. Discovering the species-specific TE repertoire in nematode genomes is a crucial element in understanding the evolution of genomes in plant-parasitic nematodes. It allows us to gain a deeper insight into the intricate workings of these organisms and their genetic makeup. With this knowledge, we are gaining a fundamental piece in the puzzle of understanding the evolution of these parasites. Moreover, recent transpositions have led to the acquisition of new TE superfamilies, especially Gypsy and Pao retrotransposons, further expanding the diversity of TEs in these nematodes. Significantly, the widely distributed Gypsy superfamily possesses proteases that are exclusively associated with parasitism during nematode-host interactions. These discoveries provide a deeper understanding of the TE landscape within plant-parasitic nematodes.

Transposable Element Expression Profiles in Premalignant Pigment Cell Lesions and Melanoma of Xiphophorus

Transposable elements (TEs) are characterized by their ability to change their genomic position. Through insertion or recombination leading to deletions and other chromosomal aberrations, they can cause genetic instability. The extent to which they thereby exert regulatory influence on cellular functions is unclear. To better characterize TEs in processes such as carcinogenesis, we used the well-established Xiphophorus melanoma model. By transcriptome sequencing, we show that an increasing total number in transposons correlates with progression of malignancy in melanoma samples from Xiphophorus interspecific hybrids. Further, by comparing the presence of TEs in the parental genomes of Xiphophorus maculatus and Xiphophorus hellerii, we could show that even in closely related species, genomic location and spectrum of TEs are considerably different.

Immune system modulation & virus transmission during parasitism identified by multi-species transcriptomics of a declining insect biocontrol system

BACKGROUND

The Argentine stem weevil (ASW, Listronotus bonariensis) is a significant pasture pest in Aotearoa New Zealand, primarily controlled by the parasitoid biocontrol agent Microctonus hyperodae. Despite providing effective control of ASW soon after release, M. hyperodae parasitism rates have since declined significantly, with ASW hypothesised to have evolved resistance to its biocontrol agent. While the parasitism arsenal of M. hyperodae has previously been investigated, revealing many venom components and an exogenous novel DNA virus Microctonus hyperodae filamentous virus (MhFV), the effects of said arsenal on gene expression in ASW during parasitism have not been examined. In this study, we performed a multi-species transcriptomic analysis to investigate the biology of ASW parasitism by M. hyperodae, as well as the decline in efficacy of this biocontrol system.

RESULTS

The transcriptomic response of ASW to parasitism by M. hyperodae involves modulation of the weevil's innate immune system, flight muscle components, and lipid and glucose metabolism. The multispecies approach also revealed continued expression of venom components in parasitised ASW, as well as the transmission of MhFV to weevils during parasitism and some interrupted parasitism attempts. Transcriptomics did not detect a clear indication of parasitoid avoidance or other mechanisms to explain biocontrol decline.

CONCLUSIONS

This study has expanded our understanding of interactions between M. hyperodae and ASW in a biocontrol system of critical importance to Aotearoa-New Zealand's agricultural economy. Transmission of MhFV to ASW during successful and interrupted parasitism attempts may link to a premature mortality phenomenon in ASW, hypothesised to be a result of a toxin-antitoxin system. Further research into MhFV and its potential role in ASW premature mortality is required to explore whether manipulation of this viral infection has the potential to increase biocontrol efficacy in future.

Evolving Together: Cassandra Retrotransposons Gradually Mirror Promoter Mutations of the 5S rRNA Genes

The 5S rRNA genes are among the most conserved nucleotide sequences across all species. Similar to the 5S preservation we observe the occurrence of 5S-related nonautonomous retrotransposons, so-called Cassandras. Cassandras harbor highly conserved 5S rDNA-related sequences within their long terminal repeats, advantageously providing them with the 5S internal promoter. However, the dynamics of Cassandra retrotransposon evolution in the context of 5S rRNA gene sequence information and structural arrangement are still unclear, especially: (1) do we observe repeated or gradual domestication of the highly conserved 5S promoter by Cassandras and (2) do changes in 5S organization such as in the linked 35S-5S rDNA arrangements impact Cassandra evolution? Here, we show evidence for gradual co-evolution of Cassandra sequences with their corresponding 5S rDNAs. To follow the impact of 5S rDNA variability on Cassandra TEs, we investigate the Asteraceae family where highly variable 5S rDNAs, including 5S promoter shifts and both linked and separated 35S-5S rDNA arrangements have been reported. Cassandras within the Asteraceae mirror 5S rDNA promoter mutations of their host genome, likely as an adaptation to the host's specific 5S transcription factors and hence compensating for evolutionary changes in the 5S rDNA sequence. Changes in the 5S rDNA sequence and in Cassandras seem uncorrelated with linked/separated rDNA arrangements. We place all these observations into the context of angiosperm 5S rDNA-Cassandra evolution, discuss Cassandra's origin hypotheses (single or multiple) and Cassandra's possible impact on rDNA and plant genome organization, giving new insights into the interplay of ribosomal genes and transposable elements.

Large Deletions, Cleavage of the Telomeric Repeat Sequence, and Reverse Transcriptase-Mediated DNA Damage Response Associated with Long Interspersed Element-1 ORF2p Enzymatic Activities

L1 elements can cause DNA damage and genomic variation via retrotransposition and the generation of endonuclease-dependent DNA breaks. These processes require L1 ORF2p protein that contains an endonuclease domain, which cuts genomic DNA, and a reverse transcriptase domain, which synthesizes cDNA. The complete impact of L1 enzymatic activities on genome stability and cellular function remains understudied, and the spectrum of L1-induced mutations, other than L1 insertions, is mostly unknown. Using an inducible system, we demonstrate that an ORF2p containing functional reverse transcriptase is sufficient to elicit DNA damage response even in the absence of the functional endonuclease. Using a TK/Neo reporter system that captures misrepaired DNA breaks, we demonstrate that L1 expression results in large genomic deletions that lack any signatures of L1 involvement. Using an in vitro cleavage assay, we demonstrate that L1 endonuclease efficiently cuts telomeric repeat sequences. These findings support that L1 could be an unrecognized source of disease-promoting genomic deletions, telomere dysfunction, and an underappreciated source of chronic RT-mediated DNA damage response in mammalian cells. Our findings expand the spectrum of biological processes that can be triggered by functional and nonfunctional L1s, which have impactful evolutionary- and health-relevant consequences.

MegaLTR: a web server and standalone pipeline for detecting and annotating LTR-retrotransposons in plant genomes

LTR-retrotransposons (LTR-RTs) are a class of RNA-replicating transposon elements (TEs) that can alter genome structure and function by moving positions, repositioning genes, shifting exons, and causing chromosomal rearrangements. LTR-RTs are widespread in many plant genomes and constitute a significant portion of the genome. Their movement and activity in eukaryotic genomes can provide insight into genome evolution and gene function, especially when LTR-RTs are located near or within genes. Building the redundant and non-redundant LTR-RTs libraries and their annotations for species lacking this resource requires extensive bioinformatics pipelines and expensive computing power to analyze large amounts of genomic data. This increases the need for online services that provide computational resources with minimal overhead and maximum efficiency. Here, we present MegaLTR as a web server and standalone pipeline that detects intact LTR-RTs at the whole-genome level and integrates multiple tools for structure-based, homologybased, and de novo identification, classification, annotation, insertion time determination, and LTR-RT gene chimera analysis. MegaLTR also provides statistical analysis and visualization with multiple tools and can be used to accelerate plant species discovery and assist breeding programs in their efforts to improve genomic resources. We hope that the development of online services such as MegaLTR, which can analyze large amounts of genomic data, will become increasingly important for the automated detection and annotation of LTR-RT elements.

Origin and Evolution of Plant Long Terminal Repeat Retrotransposons with Additional Ribonuclease H

Retroviruses originated from long terminal repeat retrotransposons (LTR-RTs) through several structural adaptations. One such modification was the arrangement of an additional ribonuclease H (aRH) domain next to native RH, followed by degradation and subfunctionalization of the latter. We previously showed that this retrovirus-like structure independently evolved in Tat LTR-RTs in flowering plants, proposing its origin from sequential rearrangements of ancestral Tat structures identified in lycophytes and conifers. However, most nonflowering plant genome assemblies were not available at that time, therefore masking the history of aRH acquisition by Tat and challenging our hypothesis. Here, we revisited Tat's evolution scenario upon the aRH acquisition by covering most of the extant plant phyla. We show that Tat evolved and obtained aRH in an ancestor of land plants. Importantly, we found the retrovirus-like structure in clubmosses, hornworts, ferns, and gymnosperms, suggesting its ancient origin, broad propagation, and yet-to-be-understood benefit for the LTR-RTs' adaptation.

From parasites to partners: exploring the intricacies of host-transposon dynamics and coevolution

Transposable elements, often referred to as "jumping genes," have long been recognized as genomic parasites due to their ability to integrate and disrupt normal gene function and induce extensive genomic alterations, thereby compromising the host's fitness. To counteract this, the host has evolved a plethora of mechanisms to suppress the activity of the transposons. Recent research has unveiled the host-transposon relationships to be nuanced and complex phenomena, resulting in the coevolution of both entities. Transposition increases the mutational rate in the host genome, often triggering physiological pathways such as immune and stress responses. Current gene transfer technologies utilizing transposable elements have potential drawbacks, including off-target integration, induction of mutations, and modifications of cellular machinery, which makes an in-depth understanding of the host-transposon relationship imperative. This review highlights the dynamic interplay between the host and transposable elements, encompassing various factors and components of the cellular machinery. We provide a comprehensive discussion of the strategies employed by transposable elements for their propagation, as well as the mechanisms utilized by the host to mitigate their parasitic effects. Additionally, we present an overview of recent research identifying host proteins that act as facilitators or inhibitors of transposition. We further discuss the evolutionary outcomes resulting from the genetic interactions between the host and the transposable elements. Finally, we pose open questions in this field and suggest potential avenues for future research.

Reverse transcriptases prime DNA synthesis

The discovery of reverse transcriptases (RTs) challenged the central dogma by establishing that genetic information can also flow from RNA to DNA. Although they act as DNA polymerases, RTs are distantly related to replicases that also possess de novo primase activity. Here we identify that CRISPR associated RTs (CARTs) directly prime DNA synthesis on both RNA and DNA. We demonstrate that RT-dependent priming is utilized by some CRISPR-Cas complexes to synthesise new spacers and integrate these into CRISPR arrays. Expanding our analyses, we show that primer synthesis activity is conserved in representatives of other major RT classes, including group II intron RT, telomerase and retroviruses. Together, these findings establish a conserved innate ability of RTs to catalyse de novo DNA primer synthesis, independently of accessory domains or alternative priming mechanisms, which likely plays important roles in a wide variety of biological pathways.

A novel role of TRIM28 B box domain in L1 retrotransposition and ORF2p-mediated cDNA synthesis

The long interspersed element 1 (LINE-1 or L1) integration is affected by many cellular factors through various mechanisms. Some of these factors are required for L1 amplification, while others either suppress or enhance specific steps during L1 propagation. Previously, TRIM28 has been identified to suppress transposable elements, including L1 expression via its canonical role in chromatin remodeling. Here, we report that TRIM28 through its B box domain increases L1 retrotransposition and facilitates shorter cDNA and L1 insert generation in cultured cells. Consistent with the latter, we observe that tumor specific L1 inserts are shorter in endometrial, ovarian, and prostate tumors with higher TRIM28 mRNA expression than in those with lower TRIM28 expression. We determine that three amino acids in the B box domain that are involved in TRIM28 multimerization are critical for its effect on both L1 retrotransposition and cDNA synthesis. We provide evidence that B boxes from the other two members in the Class VI TRIM proteins, TRIM24 and TRIM33, also increase L1 retrotransposition. Our findings could lead to a better understanding of the host/L1 evolutionary arms race in the germline and their interplay during tumorigenesis.

Insights into mammalian TE diversity through the curation of 248 genome assemblies

We examined transposable element (TE) content of 248 placental mammal genome assemblies, the largest de novo TE curation effort in eukaryotes to date. We found that although mammals resemble one another in total TE content and diversity, they show substantial differences with regard to recent TE accumulation. This includes multiple recent expansion and quiescence events across the mammalian tree. Young TEs, particularly long interspersed elements, drive increases in genome size, whereas DNA transposons are associated with smaller genomes. Mammals tend to accumulate only a few types of TEs at any given time, with one TE type dominating. We also found association between dietary habit and the presence of DNA transposon invasions. These detailed annotations will serve as a benchmark for future comparative TE analyses among placental mammals.

PlantLTRdb: An interactive database for 195 plant species LTR-retrotransposons

LTR-retrotransposons (LTR-RTs) are a large group of transposable elements that replicate through an RNA intermediate and alter genome structure. The activities of LTR-RTs in plant genomes provide helpful information about genome evolution and gene function. LTR-RTs near or within genes can directly alter gene function. This work introduces PlantLTRdb, an intact LTR-RT database for 195 plant species. Using homology- and de novo structure-based methods, a total of 150.18 Gbp representing 3,079,469 pseudomolecules/scaffolds were analyzed to identify, characterize, annotate LTR-RTs, estimate insertion ages, detect LTR-RT-gene chimeras, and determine nearby genes. Accordingly, 520,194 intact LTR-RTs were discovered, including 29,462 autonomous and 490,732 nonautonomous LTR-RTs. The autonomous LTR-RTs included 10,286 Gypsy and 19,176 Copia, while the nonautonomous were divided into 224,906 Gypsy, 218,414 Copia, 1,768 BARE-2, 3,147 TR-GAG and 4,2497 unknown. Analysis of the identified LTR-RTs located within genes showed that a total of 36,236 LTR-RTs were LTR-RT-gene chimeras and 11,619 LTR-RTs were within pseudo-genes. In addition, 50,026 genes are within 1 kbp of LTR-RTs, and 250,587 had a distance of 1 to 10 kbp from LTR-RTs. PlantLTRdb allows researchers to search, visualize, BLAST and analyze plant LTR-RTs. PlantLTRdb can contribute to the understanding of structural variations, genome organization, functional genomics, and the development of LTR-RT target markers for molecular plant breeding. PlantLTRdb is available at https://bioinformatics.um6p.ma/PlantLTRdb.

An update on post-transcriptional regulation of retrotransposons

Retrotransposons, including LINE-1, Alu, SVA, and endogenous retroviruses, are one of the major constituents of human genomic repetitive sequences. Through the process of retrotransposition, some of them occasionally insert into new genomic locations by a copy-paste mechanism involving RNA intermediates. Irrespective of de novo genomic insertions, retrotransposon expression can lead to DNA double-strand breaks and stimulate cellular innate immunity through endogenous patterns. As a result, retrotransposons are tightly regulated by multi-layered regulatory processes to prevent the dangerous effects of their expression. In recent years, significant progress was made in revealing how retrotransposon biology intertwines with general post-transcriptional RNA metabolism. Here, I summarize current knowledge on the involvement of post-transcriptional factors in the biology of retrotransposons, focusing on LINE-1. I emphasize general RNA metabolisms such as methylation of adenine (m⁶ A), RNA 3'-end polyadenylation and uridylation, RNA decay and translation regulation. I discuss the effects of retrotransposon RNP sequestration in cytoplasmic bodies and autophagy. Finally, I summarize how innate immunity restricts retrotransposons and how retrotransposons make use of cellular enzymes, including the DNA repair machinery, to complete their replication cycles.

Alternative and aberrant splicing of human endogenous retroviruses in cancer. What about head and neck? —A mini review

Human endogenous retroviruses (HERVs) are transcribed in many cancer types, including head and neck cancer. Because of accumulating mutations at proviral loci over evolutionary time, HERVs are functionally defective and cannot complete their viral life cycle. Despite that, HERV transcripts, including full-length viral RNAs and viral RNAs spliced as expected at the conventional viral splice sites, can be detected in particular conditions, such as cancer. Interestingly, non-viral–related transcription, including aberrant, non-conventionally spliced RNAs, has been reported as well. The role of HERV transcription in cancer and its contribution to oncogenesis or progression are still debated. Nonetheless, HERVs may constitute a suitable cancer biomarker or a target for therapy. Thus, ongoing research aims both to clarify the basic mechanisms underlying HERV transcription in cancer and to exploit its potential toward clinical application. In this mini-review, we summarize the current knowledge, the most recent findings, and the future perspectives of research on HERV transcription and splicing, with particular focus on head and neck cancer.

Identification of an integrase-independent pathway of retrotransposition

Retroviruses and long terminal repeat retrotransposons rely on integrase (IN) to insert their complementary DNA (cDNA) into the genome of host cells. Nevertheless, in the absence of IN, retroelements can retain notable levels of insertion activity. We have characterized the IN-independent pathway of Tf1 and found that insertion sites had homology to the primers of reverse transcription, which form single-stranded DNAs at the termini of the cDNA. In the absence of IN activity, a similar bias was observed with HIV-1. Our studies showed that the Tf1 insertions result from single-strand annealing, a noncanonical form of homologous recombination mediated by Rad52. By expanding our analysis of insertions to include repeat sequences, we found most formed tandem elements by inserting at preexisting copies of a related transposable element. Unexpectedly, we found that wild-type Tf1 uses the IN-independent pathway as an alternative mode of insertion.

SubPhaser: a robust allopolyploid subgenome phasing method based on subgenome-specific k-mers

With advanced sequencing technology, dozens of complex polyploid plant genomes have been characterized. However, for many polyploid species, their diploid ancestors are unknown or extinct, making it impossible to unravel the subgenomes and genome evolution directly. We developed a novel subgenome-phasing algorithm, SubPhaser, specifically designed for a neoallopolyploid or a homoploid hybrid. SubPhaser first searches for the subgenome-specific sequence (k-mer), then assigns homoeologous chromosomes into subgenomes, and further provides tools to annotate and investigate specific sequences. SubPhaser works well on neoallopolyploids and homoploid hybrids containing subgenome-specific sequences like wheat, but fails on autopolyploids lacking subgenome-specific sequences like alfalfa, indicating that SubPhaser can phase neoallopolyploid/homoploid hybrids with high accuracy, sensitivity and performance. This highly accurate, highly sensitive, ancestral data free chromosome phasing algorithm, SubPhaser, offers significant application value for subgenome phasing in neoallopolyploids and homoploid hybrids, and for the subsequent exploration of genome evolution and related genetic/epigenetic mechanisms.

The Role of Transposable Elements of the Human Genome in Neuronal Function and Pathology

Transposable elements (TEs) have been extensively studied for decades. In recent years, the introduction of whole-genome and whole-transcriptome approaches, as well as single-cell resolution techniques, provided a breakthrough that uncovered TE involvement in host gene expression regulation underlying multiple normal and pathological processes. Of particular interest is increased TE activity in neuronal tissue, and specifically in the hippocampus, that was repeatedly demonstrated in multiple experiments. On the other hand, numerous neuropathologies are associated with TE dysregulation. Here, we provide a comprehensive review of literature about the role of TEs in neurons published over the last three decades. The first chapter of the present review describes known mechanisms of TE interaction with host genomes in general, with the focus on mammalian and human TEs; the second chapter provides examples of TE exaptation in normal neuronal tissue, including TE involvement in neuronal differentiation and plasticity; and the last chapter lists TE-related neuropathologies. We sought to provide specific molecular mechanisms of TE involvement in neuron-specific processes whenever possible; however, in many cases, only phenomenological reports were available. This underscores the importance of further studies in this area.

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.

A Missing Link between Retrotransposons and Retroviruses

The origin and deep evolution of retroviruses remain largely unclear. It has been proposed that retroviruses might have originated from a Ty3/Gypsy retrotransposon, but all known Ty3/Gypsy retrotransposons are only distantly related to retroviruses. Retroviruses and some plant Athila/Tat elements (within Ty3/Gypsy retrotransposons) independently evolved a dual RNase H domain and an env/env-like gene. Here, we reported the discovery of a novel lineage of retrotransposons, designated Odin retrotransposons, in the genomes of eight sea anemones (order Actinaria) within the Cnidaria phylum. Odin retrotransposons exhibited unique genome features, encoding a dual RNase H domain (like retroviruses) but no env gene (like most Ty3/Gypsy retrotransposons). Phylogenetic analyses based on reverse transcriptase showed that Odin retrotransposons formed a sister group to lokiretroviruses, and lokiretroviruses and Odin retrotransposons together were sister to canonical retroviruses. Moreover, phylogenetic analyses based on RNase H and integrase also supported the hypothesis that Odin retrotransposons were sisters to lokiretroviruses. Lokiretroviruses and canonical retroviruses did not form a monophyletic group, indicating that lokiretroviruses and canonical retroviruses might represent two distinct virus families. Taken together, the discovery of Odin retrotransposons narrowed down the evolutionary gaps between retrotransposons and canonical retroviruses and lokiretroviruses.

An endogenous retrovirus presumed to have been endogenized or relocated recently in a marsupial, the red-necked wallaby

An albino infant wallaby was born to a mother with the wild-type body color. PCR and sequencing analyses of TYR (encoding tyrosinase, which is essential for melanin biosynthesis) of this albino wallaby revealed a 7.1-kb-long DNA fragment inserted in the first exon. Because the fragment carried long terminal repeats, we assumed it to be a copy of an endogenous retrovirus, which we named walb. We cloned other walb copies residing in the genomes of this species and another wallaby species. The copies exhibited length variation, and the longest copy (>8.0 kb) contained open reading frames whose deduced amino acid sequences were well aligned with those of gag, pol, and env of retroviruses. It is not known through which of the following likely processes the walb copy was inserted into TYR: endogenization (infection of a germline cell by an exogenous virus), reinfection (infection by a virus produced from a previously endogenized provirus), or retrotransposition (intracellular relocation of a provirus). In any case, the insertion into TYR is considered to have been a recent event on an evolutionary timescale because albino mutant alleles generally do not persist for long because of their deleterious effects in wild circumstances. .

DARTS: An Algorithm for Domain-Associated Retrotransposon Search in Genome Assemblies

Retrotransposons comprise a substantial fraction of eukaryotic genomes, reaching the highest proportions in plants. Therefore, identification and annotation of retrotransposons is an important task in studying the regulation and evolution of plant genomes. The majority of computational tools for mining transposable elements (TEs) are designed for subsequent genome repeat masking, often leaving aside the element lineage classification and its protein domain composition. Additionally, studies focused on the diversity and evolution of a particular group of retrotransposons often require substantial customization efforts from researchers to adapt existing software to their needs. Here, we developed a computational pipeline to mine sequences of protein-coding retrotransposons based on the sequences of their conserved protein domains—DARTS (Domain-Associated Retrotransposon Search). Using the most abundant group of TEs in plants—long terminal repeat (LTR) retrotransposons (LTR-RTs)—we show that DARTS has radically higher sensitivity for LTR-RT identification compared to the widely accepted tool LTRharvest. DARTS can be easily customized for specific user needs. As a result, DARTS returns a set of structurally annotated nucleotide and amino acid sequences which can be readily used in subsequent comparative and phylogenetic analyses. DARTS may facilitate researchers interested in the discovery and detailed analysis of the diversity and evolution of retrotransposons, LTR-RTs, and other protein-coding TEs.

Identification and classification of reverse transcriptases in bacterial genomes and metagenomes

Reverse transcriptases (RTs) are found in different systems including group II introns, Diversity Generating Retroelements (DGRs), retrons, CRISPR-Cas systems, and Abortive Infection (Abi) systems in prokaryotes. Different classes of RTs can play different roles, such as template switching and mobility in group II introns, spacer acquisition in CRISPR-Cas systems, mutagenic retrohoming in DGRs, programmed cell suicide in Abi systems, and recently discovered phage defense in retrons. While some classes of RTs have been studied extensively, others remain to be characterized. There is a lack of computational tools for identifying and characterizing various classes of RTs. In this study, we built a tool (called myRT) for identification and classification of prokaryotic RTs. In addition, our tool provides information about the genomic neighborhood of each RT, providing potential functional clues. We applied our tool to predict RTs in all complete and draft bacterial genomes, and created a collection that can be used for exploration of putative RTs and their associated protein domains. Application of myRT to metagenomes showed that gut metagenomes encode proportionally more RTs related to DGRs, outnumbering retron-related RTs, as compared to the collection of reference genomes. MyRT is both available as a standalone software (https://github.com/mgtools/myRT) and also through a website (https://omics.informatics.indiana.edu/myRT/).

Long terminal repeats (LTR) and transcription factors regulate PHRE1 and PHRE2 activity in Moso bamboo under heat stress

BACKGROUND

LTR retrotransposons play a significant role in plant growth, genome evolution, and environmental stress response, but their regulatory response to heat stress remains unclear. We have investigated the activities of two LTR retrotransposons, PHRE1 and PHRE2, of moso bamboo (Phyllostachys edulis) in response to heat stress.

RESULTS

The differential overexpression of PHRE1 and PHRE2 with or without CaMV35s promoter showed enhanced expression under heat stress in transgenic plants. The transcriptional activity studies showed an increase in transposition activity and copy number among moso bamboo wild type and Arabidopsis transgenic plants under heat stress. Comparison of promoter activity in transgenic plants indicated that 5'LTR promoter activity was higher than CaMV35s promoter. Additionally, yeast one-hybrid (Y1H) system and in planta biomolecular fluorescence complementation (BiFC) assay revealed interactions of heat-dependent transcription factors (TFs) with 5'LTR sequence and direct interactions of TFs with pol and gag.

CONCLUSIONS

Our results conclude that the 5'LTR acts as a promoter and could regulate the LTR retrotransposons in moso bamboo under heat stress.

Prokaryotic reverse transcriptases: from retroelements to specialized defense systems

Reverse transcriptases (RTs) catalyze the polymerization of DNA from an RNA template. These enzymes were first discovered in RNA tumor viruses in 1970, but it was not until 1989 that they were found in prokaryotes as a key component of retrons. Apart from RTs encoded by the 'selfish' mobile retroelements known as group II introns, prokaryotic RTs are extraordinarily diverse, but their function has remained elusive. However, recent studies have revealed that different lineages of prokaryotic RTs, including retrons, those associated with CRISPR-Cas systems, Abi-like RTs and other yet uncharacterized RTs, are key components of different lines of defense against phages and other mobile genetic elements. Prokaryotic RTs participate in various antiviral strategies, including abortive infection (Abi), in which the infected cell is induced to commit suicide to protect the host population, adaptive immunity, in which a memory of previous infection is used to build an efficient defense, and other as yet unidentified mechanisms. These prokaryotic enzymes are attracting considerable attention, both for use in cutting-edge technologies, such as genome editing, and as an emerging research topic. In this review, we discuss what is known about prokaryotic RTs, and the exciting evidence for their domestication from retroelements to create specialized defense systems.

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

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

Comparative genomic and transcriptomic analyses of transposable elements in polychaetous annelids highlight LTR retrotransposon diversity and evolution

Background

With the expansion of high throughput sequencing, we now have access to a larger number of genome-wide studies analyzing the Transposable elements (TEs) composition in a wide variety of organisms. However, genomic analyses often remain too limited in number and diversity of species investigated to study in depth the dynamics and evolutionary success of the different types of TEs among metazoans. Therefore, we chose to investigate the use of transcriptomes to describe the diversity of TEs in phylogenetically related species by conducting the first comparative analysis of TEs in two groups of polychaetes and evaluate the diversity of TEs that might impact genomic evolution as a result of their mobility.

Results

We present a detailed analysis of TEs distribution in transcriptomes extracted from 15 polychaetes depending on the number of reads used during assembly, and also compare these results with additional TE scans on associated low-coverage genomes. We then characterized the clades defined by 1021 LTR-retrotransposon families identified in 26 species. Clade richness was highly dependent on the considered superfamily. Copia elements appear rare and are equally distributed in only three clades, GalEa, Hydra and CoMol. Among the eight BEL/Pao clades identified in annelids, two small clades within the Sailor lineage are new for science. We characterized 17 Gypsy clades of which only 4 are new; the C-clade largely dominates with a quarter of the families. Finally, all species also expressed for the majority two distinct transcripts encoding PIWI proteins, known to be involved in control of TEs mobilities.

Conclusions

This study shows that the use of transcriptomes assembled from 40 million reads was sufficient to access to the diversity and proportion of the transposable elements compared to those obtained by low coverage sequencing. Among LTR-retrotransposons Gypsy elements were unequivocally dominant but results suggest that the number of Gypsy clades, although high, may be more limited than previously thought in metazoans. For BEL/Pao elements, the organization of clades within the Sailor lineage appears more difficult to establish clearly. The Copia elements remain rare and result from the evolutionary consistent success of the same three clades.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00252-0.

Transposable Element Dynamics and Regulation during Zygotic Genome Activation in Mammalian Embryos and Embryonic Stem Cell Model Systems

Transposable elements (TEs) are mobile genetic sequences capable of duplicating and reintegrating at new regions within the genome. A growing body of evidence has demonstrated that these elements play important roles in host genome evolution, despite being traditionally viewed as parasitic elements. To prevent ectopic activation of TE transposition and transcription, they are epigenetically silenced in most somatic tissues. Intriguingly, a specific class of TEs-retrotransposons-is transiently expressed at discrete phases during mammalian development and has been linked to the establishment of totipotency during zygotic genome activation (ZGA). While mechanisms controlling TE regulation in somatic tissues have been extensively studied, the significance underlying the unique transcriptional reactivation of retrotransposons during ZGA is only beginning to be uncovered. In this review, we summarize the expression dynamics of key retrotransposons during ZGA, focusing on findings from in vivo totipotent embryos and in vitro totipotent-like embryonic stem cells (ESCs). We then dissect the functions of retrotransposons and discuss how their transcriptional activities are finetuned during early stages of mammalian development.

Low-bias ncRNA libraries using ordered two-template relay: Serial template jumping by a modified retroelement reverse transcriptase

Retrotransposons are noninfectious, mobile genetic elements that proliferate in host genomes via an RNA intermediate that is copied into DNA by a reverse transcriptase (RT) enzyme. RTs are important for biotechnological applications involving information capture from RNA since RNA is first converted into complementary DNA for detection or sequencing. Here, we biochemically characterized RTs from two retroelements and uncovered several activities that allowed us to design a streamlined, efficient workflow for determining the inventory of RNA sequences in processed RNA pools. The unique properties of nonretroviral RT activities obviate many technical issues associated with current methods of RNA sequence analysis, with wide applications in research, biotechnology, and diagnostics.

Selfish, non-long terminal repeat (non-LTR) retroelements and mobile group II introns encode reverse transcriptases (RTs) that can initiate DNA synthesis without substantial base pairing of primer and template. Biochemical characterization of these enzymes has been limited by recombinant expression challenges, hampering understanding of their properties and the possible exploitation of their properties for research and biotechnology. We investigated the activities of representative RTs using a modified non-LTR RT from Bombyx mori and a group II intron RT from Eubacterium rectale. Only the non-LTR RT supported robust and serial template jumping, producing one complementary DNA (cDNA) from several templates each copied end to end. We also discovered an unexpected terminal deoxynucleotidyl transferase activity of the RTs that adds nucleotide(s) of choice to 3′ ends of single- and/or double-stranded RNA or DNA. Combining these two types of activity with additional insights about nontemplated nucleotide additions to duplexed cDNA product, we developed a streamlined protocol for fusion of next-generation sequencing adaptors to both cDNA ends in a single RT reaction. When benchmarked using a reference pool of microRNAs (miRNAs), library production by Ordered Two-Template Relay (OTTR) using recombinant non-LTR retroelement RT outperformed all commercially available kits and rivaled the low bias of technically demanding home-brew protocols. We applied OTTR to inventory RNAs purified from extracellular vesicles, identifying miRNAs as well as myriad other noncoding RNAs (ncRNAs) and ncRNA fragments. Our results establish the utility of OTTR for automation-friendly, low-bias, end-to-end RNA sequence inventories of complex ncRNA samples.

The genomic ecosystem of transposable elements in maize

Transposable elements (TEs) constitute the majority of flowering plant DNA, reflecting their tremendous success in subverting, avoiding, and surviving the defenses of their host genomes to ensure their selfish replication. More than 85% of the sequence of the maize genome can be ascribed to past transposition, providing a major contribution to the structure of the genome. Evidence from individual loci has informed our understanding of how transposition has shaped the genome, and a number of individual TE insertions have been causally linked to dramatic phenotypic changes. Genome-wide analyses in maize and other taxa have frequently represented TEs as a relatively homogeneous class of fragmentary relics of past transposition, obscuring their evolutionary history and interaction with their host genome. Using an updated annotation of structurally intact TEs in the maize reference genome, we investigate the family-level dynamics of TEs in maize. Integrating a variety of data, from descriptors of individual TEs like coding capacity, expression, and methylation, as well as similar features of the sequence they inserted into, we model the relationship between attributes of the genomic environment and the survival of TE copies and families. In contrast to the wholesale relegation of all TEs to a single category of junk DNA, these differences reveal a diversity of survival strategies of TE families. Together these generate a rich ecology of the genome, with each TE family representing the evolution of a distinct ecological niche. We conclude that while the impact of transposition is highly family- and context-dependent, a family-level understanding of the ecology of TEs in the genome can refine our ability to predict the role of TEs in generating genetic and phenotypic diversity.

A comprehensive annotation dataset of intact LTR retrotransposons of 300 plant genomes

LTR retrotransposons (LTR-RTs) are ubiquitous and represent the dominant repeat element in plant genomes, playing important roles in functional variation, genome plasticity and evolution. With the advent of new sequencing technologies, a growing number of whole-genome sequences have been made publicly available, making it possible to carry out systematic analyses of LTR-RTs. However, a comprehensive and unified annotation of LTR-RTs in plant groups is still lacking. Here, we constructed a plant intact LTR-RTs dataset, which is designed to classify and annotate intact LTR-RTs with a standardized procedure. The dataset currently comprises a total of 2,593,685 intact LTR-RTs from genomes of 300 plant species representing 93 families of 46 orders. The dataset is accompanied by sequence, diverse structural and functional annotation, age determination and classification information associated with the LTR-RTs. This dataset will contribute valuable resources for investigating the evolutionary dynamics and functional implications of LTR-RTs in plant genomes.

A Sister Lineage of Sampled Retroviruses Corroborates the Complex Evolution of Retroviruses

The origin and deep history of retroviruses remain mysterious and contentious, largely because the diversity of retroviruses is incompletely understood. Here, we report the discovery of lokiretroviruses, a novel major lineage of retroviruses, within the genomes of a wide range of vertebrates (at least 137 species), including lampreys, ray-finned fishes, lobe-finned fishes, amphibians, and reptiles. Lokiretroviruses share a similar genome architecture with known retroviruses, but display some unique features. Interestingly, lokiretrovirus Env proteins share detectable similarity with fusion glycoproteins of viruses within the Mononegavirales order, blurring the boundary between retroviruses and negative sense single-stranded RNA viruses. Phylogenetic analyses based on reverse transcriptase demonstrate that lokiretroviruses are sister to all the retroviruses sampled to date, providing a crucial nexus for studying the deep history of retroviruses. Comparing congruence between host and virus phylogenies suggests lokiretroviruses mainly underwent cross-species transmission. Moreover, we find that retroviruses replaced their ribonuclease H and integrase domains multiple times during their evolutionary course, revealing the importance of domain shuffling in the evolution of retroviruses. Overall, our findings greatly expand our views of the diversity of retroviruses, and provide novel insights into the origin and complex evolutionary history of retroviruses.

Characterizing Repeats in Two Whole-Genome Amplification Methods in the Reniform Nematode Genome

One of the major problems in the U.S. and global cotton production is the damage caused by the reniform nematode, Rotylenchulus reniformis. Amplification of DNA from single nematodes for further molecular analysis can be challenging sometimes. In this research, two whole-genome amplification (WGA) methods were evaluated for their efficiencies in DNA amplification from a single reniform nematode. The WGA was carried out using both REPLI-g Mini and Midi kits, and the GenomePlex single cell whole-genome amplification kit. Sequence analysis produced 4 Mb and 12 Mb of genomic sequences for the reniform nematode using REPLI-g and SIGMA libraries. These sequences were assembled into 28,784 and 24,508 contigs, respectively, for REPLI-g and SIGMA libraries. The highest repeats in both libraries were of low complexity, and the lowest for the REPLI-g library were for satellites and for the SIGMA library, RTE/BOV-B. The same kind of repeats were observed for both libraries; however, the SIGMA library had four other repeat elements (Penelope (long interspersed nucleotide element (LINE)), RTE/BOV-B (LINE), PiggyBac, and Mirage/P-element/Transib), which were not seen in the REPLI-g library. DNA transposons were also found in both libraries. Both reniform nematode 18S rRNA variants (RN_VAR1 and RN_VAR2) could easily be identified in both libraries. This research has therefore demonstrated the ability of using both WGA methods, in amplification of gDNA isolated from single reniform nematodes.

Large Multidomain Protein NMR: HIV-1 Reverse Transcriptase Precursor in Solution

NMR studies of large proteins, over 100 kDa, in solution are technically challenging and, therefore, of considerable interest in the biophysics field. The challenge arises because the molecular tumbling of a protein in solution considerably slows as molecular mass increases, reducing the ability to detect resonances. In fact, the typical ¹H-¹³C or ¹H-¹⁵N correlation spectrum of a large protein, using a ¹³C- or ¹⁵N-uniformly labeled protein, shows severe line-broadening and signal overlap. Selective isotope labeling of methyl groups is a useful strategy to reduce these issues, however, the reduction in the number of signals that goes hand-in-hand with such a strategy is, in turn, disadvantageous for characterizing the overall features of the protein. When domain motion exists in large proteins, the domain motion differently affects backbone amide signals and methyl groups. Thus, the use of multiple NMR probes, such as ¹H, ¹⁹F, ¹³C, and ¹⁵N, is ideal to gain overall structural or dynamical information for large proteins. We discuss the utility of observing different NMR nuclei when characterizing a large protein, namely, the 66 kDa multi-domain HIV-1 reverse transcriptase that forms a homodimer in solution. Importantly, we present a biophysical approach, complemented by biochemical assays, to understand not only the homodimer, p66/p66, but also the conformational changes that contribute to its maturation to a heterodimer, p66/p51, upon HIV-1 protease cleavage.

Stress, Adaptation, and the Deep Genome: Why Transposons Matter

Stress is a common, if often unpredictable life event. It can be defined from an evolutionary perspective as a force an organism perceives it must adapt to. Thus stress is a useful tool to study adaptation and the adaptive capacity of organisms. The deep genome, long neglected as a pile of “junk” has emerged as a source of regulatory DNA and RNA as well as a potential stockpile of adaptive capacity at the organismal and species levels. Recent work on the regulation of transposable elements (TEs), the principle constituents of the deep genome, by stress has shown that these elements are responsive to host stress and other environmental cues. Further, we have shown that some are likely directly regulated by the glucocorticoid receptor (GR), one of the two major vertebrate stress steroid receptors in a fashion that appears adaptive. On the basis of this and other emerging evidence I argue that the deep genome may represent an adaptive toolkit for organisms to respond to their environments at both individual and evolutionary scales. This argues that genomes may be adapted for what Waddington called “trait adaptability” rather than being purely passive objects of natural selection and single nucleotide level mutation.

Measuring Performance Metrics of Machine Learning Algorithms for Detecting and Classifying Transposable Elements

Because of the promising results obtained by machine learning (ML) approaches in several fields, every day is more common, the utilization of ML to solve problems in bioinformatics. In genomics, a current issue is to detect and classify transposable elements (TEs) because of the tedious tasks involved in bioinformatics methods. Thus, ML was recently evaluated for TE datasets, demonstrating better results than bioinformatics applications. A crucial step for ML approaches is the selection of metrics that measure the realistic performance of algorithms. Each metric has specific characteristics and measures properties that may be different from the predicted results. Although the most commonly used way to compare measures is by using empirical analysis, a non-result-based methodology has been proposed, called measure invariance properties. These properties are calculated on the basis of whether a given measure changes its value under certain modifications in the confusion matrix, giving comparative parameters independent of the datasets. Measure invariance properties make metrics more or less informative, particularly on unbalanced, monomodal, or multimodal negative class datasets and for real or simulated datasets. Although several studies applied ML to detect and classify TEs, there are no works evaluating performance metrics in TE tasks. Here, we analyzed 26 different metrics utilized in binary, multiclass, and hierarchical classifications, through bibliographic sources, and their invariance properties. Then, we corroborated our findings utilizing freely available TE datasets and commonly used ML algorithms. Based on our analysis, the most suitable metrics for TE tasks must be stable, even using highly unbalanced datasets, multimodal negative class, and training datasets with errors or outliers. Based on these parameters, we conclude that the F1-score and the area under the precision-recall curve are the most informative metrics since they are calculated based on other metrics, providing insight into the development of an ML application.

DIRS retrotransposons amplify via linear, single-stranded cDNA intermediates

The Dictyostelium Intermediate Repeat Sequence 1 (DIRS-1) is the name-giving member of the DIRS order of tyrosine recombinase retrotransposons. In Dictyostelium discoideum, DIRS-1 is highly amplified and enriched in heterochromatic centromers of the D. discoideum genome. We show here that DIRS-1 it tightly controlled by the D. discoideum RNA interference machinery and is only mobilized in mutants lacking either the RNA dependent RNA polymerase RrpC or the Argonaute protein AgnA. DIRS retrotransposons contain an internal complementary region (ICR) that is thought to be required to reconstitute a full-length element from incomplete RNA transcripts. Using different versions of D. discoideum DIRS-1 equipped with retrotransposition marker genes, we show experimentally that the ICR is in fact essential to complete retrotransposition. We further show that DIRS-1 produces a mixture of single-stranded, mostly linear extrachromosomal cDNA intermediates. If this cDNA is isolated and transformed into D. discoideum cells, it can be used by DIRS-1 proteins to complete productive retrotransposition. This work provides the first experimental evidence to propose a general retrotransposition mechanism of the class of DIRS like tyrosine recombinase retrotransposons.

Becoming a Selfish Clan: Recombination Associated to Reverse-Transcription in LTR Retrotransposons

Transposable elements (TEs) are parasitic DNA bits capable of mobilization and mutagenesis, typically suppressed by host’s epigenetic silencing. Since the selfish DNA concept, it is appreciated that genomes are also molded by arms-races against natural TE inhabitants. However, our understanding of evolutionary processes shaping TEs adaptive populations is scarce. Here, we review the events of recombination associated to reverse-transcription in LTR retrotransposons, a process shuffling their genetic variants during replicative mobilization. Current evidence may suggest that recombinogenic retrotransposons could beneficially exploit host suppression, where clan behavior facilitates their speciation and diversification. Novel refinements to retrotransposons life-cycle and evolution models thus emerge.

Interactions among eukaryotes, retrotransposons and riboviruses: endogenous riboviral elements in eukaryotic genomes

Riboviruses are viruses that have RNA genomes and replicate only via RNA intermediates. Although they do not require a DNA phase for replication and do not encode reverse transcriptase, the presence of DNA forms of riboviral sequences in ribovirus-infected cells has been reported since the 1970s. Additionally, heritable ribovirus-derived sequences, called riboviral endogenous viral elements (EVEs), have been found in the genomes of many eukaryotes. These are now thought to be formed by the reverse transcription machineries of retrotransposons within eukaryotic genomes sometimes referred to as selfish elements. Surprisingly, some reverse-transcribed riboviral DNAs (including EVEs) provide physiological functions for their hosts, suggesting the occurrence of novel interactions among eukaryotic genomes, retrotransposons and riboviruses, and opening the door to new avenues of investigation. Here I review current knowledge on these triangular interactions, and discuss future directions in this field.

Benchmarking transposable element annotation methods for creation of a streamlined, comprehensive pipeline

BACKGROUND

Sequencing technology and assembly algorithms have matured to the point that high-quality de novo assembly is possible for large, repetitive genomes. Current assemblies traverse transposable elements (TEs) and provide an opportunity for comprehensive annotation of TEs. Numerous methods exist for annotation of each class of TEs, but their relative performances have not been systematically compared. Moreover, a comprehensive pipeline is needed to produce a non-redundant library of TEs for species lacking this resource to generate whole-genome TE annotations.

RESULTS

We benchmark existing programs based on a carefully curated library of rice TEs. We evaluate the performance of methods annotating long terminal repeat (LTR) retrotransposons, terminal inverted repeat (TIR) transposons, short TIR transposons known as miniature inverted transposable elements (MITEs), and Helitrons. Performance metrics include sensitivity, specificity, accuracy, precision, FDR, and F₁. Using the most robust programs, we create a comprehensive pipeline called Extensive de-novo TE Annotator (EDTA) that produces a filtered non-redundant TE library for annotation of structurally intact and fragmented elements. EDTA also deconvolutes nested TE insertions frequently found in highly repetitive genomic regions. Using other model species with curated TE libraries (maize and Drosophila), EDTA is shown to be robust across both plant and animal species.

CONCLUSIONS

The benchmarking results and pipeline developed here will greatly facilitate TE annotation in eukaryotic genomes. These annotations will promote a much more in-depth understanding of the diversity and evolution of TEs at both intra- and inter-species levels. EDTA is open-source and freely available: https://github.com/oushujun/EDTA.

Benchmarking transposable element annotation methods for creation of a streamlined, comprehensive pipeline

BACKGROUND

Sequencing technology and assembly algorithms have matured to the point that high-quality de novo assembly is possible for large, repetitive genomes. Current assemblies traverse transposable elements (TEs) and provide an opportunity for comprehensive annotation of TEs. Numerous methods exist for annotation of each class of TEs, but their relative performances have not been systematically compared. Moreover, a comprehensive pipeline is needed to produce a non-redundant library of TEs for species lacking this resource to generate whole-genome TE annotations.

RESULTS

We benchmark existing programs based on a carefully curated library of rice TEs. We evaluate the performance of methods annotating long terminal repeat (LTR) retrotransposons, terminal inverted repeat (TIR) transposons, short TIR transposons known as miniature inverted transposable elements (MITEs), and Helitrons. Performance metrics include sensitivity, specificity, accuracy, precision, FDR, and F1. Using the most robust programs, we create a comprehensive pipeline called Extensive de-novo TE Annotator (EDTA) that produces a filtered non-redundant TE library for annotation of structurally intact and fragmented elements. EDTA also deconvolutes nested TE insertions frequently found in highly repetitive genomic regions. Using other model species with curated TE libraries (maize and Drosophila), EDTA is shown to be robust across both plant and animal species.

CONCLUSIONS

The benchmarking results and pipeline developed here will greatly facilitate TE annotation in eukaryotic genomes. These annotations will promote a much more in-depth understanding of the diversity and evolution of TEs at both intra- and inter-species levels. EDTA is open-source and freely available: https://github.com/oushujun/EDTA.

Evolutionary convergence or homology? Comparative cytogenomics of Caesalpinia group species (Leguminosae) reveals diversification in the pericentromeric heterochromatic composition

We demonstrated by cytogenomic analysis that the proximal heterochromatin of the Northeast Brazilian species of Caesalpinia group is enriched with phylogenetically conserved Ty3/Gypsy-Tekay RT, but diverge in the presence of Ty3/Gypsy-Athila RT and satDNA. The Caesalpinia Group includes 225 species and 27 monophyletic genera of which four occur in Northeastern Brazil: Erythrostemon (1 sp.), Cenostigma (7 spp.), Libidibia (1 sp.), and Paubrasilia (1 sp.). The last three genera are placed in different clades in the Caesalpinia Group phylogeny, and yet they are characterized by having a numerically stable karyotype 2n = 24 (16 M+8A) and GC-rich heterochromatic bands (chromomycin A₃ positive/CMA⁺ bands) in the proximal chromosome regions. To characterize the composition of their heterochromatin and test for the homology of these chromosomal regions, genomic DNA was extracted from Cenostigma microphyllum, Libidibia ferrea, and Paubrasilia echinata, and sequenced at low coverage using the Illumina platform. The genomic repetitive fractions were characterized using a Galaxy/RepeatExplorer-Elixir platform. The most abundant elements of each genome were chromosomally located by fluorescent in situ hybridization (FISH) and compared to the CMA⁺ heterochromatin distribution. The repetitive fraction of the genomes of C. microphyllum, L. ferrea, and P. echinata were estimated to be 41.70%, 38.44%, and 72.51%, respectively. Ty3/Gypsy retrotransposons (RT), specifically the Tekay lineage, were the most abundant repeats in each of the three genomes. FISH mapping revealed species-specific patterns for the Tekay elements in the proximal regions of the chromosomes, co-localized with CMA⁺ bands. Other species-specific patterns were observed, e.g., for the Ty3/Gypsy RT Athila elements which were found in all the proximal heterochromatin of L. ferrea or restricted to the acrocentric chromosomes of C. microphyllum. This Athila labeling co-localized with satellite DNAs (satDNAs). Although the Caesalpinia Group diverged around 55 Mya, our results suggest an ancestral colonization of Tekay RT in the proximal heterochromatin. Thus, the present-day composition of the pericentromeric heterochromatin in these Northeast Brazilian species is a combination of the maintenance of an ancestral Tekay distribution with a species-specific accumulation of other repeats.

Integrative genetic map of repetitive DNA in the sole Solea senegalensis genome shows a Rex transposon located in a proto-sex chromosome

Repetitive sequences play an essential role in the structural and functional evolution of the genome, particularly in the sexual chromosomes. The Senegalese sole (Solea senegalensis) is a valuable flatfish in aquaculture albeit few studies have addressed the mapping and characterization of repetitive DNA families. Here we analyzed the Simple Sequence Repeats (SSRs) and Transposable elements (TEs) content from fifty-seven BAC clones (spanning 7.9 Mb) of this species, located in chromosomes by multiple fluorescence in situ hybridization (m-BAC-FISH) technique. The SSR analysis revealed an average density of 675.1 loci per Mb and a high abundance (59.69%) of dinucleotide coverage was observed, being ‘AC’ the most abundant. An SSR-FISH analysis using eleven probes was also carried out and seven of the 11 probes yielded positive signals. ‘AC’ probes were present as large clusters in almost all chromosomes, supporting the bioinformatic analysis. Regarding TEs, DNA transposons (Class II) were the most abundant. In Class I, LINE elements were the most abundant and the hAT family was the most represented in Class II. Rex/Babar subfamily, observed in two BAC clones mapping to chromosome pair 1, showed the longest match. This chromosome pair has been recently reported as a putative sexual proto-chromosome in this species, highlighting the possible role of the Rex element in the evolution of this chromosome. In the Rex1 phylogenetic tree, the Senegalese sole Rex1 retrotransposon could be associated with one of the four major ancient lineages in fish genomes, in which it is included O. latipes.

Template-switching mechanism of a group II intron-encoded reverse transcriptase and its implications for biological function and RNA-Seq

The reverse transcriptases (RTs) encoded by mobile group II introns and other non-LTR retroelements differ from retroviral RTs in being able to template-switch efficiently from the 5' end of one template to the 3' end of another with little or no complementarity between the donor and acceptor templates. Here, to establish a complete kinetic framework for the reaction and to identify conditions that more efficiently capture acceptor RNAs or DNAs, we used a thermostable group II intron RT (TGIRT; GsI-IIC RT) that can template switch directly from synthetic RNA template/DNA primer duplexes having either a blunt end or a 3'-DNA overhang end. We found that the rate and amplitude of template switching are optimal from starter duplexes with a single nucleotide 3'-DNA overhang complementary to the 3' nucleotide of the acceptor RNA, suggesting a role for nontemplated nucleotide addition of a complementary nucleotide to the 3' end of cDNAs synthesized from natural templates. Longer 3'-DNA overhangs progressively decreased the template-switching rate, even when complementary to the 3' end of the acceptor template. The reliance on only a single bp with the 3' nucleotide of the acceptor together with discrimination against mismatches and the high processivity of group II intron RTs enable synthesis of full-length DNA copies of nucleic acids beginning directly at their 3' end. We discuss the possible biological functions of the template-switching activity of group II intron- and other non-LTR retroelement-encoded RTs, as well as the optimization of this activity for adapter addition in RNA- and DNA-Seq protocols.

Determinants of Active-Site Inhibitor Interaction with HIV-1 RNase H

The ribonuclease H (RNH) activity of HIV-1 reverse transcriptase (RT) is essential for viral replication and can be a target for drug development. Yet, no RNH inhibitor to date has substantial antiviral activity to allow advancement into clinical development. Herein, we describe our characterization of the detailed binding mechanisms of RNH active-site inhibitors, YLC2-155 and ZW566, that bind to the RNH domain through divalent metal ions, using NMR, molecular docking, and quantum mechanical calculations. In the presence of Mg²⁺, NMR spectra of RNH exhibited split (two) resonances for some residues upon inhibitor binding, suggesting two binding modes, an observation consistent with the docking results. The relative populations of the two binding conformers were independent of inhibitor or Mg²⁺ concentration, with one conformation consistently more favored. In our docking study, one distinctive pose of ZW566 showed more interactions with surrounding residues of RNH compared to the analogous binding pose of YLC2-155. Inhibitor titration experiments revealed a lower dissociation constant for ZW566 compared to YLC2-155, in agreement with its higher inhibitory activity. Mg²⁺ titration data also indicated a stronger dependence on Mg²⁺ for the RNH interaction with ZW566 compared to YLC2-155. Combined docking and quantum mechanical calculation results suggest that stronger metal coordination as well as more protein-inhibitor interactions may account for the higher binding affinity of ZW566. These findings support the idea that strategies for the development of potent competitive active site RNH inhibitors should take into account not only metal-inhibitor coordination but also protein-inhibitor interaction and conformational selectivity.

Retrotransposons in Plant Genomes: Structure, Identification, and Classification through Bioinformatics and Machine Learning

Transposable elements (TEs) are genomic units able to move within the genome of virtually all organisms. Due to their natural repetitive numbers and their high structural diversity, the identification and classification of TEs remain a challenge in sequenced genomes. Although TEs were initially regarded as "junk DNA", it has been demonstrated that they play key roles in chromosome structures, gene expression, and regulation, as well as adaptation and evolution. A highly reliable annotation of these elements is, therefore, crucial to better understand genome functions and their evolution. To date, much bioinformatics software has been developed to address TE detection and classification processes, but many problematic aspects remain, such as the reliability, precision, and speed of the analyses. Machine learning and deep learning are algorithms that can make automatic predictions and decisions in a wide variety of scientific applications. They have been tested in bioinformatics and, more specifically for TEs, classification with encouraging results. In this review, we will discuss important aspects of TEs, such as their structure, importance in the evolution and architecture of the host, and their current classifications and nomenclatures. We will also address current methods and their limitations in identifying and classifying TEs.

Study of VIPER and TATE in kinetoplastids and the evolution of tyrosine recombinase retrotransposons

Background

Kinetoplastids are a flagellated group of protists, including some parasites, such as Trypanosoma and Leishmania species, that can cause diseases in humans and other animals. The genomes of these species enclose a fraction of retrotransposons including VIPER and TATE, two poorly studied transposable elements that encode a tyrosine recombinase (YR) and were previously classified as DIRS elements. This study investigated the distribution and evolution of VIPER and TATE in kinetoplastids to understand the relationships of these elements with other retrotransposons.

Results

We observed that VIPER and TATE have a discontinuous distribution among Trypanosomatidae, with several events of loss and degeneration occurring during a vertical transfer evolution. We were able to identify the terminal repeats of these elements for the first time, and we showed that these elements are potentially active in some species, including T. cruzi copies of VIPER. We found that VIPER and TATE are strictly related elements, which were named in this study as VIPER-like. The reverse transcriptase (RT) tree presented a low resolution, and the origin and relationships among YR groups remain uncertain. Conversely, for RH, VIPER-like grouped with Hepadnavirus, whereas for YR, VIPER-like sequences constituted two different clades that are closely allied to Crypton. Distinct topologies among RT, RH and YR trees suggest ancient rearrangements/exchanges in domains and a modular pattern of evolution with putative independent origins for each ORF.

Conclusions

Due to the presence of both elements in Bodo saltans, a nontrypanosomatid species, we suggested that VIPER and TATE have survived and remained active for more than 400 million years or were reactivated during the evolution of the host species. We did not find clear evidence of independent origins of VIPER-like from the other YR retroelements, supporting the maintenance of the DIRS group of retrotransposons. Nevertheless, according to phylogenetic findings and sequence structure obtained by this study and other works, we proposed separating DIRS elements into four subgroups: DIRS-like, PAT-like, Ngaro-like, and VIPER-like.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0175-2) contains supplementary material, which is available to authorized users.

MGERT: a pipeline to retrieve coding sequences of mobile genetic elements from genome assemblies

Background

Genomes of eukaryotes are inhabited by myriads of mobile genetic elements (MGEs) – transposons and retrotransposons - which play a great role in genome plasticity and evolution. A lot of computational tools were developed to annotate them either in genomic assemblies or raw reads using de novo or homology-based approaches. But there has been no pipeline enabling users to get coding and flanking sequences of MGEs suitable for a downstream analysis from genome assemblies.

Results

We developed a new pipeline, MGERT (Mobile Genetic Elements Retrieving Tool), that automates all the steps necessary to obtain protein-coding sequences of mobile genetic elements from genomic assemblies even if no previous knowledge on MGE content of a particular genome is available.

Conclusions

Using MGERT, researchers can easily find MGEs, their coding and flanking sequences in the genome of interest. Thus, this pipeline helps researchers to focus on the biological analysis of MGEs rather than excessive scripting and pipelining.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0163-6) contains supplementary material, which is available to authorized users.

Low coverage sequencing for repetitive DNA analysis in Passiflora edulis Sims: citogenomic characterization of transposable elements and satellite DNA

BACKGROUND

The cytogenomic study of repetitive regions is fundamental for the understanding of morphofunctional mechanisms and genome evolution. Passiflora edulis a species of relevant agronomic value, this work had its genome sequenced by next generation sequencing and bioinformatics analysis performed by RepeatExplorer pipeline. The clusters allowed the identification and characterization of repetitive elements (predominant contributors to most plant genomes). The aim of this study was to identify, characterize and map the repetitive DNA of P. edulis, providing important cytogenomic markers, especially sequences associated with the centromere.

RESULTS

Three clusters of satellite DNAs (69, 118 and 207) and seven clusters of Long Terminal Repeat (LTR) retrotransposons of the superfamilies Ty1/Copy and Ty3/Gypsy and families Angela, Athila, Chromovirus and Maximus-Sire (6, 11, 36, 43, 86, 94 and 135) were characterized and analyzed. The chromosome mapping of satellite DNAs showed two hybridization sites co-located in the 5S rDNA region (PeSat_1), subterminal hybridizations (PeSat_3) and hybridization in four sites, co-located in the 45S rDNA region (PeSat_2). Most of the retroelements hybridizations showed signals scattered in the chromosomes, diverging in abundance, and only the cluster 6 presented pericentromeric regions marking. No satellite DNAs and retroelement associated with centromere was observed.

CONCLUSION

P. edulis has a highly repetitive genome, with the predominance of Ty3/Gypsy LTR retrotransposon. The satellite DNAs and LTR retrotransposon characterized are promising markers for investigation of the evolutionary patterns and genetic distinction of species and hybrids of Passiflora.