Retrotransposon populations of Vicia species with varying genome size

Genome features of common vetch (Vicia sativa) in natural habitats

Wild plants are often tolerant to biotic and abiotic stresses in their natural environments, whereas domesticated plants such as crops frequently lack such resilience. This difference is thought to be due to the high levels of genome heterozygosity in wild plant populations and the low levels of heterozygosity in domesticated crop species. In this study, common vetch (Vicia sativa) was used as a model to examine this hypothesis. The common vetch genome (2n = 14) was estimated as 1.8 Gb in size. Genome sequencing produced a reference assembly that spanned 1.5 Gb, from which 31,146 genes were predicted. Using this sequence as a reference, 24,118 single nucleotide polymorphisms were discovered in 1243 plants from 12 natural common vetch populations in Japan. Common vetch genomes exhibited high heterozygosity at the population level, with lower levels of heterozygosity observed at specific genome regions. Such patterns of heterozygosity are thought to be essential for adaptation to different environments. The resources generated in this study will provide insights into de novo domestication of wild plants and agricultural enhancement.

Retrotransposons and the Evolution of Genome Size in Pisum

Here we investigate the plant population genetics of retrotransposon insertion sites in pea to find out whether genetic drift and the neutral theory of molecular evolution can account for their abundance in the pea genome. (1) We asked whether two contrasting types of pea LTR-containing retrotransposons have the frequency and age distributions consistent with the behavior of neutral alleles and whether these parameters can explain the rate of change of genome size in legumes. (2) We used the recently assembled v1a pea genome sequence to obtain data on LTR-LTR divergence from which their age can be estimated. We coupled these data to prior information on the distribution of insertion site alleles. (3) We found that the age and frequency distribution data are consistent with the neutral theory. (4) We concluded that demographic processes are the underlying cause of genome size variation in legumes.

Obtaining retrotransposon sequences, analysis of their genomic distribution and use of retrotransposon-derived genetic markers in lentil (Lens culinaris Medik.)

Retrotransposons with long terminal repeats (LTR-RTs) are widespread mobile elements in eukaryotic genomes. We obtained a total of 81 partial LTR-RT sequences from lentil corresponding to internal retrotransposon components and LTRs. Sequences were obtained by PCR from genomic DNA. Approximately 37% of the LTR-RT internal sequences presented premature stop codons, pointing out that these elements must be non-autonomous. LTR sequences were obtained using the iPBS technique which amplifies sequences between LTR-RTs. A total of 193 retrotransposon-derived genetic markers, mainly iPBS, were used to obtain a genetic linkage map from 94 F₇ inbred recombinant lines derived from the cross between the cultivar Lupa and the wild ancestor L. culinaris subsp. orientalis. The genetic map included 136 markers located in eight linkage groups. Clusters of tightly linked retrotransposon-derived markers were detected in linkage groups LG1, LG2, and LG6, hence denoting a non-random genomic distribution. Phylogenetic analyses identified the LTR-RT families in which internal and LTR sequences are included. Ty3-gypsy elements were more frequent than Ty1-copia, mainly due to the high Ogre element frequency in lentil, as also occurs in other species of the tribe Vicieae. LTR and internal sequences were used to analyze in silico their distribution among the contigs of the lentil draft genome. Up to 8.8% of the lentil contigs evidenced the presence of at least one LTR-RT similar sequence. A statistical analysis suggested a non-random distribution of these elements within of the lentil genome. In most cases (between 97% and 72%, depending on the LTR-RT type) none of the internal sequences flanked by the LTR sequence pair was detected, suggesting that defective and non-autonomous LTR-RTs are very frequent in lentil. Results support that LTR-RTs are abundant and widespread throughout of the lentil genome and that they are a suitable source of genetic markers useful to carry out further genetic analyses.

The Small Nuclear Genomes of Selaginella Are Associated with a Low Rate of Genome Size Evolution

