PIGY, a new plant envelope-class LTR retrotransposon

A widespread occurrence of extra open reading frames in plant Ty3/gypsy retrotransposons

Long terminal repeat (LTR) retrotransposons make up substantial parts of most higher plant genomes where they accumulate due to their replicative mode of transposition. Although the transposition is facilitated by proteins encoded within the gag-pol region which is common to all autonomous elements, some LTR retrotransposons were found to potentially carry an additional protein coding capacity represented by extra open reading frames located upstream or downstream of gag-pol. In this study, we performed a comprehensive in silico survey and comparative analysis of these extra open reading frames (ORFs) in the group of Ty3/gypsy LTR retrotransposons as the first step towards our understanding of their origin and function. We found that extra ORFs occur in all three major lineages of plant Ty3/gypsy elements, being the most frequent in the Tat lineage where most (77 %) of identified elements contained extra ORFs. This lineage was also characterized by the highest diversity of extra ORF arrangement (position and orientation) within the elements. On the other hand, all of these ORFs could be classified into only two broad groups based on their mutual similarities or the presence of short conserved motifs in their inferred protein sequences. In the Athila lineage, the extra ORFs were confined to the element 3' regions but they displayed much higher sequence diversity compared to those found in Tat. In the lineage of Chromoviruses the extra ORFs were relatively rare, occurring only in 5' regions of a group of elements present in a single plant family (Poaceae). In all three lineages, most extra ORFs lacked sequence similarities to characterized gene sequences or functional protein domains, except for two Athila-like elements with similarities to LOGL4 gene and part of the Chromoviruses extra ORFs that displayed partial similarity to histone H3 gene. Thus, in these cases the extra ORFs most likely originated by transduction or recombination of cellular gene sequences. In addition, the protein domain which is otherwise associated with DNA transposons have been detected in part of the Tat-like extra ORFs, pointing to their origin from an insertion event of a mobile element.

Plant centromeric retrotransposons: a structural and cytogenetic perspective

Background

The centromeric and pericentromeric regions of plant chromosomes are colonized by Ty3/gypsy retrotransposons, which, on the basis of their reverse transcriptase sequences, form the chromovirus CRM clade. Despite their potential importance for centromere evolution and function, they have remained poorly characterized. In this work, we aimed to carry out a comprehensive survey of CRM clade elements with an emphasis on their diversity, structure, chromosomal distribution and transcriptional activity.

Results

We have surveyed a set of 190 CRM elements belonging to 81 different retrotransposon families, derived from 33 host species and falling into 12 plant families. The sequences at the C-terminus of their integrases were unexpectedly heterogeneous, despite the understanding that they are responsible for targeting to the centromere. This variation allowed the division of the CRM clade into the three groups A, B and C, and the members of each differed considerably with respect to their chromosomal distribution. The differences in chromosomal distribution coincided with variation in the integrase C-terminus sequences possessing a putative targeting domain (PTD). A majority of the group A elements possess the CR motif and are concentrated in the centromeric region, while members of group C have the type II chromodomain and are dispersed throughout the genome. Although representatives of the group B lack a PTD of any type, they appeared to be localized preferentially in the centromeres of tested species. All tested elements were found to be transcriptionally active.

Conclusions

Comprehensive analysis of the CRM clade elements showed that genuinely centromeric retrotransposons represent only a fraction of the CRM clade (group A). These centromeric retrotransposons represent an active component of centromeres of a wide range of angiosperm species, implying that they play an important role in plant centromere evolution. In addition, their transcriptional activity is consistent with the notion that the transcription of centromeric retrotransposons has a role in normal centromere function.

Natural history of eukaryotic DNA methylation systems

