Genomic analysis of Caenorhabditis elegans reveals ancient families of retroviral-like elements

The CERV protein of Cer1, a C. elegans LTR retrotransposon, is required for nuclear export of viral genomic RNA and can form giant nuclear rods

Retroviruses and closely related LTR retrotransposons export full-length, unspliced genomic RNA (gRNA) for packaging into virions and to serve as the mRNA encoding GAG and POL polyproteins. Because gRNA often includes splice acceptor and donor sequences used to splice viral mRNAs, retroelements must overcome host mechanisms that retain intron-containing RNAs in the nucleus. Here we examine gRNA expression in Cer1, an LTR retrotransposon in C. elegans which somehow avoids silencing and is highly expressed in germ cells. Newly exported Cer1 gRNA associates rapidly with the Cer1 GAG protein, which has structural similarity with retroviral GAG proteins. gRNA export requires CERV (C. elegans regulator of viral expression), a novel protein encoded by a spliced Cer1 mRNA. CERV phosphorylation at S214 is essential for gRNA export, and phosphorylated CERV colocalizes with nuclear gRNA at presumptive sites of transcription. By electron microscopy, tagged CERV proteins surround clusters of distinct, linear fibrils that likely represent gRNA molecules. Single fibrils, or groups of aligned fibrils, also localize near nuclear pores. During the C. elegans self-fertile period, when hermaphrodites fertilize oocytes with their own sperm, CERV concentrates in two nuclear foci that are coincident with gRNA. However, as hermaphrodites cease self-fertilization, and can only produce cross-progeny, CERV undergoes a remarkable transition to form giant nuclear rods or cylinders that can be up to 5 microns in length. We propose a novel mechanism of rod formation, in which stage-specific changes in the nucleolus induce CERV to localize to the nucleolar periphery in flattened streaks of protein and gRNA; these streaks then roll up into cylinders. The rods are a widespread feature of Cer1 in wild strains of C. elegans, but their function is not known and might be limited to cross-progeny. We speculate that the adaptive strategy Cer1 uses for the identical self-progeny of a host hermaphrodite might differ for heterozygous cross-progeny sired by males. For example, mating introduces male chromosomes which can have different, or no, Cer1 elements.

Non-canonical nematode endogenous retroviruses resulting from RNA virus glycoprotein gene capture by a metavirus

Reverse-transcribing retroviruses exist as horizontally transmitted infectious agents or vertically transmitted endogenous retroviruses (ERVs) resident in eukaryotic genomes, and they are phylogenetically related to the long terminal repeat (LTR) class of retrotransposons. ERVs and retrotransposons are often distinguished only by the presence or absence of a gene encoding the envelope glycoprotein (env). Endogenous elements of the virus family Metaviridae include the insect-restricted Errantivirus genus of ERVs, for which some members possess env, and the pan-eukaryotic Metavirus genus that lacks an envelope glycoprotein gene. Here we report a novel Nematoda endogenous retrovirus (NERV) clade with core retroviral genes arranged uniquely as a continuous gag-env-pro-pol ORF. Reverse transcriptase sequences were phylogenetically related to metaviruses, but envelope glycoprotein sequences resembled those of the Nyamiviridae and Chrysoviridae RNA virus families, suggesting env gene capture during host cell infection by an RNA virus. NERVs were monophyletic, restricted to the nematode subclass Chromadoria, and included additional ORFs for a small hypothetical protein or a large Upf1-like RNA-dependent AAA-ATPase/helicase indicative of viral transduction of a host gene. Provirus LTR identity, low copy number, ORF integrity and segregation of three loci in Meloidogyne incognita, taken together with detection of NERV transcriptional activity, support potential infectivity of NERVs, along with their recent emergence and integration. Altogether, NERVs constitute a new and distinct Metaviridae lineage demonstrating retroviral evolution through sequential heterologous gene capture events.

A bioactive phlebovirus-like envelope protein in a hookworm endogenous virus

Endogenous viral elements (EVEs), accounting for 15% of our genome, serve as a genetic reservoir from which new genes can emerge. Nematode EVEs are particularly diverse and informative of virus evolution. We identify Atlas virus-an intact retrovirus-like EVE in the human hookworm Ancylostoma ceylanicum, with an envelope protein genetically related to GN-GC glycoproteins from the family Phenuiviridae. A cryo-EM structure of Atlas GC reveals a class II viral membrane fusion protein fold not previously seen in retroviruses. Atlas GC has the structural hallmarks of an active fusogen. Atlas GC trimers insert into membranes with endosomal lipid compositions and low pH. When expressed on the plasma membrane, Atlas GC has cell-cell fusion activity. With its preserved biological activities, Atlas GC has the potential to acquire a cellular function. Our work reveals structural plasticity in reverse-transcribing RNA viruses.

Genome-wide characterization of LTR retrotransposons in the non-model deep-sea annelid Lamellibrachia luymesi

Background

Long Terminal Repeat retrotransposons (LTR retrotransposons) are mobile genetic elements composed of a few genes between terminal repeats and, in some cases, can comprise over half of a genome’s content. Available data on LTR retrotransposons have facilitated comparative studies and provided insight on genome evolution. However, data are biased to model systems and marine organisms, including annelids, have been underrepresented in transposable elements studies. Here, we focus on genome of Lamellibrachia luymesi, a vestimentiferan tubeworm from deep-sea hydrocarbon seeps, to gain knowledge of LTR retrotransposons in a deep-sea annelid.

Results

We characterized LTR retrotransposons present in the genome of L. luymesi using bioinformatic approaches and found that intact LTR retrotransposons makes up about 0.1% of L. luymesi genome. Previous characterization of the genome has shown that this tubeworm hosts several known LTR-retrotransposons. Here we describe and classify LTR retrotransposons in L. luymesi as within the Gypsy, Copia and Bel-pao superfamilies. Although, many elements fell within already recognized families (e.g., Mag, CSRN1), others formed clades distinct from previously recognized families within these superfamilies. However, approximately 19% (41) of recovered elements could not be classified. Gypsy elements were the most abundant while only 2 Copia and 2 Bel-pao elements were present. In addition, analysis of insertion times indicated that several LTR-retrotransposons were recently transposed into the genome of L. luymesi, these elements had identical LTR’s raising possibility of recent or ongoing retrotransposon activity.

Conclusions

Our analysis contributes to knowledge on diversity of LTR-retrotransposons in marine settings and also serves as an important step to assist our understanding of the potential role of retroelements in marine organisms. We find that many LTR retrotransposons, which have been inserted in the last few million years, are similar to those found in terrestrial model species. However, several new groups of LTR retrotransposons were discovered suggesting that the representation of LTR retrotransposons may be different in marine settings. Further study would improve understanding of the diversity of retrotransposons across animal groups and environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07749-1.

Mollusc genomes reveal variability in patterns of LTR-retrotransposons dynamics

Background

The three superfamilies of Long Terminal Repeat (LTR) retrotransposons are a widespread kind of transposable element and a major factor in eukaryotic genome evolution. In metazoans, recent studies suggested that Copia LTR-retrotransposons display specific dynamic compared to the more abundant and diverse Gypsy elements. Indeed, Copia elements show a relative scarcity and the prevalence of only a few clades in specific hosts. Thus, BEL/Pao seems to be the second most abundant superfamily. However, the generality of these assumptions remains to be assessed. Therefore, we carried out the first large-scale comparative genomic analysis of LTR-retrotransposons in molluscs. The aim of this study was to analyse the diversity, copy numbers, genomic proportions and distribution of LTR-retrotransposons in a large host phylum.

Results

We compare nine genomes of molluscs and further added LTR-retrotransposons sequences detected in databases for 47 additional species. We identified 1709 families, which enabled us to define 31 clades. We show that clade richness was highly dependent on the considered superfamily. We found only three Copia clades, including GalEa and Hydra which appear to be widely distributed and highly dominant as they account for 96% of the characterised Copia elements. Among the seven BEL/Pao clades identified, Sparrow and Surcouf are characterised for the first time. We find no BEL or Pao elements, but the rare clades Dan and Flow are present in molluscs. Finally, we characterised 21 Gypsy clades, only five of which had been previously described, the C-clade being the most abundant one. Even if they are found in the same number of host species, Copia elements are clearly less abundant than BEL/Pao elements in copy number or genomic proportions, while Gypsy elements are always the most abundant ones whatever the parameter considered.

