Tyrosine Recombinase Retrotransposons and Transposons

Bioinformatic Annotation of Transposon DNA Processing Genes on the Long-Read Genome Assembly of Caenorhabditis elegans

Transposable elements (TEs) or transposons are thought to play roles in animal physiological processes, such as germline, early embryonic, and brain development, as well as aging. However, their roles have not been systematically investigated through experimental studies. In this study, we created a catalog of genes directly involved in replication, excision, or integration of transposon-coding DNA, which we refer to as transposon DNA processing genes (TDPGs). Specifically, to bridge the gap to experimental studies, we sought potentially functional TDPGs which maintain intact open reading frames and the amino acids at their catalytic cores on the latest long-read genome assembly of Caenorhabditis elegans, VC2010. Among 52 519 TE loci, we identified 145 potentially functional TDPGs encoded in long terminal repeat elements, long interspersed nuclear elements, terminal inverted repeat elements, Helitrons, and Mavericks/Polintons. Our TDPG catalog, which contains a feasible number of genes, allows for the experimental manipulation of TE mobility in vivo, regardless of whether the TEs are autonomous or non-autonomous, thereby potentially promoting the study of the physiological functions of TE mobility.

Starships: a new frontier for fungal biology

Transposable elements (TEs) are semiautonomous genetic entities that proliferate in genomes. We recently discovered the Starships, a previously hidden superfamily of giant TEs found in a diverse subphylum of filamentous fungi, the Pezizomycotina. Starships are unlike other eukaryotic TEs because they have evolved mechanisms for both mobilizing entire genes, including those encoding conditionally beneficial phenotypes, and for horizontally transferring between individuals. We argue that Starships have unrivaled capacity to engage their fungal hosts as genetic parasites and mutualists, revealing unexplored terrain for investigating the ecoevolutionary dynamics of TE-eukaryote interactions. We build on existing models of fungal genome evolution by conceptualizing Starships as a distinct genomic compartment whose dynamics profoundly shape fungal biology.

DNA methylation enables recurrent endogenization of giant viruses in an animal relative

5-Methylcytosine (5mC) is a widespread silencing mechanism that controls genomic parasites. In eukaryotes, 5mC has gained complex roles in gene regulation beyond parasite control, yet 5mC has also been lost in many lineages. The causes for 5mC retention and its genomic consequences are still poorly understood. Here, we show that the protist closely related to animals Amoebidium appalachense features both transposon and gene body methylation, a pattern reminiscent of invertebrates and plants. Unexpectedly, hypermethylated genomic regions in Amoebidium derive from viral insertions, including hundreds of endogenized giant viruses, contributing 14% of the proteome. Using a combination of inhibitors and genomic assays, we demonstrate that 5mC silences these giant virus insertions. Moreover, alternative Amoebidium isolates show polymorphic giant virus insertions, highlighting a dynamic process of infection, endogenization, and purging. Our results indicate that 5mC is critical for the controlled coexistence of newly acquired viral DNA into eukaryotic genomes, making Amoebidium a unique model to understand the hybrid origins of eukaryotic DNA.

A beginner's guide to manual curation of transposable elements

BACKGROUND

In the study of transposable elements (TEs), the generation of a high confidence set of consensus sequences that represent the diversity of TEs found in a given genome is a key step in the path to investigate these fascinating genomic elements. Many algorithms and pipelines are available to automatically identify putative TE families present in a genome. Despite the availability of these valuable resources, producing a library of high-quality full-length TE consensus sequences largely remains a process of manual curation. This know-how is often passed on from mentor-to-mentee within research groups, making it difficult for those outside the field to access this highly specialised skill.

RESULTS

Our manuscript attempts to fill this gap by providing a set of detailed computer protocols, software recommendations and video tutorials for those aiming to manually curate TEs. Detailed step-by-step protocols, aimed at the complete beginner, are presented in the Supplementary Methods.

CONCLUSIONS

The proposed set of programs and tools presented here will make the process of manual curation achievable and amenable to all researchers and in special to those new to the field of TEs.

Evolutionary dynamics of DIRS-like and Ngaro-like retrotransposons in Xenopus laevis and Xenopus tropicalis genomes