Methylation of cytosines and adenines in DNA is a widespread epigenetic mark in both prokaryotes and eukaryotes. In eukaryotes, it has a profound influence on chromatin structure and dynamics. Recent advances in genomics and biochemistry have considerably elucidated the functions and provenance of these DNA modifications. DNA methylases appear to have emerged first in bacterial restriction-modification (R-M) systems from ancient RNA-modifying enzymes, in transitions that involved acquisition of novel catalytic residues and DNA-recognition features. DNA adenine methylases appear to have been acquired by ciliates, heterolobosean amoeboflagellates, and certain chlorophyte algae. Six distinct clades of cytosine methylases, including the DNMT1, DNMT2, and DNMT3 clades, were acquired by eukaryotes through independent lateral transfer of their precursors from bacteria or bacteriophages. In addition to these, multiple adenine and cytosine methylases were acquired by several families of eukaryotic transposons. In eukaryotes, the DNA-methylase module was often combined with distinct modified and unmodified peptide recognition domains and other modules mediating specialized interactions, for example, the RFD module of DNMT1 which contains a permuted Sm domain linked to a helix-turn-helix domain. In eukaryotes, the evolution of DNA methylases appears to have proceeded in parallel to the elaboration of histone-modifying enzymes and the RNAi system, with functions related to counter-viral and counter-transposon defense, and regulation of DNA repair and differential gene expression being their primary ancestral functions. Diverse DNA demethylation systems that utilize base-excision repair via DNA glycosylases and cytosine deaminases appear to have emerged in multiple eukaryotic lineages. Comparative genomics suggests that the link between cytosine methylation and DNA glycosylases probably emerged first in a novel R-M system in bacteria. Recent studies suggest that the 5mC is not a terminal DNA modification, with enzymes of the Tet/JBP family of 2-oxoglutarate- and iron-dependent dioxygenases further hydroxylating it to form 5-hydroxymethylcytosine (5hmC). These enzymes emerged first in bacteriophages and appear to have been transferred to eukaryotes on one or more occasions. Eukaryotes appear to have recruited three major types of DNA-binding domains (SRA/SAD, TAM/MBD, and CXXC) in discriminating DNA with methylated or unmethylated cytosines. Analysis of the domain architectures of these domains and the DNA methylases suggests that early in eukaryotic evolution they developed a close functional link with SET-domain methylases and Jumonji-related demethylases that operate on peptides in chromatin proteins. In several eukaryotes, other functional connections were elaborated in the form of various combinations between domains related to DNA methylation and those involved in ATP-dependent chromatin remodeling and RNAi. In certain eukaryotes, such as mammals and angiosperms, novel dependencies on the DNA methylation system emerged, which resulted in it affecting unexpected aspects of the biology of these organisms such as parent-offspring interactions. In genomic terms, this was reflected in the emergence of new proteins related to methylation, such as Stella. The well-developed methylation systems of certain heteroloboseans, stramenopiles, chlorophytes, and haptophyte indicate that these might be new model systems to explore the relevance of DNA modifications in eukaryotes.

FIDEL-a retrovirus-like retrotransposon and its distinct evolutionary histories in the A- and B-genome components of cultivated peanut

In this paper, we describe a Ty3-gypsy retrotransposon from allotetraploid peanut (Arachis hypogaea) and its putative diploid ancestors Arachis duranensis (A-genome) and Arachis ipaënsis (B-genome). The consensus sequence is 11,223 bp. The element, named FIDEL (Fairly long Inter-Dispersed Euchromatic LTR retrotransposon), is more frequent in the A- than in the B-genome, with copy numbers of about 3,000 (±950, A. duranensis), 820 (±480, A. ipaënsis), and 3,900 (±1,500, A. hypogaea) per haploid genome. Phylogenetic analysis of reverse transcriptase sequences showed distinct evolution of FIDEL in the ancestor species. Fluorescent in situ hybridization revealed disperse distribution in euchromatin and absence from centromeres, telomeric regions, and the nucleolar organizer region. Using paired sequences from bacterial artificial chromosomes, we showed that elements appear less likely to insert near conserved ancestral genes than near the fast evolving disease resistance gene homologs. Within the Ty3-gypsy elements, FIDEL is most closely related with the Athila/Calypso group of retrovirus-like retrotransposons. Putative transmembrane domains were identified, supporting the presence of a vestigial envelope gene. The results emphasize the importance of FIDEL in the evolution and divergence of different Arachis genomes and also may serve as an example of the role of retrotransposons in the evolution of legume genomes in general.

