Transposable Element Domestication As an Adaptation to Evolutionary Conflicts

Integrated Analysis of Large-Scale Omics Data Revealed Relationship Between Tissue Specificity and Evolutionary Dynamics of Small RNAs in Maize (Zea mays)

The evolutionary dynamics and tissue specificity of protein-coding genes are well documented in plants. However, the evolutionary consequences of small RNAs (sRNAs) on tissue-specific functions remain poorly understood. Here, we performed integrated analysis of 195 deeply sequenced sRNA libraries of maize B73, representing more than 10 tissues, and identified a comprehensive list of 419 maize microRNA (miRNA) genes, 271 of which were newly discovered in this study. We further characterized the evolutionary dynamics and tissue specificity of miRNA genes and corresponding miRNA isoforms (isomiRs). Our analysis revealed that tissue specificity of isomiR events tends to be associated with miRNA gene abundance and suggested that the frequencies of isomiR types are affected by the local genomic regions. Moreover, genome duplication (GD) events have dramatic effect on evolutionary dynamics of maize miRNA genes, and the abundance divergence for tissue-specific miRNA genes is associated with GD events. Further study indicated that duplicate miRNA genes with tissue-specific expression patterns, such as miR2275a, a phased siRNA (phasiRNA) trigger, contribute to phenotypic traits in maize. Additionally, our study revealed the expression preference of 21- and 24-nt phasiRNAs in relation to tissue specificity. This large-scale sRNAomic study depicted evolutionary implications of tissue-specific maize sRNAs, which coordinate genome duplication, isomiR modification, phenotypic traits and phasiRNAs differentiation.

Complete Genome Sequence Reveals Evolutionary and Comparative Genomic Features of Xanthomonas albilineans Causing Sugarcane Leaf Scald

Leaf scald (caused by Xanthomonas albilineans) is an important bacterial disease affecting sugarcane in most sugarcane growing countries, including China. High genetic diversity exists among strains of X. albilineans from diverse geographic regions. To highlight the genomic features associated with X. albilineans from China, we sequenced the complete genome of a representative strain (Xa-FJ1) of this pathogen using the PacBio and Illumina platforms. The complete genome of strain Xa-FJ1 consists of a circular chromosome of 3,724,581 bp and a plasmid of 31,536 bp. Average nucleotide identity analysis revealed that Xa-FJ1 was closest to five strains from the French West Indies and the USA, particularly to the strain GPE PC73 from Guadeloupe. Comparative genomic analysis between Xa-FJ1 and GPE PC73 revealed prophage integration, homologous recombination, transposable elements, and a clustered regulatory interspaced short palindromic repeats (CRISPR) system that were linked with 16 insertions/deletions (InDels). Ten and 82 specific genes were found in Xa-FJ1 and GPE PC73, respectively, and some of these genes were subjected to phage-related proteins, zona occludens toxin, and DNA methyltransferases. Our findings highlight intra-species genetic variability of the leaf scald pathogen and provide additional genomic resources to investigate its fitness and virulence.

Compensating for over-production inhibition of the Hsmar1 transposon in Escherichia coli using a series of constitutive promoters

Background

Transposable elements (TEs) are a diverse group of self-mobilizing DNA elements. Transposition has been exploited as a powerful tool for molecular biology and genomics. However, transposition is sometimes limited because of auto-regulatory mechanisms that presumably allow them to cohabit within their hosts without causing excessive genomic damage. The papillation assay provides a powerful visual screen for hyperactive transposases. Transposition is revealed by the activation of a promoter-less lacZ gene when the transposon integrates into a non-essential gene on the host chromosome. Transposition events are detected as small blue speckles, or papillae, on the white background of the main Escherichia coli colony.

Results

We analysed the parameters of the papillation assay including the strength of the transposase transcriptional and translational signals. To overcome certain limitations of inducible promoters, we constructed a set of vectors based on constitutive promoters of different strengths to widen the range of transposase expression. We characterized and validated our expression vectors with Hsmar1, a member of the mariner transposon family. The highest rate of transposition was observed with the weakest promoters. We then took advantage of our approach to investigate how the level of transposition responds to selected point mutations and the effect of joining the transposase monomers into a single-chain dimer.

Conclusions

We generated a set of vectors to provide a wide range of transposase expression which will be useful for screening libraries of transposase mutants. The use of weak promoters should allow screening for truly hyperactive transposases rather than those that are simply resistant to auto-regulatory mechanisms, such as overproduction inhibition (OPI). We also found that mutations in the Hsmar1 dimer interface provide resistance to OPI in bacteria, which could be valuable for improving bacterial transposon mutagenesis techniques.

Evolution of the RNA N 6-Methyladenosine Methylome Mediated by Genomic Duplication

RNA N ⁶-methyladenosine (m⁶A) modification is the most abundant form of RNA epigenetic modification in eukaryotes. Given that m⁶A evolution is associated with the selective constraints of nucleotide sequences in mammalian genomes, we hypothesize that m⁶A evolution can be linked, at least in part, to genomic duplication events in complex polyploid plant genomes. To test this hypothesis, we presented the maize (Zea mays) m⁶A modification landscape in a transcriptome-wide manner and identified 11,968 m⁶A peaks carried by 5,893 and 3,811 genes from two subgenomes (maize1 and maize2, respectively). Each of these subgenomes covered over 2,200 duplicate genes. Within these duplicate genes, those carrying m⁶A peaks exhibited significant differences in retention rate. This biased subgenome fractionation of m⁶A-methylated genes is associated with multiple sequence features and is influenced by asymmetric evolutionary rates. We also characterized the coevolutionary patterns of m⁶A-methylated genes and transposable elements, which can be mediated by whole genome duplication and tandem duplication. We revealed the evolutionary conservation and divergence of duplicated m⁶A functional factors and the potential role of m⁶A modification in maize responses to drought stress. This study highlights complex interplays between m⁶A modification and gene duplication, providing a reference for understanding the mechanisms underlying m⁶A evolution mediated by genome duplication events.

