Target Site Duplications Research Articles

How Well Does Evolution Explain Endogenous Retroviruses?-A Lakatosian Assessment.

One of the most sophisticated philosophies of science is the methodology of scientific research programmes (MSRP), developed by Imre Lakatos. According to MSRP, scientists are working within so-called research programmes, consisting of a hard core of fixed convictions and a flexible protective belt of auxiliary hypotheses. Anomalies are accommodated by changes to the protective belt that do not affect the hard core. Under MSRP, research programmes are appraised as ‘progressive’ if they successfully predict novel facts but are judged as ‘degenerative’ if they merely offer ad hoc solutions to anomalies. This paper applies these criteria to the evolutionary research programme as it has performed during half a century of ERV research. It describes the early history of the field and the emergence of the endogenization-amplification theory on the origins of retroviral-like sequences. It then discusses various predictions and postdictions that were generated by the programme, regarding orthologous ERVs in different species, the presence of target site duplications and the divergence of long terminal repeats, and appraises how the programme has dealt with data that did not conform to initial expectations. It is concluded that the evolutionary research programme has been progressive with regard to the issues here examined.

The genome of the diatom Chaetoceros tenuissimus carries an ancient integrated fragment of an extant virus

Diatoms are one of the most prominent oceanic primary producers and are now recognized to be distributed throughout the world. They maintain their population despite predators, infections, and unfavourable environmental conditions. One of the smallest diatoms, Chaetoceros tenuissimus, can coexist with infectious viruses during blooms. To further understand this relationship, we sequenced the C. tenuissimus strain NIES-3715 genome. A gene fragment of a replication-associated gene from the infectious ssDNA virus (designated endogenous virus-like fragment, EVLF) was found to be integrated into each 41 Mb of haploid assembly. In addition, the EVLF was transcriptionally active and conserved in nine other C. tenuissimus strains from different geographical areas, although the primary structures of their proteins varied. The phylogenetic tree further suggested that the EVLF was acquired by the ancestor of C. tenuissimus. Additionally, retrotransposon genes possessing a reverse transcriptase function were more abundant in C. tenuissimus than in Thalassiosira pseudonana and Phaeodactylum tricornutum. Moreover, a target site duplication, a hallmark for long interspersed nuclear element retrotransposons, flanked the EVLF. Therefore, the EVLF was likely integrated by a retrotransposon during viral infection. The present study provides further insights into the diatom-virus evolutionary relationship.

Prokaryotic and Eukaryotic Horizontal Transfer of Sailor (DD82E), a New Superfamily of IS630-Tc1-Mariner DNA Transposons.

Simple SummaryTransposable elements, including DNA transposons, play a significant role in genetic material exchanges between prokaryotes and eukaryotes. Comparative profiling of the evolution pattern of DNA transposons between prokaryotes and eukaryotes may identify potential genetic material exchanges between them and provide insights into the evolutionary history of prokaryotic and eukaryotic genomes. The members of the IS630-Tc1-mariner (ITm) group may represent the most diverse and widely distributed DNA transposons in nature, and the discovery of new members of this group is highly expected based on the increasing availability of genome sequencing data. We discovered a new superfamily (termed Sailor) belonging to the ITm hyperfamily, which differed from the known superfamilies of Tc1/mariner, DDxD/pogo and DD34E/Gambol, regarding phylogenetic position and catalytic domain. Our data revealed that Sailor was distributed in both prokaryotes and eukaryotes and suggested that horizontal transfer (HT) events of Sailor may occur from prokaryotic to eukaryotic genomes. Finally, internal transmissions of Sailor in prokaryotes and eukaryotes were also detected.Here, a new superfamily of IS630-Tc1-mariner (ITm) DNA transposons, termed Sailor, is identified, that is characterized by a DD82E catalytic domain and is distinct from all previously known superfamilies of the ITm group. Phylogenetic analyses revealed that Sailor forms a monophyletic clade with a more intimate link to the clades of Tc1/mariner and DD34E/Gambol. Sailor was detected in both prokaryotes and eukaryotes and invaded a total of 256 species across six kingdoms. Sailor is present in nine species of bacteria, two species of plantae, four species of protozoa, 23 species of Chromista, 12 species of Fungi and 206 species of animals. Moreover, Sailor is extensively distributed in invertebrates (a total of 206 species from six phyla) but is absent in vertebrates. Sailor transposons are 1.38–6.98 kb in total length and encoded transposases of ~676 aa flanked by TIRs with lengths between 18, 1362 and 4 bp (TATA) target-site duplications. Furthermore, our analysis provided strong evidence of Sailor transmissions from prokaryotes to eukaryotes and internal transmissions in both. These data update the classification of the ITm group and will contribute to the understanding of the evolution of ITm transposons and that of their hosts.