Sequence heterogeneity of the envelope-like domain in cultivated allotetraploid Gossypium species and their diploid progenitors

Retroviral envelope (env)-like sequences in 2 cultivated allotetraploid cottons and their diploid progenitors have been identified and characterized in this study. DNA sequence analysis reveals that these sequences are heterogeneous. The observed sequence diversity, however, seems to preserve coding information. This is evidenced by the detection of the transmembrane domain (TM), which is the most conserved feature of the divergent retroviral env genes. The high ratio of synonymous to nonsynonymous changes suggests that these sequences are evolving under purifying selection. Phylogenetic analysis shows that Gossypium sequences closely cluster with a lineage of plant endogenous retroviruses that have an env-like gene. These results provide evidence for the antiquity and the wide diversity of env-like sequences in the Gossypium genome.

Envelope-like retrotransposons in the plant kingdom: evidence of their presence in gymnosperms (Pinus pinaster)

Retroviruses differ from retrotransposons due to their infective capacity, which depends critically on the encoded envelope. Some plant retroelements contain domains reminiscent of the env of animal retroviruses but the number of such elements described to date is restricted to angiosperms. We show here the first evidence of the presence of putative env-like gene sequences in a gymnosperm species, Pinus pinaster (maritime pine). Using a degenerate primer approach for conserved domains of RNaseH gene, three clones from putative envelope-like retrotransposons (PpRT2, PpRT3, and PpRT4) were identified. The env-like sequences of P. pinaster clones are predicted to encode proteins with transmembrane domains. These sequences showed identity scores of up to 30% with env-like sequences belonging to different organisms. A phylogenetic analysis based on protein alignment of deduced aminoacid sequences revealed that these clones clustered with env-containing plant retrotransposons, as well as with retrotransposons from invertebrate organisms. The differences found among the sequences of maritime pine clones isolated here suggest the existence of different putative classes of env-like retroelements. The identification for the first time of env-like genes in a gymnosperm species may support the ancestrality of retroviruses among plants shedding light on their role in plant evolution.

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.

Repetitive DNA in the pea (Pisum sativum L.) genome: comprehensive characterization using 454 sequencing and comparison to soybean and Medicago truncatula

BACKGROUND

Extraordinary size variation of higher plant nuclear genomes is in large part caused by differences in accumulation of repetitive DNA. This makes repetitive DNA of great interest for studying the molecular mechanisms shaping architecture and function of complex plant genomes. However, due to methodological constraints of conventional cloning and sequencing, a global description of repeat composition is available for only a very limited number of higher plants. In order to provide further data required for investigating evolutionary patterns of repeated DNA within and between species, we used a novel approach based on massive parallel sequencing which allowed a comprehensive repeat characterization in our model species, garden pea (Pisum sativum).

RESULTS

Analysis of 33.3 Mb sequence data resulted in quantification and partial sequence reconstruction of major repeat families occurring in the pea genome with at least thousands of copies. Our results showed that the pea genome is dominated by LTR-retrotransposons, estimated at 140,000 copies/1C. Ty3/gypsy elements are less diverse and accumulated to higher copy numbers than Ty1/copia. This is in part due to a large population of Ogre-like retrotransposons which alone make up over 20% of the genome. In addition to numerous types of mobile elements, we have discovered a set of novel satellite repeats and two additional variants of telomeric sequences. Comparative genome analysis revealed that there are only a few repeat sequences conserved between pea and soybean genomes. On the other hand, all major families of pea mobile elements are well represented in M. truncatula.

CONCLUSION

We have demonstrated that even in a species with a relatively large genome like pea, where a single 454-sequencing run provided only 0.77% coverage, the generated sequences were sufficient to reconstruct and analyze major repeat families corresponding to a total of 35-48% of the genome. These data provide a starting point for further investigations of legume plant genomes based on their global comparative analysis and for the development of more sophisticated approaches for data mining.

Repetitive DNA in the pea (Pisum sativum L.) genome: comprehensive characterization using 454 sequencing and comparison to soybean and Medicago truncatula

