Evolutionary history of Cer elements and their impact on the C. elegans genome

From parasites to partners: exploring the intricacies of host-transposon dynamics and coevolution

Transposable elements, often referred to as "jumping genes," have long been recognized as genomic parasites due to their ability to integrate and disrupt normal gene function and induce extensive genomic alterations, thereby compromising the host's fitness. To counteract this, the host has evolved a plethora of mechanisms to suppress the activity of the transposons. Recent research has unveiled the host-transposon relationships to be nuanced and complex phenomena, resulting in the coevolution of both entities. Transposition increases the mutational rate in the host genome, often triggering physiological pathways such as immune and stress responses. Current gene transfer technologies utilizing transposable elements have potential drawbacks, including off-target integration, induction of mutations, and modifications of cellular machinery, which makes an in-depth understanding of the host-transposon relationship imperative. This review highlights the dynamic interplay between the host and transposable elements, encompassing various factors and components of the cellular machinery. We provide a comprehensive discussion of the strategies employed by transposable elements for their propagation, as well as the mechanisms utilized by the host to mitigate their parasitic effects. Additionally, we present an overview of recent research identifying host proteins that act as facilitators or inhibitors of transposition. We further discuss the evolutionary outcomes resulting from the genetic interactions between the host and the transposable elements. Finally, we pose open questions in this field and suggest potential avenues for future research.

Virus-like transposons cross the species barrier and drive the evolution of genetic incompatibilities

Horizontal gene transfer, the movement of genetic material between species, has been reported across all major eukaryotic lineages. However, the underlying mechanisms of transfer and their impact on genome evolution are still poorly understood. While studying the evolutionary origin of a selfish element in the nematode Caenorhabditis briggsae, we discovered that Mavericks, ancient virus-like transposons related to giant viruses and virophages, are one of the long-sought vectors of horizontal gene transfer. We found that Mavericks gained a novel herpesvirus-like fusogen in nematodes, leading to the widespread exchange of cargo genes between extremely divergent species, bypassing sexual and genetic barriers spanning hundreds of millions of years. Our results show how the union between viruses and transposons causes horizontal gene transfer and ultimately genetic incompatibilities in natural populations.

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.

Long-term experimental evolution reveals purifying selection on piRNA-mediated control of transposable element expression

Background

Transposable elements (TEs) are an almost universal constituent of eukaryotic genomes. In animals, Piwi-interacting small RNAs (piRNAs) and repressive chromatin often play crucial roles in preventing TE transcription and thus restricting TE activity. Nevertheless, TE content varies widely across eukaryotes and the dynamics of TE activity and TE silencing across evolutionary time is poorly understood.

Results

Here, we used experimentally evolved populations of C. elegans to study the dynamics of TE expression over 409 generations. The experimental populations were evolved at population sizes of 1, 10 and 100 individuals to manipulate the efficiency of natural selection versus genetic drift. We demonstrate increased TE expression relative to the ancestral population, with the largest increases occurring in the smallest populations. We show that the transcriptional activation of TEs within active regions of the genome is associated with failure of piRNA-mediated silencing, whilst desilenced TEs in repressed chromatin domains retain small RNAs. Additionally, we find that the sequence context of the surrounding region influences the propensity of TEs to lose silencing through failure of small RNA-mediated silencing.

Conclusions

Our results show that natural selection in C. elegans is responsible for maintaining low levels of TE expression, and provide new insights into the epigenomic features responsible.

Human Endogenous Retrovirus K (HERV-K) can drive gene expression as a promoter in Caenorhabditis elegans

Endogenous retroviruses (ERVs) are retrotransposons present in various metazoan genomes and have been implicated in metazoan evolution as well as in nematodes and humans. The long terminal repeat (LTR) retrotransposons contain several regulatory sequences including promoters and enhancers that regulate endogenous gene expression and thereby control organismal development and response to environmental change. ERVs including the LTR retrotransposons constitute 8% of the human genome and less than 0.6% of the Caenorhabditis elegans (C.elegans) genome, a nematode genetic model system. To investigate the evolutionarily conserved mechanism behind the transcriptional activity of retrotransposons, we generated a transgenic worm model driving green fluorescent protein (GFP) expression using Human endogenous retroviruses (HERV)-K LTR as a promoter. The promoter activity of HERV-K LTR was robust and fluorescence was observed in various tissues throughout the developmental process. Interestingly, persistent GFP expression was specifically detected in the adult vulva muscle. Using deletion constructs, we found that the region from positions 675 to 868 containing the TATA box was necessary for promoter activity driving gene expression in the vulva. Interestingly, we found that the promoter activity of the LTR was dependent on che-1 transcription factor, a sensory neuron driver, and lin-15b, a negative regulator of RNAi and germline gene expression. These results suggest evolutionary conservation of the LTR retrotransposon activity in transcriptional regulation as well as the possibility of che-1 function in non-neuronal tissues.

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 Importance of Satellite Sequence Repression for Genome Stability