The function of twister ribozyme variants in non-LTR retrotransposition in Schistosoma mansoni.

The twister ribozyme is widely distributed over numerous organisms and is especially abundant in Schistosoma mansoni, but has no confirmed biological function. Of the 17 non-LTR retrotransposons known in S. mansoni, none have thus far been associated with ribozymes. Here we report the identification of novel twister variant (T-variant) ribozymes and their function in S. mansoni non-LTR retrotransposition. We show that T-variant ribozymes are located at the 5′ end of Perere-3 non-LTR retrotransposons in the S. mansoni genome. T-variant ribozymes were demonstrated to be catalytically active in vitro. In reporter constructs, T-variants were shown to cleave in vivo, and cleavage of T-variants was sufficient for the translation of downstream reporter genes. Our analysis shows that the T-variants and Perere-3 are transcribed together. Target site duplications (TSDs); markers of target-primed reverse transcription (TPRT) and footmarks of retrotransposition, are located adjacent to the T-variant cleavage site and suggest that T-variant cleavage has taken place inS. mansoni. Sequence heterogeneity in the TSDs indicates that Perere-3 retrotransposition is not site-specific. The TSD sequences contribute to the 5′ end of the terminal ribozyme helix (P1 stem). Based on these results we conclude that T-variants have a functional role in Perere-3 retrotransposition.

Sequence specific integration by the family 1 casposase from CandidatusNitrosopumilus koreensis AR1

Casposase, a homolog of Cas1 integrase, is encoded by a superfamily of mobile genetic elements known as casposons. While family 2 casposase has been well documented in both function and structure, little is known about the other three casposase families. Here, we studied the family 1 casposase lacking the helix-turn-helix (HTH) domain from Candidatus Nitrosopumilus koreensis AR1 (Ca. N. koreensis). The determinants for integration by Ca. N. koreensis casposase were extensively investigated, and it was found that a 13-bp target site duplication (TSD) sequence, a minimal 3-bp leader and three different nucleotides of the TSD sequences are indispensable for target specific integration. Significantly, the casposase can site-specifically integrate a broad range of terminal inverted repeat (TIR)-derived oligonucleotides ranging from 7-nt to ∼4000-bp, and various oligonucleotides lacking the 5′-TTCTA-3′ motif at the 3′ end of TIR sequence can be integrated efficiently. Furthermore, similar to some Cas1 homologs, the casposase utilizes a 5′-ATAA-3′ motif in the TSD as a molecular ruler to dictate nucleophilic attack at 9-bp downstream of the end of the ruler during the spacer-side integration. By characterizing the family 1 Ca. N. koreensis casposase, we have extended our understanding on mechanistic similarities and evolutionary connections between casposons and the adaptation elements of CRISPR-Cas immunity.

In and out: Evolution of viral sequences in the mitochondrial genomes of legumes (Fabaceae)