Telomere-Specialized Retroelements in Drosophila: Adaptive Symbionts of the Genome, Neutral, or in Conflict?

Linear chromosomes shorten in every round of replication. In Drosophila, telomere-specialized long interspersed retrotransposable elements (LINEs) belonging to the jockey clade offset this shortening by forming head-to-tail arrays at Drosophila telomere ends. As such, these telomeric LINEs have been considered adaptive symbionts of the genome, protecting it from premature decay, particularly as Drosophila lacks a conventional telomerase holoenzyme. However, as reviewed here, recent work reveals a high degree of variation and turnover in the telomere-specialized LINE lineages across Drosophila. There appears to be no absolute requirement for LINE activity to maintain telomeres in flies, hence the suggestion that the telomere-specialized LINEs may instead be neutral or in conflict with the host, rather than adaptive.

HIV‑1 integrase inhibitors targeting various DDE transposases: Retroviral integration versus RAG‑mediated recombination (Review)

Transposases are ubiquitous mobile genetic elements responsible for genome development, driving rearrangements, such as insertions, deletions and translocations. Across species evolution, some transposases are tamed by their host and are made part of complex cellular systems. The proliferation of retroviruses is also dependent on transposase related enzymes termed integrases. Recombination‑activating gene protein (RAG)1 and metnase are just two examples of transposase domestication and together with retroviral integrases (INs), they belong to the DDE polynucleotidyl transferases superfamily. They share mechanistic and structural features linked to the RNase H‑like fold, harboring a DDE(D) metal dependent catalytic motif. Recent antiretroviral compounds target the catalytic domain of integrase, but they also have the potential of inhibiting other related enzymes. In this review, we report the activity of different classes of integrase inhibitors on various DDE transposases. Computational simulations are useful to predict the extent of off‑target activity and have been employed to study the interactions between RAG1 recombinase and compounds from three different pharmacologic classes. We demonstrate that strand‑transfer inhibitors display a higher affinity towards the RAG1 RNase H domain, as suggested by experimental data compared to allosteric inhibitors. While interference with RAG1 and 2 recombination is associated with a negative impact on immune function, the inhibition of metnase or HTLV‑1 integrase opens the way for the development of novel therapies for refractory cancers.

Has gene expression neofunctionalization in the fire ant antennae contributed to queen discrimination behavior?

Queen discrimination behavior in the fire ant Solenopsis invicta maintains its two types of societies: colonies with one (monogyne) or many (polygyne) queens, yet the underlying genetic mechanism is poorly understood. This behavior is controlled by two supergene alleles, SB and Sb, with ~600 genes. Polygyne workers, having either the SB/SB or SB/Sb genotype, accept additional SB/Sb queens into their colonies but kill SB/SB queens. In contrast, monogyne workers, all SB/SB, reject all additional queens regardless of genotype. Because the SB and Sb alleles have suppressed recombination, determining which genes within the supergene mediate this differential worker behavior is difficult. We hypothesized that the alternate worker genotypes sense queens differently because of the evolution of differential expression of key genes in their main sensory organ, the antennae. To identify such genes, we sequenced RNA from four replicates of pooled antennae from three classes of workers: monogyne SB/SB, polygyne SB/SB, and polygyne SB/Sb. We identified 81 differentially expressed protein-coding genes with 13 encoding potential chemical metabolism or perception proteins. We focused on the two odorant perception genes: an odorant receptor SiOR463 and an odorant-binding protein SiOBP12. We found that SiOR463 has been lost in the Sb genome. In contrast, SiOBP12 has an Sb-specific duplication, SiOBP12b', which is expressed in the SB/Sb worker antennae, while both paralogs are expressed in the body. Comparisons with another fire ant species revealed that SiOBP12b' antennal expression is specific to S. invicta and suggests that queen discrimination may have evolved, in part, through expression neofunctionalization.

Evolutionary history and classification of Micropia retroelements in Drosophilidae species

Transposable elements (TEs) have the main role in shaping the evolution of genomes and host species, contributing to the creation of new genes and promoting rearrangements frequently associated with new regulatory networks. Support for these hypotheses frequently results from studies with model species, and Drosophila provides a great model organism to the study of TEs. Micropia belongs to the Ty3/Gypsy group of long terminal repeats (LTR) retroelements and comprises one of the least studied Drosophila transposable elements. In this study, we assessed the evolutionary history of Micropia within Drosophilidae, while trying to assist in the classification of this TE. At first, we performed searches of Micropia presence in the genome of natural populations from several species. Then, based on searches within online genomic databases, we retrieved Micropia-like sequences from the genomes of distinct Drosophilidae species. We expanded the knowledge of Micropia distribution within Drosophila species. The Micropia retroelements we detected consist of an array of divergent sequences, which we subdivided into 20 subfamilies. Even so, a patchy distribution of Micropia sequences within the Drosophilidae phylogeny could be identified, with incongruences between the species phylogeny and the Micropia phylogeny. Comparing the pairwise synonymous distance (dS) values between Micropia and three host nuclear sequences, we found several cases of unexpectedly high levels of similarity between Micropia sequences in divergent species. All these findings provide a hypothesis to the evolution of Micropia within Drosophilidae, which include several events of vertical and horizontal transposon transmission, associated with ancestral polymorphisms and recurrent Micropia sequences diversification.

Defensive Function of Transposable Elements in Bacteria