Conclusions

Our analysis confirms the contrasting dynamics of Copia and Gypsy elements in metazoans and indicates that BEL/Pao represents the second most abundant superfamily, probably reflecting an intermediate dynamic. Altogether, the data obtained in several taxa highly suggest that these patterns can be generalised for most metazoans. Finally, we highlight the importance of using database information in complement of genome analyses when analyzing transposable element diversity.

Electronic supplementary material

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

The spatial and temporal dynamics of nuclear RNAi-targeted retrotransposon transcripts in Caenorhabditis elegans

Nuclear RNA interference provides a unique approach to the study of RNA-mediated transgenerational epigenetic inheritance. A paradox in the field is that expression of target loci is necessary for the initiation and maintenance of their silencing. How expression and repression are coordinated during animal development is poorly understood. To resolve this gap, we took imaging, deep-sequencing and genetic approaches towards delineating the developmental regulation and subcellular localization of RNA transcripts of two representative endogenous targets, the LTR retrotransposons Cer3 and Cer8. By examining wild-type worms and a collection of mutant strains, we found that the expression and silencing cycle of Cer3 and Cer8 is coupled with embryonic and germline development. Strikingly, endogenous targets exhibit a hallmark of nuclear enrichment of their RNA transcripts. In addition, germline and somatic repressions of Cer3 have different genetic requirements for three heterochromatin enzymes, MET-2, SET-25 and SET-32, in conjunction with the nuclear Argonaute protein HRDE-1. These results provide the first comprehensive cellular and developmental characterization of nuclear RNAi activities throughout the animal reproductive cycle.

The first complete Mag family retrotransposons discovered in Drosophila

A retrotransposon of the Mag family was found in the Drosophila simulans genome for the first time. We also identified novel transposable elements representing the Mag family in seven Drosophila species. The high similarity between the 3' and 5' long terminal repeats in the found copies of transposable elements indicates that their retrotransposition has occurred relatively recently. Thus, the Mag family of retrotransposons is quite common for the genus Drosophila.

Genome Filtering for New DNA Biomarkers of Loa loa Infection Suitable for Loop-Mediated Isothermal Amplification

Loa loa infections have emerged as a serious public health problem in patients co-infected with Onchocerca volvulus or Wuchereria bancrofti because of severe adverse neurological reactions after treatment with ivermectin. Accurate diagnostic tests are needed for careful mapping in regions where mass drug administration is underway. Loop-mediated isothermal amplification (LAMP) has become a widely adopted screening method because of its operational simplicity, rapidity and versatility of visual detection readout options. Here, we present a multi-step bioinformatic pipeline to generate diagnostic candidates suitable for LAMP and experimentally validate this approach using one of the identified candidates to develop a species-specific LAMP assay for L. loa. The pipeline identified ~140 new L. loa specific DNA repeat families as putative biomarkers of infection. The consensus sequence of one family, repeat family 4 (RF4), was compiled from ~ 350 sequences dispersed throughout the L. loa genome and maps to a L. loa-specific region of the long terminal repeats found at the boundaries of Bel/Pao retrotransposons. PCR and LAMP primer sets targeting RF4 specifically amplified L. loa but not W. bancrofti, O. volvulus, Brugia malayi, human or mosquito DNA. RF4 LAMP detects the DNA equivalent of one microfilaria (100 pg) in 25-30 minutes and as little as 0.060 pg of L. loa DNA (~1/1600th of a microfilaria) purified from spiked blood samples in approximately 50 minutes. In summary, we have successfully employed a bioinformatic approach to mine the L. loa genome for species-specific repeat families that can serve as new DNA biomarkers for LAMP. The RF4 LAMP assay shows promise as a field tool for the implementation and management of mass drug administration programs and warrants further testing on clinical samples as the next stage in development towards this goal.

Complex coding of endogenous siRNA, transcriptional silencing and H3K9 methylation on native targets of germline nuclear RNAi in C. elegans

Background

Small RNA-guided transcriptional silencing (nuclear RNAi) is fundamental to genome integrity and epigenetic inheritance. Despite recent progress in identifying the capability and genetic requirements for nuclear RNAi in Caenorhabditis elegans, the natural targets and cellular functions of nuclear RNAi remain elusive.

Methods

To resolve this gap, we coordinately examined the genome-wide profiles of transcription, histone H3 lysine 9 methylation (H3K9me) and endogenous siRNAs of a germline nuclear Argonaute (hrde-1/wago-9) mutant and identified regions on which transcription activity is markedly increased and/or H3K9me level is markedly decreased relative to wild type animals.

Results

Our data revealed a distinct set of native targets of germline nuclear RNAi, with the H3K9me response exhibiting both overlapping and non-overlapping distribution with the transcriptional silencing response. Interestingly LTR retrotransposons, but not DNA transposons, are highly enriched in the targets of germline nuclear RNAi. The genomic distribution of the native targets is highly constrained, with >99% of the identified targets present in five autosomes but not in the sex chromosome. By contrast, HRDE-1-associated small RNAs correspond to all chromosomes. In addition, we found that the piRNA pathway is not required for germline nuclear RNAi activity on native targets.

Conclusion

Germline nuclear RNAi in C. elegans is required to silence retrotransposons but not DNA transposon. Transcriptional silencing and H3K9me can occur independently of each other on the native targets of nuclear RNAi in C. elegans. Our results rule out a simple model in which nuclear Argonaute protein-associated-small RNAs are sufficient to trigger germline nuclear RNAi responses. In addition, the piRNA pathway and germline nuclear RNAi are specialized to target different types of foreign genetic elements for genome surveillance in C. elegans.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1157) contains supplementary material, which is available to authorized users.

LTR retroelements are intrinsic components of transcriptional networks in frogs

Background

LTR retroelements (LTR REs) constitute a major group of transposable elements widely distributed in eukaryotic genomes. Through their own mechanism of retrotranscription LTR REs enrich the genomic landscape by providing genetic variability, thus contributing to genome structure and organization. Nonetheless, transcriptomic activity of LTR REs still remains an obscure domain within cell, developmental, and organism biology.

Results

Here we present a first comparative analysis of LTR REs for anuran amphibians based on a full depth coverage transcriptome of the European pool frog, Pelophylax lessonae, the genome of the African clawed frog, Silurana tropicalis (release v7.1), and additional transcriptomes of S. tropicalis and Cyclorana alboguttata. We identified over 1000 copies of LTR REs from all four families (Bel/Pao, Ty1/Copia, Ty3/Gypsy, Retroviridae) in the genome of S. tropicalis and discovered transcripts of several of these elements in all RNA-seq datasets analyzed. Elements of the Ty3/Gypsy family were most active, especially Amn-san elements, which accounted for approximately 0.27% of the genome in Silurana. Some elements exhibited tissue specific expression patterns, for example Hydra1.1 and MuERV-like elements in Pelophylax. In S. tropicalis considerable transcription of LTR REs was observed during embryogenesis as soon as the embryonic genome became activated, i.e. at midblastula transition. In the course of embryonic development the spectrum of transcribed LTR REs changed; during gastrulation and neurulation MuERV-like and SnRV like retroviruses were abundantly transcribed while during organogenesis transcripts of the XEN1 retroviruses became much more active.

Conclusions

The differential expression of LTR REs during embryogenesis in concert with their tissue-specificity and the protein domains they encode are evidence for the functional roles these elements play as integrative parts of complex regulatory networks. Our results support the meanwhile widely accepted concept that retroelements are not simple “junk DNA” or “harmful genomic parasites” but essential components of the transcriptomic machinery in vertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-626) contains supplementary material, which is available to authorized users.