Up to two-thirds of eukaryotic genomes consist of repetitive sequences, which include both transposable elements and tandemly arranged simple or satellite repeats. Whereas extensive progress has been made toward understanding the danger of and control over transposon expression, only recently has it been recognized that DNA damage can arise from satellite sequence transcription. Although the structural role of satellite repeats in kinetochore function and end protection has long been appreciated, it has now become clear that it is not only these functions that are compromised by elevated levels of transcription. RNA from simple repeat sequences can compromise replication fork stability and genome integrity, thus compromising germline viability. Here we summarize recent discoveries on how cells control the transcription of repeat sequence and the dangers that arise from their expression. We propose that the link between the DNA damage response and the transcriptional silencing machinery may help a cell or organism recognize foreign DNA insertions into an evolving genome.

Alternative Oxidase Gene Family in Hypericum perforatum L.: Characterization and Expression at the Post-germinative Phase

Alternative oxidase (AOX) protein is located in the inner mitochondrial membrane and is encoded in the nuclear genome being involved in plant response upon a diversity of environmental stresses and also in normal plant growth and development. Here we report the characterization of the AOX gene family of Hypericum perforatum L. Two AOX genes were identified, both with a structure of four exons (HpAOX1, acc. KU674355 and HpAOX2, acc. KU674356). High variability was found at the N-terminal region of the protein coincident with the high variability identified at the mitochondrial transit peptide. In silico analysis of regulatory elements located at intronic regions identified putative sequences coding for miRNA precursors and trace elements of a transposon. Simple sequence repeats were also identified. Additionally, the mRNA levels for the HpAOX1 and HpAOX2, along with the ones for the HpGAPA (glyceraldehyde-3-phosphate dehydrogenase A subunit) and the HpCAT1 (catalase 1), were evaluated during the post-germinative development. Gene expression analysis was performed by RT-qPCR with accurate data normalization, pointing out HpHYP1 (chamba phenolic oxidative coupling protein 1) and HpH2A (histone 2A) as the most suitable reference genes (RGs) according to GeNorm algorithm. The HpAOX2 transcript demonstrated larger stability during the process with a slight down-regulation in its expression. Contrarily, HpAOX1 and HpGAPA (the corresponding protein is homolog to the chloroplast isoform involved in the photosynthetic carbon assimilation in other plant species) transcripts showed a marked increase, with a similar expression pattern between them, during the post-germinative development. On the other hand, the HpCAT1 (the corresponding protein is homolog to the major H2O2-scavenging enzyme in other plant species) transcripts showed an opposite behavior with a down-regulation during the process. In summary, our findings, although preliminary, highlight the importance to investigate in more detail the participation of AOX genes during the post-germinative development in H. perforatum, in order to explore their functional role in optimizing photosynthesis and in the control of reactive oxygen species (ROS) levels during the process.

Analysis of the repetitive component and retrotransposon population in the genome of a marine angiosperm, Posidonia oceanica (L.) Delile

Posidonia oceanica is a monocotyledonous marine plant that plays a crucial role in maintaining the Mediterranean environment. Despite its ecological importance, basic knowledge of the functional and structural genomics of this species is still limited, as it is for the other seagrasses. Here, for the first time, we report data on the repetitive component of the genome of this seagrass using a low coverage of Illumina sequences and different assembly approaches. A dataset of 19,760 assembled sequences, mostly belonging to the repetitive fraction of the genome, was produced and annotated. Based on mapping Illumina reads onto this dataset, the genome structure of P. oceanica and its repetitive component was inferred. A very large proportion of the genome is represented by long-terminal-repeat (LTR) retrotransposons of both the Copia and Gypsy superfamilies. Posidonia LTR-retrotransposons were classified and their sequences analysed. Gypsy elements belong to three main lineages, while Copia ones belong to seven lineages. Gypsy elements were more represented than Copia ones in the set of assembled sequences and in the genome. Analysis of sequence variability indicated that Gypsy lineages have experienced amplification in more recent times compared to Copia ones.

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.

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.