It has been widely debated whether transposable elements have a positive or a negative effect on their host cells. This study demonstrated that transposable elements, specifically insertion sequences (ISs), can adopt a defensive role in Escherichia coli. In three different E. coli strains (S17, DH5α, and Nissle 1917), IS1 and IS10 rapidly disrupted the I-CeuI gene (encoding I-CeuI endonuclease) on the plasmid pLO11-ICeuI as early as the first generation, despite the gene-circuit being under control of an arabinose promoter. Proteomics analysis showed that the protein abundance profile of E. coli DH5α with pLO11-ICeuI in the fifth generation was nearly opposite to that of control strain (E. coli with pLO11, no I-CeuI). The DNA damage caused by the leaky expression of I-CeuI was enough to trigger a SOS response and alter lipid synthesis, ribosomal activity, RNA/DNA metabolism, central dogma and cell cycle processes in E. coli DH5α. After the ISs disrupted the expression of I-CeuI, cells fully recovered by the 31st generation had a protein abundance profile similar to that of the control strain. This study showed that ISs readily mutated a harmful gene which subsequently restored host fitness. These observations have implications for the stability of designed gene circuits in synthetic biology.

Environmentally-Induced Transgenerational Epigenetic Inheritance: Implication of PIWI Interacting RNAs

Environmentally-induced transgenerational epigenetic inheritance is an emerging field. The understanding of associated epigenetic mechanisms is currently in progress with open questions still remaining. In this review, we present an overview of the knowledge of environmentally-induced transgenerational inheritance and associated epigenetic mechanisms, mainly in animals. The second part focuses on the role of PIWI-interacting RNAs (piRNAs), a class of small RNAs involved in the maintenance of the germline genome, in epigenetic memory to put into perspective cases of environmentally-induced transgenerational inheritance involving piRNA production. Finally, the last part addresses how genomes are facing production of new piRNAs, and from a broader perspective, how this process might have consequences on evolution and on sporadic disease development.

Host-transposon interactions: conflict, cooperation, and cooption

Transposable elements (TEs) are mobile DNA sequences that colonize genomes and threaten genome integrity. As a result, several mechanisms appear to have emerged during eukaryotic evolution to suppress TE activity. However, TEs are ubiquitous and account for a prominent fraction of most eukaryotic genomes. We argue that the evolutionary success of TEs cannot be explained solely by evasion from host control mechanisms. Rather, some TEs have evolved commensal and even mutualistic strategies that mitigate the cost of their propagation. These coevolutionary processes promote the emergence of complex cellular activities, which in turn pave the way for cooption of TE sequences for organismal function.

Human SETMAR is a DNA sequence-specific histone-methylase with a broad effect on the transcriptome

Transposons impart dynamism to the genomes they inhabit and their movements frequently rewire the control of nearby genes. Occasionally, their proteins are domesticated when they evolve a new function. SETMAR is a protein methylase with a sequence-specific DNA binding domain. It began to evolve about 50 million years ago when an Hsmar1 transposon integrated downstream of a SET-domain methylase gene. Here we show that the DNA-binding domain of the transposase targets the enzyme to transposon-end remnants and that this is capable of regulating gene expression, dependent on the methylase activity. When SETMAR was modestly overexpressed in human cells, almost 1500 genes changed expression by more than 2-fold (65% up- and 35% down-regulated). These genes were enriched for the KEGG Pathways in Cancer and include several transcription factors important for development and differentiation. Expression of a similar level of a methylase-deficient SETMAR changed the expression of many fewer genes, 77% of which were down-regulated with no significant enrichment of KEGG Pathways. Our data is consistent with a model in which SETMAR is part of an anthropoid primate-specific regulatory network centered on the subset of genes containing a transposon end.

Transposition favors the generation of large effect mutations that may facilitate rapid adaption

Transposable elements (TEs) are mobile parasitic sequences that have been repeatedly coopted during evolution to generate new functions and rewire gene regulatory networks. Yet, the contribution of active TEs to the creation of heritable mutations remains unknown. Using TE accumulation lines in Arabidopsis thaliana we show that once initiated, transposition produces an exponential spread of TE copies, which rapidly leads to high mutation rates. Most insertions occur near or within genes and targets differ between TE families. Furthermore, we uncover an essential role of the histone variant H2A.Z in the preferential integration of Ty1/copia retrotransposons within environmentally responsive genes and away from essential genes. We also show that epigenetic silencing of new Ty1/copia copies can affect their impact on major fitness-related traits, including flowering time. Our findings demonstrate that TEs are potent episodic (epi)mutagens that, thanks to marked chromatin tropisms, limit the mutation load and increase the potential for rapid adaptation.

Impaired plasma membrane localization of ubiquitin ligase complex underlies 3-M syndrome development

3-M primordial dwarfism is an inherited disease characterized by severe pre- and postnatal growth retardation and by mutually exclusive mutations in three genes, CUL7, OBSL1, and CCDC8. The mechanism underlying 3-M dwarfism is not clear. We showed here that CCDC8, derived from a retrotransposon Gag protein in placental mammals, exclusively localized on the plasma membrane and was phosphorylated by CK2 and GSK3. Phosphorylation of CCDC8 resulted in its binding first with OBSL1, and then CUL7, leading to the membrane assembly of the 3-M E3 ubiquitin ligase complex. We identified LL5β, a plasma membrane protein that regulates cell migration, as a substrate of 3-M ligase. Wnt inhibition of CCDC8 phosphorylation or patient-derived mutations in 3-M genes disrupted membrane localization of the 3-M complex and accumulated LL5β. Deletion of Ccdc8 in mice impaired trophoblast migration and placental development, resulting in intrauterine growth restriction and perinatal lethality. These results identified a mechanism regulating cell migration and placental development that underlies the development of 3-M dwarfism.

Heterogeneous transposable elements as silencers, enhancers and targets of meiotic recombination

During meiosis, DNA double-strand breaks are initiated by the topoisomerase-like enzyme SPO11 and are repaired by inter-sister chromatid and inter-homologue DNA repair pathways. Genome-wide maps of initiating DNA double-strand breaks and inter-homologue repair events are now available for a number of mammalian, fungal and plant species. In mammals, PRDM9 specifies the location of meiotic recombination initiation via recognition of specific DNA sequence motifs by its C2H2 zinc finger array. In fungi and plants, meiotic recombination appears to be initiated less discriminately in accessible chromatin, including at gene promoters. Generally, meiotic crossover is suppressed in highly repetitive genomic regions that are made up of transposable elements (TEs), to prevent deleterious non-allelic homologous recombination events. However, recent and older studies have revealed intriguing relationships between meiotic recombination initiation and repair, and transposable elements. For instance, gene conversion events have been detected in maize centromeric retroelements, mouse MULE-MuDR DNA transposons undergo substantial meiotic recombination initiation, Arabidopsis Helitron TEs are among the hottest of recombination initiation hotspots, and human TE sequences can modify the crossover rate at adjacent PRDM9 motifs in cis. Here, we summarize the relationship between meiotic recombination and TEs, discuss recent insights from highly divergent eukaryotes and highlight outstanding questions in the field.

