Evolutionary dynamics of an ancient retrotransposon family provides insights into evolution of genome size in the genus Oryza

A lineage-specific centromere retrotransposon in Oryza brachyantha

Most eukaryotic centromeres contain large quantities of repetitive DNA, such as satellite repeats and retrotransposons. Unlike most transposons in plant genomes, the centromeric retrotransposon (CR) family is conserved over long evolutionary periods among a majority of the grass species. CR elements are highly concentrated in centromeres, and are likely to play a role in centromere function. In order to study centromere evolution in the Oryza (rice) genus, we sequenced the orthologous region to centromere 8 of Oryza sativa from a related species, Oryza brachyantha. We found that O. brachyantha does not have the canonical CRR (CR of rice) found in the centromeres of all other Oryza species. Instead, a new Ty3-gypsy (Metaviridae) retroelement (FRetro3) was found to colonize the centromeres of this species. This retroelement is found in high copy numbers in the O. brachyantha genome, but not in other Oryza genomes, and based on the dating of long terminal repeats (LTRs) of FRetro3 it was amplified in the genome in the last few million years. Interestingly, there is a high level of removal of FRetro3 based on solo-LTRs to full-length elements, and this rapid turnover may have played a role in the replacement of the canonical CRR with the new element by active deletion. Comparison with previously described ChIP cloning data revealed that FRetro3 is found in CENH3-associated chromatin sequences. Thus, within a single lineage of the Oryza genus, the canonical component of grass centromeres has been replaced with a new retrotransposon that has all the hallmarks of a centromeric retroelement.

Diversity of the Ty-1 copia retrotransposon Tos17 in rice (Oryza sativa L.) and the AA genome of the Oryza genus

Retrotransposons are mobile genetic elements, ubiquitous in Eukaryotic genomes, which have proven to be major genetic tools in determining phylogeny and structuring genetic diversity, notably in plants. We investigate here the diversity of the Ty1-copia retrotransposon Tos17 in the cultivated rice of Asian origin (Oryza sativa L.) and related AA genome species of the Oryza genus, to contribute understanding of the complex evolutionary history in this group of species through that of the element in the lineages. In that aim, we used a combination of Southern hybridization with a reverse transcriptase (RT) probe and an adapter-PCR mediated amplification, which allowed the sequencing of the genomic regions flanking Tos17 insertions. This analysis was carried out in a collection of 47 A-genome Oryza species accessions and 202 accessions of a core collection of Oryza sativa L. representative of the diversity of the species. Our Southern hybridization results show that Tos17 is present in all the accessions of the A-genome Oryza species, except for the South American species O. glumaepatula and the African species O. glaberrima and O. breviligulata. In O. sativa, the number of putative copies of Tos17 per accession ranged from 1 to 11 and multivariate analysis based on presence/absence of putative copies yielded a varietal clustering which is consistent with the isozyme classification of rice. Adapter PCR amplification and sequencing of flanking regions of Tos17 insertions in A-genome species other than O. sativa, followed by anchoring on the Nipponbare genome sequence, revealed 13 insertion sites of Tos17 in the surveyed O. rufipogon and O. longistaminata accessions, including one shared by both species. In O. sativa, the same approach revealed 25 insertions in the 6 varietal groups. Four insertion sites located on chromosomes 1, 2, 10, and 11 were found orthologous in O. rufipogon and O. sativa. The chromosome 1 insertion was also shared between O. rufipogon and O. longistaminata. The presence of Tos17 at three insertion sites was confirmed by retrotransposon-based insertion polymorphism (RBIP) in a sample of O. sativa accessions. The first insertion, located on chromosome 3 was only found in two japonica accessions from the Bhutan region while the second insertion, located on chromosome 10 was specific to the varietal groups 1, 2, and 5. The third insertion located on chromosome 7 corresponds to the only insertion shown active in rice so far, notably in cv. Nipponbare, where it has been extensively used for insertion mutagenesis. This copy was only found in a few varieties of the japonica group 6 and in one group 5 accession. Taken together, these results confirm that Tos17 was probably present in the ancestor of A-genome species and that some copies of the element remained active in some Oryza lineages--notably in O. rufipogon and O. longistaminata--as well as in the indica and japonica O. sativa L. lineages.

Genetic variability in sunflower (Helianthus annuus L.) and in the Helianthus genus as assessed by retrotransposon-based molecular markers