The Tvv1 retrotransposon family is conserved between plant genomes separated by over 100 million years

Combining several different approaches, we have examined the structure, variability, and distribution of Tvv1 retrotransposons. Tvv1 is an unusual example of a low-copy retrotransposon metapopulation dispersed unevenly among very distant species and is promising for the development of molecular markers. Retrotransposons are ubiquitous throughout the genomes of the vascular plants, but individual retrotransposon families tend to be confined to the level of plant genus or at most family. This restricts the general applicability of a family as molecular markers. Here, we characterize a new plant retrotransposon named Tvv1_Sdem, a member of the Copia superfamily of LTR retrotransposons, from the genome of the wild potato Solanum demissum. Comparative analyses based on structure and sequence showed a high level of similarity of Tvv1_Sdem with Tvv1-VB, a retrotransposon previously described in the grapevine genome Vitis vinifera. Extending the analysis to other species by in silico and in vitro approaches revealed the presence of Tvv1 family members in potato, tomato, and poplar genomes, and led to the identification of full-length copies of Tvv1 in these species. We were also able to identify polymorphism in UTL sequences between Tvv1_Sdem copies from wild and cultivated potatoes that are useful as molecular markers. Combining different approaches, our results suggest that the Tvv1 family of retrotransposons has a monophyletic origin and has been maintained in both the rosids and the asterids, the major clades of dicotyledonous plants, since their divergence about 100 MYA. To our knowledge, Tvv1 represents an unusual plant retrotransposon metapopulation comprising highly similar members disjointedly dispersed among very distant species. The twin features of Tvv1 presence in evolutionarily distant genomes and the diversity of its UTL region in each species make it useful as a source of robust molecular markers for diversity studies and breeding.

Development of a polymerase chain reaction applicable to rapid and sensitive detection of Clonorchis sinensis eggs in human stool samples

Microscopic examination of eggs of parasitic helminths in stool samples has been the most widely used classical diagnostic method for infections, but tiny and low numbers of eggs in stool samples often hamper diagnosis of helminthic infections with classical microscopic examination. Moreover, it is also difficult to differentiate parasite eggs by the classical method, if they have similar morphological characteristics. In this study, we developed a rapid and sensitive polymerase chain reaction (PCR)-based molecular diagnostic method for detection of Clonorchis sinensis eggs in stool samples. Nine primers were designed based on the long-terminal repeat (LTR) of C. sinensis retrotransposon1 (CsRn1) gene, and seven PCR primer sets were paired. Polymerase chain reaction with each primer pair produced specific amplicons for C. sinensis, but not for other trematodes including Metagonimus yokogawai and Paragonimus westermani. Particularly, three primer sets were able to detect 10 C. sinensis eggs and were applicable to amplify specific amplicons from DNA samples purified from stool of C. sinensis-infected patients. This PCR method could be useful for diagnosis of C. sinensis infections in human stool samples with a high level of specificity and sensitivity.

C. elegans germ cells show temperature and age-dependent expression of Cer1, a Gypsy/Ty3-related retrotransposon

Virus-like particles (VLPs) have not been observed in Caenorhabditis germ cells, although nematode genomes contain low numbers of retrotransposon and retroviral sequences. We used electron microscopy to search for VLPs in various wild strains of Caenorhabditis, and observed very rare candidate VLPs in some strains, including the standard laboratory strain of C. elegans, N2. We identified the N2 VLPs as capsids produced by Cer1, a retrotransposon in the Gypsy/Ty3 family of retroviruses/retrotransposons. Cer1 expression is age and temperature dependent, with abundant expression at 15°C and no detectable expression at 25°C, explaining how VLPs escaped detection in previous studies. Similar age and temperature-dependent expression of Cer1 retrotransposons was observed for several other wild strains, indicating that these properties are common, if not integral, features of this retroelement. Retrotransposons, in contrast to DNA transposons, have a cytoplasmic stage in replication, and those that infect non-dividing cells must pass their genomic material through nuclear pores. In most C. elegans germ cells, nuclear pores are largely covered by germline-specific organelles called P granules. Our results suggest that Cer1 capsids target meiotic germ cells exiting pachytene, when free nuclear pores are added to the nuclear envelope and existing P granules begin to be removed. In pachytene germ cells, Cer1 capsids concentrate away from nuclei on a subset of microtubules that are exceptionally resistant to microtubule inhibitors; the capsids can aggregate these stable microtubules in older adults, which exhibit a temperature-dependent decrease in egg viability. When germ cells exit pachytene, the stable microtubules disappear and capsids redistribute close to nuclei that have P granule-free nuclear pores. This redistribution is microtubule dependent, suggesting that capsids that are released from stable microtubules transfer onto new, dynamic microtubules to track toward nuclei. These studies introduce C. elegans as a model to study the interplay between retroelements and germ cell biology.

Retrotransposons and retroviruses pose enormous threats to animal and plants because of their ability to insert into host genes. Retroelements that replicate in germ cells can, if left unchecked, expand exponentially in the host genome. C. elegans has proven to be an exceptional model system for studying many facets of cell and molecular biology, and the genome contains both retrotransposon and retroviral sequences. However, no virus-like particles have been observed in C. elegans germ cells. We show here that Cer1, an endogenous Gypsy/Ty3 class retrotransposon, is expressed at very high levels in C. elegans germ cells, but escaped detection in previous studies because its expression is both temperature and age dependent. These studies reveal new aspects of microtubule regulation in C. elegans that the retroelement appears to exploit to navigate the germ cell cytoplasm, and demonstrate the power of C. elegans for studying host/pathogen interactions in germ cell biology.

Viral ancestors of antiviral systems

All life must survive their corresponding viruses. Thus antiviral systems are essential in all living organisms. Remnants of virus derived information are also found in all life forms but have historically been considered mostly as junk DNA. However, such virus derived information can strongly affect host susceptibility to viruses. In this review, I evaluate the role viruses have had in the origin and evolution of host antiviral systems. From Archaea through bacteria and from simple to complex eukaryotes I trace the viral components that became essential elements of antiviral immunity. I conclude with a reexamination of the 'Big Bang' theory for the emergence of the adaptive immune system in vertebrates by horizontal transfer and note how viruses could have and did provide crucial and coordinated features.

Viral ancestors of antiviral systems

All life must survive their corresponding viruses. Thus antiviral systems are essential in all living organisms. Remnants of virus derived information are also found in all life forms but have historically been considered mostly as junk DNA. However, such virus derived information can strongly affect host susceptibility to viruses. In this review, I evaluate the role viruses have had in the origin and evolution of host antiviral systems. From Archaea through bacteria and from simple to complex eukaryotes I trace the viral components that became essential elements of antiviral immunity. I conclude with a reexamination of the ‘Big Bang’ theory for the emergence of the adaptive immune system in vertebrates by horizontal transfer and note how viruses could have and did provide crucial and coordinated features.

Novel transposable elements from Anopheles gambiae

Background

Transposable elements (TEs) are DNA sequences, present in the genome of most eukaryotic organisms that hold the key characteristic of being able to mobilize and increase their copy number within chromosomes. These elements are important for eukaryotic genome structure and evolution and lately have been considered as potential drivers for introducing transgenes into pathogen-transmitting insects as a means to control vector-borne diseases. The aim of this work was to catalog the diversity and abundance of TEs within the Anopheles gambiae genome using the PILER tool and to consolidate a database in the form of a hyperlinked spreadsheet containing detailed and readily available information about the TEs present in the genome of An. gambiae.

Results

Here we present the spreadsheet named AnoTExcel that constitutes a database with detailed information on most of the repetitive elements present in the genome of the mosquito. Despite previous work on this topic, our approach permitted the identification and characterization both of previously described and novel TEs that are further described in detailed.

Conclusions

Identification and characterization of TEs in a given genome is important as a way to understand the diversity and evolution of the whole set of TEs present in a given species. This work contributes to a better understanding of the landscape of TEs present in the mosquito genome. It also presents a novel platform for the identification, analysis, and characterization of TEs on sequenced genomes.