Regulation of transposable elements by DNA modifications

Maintenance of genome stability requires control over the expression of transposable elements (TEs), whose activity can have substantial deleterious effects on the host. Chemical modification of DNA is a commonly used strategy to achieve this, and it has long been argued that the emergence of 5-methylcytosine (5mC) in many species was driven by the requirement to silence TEs. Potential roles in TE regulation have also been suggested for other DNA modifications, such as N6-methyladenine and oxidation derivatives of 5mC, although the underlying mechanistic relationships are poorly understood. Here, we discuss current evidence implicating DNA modifications and DNA-modifying enzymes in TE regulation across different species.

Transposable Elements Adaptive Role in Genome Plasticity, Pathogenicity and Evolution in Fungal Phytopathogens

Transposable elements (TEs) are agents of genetic variability in phytopathogens as they are a source of adaptive evolution through genome diversification. Although many studies have uncovered information on TEs, the exact mechanism behind TE-induced changes within the genome remains poorly understood. Furthermore, convergent trends towards bigger genomes, emergence of novel genes and gain or loss of genes implicate a TE-regulated genome plasticity of fungal phytopathogens. TEs are able to alter gene expression by revamping the cis-regulatory elements or recruiting epigenetic control. Recent findings show that TEs recruit epigenetic control on the expression of effector genes as part of the coordinated infection strategy. In addition to genome plasticity and diversity, fungal pathogenicity is an area of economic concern. A survey of TE distribution suggests that their proximity to pathogenicity genes TEs may act as sites for emergence of novel pathogenicity factors via nucleotide changes and expansion or reduction of the gene family. Through a systematic survey of literature, we were able to conclude that the role of TEs in fungi is wide: ranging from genome plasticity, pathogenicity to adaptive behavior in evolution. This review also identifies the gaps in knowledge that requires further elucidation for a better understanding of TEs' contribution to genome architecture and versatility.

Genomic Takeover by Transposable Elements in the Strawberry Poison Frog

We sequenced the genome of the strawberry poison frog, Oophaga pumilio, at a depth of 127.5× using variable insert size libraries. The total genome size is estimated to be 6.76 Gb, of which 4.76 Gb are from high copy number repetitive elements with low differentiation across copies. These repeats encompass DNA transposons, RNA transposons, and LTR retrotransposons, including at least 0.4 and 1.0 Gb of Mariner/Tc1 and Gypsy elements, respectively. Expression data indicate high levels of gypsy and Mariner/Tc1 expression in ova of O. pumilio compared with Xenopus laevis. We further observe phylogenetic evidence for horizontal transfer (HT) of Mariner elements, possibly between fish and frogs. The elements affected by HT are present in high copy number and are highly expressed, suggesting ongoing proliferation after HT. Our results suggest that the large amphibian genome sizes, at least partially, can be explained by a process of repeated invasion of new transposable elements that are not yet suppressed in the germline. We also find changes in the spliceosome that we hypothesize are related to permissiveness of O. pumilio to increases in intron length due to transposon proliferation. Finally, we identify the complement of ion channels in the first genomic sequenced poison frog and discuss its relation to the evolution of autoresistance to toxins sequestered in the skin.

Capture of a functionally active methyl-CpG binding domain by an arthropod retrotransposon family

The repressive capacity of cytosine DNA methylation is mediated by recruitment of silencing complexes by methyl-CpG binding domain (MBD) proteins. Despite MBD proteins being associated with silencing, we discovered that a family of arthropod Copia retrotransposons have incorporated a host-derived MBD. We functionally show how retrotransposon-encoded MBDs preferentially bind to CpG-dense methylated regions, which correspond to transposable element regions of the host genome, in the myriapod Strigamia maritima Consistently, young MBD-encoding Copia retrotransposons (CopiaMBD) accumulate in regions with higher CpG densities than other LTR-retrotransposons also present in the genome. This would suggest that retrotransposons use MBDs to integrate into heterochromatic regions in Strigamia, avoiding potentially harmful insertions into host genes. In contrast, CopiaMBD insertions in the spider Stegodyphus dumicola genome disproportionately accumulate in methylated gene bodies compared with other spider LTR-retrotransposons. Given that transposons are not actively targeted by DNA methylation in the spider genome, this distribution bias would also support a role for MBDs in the integration process. Together, these data show that retrotransposons can co-opt host-derived epigenome readers, potentially harnessing the host epigenome landscape to advantageously tune the retrotransposition process.

Characterization of a novel Helitron family in insect genomes: insights into classification, evolution and horizontal transfer

Background

Helitrons play an important role in shaping eukaryotic genomes due to their ability to transfer horizontally between distantly related species and capture gene fragments during the transposition. However, the mechanisms of horizontal transfer (HT) and the process of gene fragment capturing of Helitrons still remain to be further clarified.

Results

Here, we characterized a novel Helitron family discontinuously distributed in 27 out of 256 insect genomes. The most prominent characteristic of Hel1 family is its high sequence similarity among species of different insect orders. Related elements were also identified in two spiders, representing the first report of spider Helitrons. All these elements were classified into 2 families, 9 subfamilies and 35 exemplars based on our new classification criteria. Autonomous partners of Helitron were reconstructed in the genomes of three insects and one spider. Integration pattern analysis showed that majority of Hel1A elements in Papilio xuthus and Pieris rapae inserted into introns. Consistent with filler DNA model, stepwise sequence acquisition was observed in Sfru_Hel1Aa, Sfru_Hel1Ab and Sfru_Hel1Ac in Spodoptera frugiperda. Remarkably, the evidence that Prap_Hel1Aa in a Lepdidoptera insect, Pieris rapae, was derived from Cves_Hel1Aa in a parasitoid wasp, Cotesia vestalis, suggested the role of nonregular host-parasite interactions in HT of Helitrons.

