Search

Background: SINE-VNTR-Alu (SVA) retrotransposons move from one genomic location to another in a 'copy-and-paste' manner. They continue to move actively and cause monogenic diseases through various mechanisms. Currently, disease-causing SVA retrotransposons are classified into human-specific young SVA_E or SVA_F subfamilies. In this study, we identified an evolutionarily old SVA_D retrotransposon as a novel cause of occipital horn syndrome (OHS). OHS is an X-linked, copper metabolism disorder caused by dysfunction of the copper transporter, ATP7A., Methods: We investigated a 16-year-old boy with OHS whose pathogenic variant could not be detected via routine molecular genetic analyses., Results: A 2.8 kb insertion was detected deep within the intron of the patient's ATP7A gene. This insertion caused aberrant mRNA splicing activated by a new donor splice site located within it. Long-read circular consensus sequencing enabled us to accurately read the entire insertion sequence, which contained highly repetitive and GC-rich segments. Consequently, the insertion was identified as an SVA_D retrotransposon. Antisense oligonucleotides (AOs) targeting the new splice site restored the expression of normal transcripts and functional ATP7A proteins. AO treatment alleviated excessive accumulation of copper in patient fibroblasts in a dose-dependent manner. Pedigree analysis revealed that the retrotransposon had moved into the OHS-causing position two generations ago., Conclusion: This is the first report of a human monogenic disease caused by the SVA_D retrotransposon. The fact that the evolutionarily old SVA_D is still actively transposed, leading to increased copy numbers may make a notable impact on rare genetic disease research., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2024. No commercial re-use. See rights and permissions. Published by BMJ.)

Terminal repeat retrotransposons in miniature (TRIMs) are short non-autonomous long terminal repeat (LTR) retrotransposons found from various eukaryotes. Cassandra is a unique TRIM lineage which contains a 5S rRNA-derived sequence in its LTRs. Here, two new groups of TRIMs, designated Helenus and Ajax, are reported based on bioinformatics analysis and the usage of Repbase. Helenus is found from fungi, animals, and plants, and its LTRs contain a tRNA-like sequence. It includes two LTRs and between them, a primer-binding site (PBS) and polypurine tract (PPT) exist. Fungal and plant Helenus generate 5 bp target site duplications (TSDs) upon integration, while animal Helenus generates 4 bp TSDs. Ajax includes a 5S rRNA-derived sequence in its LTR and is found from two nemertean genomes. Ajax generates 5 bp TSDs upon integration. These results suggest that despite their unique promoters, Helenus and Ajax are TRIMs whose transposition is dependent on autonomous LTR retrotransposon. These TRIMs can originate through an insertion of SINE in an LTR of TRIM. The discovery of Helenus and Ajax suggests the presence of TRIMs with a promoter for RNA polymerase III derived from a small RNA gene, which is here collectively termed TRIMp3.

Long terminal repeat (LTR) retrotransposons are the major contributor to genome size expansion, as in the cases of the maize genome or the axolotl genome. Despite their impact on the genome size, the length of each retrotransposon is limited, compared to DNA transposons, which sometimes exceed over 100 kb. The longest LTR retrotransposon known to date is Burro-1 from the planarian Schmidtea medierranea, which is around 35.7 kb long. Here through bioinformatics analysis, a new lineage of gigantic LTR retrotransposons, designated Daidara, is reported from the springtail Allacma fusca genome. Their entire length (25-33 kb) rivals Burro families, while their LTRs are shorter than 1.5 kb, in contrast to other gigantic LTR retrotransposon lineages Burro and Ogre, whose LTRs are around 5 kb long. Daidara encodes three core proteins corresponding to gag, pol, and an additional protein of unknown function. The phylogenetic analysis supports the independent gigantification of Daidara from Burro or Ogre., (© 2023 Molecular Biology Society of Japan and John Wiley & Sons Australia, Ltd.)

The restriction enzymes examined so far are phosphodiesterases, which cleave DNA strands by hydrolysing phosphodiester bonds. Based on the mobility of restriction-modification systems, recent studies have identified a family of restriction enzymes that excise a base in their recognition sequence to generate an abasic (AP) site unless the base is properly methylated. These restriction glycosylases also show intrinsic but uncoupled AP lyase activity at the AP site, generating an atypical strand break. Action of an AP endonuclease at the AP site may generate another atypical break, rejoining/repairing of which is difficult. This PabI family of restriction enzymes contain a novel fold (HALFPIPE) and show unusual properties, such as non-requirement of divalent cations for cleavage. These enzymes are present in Helicobacteraceae/Campylobacteraceae and in few hyperthermophilic archaeal species. In Helicobacter genomes, their recognition sites are strongly avoided, and the encoding genes are often inactivated by mutations or replacement, indicating that their expression is toxic for the cells. The discovery of restriction glycosylases generalizes the concept of restriction-modification systems to epigenetic immune systems, which may use any mode of damage to DNA that are considered 'non-self' based on epigenetic modifications. This concept will add to our understanding of immunity and epigenetics., (© The Author(s) 2023. Published by Oxford University Press on behalf of Kazusa DNA Research Institute.)

DDD/E transposase gene is the most abundant gene in nature and many DNA transposons in all three domains of life use it for their transposition. A substantial number of eukaryotic DNA transposons show similarity to prokaryotic insertion sequences (ISs). The presence of IS 481 -like DNA transposons was indicated in the genome of Trichomonas vaginalis . Here, we surveyed IS 481 -like eukaryotic sequences using a bioinformatics approach and report a group of eukaryotic IS 481 -like DNA transposons, designated IS481EU , from parabasalids including T. vaginalis . The lengths of target site duplications (TSDs) of IS481EU are around 4 bps, around 15 bps, or around 25 bps, and strikingly, these discrete lengths of TSDs can be observed even in a single IS481EU family. Phylogenetic analysis indicated the close relationships of IS481EU with some of the prokaryotic IS 481 family members. IS481EU was not well separated from IS3EU / GingerRoot in the phylogenetic analysis, but was distinct from other eukaryotic DNA transposons including Ginger1 and Ginger2 . The unique characteristics of IS481EU in protein sequences and the distribution of TSD lengths support its placement as a new superfamily of eukaryotic DNA transposons.

Background: DNA transposons are ubiquitous components of eukaryotic genomes. A major group of them encode a DDD/E transposase and contain terminal inverted repeats (TIRs) of varying lengths. The Kolobok superfamily of DNA transposons has been found in a wide spectrum of organisms., Results: Here we report a new Kolobok lineage, designated KolobokP. They were identified in 7 animal phyla (Mollusca, Phoronida, Annelida, Nemertea, Bryozoa, Chordata, and Echinodermata), and are especially rich in bivalves. Unlike other Kolobok families, KolobokP adopts a composite-like architecture: an internal region (INT) flanked by two long terminal direct repeats (LTDRs), which exhibit their own short terminal inverted repeats ranging up to 18 bps. The excision of LTDRs was strongly suggested. The LTDR lengths seem to be constrained to be either around 450-bp or around 660-bp. The internal region encodes a DDD/E transposase and a small His-Me finger nuclease, which likely originated from the homing endonuclease encoded by a group I intron from a eukaryotic species. The architecture of KolobokP resembles composite DNA transposons, usually observed in bacterial genomes, and long terminal repeat (LTR) retrotransposons. In addition to this monomeric LTDR-INT-LTDR structure, plenty of solo LTDRs and multimers represented as (LTDR-INT) n -LTDR are also observed. Our structural and phylogenetic analysis supported the birth of KolobokP in the late stage of the Kolobok evolution. We propose KolobokP families propagate themselves in two ways: the canonical transposition catalyzed by their transposase and the sequence-specific cleavage by their endonuclease followed by the multimerization through the unequal crossover., Conclusions: The presence of homing endonuclease and long terminal direct repeats of KolobokP families suggest their unique dual replication mechanisms: transposition and induced unequal crossover., (© 2022. The Author(s).)