The source of self: genetic parasites and the origin of adaptive immunity

Stable colonization of the host by viruses (genetic parasites) can alter the systems of host identity and provide immunity against related viruses. To attain the needed stability, some viruses of prokaryotes (P1 phage) use a strategy called an addiction module. The linked protective and destructive gene functions of an addiction module insures both virus persistence but will also destroy cells that interrupt this module and thereby prevent infection by competitors. Previously, I have generalized this concept to also include persistent and lytic states of virus infection, which can be considered as a virus addiction module. Such states often involve defective viruses. In this report, I examine the origin of the adaptive immune system from the perspective of a virus addiction module. The likely role of both endogenous and exogenous retroviruses, DNA viruses, and their defective elements is considered in the origin of all the basal components of adaptive immunity (T-cell receptor, RAG-mediated gene rearrangement, clonal lymphocyte proliferation, antigen surface presentation, apoptosis, and education of immune cells). It is concluded that colonization by viruses and their defectives provides a more coherent explanation for the origin of adaptive immunity.

Network dynamics of eukaryotic LTR retroelements beyond phylogenetic trees

BACKGROUND

Sequencing projects have allowed diverse retroviruses and LTR retrotransposons from different eukaryotic organisms to be characterized. It is known that retroviruses and other retro-transcribing viruses evolve from LTR retrotransposons and that this whole system clusters into five families: Ty3/Gypsy, Retroviridae, Ty1/Copia, Bel/Pao and Caulimoviridae. Phylogenetic analyses usually show that these split into multiple distinct lineages but what is yet to be understood is how deep evolution occurred in this system.

RESULTS

We combined phylogenetic and graph analyses to investigate the history of LTR retroelements both as a tree and as a network. We used 268 non-redundant LTR retroelements, many of them introduced for the first time in this work, to elucidate all possible LTR retroelement phylogenetic patterns. These were superimposed over the tree of eukaryotes to investigate the dynamics of the system, at distinct evolutionary times. Next, we investigated phenotypic features such as duplication and variability of amino acid motifs, and several differences in genomic ORF organization. Using this information we characterized eight reticulate evolution markers to construct phenotypic network models.

CONCLUSION

The evolutionary history of LTR retroelements can be traced as a time-evolving network that depends on phylogenetic patterns, epigenetic host-factors and phenotypic plasticity. The Ty1/Copia and the Ty3/Gypsy families represent the oldest patterns in this network that we found mimics eukaryotic macroevolution. The emergence of the Bel/Pao, Retroviridae and Caulimoviridae families in this network can be related with distinct inflations of the Ty3/Gypsy family, at distinct evolutionary times. This suggests that Ty3/Gypsy ancestors diversified much more than their Ty1/Copia counterparts, at distinct geological eras. Consistent with the principle of preferential attachment, the connectivities among phenotypic markers, taken as network-represented combinations, are power-law distributed. This evidences an inflationary mode of evolution where the system diversity; 1) expands continuously alternating vertical and gradual processes of phylogenetic divergence with episodes of modular, saltatory and reticulate evolution; 2) is governed by the intrinsic capability of distinct LTR retroelement host-communities to self-organize their phenotypes according to emergent laws characteristic of complex systems.

REVIEWERS

This article was reviewed by Eugene V. Koonin, Eric Bapteste, and Enmanuelle Lerat (nominated by King Jordan).

Bioinformatic flowchart and database to investigate the origins and diversity of clan AA peptidases

BACKGROUND

Clan AA of aspartic peptidases relates the family of pepsin monomers evolutionarily with all dimeric peptidases encoded by eukaryotic LTR retroelements. Recent findings describing various pools of single-domain nonviral host peptidases, in prokaryotes and eukaryotes, indicate that the diversity of clan AA is larger than previously thought. The ensuing approach to investigate this enzyme group is by studying its phylogeny. However, clan AA is a difficult case to study due to the low similarity and different rates of evolution. This work is an ongoing attempt to investigate the different clan AA families to understand the cause of their diversity.

RESULTS

In this paper, we describe in-progress database and bioinformatic flowchart designed to characterize the clan AA protein domain based on all possible protein families through ancestral reconstructions, sequence logos, and hidden markov models (HMMs). The flowchart includes the characterization of a major consensus sequence based on 6 amino acid patterns with correspondence with Andreeva's model, the structural template describing the clan AA peptidase fold. The set of tools is work in progress we have organized in a database within the GyDB project, referred to as Clan AA Reference Database http://gydb.uv.es/gydb/phylogeny.php?tree=caard.

CONCLUSION

The pre-existing classification combined with the evolutionary history of LTR retroelements permits a consistent taxonomical collection of sequence logos and HMMs. This set is useful for gene annotation but also a reference to evaluate the diversity of, and the relationships among, the different families. Comparisons among HMMs suggest a common ancestor for all dimeric clan AA peptidases that is halfway between single-domain nonviral peptidases and those coded by Ty3/Gypsy LTR retroelements. Sequence logos reveal how all clan AA families follow similar protein domain architecture related to the peptidase fold. In particular, each family nucleates a particular consensus motif in the sequence position related to the flap. The different motifs constitute a network where an alanine-asparagine-like variable motif predominates, instead of the canonical flap of the HIV-1 peptidase and closer relatives.

Characterization of ten novel Ty1/copia-like retrotransposon families of the grapevine genome

Background

Retrotransposons make a significant contribution to the size, organization and genetic diversity of their host genomes. To characterize retrotransposon families in the grapevine genome (the fourth crop plant genome sequenced) we have combined two approaches: a PCR-based method for the isolation of RnaseH-LTR sequences with a computer-based sequence similarity search in the whole-genome sequence of PN40024.

Results

Supported by a phylogenic analysis, ten novel Ty1/copia families were distinguished in this study. To select a canonical reference element sequence from amongst the various insertions in the genome belonging to each retroelement family, the following screening criteria were adopted to identify the element sequence with: (1) perfect 5 bp-duplication of target sites, (2) the highest level of identity between 5' and 3'-LTR within a single insertion sequence, and (3) longest, un-interrupted coding capacity within the gag-pol ORF. One to eight copies encoding a single putatively functional gag-pol polyprotein were found for three families, indicating that these families could be still autonomous and active. For the others, no autonomous copies were identified. However, a subset of copies within the presumably non-autonomous families had perfect identity between their 5' and 3' LTRs, indicating a recent insertion event. A phylogenic study based on the sequence alignment of the region located between reverse transcriptase domains I and VII distinguished these 10 families from other plant retrotransposons. Including the previously characterized Ty1/copia-like grapevine retrotransposons Tvv1 and Vine 1 and the Ty3/gypsy-like Gret1 in this assessment, a total of 1709 copies were identified for the 13 retrotransposon families, representing 1.24% of the sequenced genome. The copy number per family ranged from 91–212 copies. We performed insertion site profiling for 8 out of the 13 retrotransposon families and confirmed multiple insertions of these elements across the Vitis genus. Insertional polymorphism analysis and dating of full-length copies based on their LTR divergence demonstrated that each family has a particular amplification history, with 71% of the identified copies being inserted within the last 2 million years.

Conclusion

The strategy we used efficiently delivered new Ty1/copia-like retrotransposon sequences, increasing the total number of characterized grapevine retrotrotransposons from 3 to 13. We provide insights into the representation and dynamics of the 13 families in the genome. Our data demonstrated that each family has a particular amplification pattern, with 7 families having copies recently inserted within the last 0.2 million year. Among those 7 families with recent insertions, three retain the capacity for activity in the grape genome today.

Relationships of gag-pol diversity between Ty3/Gypsy and Retroviridae LTR retroelements and the three kings hypothesis

BACKGROUND