Mosquitoes LTR retrotransposons: a deeper view into the genomic sequence of Culex quinquefasciatus

A set of 67 novel LTR-retrotransposon has been identified by in silico analyses of the Culex quinquefasciatus genome using the LTR_STRUC program. The phylogenetic analysis shows that 29 novel and putatively functional LTR-retrotransposons detected belong to the Ty3/gypsy group. Our results demonstrate that, by considering only families containing potentially autonomous LTR-retrotransposons, they account for about 1% of the genome of C. quinquefasciatus. In previous studies it has been estimated that 29% of the genome of C. quinquefasciatus is occupied by mobile genetic elements.The potential role of retrotransposon insertions strictly associated with host genes is described and discussed along with the possible origin of a retrotransposon with peculiar Primer Binding Site region. Finally, we report the presence of a group of 38 retrotransposons, carrying tandem repeated sequences but lacking coding potential, and apparently lacking "master copy" elements from which they could have originated. The features of the repetitive sequences found in these non-autonomous LTR retrotransposons are described, and their possible role discussed.These results integrate the existing data on the genomics of an important virus-borne disease vector.

BEL/Pao retrotransposons in metazoan genomes

BACKGROUND

Long terminal repeat (LTR) retrotransposons are a widespread kind of transposable element present in eukaryotic genomes. They are a major factor in genome evolution due to their ability to create large scale mutations and genome rearrangements. Compared to other transposable elements, little attention has been paid to elements belonging to the metazoan BEL/Pao subclass of LTR retrotransposons. No comprehensive characterization of these elements is available so far. The aim of this study was to describe all BEL/Pao elements in a set of 62 sequenced metazoan genomes, and to analyze their phylogenetic relationship.

RESULTS

We identified a total of 7,861 BEL/Pao elements in 53 of our 62 study genomes. We identified BEL/Pao elements in 20 genomes where such elements had not been found so far. Our analysis shows that BEL/Pao elements are the second-most abundant class of LTR retrotransposons in the genomes we study, more abundant than Ty1/Copia elements, and second only to Ty3/Gypsy elements. They occur in multiple phyla, including basal metazoan phyla, suggesting that BEL/Pao elements arose early in animal evolution. We confirm findings from previous studies that BEL/Pao elements do not occur in mammals. The elements we found can be grouped into more than 1725 families, 1623 of which are new, previously unknown families. These families fall into seven superfamilies, only five of which have been characterized so far. One new superfamily is a major subdivision of the Pao superfamily which we propose to call Dan, because it is restricted to the genome of the zebrafish Danio rerio. The other new superfamily comprises 83 elements and is restricted to lower aquatic eumetazoans. We propose to call this superfamily Flow. BEL/Pao elements do not show any signs of recent horizontal gene transfer between distantly related species.

CONCLUSIONS

In sum, our analysis identifies thousands of new BEL/Pao elements and provides new insights into their distribution, abundance, and evolution.

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.

Evolution of the Caenorhabditis elegans genome

A fundamental problem in genome biology is to elucidate the evolutionary forces responsible for generating nonrandom patterns of genome organization. As the first metazoan to benefit from full-genome sequencing, Caenorhabditis elegans has been at the forefront of research in this area. Studies of genomic patterns, and their evolutionary underpinnings, continue to be augmented by the recent push to obtain additional full-genome sequences of related Caenorhabditis taxa. In the near future, we expect to see major advances with the onset of whole-genome resequencing of multiple wild individuals of the same species. In this review, we synthesize many of the important insights to date in our understanding of genome organization and function that derive from the evolutionary principles made explicit by theoretical population genetics and molecular evolution and highlight fertile areas for future research on unanswered questions in C. elegans genome evolution. We call attention to the need for C. elegans researchers to generate and critically assess nonadaptive hypotheses for genomic and developmental patterns, in addition to adaptive scenarios. We also emphasize the potential importance of evolution in the gonochoristic (female and male) ancestors of the androdioecious (hermaphrodite and male) C. elegans as the source for many of its genomic and developmental patterns.