Plant specific mitoviruses in the ‘genus’ Mitovirus (Narnaviridae) and their integrated sequences (non-retroviral endogenous RNA viral elements or NERVEs) have been recently identified in various plant lineages. However, the sparse phylogenetic coverage of complete plant mitochondrial genome (mitogenome) sequences and the non-conserved nature of mitochondrial intergenic regions have hindered comparative studies on mitovirus NERVEs in plants. In this study, 10 new mitogenomes were sequenced from legumes (Fabaceae). Based on comparative genomic analysis of 27 total mitogenomes, we identified mitovirus NERVEs and transposable elements across the family. All legume mitogenomes included NERVEs and total NERVE length varied from ca. 2 kb in the papilionoid Trifolium to 35 kb in the mimosoid Acacia. Most of the NERVE integration sites were in highly variable intergenic regions, however, some were positioned in six cis-spliced mitochondrial introns. In the Acacia mitogenome, there were L1-like transposon sequences including an almost full-length copy with target site duplications (TSDs). The integration sites of NERVEs in four introns showed evidence of L1-like retrotransposition events. Phylogenetic analysis revealed that there were multiple instances of precise deletion of NERVEs between TSDs. This study provides clear evidence that a L1-like retrotransposition mechanism has a long history of contributing to the integration of viral RNA into plant mitogenomes while microhomology-mediated deletion can restore the integration site.

Nuku, a family of primate retrocopies derived from KU70.

The gene encoding the ubiquitous DNA repair protein, Ku70p, has undergone extensive copy number expansion during primate evolution. Gene duplications of KU70 have the hallmark of long interspersed element-1 mediated retrotransposition with evidence of target-site duplications, the poly-A tails, and the absence of introns. Evolutionary analysis of this expanded family of KU70-derived “NUKU” retrocopies reveals that these genes are both ancient and also actively being created in extant primate species. NUKU retrocopies show evidence of functional divergence away from KU70, as evinced by their altered pattern of tissue expression and possible tissue-specific translation. Molecular modeling predicted that amino acid changes in Nuku2p at the interaction interface with Ku80p would prevent the assembly of the Ku heterodimer. The lack of Nuku2p-Ku80p interaction was confirmed by yeast two-hybrid assay, which contrasts the robust interaction of Ku70p-Ku80p. While several NUKU retrocopies appear to have been degraded by mutation, NUKU2 shows evidence of positive natural selection, suggesting that this retrocopy is undergoing neofunctionalization. Although Nuku proteins do not appear to antagonize retrovirus transduction in cell culture, the observed expansion and rapid evolution of NUKUs could be being driven by alternative selective pressures related to infectious disease or an undefined role in primate physiology.

Reverse-transcribed SARS-CoV-2 RNA can integrate into the genome of cultured human cells and can be expressed in patient-derived tissues

Prolonged detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA and recurrence of PCR-positive tests have been widely reported in patients after recovery from COVID-19, but some of these patients do not appear to shed infectious virus. We investigated the possibility that SARS-CoV-2 RNAs can be reverse-transcribed and integrated into the DNA of human cells in culture and that transcription of the integrated sequences might account for some of the positive PCR tests seen in patients. In support of this hypothesis, we found that DNA copies of SARS-CoV-2 sequences can be integrated into the genome of infected human cells. We found target site duplications flanking the viral sequences and consensus LINE1 endonuclease recognition sequences at the integration sites, consistent with a LINE1 retrotransposon-mediated, target-primed reverse transcription and retroposition mechanism. We also found, in some patient-derived tissues, evidence suggesting that a large fraction of the viral sequences is transcribed from integrated DNA copies of viral sequences, generating viral-host chimeric transcripts. The integration and transcription of viral sequences may thus contribute to the detection of viral RNA by PCR in patients after infection and clinical recovery. Because we have detected only subgenomic sequences derived mainly from the 3' end of the viral genome integrated into the DNA of the host cell, infectious virus cannot be produced from the integrated subgenomic SARS-CoV-2 sequences.

A Cas12a-based gene editing system for Phytophthora infestans reveals monoallelic expression of an elicitor.