Objectives: The aim of the present study was to clarify the current state of outpatient cancer rehabilitation and coordination systems provided by designated cancer hospitals in Japan., Methods: A questionnaire was sent to 427 designated cancer hospitals in Japan to investigate the status of outpatient cancer rehabilitation and whether it was sufficiently conducted. The status of regional coordination with post-discharge rehabilitation facilities was surveyed., Results: Responses were received from 235/427 facilities (55.0%). Outpatient cancer rehabilitation was implemented in 92 (39.1% of responding facilities), and of these facilities, 83.7% answered that the provision of rehabilitation was insufficient. The reasons were ineligibility for reimbursement of medical fees, a lack of human resources, a lack of awareness of the need, and a lack of education. Regional coordination was conducted by 39.1% of responding facilities, yet a regional alliance path had been established in only 9.8% of centers. The absence of coordination was associated with large facility size, the absence of physiatrists, and few rehabilitation professionals who had completed the training program; an insufficient framework for regional coordination was also given as a reason., Conclusions: To provide adequate outpatient cancer rehabilitation, sufficient human resources, the reimbursement of medical fees in the outpatient setting, and education and a framework to promote regional coordination are necessary., Competing Interests: CONFLICTS OF INTEREST: The authors declare that they have no conflicts of interest., (2022 The Japanese Association of Rehabilitation Medicine.)

Dada is a unique superfamily of DNA transposons, inserted specifically in multicopy RNA genes. The zebrafish genome harbors five families of Dada transposons, whose targets are U6 and U1 snRNA genes, and tRNA-Ala and tRNA-Leu genes. Dada-U6 , which is inserted specifically in U6 snRNA genes, is found in four animal phyla, but other target-specific lineages have been reported only from one or two species. Here, vertebrate genomes and transcriptomes were surveyed to characterize Dada families with new target specificities, and over 120 Dada families were characterized from the genomes of actinopterygian fish. They were classified into 12 groups with confirmed target specificities. Newly characterized Dada families target tRNA genes for Asp, Asn, Arg, Gly, Lys, Ser, Tyr, and Val, and 5S rRNA genes. Targeted positions inside of tRNA genes are concentrated in two regions: around the anticodon and the A box of RNA polymerase III promoter. Phylogenetic analysis revealed the relationships among actinopterygian Dada families, and one domestication event in the common ancestor of carps and minnows belonging to Cyprinoidei, Cypriniformes. Sequences targeted by phylogenetically related Dada families show sequence similarities, indicating that the target specificity of Dada is accomplished through the recognition of primary nucleotide sequences.

Background: Short interspersed elements (SINEs) are ubiquitous components of eukaryotic genomes. SINEs are composite transposable elements that are mobilized by non-long terminal repeat (non-LTR) retrotransposons, also called long interspersed elements (LINEs). The 3' part of SINEs usually originated from that of counterpart non-LTR retrotransposons. The 5' part of SINEs mostly originated from small RNA genes. SINE1 is a group of SINEs whose 5' part originated from 7SL RNA, and is represented by primate Alu and murine B1 . Well-defined SINE1 has been found only from Euarchontoglires, a group of mammals, in contrast to the wide distribution of SINE2, which has a tRNA-derived sequence, from animals to plants to protists. Both Alu and B1 are mobilized by L1 -type non-LTR retrotransposons, which are the only lineage of autonomous non-LTR retrotransposons active in these mammalian lineages., Results: Here a new lineage of SINE1 is characterized from the seashore hagfish Eptatretus burgeri genome. This SINE1 family, designated SINE1 - 1_EBu , is young, and is transposed by RTE -type non-LTR retrotransposon, not L1 -type. Comparison with other SINE families from hagfish indicated the birth of SINE1 - 1_EBu through chimera formation of a 7SL RNA-derived sequence and an older tRNA-derived SINE family. It reveals parallel evolution of SINE1 in two vertebrate lineages with different autonomous non-LTR retrotransposon partners. The comparison between two SINE1 lineages supports that the RNA secondary structure of the Alu domain of 7SL RNA is required for the efficient retrotransposition., Conclusions: The hagfish SINE1 is the first evident SINE1 family found outside of Euarchontoglires. Independent evolution of SINE1 with similar RNA secondary structure originated in 7SL RNA indicates the functional importance of 7SL RNA-derived sequence in the proliferation of SINEs., Competing Interests: Competing interestsThe authors declare that they have no competing interests., (© The Author(s) 2020.)

Bioinformatic and molecular characterization of satellite repeats was performed to understand the impact of their diversification on Vaccinium genome evolution. Satellite repeat diversity was evaluated in four cultivated and wild species, including the diploid species Vaccinium myrtillus and Vaccinium uliginosum , as well as the tetraploid species Vaccinium corymbosum and Vaccinium arctostaphylos . We comparatively characterized six satellite repeat families using in total 76 clones with 180 monomers. We observed that the monomer units of VaccSat1, VaccSat2, VaccSat5, and VaccSat6 showed a higher order repeat (HOR) structure, likely originating from the organization of two adjacent subunits with differing similarity, length and size. Moreover, VaccSat1, VaccSat3, VaccSat6, and VaccSat7 were found to have sequence similarity to parts of transposable elements. We detected satellite-typical tandem organization for VaccSat1 and VaccSat2 in long arrays, while VaccSat5 and VaccSat6 distributed in multiple sites over all chromosomes of tetraploid V. corymbosum , presumably in long arrays. In contrast, very short arrays of VaccSat3 and VaccSat7 are dispersedly distributed over all chromosomes in the same species, likely as internal parts of transposable elements. We provide a comprehensive overview on satellite species specificity in Vaccinium , which are potentially useful as molecular markers to address the taxonomic complexity of the genus, and provide information for genome studies of this genus.

Background: DNA transposons are ubiquitous components of eukaryotic genomes. Academ superfamily of DNA transposons is one of the least characterized DNA transposon superfamilies in eukaryotes. DNA transposons belonging to the Academ superfamily have been reported from various animals, one red algal species Chondrus crispus , and one fungal species Puccinia graminis . Six Academ families from P. graminis encode a helicase in addition to putative transposase, while some other families encode a single protein which contains a putative transposase and an XPG nuclease., Results: Systematic searches on Repbase and BLAST searches against publicly available genome sequences revealed that several species of fungi and animals contain multiple Academ transposon families encoding a helicase. These AcademH families generate 9 or 10-bp target site duplications (TSDs) while Academ families lacking helicase generate 3 or 4-bp TSDs. Phylogenetic analysis clearly shows two lineages inside of Academ , designated here as AcademH and AcademX for encoding helicase or XPG nuclease, respectively. One sublineage of AcademH in animals encodes plant homeodomain (PHD) finger in its transposase, and its remnants are found in several fish genomes., Conclusions: The AcademH lineage of TEs is widely distributed in animals and fungi, and originated early in the evolution of Academ DNA transposons. This analysis highlights the structural diversity in one less studied superfamily of eukaryotic DNA transposons., Competing Interests: Competing interestsThe author declare that he has no competing interests., (© The Author(s) 2020.)