The inter-retrotransposon amplified polymorphism (IRAP) protocol was applied for the first time within the genus Helianthus to assess intraspecific variability based on retrotransposon sequences among 36 wild accessions and 26 cultivars of Helianthus annuus L., and interspecific variability among 39 species of Helianthus. Two groups of LTRs, one belonging to a Copia-like retroelement and the other to a putative retrotransposon of unknown nature (SURE) have been isolated, sequenced and primers were designed to obtain IRAP fingerprints. The number of polymorphic bands in H. annuus wild accessions is as high as in Helianthus species. If we assume that a polymorphic band can be related to a retrotransposon insertion, this result suggests that retrotransposon activity continued after Helianthus speciation. Calculation of similarity indices from binary matrices (Shannon's and Jaccard's indices) show that variability is reduced among domesticated H. annuus. On the contrary, similarity indices among Helianthus species were as large as those observed among wild H. annuus accessions, probably related to their scattered geographic distribution. Principal component analysis of IRAP fingerprints allows the distinction between perennial and annual Helianthus species especially when the SURE element is concerned.

Comparative sequence analysis of MONOCULM1-orthologous regions in 14 Oryza genomes

Comparative genomics is a powerful tool to decipher gene and genome evolution. Placing multiple genome comparisons in a phylogenetic context improves the sensitivity of evolutionary inferences. In the genus Oryza, this comparative approach can be used to investigate gene function, genome evolution, domestication, polyploidy, and ecological adaptation. A large genomic region surrounding the MONOCULM1 (MOC1) locus was chosen for study in 14 Oryza species, including 10 diploids and 4 allotetraploids. Sequencing and annotation of 18 bacterial artificial chromosome clones for these species revealed highly conserved gene colinearity and structure in the MOC1 region. Since the Oryza radiation about 14 Mya, differences in transposon amplification appear to be responsible for the different current sizes of the Oryza genomes. In the MOC1 region, transposons were only conserved between genomes of the same type (e.g., AA or BB). In addition to the conserved gene content, several apparent genes have been generated de novo or uniquely retained in the AA lineage. Two different 3-gene segments have been inserted into the MOC1 region of O. coarctata (KK) or O. sativa by unknown mechanism(s). Large and apparently noncoding sequences flanking the MOC1 gene were observed to be under strong purifying selection. The allotetraploids Oryza alta and Oryza minuta were found to be products of recent polyploidization, less than 1.6 and 0.4 Mya, respectively. In allotetraploids, pseudogenization of duplicated genes was common, caused by large deletions, small frame-shifting insertions/deletions, or nonsense mutations.

The Tnt1 family member Retrosol copy number and structure disclose retrotransposon diversification in different Solanum species

Eukaryotic genome expansion/retraction caused by LTR-retrotransposon activity is dependent on the expression of full length copies to trigger efficient transposition and recombination-driven events. The Tnt1 family of retrotransposons has served as a model to evaluate the diversity among closely related elements within Solanaceae species and found that members of the family vary mainly in their U3 region of the long terminal repeats (LTRs). Recovery of a full length genomic copy of Retrosol was performed through a PCR-based approach from wild potato, Solanum oplocense. Further characterization focusing on both LTR sequences of the amplified copy allowed estimating an approximate insertion time at 2 million years ago thus supporting the occurrence of transposition cycles after genus divergence. Copy number of Tnt1-like elements in Solanum species were determined through genomic quantitative PCR whereby results sustain that Retrosol in Solanum species is a low copy number retrotransposon (1-4 copies) while Retrolyc1 has an intermediate copy number (38 copies) in S. peruvianum. Comparative analysis of retrotransposon content revealed no correlation between genome size or ploidy level and Retrosol copy number. The tetraploid cultivated potato with a cellular genome size of 1,715 Mbp harbours similar copy number per monoploid genome than other diploid Solanum species (613-884 Mbp). Conversely, S. peruvianum genome (1,125 Mbp) has a higher copy number. These results point towards a lineage specific dynamic flux regarding the history of amplification/activity of Tnt1-like elements in the genome of Solanum species.

Replication of nonautonomous retroelements in soybean appears to be both recent and common