Phytophthora infestans is a destructive pathogen of potato and a model for investigations of oomycete biology. The successful application of a CRISPR gene editing system to P. infestans is so far unreported. We discovered that it is difficult to express CRISPR/Cas9 but not a catalytically inactive form in transformants, suggesting that the active nuclease is toxic. We were able to achieve editing with CRISPR/Cas12a using vectors in which the nuclease and its guide RNA were expressed from a single transcript. Using the elicitor gene Inf1 as a target, we observed editing of one or both alleles in up to 13% of transformants. Editing was more efficient when guide RNA processing relied on the Cas12a direct repeat instead of ribozyme sequences. INF1 protein was not made when both alleles were edited in the same transformant, but surprisingly also when only one allele was altered. We discovered that the isolate used for editing, 1306, exhibited monoallelic expression of Inf1 due to insertion of a copia‐like element in the promoter of one allele. The element exhibits features of active retrotransposons, including a target site duplication, long terminal repeats, and an intact polyprotein reading frame. Editing occurred more often on the transcribed allele, presumably due to differences in chromatin structure. The Cas12a system not only provides a tool for modifying genes in P. infestans, but also for other members of the genus by expanding the number of editable sites. Our work also highlights a natural mechanism that remodels oomycete genomes.

Divergent evolution profiles of DD37D and DD39D families of Tc1/mariner transposons in eukaryotes

DNA transposons play a significant role in shaping the size and structure of eukaryotic genomes. The Tc1/mariner transposons are the most diverse and widely distributed superfamily of DNA transposons and the structure and distribution of several Tc1/mariner families, such as DD35E/TR, DD36E/IC, DD37E/TRT, and DD41D/VS, have been well studied. Nonetheless, a greater understanding of the structure and diversity of Tc1/mariner transposons will provide insight into the evolutionary history of eukaryotic genomes. Here, we conducted further analysis of DD37D/maT and DD39D (named Guest, GT), which were identified by the specific catalytic domains DD37D and DD39D. Most transposons of the maT family have a total length of approximately 1.3 kb and harbor a single open reading frame encoding a ~ 346 amino acid (range 302–398 aa) transposase protein, flanked by short terminal inverted repeats (TIRs) (13–48 base pairs, bp). In contrast, GTs transposons were longer (2.0–5.8 kb), encoded a transposase protein of ~400 aa (range 140–592 aa), and were flanked by short TIRs (19–41 bp). Several conserved motifs, including two helix–turn–helix (HTH) motifs, a GRPR (GRKR) motif, a nuclear localization sequence, and a DDD domain, were also identified in maT and GT transposases. Phylogenetic analyses of the DDD domain showed that the maT and GT families each belong to a monophyletic clade and appear to be closely related to DD41D/VS and DD34D/mariner. In addition, maTs are mainly distributed in invertebrates (144 species), whereas GTs are mainly distributed in land plants through a small number of GTs are present in Chromista and animals. Sequence identity and phylogenetic analysis revealed that horizontal transfer (HT) events of maT and GT might occur between kingdoms and phyla of eukaryotes; however, pairwise distance comparisons between host genes and transposons indicated that HT events involving maTs might be less frequent between invertebrate species and HT events involving GTs may be less frequent between land plant species. Overall, the DD37D/maT and DD39D/GT families display significantly different distribution and tend to be identified in more ancient evolutionary families. The discovery of intact transposases, perfect TIRs, and target site duplications (TSD) of maTs and GTs illustrates that the DD37D/maT and DD39D/GT families may be active. Together, these findings improve our understanding of the diversity of Tc1/mariner transposons and their impact on eukaryotic genome evolution.

Dissemination of novel Tn7 family transposons carrying genes for synthesis and uptake of fimsbactin siderophores among Acinetobacter baumannii isolates.