Anuran genomes have a large number and diversity of transposable elements, but are little explored, mainly in relation to their molecular structure and evolutionary dynamics. Here, we investigated the retrotransposons containing tyrosine recombinase (YR) (order DIRS) in the genome of Xenopus tropicalis and Xenopus laevis. These anurans show 2n = 20 and the 2n = 36 karyotypes, respectively. They diverged about 48 million years ago (mya) and X. laevis had an allotetraploid origin (around 17–18 mya). Our investigation is based on the analysis of the molecular structure and the phylogenetic relationships of 95 DIRS families of Xenopus belonging to DIRS-like and Ngaro-like superfamilies. We were able to identify molecular signatures in the 5' and 3' noncoding terminal regions, preserved open reading frames, and conserved domains that are specific to distinguish each superfamily. We recognize two ancient amplification waves of DIRS-like elements that occurred in the ancestor of both species and a higher density of the old/degenerate copies detected in both subgenomes of X. laevis. More recent amplification waves are seen in X. tropicalis (less than 3.2 mya) and X. laevis (around 10 mya) corroborating with transcriptional activity evidence. All DIRS-like families were found in both X. laevis subgenomes, while a few were most represented in the L subgenome. Ngaro-like elements presented less diversity and quantity in X. tropicalis and X. laevis genomes, although potentially active copies were found in both species and this is consistent with a recent amplification wave seen in the evolutionary landscape. Our findings highlight a differential diversity-level and evolutionary dynamics of the YR retrotransposons in X. tropicalis and X. laevis species expanding our comprehension of the behavior of these elements in both genomes during the diversification process.

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.

Differential expression of transposable elements in the medaka melanoma model

Malignant melanoma incidence is rising worldwide. Its treatment in an advanced state is difficult, and the prognosis of this severe disease is still very poor. One major source of these difficulties is the high rate of metastasis and increased genomic instability leading to a high mutation rate and the development of resistance against therapeutic approaches. Here we investigate as one source of genomic instability the contribution of activation of transposable elements (TEs) within the tumor. We used the well-established medaka melanoma model and RNA-sequencing to investigate the differential expression of TEs in wildtype and transgenic fish carrying melanoma. We constructed a medaka-specific TE sequence library and identified TE sequences that were specifically upregulated in tumors. Validation by qRT- PCR confirmed a specific upregulation of a LINE and an LTR element in malignant melanomas of transgenic fish.

Virophages and retrotransposons colonize the genomes of a heterotrophic flagellate

Virophages can parasitize giant DNA viruses and may provide adaptive anti-giant virus defense in unicellular eukaryotes. Under laboratory conditions, the virophage mavirus integrates into the nuclear genome of the marine flagellate Cafeteria burkhardae and reactivates upon superinfection with the giant virus CroV. In natural systems, however, the prevalence and diversity of host-virophage associations has not been systematically explored. Here, we report dozens of integrated virophages in four globally sampled C. burkhardae strains that constitute up to 2% of their host genomes. These endogenous mavirus-like elements (EMALEs) separated into eight types based on GC-content, nucleotide similarity, and coding potential and carried diverse promoter motifs implicating interactions with different giant viruses. Between host strains, some EMALE insertion loci were conserved indicating ancient integration events, whereas the majority of insertion sites were unique to a given host strain suggesting that EMALEs are active and mobile. Furthermore, we uncovered a unique association between EMALEs and a group of tyrosine recombinase retrotransposons, revealing yet another layer of parasitism in this nested microbial system. Our findings show that virophages are widespread and dynamic in wild Cafeteria populations, supporting their potential role in antiviral defense in protists.

Viruses exist in all ecosystems in vast numbers and infect many organisms. Some of them are harmful but others can protect the organisms they infect. For example, a group of viruses called virophages protect microscopic sea creatures called plankton from deadly infections by so-called giant viruses. In fact, virophages need plankton infected with giant viruses to survive because they use enzymes from the giant viruses to turn on their own genes.

A virophage called mavirus integrates its genes into the DNA of a type of plankton called Cafeteria. It lays dormant in the DNA until a giant virus called CroV infects the plankton. This suggests that the mavirus may be a built-in defense against CroV infections and laboratory studies seem to confirm this. But whether wild Cafeteria also use these defenses is unknown.

Hackl et al. show that virophages are common in the DNA of wild Cafeteria and that the two appear to have a mutually beneficial relationship. In the experiments, the researchers sequenced the genomes of four Cafeteria populations from the Atlantic and Pacific Oceans and looked for virophages in their DNA. Each of the four Cafeteria genomes contained dozens of virophages, which suggests that virophages are important to these plankton. This included several relatives of the mavirus and seven new virophages. Virophage genes were often interrupted by so called jumping genes, which may take advantage of the virophages the way the virophages use giant viruses to meet their own needs.

The experiments show that virophages often co-exist with marine plankton from around the world and these relationships are likely beneficial. In fact, the experiments suggest that the virophages may have played an important role in the evolution of these plankton. Further studies may help scientists learn more about virus ecology and how viruses have shaped the evolution of other creatures.

Transposable elements expression in Rhinella marina (cane toad) specimens submitted to immune and stress challenge