The haploid nuclear genome size (1C DNA) of vascular land plants varies over several orders of magnitude. Much of this observed diversity in genome size is due to the proliferation and deletion of transposable elements. To date, all vascular land plant lineages with extremely small nuclear genomes represent recently derived states, having ancestors with much larger genome sizes. The Selaginellaceae represent an ancient lineage with extremely small genomes. It is unclear how small nuclear genomes evolved in Selaginella We compared the rates of nuclear genome size evolution in Selaginella and major vascular plant clades in a comparative phylogenetic framework. For the analyses, we collected 29 new flow cytometry estimates of haploid genome size in Selaginella to augment publicly available data. Selaginella possess some of the smallest known haploid nuclear genome sizes, as well as the lowest rate of genome size evolution observed across all vascular land plants included in our analyses. Additionally, our analyses provide strong support for a history of haploid nuclear genome size stasis in Selaginella Our results indicate that Selaginella, similar to other early diverging lineages of vascular land plants, has relatively low rates of genome size evolution. Further, our analyses highlight that a rapid transition to a small genome size is only one route to an extremely small genome.

Analysis of transposable elements in the genome of Asparagus officinalis from high coverage sequence data

Asparagus officinalis is an economically and nutritionally important vegetable crop that is widely cultivated and is used as a model dioecious species to study plant sex determination and sex chromosome evolution. To improve our understanding of its genome composition, especially with respect to transposable elements (TEs), which make up the majority of the genome, we performed Illumina HiSeq2000 sequencing of both male and female asparagus genomes followed by bioinformatics analysis. We generated 17 Gb of sequence (12×coverage) and assembled them into 163,406 scaffolds with a total cumulated length of 400 Mbp, which represent about 30% of asparagus genome. Overall, TEs masked about 53% of the A. officinalis assembly. Majority of the identified TEs belonged to LTR retrotransposons, which constitute about 28% of genomic DNA, with Ty1/copia elements being more diverse and accumulated to higher copy numbers than Ty3/gypsy. Compared with LTR retrotransposons, non-LTR retrotransposons and DNA transposons were relatively rare. In addition, comparison of the abundance of the TE groups between male and female genomes showed that the overall TE composition was highly similar, with only slight differences in the abundance of several TE groups, which is consistent with the relatively recent origin of asparagus sex chromosomes. This study greatly improves our knowledge of the repetitive sequence construction of asparagus, which facilitates the identification of TEs responsible for the early evolution of plant sex chromosomes and is helpful for further studies on this dioecious plant.

Isolation and characterization of reverse transcriptase fragments of LTR retrotransposons from the genome of Chenopodium quinoa (Amaranthaceae)

High heterogeneity was observed among conserved domains of reverse transcriptase ( rt ) isolated from quinoa. Only one Ty1- copia rt was highly amplified. Reverse transcriptase sequences were located predominantly in pericentromeric region of quinoa chromosomes. The heterogeneity, genomic abundance, and chromosomal distribution of reverse transcriptase (rt)-coding fragments of Ty1-copia and Ty3-gypsy long terminal repeat retrotransposons were analyzed in the Chenopodium quinoa genome. Conserved domains of the rt gene were amplified and characterized using degenerate oligonucleotide primer pairs. Sequence analyses indicated that half of Ty1-copia rt (51 %) and 39 % of Ty3-gypsy rt fragments contained intact reading frames. High heterogeneity among rt sequences was observed for both Ty1-copia and Ty3-gypsy rt amplicons, with Ty1-copia more heterogeneous than Ty3-gypsy. Most of the isolated rt fragments were present in quinoa genome in low copy numbers, with only one highly amplified Ty1-copia rt sequence family. The gypsy-like RNase H fragments co-amplified with Ty1-copia-degenerate primers were shown to be highly amplified in the quinoa genome indicating either higher abundance of some gypsy families of which rt domains could not be amplified, or independent evolution of this gypsy-region in quinoa. Both Ty1-copia and Ty3-gypsy retrotransposons were preferentially located in pericentromeric heterochromatin of quinoa chromosomes. Phylogenetic analyses of newly amplified rt fragments together with well-characterized retrotransposon families from other organisms allowed identification of major lineages of retroelements in the genome of quinoa and provided preliminary insight into their evolutionary dynamics.

Assembly and comparative analysis of transposable elements from low coverage genomic sequence data in Asparagales

The research field of comparative genomics is moving from a focus on genes to a more holistic view including the repetitive complement. This study aimed to characterize relative proportions of the repetitive fraction of large, complex genomes in a nonmodel system. The monocotyledonous plant order Asparagales (onion, asparagus, agave) comprises some of the largest angiosperm genomes and represents variation in both genome size and structure (karyotype). Anonymous, low coverage, single-end Illumina data from 11 exemplar Asparagales taxa were assembled using a de novo method. Resulting contigs were annotated using a reference library of available monocot repetitive sequences. Mapping reads to contigs provided rough estimates of relative proportions of each type of transposon in the nuclear genome. The results were parsed into general repeat types and synthesized with genome size estimates and a phylogenetic context to describe the pattern of transposable element evolution among these lineages. The major finding is that although some lineages in Asparagales exhibit conservation in repeat proportions, there is generally wide variation in types and frequency of repeats. This approach is an appropriate first step in characterizing repeats in evolutionary lineages with a paucity of genomic resources.