Conclusions

We proposed a modified classification criteria of Helitrons based on the important role of the 5′-end of Helitrons in transposition, and provided evidence for stepwise sequence acquisition and recurrent HT of a novel Helitron family. Our findings of the nonregular host-parasite interactions may be more conducive to the HT of transposons.

Electronic supplementary material

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

On the Population Dynamics of Junk: A Review on the Population Genomics of Transposable Elements

Transposable elements (TEs) play an important role in shaping genomic organization and structure, and may cause dramatic changes in phenotypes. Despite the genetic load they may impose on their host and their importance in microevolutionary processes such as adaptation and speciation, the number of population genetics studies focused on TEs has been rather limited so far compared to single nucleotide polymorphisms (SNPs). Here, we review the current knowledge about the dynamics of transposable elements at recent evolutionary time scales, and discuss the mechanisms that condition their abundance and frequency. We first discuss non-adaptive mechanisms such as purifying selection and the variable rates of transposition and elimination, and then focus on positive and balancing selection, to finally conclude on the potential role of TEs in causing genomic incompatibilities and eventually speciation. We also suggest possible ways to better model TEs dynamics in a population genomics context by incorporating recent advances in TEs into the rich information provided by SNPs about the demography, selection, and intrinsic properties of genomes.

Diversification and collapse of a telomere elongation mechanism

In most eukaryotes, telomerase counteracts chromosome erosion by adding repetitive sequence to terminal ends. Drosophila melanogaster instead relies on specialized retrotransposons that insert exclusively at telomeres. This exchange of goods between host and mobile element-wherein the mobile element provides an essential genome service and the host provides a hospitable niche for mobile element propagation-has been called a "genomic symbiosis." However, these telomere-specialized, jockey family retrotransposons may actually evolve to "selfishly" overreplicate in the genomes that they ostensibly serve. Under this model, we expect rapid diversification of telomere-specialized retrotransposon lineages and, possibly, the breakdown of this ostensibly symbiotic relationship. Here we report data consistent with both predictions. Searching the raw reads of the 15-Myr-old melanogaster species group, we generated de novo jockey retrotransposon consensus sequences and used phylogenetic tree-building to delineate four distinct telomere-associated lineages. Recurrent gains, losses, and replacements account for this retrotransposon lineage diversity. In Drosophila biarmipes, telomere-specialized elements have disappeared completely. De novo assembly of long reads and cytogenetics confirmed this species-specific collapse of retrotransposon-dependent telomere elongation. Instead, telomere-restricted satellite DNA and DNA transposon fragments occupy its terminal ends. We infer that D. biarmipes relies instead on a recombination-based mechanism conserved from yeast to flies to humans. Telomeric retrotransposon diversification and disappearance suggest that persistently "selfish" machinery shapes telomere elongation across Drosophila rather than completely domesticated, symbiotic mobile elements.

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.

Evolutionary origins and impacts of genome architecture in ciliates

Genome architecture is well diversified among eukaryotes in terms of size and content, with many being radically shaped by ancient and ongoing genome conflicts with transposable elements (e.g., the large transposon‐rich genomes common among plants). In ciliates, a group of microbial eukaryotes with distinct somatic and germ‐line genomes present in a single cell, the consequences of these genome conflicts are most apparent in their developmentally programmed genome rearrangements. This complicated developmental phenomenon has largely overshadowed and outpaced our understanding of how germ‐line and somatic genome architectures have influenced the evolutionary dynamism and potential in these taxa. In our review, we highlight three central concepts: how the evolution of atypical ciliate germ‐line genome architectures is linked to ancient genome conflicts; how the complex, epigenetically guided transformation of germline to soma during development can generate widespread genetic variation; and how these features, coupled with their unusual life cycle, have increased the rate of molecular evolution linked to genome architecture in these taxa.

Diversification of Transposable Elements in Arthropods and Its Impact on Genome Evolution

Transposable elements (TEs) are ubiquitous in arthropods. However, analyses of large-scale and long-term coevolution between TEs and host genomes remain scarce in arthropods. Here, we choose 14 representative Arthropoda species from eight orders spanning more than 500 million years of evolution. By developing an unbiased TE annotation pipeline, we obtained 87 to 2266 TE reference sequences in a species, which is a considerable improvement compared to the reference TEs previously annotated in Repbase. We find that TE loads are diversified among species and were previously underestimated. The highly species- and time-specific expansions and contractions, and intraspecific sequence diversification are the leading driver of long terminal repeat (LTR) dynamics in Lepidoptera. Terminal inverted repeats (TIRs) proliferated substantially in five species with large genomes. A phylogenetic comparison reveals that the loads of multiple TE subfamilies are positively correlated with genome sizes. We also identified a few horizontally transferred TE candidates across nine species. In addition, we set up the Arthropod Transposable Elements database (ArTEdb) to provide TE references and annotations. Collectively, our results provide high-quality TE references and uncover that TE loads and expansion histories vary greatly among arthropods, which implies that TEs are an important driving force shaping the evolution of genomes through gain and loss.

Transposon molecular domestication and the evolution of the RAG recombinase