Transposable elements (TEs) are important components of eukaryotic genomes and compose around 30% of the genome of Rhinella marina, an invasive toad species. Considering the possible role of TEs in the adaptation of populations, we have analyzed the expression of TEs in publicly available spleen tissue transcriptomic data generated for this species after immune and stress challenge. By analyzing the transcriptome assembly, we detected a high number of TE segments. Moreover, some distinct TE families were differentially expressed in some conditions. Our result shows that several TEs are capable of being transcribed in R. marina and they could help to generate a rapid response of specimens to the environment. Also, we can suggest that these TEs could be activated in the germinative cells as well producing variability to be selected and shaped by the evolutionary processes behind the success of this invasive species. Thus, the TEs are important targets for investigation in the context of R. marina adaptation.

Comparative genomics of Chlamydomonas

Despite its role as a reference organism in the plant sciences, the green alga Chlamydomonas reinhardtii entirely lacks genomic resources from closely related species. We present highly contiguous and well-annotated genome assemblies for three unicellular C. reinhardtii relatives: Chlamydomonas incerta, Chlamydomonas schloesseri, and the more distantly related Edaphochlamys debaryana. The three Chlamydomonas genomes are highly syntenous with similar gene contents, although the 129.2 Mb C. incerta and 130.2 Mb C. schloesseri assemblies are more repeat-rich than the 111.1 Mb C. reinhardtii genome. We identify the major centromeric repeat in C. reinhardtii as a LINE transposable element homologous to Zepp (the centromeric repeat in Coccomyxa subellipsoidea) and infer that centromere locations and structure are likely conserved in C. incerta and C. schloesseri. We report extensive rearrangements, but limited gene turnover, between the minus mating type loci of these Chlamydomonas species. We produce an eight-species core-Reinhardtinia whole-genome alignment, which we use to identify several hundred false positive and missing genes in the C. reinhardtii annotation and >260,000 evolutionarily conserved elements in the C. reinhardtii genome. In summary, these resources will enable comparative genomics analyses for C. reinhardtii, significantly extending the analytical toolkit for this emerging model system.

Cell-cycle-dependent EBNA1-DNA crosslinking promotes replication termination at oriP and viral episome maintenance

Epstein-Barr virus (EBV) is an oncogenic human herpesvirus that persists as a multicopy episome in proliferating host cells. Episome maintenance is strictly dependent on EBNA1, a sequence-specific DNA-binding protein with no known enzymatic activities. Here, we show that EBNA1 forms a cell cycle-dependent DNA crosslink with the EBV origin of plasmid replication oriP. EBNA1 tyrosine 518 (Y518) is essential for crosslinking to oriP and functionally required for episome maintenance and generation of EBV-transformed lymphoblastoid cell lines (LCLs). Mechanistically, Y518 is required for replication fork termination at oriP in vivo and for formation of SDS-resistant complexes in vitro. EBNA1-DNA crosslinking corresponds to single-strand endonuclease activity specific to DNA structures enriched at replication-termination sites, such as 4-way junctions. These findings reveal that EBNA1 forms tyrosine-dependent DNA-protein crosslinks and single-strand cleavage at oriP required for replication termination and viral episome maintenance.

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.

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.

Mobilization of retrotransposons as a cause of chromosomal diversification and rapid speciation: the case for the Antarctic teleost genus Trematomus

Background

The importance of transposable elements (TEs) in the genomic remodeling and chromosomal rearrangements that accompany lineage diversification in vertebrates remains the subject of debate. The major impediment to understanding the roles of TEs in genome evolution is the lack of comparative and integrative analyses on complete taxonomic groups. To help overcome this problem, we have focused on the Antarctic teleost genus Trematomus (Notothenioidei: Nototheniidae), as they experienced rapid speciation accompanied by dramatic chromosomal diversity. Here we apply a multi-strategy approach to determine the role of large-scale TE mobilization in chromosomal diversification within Trematomus species.

Results

Despite the extensive chromosomal rearrangements observed in Trematomus species, our measurements revealed strong interspecific genome size conservation. After identifying the DIRS1, Gypsy and Copia retrotransposon superfamilies in genomes of 13 nototheniid species, we evaluated their diversity, abundance (copy numbers) and chromosomal distribution. Four families of DIRS1, nine of Gypsy, and two of Copia were highly conserved in these genomes; DIRS1 being the most represented within Trematomus genomes. Fluorescence in situ hybridization mapping showed preferential accumulation of DIRS1 in centromeric and pericentromeric regions, both in Trematomus and other nototheniid species, but not in outgroups: species of the Sub-Antarctic notothenioid families Bovichtidae and Eleginopsidae, and the non-notothenioid family Percidae.

Conclusions