Identification and characterization of jute LTR retrotransposons:: Their abundance, heterogeneity and transcriptional activity

Long Terminal Repeat (LTR) retrotransposons constitute a significant part of eukaryotic genomes and play an important role in genome evolution especially in plants. Jute is an important fiber crop with a large genome of 1,250 Mbps. This genome is still mostly unexplored. In this study we aimed at identifying and characterizing the LTR retrotransposons of jute with a view to understanding the jute genome better. In this study, the Reverse Transcriptase domain of Ty1-copia and Ty3-gypsy LTR retrotransposons of jute were amplified by degenerate primers and their expressions were examined by reverse transcription PCR. Copy numbers of reverse transcriptase (RT) genes of Ty1-copia and Ty3-gypsy elements were determined by dot blot analysis. Sequence analysis revealed higher heterogeneity among Ty1-copia retrotransposons than Ty3-gypsy and clustered each of them in three groups. Copy number of RT genes in Ty1-copia was found to be higher than that of Ty3-gypsy elements from dot blot hybridization. Cumulatively Ty1-copia and Ty3-gypsy may constitute around 19% of the jute genome where two groups of Ty1-copia were found to be transcriptionally active. Since the LTR retrotransposons constitute a large portion of jute genome, these findings imply the importance of these elements in the evolution of jute genome.

Two novel Ty1-copia retrotransposons isolated from coffee trees can effectively reveal evolutionary relationships in the Coffea genus (Rubiaceae)

In the study, we developed new markers for phylogenetic relationships and intraspecies differentiation in Coffea. Nana and Divo, two novel Ty1-copia LTR-retrotransposon families, were isolated through C. canephora BAC clone sequencing. Nana- and Divo-based markers were used to test their: (1) ability to resolve recent phylogenetic relationships; (2) efficiency in detecting intra-species differentiation. Sequence-specific amplification polymorphism (SSAP), retrotransposon-microsatellite amplified polymorphism (REMAP) and retrotransposon-based insertion polymorphism (RBIP) approaches were applied to 182 accessions (31 Coffea species and one Psilanthus accession). Nana- and Divo-based markers revealed contrasted transpositional histories. At the BAC clone locus, RBIP results on C. canephora demonstrated that Nana insertion took place prior to C. canephora differentiation, while Divo insertion occurred after differentiation. Combined SSAP and REMAP data showed that Nana could resolve Coffea lineages, while Divo was efficient at a lower taxonomic level. The combined results indicated that the retrotransposon-based markers were useful in highlighting Coffea genetic diversity and the chronological pattern of speciation/differentiation events. Ongoing complete sequencing of the C. canephora genome will soon enable exhaustive identification of LTR-RTN families, as well as more precise in-depth analyses on contributions to genome size variation and Coffea evolution.

Analysis of plant diversity with retrotransposon-based molecular markers

Retrotransposons are both major generators of genetic diversity and tools for detecting the genomic changes associated with their activity because they create large and stable insertions in the genome. After the demonstration that retrotransposons are ubiquitous, active and abundant in plant genomes, various marker systems were developed to exploit polymorphisms in retrotransposon insertion patterns. These have found applications ranging from the mapping of genes responsible for particular traits and the management of backcrossing programs to analysis of population structure and diversity of wild species. This review provides an insight into the spectrum of retrotransposon-based marker systems developed for plant species and evaluates the contributions of retrotransposon markers to the analysis of population diversity in plants.

Analysis of transposons and repeat composition of the sunflower (Helianthus annuus L.) genome