The origin of vertebrate retroviruses (Retroviridae) is yet to be thoroughly investigated, but due to their similarity and identical gag-pol (and env) genome structure, it is accepted that they evolve from Ty3/Gypsy LTR retroelements the retrotransposons and retroviruses of plants, fungi and animals. These 2 groups of LTR retroelements code for 3 proteins rarely studied due to the high variability - gag polyprotein, protease and GPY/F module. In relation to 3 previously proposed Retroviridae classes I, II and II, investigation of the above proteins conclusively uncovers important insights regarding the ancient history of Ty3/Gypsy and Retroviridae LTR retroelements.

RESULTS

We performed a comprehensive study of 120 non-redundant Ty3/Gypsy and Retroviridae LTR retroelements. Phylogenetic reconstruction inferred based on the concatenated analysis of the gag and pol polyproteins shows a robust phylogenetic signal regarding the clustering of OTUs. Evaluation of gag and pol polyproteins separately yields discordant information. While pol signal supports the traditional perspective (2 monophyletic groups), gag polyprotein describes an alternative scenario where each Retroviridae class can be distantly related with one or more Ty3/Gypsy lineages. We investigated more in depth this evidence through comparative analyses performed based on the gag polyprotein, the protease and the GPY/F module. Our results indicate that contrary to the traditional monophyletic view of the origin of vertebrate retroviruses, the Retroviridae class I is a molecular fossil, preserving features that were probably predominant among Ty3/Gypsy ancestors predating the split of plants, fungi and animals. In contrast, classes II and III maintain other phenotypes that emerged more recently during Ty3/Gypsy evolution.

CONCLUSION

The 3 Retroviridae classes I, II and III exhibit phenotypic differences that delineate a network never before reported between Ty3/Gypsy and Retroviridae LTR retroelements. This new scenario reveals how the diversity of vertebrate retroviruses is polyphyletically recurrent into the Ty3/Gypsy evolution, i.e. older than previously thought. The simplest hypothesis to explain this finding is that classes I, II and III trace back to at least 3 Ty3/Gypsy ancestors that emerged at different evolutionary times prior to protostomes-deuterostomes divergence. We have called this "the three kings hypothesis" concerning the origin of vertebrate retroviruses.

How Athila retrotransposons survive in the Arabidopsis genome

BACKGROUND

Transposable elements are selfish genetic sequences which only occasionally provide useful functions to their host species. In addition, models of mobile element evolution assume a second type of selfishness: elements of different families do not cooperate, but they independently fight for their survival in the host genome.

RESULTS

We show that recombination events among distantly related Athila retrotransposons have led to the generation of new Athila lineages. Their pattern of diversification suggests that Athila elements survive in Arabidopsis by a combination of selfish replication and of amplification of highly diverged copies with coding potential. Many Athila elements are non-autonomous but still conserve intact open reading frames which are under the effect of negative, purifying natural selection.

CONCLUSION

The evolution of these mobile elements is far more complex than hitherto assumed. Strict selfish replication does not explain all the patterns observed.

The diversity of retrotransposons and the properties of their reverse transcriptases

A number of abundant mobile genetic elements called retrotransposons reverse transcribe RNA to generate DNA for insertion into eukaryotic genomes. Four major classes of retrotransposons are described here. First, the long-terminal-repeat (LTR) retrotransposons have similar structures and mechanisms to those of the vertebrate retroviruses. Genes that may enable these retrotransposons to leave a cell have been acquired by these elements in a number of animal and plant lineages. Second, the tyrosine recombinase retrotransposons are similar to the LTR retrotransposons except that they have substituted a recombinase for the integrase and recombine into the host chromosomes. Third, the non-LTR retrotransposons use a cleaved chromosomal target site generated by an encoded endonuclease to prime reverse transcription. Finally, the Penelope-like retrotransposons are not well understood but appear to also use cleaved DNA or the ends of chromosomes as primer for reverse transcription. Described in the second part of this review are the enzymatic properties of the reverse transcriptases (RTs) encoded by retrotransposons. The RTs of the LTR retrotransposons are highly divergent in sequence but have similar enzymatic activities to those of retroviruses. The RTs of the non-LTR retrotransposons have several unique properties reflecting their adaptation to a different mechanism of retrotransposition.

The Gypsy Database (GyDB) of mobile genetic elements

In this article, we introduce the Gypsy Database (GyDB) of mobile genetic elements, an in-progress database devoted to the non-redundant analysis and evolutionary-based classification of mobile genetic elements. In this first version, we contemplate eukaryotic Ty3/Gypsy and Retroviridae long terminal repeats (LTR) retroelements. Phylogenetic analyses based on the gag-pro-pol internal region commonly presented by these two groups strongly support a certain number of previously described Ty3/Gypsy lineages originally reported from reverse-transcriptase (RT) analyses. Vertebrate retroviruses (Retroviridae) are also constituted in several monophyletic groups consistent with genera proposed by the ICTV nomenclature, as well as with the current tendency to classify both endogenous and exogenous retroviruses by three major classes (I, II and III). Our inference indicates that all protein domains codified by the gag-pro-pol internal region of these two groups agree in a collective presentation of a particular evolutionary history, which may be used as a main criterion to differentiate their molecular diversity in a comprehensive collection of phylogenies and non-redundant molecular profiles useful in the identification of new Ty3/Gypsy and Retroviridae species. The GyDB project is available at http://gydb.uv.es.

Two novel elements (CFG1 and PYG1) of Mag lineage of Ty3/Gypsy retrotransposons from Zhikong scallop (Chlamys farreri) and Japanese scallop (Patinopecten yessoensis)

Two novel elements (CFG1 and PYG1) of Mag lineage of Ty3/Gypsy retrotransposons were cloned from Zhikong scallop (Chlamys farreri) and Japanese scallop (Patinopecten yessoensis). The total length of the CFG1 element is 4826 bp, including 5'-LTR (192 bp), the entire ORF (4047 bp) and 3'-LTR (189 bp). The entire ORFs of both CFG1 and PYG1 elements are composed of 1348 aa and do not have any frameshifts. Their closest relative is Jule element from the poeciliid fish (Xiphophorus maculatus). On average, the diploid genome of C. farreri contains approximately 84 copies of CFG1 elements. We summarize the major features of CFG1, PYG1 and other elements of Mag lineage of the Ty3/Gypsy group. mRNA expression of CFG1 element in larvae increases gradually before the gastrulae stage and decreases gradually afterward, whereas in adductor such expression in adductor muscle and digestive gland are lower than those in other tissues. Overall, mRNA expression of CFG1 element in the early larvae is significantly higher than that in adult tissues. In muscle tissue, while the promoter and partial GAG domain of CFG1 element are unmethylated, the partial RT domain is highly methylated. These results suggest that CFG1 expression may be controlled by a post-transcriptional gene silencing mechanism that is associated with coding-region (RT domain) methylation.

Both sense and antisense strands of the LTR of the Schistosoma mansoni Pao-like retrotransposon Sinbad drive luciferase expression

Long terminal repeat (LTR) retrotransposons, mobile genetic elements comprising substantial proportions of many eukaryotic genomes, are so named for the presence of LTRs, direct repeats about 250-600 bp in length flanking the open reading frames that encode the retrotransposon enzymes and structural proteins. LTRs include promotor functions as well as other roles in retrotransposition. LTR retrotransposons, including the Gypsy-like Boudicca and the Pao/BEL-like Sinbad elements, comprise a substantial proportion of the genome of the human blood fluke, Schistosoma mansoni. In order to deduce the capability of specific copies of Boudicca and Sinbad LTRs to function as promotors, these LTRs were investigated analytically and experimentally. Sequence analysis revealed the presence of TATA boxes, canonical polyadenylation signals, and direct inverted repeats within the LTRs of both the Boudicca and Sinbad retrotransposons. Inserted in the reporter plasmid pGL3, the LTR of Sinbad drove firefly luciferase activity in HeLa cells in its forward and inverted orientation. In contrast, the LTR of Boudicca did not drive luciferase activity in HeLa cells. The ability of the Sinbad LTR to transcribe in both its forward and inverted orientation represents one of few documented examples of bidirectional promotor function.

Ongoing loss of the tirant transposable element in natural populations of Drosophila simulans