BACKGROUND

Extraordinary size variation of higher plant nuclear genomes is in large part caused by differences in accumulation of repetitive DNA. This makes repetitive DNA of great interest for studying the molecular mechanisms shaping architecture and function of complex plant genomes. However, due to methodological constraints of conventional cloning and sequencing, a global description of repeat composition is available for only a very limited number of higher plants. In order to provide further data required for investigating evolutionary patterns of repeated DNA within and between species, we used a novel approach based on massive parallel sequencing which allowed a comprehensive repeat characterization in our model species, garden pea (Pisum sativum).

RESULTS

Analysis of 33.3 Mb sequence data resulted in quantification and partial sequence reconstruction of major repeat families occurring in the pea genome with at least thousands of copies. Our results showed that the pea genome is dominated by LTR-retrotransposons, estimated at 140,000 copies/1C. Ty3/gypsy elements are less diverse and accumulated to higher copy numbers than Ty1/copia. This is in part due to a large population of Ogre-like retrotransposons which alone make up over 20% of the genome. In addition to numerous types of mobile elements, we have discovered a set of novel satellite repeats and two additional variants of telomeric sequences. Comparative genome analysis revealed that there are only a few repeat sequences conserved between pea and soybean genomes. On the other hand, all major families of pea mobile elements are well represented in M. truncatula.

CONCLUSION

We have demonstrated that even in a species with a relatively large genome like pea, where a single 454-sequencing run provided only 0.77% coverage, the generated sequences were sufficient to reconstruct and analyze major repeat families corresponding to a total of 35-48% of the genome. These data provide a starting point for further investigations of legume plant genomes based on their global comparative analysis and for the development of more sophisticated approaches for data mining.

The centromeric retrotransposons of rice are transcribed and differentially processed by RNA interference

Retrotransposons consist of significant portions of many complex eukaryotic genomes and are often enriched in heterochromatin. The centromeric retrotransposon (CR) family in grass species is colonized in the centromeres and highly conserved among species that have been diverged for >50 MY. These unique characteristics have inspired scientists to speculate about the roles of CR elements in organization and function of centromeric chromatin. Here we report that the CRR (CR of rice) elements in rice are highly enriched in chromatin associated with H3K9me2, a hallmark for heterochromatin. CRR elements were transcribed in root, leaf, and panicle tissues, suggesting a constitutive transcription of this retrotransposon family. However, the overall transcription level was low and the CRR transcripts appeared to be derived from relatively few loci. The majority of the CRR transcripts had chimerical structures and contained only partial CRR sequences. We detected small RNAs (smRNAs) cognate to nonautonomous CRR1 (noaCRR1) and CRR1, but not CRR2 elements. This result was also confirmed by in silico analysis of rice smRNA sequences. These results suggest that different CRR subfamilies may play different roles in the RNAi-mediated pathway for formation and maintenance of centromeric heterochromatin.

SIRE-1, a putative plant retrovirus is closely related to a legume TY1-copia retrotransposon family

SIRE-1 is a potential soybean retrovirus which has a gene order similar to Ty1-copia retrotransposons but also contains an envelope-like open reading frame (ORF), which is characteristic of retroviruses. PCR and Southern analysis reveals that SIRE-1 is closely related to a legume-wide family of envelope-lacking Ty1-copia group retrotransposons which suggests that SIRE-1 was formed by the recent acquisition of an envelope gene by a Ty1-copia retrotransposon.

Ogre elements--a distinct group of plant Ty3/gypsy-like retrotransposons