In contrast to the outgroups, High-Antarctic notothenioid species, including the genus Trematomus, were subjected to strong environmental stresses involving repeated bouts of warming above the freezing point of seawater and cooling to sub-zero temperatures on the Antarctic continental shelf during the past 40 millions of years (My). As a consequence of these repetitive environmental changes, including thermal shocks; a breakdown of epigenetic regulation that normally represses TE activity may have led to sequential waves of TE activation within their genomes. The predominance of DIRS1 in Trematomus species, their transposition mechanism, and their strategic location in “hot spots” of insertion on chromosomes are likely to have facilitated nonhomologous recombination, thereby increasing genomic rearrangements. The resulting centric and tandem fusions and fissions would favor the rapid lineage diversification, characteristic of the nototheniid adaptive radiation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4714-x) contains supplementary material, which is available to authorized users.

Using bioinformatic and phylogenetic approaches to classify transposable elements and understand their complex evolutionary histories

In recent years, much attention has been paid to comparative genomic studies of transposable elements (TEs) and the ensuing problems of their identification, classification, and annotation. Different approaches and diverse automated pipelines are being used to catalogue and categorize mobile genetic elements in the ever-increasing number of prokaryotic and eukaryotic genomes, with little or no connectivity between different domains of life. Here, an overview of the current picture of TE classification and evolutionary relationships is presented, updating the diversity of TE types uncovered in sequenced genomes. A tripartite TE classification scheme is proposed to account for their replicative, integrative, and structural components, and the need to expand in vitro and in vivo studies of their structural and biological properties is emphasized. Bioinformatic studies have now become front and center of novel TE discovery, and experimental pursuits of these discoveries hold great promise for both basic and applied science.

Sub-Minimum Inhibitory Concentrations of Rhubarb Water Extracts Inhibit Streptococcus suis Biofilm Formation

Streptococcus suis is one of the most important swine pathogens, which can cause persistent infection by forming biofilms. In this study, sub-minimum inhibitory concentration (sub-MIC) of rhubarb water extracts were found to inhibit biofilm formation. Two-component signal transduction systems (TCSs), transcriptional regulators, and DNA binding proteins were compared under two conditions: (1) cells treated with sub-MIC rhubarb water extracts and (2) untreated cells. Using an isobaric tags for relative and absolute quantitation (iTRAQ) strategy, we found that TCSs constituent proteins of histidine kinase and response regulator were significantly down-regulated. This down-regulation can affect the transfer of information during biofilm formation. The transcriptional regulators and DNA binding proteins that can interact with TCSs and interrupt gene transcription were also significantly altered. For these reasons, the levels of protein expressions varied in different parts of the treated vs. untreated cells. In summary, rhubarb water extracts might serve as potential inhibitor for the control of S. suis biofilm formation. The change in TCSs, transcriptional regulators, and DNA binding proteins may be important factors in S. suis biofilm inhibition.

Identification of the agg1 mutation responsible for negative phototaxis in a "wild-type" strain of Chlamydomonas reinhardtii

The unicellular green alga Chlamydomonas reinhardtii is a model organism for various studies in biology. CC-124 is a laboratory strain widely used as a wild type. However, this strain is known to carry agg1 mutation, which causes cells to swim away from the light source (negative phototaxis), in contrast to the cells of other wild-type strains, which swim toward the light source (positive phototaxis). Here we identified the causative gene of agg1 (AGG1) using AFLP-based gene mapping and whole genome next-generation sequencing. This gene encodes a 36-kDa protein containing a Fibronectin type III domain and a CHORD-Sgt1 (CS) domain. The gene product is localized to the cell body and not to flagella or basal body.

•

agg1, a mutation harbored in CC-124 (a widely used wild-type strain), was identified.

•

A retrotransposon, TOC1, is inserted in the 5′ UTR of AGG1 gene (Cre13.g590400).

•

AGG1 encodes a 36-kDa protein of unknown function.

•

This protein localizes to the cell body but not to the eyespot, flagella or basal body.

Not so bad after all: retroviruses and long terminal repeat retrotransposons as a source of new genes in vertebrates

Viruses and transposable elements, once considered as purely junk and selfish sequences, have repeatedly been used as a source of novel protein-coding genes during the evolution of most eukaryotic lineages, a phenomenon called 'molecular domestication'. This is exemplified perfectly in mammals and other vertebrates, where many genes derived from long terminal repeat (LTR) retroelements (retroviruses and LTR retrotransposons) have been identified through comparative genomics and functional analyses. In particular, genes derived from gag structural protein and envelope (env) genes, as well as from the integrase-coding and protease-coding sequences, have been identified in humans and other vertebrates. Retroelement-derived genes are involved in many important biological processes including placenta formation, cognitive functions in the brain and immunity against retroelements, as well as in cell proliferation, apoptosis and cancer. These observations support an important role of retroelement-derived genes in the evolution and diversification of the vertebrate lineage.