A sample-sequencing strategy combined with slot-blot hybridization and FISH was used to study the composition of the repetitive component of the sunflower genome. One thousand six hundred thirty-eight sequences for a total of 954,517 bp were analyzed. The fraction of sequences that can be classified as repetitive using computational and hybridization approaches amounts to 62% in total. Almost two thirds remain as yet uncharacterized in nature. Of those characterized, most belong to the gypsy superfamily of LTR-retrotransposons. Unlike in other species, where single families can account for large fractions of the genome, it appears that no transposon family has been amplified to very high levels in sunflower. All other known classes of transposable elements were also found. One family of unknown nature (contig 61) was the most repeated in the sunflower genome. The evolution of the repetitive component in the Helianthus genus and in other Asteraceae was studied by comparative analysis of the hybridization of total genomic DNAs from these species to the sunflower small-insert library and compared to gene-based phylogeny. Very little similarity is observed between Helianthus species and two related Asteraceae species outside of the genus. Most repetitive elements are similar in annual and perennial Helianthus species indicating that sequence amplification largely predates such divergence. Gypsy-like elements are more represented in the annuals than in the perennials, while copia-like elements are similarly represented, attesting a different amplification history of the two superfamilies of LTR-retrotransposons in the Helianthus genus.

An abundant and heavily truncated non-LTR retrotransposon (LINE) family in Beta vulgaris

We describe a non-LTR retrotransposon family,BvL, of the long interspersed nuclear elements L1 clade isolated from sugar beet (Beta vulgaris). Characteristic molecular domains of three full-length BvL elements were determined in detail, showing that coding sequences are interrupted and most likely non-functionally. In addition,eight highly conserved endonuclease regions were defined by comparison with other plant LINEs. The abundant BvL family is widespread within the genus Beta, however, the vast majority of BvL copies are extremely 50 truncated indicating an error-prone reverse transcriptase activity. The dispersed distribution of BvL copies on all sugar beet chromosomes with exclusion of most heterochromatic regions was shown by fluorescent in situ hybridization. The analysis of BvL 30 end sequences and corresponding flanking regions, respectively, revealed the preferred integration of BvL into A/T-rich regions of the sugar beet genome, but no specific target sequences.

Rapid DNA loss as a counterbalance to genome expansion through retrotransposon proliferation in plants

Transposable elements, particularly LTR-retrotransposons, comprise the primary vehicle for genome size expansion in plants, while DNA removal through illegitimate recombination and intrastrand homologous recombination serve as the most important counteracting forces to plant genomic obesity. Despite extensive research, the relative impact of these opposing forces and hence the directionality of genome size change remains unknown. In Gossypium (cotton), the 3-fold genome size variation among diploids is due largely to copy number variation of the gypsy-like retrotransposon Gorge3. Here we combine comparative sequence analysis with a modeling approach to study the directionality of genome size change in Gossypium. We demonstrate that the rate of DNA removal in the smaller genomes is sufficient to reverse genome expansion through Gorge3 proliferation. These data indicate that rates of DNA loss can be highly variable even within a single plant genus, and that the known mechanisms of DNA loss can indeed reverse the march toward genomic obesity.

Evolution of genome size and complexity in Pinus

BACKGROUND

Genome evolution in the gymnosperm lineage of seed plants has given rise to many of the most complex and largest plant genomes, however the elements involved are poorly understood.

METHODOLOGY/PRINCIPAL FINDINGS

Gymny is a previously undescribed retrotransposon family in Pinus that is related to Athila elements in Arabidopsis. Gymny elements are dispersed throughout the modern Pinus genome and occupy a physical space at least the size of the Arabidopsis thaliana genome. In contrast to previously described retroelements in Pinus, the Gymny family was amplified or introduced after the divergence of pine and spruce (Picea). If retrotransposon expansions are responsible for genome size differences within the Pinaceae, as they are in angiosperms, then they have yet to be identified. In contrast, molecular divergence of Gymny retrotransposons together with other families of retrotransposons can account for the large genome complexity of pines along with protein-coding genic DNA, as revealed by massively parallel DNA sequence analysis of Cot fractionated genomic DNA.

CONCLUSIONS/SIGNIFICANCE

Most of the enormous genome complexity of pines can be explained by divergence of retrotransposons, however the elements responsible for genome size variation are yet to be identified. Genomic resources for Pinus including those reported here should assist in further defining whether and how the roles of retrotransposons differ in the evolution of angiosperm and gymnosperm genomes.

Genome-wide characterization of long terminal repeat -retrotransposons in apple reveals the differences in heterogeneity and copy number between Ty1-copia and Ty3-gypsy retrotransposons