LTR retrotransposons reveal recent extensive inter-subspecies nonreciprocal recombination in Asian cultivated rice

Background

Long Terminal Repeats retrotransposons (LTR elements) are ubiquitous Eukaryotic transposable elements (TEs). They are considered to be one of the major forces underlying plant genome evolution. Because of relatively high evolutionary speed, active transposition of LTR elements in the host genomes provides rich information on their short-term history. As more and more genomes, especially those of closely related organisms, have been sequenced, it is possible to perform global comparative study of their LTR retrotransposons to reveal events in the history.

Results

The present research is designed to investigate important evolutionary events in the origin of Asian cultivated rice through the comparison of LTR elements. We have developed LTR_INSERT, a new method for LTR elements discovery in two closely related genomes. Our method has a distinctive feature that it is capable of judging whether an insertion occurs prior or posterior to the divergence of genomes. LTR_INSERT identifies 993 full-length LTR elements, annotates 15916 copies related with them, and discovers at least 16 novel LTR families in the whole-genome comparative map of two cultivated rice subspecies. From the full-length LTR elements, we estimate that a significant proportion of the rice genome has experienced inter-subspecies nonreciprocal recombination (ISNR) in as recent as 53,000 years. Large-scale samplings further support that more than 15% of the rice genome has been involved in such recombination. In addition, LTR elements confirm that the genome of O. sativa ssp. indica and that of japonica diverged about 600,000 years ago.

Conclusion

A new LTR retrotransposon identification method integrating both comparative genomics and ab initio algorithm is introduced and applied to Asian cultivated rice genomes. At whole-genome level, this work confirms that recent ISNR is an important factor that molds modern cultivated rice genome.

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.

Molecular basis of the copulatory plug polymorphism in Caenorhabditis elegans

Heritable variation is the raw material for evolutionary change, and understanding its genetic basis is one of the central problems in modern biology. We investigated the genetic basis of a classic phenotypic dimorphism in the nematode Caenorhabditis elegans. Males from many natural isolates deposit a copulatory plug after mating, whereas males from other natural isolates?including the standard wild-type strain (N2 Bristol) that is used in most research laboratories?do not deposit plugs. The copulatory plug is a gelatinous mass that covers the hermaphrodite vulva, and its deposition decreases the mating success of subsequent males. We show that the plugging polymorphism results from the insertion of a retrotransposon into an exon of a novel mucin-like gene, plg-1, whose product is a major structural component of the copulatory plug. The gene is expressed in a subset of secretory cells of the male somatic gonad, and its loss has no evident effects beyond the loss of male mate-guarding. Although C. elegans descends from an obligate-outcrossing, male?female ancestor, it occurs primarily as self-fertilizing hermaphrodites. The reduced selection on male?male competition associated with the origin of hermaphroditism may have permitted the global spread of a loss-of-function mutation with restricted pleiotropy.

LTR retrotransposon landscape in Medicago truncatula: more rapid removal than in rice

Background

Long terminal repeat retrotransposons (LTR elements) are ubiquitous Eukaryotic TEs that transpose through RNA intermediates. Accounting for significant proportion of many plant genomes, LTR elements have been well established as one of the major forces underlying the evolution of plant genome size, structure and function. The accessibility of more than 40% of genomic sequences of the model legume Medicago truncatula (Mt) has made the comprehensive study of its LTR elements possible.

Results

We use a newly developed tool LTR_FINDER to identify LTR retrotransposons in the Mt genome and detect 526 full-length elements as well as a great number of copies related to them. These elements constitute about 9.6% of currently available genomic sequences. They are classified into 85 families of which 64 are reported for the first time. The majority of the LTR retrotransposons belong to either Copia or Gypsy superfamily and the others are categorized as TRIMs or LARDs by their length. We find that the copy-number of Copia-like families is 3 times more than that of Gypsy-like ones but the latter contribute more to the genome. The analysis of PBS and protein-coding domain structure of the LTR families reveals that they tend to use only 4–5 types of tRNAs and many families have quite conservative ORFs besides known TE domains. For several important families, we describe in detail their abundance, conservation, insertion time and structure. We investigate the amplification-deletion pattern of the elements and find that the detectable full-length elements are relatively young and most of them were inserted within the last 0.52 MY. We also estimate that more than ten million bp of the Mt genomic sequences have been removed by the deletion of LTR elements and the removal of the full-length structures in Mt has been more rapid than in rice.