Domestication of a transposon (a DNA sequence that can change its position in a genome) to give rise to the RAG1-RAG2 recombinase (RAG) and V(D)J recombination, which produces the diverse repertoire of antibodies and T cell receptors, was a pivotal event in the evolution of the adaptive immune system of jawed vertebrates. The evolutionary adaptations that transformed the ancestral RAG transposase into a RAG recombinase with appropriately regulated DNA cleavage and transposition activities are not understood. Here, beginning with cryo-electron microscopy structures of the amphioxus ProtoRAG transposase (an evolutionary relative of RAG), we identify amino acid residues and domains the acquisition or loss of which underpins the propensity of RAG for coupled cleavage, its preference for asymmetric DNA substrates and its inability to perform transposition in cells. In particular, we identify two adaptations specific to jawed-vertebrates-arginine 848 in RAG1 and an acidic region in RAG2-that together suppress RAG-mediated transposition more than 1,000-fold. Our findings reveal a two-tiered mechanism for the suppression of RAG-mediated transposition, illuminate the evolution of V(D)J recombination and provide insight into the principles that govern the molecular domestication of transposons.

Transcriptionally promiscuous "blurry" promoters in Tc1/mariner transposons allow transcription in distantly related genomes

Background

We have recently described a peculiar feature of the promoters in two Drosophila Tc1-like elements, Bari1 and Bari3. The AT-richness and the presence of weak core-promoter motifs make these promoters, that we have defined “blurry”, able to activate transcription of a reporter gene in cellular systems as diverse as fly, human, yeast and bacteria. In order to clarify whether the blurry promoter is a specific feature of the Bari transposon family, we have extended this study to promoters isolated from three additional DNA transposon and from two additional LTR retrotransposons.

Results

Here we show that the blurry promoter is also a feature of two vertebrate transposable elements, Sleeping Beauty and Hsmar1, belonging to the Tc1/mariner superfamily. In contrast, this feature is not shared by the promoter of the hobo transposon, which belongs to the hAT superfamily, nor by LTR retrotransposon-derived promoters, which, in general, do not activate transcription when introduced into non-related genomes.

Conclusions

Our results suggest that the blurry promoter could be a shared feature of the members of the Tc1/mariner superfamily with possible evolutionary and biotechnological implications.

Electronic supplementary material

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

Integrated mobile genetic elements in Thaumarchaeota

To explore the diversity of mobile genetic elements (MGE) associated with archaea of the phylum Thaumarchaeota, we exploited the property of most MGE to integrate into the genomes of their hosts. Integrated MGE (iMGE) were identified in 20 thaumarchaeal genomes amounting to 2 Mbp of mobile thaumarchaeal DNA. These iMGE group into five major classes: (i) proviruses, (ii) casposons, (iii) insertion sequence-like transposons, (iv) integrative-conjugative elements and (v) cryptic integrated elements. The majority of the iMGE belong to the latter category and might represent novel families of viruses or plasmids. The identified proviruses are related to tailed viruses of the order Caudovirales and to tailless icosahedral viruses with the double jelly-roll capsid proteins. The thaumarchaeal iMGE are all connected within a gene sharing network, highlighting pervasive gene exchange between MGE occupying the same ecological niche. The thaumarchaeal mobilome carries multiple auxiliary metabolic genes, including multicopper oxidases and ammonia monooxygenase subunit C (AmoC), and stress response genes, such as those for universal stress response proteins (UspA). Thus, iMGE might make important contributions to the fitness and adaptation of their hosts. We identified several iMGE carrying type I-B CRISPR-Cas systems and spacers matching other thaumarchaeal iMGE, suggesting antagonistic interactions between coexisting MGE and symbiotic relationships with the ir archaeal hosts.

Dynamic Interactions Between the Genome and an Endogenous Retrovirus: Tirant in Drosophila simulans Wild-Type Strains

All genomes contain repeated sequences that are known as transposable elements (TEs). Among these are endogenous retroviruses (ERVs), which are sequences similar to retroviruses and are transmitted across generations from parent to progeny. These sequences are controlled in genomes through epigenetic mechanisms. At the center of the epigenetic control of TEs are small interfering RNAs of the piRNA class, which trigger heterochromatinization of TE sequences. The tirant ERV of Drosophila simulans displays intra-specific variability in copy numbers, insertion sites, and transcription levels, providing us with a well-suited model to study the dynamic relationship between a TE family and the host genome through epigenetic mechanisms. We show that tirant transcript amounts and piRNA amounts are positively correlated in ovaries in normal conditions, unlike what was previously described following divergent crosses. In addition, we describe tirant insertion polymorphism in the genomes of three D. simulans wild-type strains, which reveals a limited number of insertions that may be associated with gene transcript level changes through heterochromatin spreading and have phenotypic impacts. Taken together, our results participate in the understanding of the equilibrium between the host genome and its TEs.

Hitchhiking on chromosomes: A persistence strategy shared by diverse selfish DNA elements

An underlying theme in the segregation of low-copy bacterial plasmids is the assembly of a 'segrosome' by DNA-protein and protein-protein interactions, followed by energy-driven directed movement. Analogous partitioning mechanisms drive the segregation of host chromosomes as well. Eukaryotic extra-chromosomal elements, exemplified by budding yeast plasmids and episomes of certain mammalian viruses, harbor partitioning systems that promote their physical association with chromosomes. In doing so, they indirectly take advantage of the spindle force that directs chromosome movement to opposite cell poles. Molecular-genetic, biochemical and cell biological studies have revealed several unsuspected aspects of 'chromosome hitchhiking' by the yeast 2-micron plasmid, including the ability of plasmid sisters to associate symmetrically with sister chromatids. As a result, the plasmid overcomes the 'mother bias' experienced by plasmids lacking a partitioning system, and elevates itself to near chromosome status in equal segregation. Chromosome association for stable propagation, without direct energy expenditure, may also be utilized by a small minority of bacterial plasmids-at least one case has been reported. Given the near perfect accuracy of chromosome segregation, it is not surprising that elements residing in evolutionarily distant host organisms have converged upon the common strategy of gaining passage to daughter cells as passengers on chromosomes.

Retrotranspositional landscape of Asian rice revealed by 3000 genomes