Tirant is a long terminal repeat (LTR) retrotransposon with an average of 11 insertion sites on the chromosome arms of Drosophila melanogaster flies collected from natural populations worldwide. In the sibling species Drosophila simulans, tirant is found only in African populations, which harbor a few insertion sites (1 to 5) on the chromosome arms, although some tirant sequences are present in the heterochromatin of most populations. This distribution in D. simulans reflects either the recent genomic invasion of African populations by a new variant of tirant, or a loss of tirant from the entire species apart from some sequence relics still present in Africa. In an attempt to clarify the situation, we focused on the LTR-UTR region of tirant copies from various populations of both D. melanogaster and D. simulans. We found two distinct types of regulatory region: one type was present in both D. melanogaster and D. simulans, and the other was present only in D. simulans. Copies of this latter type of tirant were transcriptionally inactive in gonads. Here we propose that the present day distribution of tirant in D. simulans populations reflects an ancient invasion of D. simulans by tirant copies followed by the loss of active copies from most populations, apart from the African ones, suggesting that this loss is still ongoing in this species.

Identification, characterization and comparative genomics of chimpanzee endogenous retroviruses

BACKGROUND

Retrotransposons, the most abundant and widespread class of eukaryotic transposable elements, are believed to play a significant role in mutation and disease and to have contributed significantly to the evolution of genome structure and function. The recent sequencing of the chimpanzee genome is providing an unprecedented opportunity to study the functional significance of these elements in two closely related primate species and to better evaluate their role in primate evolution.

RESULTS

We report here that the chimpanzee genome contains at least 42 separate families of endogenous retroviruses, nine of which were not previously identified. All but two (CERV 1/PTERV1 and CERV 2) of the 42 families of chimpanzee endogenous retroviruses were found to have orthologs in humans. Molecular analysis (PCR and Southern hybridization) of CERV 2 elements demonstrates that this family is present in chimpanzee, bonobo, gorilla and old-world monkeys but absent in human, orangutan and new-world monkeys. A survey of endogenous retroviral positional variation between chimpanzees and humans determined that approximately 7% of all chimpanzee-human INDEL variation is associated with endogenous retroviral sequences.

CONCLUSION

Nine families of chimpanzee endogenous retroviruses have been transpositionally active since chimpanzees and humans diverged from a common ancestor. Seven of these transpositionally active families have orthologs in humans, one of which has also been transpositionally active in humans since the human-chimpanzee divergence about six million years ago. Comparative analyses of orthologous regions of the human and chimpanzee genomes have revealed that a significant portion of INDEL variation between chimpanzees and humans is attributable to endogenous retroviruses and may be of evolutionary significance.

A genome-wide screening of BEL-Pao like retrotransposons in Anopheles gambiae by the LTR_STRUC program

The advanced status of assembly of the nematoceran Anopheles gambiae genomic sequence allowed us to perform a wide genome analysis to looking at the presence of Long Terminal Repeats (LTRs) in the range of 10 kb by means of the LTR_STRUC tool. More than three hundred sequences were retrieved and 210 were treated as putative complete retrotransposons that were individually analysed with respect to known retrotransposons of A. gambiae and D. melanogaster. The results show that the vast majority of the retrotransposons analysed belong to the Ty3/gypsy class and only 8% to the Ty1/copia class. In addition, phylogenetic analysis allowed us to characterize in more detail the relationship of a large BEL-Pao lineage in which a single family was shown to harbour an additional env gene.

The Sinbad retrotransposon from the genome of the human blood fluke, Schistosoma mansoni, and the distribution of related Pao-like elements

Background

Of the major families of long terminal repeat (LTR) retrotransposons, the Pao/BEL family is probably the least well studied. It is becoming apparent that numerous LTR retrotransposons and other mobile genetic elements have colonized the genome of the human blood fluke, Schistosoma mansoni.

Results

A proviral form of Sinbad, a new LTR retrotransposon, was identified in the genome of S. mansoni. Phylogenetic analysis indicated that Sinbad belongs to one of five discreet subfamilies of Pao/BEL like elements. BLAST searches of whole genomes and EST databases indicated that members of this clade occurred in species of the Insecta, Nematoda, Echinodermata and Chordata, as well as Platyhelminthes, but were absent from all plants, fungi and lower eukaryotes examined. Among the deuterostomes examined, only aquatic species harbored these types of elements. All four species of nematode examined were positive for Sinbad sequences, although among insect and vertebrate genomes, some were positive and some negative. The full length, consensus Sinbad retrotransposon was 6,287 bp long and was flanked at its 5'- and 3'-ends by identical LTRs of 386 bp. Sinbad displayed a triple Cys-His RNA binding motif characteristic of Gag of Pao/BEL-like elements, followed by the enzymatic domains of protease, reverse transcriptase (RT), RNAseH, and integrase, in that order. A phylogenetic tree of deduced RT sequences from 26 elements revealed that Sinbad was most closely related to an unnamed element from the zebrafish Danio rerio and to Saci-1, also from S. mansoni. It was also closely related to Pao from Bombyx mori and to Ninja of Drosophila simulans. Sinbad was only distantly related to the other schistosome LTR retrotransposons Boudicca, Gulliver, Saci-2, Saci-3, and Fugitive, which are gypsy-like. Southern hybridization and bioinformatics analyses indicated that there were about 50 copies of Sinbad in the S. mansoni genome. The presence of ESTs representing transcripts of Sinbad in numerous developmental stages of S. mansoni along with the identical 5'- and 3'-LTR sequences suggests that Sinbad is an active retrotransposon.

Conclusion

Sinbad is a Pao/BEL type retrotransposon from the genome of S. mansoni. The Pao/BEL group appears to be comprised of at least five discrete subfamilies, which tend to cluster with host species phylogeny. Pao/BEL type elements appear to have colonized only the genomes of the Animalia. The distribution of these elements in the Ecdysozoa, Deuterostomia, and Lophotrochozoa is discontinuous, suggesting horizontal transmission and/or efficient elimination of Pao-like mobile genetic elements from some genomes.

New regulatory regions of Drosophila 412 retrotransposable element generated by recombination

There are no doubts that transposable elements (TEs) have greatly influenced genomes evolution. They have, however, evolved in different ways throughout mammals, plants, and invertebrates. In mammals they have been shown to be widely present but with low transposition activity; in plants they are responsible for large increases in genome size. In Drosophila, despite their low amount, transposition seems to be higher. Therefore, to understand how these elements have evolved in different genomes and how host genomes have proposed to go around them, are major questions on genome evolution. We analyzed sequences of the retrotransposable elements 412 in natural populations of the Drosophila simulans and D. melanogaster species that greatly differ in their amount of TEs. We identified new subfamilies of this element that were the result of mutation or insertion-deletion process, but also of interfamily recombinations. These new elements were well conserved in the D. simulans natural populations. The new regulatory regions produced by recombination could give rise to new elements able to overcome host control of transposition and, thus, become potential genome invaders.

The Evolution of Alloimmunity and the Genesis of Adaptive Immunity

Infectious agents select for host immune responses that destroy infectious nonself yet maintain tolerance to self. Here we propose that retroviruses and other host-antigen associated pathogens (HAAPs) select for the genetic, biochemical, and cell biological properties of alloimmunity, also known as the histocompatibility or tissue rejection response. This hypothesis predicts the major observations regarding histocompatibility responses, including: (i) their existence in animals as diverse as sponges and humans; (ii) extreme polymorphism and balanced allele frequencies at histocompatibility loci, including the human MHC and blood group loci; (iii) the frequency dependent selection of histocompatibility alleles; (iv) the ancient age of many alloantigenic polymorphisms; (v) the high ratio of nonsynonymous mutations to synonymous mutations at histocompatibility loci; (vi) disassortative mating based on MHC alleles; (vii) the inability to explain the existence and continuing selection of histocompatibility alleles by other more conventional biochemical and genetic paradigms; and (viii) the susceptibility of HAAPs, particularly retroviruses such as HIV (human immunodeficiency virus), to histocompatibility reactions. In addition, the hypothesis that HAAPs select the forms and molecules of alloimmunity offers simple explanations for the evolution of histocompatibility systems over time, the initial selection of hypervariable immune mechanisms, and the genesis of adaptive immunity.