Conclusion

This report is the first comprehensive description and analysis of LTR retrotransposons in the Mt genome. Many important novel LTR families were discovered and their evolution is elucidated. Our results may outline the LTR retrotransposon landscape of the model legume.

LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons

Long terminal repeat retrotransposons (LTR elements) are ubiquitous eukaryotic transposable elements. They play important roles in the evolution of genes and genomes. Ever-growing amount of genomic sequences of many organisms present a great challenge to fast identifying them. That is the first and indispensable step to study their structure, distribution, functions and other biological impacts. However, until today, tools for efficient LTR retrotransposon discovery are very limited. Thus, we developed LTR_FINDER web server. Given DNA sequences, it predicts locations and structure of full-length LTR retrotransposons accurately by considering common structural features. LTR_FINDER is a system capable of scanning large-scale sequences rapidly and the first web server for ab initio LTR retrotransposon finding. We illustrate its usage and performance on the genome of Saccharomyces cerevisiae. The web server is freely accessible at http://tlife.fudan.edu.cn/ltr_finder/.

De novo identification of LTR retrotransposons in eukaryotic genomes

Background

LTR retrotransposons are a class of mobile genetic elements containing two similar long terminal repeats (LTRs). Currently, LTR retrotransposons are annotated in eukaryotic genomes mainly through the conventional homology searching approach. Hence, it is limited to annotating known elements.

Results

In this paper, we report a de novo computational method that can identify new LTR retrotransposons without relying on a library of known elements. Specifically, our method identifies intact LTR retrotransposons by using an approximate string matching technique and protein domain analysis. In addition, it identifies partially deleted or solo LTRs using profile Hidden Markov Models (pHMMs). As a result, this method can de novo identify all types of LTR retrotransposons. We tested this method on the two pairs of eukaryotic genomes, C. elegans vs. C. briggsae and D. melanogaster vs. D. pseudoobscura. LTR retrotransposons in C. elegans and D. melanogaster have been intensively studied using conventional annotation methods. Comparing with previous work, we identified new intact LTR retroelements and new putative families, which may imply that there may still be new retroelements that are left to be discovered even in well-studied organisms. To assess the sensitivity and accuracy of our method, we compared our results with a previously published method, LTR_STRUC, which predominantly identifies full-length LTR retrotransposons. In summary, both methods identified comparable number of intact LTR retroelements. But our method can identify nearly all known elements in C. elegans, while LTR_STRUCT missed about 1/3 of them. Our method also identified more known LTR retroelements than LTR_STRUCT in the D. melanogaster genome. We also identified some LTR retroelements in the other two genomes, C. briggsae and D. pseudoobscura, which have not been completely finished. In contrast, the conventional method failed to identify those elements. Finally, the phylogenetic and chromosomal distributions of the identified elements are discussed.

Conclusion

We report a novel method for de novo identification of LTR retrotransposons in eukaryotic genomes with favorable performance over the existing methods.

A novel non-coding DNA family in Caenorhabditis elegans

Many repetitive elements, for example, SINEs, LINEs, LTR-retrotransposons and other SSRs are dispersed throughout eukaryotic genomes. To understand the biological function of these repetitive elements is of great current research interest. In this study, we report on the identification of a novel non-coding DNA family, designated CE1 family, in the nematode C. elegans genome. Some CE1 elements constituted a large palindrome sequence. The CE1 elements were interspersed at 95 sites in the C. elegans genome. Most of the CE1 elements were associated with, or were within, protein-coding genes. The sequence of the CE1 elements indicated that some could form a hairpin structure. One of the CE1 family, CE1(bs258), is located in the first intron of a novel gene, C46H11.6 which encodes a PDZ/DHR/GLGF domain protein. In gfp and lacZ reporter gene assays the CE1(bs258) element appeared to behave as an enhancer element for the expression of C46H11.6 but no effect on the expression of the opposite direction gene, pat-10 which encodes the body-wall muscle troponin C. The CE1(bs258) RNA transcript was detected by RT-PCR even when CE1(bs258) was located in an intron. We conclude that CE1 elements are involved in the expression of adjacent genes and are therefore selectively retained in the C. elegans genome. We discussed a biological function of the CE1(bs258) having many transcription factor-binding sites.

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.

Abundance, distribution and dynamics of retrotransposable elements and transposons: similarities and differences