Acinetobacter baumannii is a successful opportunistic pathogen that can compete for iron under iron-limiting conditions. Here, large novel transposons that carry genes for synthesis and transport of the fimsbactin siderophores present in some A. baumannii strains were examined. Tn6171, originally found in the A. baumannii global clone 1 (GC1) lineage 2 isolate D36, includes tns genes encoding proteins related to the TnsA, TnsB, TnsC transposition proteins (50–59 % identity), TnsD targeting protein (43 % identity) and TnsE (31 % identity) of Tn7, and is found in the chromosome downstream of the glmS gene, the preferred location for Tn7, flanked by a 5 bp target site duplication. Tn6171 is bounded by 29 bp inverted repeats and, like Tn7, includes additional TnsB binding sites at each end. Tn6171 or minor variants were detected in the equivalent location in complete or draft genomes of several further A. baumannii isolates belonging to GC1 [sequence type (ST) 1, ST81, ST94, ST328, ST623, ST717], GC2 (ST2) and ST10. However, in some of these isolates the surrounding glmS region was clearly derived from a different A. baumannii lineage, indicating that the transposon may have been acquired by replacement of a segment of the chromosome. A recombination-free phylogeny revealed that there were several transposon acquisition events in GC1. The GC1 isolates were mainly lineage 2, but a potential third lineage was also detected. A related transposon, designated Tn6552, was detected in ATCC 17978 (ST437) and other ST437 isolates. However, the Tn6552 tnsD targeting gene was interrupted by an ISAba12, and Tn6552 is not downstream of glmS.

Detection of Transposition Events from Next-Generation Sequencing Data.

Detection of transposition events of a transposon from short reads of next-generation sequencing (NGS) is challenging because transposons are repetitive and difficult to be distinguished from already existing transposons in the genome. Many transposons generate target site duplication (TSD) as the result of chromosomal integration. Since TSDs flanking the 5'-end (head) and 3'-end (tail) of a transposon has the identical sequences which are absent from the reference copy, the short reads containing the head or tail sequences of the transposon following the same TSD sequence may reveal the evidence of transposition. Transposon Insertion Finder (TIF) focuses on the TSD with flanking sequence of transposon and detects transposition events from NGS data. TIF software is available at https://github.com/akiomiyao/tif .

Growth by Insertion: The Family of Bacterial DDxP Proteins.

We have identified a variety of proteins in species of the Legionella, Aeromonas, Pseudomonas, Vibrio, Nitrosomonas, Nitrosospira, Variovorax, Halomonas, and Rhizobia genera, which feature repetitive modules of different length and composition, invariably ending at the COOH side with Asp–Asp–x–Pro (DDxP) motifs. DDxP proteins range in size from 900 to 6200 aa (amino acids), and contain 1 to 5 different module types, present in one or multiple copies. We hypothesize that DDxP proteins were modeled by the action of specific endonucleases inserting DNA segments into genes encoding DDxP motifs. Target site duplications (TSDs) formed upon repair of staggered ends generated by endonuclease cleavage would explain the DDxP motifs at repeat ends. TSDs acted eventually as targets for the insertion of more modules of the same or different types. Repeat clusters plausibly resulted from amplification of both repeat and flanking TSDs. The proposed growth shown by the insertion model is supported by the identification of homologous proteins lacking repeats in Pseudomonas and Rhizobium. The 85 DDxP repeats identified in this work vary in length, and can be sorted into short (136–215 aa) and long (243–304 aa) types. Conserved Asp–Gly–Asp–Gly–Asp motifs are located 11–19 aa from the terminal DDxP motifs in all repeats, and far upstream in most long repeats.

Novel insertion sequence ISChh1-like mediating acquisition of optrA gene in foodborne pathogen Campylobacter coli of swine origin