The recent release of genomic sequences for 3000 rice varieties provides access to the genetic diversity at species level for this crop. We take advantage of this resource to unravel some features of the retrotranspositional landscape of rice. We develop software TRACKPOSON specifically for the detection of transposable elements insertion polymorphisms (TIPs) from large datasets. We apply this tool to 32 families of retrotransposons and identify more than 50,000 TIPs in the 3000 rice genomes. Most polymorphisms are found at very low frequency, suggesting that they may have occurred recently in agro. A genome-wide association study shows that these activations in rice may be triggered by external stimuli, rather than by the alteration of genetic factors involved in transposable element silencing pathways. Finally, the TIPs dataset is used to trace the origin of rice domestication. Our results suggest that rice originated from three distinct domestication events.

Ten things you should know about transposable elements

Transposable elements (TEs) are major components of eukaryotic genomes. However, the extent of their impact on genome evolution, function, and disease remain a matter of intense interrogation. The rise of genomics and large-scale functional assays has shed new light on the multi-faceted activities of TEs and implies that they should no longer be marginalized. Here, we introduce the fundamental properties of TEs and their complex interactions with their cellular environment, which are crucial to understanding their impact and manifold consequences for organismal biology. While we draw examples primarily from mammalian systems, the core concepts outlined here are relevant to a broad range of organisms.

Firmicutes-enriched IS1447 represents a group of IS3-family insertion sequences exhibiting unique + 1 transcriptional slippage

BACKGROUND

Bacterial insertion sequences (ISs) are ubiquitous mobile genetic elements that play important roles in genome plasticity, cell adaptability, and function evolution. ISs of various families and subgroups contain significantly diverse molecular features and functional mechanisms that are not fully understood.

RESULTS

IS1447 is a member of the widespread IS3 family and was previously detected to have transposing activity in a typical thermophilic and cellulolytic microorganism Clostridium thermocellum. Phylogenetic analysis showed that IS1447-like elements are widely distributed in Firmicutes and possess unique features in the IS3 family. Therefore, IS1447 may represent a novel subgroup of the IS3 family. Unlike other well-known IS3 subgroups performing programmed - 1 translational frameshifting for the expression of the transposase, IS1447 exhibits transcriptional slippage in both the + 1 and - 1 directions, each with a frequency of ~ 16%, and only + 1 slippage results in full-length and functional transposase. The slippage-prone region of IS1447 contains a run of nine A nucleotides following a stem-loop structure in mRNA, but mutagenesis analysis indicated that seven of them are sufficient for the observed slippage. Western blot analysis indicated that IS1447 produces three types of transposases with alternative initiations. Furthermore, the IS1447-subgroup elements are abundant in the genomes of several cellulolytic bacteria.

CONCLUSION

Our result indicated that IS1447 represents a new Firmicutes-enriched subgroup of the IS3 family. The characterization of the novel IS3-family member will enrich our understanding of the transposition behavior of IS elements and may provide insight into developing IS-based mutagenesis tools for thermophiles.

Six domesticated PiggyBac transposases together carry out programmed DNA elimination in Paramecium

The domestication of transposable elements has repeatedly occurred during evolution and domesticated transposases have often been implicated in programmed genome rearrangements, as remarkably illustrated in ciliates. In Paramecium, PiggyMac (Pgm), a domesticated PiggyBac transposase, carries out developmentally programmed DNA elimination, including the precise excision of tens of thousands of gene-interrupting germline Internal Eliminated Sequences (IESs). Here, we report the discovery of five groups of distant Pgm-like proteins (PgmLs), all able to interact with Pgm and essential for its nuclear localization and IES excision genome-wide. Unlike Pgm, PgmLs lack a conserved catalytic site, suggesting that they rather have an architectural function within a multi-component excision complex embedding Pgm. PgmL depletion can increase erroneous targeting of residual Pgm-mediated DNA cleavage, indicating that PgmLs contribute to accurately position the complex on IES ends. DNA rearrangements in Paramecium constitute a rare example of a biological process jointly managed by six distinct domesticated transposases.

The impact of transposable elements in adaptive evolution

The growing knowledge about the influence of transposable elements (TEs) on (a) long-term genome and transcriptome evolution; (b) genomic, transcriptomic and epigenetic variation within populations; and (c) patterns of somatic genetic differences in individuals continues to spur the interest of evolutionary biologists in the role of TEs in adaptive evolution. As TEs can trigger a broad range of molecular variation in a population with potentially severe fitness and phenotypic consequences for individuals, different mechanisms evolved to keep TE activity in check, allowing for a dynamic interplay between the host, its TEs and the environment in evolution. Here, we review evidence for adaptive phenotypic changes associated with TEs and the basic molecular mechanisms by which the underlying genetic changes arise: (a) domestication, (b) exaptation, (c) host gene regulation, (d) TE-mediated formation of intronless gene copies-so-called retrogenes and (e) overall increased genome plasticity. Furthermore, we review and discuss how the stress-dependent incapacitation of defence mechanisms against the activity of TEs might facilitate adaptive responses to environmental challenges and how such mechanisms might be particularly relevant in species frequently facing novel environments, such as invasive, pathogenic or parasitic species.

Border collies of the genome: domestication of an autonomous retrovirus-like transposon

Retrotransposons often spread rapidly through eukaryotic genomes until they are neutralized by host-mediated silencing mechanisms, reduced by recombination and mutation, and lost or transformed into benevolent entities. But the Ty1 retrotransposon appears to have been domesticated to guard the genome of Saccharomyces cerevisiae.

Multiple massive domestication and recent amplification of Kolobok superfamily transposons in the clawed frog Xenopus

BACKGROUND

DNA transposons are generally destroyed by mutations and have short lifespans in hosts, as they are neutral or harmful to the host and therefore not conserved by natural selection. The clawed frog Xenopus harbors many DNA transposons and certain families, such as T2-MITE, have extremely long lives. These have ancient origins, but have shown recent transposition activity. In addition, certain transposase genes may have been "domesticated" by Xenopus and conserved over long time periods by natural selection. The aim of this study was to elucidate the evolutionary interactions between the host and the long-lived DNA transposon family it contains. Here, we investigated the molecular evolution of the Kolobok DNA transposon superfamily. Kolobok is thought to contribute to T2-MITE transposition.