Retrotransposable elements and transposons are generally both found in most eukaryotes. These two classes of elements are usually distinguished on the basis of their differing mechanisms of transposition. However, their respective frequencies, their intragenomic dynamics and distributions, and the frequencies of their horizontal transfer from one species to another can also differ. The main objective of this review is to compare these two types of elements from a new perspective, using data provided by genome sequencing projects and relating this to the theoretical and observed dynamics. It is shown that the traditional division into two classes, based on the transposition mechanisms, becomes less obvious when other factors are taken into consideration. A great diversity in distribution and dynamics within each class is observed. In contrast, the impact on and the interactions with the genome can show striking similarities between families of the two classes.

Diversity and clustered distribution of retrotransposable elements in the compact genome of the pufferfish Tetraodon nigroviridis

We report the characterization and chromosomal distribution of retroelements in the compact genome of the pufferfish Tetraodon nigroviridis. We have reconstructed partial/complete retroelement sequences, established their phylogenetic relationship to other known eukaryotic retrotransposons, and performed double-color FISH analyses to gain new insights into their patterns of chromosomal distribution. We could identify 43 different reverse transcriptase retrotransposons belonging to the three major known subclasses (14 non-LTR retrotransposons from seven clades, 25 LTR retrotransposons representing the five major known groups, and four Penelope-like elements), and well as two SINEs (non-autonomous retroelements). Such a diversity of retrotransposable elements, which seems to be relatively common in fish but not in mammals, is astonishing in such a compact genome. The total number of retroelements was approximately 3000, roughly representing only 2.6% of the genome of T. nigroviridis. This is much less than in other vertebrate genomes, reflecting the compact nature of the genome of this pufferfish. Major differences in copy number were observed between different clades, indicating differential success in invading and persisting in the genome. Some retroelements displayed evidence of recent activity. Finally, FISH analysis showed that retrotransposable elements preferentially accumulate in specific heterochromatic regions of the genome of T. nigroviridis, revealing a degree of genomic compartmentalization not observed in the human genome.

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.

Retrotransposon-gene associations are widespread among D. melanogaster populations

We have surveyed 18 natural populations of Drosophila melanogaster for the presence of 23 retrotransposon-gene-association alleles (i.e., the presence of an LTR retrotransposon sequence in or within 1,000 bp of a gene) recently identified in the sequenced D. melanogaster genome. The identified associations were detected only in the D. melanogaster populations. The majority (61%) of the identified retrotransposon-gene associations were present only in the sequenced strain in which they were first identified. Thirty percent of the associations were detected in at least one of the natural populations, and 9% of the associations were detected in all of the D. melanogaster populations surveyed. Sequence analysis of an association allele present in all populations indicates that selection is a significant factor in the spread and/or maintenance of at least some of retroelement-gene associations in D. melanogaster.

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.

Translational recoding signals between gag and pol in diverse LTR retrotransposons

Because of their compact genomes, retroelements (including retrotransposons and retroviruses) employ a variety of translational recoding mechanisms to express Gag and Pol. To assess the diversity of recoding strategies, we surveyed gag/pol gene organization among retroelements from diverse host species, including elements exhaustively recovered from the genome sequences of Caenorhabditis elegans, Drosophila melanogaster, Schizosaccharomyces pombe, Candida albicans, and Arabidopsis thaliana. In contrast to the retroviruses, which typically encode pol in the -1 frame relative to gag, nearly half of the retroelements surveyed encode a single gag-pol open reading frame. This was particularly true for the Ty1/copia group retroelements. Most animal Ty3/gypsy retroelements, on the other hand, encode gag and pol in separate reading frames, and likely express Pol through +1 or -1 frameshifting. Conserved sequences conforming to slippery sites that specify viral ribosomal frameshifting were identified among retroelements with pol in the -1 frame. None of the plant retroelements encoded pol in the -1 frame relative to gag; however, two closely related plant Ty3/gypsy elements encode pol in the +1 frame. Interestingly, a group of plant Ty1/copia retroelements encode pol either in a +1 frame relative to gag or in two nonoverlapping reading frames. These retroelements have a conserved stem-loop at the end of gag, and likely express pol either by a novel means of internal ribosomal entry or by a bypass mechanism.

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.

Continuous exchange of sequence information between dispersed Tc1 transposons in the Caenorhabditis elegans genome