The optrA gene encodes an ABC-F protein which confers cross-resistance to oxazolidinones and phenicols, and its transmission has so far been associated with multiple transposable elements, including IS1216E, prophages, ICEs and Tn558. Here, we identified an optrA gene flanked by two copies of a novel insertion sequence ISChh1-like in the same direction in Campylobacter coli. Seven optrA-positive C. coli were identified from 263 Campylobacter isolates obtained from one swine farm and two slaughterhouses in Sichuan province of China. The optrA genes in 6 isolates were directly flanked by two copies of ISChh1-like elements in the same orientation, in the remaining one isolate, optrA was co-located with fexA and flanked by two copies of IS1216E and inserted in the downstream of the aadE-sat4-aphA3 gene cluster on chromosome. Cloning of optrA into C. jejuni NCTC 11,168 confirmed its role in elevated MICs to oxazolidinones and phenicols. Translocatable units (TUs) and natural transformants were only detected from the isolate with optrA bracketed by IS1216E, not from these with optrA bracketed by ISChh1-like. ISChh1-like in all isolates inserts specifically either next to or between adjacent GG nucleotides, neither have terminal inverted repeats nor generate target site duplications. Interestingly, ISChh1-like were also found mediating integration of other antibiotic resistance genes in Campylobacter, including tet(O), aphA3 and aadE-sat4-aphA3 gene cluster. Taken together, these results identify ISChh1-like as a novel transposon mediating acquisition of multiple antibiotic resistance genes in Campylobacter, including the very important optrA gene, suggesting that it plays an essential role in the transmission of antibiotic resistance genes to Campylobacter.

Pooled Plasmid Sequencing Reveals the Relationship Between Mobile Genetic Elements and Antimicrobial Resistance Genes in Clinically Isolated Klebsiella pneumoniae

Plasmids remain important microbial components mediating the horizontal gene transfer (HGT) and dissemination of antimicrobial resistance. To systematically explore the relationship between mobile genetic elements (MGEs) and antimicrobial resistance genes (ARGs), a novel strategy using single-molecule real-time (SMRT) sequencing was developed. This approach was applied to pooled conjugative plasmids from clinically isolated multidrug-resistant (MDR) Klebsiella pneumoniae from a tertiary referral hospital over a 9-month period. The conjugative plasmid pool was obtained from transconjugants that acquired antimicrobial resistance after plasmid conjugation with 53 clinical isolates. The plasmid pool was then subjected to SMRT sequencing, and 82 assembled plasmid fragments were obtained. In total, 124 ARGs (responsible for resistance to β-lactam, fluoroquinolone, and aminoglycoside, among others) and 317 MGEs [including transposons (Tns), insertion sequences (ISs), and integrons] were derived from these fragments. Most of these ARGs were linked to MGEs, allowing for the establishment of a relationship network between MGEs and/or ARGs that can be used to describe the dissemination of resistance by mobile elements. Key elements involved in resistance transposition were identified, including IS26, Tn3, IS903B, ISEcp1, and ISKpn19. As the most predominant IS in the network, a typical IS26-mediated multicopy composite transposition event was illustrated by tracing its flanking 8-bp target site duplications (TSDs). The landscape of the pooled plasmid sequences highlights the diversity and complexity of the relationship between MGEs and ARGs, underpinning the clinical value of dominant HGT profiles.

Multiple insertions of COIN, a novel maize Foldback transposable element, in the Conring gene cause a spontaneous progressive cell death phenotype.

Similar progressive leaf lesion phenotypes, named conring for "concentric ring," were identified in 10 independently derived maize lines. Complementation and mapping experiments indicated that the phenotype had the same genetic basis in each line - a single recessive gene located in a 1.1-Mb region on chromosome 2. Among the 15 predicted genes in this interval, Zm00001d003866 (subsequently renamed Conring or Cnr) had insertions of four related 138bp transposable element (TE) sequences at precisely the same site in exon 4 in nine of the 10 cnr alleles. The 10th cnr allele had a distinct insertion of 226bp of in exon 3. Genetic evidence suggested that the 10 cnr alleles were independently derived, and arose during the derivation of each line. The four TEs, named COINa (for COnring INsertion) through COINd, have not been previously characterized and consist entirely of imperfect 69-bp terminal inverted repeats characteristic of the Foldback class of TEs. They belong to three clades of a family of maize TEs comprising hundreds of sequences in the genome of the B73 maize line. COIN elements preferentially insert at TNA sequences with a preference for C and G nucleotides in the immediately flanking 5' and 3' regions, respectively. They produce a three-base target site duplication and do not have homology to other characterized TEs. We propose that Cnr is an unstable gene that is mutated insertionally at high frequency, most commonly due to COIN element insertions at a specific site in the gene.