The conserved domains of reverse transcriptase (RT) genes of Ty1-copia and Ty3-gypsy groups of long terminal repeat (LTR) retrotransposons were isolated from the Malus domestica genome using degenerate oligonucleotide primers. Sequence analysis showed that 45% of Ty1-copia and 63% of Ty3-gypsy RT sequences contained premature stop codons and/or indels disrupting the reading frame. High heterogeneity among RT sequences of both Ty1-copia and Ty3-gypsy group retrotransposons was observed, but Ty3-gypsy group retrotransposons in the apple genome are less heterogeneous than Ty1-copia elements. Retrotransposon copy number was estimated by dot blot hybridizations for Ty1-copia (approximately 5,000) and Ty3-gypsy (approximately 26,000). All elements of the two types of LTR retrotransposons comprise approximately 38% of the M. domestica genome, with the Ty3-gypsy group contribution being higher (33.5%) than the Ty1-copia one (4.6%). Transcription was not detected by reverse transcription-polymerase chain reaction for either Ty1-copia or Ty3-gypsy retrotransposons in the leaves of plants in vitro or in leaf explants cultured on medium supplemented with high concentration benzylaminopurine. This research reveals the differences in heterogeneity and copy number between Ty1-copia and Ty3-gypsy retrotransposons in the apple genome. Ty1-copia retrotransposon has higher heterogeneity than Ty3-gypsy retrotransposon, but the latter has a higher copy number, which implies that Ty3-gypsy retrotransposons may play a more important role in the apple genome evolution.

The pattern of amplification and differentiation of Ty1-copia and Ty3-gypsy retrotransposons in Brassicaceae species

One of the causes of genome size expansion is considered to be amplification of retrotransposons. We determined nucleotide sequences of 24 PCR products for each of six retrotransposons in Brassica rapa and Brassica oleracea. Phylogenetic trees of these sequences showed species-specific clades. We also sequenced STF7a homologs and Tto1 homologs, 24 PCR products each, in nine diploids and three allopolyploids, and constructed phylogenetic trees. In these phylogenetic trees, species-specific clades of diploid species were also formed, but retrotransposons of allopolyploids were clustered into the clades of their original genomes, indicating that these two retrotransposons amplified after speciation of the nine diploids. Genetic variation in these retrotransposons may have arisen before emergence of allopolyploid species. There was a positive correlation between the genome size and the average number of substitutions of STF7a and Tto1 homologs in at least seven diploids. The implications of these results in the genome evolution of Brassicaceae are herein discussed.

A phylogenetic analysis of indel dynamics in the cotton genus

Genome size evolution is a dynamic process involving counterbalancing mechanisms whose actions vary across lineages and over time. Whereas the primary mechanism of expansion, transposable element (TE) amplification, has been widely documented, the evolutionary dynamics of genome contraction have been less thoroughly explored. To evaluate the relative impact and evolutionary stability of the mechanisms that affect genome size, we conducted a phylogenetic analysis of indel rates for 2 genomic regions in 4 Gossypium genomes: the 2 coresident genomes (A(T) and D(T)) of tetraploid cotton and its model diploid progenitors, Gossypium arboreum (A) and Gossypium raimondii (D). We determined the rates of sequence gain or loss along each branch, partitioned by mechanism, and how these changed during species divergence. In general, there has been a propensity toward growth of the diploid genomes and contraction in the polyploid. Most of the size difference between the diploid species occurred prior to polyploid divergence and was largely attributable to TE amplification in the A/A(T) genome. After separating from the true parents of the polyploid genomes, both diploid genomes experienced slower sequence gain than in the ancestor, due to fewer TE insertions in the A genome and a combination of increased deletions and decreased TE insertions in the D genome. Both genomes of the polyploid displayed increased rates of deletion and decreased rates of insertion, leading to a rate of near stasis in D(T) and overall contraction in A(T) resulting in polyploid genome contraction. As expected, TE insertions contributed significantly to the genome size differences; however, intrastrand homologous recombination, although rare, had the most significant impact on the rate of deletion. Small indel data for the diploids suggest the possibility of a bias as the smaller genomes add less or delete more sequence through small indels than do the larger genomes, whereas data for the polyploid suggest increased sequence turnover in general (both as small deletions and small insertions). Illegitimate recombination, although not demonstrated to be a dominant mechanism of genome size change, was biased in the polyploid toward deletions, which may provide a partial explanation of polyploid genomic downsizing.

Phylogenetic determination of the pace of transposable element proliferation in plants: copia and LINE-like elements in Gossypium