In a genome-wide analysis of the active transposons in Caenorhabditis elegans we determined the localization and sequence of all copies of each of the six active transposon families. Most copies of the most active transposons, Tc1 and Tc3, are intact but individually have a unique sequence, because of unique patterns of single-nucleotide polymorphisms. The sequence of each of the 32 Tc1 elements is invariant in the C. elegans strain N2, which has no germline transposition. However, at the same 32 Tc1 loci in strains with germline transposition, Tc1 elements can acquire the sequence of Tc1 elements elsewhere in the N2 genome or a chimeric sequence derived from two dispersed Tc1 elements. We hypothesize that during double-strand-break repair after Tc1 excision, the template for repair can switch from the Tc1 element on the sister chromatid or homologous chromosome to a Tc1 copy elsewhere in the genome. Thus, the population of active transposable elements in C. elegans is highly dynamic because of a continuous exchange of sequence information between individual copies, potentially allowing a higher evolution rate than that found in endogenous genes.

Evidence for the adaptive significance of an LTR retrotransposon sequence in a Drosophila heterochromatic gene

BACKGROUND

The potential adaptive significance of transposable elements (TEs) to the host genomes in which they reside is a topic that has been hotly debated by molecular evolutionists for more than two decades. Recent genomic analyses have demonstrated that TE fragments are associated with functional genes in plants and animals. These findings suggest that TEs may contribute significantly to gene evolution.

RESULTS

We have analyzed two transposable elements associated with genes in the sequenced Drosophila melanogaster y; cn bw sp strain. A fragment of the Antonia long terminal repeat (LTR) retrotransposon is present in the intron of Chitinase 3 (Cht3), a gene located within the constitutive heterochromatin of chromosome 2L. Within the euchromatin of chromosome 2R a full-length Burdock LTR retrotransposon is located immediately 3' to cathD, a gene encoding cathepsin D. We tested for the presence of these two TE/gene associations in strains representing 12 geographically diverse populations of D. melanogaster. While the cathD insertion variant was detected only in the sequenced y; cn bw sp strain, the insertion variant present in the heterochromatic Cht3 gene was found to be fixed throughout twelve D. melanogaster populations and in a D. mauritiana strain suggesting that it maybe of adaptive significance. To further test this hypothesis, we sequenced a 685bp region spanning the LTR fragment in the intron of Cht3 in strains representative of the two sibling species D. melanogaster and D. mauritiana (approximately 2.7 million years divergent). The level of sequence divergence between the two species within this region was significantly lower than expected from the neutral substitution rate and lower than the divergence observed between a randomly selected intron of the Drosophila Alcohol dehydrogenase gene (Adh).

CONCLUSIONS

Our results suggest that a 359 bp fragment of an Antonia retrotransposon (complete LTR is 659 bp) located within the intron of the Drosophila melanogaster Cht3 gene is of adaptive evolutionary significance. Our results are consistent with previous suggestions that the presence of TEs in constitutive heterochromatin may be of significance to the expression of heterochromatic genes.

The transposable elements of the Drosophila melanogaster euchromatin: a genomics perspective

BACKGROUND

Transposable elements are found in the genomes of nearly all eukaryotes. The recent completion of the Release 3 euchromatic genomic sequence of Drosophila melanogaster by the Berkeley Drosophila Genome Project has provided precise sequence for the repetitive elements in the Drosophila euchromatin. We have used this genomic sequence to describe the euchromatic transposable elements in the sequenced strain of this species.

RESULTS

We identified 85 known and eight novel families of transposable element varying in copy number from one to 146. A total of 1,572 full and partial transposable elements were identified, comprising 3.86% of the sequence. More than two-thirds of the transposable elements are partial. The density of transposable elements increases an average of 4.7 times in the centromere-proximal regions of each of the major chromosome arms. We found that transposable elements are preferentially found outside genes; only 436 of 1,572 transposable elements are contained within the 61.4 Mb of sequence that is annotated as being transcribed. A large proportion of transposable elements is found nested within other elements of the same or different classes. Lastly, an analysis of structural variation from different families reveals distinct patterns of deletion for elements belonging to different classes.

CONCLUSIONS

This analysis represents an initial characterization of the transposable elements in the Release 3 euchromatic genomic sequence of D. melanogaster for which comparison to the transposable elements of other organisms can begin to be made. These data have been made available on the Berkeley Drosophila Genome Project website for future analyses.