The majority of eukaryotic genomes contain a large fraction of repetitive sequences that primarily originate from transpositional bursts of transposable elements (TEs). Repbase serves as a database for eukaryotic repetitive sequences and has now become the largest collection of eukaryotic TEs. During the development of Repbase, many new superfamilies/lineages of TEs, which include Helitron, Polinton, Ginger and SINEU, were reported. The unique composition of protein domains and DNA motifs in TEs sometimes indicates novel mechanisms of transposition, replication, anti-suppression or proliferation. In this review, our current understanding regarding the diversity of eukaryotic TEs in sequence, protein domain composition and structural hallmarks is introduced and summarized, based on the classification system implemented in Repbase. Autonomous eukaryotic TEs can be divided into two groups: Class I TEs, also called retrotransposons, and Class II TEs, or DNA transposons. Long terminal repeat (LTR) retrotransposons, including endogenous retroviruses, non-LTR retrotransposons, tyrosine recombinase retrotransposons and Penelope-like elements, are well accepted groups of autonomous retrotransposons. They share reverse transcriptase for replication but are distinct in the catalytic components responsible for integration into the host genome. Similarly, at least three transposition machineries have been reported in eukaryotic DNA transposons: DDD/E transposase, tyrosine recombinase and HUH endonuclease combined with helicase. Among these, TEs with DDD/E transposase are dominant and are classified into 21 superfamilies in Repbase. Non-autonomous TEs are either simple derivatives generated by internal deletion, or are composed of several units that originated independently.