Ogre elements are a group of LTR retrotransposons recently discovered in legume plants, where they constitute almost 40% of the genome in some species. They are exceptional in their size (reaching 25 kb) and possess several specific features, including an intron within a polyprotein-coding region, and an extra open reading frame (ORF1) encoding a protein of unknown function located upstream of the gag gene. Although these features make Ogres interesting for further research, identification of additional elements from a broader range of plant taxa has been complicated by the divergence of their sequences, preventing their detection using similarity-based searches. Here we report the results of structure-based computational searches for Ogre elements in available plant genomic sequences, which proved to be more efficient and revealed occurrences of Ogres in three families of dicot plants (Leguminosae, Solanaceae and Salicaceae). In addition, a representative set of 85 elements was retrieved from a model legume species Medicago truncatula. All identified full-length elements were used for comparative analysis, which showed that in spite of only little conservation of their nucleotide sequences, their protein domains were highly conserved, including several regions within ORF1. Further, the elements shared the same functional regions, including a primer binding site complementary to tRNA(arg), a conserved motif within a polypurine tract, and a putative intron between the pro and rt/rh coding domains. These findings, together with analysis of their phylogenetic relationship to other retrotransposons based on similarities of rt domains suggest that Ogre elements from different plant taxa have a common origin and thus constitute a distinct group of Ty3/gypsy retrotransposons.

Retand: a novel family of gypsy-like retrotransposons harboring an amplified tandem repeat

In this paper we describe a pair of novel Ty3/gypsy retrotransposons isolated from the dioecious plant Silene latifolia, consisting of a non-autonomous element Retand-1 (3.7 kb) and its autonomous partner Retand-2 (11.1 kb). These two elements have highly similar long terminal repeat (LTR) sequences but differ in the presence of the typical retroelement coding regions (gag-pol genes), most of which are missing in Retand-1. Moreover, Retand-2 contains two additional open reading frames in antisense orientation localized between the pol gene and right LTR. Retand transcripts were detected in all organs tested (leaves, flower buds and roots) which, together with the high sequence similarity of LTRs in individual elements, indicates their recent transpositional activity. The autonomous elements are similarly abundant (2,700 copies) as non-autonomous ones (2,100 copies) in S. latifolia genome. Retand elements are also present in other Silene species, mostly in subtelomeric heterochromatin regions of all chromosomes. The only exception is the subtelomere of the short arm of the Y chromosome in S. latifolia which is known to lack the terminal heterochromatin. An interesting feature of the Retand elements is the presence of a tandem repeat sequence, which is more amplified in the non-autonomous Retand-1.

Significant expansion of Vicia pannonica genome size mediated by amplification of a single type of giant retroelement

Amplification and eventual elimination of dispersed repeats, especially those of the retroelement origin, account for most of the profound size variability observed among plant genomes. In most higher plants investigated so far, differential accumulation of various families of elements contributes to these differences. Here we report the identification of giant Ty3/gypsy-like retrotransposons from the legume plant Vicia pannonica, which alone make up approximately 38% of the genome of this species. These retrotransposons have structural features of the Ogre elements previously identified in the genomes of pea and Medicago. These features include extreme size (25 kb), the presence of an extra ORF upstream of the gag-pol region, and a putative intron dividing the prot and rt coding sequences. The Ogre elements are evenly dispersed on V. pannonica chromosomes except for terminal regions containing satellite repeats, their individual copies show extraordinary sequence similarity, and at least part of them are transcriptionally active, which suggests their recent amplification. Similar elements were also detected in several other Vicia species but in most cases in significantly lower numbers. However, there was no obvious correlation of the abundance of Ogre sequences with the genome size of these species.

Development of an efficient retrotransposon-based fingerprinting method for rapid pea variety identification

Fast and efficient DNA fingerprinting of crop cultivars and individuals is frequently used in both theoretical population genetics and in practical breeding. Numerous DNA marker technologies exist and the ratio of speed, cost and accuracy are of importance. Therefore even in species where highly accurate and polymorphic marker systems are available, such as microsatellite SSR (simple sequence repeats), also alternative methods may be of interest. Thanks to their high abundance and ubiquity, temporary mobile retrotransposable elements come into recent focus. Their properties, such as genome wide distribution and well-defined origin of individual insertions by descent, predetermine them for use as molecular markers. In this study, several Ty3-gypsy type retrotransposons have been developed and adopted for the inter-retrotransposon amplified polymorphism (IRAP) method, which is suitable for fast and efficient pea cultivar fingerprinting. The method can easily distinguish even between genetically closely related pea cultivars and provide high polymorphic information content (PIC) in a single PCR analysis.