Transposable elements contribute significantly to plant genome evolution in myriad ways, ranging from local insertional mutations to global effects exerted on genome size through accumulation. Differential accumulation and deletion of transposable elements may profoundly affect genome size, even among members of the same genus. One example is that of Gossypium (cotton), where much of the 3-fold genome size variation is due to differential accumulation of one gypsy-like LTR retrotransposon, Gorge3. Copia and non-LTR LINE retrotransposons are also major components of the Gossypium genome, but unlike Gorge3, their extant copy numbers do not correlate with genome size. In the present study, we describe the nature and timing of transposition for copia and LINE retrotransposons in Gossypium. Our findings indicate that copia retrotransposons have been active in each lineage since divergence from a common ancestor, and that they have proliferated in a punctuated manner. However, the evolutionary history of LINEs contrasts markedly with that of the copia retrotransposons. Although LINEs have also been active in each lineage, they have accumulated in a stochastically regular manner, and phylogenetic analysis suggests that extant LINE populations in Gossypium are dominated by ancient insertions. Interestingly, the magnitude of transpositional bursts in each lineage corresponds directly with extant estimated copy number.

Isolation and characterization of genomic retrotransposon sequences from octoploid strawberry (Fragaria x ananassa Duch.)

Strawberry (Fragaria spp.) is a kind of herbaceous perennial plant that propagates vegetatively. The conserved domains of reverse transcriptase (RT) genes of Ty1-copia and Ty3-gypsy groups of LTR retrotransposons were amplified from the cultivated strawberry (Fragaria x ananassa Duch.). Sequence analysis of clones demonstrated that 5 of 19 Ty1-copia group unique sequences and 2 of 10 Ty3-gypsy unique sequences in F. x ananassa genome possessed either stop codon or frameshift. Ty1-copia group sequences are highly heterogeneous (divergence ranged from 1 to 69.8%), but the Ty3-gypsy group sequences are less (divergence ranged from 1 to 10%). Southern dot blot hybridization result suggested that both of the LTR retrotransposons are present in the genome of cultivated strawberry with high copy number (Ty1-copia group 2,875 Ty3-gypsy group 348). RT-PCR amplification from total RNA, which was extracted from leaves of micropropagated strawberry plants, did not yield either of the RT fragments. This is the first report on the presence of RT sequences of Ty1-copia and Ty3-gypsy group retrotransposons in F. x ananassa genome.

The heterochromatic copies of the LTR retrotransposons as a record of the genomic events that have shaped the Drosophila melanogaster genome

Transposable elements, which are major components of most genomes, are known to accumulate in heterochromatic regions in which they have progressively diverged in sequence by mutations and internal deletions and insertions (indels) during the course of evolution. They therefore provide a record of the genomic events that have shaped the genomes, some of which could correspond to speciation events. Using the sequence divergence between the long terminal repeats (LTRs), we estimated the date of the insertion events of the LTR retrotransposon copies embedded within the heterochromatin regions of the Drosophila melanogaster genome. We did not detect traces of any specific waves of mobilization of retrotransposons within heterochromatin, apart from a very recent wave, which corresponds to the numerous LTR retrotransposon copies found in euchromatin.

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.

Transposable element distribution, abundance and role in genome size variation in the genus Oryza

Background

The genus Oryza is composed of 10 distinct genome types, 6 diploid and 4 polyploid, and includes the world's most important food crop – rice (Oryza sativa [AA]). Genome size variation in the Oryza is more than 3-fold and ranges from 357 Mbp in Oryza glaberrima [AA] to 1283 Mbp in the polyploid Oryza ridleyi [HHJJ]. Because repetitive elements are known to play a significant role in genome size variation, we constructed random sheared small insert genomic libraries from 12 representative Oryza species and conducted a comprehensive study of the repetitive element composition, distribution and phylogeny in this genus. Particular attention was paid to the role played by the most important classes of transposable elements (Long Terminal Repeats Retrotransposons, Long interspersed Nuclear Elements, helitrons, DNA transposable elements) in shaping these genomes and in their contributing to genome size variation.

Results

We identified the elements primarily responsible for the most strikingly genome size variation in Oryza. We demonstrated how Long Terminal Repeat retrotransposons belonging to the same families have proliferated to very different extents in various species. We also showed that the pool of Long Terminal Repeat Retrotransposons is substantially conserved and ubiquitous throughout the Oryza and so its origin is ancient and its existence predates the speciation events that originated the genus. Finally we described the peculiar behavior of repeats in the species Oryza coarctata [HHKK] whose placement in the Oryza genus is controversial.