Retrotransposons and their remnants often constitute more than 50% of higher plant genomes. Although extensively studied in monocot crops such as maize (Zea mays) and rice (Oryza sativa), the impact of retrotransposons on dicot crop genomes is not well documented. Here, we present an analysis of retrotransposons in soybean (Glycine max). Analysis of approximately 3.7 megabases (Mb) of genomic sequence, including 0.87 Mb of pericentromeric sequence, uncovered 45 intact long terminal repeat (LTR)-retrotransposons. The ratio of intact elements to solo LTRs was 8:1, one of the highest reported to date in plants, suggesting that removal of retrotransposons by homologous recombination between LTRs is occurring more slowly in soybean than in previously characterized plant species. Analysis of paired LTR sequences uncovered a low frequency of deletions relative to base substitutions, indicating that removal of retrotransposon sequences by illegitimate recombination is also operating more slowly. Significantly, we identified three subfamilies of nonautonomous elements that have replicated in the recent past, suggesting that retrotransposition can be catalyzed in trans by autonomous elements elsewhere in the genome. Analysis of 1.6 Mb of sequence from Glycine tomentella, a wild perennial relative of soybean, uncovered 23 intact retroelements, two of which had accumulated no mutations in their LTRs, indicating very recent insertion. A similar pattern was found in 0.94 Mb of sequence from Phaseolus vulgaris (common bean). Thus, autonomous and nonautonomous retrotransposons appear to be both abundant and active in Glycine and Phaseolus. The impact of nonautonomous retrotransposon replication on genome size appears to be much greater than previously appreciated.

Dynamic evolution of oryza genomes is revealed by comparative genomic analysis of a genus-wide vertical data set

Oryza (23 species; 10 genome types) contains the world's most important food crop - rice. Although the rice genome serves as an essential tool for biological research, little is known about the evolution of the other Oryza genome types. They contain a historical record of genomic changes that led to diversification of this genus around the world as well as an untapped reservoir of agriculturally important traits. To investigate the evolution of the collective Oryza genome, we sequenced and compared nine orthologous genomic regions encompassing the Adh1-Adh2 genes (from six diploid genome types) with the rice reference sequence. Our analysis revealed the architectural complexities and dynamic evolution of this region that have occurred over the past approximately 15 million years. Of the 46 intact genes and four pseudogenes in the japonica genome, 38 (76%) fell into eight multigene families. Analysis of the evolutionary history of each family revealed independent and lineage-specific gain and loss of gene family members as frequent causes of synteny disruption. Transposable elements were shown to mediate massive replacement of intergenic space (>95%), gene disruption, and gene/gene fragment movement. Three cases of long-range structural variation (inversions/deletions) spanning several hundred kilobases were identified that contributed significantly to genome diversification.

Natural selection on gene function drives the evolution of LTR retrotransposon families in the rice genome

Although the proliferation of LTR retrotransposons can cause major genomic modification and reorganization, the evolutionary dynamics that affect their frequency in host genomes are poorly understood. We analyzed patterns of genetic variation among LTR retrotransposons from Oryza sativa to investigate the type of selective forces that potentially limit their amplification and subsequent population of a nuclear genome. We performed both intra- and interfamily analyses of patterns of molecular sequence variation across multiple LTR retrotransposon genes. This analysis involved more than 1000 LTR retrotransposon sequences from 14 separate families that varied in both their insertion dates and full-length copy numbers. We uncovered evidence of strong purifying selection across all gene regions, but also indications that rare episodes of positive selection and adaptation to the host genome occur. Furthermore, our results indicate that LTR retrotransposons exhibit different but predictable patterns of sequence variation depending on their date of transposition, suggesting that LTR retrotransposons, regardless of superfamily and family classifications, show similar "life-histories."

Evidence of multiple horizontal transfers of the long terminal repeat retrotransposon RIRE1 within the genus Oryza

Horizontal gene transfer, defined as the transmission of genetic material between reproductively isolated species, has been considered for a long time to be a rare phenomenon. Most well-documented cases of horizontal gene transfer have been described in prokaryotes or in animals and they often involve transposable elements. The most abundant class of transposable elements in plant genomes are the long terminal repeat (LTR) retrotransposons. Because of their propensity to increase their copy number while active, LTR retrotransposons can have a significant impact on genomics changes during evolution. In a previous study, we showed that in the wild rice species Oryza australiensis, 60% of the genome is composed of only three families of LTR retrotransposons named RIRE1, Wallabi and Kangourou. In the present study, using both in silico and experimental approaches, we show that one of these three families, RIRE1, has been transferred horizontally between O. australiensis and seven other reproductively isolated Oryza species. This constitutes a new case of horizontal transfer in plants.

Construction, alignment and analysis of twelve framework physical maps that represent the ten genome types of the genus Oryza

Bacterial artificial chromosome (BAC) fingerprint and end-sequenced physical maps representing the ten genome types of Oryza are presented

We describe the establishment and analysis of a genus-wide comparative framework composed of 12 bacterial artificial chromosome fingerprint and end-sequenced physical maps representing the 10 genome types of Oryza aligned to the O. sativa ssp. japonica reference genome sequence. Over 932 Mb of end sequence was analyzed for repeats, simple sequence repeats, miRNA and single nucleotide variations, providing the most extensive analysis of Oryza sequence to date.