Repbase is a comprehensive database of eukaryotic transposable elements (TEs) and repeat sequences, containing over 1300 human repeat sequences. Recent analyses of these repeat sequences have accumulated evidences for their contribution to human evolution through becoming functional elements, such as protein-coding regions or binding sites of transcriptional regulators. However, resolving the origins of repeat sequences is a challenge, due to their age, divergence, and degradation. Ancient repeats have been continuously classified as TEs by finding similar TEs from other organisms. Here, the most comprehensive picture of human repeat sequences is presented. The human genome contains traces of 10 clades ( L1, CR1, L2, Crack, RTE, RTEX, R4, Vingi, Tx1 and Penelope ) of non-long terminal repeat (non-LTR) retrotransposons (long interspersed elements, LINEs), 3 types ( SINE1/7SL, SINE2/tRNA , and SINE3/5S ) of short interspersed elements (SINEs), 1 composite retrotransposon ( SVA ) family, 5 classes ( ERV1, ERV2, ERV3 , Gypsy and DIRS ) of LTR retrotransposons, and 12 superfamilies ( Crypton , Ginger1 , Harbinger , hAT, Helitron, Kolobok, Mariner, Merlin, MuDR, P, piggyBac and Transib ) of DNA transposons. These TE footprints demonstrate an evolutionary continuum of the human genome., Competing Interests: Not applicable.Not applicable.The author declares that he has no competing interests.Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Short interspersed elements (SINEs), which are nonautonomous transposable elements, require the transposition machinery of long interspersed elements (LINEs) to mobilize. SINEs are composed of two or more independently originating parts. The 5' region is called the "head" and is derived mainly from small RNAs, and the 3' region ("tail") originates from the 3' region of LINEs and is responsible for being recognized by counterpart LINE proteins. The origin of the middle "body" of SINEs is enigmatic, although significant sequence similarities among SINEs from very diverse species have been observed. Here, a systematic analysis of the similarities among SINEs and LINEs deposited on Repbase, a comprehensive database of eukaryotic repeat sequences was performed. Three primary findings are described: 1) The 5' regions of only two clades of LINEs, RTE and Vingi, were revealed to have contributed to the middle parts of SINEs; 2) The linkage of the 5' and 3' parts of LINEs can be lost due to occasional tail exchange of SINEs; and 3) The previously proposed Ceph-domain was revealed to be a fusion of a CORE-domain and a 5' part of RTE clade of LINE. Based on these findings, a hypothesis that the 5' parts of bipartite nonautonomous LINEs, which possess only the 5' and 3' regions of the original LINEs, can contribute to the undefined middle part of SINEs is proposed., (© The Author(s) 2018. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

A bacterial insertion sequence (IS) is a mobile DNA sequence carrying only the transposase gene (tnp) that acts as a mutator to disrupt genes, alter gene expressions, and cause genomic rearrangements. "Canonical" ISs have historically been characterized by their terminal inverted repeats (IRs), which may form a stem-loop structure, and duplications of a short (non-IR) target sequence at both ends, called target site duplications (TSDs). The IS distributions and virulence potentials of Staphylococcus aureus genomes in familial infection cases are unclear. Here, we determined the complete circular genome sequences of familial strains from a Panton-Valentine leukocidin (PVL)-positive ST50/agr4 S. aureus (GN) infection of a 4-year old boy with skin abscesses. The genomes of the patient strain (GN1) and parent strain (GN3) were rich for "canonical" IS1272 with terminal IRs, both having 13 commonly-existing copies (ce-IS1272). Moreover, GN1 had a newly-inserted IS1272 (ni-IS1272) on the PVL-converting prophage, while GN3 had two copies of ni-IS1272 within the DNA helicase gene and near rot. The GN3 genome also had a small deletion. The targets of ni-IS1272 transposition were IR structures, in contrast with previous "canonical" ISs. There were no TSDs. Based on a database search, the targets for ce-IS1272 were IRs or "non-IRs". IS1272 included a larger structure with tandem duplications of the left (IRL) side sequence; tnp included minor cases of a long fusion form and truncated form. One ce-IS1272 was associated with the segments responsible for immune evasion and drug resistance. Regarding virulence, GN1 expressed cytolytic peptides (phenol-soluble modulin α and δ-hemolysin) and PVL more strongly than some other familial strains. These results suggest that IS1272 transposes through an IR-replacing mechanism, with an irreversible process unlike that of "canonical" transpositions, resulting in genomic variations, and that, among the familial strains, the patient strain has strong virulence potential based on community-associated virulence factors.

Decoding of closely related genomes is now revealing the process of population evolution. In bacteria, population divergence appears associated with a unique set of sequence-specific epigenetic DNA methylation systems, often within restriction-modification (RM) systems. They might define a unique gene expression pattern and limit genetic flux between lineages in population divergence. We addressed the contribution of methylation systems to population diversification in panmictic bacterial species, Helicobacter pylori, which shows an interconnected population structure through frequent mutual recombination. We analyzed complete genome sequences of 28 strains collected in Fukui, Japan. Their nucleotide sequences are closely related although fine-scale analyses revealed two subgroups likely reflecting human subpopulations. Their sequences are tightly connected by homologous recombination. Our extensive analysis of RM systems revealed an extreme variability in DNA methyltransferases, especially in their target recognition domains. Their diversity was, however, not immediately related to the genome sequence diversity, except for very closely related strains. An interesting exception is a hybrid strain, which likely has conserved the methylation gene repertoire from one parent but diversified in sequence by massive acquisition of fragmentary DNA sequences from the other parent. Our results demonstrate how a bacterial population can be extremely divergent in epigenetics and yet homogenized in sequence., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Transposons, or transposable elements, are the major components of genomes in most eukaryotes. Some groups of transposons have developed target specificity that limits the integration sites to a specific nonessential sequence or a genomic region to avoid gene disruption caused by insertion into an essential gene. R2 is one of the most intensively investigated groups of sequence-specific non-LTR retrotransposons and is inserted at a specific site inside of 28S ribosomal RNA (rRNA) genes. R2 is known to be distributed among at least six animal phyla even though its occurrence is reported to be patchy. Here, in order to obtain a more detailed picture of the distribution of R2, we surveyed R2 using both in silico screening and degenerate PCR, particularly focusing on actinopterygian fish. We found two families of the R2C lineage from vertebrates, although it has previously only been found in platyhelminthes. We also revealed the apparent movement of insertion sites of a lineage of actinopterygian R2, which was likely concurrent with the acquisition of a 28S rRNA-derived sequence in their 3' UTR. Outside of actinopterygian fish, we revealed the maintenance of a single R2 lineage in birds; the co-existence of four lineages of R2 in the leafcutter bee Megachile rotundata; the first examples of R2 in Ctenophora, Mollusca, and Hemichordata; and two families of R2 showing no target specificity. These findings indicate that R2 is relatively stable and universal, while differences in the distribution and maintenance of R2 lineages probably reflect characteristics of some combination of both R2 lineages and host organisms., Competing Interests: The authors have declared that no competing interests exist.

Genomes of a given bacterial species can show great variation in gene content and thus systematic analysis of the entire gene repertoire, termed the pan-genome, is important for understanding bacterial intra-species diversity, population genetics, and evolution. Here, we analyzed the pan-genome from 30 completely sequenced strains of the human gastric pathogen Helicobacter pylori belonging to various phylogeographic groups, focusing on 991 accessory (not fully conserved) orthologous groups (OGs). We developed a method to evaluate the mobility of genes within a genome, using the gene order in the syntenically conserved regions as a reference, and classified the 991 accessory OGs into five classes: Core, Stable, Intermediate, Mobile, and Unique. Phylogenetic networks based on the gene content of Core and Stable classes are highly congruent with that created from the concatenated alignment of fully conserved core genes, in contrast to those of Intermediate and Mobile classes, which show quite different topologies. By clustering the accessory OGs on the basis of phylogenetic pattern similarity and chromosomal proximity, we identified 60 co-occurring gene clusters (CGCs). In addition to known genomic islands, including cag pathogenicity island, bacteriophages, and integrating conjugative elements, we identified some novel ones. One island encodes TerY-phosphorylation triad, which includes the eukaryote-type protein kinase/phosphatase gene pair, and components of type VII secretion system. Another one contains a reverse-transcriptase homolog, which may be involved in the defense against phage infection through altruistic suicide. Many of the CGCs contained restriction-modification (RM) genes. Different RM systems sometimes occupied the same (orthologous) locus in the strains. We anticipate that our method will facilitate pan-genome studies in general and help identify novel genomic islands in various bacterial species.

Genes encoding reverse transcriptases (RTs) are found in most eukaryotes, often as a component of retrotransposons, as well as in retroviruses and in prokaryotic retroelements. We investigated the abundance, classification and transcriptional status of RTs based on Tara Oceans marine metagenomes and metatranscriptomes encompassing a wide organism size range. Our analyses revealed that RTs predominate large-size fraction metagenomes (>5 μm), where they reached a maximum of 13.5% of the total gene abundance. Metagenomic RTs were widely distributed across the phylogeny of known RTs, but many belonged to previously uncharacterized clades. Metatranscriptomic RTs showed distinct abundance patterns across samples compared with metagenomic RTs. The relative abundances of viral and bacterial RTs among identified RT sequences were higher in metatranscriptomes than in metagenomes and these sequences were detected in all metatranscriptome size fractions. Overall, these observations suggest an active proliferation of various RT-assisted elements, which could be involved in genome evolution or adaptive processes of plankton assemblage.

Most non-long terminal repeat (non-LTR) retrotransposons encoding a restriction-like endonuclease show target-specific integration into repetitive sequences such as ribosomal RNA genes and microsatellites. However, only a few target-specific lineages of non-LTR retrotransposons are distributed widely and no lineage is found across the eukaryotic kingdoms. Here we report the most widely distributed lineage of target sequence-specific non-LTR retrotransposons, designated Utopia. Utopia is found in three supergroups of eukaryotes: Amoebozoa, SAR, and Opisthokonta. Utopia is inserted into a specific site of U2 small nuclear RNA genes with different strength of specificity for each family. Utopia families from oomycetes and wasps show strong target specificity while only a small number of Utopia copies from reptiles are flanked with U2 snRNA genes. Oomycete Utopia families contain an "archaeal" RNase H domain upstream of reverse transcriptase (RT), which likely originated from a plant RNase H gene. Analysis of Utopia from oomycetes indicates that multiple lineages of Utopia have been maintained inside of U2 genes with few copy numbers. Phylogenetic analysis of RT suggests the monophyly of Utopia, and it likely dates back to the early evolution of eukaryotes.

Background: R.PabI is an exceptional restriction enzyme that functions as a DNA glycosylase. The enzyme excises an unmethylated base from its recognition sequence to generate apurinic/apyrimidinic (AP) sites, and also displays AP lyase activity, cleaving the DNA backbone at the AP site to generate the 3'-phospho alpha, beta-unsaturated aldehyde end in addition to the 5'-phosphate end. The resulting ends are difficult to religate with DNA ligase. The enzyme was originally isolated in Pyrococcus, a hyperthermophilic archaeon, and additional homologs subsequently identified in the epsilon class of the Gram-negative bacterial phylum Proteobacteria, such as Helicobacter pylori., Results: Systematic analysis of R.PabI homologs and their neighboring genes in sequenced genomes revealed co-occurrence of R.PabI with M.PabI homolog methyltransferase genes. R.PabI and M.PabI homolog genes are occasionally found at corresponding (orthologous) loci in different species, such as Helicobacter pylori, Helicobacter acinonychis and Helicobacter cetorum, indicating long-term maintenance of the gene pair. One R.PabI and M.PabI homolog gene pair is observed immediately after the GMP synthase gene in both Campylobacter and Helicobacter, representing orthologs beyond genera. The mobility of the PabI family of restriction-modification (RM) system between genomes is evident upon comparison of genomes of sibling strains/species. Analysis of R.PabI and M.PabI homologs in H. pylori revealed an insertion of integrative and conjugative elements (ICE), and replacement with a gene of unknown function that may specify a membrane-associated toxin (hrgC). In view of the similarity of HrgC with toxins in type I toxin-antitoxin systems, we addressed the biological significance of this substitution. Our data indicate that replacement with hrgC occurred in the common ancestor of hspAmerind and hspEAsia. Subsequently, H. pylori with and without hrgC were intermixed at this locus, leading to complex distribution of hrgC in East Asia and the Americas. In Malaysia, hrgC was horizontally transferred from hspEAsia to hpAsia2 strains., Conclusions: The PabI family of RM system behaves as a mobile, selfish genetic element, similar to the other families of Type II RM systems. Our analysis additionally revealed some cases of long-term inheritance. The distribution of the hrgC gene replacing the PabI family in the subpopulations of H. pylori, hspAmerind, hspEAsia and hpAsia2, corresponds to the two human migration events, one from East Asia to Americas and the other from China to Malaysia.

Repbase Update (RU) is a database of representative repeat sequences in eukaryotic genomes. Since its first development as a database of human repetitive sequences in 1992, RU has been serving as a well-curated reference database fundamental for almost all eukaryotic genome sequence analyses. Here, we introduce recent updates of RU, focusing on technical issues concerning the submission and updating of Repbase entries and will give short examples of using RU data. RU sincerely invites a broader submission of repeat sequences from the research community.

Eukaryotic genomes are colonized by various transposons including short interspersed elements (SINEs). The 5' region (head) of the majority of SINEs is derived from one of the three types of RNA genes--7SL RNA, transfer RNA (tRNA), or 5S ribosomal RNA (rRNA)--and the internal promoter inside the head promotes the transcription of the entire SINEs. Here I report a new group of SINEs whose heads originate from either the U1 or U2 small nuclear RNA gene. These SINEs, named SINEU, are distributed among crocodilians and classified into three families. The structures of the SINEU-1 subfamilies indicate the recurrent addition of a U1- or U2-derived sequence onto the 5' end of SINEU-1 elements. SINEU-1 and SINEU-3 are ancient and shared among alligators, crocodiles, and gharials, while SINEU-2 is absent in the alligator genome. SINEU-2 is the only SINE family that was active after the split of crocodiles and gharials. All SINEU families, especially SINEU-3, are preferentially inserted into a family of Mariner DNA transposon, Mariner-N4_AMi. A group of Tx1 non-long terminal repeat retrotransposons designated Tx1-Mar also show target preference for Mariner-N4_AMi, indicating that SINEU was mobilized by Tx1-Mar., (© The Author(s) 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background: Crocodilians are thought to be hosts to a diverse and divergent complement of endogenous retroviruses (ERVs) but a comprehensive investigation is yet to be performed. The recent sequencing of three crocodilian genomes provides an opportunity for a more detailed and accurate representation of the ERV diversity that is present in these species. Here we investigate the diversity, distribution and evolution of ERVs from the genomes of three key crocodilian species, and outline the key processes driving crocodilian ERV proliferation and evolution., Results: ERVs and ERV related sequences make up less than 2% of crocodilian genomes. We recovered and described 45 ERV groups within the three crocodilian genomes, many of which are species specific. We have also revealed a new class of ERV, ERV4, which appears to be common to crocodilians and turtles, and currently has no characterised exogenous counterpart. For the first time, we formally describe the characteristics of this ERV class and its classification relative to other recognised ERV and retroviral classes. This class shares some sequence similarity and sequence characteristics with ERV3, although it is phylogenetically distinct from the other ERV classes. We have also identified two instances of gene capture by crocodilian ERVs, one of which, the capture of a host KIT-ligand mRNA has occurred without the loss of an ERV domain., Conclusions: This study indicates that crocodilian ERVs comprise a wide variety of lineages, many of which appear to reflect ancient infections. In particular, ERV4 appears to have a limited host range, with current data suggesting that it is confined to crocodilians and some lineages of turtles. Also of interest are two ERV groups that demonstrate evidence of host gene capture. This study provides a framework to facilitate further studies into non-mammalian vertebrates and highlights the need for further studies into such species.

To provide context for the diversification of archosaurs--the group that includes crocodilians, dinosaurs, and birds--we generated draft genomes of three crocodilians: Alligator mississippiensis (the American alligator), Crocodylus porosus (the saltwater crocodile), and Gavialis gangeticus (the Indian gharial). We observed an exceptionally slow rate of genome evolution within crocodilians at all levels, including nucleotide substitutions, indels, transposable element content and movement, gene family evolution, and chromosomal synteny. When placed within the context of related taxa including birds and turtles, this suggests that the common ancestor of all of these taxa also exhibited slow genome evolution and that the comparatively rapid evolution is derived in birds. The data also provided the opportunity to analyze heterozygosity in crocodilians, which indicates a likely reduction in population size for all three taxa through the Pleistocene. Finally, these data combined with newly published bird genomes allowed us to reconstruct the partial genome of the common ancestor of archosaurs, thereby providing a tool to investigate the genetic starting material of crocodilians, birds, and dinosaurs., (Copyright © 2014, American Association for the Advancement of Science.)

Chicken repeat 1 (CR1) retroposons are long interspersed elements (LINEs) that are ubiquitous within amniote genomes and constitute the most abundant family of transposed elements in birds, crocodilians, turtles, and snakes. They are also present in mammalian genomes, where they reside as numerous relics of ancient retroposition events. Yet, despite their relevance for understanding amniote genome evolution, the diversity and evolution of CR1 elements has never been studied on an amniote-wide level. We reconstruct the temporal and quantitative activity of CR1 subfamilies via presence/absence analyses across crocodilian phylogeny and comparative analyses of 12 crocodilian genomes, revealing relative genomic stasis of retroposition during genome evolution of extant Crocodylia. Our large-scale phylogenetic analysis of amniote CR1 subfamilies suggests the presence of at least seven ancient CR1 lineages in the amniote ancestor; and amniote-wide analyses of CR1 successions and quantities reveal differential retention (presence of ancient relics or recent activity) of these CR1 lineages across amniote genome evolution. Interestingly, birds and lepidosaurs retained the fewest ancient CR1 lineages among amniotes and also exhibit smaller genome sizes. Our study is the first to analyze CR1 evolution in a genome-wide and amniote-wide context and the data strongly suggest that the ancestral amniote genome contained myriad CR1 elements from multiple ancient lineages, and remnants of these are still detectable in the relatively stable genomes of crocodilians and turtles. Early mammalian genome evolution was thus characterized by a drastic shift from CR1 prevalence to dominance and hyperactivity of L2 LINEs in monotremes and L1 LINEs in therians., (© The Author(s) 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Target-specific integration of transposable elements for multicopy genes, such as ribosomal RNA and small nuclear RNA (snRNA) genes, is of great interest because of the relatively harmless nature, stable inheritance and possible application for targeted gene delivery of target-specific transposable elements. To date, such strict target specificity has been observed only among non-LTR retrotransposons. We here report a new superfamily of sequence-specific DNA transposons, designated Dada. Dada encodes a DDE-type transposase that shows a distant similarity to transposases encoded by eukaryotic MuDR, hAT, P and Kolobok transposons, as well as the prokaryotic IS256 insertion element. Dada generates 6-7 bp target site duplications upon insertion. One family of Dada DNA transposons targets a specific site inside the U6 snRNA genes and are found in various fish species, water flea, oyster and polycheate worm. Other target sequences of the Dada transposons are U1 snRNA genes and different tRNA genes. The targets are well conserved in multicopy genes, indicating that copy number and sequence conservation are the primary constraints on the target choice of Dada transposons. Dada also opens a new frontier for target-specific gene delivery application.

Background: Mammalian genomes are repositories of repetitive DNA sequences derived from transposable elements (TEs). Typically, TEs generate multiple, mostly inactive copies of themselves, commonly known as repetitive families or families of repeats. Recently, we proposed that families of TEs originate in small populations by genetic drift and that the origin of small subpopulations from larger populations can be fueled by biological innovations., Results: We report three distinct groups of repetitive families preserved in the human genome that expanded and declined during the three previously described periods of regulatory innovations in vertebrate genomes. The first group originated prior to the evolutionary separation of the mammalian and bird lineages and the second one during subsequent diversification of the mammalian lineages prior to the origin of eutherian lineages. The third group of families is primate-specific., Conclusions: The observed correlation implies a relationship between regulatory innovations and the origin of repetitive families. Consistent with our previous hypothesis, it is proposed that regulatory innovations fueled the origin of new subpopulations in which new repetitive families became fixed by genetic drift.

Background: "Domestication" of transposable elements (TEs) led to evolutionary breakthroughs such as the origin of telomerase and the vertebrate adaptive immune system. These breakthroughs were accomplished by the adaptation of molecular functions essential for TEs, such as reverse transcription, DNA cutting and ligation or DNA binding. Cryptons represent a unique class of DNA transposons using tyrosine recombinase (YR) to cut and rejoin the recombining DNA molecules. Cryptons were originally identified in fungi and later in the sea anemone, sea urchin and insects., Results: Herein we report new Cryptons from animals, fungi, oomycetes and diatom, as well as widely conserved genes derived from ancient Crypton domestication events. Phylogenetic analysis based on the YR sequences supports four deep divisions of Crypton elements. We found that the domain of unknown function 3504 (DUF3504) in eukaryotes is derived from Crypton YR. DUF3504 is similar to YR but lacks most of the residues of the catalytic tetrad (R-H-R-Y). Genes containing the DUF3504 domain are potassium channel tetramerization domain containing 1 (KCTD1), KIAA1958, zinc finger MYM type 2 (ZMYM2), ZMYM3, ZMYM4, glutamine-rich protein 1 (QRICH1) and "without children" (WOC). The DUF3504 genes are highly conserved and are found in almost all jawed vertebrates. The sequence, domain structure, intron positions and synteny blocks support the view that ZMYM2, ZMYM3, ZMYM4, and possibly QRICH1, were derived from WOC through two rounds of genome duplication in early vertebrate evolution. WOC is observed widely among bilaterians. There could be four independent events of Crypton domestication, and one of them, generating WOC/ZMYM, predated the birth of bilaterian animals. This is the third-oldest domestication event known to date, following the domestication generating telomerase reverse transcriptase (TERT) and Prp8. Many Crypton-derived genes are transcriptional regulators with additional DNA-binding domains, and the acquisition of the DUF3504 domain could have added new regulatory pathways via protein-DNA or protein-protein interactions., Conclusions: Cryptons have contributed to animal evolution through domestication of their YR sequences. The DUF3504 domains are domesticated YRs of animal Crypton elements.

Background: Eukaryotic genomes harbor diverse families of repetitive DNA derived from transposable elements (TEs) that are able to replicate and insert into genomic DNA. The biological role of TEs remains unclear, although they have profound mutagenic impact on eukaryotic genomes and the origin of repetitive families often correlates with speciation events. We present a new hypothesis to explain the observed correlations based on classical concepts of population genetics., Presentation of the Hypothesis: The main thesis presented in this paper is that the TE-derived repetitive families originate primarily by genetic drift in small populations derived mostly by subdivisions of large populations into subpopulations. We outline the potential impact of the emerging repetitive families on genetic diversification of different subpopulations, and discuss implications of such diversification for the origin of new species., Testing the Hypothesis: Several testable predictions of the hypothesis are examined. First, we focus on the prediction that the number of diverse families of TEs fixed in a representative genome of a particular species positively correlates with the cumulative number of subpopulations (demes) in the historical metapopulation from which the species has emerged. Furthermore, we present evidence indicating that human AluYa5 and AluYb8 families might have originated in separate proto-human subpopulations. We also revisit prior evidence linking the origin of repetitive families to mammalian phylogeny and present additional evidence linking repetitive families to speciation based on mammalian taxonomy. Finally, we discuss evidence that mammalian orders represented by the largest numbers of species may be subject to relatively recent population subdivisions and speciation events., Implications of the Hypothesis: The hypothesis implies that subdivision of a population into small subpopulations is the major step in the origin of new families of TEs as well as of new species. The origin of new subpopulations is likely to be driven by the availability of new biological niches, consistent with the hypothesis of punctuated equilibria. The hypothesis also has implications for the ongoing debate on the role of genetic drift in genome evolution.

Background: The wheat stripe rust fungus (Puccinia striiformis f. sp. tritici, PST) is responsible for significant yield losses in wheat production worldwide. In spite of its economic importance, the PST genomic sequence is not currently available. Fortunately Next Generation Sequencing (NGS) has radically improved sequencing speed and efficiency with a great reduction in costs compared to traditional sequencing technologies. We used Illumina sequencing to rapidly access the genomic sequence of the highly virulent PST race 130 (PST-130)., Methodology/principal Findings: We obtained nearly 80 million high quality paired-end reads (>50x coverage) that were assembled into 29,178 contigs (64.8 Mb), which provide an estimated coverage of at least 88% of the PST genes and are available through GenBank. Extensive micro-synteny with the Puccinia graminis f. sp. tritici (PGTG) genome and high sequence similarity with annotated PGTG genes support the quality of the PST-130 contigs. We characterized the transposable elements present in the PST-130 contigs and using an ab initio gene prediction program we identified and tentatively annotated 22,815 putative coding sequences. We provide examples on the use of comparative approaches to improve gene annotation for both PST and PGTG and to identify candidate effectors. Finally, the assembled contigs provided an inventory of PST repetitive elements, which were annotated and deposited in Repbase., Conclusions/significance: The assembly of the PST-130 genome and the predicted proteins provide useful resources to rapidly identify and clone PST genes and their regulatory regions. Although the automatic gene prediction has limitations, we show that a comparative genomics approach using multiple rust species can greatly improve the quality of gene annotation in these species. The PST-130 sequence will also be useful for comparative studies within PST as more races are sequenced. This study illustrates the power of NGS for rapid and efficient access to genomic sequence in non-model organisms.

Background: Long interspersed nuclear element-1 (LINE-1 or L1) is a dominant repetitive sequence in the human genome. Besides mediating its own retrotransposition, L1 can mobilize Alu and messenger RNA (mRNA) in trans, and probably also SVA and non-coding RNA. The structures of L1 copies and trans-mobilized retrocopies are variable and can be classified into three categories: full-length; 5'-truncated; and 5'-inverted insertions. These structures may be generated by different 5' integration mechanisms., Results: In this study, a method to correctly characterize insertions with short target site duplications (TSDs) is developed and extranucleotides, TSDs and microhomologies (MHs) at junctions were analysed for the three types of insertions. Only 5'-truncated L1 insertions were found to be associated with short TSDs. Both full-length and 5'-truncated retrotransposed sequences in trans, including Alu, SVA and mRNA retrocopies and also full-length and 5'-inverted L1, were not associated with short TSDs, indicating the difference of 5' attachment between retrotransposition in cis and retrotransposition in trans. Target sequence analysis suggested that short TSDs were generated in an L1 endonuclease-dependent manner. The MHs were longer for 5'-inverted L1 than for 5'-truncated L1, indicating less dependence on annealing in 5'-truncated L1 insertions., Conclusions: The results suggest that insertions flanked by short TSDs occur more often coupled with the insertion of 5'-truncated L1 than with those of other types of insertions in vivo. The method used in this study can be used to characterize elements without any apparent boundary structures.

Chromosome elimination is a process in which some chromatins are discarded from the presumptive somatic cells during early embryogenesis. Eliminated chromatins in hagfish generally consist of repetitive sequences, and they are highly heterochromatinized in germ cells. In this study, we characterized four novel eliminated DNA families, EEPs1-4, from the Taiwanese hagfish Paramyxine sheni. Sequences of these four elements occupied 20-27% of eliminated DNA in total, and each family was arranged mainly in tandem in the germline genome with high copy numbers. Although most of these elements were eliminated, a minor fraction remained in somatic cells. Some eliminated DNA families are shared as eliminated sequences between Eptatretidae and Myxinidae. Fluorescence in situ hybridization (FISH) of these elements showed that not only heterochromatic chromosomes but also both ends of euchromatic chromosomes in germ cells are absent in somatic cells of P. sheni. It strongly suggests that chromosome terminus elimination, in addition to whole chromosome elimination, contributes to somatic chromosome differentiation. Telomere-FISH further showed that chromosome fragmentation and the subsequent de novo addition of telomeric repeats are the likely mechanisms underlying chromosome terminus elimination. These characteristics make it indispensable to study the evolution and mechanisms underlying chromosome elimination in hagfish.

Human long interspersed nuclear element-1 (L1) occupies one-sixth of our genome and has contributed to genome evolution in various ways. Approximately 10% of human L1 copies are composed of two L1 segments; the 5' segment and 3' segment are in head-to-head (i.e., 5'-inverted) orientation. Besides mediating their own retrotransposition, L1 has the ability to mobilize mRNA "in trans," and the number of retrotransposed mRNA sequences (retrocopies) is estimated to be >6,000. In this study, we identified 48 human-specific retrocopies and 95 chimpanzee-specific retrocopies by comparing the human and chimpanzee genomes. Among these retrocopies, 12 were 5'-inverted. The characteristics of these 5'-inverted retrocopies were similar to those of 5'-inverted L1 copies, indicating that the 5' inversion is generated by the same mechanism. With these findings, we examined the possibility that 5' inversion of the retrocopy generates a new gene that codes for a peptide with a different N terminus. We identified several potential 5'-inverted retrogenes, including those of thymopoietin beta (TMPO) and eukaryotic translation initiation factor 3 subunit 5 (EIF3F). The most interesting candidate was the 5'-inverted retrocopy of small nuclear ribonucleoprotein polypeptide N (SNRPN). This retrocopy was transcribed in the reverse orientation in several organs, had multiple transcript variants, and encoded a protein containing a peptide fragment derived from the N-terminal portion of SNRPN. Our results suggest that mRNA retrotransposition coupled with 5' inversion may be a mechanism to generate new genes distinct from parental genes.

Retroelements, elements encoding reverse transcriptase (RT), are ubiquitous in eukaryotes and have a great influence on the evolution of our genome. Detailed information is available on eukaryotic retroelements; however, prokaryotic retroelements are poorly understood. Recently, new types of eukaryotic retroelements were characterized on the basis of their gene composition and their phylogenetic positions. Here we performed a systematic survey to identify novel types of prokaryotic retroelements by analyzing gene neighborhood and protein architecture. We found novel types of gene combination and examined whether they represent actual retroelements. Five monophyletic groups were identified that were distinct from characterized prokaryotic retroelements, showed specific gene combination, were distributed patchily, and included at least 1 example of recent integration. These results strongly indicated the frequent horizontal transfer of these elements. One group encoded DNA polymerase A. A possible function of DNA polymerase A in the life cycle of retroelements is catalyzing second-strand cDNA synthesis, which is DNA polymerization performed using a DNA template not an RNA template. Another group encoded both bacterial primase and carbon-nitrogen hydrolase. Primase is likely to synthesize primers to initiate reverse transcription. Two other groups also encoded carbon-nitrogen hydrolase as a fusion protein with RT. It is difficult to speculate on the function of hydrolase in the life cycle of retroelements. The last group encoded dual RT proteins, which are likely to form heterodimers during replication. The protein sets of these 5 groups of prokaryotic retroelements were completely different from those of eukaryotic retroelements, indicating that the survival constraints of prokaryotic elements were distinct from those of eukaryotic elements. It is likely that these prokaryotic retroelements are maintained as extrachromosomal DNA or RNA or are accidentally integrated into genomes. Our findings presented the possibility that many types of extrachromosomal prokaryotic retroelements remain to be characterized. In addition, we found 8 RT genes were associated with clustered regularly interspaced short palindrome repeats (CRISPRs) of the CRISPR-Cas system. These RT genes are likely to work in immunity against RNA phages via cDNA synthesis.

Ribosomal RNA genes are abundant repetitive sequences in most eukaryotes. Ribosomal DNA (rDNA) contains many insertions derived from mobile elements including non-long terminal repeat (non-LTR) retrotransposons. R2 is the well-characterized 28S rDNA-specific non-LTR retrotransposon family that is distributed over at least 4 bilaterian phyla. R2 is a large family sharing the same insertion specificity and classified into 4 clades (R2-A, -B, -C, and -D) based on the N-terminal domain structure and the phylogeny. There is no observation of horizontal transfer of R2; therefore, the origin of R2 dates back to before the split between protostomes and deuterostomes. Here, we in silico identified 1 R2 element from the sea anemone Nematostella vectensis and 2 R2-like retrotransposons from the hydrozoan Hydra magnipapillata. R2 from N. vectensis was inserted into the 28S rDNA like other R2, but the R2-like elements from H. magnipapillata were inserted into the specific sequence in the highly conserved region of the 18S rDNA. We designated the Hydra R2-like elements R8. R8 is inserted at 37 bp upstream from R7, another 18S rDNA-specific retrotransposon family. There is no obvious sequence similarity between targets of R2 and R8, probably because they recognize long DNA sequences. Domain structure and phylogeny indicate that R2 from N. vectensis is the member of the R2-D clade, and R8 from H. magnipapillata belongs to the R2-A clade despite its different sequence specificity. These results suggest that R2 had been generated before the split between cnidarians and bilaterians and that R8 is a retrotransposon family that changed its target from the 28S rDNA to the 18S rDNA.

Chromosomal ends of most eukaryotes are composed of simple telomeric repeats. Arthropod telomeres are generally constituted by TTAGG pentanucleotide repeats; however, some insect species including Drosophila melanogaster do not have telomeric repeats. In contrast, the domestic silkworm Bombyx mori contains TTAGG-type telomeric repeats, but the telomerase activity has not been detected in all investigated tissues. To search for a cause of unusual telomere structure in insects, we here identified telomerase reverse transcriptase (TERT) subunit from the domestic silkworm B. mori and the flour beetle Tribolium castaneum. This is the first report of telomerase genes from arthropods. The domestic silkworm TERT gene (BmoTERT) and the flour beetle TERT gene (TcasTERT) both did not have the N-terminal GQ motif. Comparison between cDNA and genomic DNA of BmoTERT revealed that it includes no introns. BmoTERT contains five ATG codons in its 5'UTR, which could reduce the translation of BmoTERT proteins. Also, Northern hybridization indicated that BmoTERT is transcribed at a very low level. These unique features of BmoTERT possibly explain the undetectable Bombyx telomerase activity.

R2 is a non-long-terminal-repeat (LTR) retrotransposon that inserts specifically into 28S rDNA. R2 has been identified in many species of arthropods and three species of chordates. R2 may be even more widely distributed in animals, and its origin may be traceable to early animal evolution. In this study, we identified R2 elements in medaka fish, White Cloud Mountain minnow, Reeves' turtle, hagfish, sea lilies, and some arthropod species, using degenerate polymerase chain reaction methods. We also identified two R2 elements from the public genomic sequence database of the bloodfluke Schistosoma mansoni. One of the two bloodfluke R2 elements has two zinc-finger motifs at the N-terminus; this differs from other known R2 elements, which have one or three zinc-finger motifs. Phylogenetic analysis revealed that the whole phylogeny of R2 can be divided into 11 parts (subclades), in which the local R2 phylogeny and the corresponding host phylogeny are consistent. Divergence-versus-age analysis revealed that there is no reliable evidence for the horizontal transfer of R2 but supports the proposition that R2 has been vertically transferred since before the divergence of the deuterostomes and protostomes. The seeming inconsistency between the R2 phylogeny and the phylogeny of their hosts is due to the existence of paralogous lineages. The number of N-terminal zinc-finger motifs is consistent with the deep phylogeny of R2 and indicates that the common ancestor of R2 had three zinc-finger motifs at the N-terminus. This study revealed the long-term vertical inheritance and the ancient origin of sequence specificity of R2, both of which seem applicable to some other non-LTR retrotransposons.

Most eukaryotic cellular mRNAs are monocistronic; however, many retroviruses and long terminal repeat (LTR) retrotransposons encode multiple proteins on a single RNA transcript using ribosomal frameshifting. Non-long terminal repeat (non-LTR) retrotransposons are considered the ancestor of LTR retrotransposons and retroviruses, but their translational mechanism of bicistronic RNA remains unknown. We used a baculovirus expression system to produce a large amount of the bicistronic RNA of SART1, a non-LTR retrotransposon of the silkworm, and were able to detect the second open reading frame protein (ORF2) by Western blotting. The ORF2 protein was translated as an independent protein, not as an ORF1-ORF2 fusion protein. We revealed by mutagenesis that the UAAUG overlapping stop-start codon and the downstream RNA secondary structure are necessary for efficient ORF2 translation. Increasing the distance between the ORF1 stop codon and the ORF2 start codon decreased translation efficiency. These results are different from the eukaryotic translation reinitiation mechanism represented by the yeast GCN4 gene, in which the probability of reinitiation increases as the distance between the two ORFs increases. The translational mechanism of SART1 ORF2 is analogous to translational coupling observed in prokaryotes and viruses. Our results indicate that translational coupling is a general mechanism for bicistronic RNA translation.

Retrotransposons commonly encode a reverse transcriptase (RT), but other functional domains are variable. The acquisition of new domains is the dominant evolutionary force that brings structural variety to retrotransposons. Non-long-terminal-repeat (non-LTR) retrotransposons are classified into two groups by their structure. Early branched non-LTR retrotransposons encode a restriction-like endonuclease (RLE), and recently branched non-LTR retrotransposons encode an apurinic/apyrimidinic endonuclease-like endonuclease (APE). In this study, we report a novel non-LTR retrotransposon family Dualen, identified from the Chlamydomonas reinhardtii genome. Dualen encodes two endonucleases, RLE and APE, with RT, ribonuclease H, and cysteine protease. Phylogenetic analyses of the RT domains revealed that Dualen is positioned at the midpoint between the early-branched and the recently branched groups. In the APE tree, Dualen was branched earlier than the I group and the Jockey group. The ribonuclease H domains among the Dualen family and other non-LTR retrotransposons are monophyletic. Phylogenies of three domains revealed the monophyly of the Dualen family members. The domain structure and the phylogeny of each domain imply that Dualen is a retrotransposon conserving the domain structure just after the acquisition of APE. From these observations, we discuss the evolution of domain structure of non-LTR retrotransposons.

Most insects have telomeres that consist of pentanucleotide (TTAGG) telomeric repeats, which are synthesized by telomerase. However, all species in Diptera so far examined and several species in other orders of insect have lost the (TTAGG)n repeats, suggesting that some of them recruit telomerase-independent telomere maintenance. The silkworm, Bombyx mori, retains the TTAGG motifs in the chromosomal ends but expresses quite a low level of telomerase activity in all stages of various tissues. Just proximal to a 6-8-kb stretch of the TTAGG repeats in B. mori, more than 1000 copies of non-LTR retrotransposons, designated TRAS and SART families, occur among the telomeric repeats and accumulate. TRAS and SART are abundantly transcribed and actively retrotransposed into TTAGG telomeric repeats in a highly sequence-specific manner. They have three possible mechanisms to ensure specific integration into the telomeric repeats. This article focuses on the telomere structure and telomere-specific non-LTR retrotransposons in B. mori and discusses the mechanisms for telomere maintenance in this insect.

Non-long-terminal-repeat (non-LTR) retrotransposons amplify their copies by reverse transcribing mRNA from the 3' end, but the initial processes of reverse transcription are still unclear. We have shown that a telomere-specific non-LTR retrotransposon of the silkworm, SART1, requires the 3' untranslated region (3' UTR) for retrotransposition. With an in vivo retrotransposition assay, we identified several novel motifs within the 3' UTR involved in precise and efficient reverse transcription. Of 461 nucleotides (nt) of the 3' UTR, the central region, from nt 163 to nt 295, was essential for SART1 retrotransposition. Of five putative stem-loops formed in RNA for the SART1 3' UTR, the second stem-loop (nt 159 to 221) is included in this region. Loss of the 3' region (nt 296 to 461) in the 3' UTR and the poly(A) tract resulted in decreased and inaccurate reverse transcription, which starts mostly from several telomeric repeat-like GGUU sequences just downstream of the second stem-loop. These results suggest that short telomeric repeat-like sequences in the 3' UTR anneal to the bottom strand of (TTAGG)(n) repeats. We also demonstrated that the mRNA for green fluorescent protein (GFP) could be retrotransposed into telomeric repeats when the GFP coding region is fused with the SART1 3' UTR and SART1 open reading frame proteins are supplied in trans.

Although most LINEs (long interspersed nuclear elements), which are autonomous non-long-terminal-repeat retrotransposons, are inserted throughout the host genome, three groups of LINEs, the early-branched group, the Tx group, and the R1 clade, are inserted into specific sites within the target sequence. We previously characterized the sequence specificity of the R1 clade elements. In this study, we screened the other two groups of sequence-specific LINEs from public DNA databases, reconstructed elements from fragmented sequences, identified their target sequences, and analyzed them phylogenetically. We characterized 13 elements in the early-branched group and 13 in the Tx group. In the early-branched group, we identified R2 elements from sea squirts and zebrafish in this study, although R2 has not been characterized outside the arthropod group to date. This is the first evidence of cross-phylum distribution of sequence-specific LINEs. The Dong element also occurs across phyla, among arthropods and mollusks. In the Tx group, we characterized five novel sequence-specific families: Kibi for TC repeats, Koshi for TTC repeats, Keno for the U2 snRNA gene, Dewa for the tRNA tandem arrays, and Mutsu for the 5S rRNA gene. Keno and Mutsu insert into the highly conserved region within small RNA genes and destroy the targets. Several copies of Dewa insert different positions of tRNA tandem array, which indicates a certain "site specifier" other than sequence-specific endonuclease. In all three groups, LINEs specific for the rRNA genes or microsatellites can occur as multiple families in one organism. This indicates that the copy number of a target sequence is the primary factor to restrict the variety of sequence specificity of LINEs.

Although most non-long terminal repeat (non-LTR) retrotransposons are inserted throughout the host genome, many non-LTR elements in the R1 clade are inserted into specific sites within the target sequence. Four R1 clade families have distinct target specificity: R1 and RT insert into specific sites of 28S rDNA, and TRAS and SART insert into different sites within the (TTAGG)(n) telomeric repeats. To study the evolutionary history of target specificity of R1-clade retrotransposons, we have screened extensively novel representatives of the clade from various insects by in silico and degenerate polymerase chain reaction (PCR) cloning. We found four novel sequence-specific elements; Waldo (WaldoAg1, 2, and WaldoFs1) inserts into ACAY repeats, Mino (MinoAg1) into AC repeats, R6 into another specific site of the 28S rDNA, and R7 into a specific site of the 18S rDNA. In contrast, several elements (HOPE, WISHBm1, HidaAg1, NotoAg1, KagaAg1, Ha1Fs1) lost target sequence specificity, although some of them have preferred target sequences. Phylogenetic trees based on the RT and EN domains of each element showed that (1) three rDNA-specific elements, RT, R6, and R7, diverged from Waldo; (2) the elements having similar target sequences are phylogenetically related; and (3) the target specificity in the R1 clade was obtained once and thereafter altered and lost several times independently. These data indicate that the target specificity in R1 clade retroelements has changed during evolution and is more divergent than has been speculated so far.

Telomeres of most insects are composed of simple (TTAGG) (n) repeats that are synthesized by telomerase. However, in some dipteran insects such as Drosophila melanogaster, (TTAGG) (n) repeats or telomerase activity has not been detected. Although telomere structure is well documented in Diptera and Lepidoptera, very limited information is available on lower insect groups. To understand general aspects of telomere function and evolution in insects, we endeavored to characterize structures of the telomeric and subtelomeric regions in a lower insect, the Taiwan cricket, Teleogryllus taiwanemma. FISH analysis of this insect's chromosomes demonstrated (TTAGG) (n) repeat elements in all distal ends. Just proximal to the telomeric repeats, the highly conserved 9-kb long terminal unit (LTU) sequences are tandemly repeated. These were observed in four of six chromosomes, three autosomal ends, and one X-chromosomal end. LTU sequences represent about 0.2% of the T. taiwanemma genome. Each LTU contains a core (TTAGG)(8)-like sequence (TRLS) and five types of conserved sequences-ST (short telomere associated), J (joint), X, SR (satellite sequence rich), and Y-which vary in length from about 150 bp to 2.7 kb. The LTU sequence is defined as ST-J-TRLS-SR-X-Y-X-Y-X. Most LTU regions may be derived from the ancestral common sequence, which is observed in ST regions six times and at many other LTU sites. We could not find the LTU-like sequence in three other crickets including the closest species, T. emma, suggesting that the LTU in T. taiwanemma has been rapidly amplified in subtelomeric regions through recent evolutional events. It is also suggested that the highly conserved structure of the LTU is maintained by recombination and may contribute to telomere elongation, as seen in dipteran insects.