The diversity of LTR retrotransposons

Although most long terminal repeat retrotransposons have common structural features and encode similar genes, there is nonetheless considerable diversity.

Eukaryotic genomes are full of long terminal repeat (LTR) retrotransposons. Although most LTR retrotransposons have common structural features and encode similar genes, there is nonetheless considerable diversity in their genomic organization, reflecting the different strategies they use to proliferate within the genomes of their hosts.

Evolutionary Genomics of Chromoviruses in Eukaryotes

The diversity, origin, and evolution of chromoviruses in Eukaryota were examined using the massive amount of genome sequence data for different eukaryotic lineages. A surprisingly large number of novel full-length chromoviral elements were found, greatly exceeding the number of the known chromoviruses. These new elements are mostly structurally intact and highly conserved. Chromoviruses in the key Amniota lineage, the reptiles, have been analyzed by PCR to explain their evolutionary dynamics in amniotes. Phylogenetic analyses provide evidence for a novel centromere-specific chromoviral clade that is widespread and highly conserved in all seed plants. Chromoviral diversity in plants, fungi, and vertebrates, as shown by phylogenetic analyses, was found to be much greater than previously expected. The age of plant chromoviruses has been significantly extended by finding their representatives in the most basal plant lineages, the green and the red algae. The evolutionary origin of chromoviruses has been found to be no earlier than in Cercozoa. The evolutionary history and dynamics of chromoviruses can be explained simply by strict vertical transmission in plants, followed by more complex evolution in fungi and in Metazoa. The currently available data clearly show that chromoviruses indeed represent the oldest and the most widespread clade of Metaviridae.

Evolutionary history of Oryza sativa LTR retrotransposons: a preliminary survey of the rice genome sequences

Background

LTR Retrotransposons transpose through reverse transcription of an RNA intermediate and are ubiquitous components of all eukaryotic genomes thus far examined. Plant genomes, in particular, have been found to be comprised of a remarkably high number of LTR retrotransposons. There is a significant body of direct and indirect evidence that LTR retrotransposons have contributed to gene and genome evolution in plants.

Results

To explore the evolutionary history of long terminal repeat (LTR) retrotransposons and their impact on the genome of Oryza sativa, we have extended an earlier computer-based survey to include all identifiable full-length, fragmented and solo LTR elements in the rice genome database as of April 2002. A total of 1,219 retroelement sequences were identified, including 217 full-length elements, 822 fragmented elements, and 180 solo LTRs. In order to gain insight into the chromosomal distribution of LTR-retrotransposons in the rice genome, a detailed examination of LTR-retrotransposon sequences on Chromosome 10 was carried out. An average of 22.3 LTR-retrotransposons per Mb were detected in Chromosome 10.

Conclusions

Gypsy-like elements were found to be >4 × more abundant than copia-like elements. Eleven of the thirty-eight investigated LTR-retrotransposon families displayed significant subfamily structure. We estimate that at least 46.5% of LTR-retrotransposons in the rice genome are older than the age of the species (< 680,000 years). LTR-retrotransposons present in the rice genome range in age from those just recently inserted up to nearly 10 million years old. Approximately 20% of LTR retrotransposon sequences lie within putative genes. The distribution of elements across chromosome 10 is non-random with the highest density (48 elements per Mb) being present in the pericentric region.

Divergent long-terminal-repeat retrotransposon families in the genome of Paragonimus westermani

To gain information on retrotransposons in the genome of Paragonimus westermani, PCR was carried out with degenerate primers, specific to protease and reverse transcriptase (rt) genes of long-terminal-repeat (LTR) retrotransposons. The PCR products were cloned and sequenced, after which 12 different retrotransposon-related sequences were isolated from the trematode genome. These showed various degrees of identity to the polyprotein of divergent retrotransposon families. A phylogenetic analysis demonstrated that these sequences could be classified into three different families of LTR retrotransposons, namely, Xena, Bel, and Gypsy families. Of these, two mRNA transcripts were detected by reverse transcriptase-PCR, showing that these two elements preserved their mobile activities. The genomic distributions of these two sequences were found to be highly repetitive. These results suggest that there are diverse retrotransposons including the ancient Xena family in the genome of P. westermani, which may have been involved in the evolution of the host genome.

Long terminal repeat retrotransposons of Mus musculus

This article describes the identification of long terminal repeat retrotransposons in the mouse genome. These elements can be classified in to 21 families, of which 13 have not been described previously.

Background

Long terminal repeat (LTR) retrotransposons make up a large fraction of the typical mammalian genome. They comprise about 8% of the human genome and approximately 10% of the mouse genome. On account of their abundance, LTR retrotransposons are believed to hold major significance for genome structure and function. Recent advances in genome sequencing of a variety of model organisms has provided an unprecedented opportunity to evaluate better the diversity of LTR retrotransposons resident in eukaryotic genomes.

Results

Using a new data-mining program, LTR_STRUC, in conjunction with conventional techniques, we have mined the GenBank mouse (Mus musculus) database and the more complete Ensembl mouse dataset for LTR retrotransposons. We report here that the M. musculus genome contains at least 21 separate families of LTR retrotransposons; 13 of these families are described here for the first time.

Conclusions

All families of mouse LTR retrotransposons are members of the gypsy-like superfamily of retroviral-like elements. Several different families of unrelated non-autonomous elements were identified, suggesting that the evolution of non-autonomy may be a common event. High sequence similarity between several LTR retrotransposons identified in this study and those found in distantly-related species suggests that horizontal transfer has been a significant factor in the evolution of mouse LTR retrotransposons.

MAX, a novel retrotransposon of the BEL-Pao family, is nested within the Bari1 cluster at the heterochromatic h39 region of chromosome 2 in Drosophila melanogaster

A homogeneous array of 80 tandem repeats of the Bari1 transposon is located in the pericentromeric h39 region of chromosome 2 of Drosophila melanogaster. Here, we report that the Bari1 cluster is interrupted by an 8556-bp insertion. DNA sequencing and database searches identified this insertion as a previously unannotated retrotransposon that we have named MAX. MAX possesses two ORFs; ORF1 putatively encodes a polyprotein comprising GAG and RT domains, while ORF2 could encode a 288-amino acid protein of unknown function. Alignment with the RT domains of known LTR retrotransposons shows that MAX belongs to the BEL-Pao family, which remarkable for its widespread presence in different taxa, including lower chordates. We have analyzed the distribution of MAX elements within representative species of the Sophophora subgroup and found that they are restricted to the species of the melanogaster complex, where they are heavily represented in the heterochromatin of all autosomes and on the Y chromosome.

Retrotransposons and their recognition of pol II promoters: a comprehensive survey of the transposable elements from the complete genome sequence of Schizosaccharomyces pombe

The complete DNA sequence of the genome of Schizosaccharomyces pombe provides the opportunity to investigate the entire complement of transposable elements (TEs), their association with specific sequences, their chromosomal distribution, and their evolution. Using homology-based sequence identification, we found that the sequenced strain of S. pombe contained only one family of full-length transposons. This family, Tf2, consisted of 13 full-length copies of a long terminal repeat (LTR) retrotransposon. We found that LTR-LTR recombination of previously existing transposons had resulted in extensive populations of solo LTRs. These included 35 solo LTRs of Tf2, as well as 139 solo LTRs from other Tf families. Phylogenetic analysis of solo Tf LTRs reveals that Tf1 and Tf2 were the most recently active elements within the genome. The solo LTRs also served as footprints for previous insertion events by the Tf retrotransposons. Analysis of 186 genomic insertion events revealed a close association with RNA polymerase II promoters. These insertions clustered in the promoter-proximal regions of genes, upstream of protein coding regions by 100 to 400 nucleotides. The association of Tf insertions with pol II promoters was very similar to the preference previously observed for Tf1 integration. We found that the recently active Tf elements were absent from centromeres and pericentromeric regions of the genome containing tandem tRNA gene clusters. In addition, our analysis revealed that chromosome III has twice the density of insertion events compared to the other two chromosomes. Finally we describe a novel repetitive sequence, wtf, which was also preferentially located on chromosome III, and was often located near solo LTRs of Tf elements.