RESULTS

In the diploid western clawed frog Xenopus tropicalis and the allotetraploid African clawed frog Xenopus laevis, we searched for transposase genes homologous to those in the Kolobok superfamily. To determine the amplification and domestication of these genes, we used molecular phylogenetics and analyses of copy numbers, conserved motifs, orthologous gene synteny, and coding sequence divergence between the orthologs of X. laevis and X. tropicalis, or between those of two distant X. tropicalis lineages. Among 38 X. tropicalis and 24 X. laevis prospective transposase genes, 10 or more in X. tropicalis and 14 or more in X. laevis were apparently domesticated. These genes may have undergone multiple independent domestications from before the divergence of X. laevis and X. tropicalis. In contrast, certain other transposases may have retained catalytic activity required for transposition and could therefore have been recently amplified.

CONCLUSION

Multiple domestication of certain transposases and prolonged conservation of the catalytic activity in others suggest that Kolobok superfamily transposons were involved in complex, mutually beneficial relationships with their Xenopus hosts. Some transposases may serve to activate long-lived T2-MITE subfamilies.

Recurrent acquisition of cytosine methyltransferases into eukaryotic retrotransposons

Transposable elements are in a constant arms race with the silencing mechanisms of their host genomes. One silencing mechanism commonly used by many eukaryotes is dependent on cytosine methylation, a covalent modification of DNA deposited by C5 cytosine methyltransferases (DNMTs). Here, we report how two distantly related eukaryotic lineages, dinoflagellates and charophytes, have independently incorporated DNMTs into the coding regions of distinct retrotransposon classes. Concomitantly, we show that dinoflagellates of the genus Symbiodinium have evolved cytosine methylation patterns unlike any other eukaryote, with most of the genome methylated at CG dinucleotides. Finally, we demonstrate the ability of retrotransposon DNMTs to methylate CGs de novo, suggesting that retrotransposons could self-methylate retrotranscribed DNA. Together, this is an example of how retrotransposons incorporate host-derived genes involved in DNA methylation. In some cases, this event could have implications for the composition and regulation of the host epigenomic environment.

Integrative analysis of large scale transcriptome data draws a comprehensive landscape of Phaeodactylum tricornutum genome and evolutionary origin of diatoms

Diatoms are one of the most successful and ecologically important groups of eukaryotic phytoplankton in the modern ocean. Deciphering their genomes is a key step towards better understanding of their biological innovations, evolutionary origins, and ecological underpinnings. Here, we have used 90 RNA-Seq datasets from different growth conditions combined with published expressed sequence tags and protein sequences from multiple taxa to explore the genome of the model diatom Phaeodactylum tricornutum, and introduce 1,489 novel genes. The new annotation additionally permitted the discovery of extensive alternative splicing in diatoms, including intron retention and exon skipping, which increase the diversity of transcripts generated in changing environments. In addition, we have used up-to-date reference sequence libraries to dissect the taxonomic origins of diatom genes. We show that the P. tricornutum genome is enriched in lineage-specific genes, with up to 47% of the gene models present only possessing orthologues in other stramenopile groups. Finally, we have performed a comprehensive de novo annotation of repetitive elements showing novel classes of transposable elements such as SINE, MITE and TRIM/LARD. This work provides a solid foundation for future studies of diatom gene function, evolution and ecology.

Exaptation of transposable element coding sequences

Transposable elements (TEs) are mobile genetic elements that were once perceived as merely selfish, but are now recognized as potent agents of adaptation. One way TEs contribute to genome evolution is through TE exaptation, a process whereby TEs, which usually persist by replicating in the genome, transform into novel host genes, which thereafter persist by conferring phenotypic benefits. Exapted TEs are known to contribute diverse and vital functions, and may facilitate punctuated equilibrium, yet we have little understanding about the process of TE exaptation. In order to facilitate our understanding of how TE coding sequences may become exapted, here we incorporate the findings of recent publications into a framework and six-step model.

Horizontal acquisition of transposable elements and viral sequences: patterns and consequences

It is becoming clear that most eukaryotic transposable elements (TEs) owe their evolutionary success in part to horizontal transfer events, which enable them to invade new species. Recent large-scale studies are beginning to unravel the mechanisms and ecological factors underlying this mode of transmission. Viruses are increasingly recognized as vectors in the process but also as a direct source of genetic material horizontally acquired by eukaryotic organisms. Because TEs and endogenous viruses are major catalysts of variation and innovation in genomes, we argue that horizontal inheritance has had a more profound impact in eukaryotic evolution than is commonly appreciated. To support this proposal, we compile a list of examples, including some previously unrecognized, whereby new host functions and phenotypes can be directly attributed to horizontally acquired TE or viral sequences. We predict that the number of examples will rapidly grow in the future as the prevalence of horizontal transfer in the life cycle of TEs becomes even more apparent, firmly establishing this form of non-Mendelian inheritance as a consequential facet of eukaryotic evolution.

Circular RNAs Biogenesis in Eukaryotes Through Self-Cleaving Hammerhead Ribozymes

Circular DNAs are frequent genomic molecules, especially among the simplest life beings, whereas circular RNAs have been regarded as weird nucleic acids in biology. Now we know that eukaryotes are able to express circRNAs, mostly derived from backsplicing mechanisms, and playing different biological roles such as regulation of RNA splicing and transcription, among others. However, a second natural and highly efficient pathway for the expression in vivo of circRNAs has been recently reported, which allows the accumulation of abundant small (100-1000 nt) non-coding RNA circles through the participation of small self-cleaving RNAs or ribozymes called hammerhead ribozymes. These genome-encoded circRNAs with ribozymes seem to be a new family of small and nonautonomous retrotransposable elements of plants and animals (so-called retrozymes), which will offer functional clues to the biology and evolution of circular RNA molecules as well as new biotechnological tools in this emerging field.