Conclusion

Long Terminal Repeat retrotransposons are the major component of the Oryza genomes analyzed and, along with polyploidization, are the most important contributors to the genome size variation across the Oryza genus. Two families of Ty3-gypsy elements (RIRE2 and Atlantys) account for a significant portion of the genome size variations present in the Oryza genus.

Diversity and evolution of Ty1-copia and Ty3-gypsy retroelements in the non-photosynthetic flowering plants Orobanche and Phelipanche (Orobanchaceae)

We present the first study on the diversity and evolution of Ty1-copia and Ty3-gypsy retroelements in a group of non-photosynthetic flowering plants. To this end partial sequences of the reverse transcriptase (rt) gene were obtained from 20 clones for each retroelement type from seven and six accessions of Orobanche and Phelipanche (Orobanchaceae), respectively. Overall sequence similarity is higher in Ty3-gypsy elements than in Ty1-copia elements in agreement with the results from other angiosperm groups. Higher sequence diversity and stronger phylogenetic structure, especially of Ty1-copia sequences, in Orobanche species compared to Phelipanche species support the previously suggested hypothesis (based on karyological and cytological data) that genomes of Orobanche species are more dynamic than those of Phelipanche species. No evidence was found for intraspecific differences of retroelement diversity nor for differences between pest taxa and their putative wild relatives, e.g., O. crenata and O. owerini. The occurrence of a few sequences from Phelipanche species in clades otherwise comprising sequences from Orobanche species might be due to horizontal gene transfer, but the alternative of vertical transmission cannot be rejected unambiguously.

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.

Phylogenetic analyses within three sections of the genus Vicia

The averaged genomic similarities based on multilocus randomly amplified polymorphic DNA (RAPD) were calculated for eight species representing three sections of the genus Vicia: faba, bithynica and narbonensis. The frequency of appearance of the sequences corresponding to 25 decamers selected at random from genomes of different Fabace species was checked, and a high correlation with the frequency observed for Vicia allowed us to assume their similar weight in typing Vicia species. The RAPD-based similarity coefficients compared with those related to whole genome hybridization with barley rDNA and those based on restriction fragment length polymorphism (RFLP) revealed similar interspecies relationships. The averaged RAPD-based similarity coefficient (Pearson's) was 0.68 for all the species, and was sectionspecific: 0.43 (bithynica), 0.50 (faba) and 0.73 (narbonensis). The averaged similarity coefficient for V. serratifolia (0.63) placed it apart from the rest (0.75) of its section. The results correspond to the interspecies relationships built upon non-genetic data. The averaged similarity coefficient for particular RAPD was related to the presence and type of tandemly repeated motif in a primer: 0.7-0.8 for heterodimers (GC, AG, CA, GT, CT), 0.5-0.6 for homodimers (CC, GG) and 0.6 for no repeat, indicating the sensitivity of diversity range to the type of target sequences.

Insertional polymorphism and antiquity of PDR1 retrotransposon insertions in pisum species

Sequences flanking 73 insertions of the retrotransposon PDR1 have been characterized, together with an additional 270 flanking regions from one side alone, from a diverse collection of Pisum germ plasm. Most of the identified flanking sequences are repetitious DNAs but more than expected (7%) lie within nuclear gene protein-coding regions. The approximate age of 52 of the PDR1 insertions has been determined by measuring sequence divergence among LTR pairs. These data show that PDR1 transpositions occurred within the last 5 MY, with a peak at 1-2.5 MYA. The insertional polymorphism of 68 insertions has been assessed across 47 selected Pisum accessions, representing the diversity of the genus. None of the insertions are fixed, showing that PDR1 insertions can persist in a polymorphic state for millions of years in Pisum. The insertional polymorphism data have been compared with the age estimations to ask what rules control the proliferation of PDR1 insertions in Pisum. Relatively recent insertions (< approximately 1.5 MYA) tend to be found in small subsets of the Pisum accessions set, "middle-aged" insertions (between approximately 1.5 and 2.5 MYA) vary greatly in their occurrence, and older insertions (> approximately 2.5 MYA) are mostly found in small subsets of Pisum. Finally, the average age estimate for PDR1 insertions, together with an existing data set for PDR1 retrotransposon SSAP markers, has been used to derive an estimate of the effective population size for Pisum of approximately 7.5 x 10(5).