A Transposon Story: From TE Content to TE Dynamic Invasion of Drosophila Genomes Using the Single-Molecule Sequencing Technology from Oxford Nanopore.

Transposable elements (TEs) are the main components of genomes. However, due to their repetitive nature, they are very difficult to study using data obtained with short-read sequencing technologies. Here, we describe an efficient pipeline to accurately recover TE insertion (TEI) sites and sequences from long reads obtained by Oxford Nanopore Technology (ONT) sequencing. With this pipeline, we could precisely describe the landscapes of the most recent TEIs in wild-type strains of Drosophila melanogaster and Drosophila simulans. Their comparison suggests that this subset of TE sequences is more similar than previously thought in these two species. The chromosome assemblies obtained using this pipeline also allowed recovering piRNA cluster sequences, which was impossible using short-read sequencing. Finally, we used our pipeline to analyze ONT sequencing data from a D. melanogaster unstable line in which LTR transposition was derepressed for 73 successive generations. We could rely on single reads to identify new insertions with intact target site duplications. Moreover, the detailed analysis of TEIs in the wild-type strains and the unstable line did not support the trap model claiming that piRNA clusters are hotspots of TE insertions.

Evolution of pogo, a separate superfamily of IS630-Tc1-mariner transposons, revealing recurrent domestication events in vertebrates

sBackgroundTc1/mariner and Zator, as two superfamilies of IS630-Tc1-mariner (ITm) group, have been well-defined. However, the molecular evolution and domestication of pogo transposons, once designated as an important family of the Tc1/mariner superfamily, are still poorly understood.ResultsHere, phylogenetic analysis show that pogo transposases, together with Tc1/mariner, DD34E/Gambol, and Zator transposases form four distinct monophyletic clades with high bootstrap supports (> = 74%), suggesting that they are separate superfamilies of ITm group. The pogo superfamily represents high diversity with six distinct families (Passer, Tigger, pogoR, Lemi, Mover, and Fot/Fot-like) and wide distribution with an expansion spanning across all the kingdoms of eukaryotes. It shows widespread occurrences in animals and fungi, but restricted taxonomic distribution in land plants. It has invaded almost all lineages of animals—even mammals—and has been domesticated repeatedly in vertebrates, with 12 genes, including centromere-associated protein B (CENPB), CENPB DNA-binding domain containing 1 (CENPBD1), Jrk helix–turn–helix protein (JRK), JRK like (JRKL), pogo transposable element derived with KRAB domain (POGK), and with ZNF domain (POGZ), and Tigger transposable element-derived 2 to 7 (TIGD2–7), deduced as originating from this superfamily. Two of them (JRKL and TIGD2) seem to have been co-domesticated, and the others represent independent domestication events. Four genes (TIGD3, TIGD4, TIGD5, and POGZ) tend to represent ancient domestications in vertebrates, while the others only emerge in mammals and seem to be domesticated recently. Significant structural variations including target site duplication (TSD) types and the DDE triad signatures (DD29–56D) were observed for pogo transposons. Most domesticated genes are derived from the complete transposase genes; but CENPB, POGK, and POGZ are chimeric genes fused with additional functional domains.ConclusionsThis is the first report to systematically reveal the evolutionary profiles of the pogo transposons, suggesting that pogo and Tc1/Mariner are two separate superfamilies of ITm group, and demonstrating the repeated domestications of pogo in vertebrates. These data indicate that pogo transposons have played important roles in shaping the genome and gene evolution of fungi and animals. This study expands our understanding of the diversity of pogo transposons and updates the classification of ITm group.