Characterization of high-copy-number retrotransposons from the large genomes of the louisiana iris species and their use as molecular markers

The Louisiana iris species Iris brevicaulis and I. fulva are morphologically and karyotypically distinct yet frequently hybridize in nature. A group of high-copy-number TY3/gypsy-like retrotransposons was characterized from these species and used to develop molecular markers that take advantage of the abundance and distribution of these elements in the large iris genome. The copy number of these IRRE elements (for iris retroelement), is approximately 1 x 10(5), accounting for approximately 6-10% of the approximately 10,000-Mb haploid Louisiana iris genome. IRRE elements are transcriptionally active in I. brevicaulis and I. fulva and their F(1) and backcross hybrids. The LTRs of the elements are more variable than the coding domains and can be used to define several distinct IRRE subfamilies. Transposon display or S-SAP markers specific to two of these subfamilies have been developed and are highly polymorphic among wild-collected individuals of each species. As IRRE elements are present in each of 11 iris species tested, the marker system has the potential to provide valuable comparative data on the dynamics of retrotransposition in large plant genomes.

Boudicca, a retrovirus-like long terminal repeat retrotransposon from the genome of the human blood fluke Schistosoma mansoni

The genome of Schistosoma mansoni contains a proviral form of a retrovirus-like long terminal repeat (LTR) retrotransposon, designated BOUDICCA: Sequence and structural characterization of the new mobile genetic element, which was found in bacterial artificial chromosomes prepared from S. mansoni genomic DNA, revealed the presence of three putative open reading frames (ORFs) bounded by direct LTRs of 328 bp in length. ORF1 encoded a retrovirus-like major homology region and a Cys/His box motif, also present in Gag polyproteins of related retrotransposons and retroviruses. ORF2 encoded enzymatic domains and motifs characteristic of a retrovirus-like polyprotein, including aspartic protease, reverse transcriptase, RNase H, and integrase, in that order, a domain order similar to that of the gypsy/Ty3 retrotransposons. An additional ORF at the 3' end of the retrotransposon may encode an envelope protein. Phylogenetic comparison based on the reverse transcriptase domain of ORF2 confirmed that Boudicca was a gypsy-like retrotransposon and showed that it was most closely related to CsRn1 from the Oriental liver fluke Clonorchis sinensis and to kabuki from Bombyx mori. Bioinformatics approaches together with Southern hybridization analysis of genomic DNA of S. mansoni and the screening of a bacterial artificial chromosome library representing approximately 8-fold coverage of the S. mansoni genome revealed that numerous copies of Boudicca were interspersed throughout the schistosome genome. By reverse transcription-PCR, mRNA transcripts were detected in the sporocyst, cercaria, and adult developmental stages of S. mansoni, indicating that Boudicca is actively transcribed in this trematode.

Multiple ribonuclease H-encoding genes in the Caenorhabditis elegans genome contrasts with the two typical ribonuclease H-encoding genes in the human genome

Database searches of the Caenorhabditis elegans and human genomic DNA sequences revealed genes encoding ribonuclease H1 (RNase H1) and RNase H2 in each genome. The human genome contains a single copy of each gene, whereas C. elegans has four genes encoding RNase H1-related proteins and one gene for RNase H2. By analyzing the mRNAs produced from the C. elegans genes, examining the amino acid sequence of the predicted protein, and expressing the proteins in Esherichia coli we have identified two active RNase H1-like proteins. One is similar to other eukaryotic RNases H1, whereas the second RNase H (rnh-1.1) is unique. The rnh-1.0 gene is transcribed as a dicistronic message with three dsRNA-binding domains; the mature mRNA is transspliced with SL2 splice leader and contains only one dsRNA-binding domain. Formation of RNase H1 is further regulated by differential cis-splicing events. A single rnh-2 gene, encoding a protein similar to several other eukaryotic RNase H2L's, also has been examined. The diversity and enzymatic properties of RNase H homologues are other examples of expansion of protein families in C. elegans. The presence of two RNases H1 in C. elegans suggests that two enzymes are required in this rather simple organism to perform the functions that are accomplished by a single enzyme in more complex organisms. Phylogenetic analysis indicates that the active C. elegans RNases H1 are distantly related to one another and that the C. elegans RNase H1 is more closely related to the human RNase H1. The database searches also suggest that RNase H domains of LTR-retrotransposons in C. elegans are quite unrelated to cellular RNases H1, but numerous RNase H domains of human endogenous retroviruses are more closely related to cellular RNases H.

Genes of the Pseudoviridae (Ty1/copia retrotransposons)

A comprehensive survey of the Pseudoviridae (Ty1/copia) retroelement family was conducted using the GenBank sequence database and completed genome sequences of several model organisms. Plant genomes were the most abundant sources of Pseudoviridae, with the Arabidopsis thaliana genome having 276 distinct elements. A reverse transcriptase amino acid sequence phylogeny indicated that the Pseudoviridae comprises highly divergent members. Coding sequences for a representative subset of elements were analyzed to identify conserved domains and differences that may underlie functional divergence. With the exception of some fungal elements (e.g., Ty1), most Pseudoviridae encode Gag and Pol on a single open reading frame. In addition to the nearly ubiquitous RNA-binding motif of nucleocapsid, three new conserved domains were identified in Gag. pol-encoded aspartic protease was similar to the retroviral enzyme and could be mapped onto the HIV-1 structure. Pol was highly conserved throughout the family. The greatest divergence among Pol sequences was seen in the C-terminus of integrase (IN). We defined a large motif (GKGY) after the IN catalytic domain that is unique to the Pseudoviridae. Additionally, the extreme C-terminus of IN is rich in simple sequence motifs. A distinct lineage of Pseudoviridae in plants have envlike genes. This lineage has undergone a large expansion of Gag characterized by an alpha-helix-rich domain containing coiled-coil motifs. In several elements, this domain is flanked on both sides by RNA-binding domains. We propose that this monophyletic lineage defines a new Pseudoviridae genus, herein referred to as the AGROVIRUS:

Sequence variation and genomic amplification of a family of Gypsy-like elements in the oomycete genus Phytophthora

A family of sequences resembling Gypsy retroelements was identified and shown to be widely distributed throughout the genus Phytophthora, a member of the algallike oomycete fungi. Polymerase chain reaction (PCR) using specific and degenerate primers detected the family in 29 of 37 species tested. DNA hybridization also failed to detect the sequences in the eight species that were negative in PCR. The element appears to have been a major force in the shaping of Phytophthora genomes because its abundance varied drastically from about 10 to more than 10,000 copies per genome within the species containing the element. Family members diverged from each other by single-base changes, insertions, and deletions, with a mean nucleotide divergence of 16.7%. By constructing phylogenies of the elements, lineages were identified that predated speciation events within Phytophthora and subfamilies that had diverged more recently. The element was studied in detail in Phytophthora infestans, in which about 30 copies are dispersed throughout the genome. Phylogenetic comparisons of the reverse transcriptases placed the family within the Ty3/Gypsy group of long terminal repeat (LTR) retrotransposons, with the closest affinities to elements from plants. However, each of 12 family members sequenced contained defects that would render their protein products inactive, including frameshift mutations within reverse transcriptase domains and truncations that appeared to eliminate gag, protease, and terminal repeat sequences.

Intragenomic conflict and cancer

Intragenomic conflict occurs when some elements within the genome produce effects that enhance their own probability of replication or transmission at the expense of other elements within the same genome. Here it is proposed that mutations involved in intragenomic conflict are particularly likely to be co-opted by evolving lineages of cancer cells, and hence should be associated with the occurrence of cancer. We discuss several types of intragenomic conflict that are associated with various forms of cancer.