Expansion of MIR482/2118 by a class-II transposable element in cotton.

Some plant microRNA (miRNA) families contain multiple members generating identical or highly similar mature miRNA variants. Mechanisms underlying the expansion of miRNA families remain elusive, although tandem and/or segmental duplications have been proposed. In this study of two tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, and their extant diploid progenitors, Gossypium arboreum and Gossypium raimondii, we investigated the gain and loss of members of the miR482/2118 superfamily, which modulates the expression of nucleotide-binding site leucine-rich repeat (NBS-LRR) disease resistance genes. We found significant expansion of MIR482/2118d in G.barbadense, G.hirsutum and G.raimondii, but not in G.arboreum. Several newly expanded MIR482/2118d loci have mutated to produce different miR482/2118 variants with altered target-gene specificity. Based on detailed analysis of sequences flanking these MIR482/2118 loci, we found that this expansion of MIR482/2118d and its derivatives resulted from an initial capture of an MIR482/2118d by a class-II DNA transposable element (TE) in G.raimondii prior to the tetraploidization event, followed by transposition to new genomic locations in G.barbadense, G.hirsutum and G.raimondii. The 'GosTE' involved in the capture and proliferation of MIR482/2118d and its derivatives belongs to the PIF/Harbinger superfamily, generating a 3-bp target site duplication upon insertion at new locations. All orthologous MIR482/2118 loci in the two diploids were retained in the two tetraploids, but mutation(s) in miR482/2118 were observed across all four species as well as in different cultivars of both G.barbadense and G.hirsutum, suggesting a dynamic co-evolution of miR482/2118 and its NBS-LRR targets. Our results provide fresh insights into the mechanisms contributing to MIRNA proliferation and enrich our knowledge on TEs.

Identification of RAG-like transposons in protostomes suggests their ancient bilaterian origin

BackgroundV(D) J recombination is essential for adaptive immunity in jawed vertebrates and is initiated by the RAG1-RAG2 endonuclease. The RAG1 and RAG2 genes are thought to have evolved from a RAGL (RAG-like) transposon containing convergently-oriented RAG1-like (RAG1L) and RAG2-like (RAG2L) genes. Elements resembling this presumptive evolutionary precursor have thus far only been detected convincingly in deuterostomes, leading to the model that the RAGL transposon first appeared in an early deuterostome.ResultsWe have identified numerous RAGL transposons in the genomes of protostomes, including oysters and mussels (phylum Mollusca) and a ribbon worm (phylum Nemertea), and in the genomes of several cnidarians. Phylogenetic analyses are consistent with vertical evolution of RAGL transposons within the Bilateria clade and with its presence in the bilaterian ancestor. Many of the RAGL transposons identified in protostomes are intact elements containing convergently oriented RAG1L and RAG2L genes flanked by terminal inverted repeats (TIRs) and target site duplications with striking similarities with the corresponding elements in deuterostomes. In addition, protostome genomes contain numerous intact RAG1L-RAG2L adjacent gene pairs that lack detectable flanking TIRs. Domains and critical active site and structural amino acids needed for endonuclease and transposase activity are present and conserved in many of the predicted RAG1L and RAG2L proteins encoded in protostome genomes.ConclusionsActive RAGL transposons were present in multiple protostome lineages and many were likely transmitted vertically during protostome evolution. It appears that RAGL transposons were broadly active during bilaterian evolution, undergoing multiple duplication and loss/fossilization events, with the RAGL genes that persist in present day protostomes perhaps constituting both active RAGL transposons and domesticated RAGL genes. Our findings raise the possibility that the RAGL transposon arose earlier in evolution than previously thought, either in an early bilaterian or prior to the divergence of bilaterians and non-bilaterians, and alter our understanding of the evolutionary history of this important group of transposons.