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

Characterization of transposable elements in the ectomycorrhizal fungus Laccaria bicolor

Background

The publicly available Laccaria bicolor genome sequence has provided a considerable genomic resource allowing systematic identification of transposable elements (TEs) in this symbiotic ectomycorrhizal fungus. Using a TE-specific annotation pipeline we have characterized and analyzed TEs in the L. bicolor S238N-H82 genome.

Methodology/Principal Findings

TEs occupy 24% of the 60 Mb L. bicolor genome and represent 25,787 full-length and partial copy elements distributed within 171 families. The most abundant elements were the Copia-like. TEs are not randomly distributed across the genome, but are tightly nested or clustered. The majority of TEs exhibits signs of ancient transposition except some intact copies of terminal inverted repeats (TIRS), long terminal repeats (LTRs) and a large retrotransposon derivative (LARD) element. There were three main periods of TE expansion in L. bicolor: the first from 57 to 10 Mya, the second from 5 to 1 Mya and the most recent from 0.5 Mya ago until now. LTR retrotransposons are closely related to retrotransposons found in another basidiomycete, Coprinopsis cinerea.

Conclusions

This analysis 1) represents an initial characterization of TEs in the L. bicolor genome, 2) contributes to improve genome annotation and a greater understanding of the role TEs played in genome organization and evolution and 3) provides a valuable resource for future research on the genome evolution within the Laccaria genus.

Recent retrotransposon insertions are methylated and phylogenetically clustered in japonica rice (Oryza sativa spp. japonica)

In plants, the genome of the host responds to the amplification of transposable elements (TEs) with DNA methylation. However, neither the factors involved in TE methylation nor the dynamics of the host-TE interaction are well resolved. Here, we identify 5,522 long terminal repeat retrotransposons (LTR-RT) in the genome of Oryza sativa ssp. japonica and then assess methylation for individual elements. Our analyses uncover three strong trends: long LTR-RTs are more highly methylated, the insertion times of LTR-RTs are negatively correlated with methylation, and young LTR-RTs tend to be closer to genes than older insertions. Additionally, a phylogenetic examination of the gypsy-like LTR-RT superfamily revealed that methylation is phylogenetically correlated. Given these observations, we present a model suggesting that the phylogenetic correlation among related LTR-RTs is a primary mechanism driving methylation. In this model, bursts of transposition produce new elements with high sequence similarity. The host machinery identifies proliferating elements as well as closely related LTR-RTs through cross-homology. In addition, our data are consistent with previous hypotheses that methylated LTR-RT elements are removed preferentially from regions near genes, explaining some of the observed age distribution.

The turbulent life of Sirevirus retrotransposons and the evolution of the maize genome: more than ten thousand elements tell the story

Sireviruses are one of the three genera of Copia long terminal repeat (LTR) retrotransposons, exclusive to and highly abundant in plants, and with a unique, among retrotransposons, genome structure. Yet, perhaps due to the few references to the Sirevirus origin of some families, compounded by the difficulty in correctly assigning retrotransposon families into genera, Sireviruses have hardly featured in recent research. As a result, analysis at this key level of classification and details of their colonization and impact on plant genomes are currently lacking. Recently, however, it became possible to accurately assign elements from diverse families to this genus in one step, based on highly conserved sequence motifs. Hence, Sirevirus dynamics in the relatively obese maize genome can now be comprehensively studied. Overall, we identified >10 600 intact and approximately 28 000 degenerate Sirevirus elements from a plethora of families, some brought into the genus for the first time. Sireviruses make up approximately 90% of the Copia population and it is the only genus that has successfully infiltrated the genome, possibly by experiencing intense amplification during the last 600 000 years, while being constantly recycled by host mechanisms. They accumulate in chromosome-distal gene-rich areas, where they insert in between gene islands, mainly in preferred zones within their own genomes. Sirevirus LTRs are heavily methylated, while there is evidence for a palindromic consensus target sequence. This work brings Sireviruses in the spotlight, elucidating their lifestyle and history, and suggesting their crucial role in the current genomic make-up of maize, and possibly other plant hosts.

Dynamics and adaptive benefits of protein domain emergence and arrangements during plant genome evolution

Plant genomes are generally very large, mostly paleopolyploid, and have numerous gene duplicates and complex genomic features such as repeats and transposable elements. Many of these features have been hypothesized to enable plants, which cannot easily escape environmental challenges, to rapidly adapt. Another mechanism, which has recently been well described as a major facilitator of rapid adaptation in bacteria, animals, and fungi but not yet for plants, is modular rearrangement of protein-coding genes. Due to the high precision of profile-based methods, rearrangements can be well captured at the protein level by characterizing the emergence, loss, and rearrangements of protein domains, their structural, functional, and evolutionary building blocks. Here, we study the dynamics of domain rearrangements and explore their adaptive benefit in 27 plant and 3 algal genomes. We use a phylogenomic approach by which we can explain the formation of 88% of all arrangements by single-step events, such as fusion, fission, and terminal loss of domains. We find many domains are lost along every lineage, but at least 500 domains are novel, that is, they are unique to green plants and emerged more or less recently. These novel domains duplicate and rearrange more readily within their genomes than ancient domains and are overproportionally involved in stress response and developmental innovations. Novel domains more often affect regulatory proteins and show a higher degree of structural disorder than ancient domains. Whereas a relatively large and well-conserved core set of single-domain proteins exists, long multi-domain arrangements tend to be species-specific. We find that duplicated genes are more often involved in rearrangements. Although fission events typically impact metabolic proteins, fusion events often create new signaling proteins essential for environmental sensing. Taken together, the high volatility of single domains and complex arrangements in plant genomes demonstrate the importance of modularity for environmental adaptability of plants.

On the Coevolution of Transposable Elements and Plant Genomes

Plant genomes are unique in an intriguing feature: the range of their size variation is unprecedented among living organisms. Although polyploidization contributes to this variability, transposable elements (TEs) seem to play the pivotal role. TEs, often considered intragenomic parasites, not only affect the genome size of the host, but also interact with other genes, disrupting and creating new functions and regulatory networks. Coevolution of plant genomes and TEs has led to tight regulation of TE activity, and growing evidence suggests their relationship became mutualistic. Although the expansions of TEs represent certain costs for the host genomes, they may also bring profits for populations, helping to overcome challenging environmental (biotic/abiotic stress) or genomic (hybridization and allopolyploidization) conditions. In this paper, we discuss the possibility that the possession of inducible TEs may provide a selective advantage for various plant populations.

Diversity of long terminal repeat retrotransposon genome distribution in natural populations of the wild diploid wheat Aegilops speltoides

The environment can have a decisive influence on the structure of the genome, changing it in a certain direction. Therefore, the genomic distribution of environmentally sensitive transposable elements may vary measurably across a species area. In the present research, we aimed to detect and evaluate the level of LTR retrotransposon intraspecific variability in Aegilops speltoides (2n = 2x = 14), a wild cross-pollinated relative of cultivated wheat. The interretrotransposon amplified polymorphism (IRAP) protocol was applied to detect and evaluate the level of retrotransposon intraspecific variability in Ae. speltoides and closely related species. IRAP analysis revealed significant diversity in TE distribution. Various genotypes from the 13 explored populations significantly differ with respect to the patterns of the four explored LTR retrotransposons (WIS2, Wilma, Daniela, and Fatima). This diversity points to a constant ongoing process of LTR retrotransposon fraction restructuring in populations of Ae. speltoides throughout the species' range and within single populations in time. Maximum changes were recorded in genotypes from small stressed populations. Principal component analysis showed that the dynamics of the Fatima element significantly differ from those of WIS2, Wilma, and Daniela. In terms of relationships between Sitopsis species, IRAP analysis revealed a grouping with Ae. sharonensis and Ae. longissima forming a separate unit, Ae. speltoides appearing as a dispersed group, and Ae. bicornis being in an intermediate position. IRAP display data revealed dynamic changes in LTR retrotransposon fractions in the genome of Ae. speltoides. The process is permanent and population specific, ultimately leading to the separation of small stressed populations from the main group.

Retrogenes in rice (Oryza sativa L. ssp. japonica) exhibit correlated expression with their source genes

Gene duplication occurs by either DNA- or RNA-based processes; the latter duplicates single genes via retroposition of messenger RNA. The expression of a retroposed gene copy (retrocopy) is expected to be uncorrelated with its source gene because upstream promoter regions are usually not part of the retroposition process. In contrast, DNA-based duplication often encompasses both the coding and the intergenic (promoter) regions; hence, expression is often correlated, at least initially, between DNA-based duplicates. In this study, we identified 150 retrocopies in rice (Oryza sativa L. ssp japonica), most of which represent ancient retroposition events. We measured their expression from high-throughput RNA sequencing (RNAseq) data generated from seven tissues. At least 66% of the retrocopies were expressed but at lower levels than their source genes. However, the tissue specificity of retrogenes was similar to their source genes, and expression between retrocopies and source genes was correlated across tissues. The level of correlation was similar between RNA- and DNA-based duplicates, and they decreased over time at statistically indistinguishable rates. We extended these observations to previously identified retrocopies in Arabidopsis thaliana, suggesting they may be general features of the process of retention of plant retrogenes.

Retrotransposon activation contributes to fragile X premutation rCGG-mediated neurodegeneration

Fragile X-associated tremor/ataxia syndrome (FXTAS) is a neurodegenerative disorder associated with fragile X premutation carriers. Previous studies have shown that fragile X rCGG repeats are sufficient to cause neurodegeneration and that the rCGG-repeat-binding proteins Pur α and heterogeneous nuclear ribonucleoprotein (hnRNP) A2/B1 could modulate rCGG-mediated neuronal toxicity. Mobile genetic elements or their remnants populate the genomes, and the activities of these elements are tightly controlled for the fitness of host genomes in different organisms. Here we provide both biochemical and genetic evidence to show that the activation of a specific retrotransposon, gypsy, can modulate rCGG-mediated neurodegeneration in an FXTAS Drosophila model. We find that one of the rCGG-repeat-binding proteins, hnRNP A2/B1, is involved in this process via interaction with heterochromatin protein 1. Knockdown of gypsy RNA by RNAi could suppress the neuronal toxicity caused by rCGG repeats. These data together point to a surprisingly active role for retrotransposition in neurodegeneration.

Temporal dynamics in the evolution of the sunflower genome as revealed by sequencing and annotation of three large genomic regions

Improved knowledge of genome composition, especially of its repetitive component, generates important informations in both theoretical and applied research. In this study, we provide the first insight into the local organization of the sunflower genome by sequencing and annotating 349,380 bp from 3 BAC clones, each including one single-copy gene. These analyses resulted in the identification of 11 putative gene sequences, 18 full-length LTR retrotransposons, 6 incomplete LTR retrotransposons, 2 non-autonomous LTR-retroelements (LINEs), 2 putative DNA transposons fragments and one putative helitron. Among LTR-retrotransposons, non-autonomous elements (the so-called LARDs), which do not carry any protein-encoding sequence, were discovered for the first time in the sunflower. The insertion time of intact retroelements was measured, based on sister LTRs divergence. All isolated elements were inserted relatively recently, especially those belonging to the Gypsy superfamily. Retrotransposon families related to those identified in the BAC clones are present also in other species of Helianthus, both annual and perennial, and even in other Asteraceae. In one of the three BAC clones, we found five copies of a lipid transfer protein (LTP) encoding gene within less than 100,000 bp, four of which are potentially functional. Two of these are interrupted by LTR retrotransposons, in the intron and in the coding sequence, respectively. The divergence between sister LTRs of the retrotransposons inserted within the genes indicates that LTP gene duplication started earlier than 1.749 MYRS ago. On the whole, the results reported in this study confirm that the sunflower is an excellent system to study transposons dynamics and evolution.

Polymorphisms and evolutionary history of retrotransposon insertions in rice promoters

Retrotransposons are ubiquitous in higher plant genomes. The presence or absence of retrotransposons in whole genome and high throughput genomic sequence (HTGS) from cultivated and wild rice was investigated to understand the organization and evolution of retrotransposon insertions in promoter regions. Approximately half of the Oryza sativa subsp. japonica ‘Nipponbare’ promoters with retrotransposons conserved in Oryza sativa subsp. indica ‘93-11’ and four wild rice species showed higher sequence conservation in retrotransposon than nonretrotransposon regions. We further investigated, in detail, the evolutionary dynamics of five retrotransposons in the promoter regions of 95 rice genotypes. Our data suggest that four of five insertions (Rp2–Rp5) occurred in the ancestor of AA genome, while the other insertion (Rp1) predates the ancestral divergence of Oryza officinalis (CC genome). Four retrotransposons (Rp2–Rp5) were present in 52% (Rp2), 29% (Rp3), 53% (Rp4), and 43% (Rp5) of the rice genotypes with AA genome type, and the fifth retrotransposon (Rp1) was present in 95% of the rice genotypes with AA, BBCC, or CC genome types. Furthermore, most of these retrotransposons were found to evolve slower than flanking promoter regions, suggesting a role in promoter function for regulating downstream genes.

Spliceosomal intron size expansion in domesticated grapevine (Vitis vinifera)

BACKGROUND

Spliceosomal introns are important components of eukaryotic genes as their structure, sizes and contents reflect the architecture of gene and genomes. Intron size, determined by both neutral evolution, repetitive elements activities and potential functional constraints, varies significantly in eukaryotes, suggesting unique dynamics and evolution in different lineages of eukaryotic organisms. However, the evolution of intron size, is rarely studied. To investigate intron size dynamics in flowering plants, in particular domesticated grapevines, a survey of intron size and content in wine grape (Vitis vinifera Pinot Noir) genes was conducted by assembling and mapping the transcriptome of V. vinifera genes from ESTs to characterize and analyze spliceosomal introns.

RESULTS

Uncommonly large size of spliceosomal intron was observed in V. vinifera genome, otherwise inconsistent with overall genome size dynamics when comparing Arabidopsis, Populus and Vitis. In domesticated grapevine, intron size is generally not related to gene function. The composition of enlarged introns in grapevines indicated extensive transposable element (TE) activity within intronic regions. TEs comprise about 80% of the expanded intron space and in particular, recent LTR retrotransposon insertions are enriched in these intronic regions, suggesting an intron size expansion in the lineage leading to domesticated grapevine, instead of size contractions in Arabidopsis and Populus. Comparative analysis of selected intronic regions in V. vinifera cultivars and wild grapevine species revealed that accelerated TE activity was associated with grapevine domestication, and in some cases with the development of specific cultivars.

CONCLUSIONS

In this study, we showed intron size expansion driven by TE activities in domesticated grapevines, likely a result of long-term vegetative propagation and intensive human care, which simultaneously promote TE proliferation and repress TE removal mechanisms such as recombination. The intron size expansion observed in domesticated grapevines provided an example of rapid plant genome evolution in response to artificial selection and propagation, and may shed light on the important genomic changes during domestication. In addition, the transcriptome approach used to gather intron size data significantly improved annotations of the V. vinifera genome.

Correlated evolution of LTR retrotransposons and genome size in the genus Eleocharis

BACKGROUND

Transposable elements (TEs) are considered to be an important source of genome size variation and genetic and phenotypic plasticity in eukaryotes. Most of our knowledge about TEs comes from large genomic projects and studies focused on model organisms. However, TE dynamics among related taxa from natural populations and the role of TEs at the species or supra-species level, where genome size and karyotype evolution are modulated in concert with polyploidy and chromosomal rearrangements, remain poorly understood. We focused on the holokinetic genus Eleocharis (Cyperaceae), which displays large variation in genome size and the occurrence of polyploidy and agmatoploidy/symploidy. We analyzed and quantified the long terminal repeat (LTR) retrotransposons Ty1-copia and Ty3-gypsy in relation to changes in both genome size and karyotype in Eleocharis. We also examined how this relationship is reflected in the phylogeny of Eleocharis.

RESULTS

Using flow cytometry, we measured the genome sizes of members of the genus Eleocharis (Cyperaceae). We found positive correlation between the independent phylogenetic contrasts of genome size and chromosome number in Eleocharis. We analyzed PCR-amplified sequences of various reverse transcriptases of the LTR retrotransposons Ty1-copia and Ty3-gypsy (762 sequences in total). Using real-time PCR and dot blot approaches, we quantified the densities of Ty1-copia and Ty3-gypsy within the genomes of the analyzed species. We detected an increasing density of Ty1-copia elements in evolutionarily younger Eleocharis species and found a positive correlation between Ty1-copia densities and C/n-values (an alternative measure of monoploid genome size) in the genus phylogeny. In addition, our analysis of Ty1-copia sequences identified a novel retrotransposon family named Helos1, which is responsible for the increasing density of Ty1-copia. The transition:transversion ratio of Helos1 sequences suggests that Helos1 recently transposed in later-diverging Eleocharis species.

CONCLUSIONS

Using several different approaches, we were able to distinguish between the roles of LTR retrotransposons, polyploidy and agmatoploidy/symploidy in shaping Eleocharis genomes and karyotypes. Our results confirm the occurrence of both polyploidy and agmatoploidy/symploidy in Eleocharis. Additionally, we introduce a new player in the process of genome evolution in holokinetic plants: LTR retrotransposons.

Stabilizing selection on genome size in a population of Festuca pallens under conditions of intensive intraspecific competition

*Stabilizing selection is a key evolutionary mechanism for which there is relatively little experimental evidence. To date, stabilizing selection has never been observed at the whole-genome level. *We tested the effect of selection on genome size in a field experiment using seeds collected in a population of Festuca pallens with a highly variable genome size. Using flow cytometry, we measured the genome size in germinating seedlings and juvenile plants grown with or without high intraspecific competition (908 individuals). Above-ground biomass and leaf number were used as measurements of individual vegetative performance. The possible confounding effect of seed weight was controlled for in a separate experiment. *Growth under high competition had a significant stabilizing effect on genome size. Because no relationship was observed between genome size and vegetative performance, we assume that the elimination of plants with extreme genome sizes was the result of decreased survival as a consequence of some unrecognized stress. *Our results indicate that genome size may be under direct selection. The equal disadvantaging of either large or small genomes indicates that the selection for optimum genome size in species may be fully context dependent. This study demonstrates the power of competition experiments for the detection of weak selection processes.

LTR retroelements in the genome of Daphnia pulex

BACKGROUND

Long terminal repeat (LTR) retroelements represent a successful group of transposable elements (TEs) that have played an important role in shaping the structure of many eukaryotic genomes. Here, we present a genome-wide analysis of LTR retroelements in Daphnia pulex, a cyclical parthenogen and the first crustacean for which the whole genomic sequence is available. In addition, we analyze transcriptional data and perform transposon display assays of lab-reared lineages and natural isolates to identify potential influences on TE mobility and differences in LTR retroelements loads among individuals reproducing with and without sex.

RESULTS

We conducted a comprehensive de novo search for LTR retroelements and identified 333 intact LTR retroelements representing 142 families in the D. pulex genome. While nearly half of the identified LTR retroelements belong to the gypsy group, we also found copia (95), BEL/Pao (66) and DIRS (19) retroelements. Phylogenetic analysis of reverse transcriptase sequences showed that LTR retroelements in the D. pulex genome form many lineages distinct from known families, suggesting that the majority are novel. Our investigation of transcriptional activity of LTR retroelements using tiling array data obtained from three different experimental conditions found that 71 LTR retroelements are actively transcribed. Transposon display assays of mutation-accumulation lines showed evidence for putative somatic insertions for two DIRS retroelement families. Losses of presumably heterozygous insertions were observed in lineages in which selfing occurred, but never in asexuals, highlighting the potential impact of reproductive mode on TE abundance and distribution over time. The same two families were also assayed across natural isolates (both cyclical parthenogens and obligate asexuals) and there were more retroelements in populations capable of reproducing sexually for one of the two families assayed.

CONCLUSIONS

Given the importance of LTR retroelements activity in the evolution of other genomes, this comprehensive survey provides insight into the potential impact of LTR retroelements on the genome of D. pulex, a cyclically parthenogenetic microcrustacean that has served as an ecological model for over a century.

Paleogenomic analysis of the short arm of chromosome 3 reveals the history of the African and Asian progenitors of cultivated rices

Rice is one of the most important crops, feeding more than half of the world population. There are two cultivated species, the African rice Oryza glaberrima and the Asian rice O. sativa. Although the African species is gradually replaced by O. sativa in most of African rice agrosystems, this species represents an important reservoir of genes of agronomical interest. Their exploitation for the development of modern African rice varieties requires a good understanding of the genetic relationships between the two cultivated species. We took advantage of the recent availability of the sequence of the chromosome 3 short arm of O. glaberrima to estimate the date of radiation between O. glaberrima and O. sativa lineages, using all the long terminal repeat (LTR)-retrotransposons as paleogenomic markers. We first demonstrated that in two distinct lineages, LTR-retrotransposons mutate at the same rate. Based on LTR-retrotransposons shared by both species in orthologous position, we then estimated that O. glaberrima and O. sativa progenitors diverged 1.2 Ma. This constitutes one of the first studies using such a large sample of transposable elements to reconstruct the phylogeny of species. Given the number of genome sequencing projects, there is no doubt that such approach will allow to resolve phylogenetic incongruities. The application of this method to other plant genomes will also facilitate further understanding of evolution of LTR-retrotransposons and eventually of the whole genome in divergent plant lineages.

Exceptional diversity, non-random distribution, and rapid evolution of retroelements in the B73 maize genome

Recent comprehensive sequence analysis of the maize genome now permits detailed discovery and description of all transposable elements (TEs) in this complex nuclear environment. Reiteratively optimized structural and homology criteria were used in the computer-assisted search for retroelements, TEs that transpose by reverse transcription of an RNA intermediate, with the final results verified by manual inspection. Retroelements were found to occupy the majority (>75%) of the nuclear genome in maize inbred B73. Unprecedented genetic diversity was discovered in the long terminal repeat (LTR) retrotransposon class of retroelements, with >400 families (>350 newly discovered) contributing >31,000 intact elements. The two other classes of retroelements, SINEs (four families) and LINEs (at least 30 families), were observed to contribute 1,991 and ∼35,000 copies, respectively, or a combined ∼1% of the B73 nuclear genome. With regard to fully intact elements, median copy numbers for all retroelement families in maize was 2 because >250 LTR retrotransposon families contained only one or two intact members that could be detected in the B73 draft sequence. The majority, perhaps all, of the investigated retroelement families exhibited non-random dispersal across the maize genome, with LINEs, SINEs, and many low-copy-number LTR retrotransposons exhibiting a bias for accumulation in gene-rich regions. In contrast, most (but not all) medium- and high-copy-number LTR retrotransposons were found to preferentially accumulate in gene-poor regions like pericentromeric heterochromatin, while a few high-copy-number families exhibited the opposite bias. Regions of the genome with the highest LTR retrotransposon density contained the lowest LTR retrotransposon diversity. These results indicate that the maize genome provides a great number of different niches for the survival and procreation of a great variety of retroelements that have evolved to differentially occupy and exploit this genomic diversity.

Although TEs are a major component of all studied plant genomes, and are the most significant contributors to genome structure and evolution in almost all eukaryotes that have been investigated, their properties and reasons for existence are not well understood in any eukaryotic genome. In order to begin a comprehensive study of TE contributions to the structure, function, and evolution of both genes and genomes, we first identified all of the TEs in maize and then investigated whether there were non-random patterns in their dispersal. We used homology and TE structure criteria in an effort to discover all of the retroelements in the recently sequenced genome from maize inbred B73. We found that the retroelements are incredibly diverse in maize, with many hundreds of families that show different insertion and/or retention specificities across the maize chromosomes. Most of these element families are present in low copy numbers and had been missed by previous searches that relied on a high-copy-number criterion. Different element families exhibited very different biases for accumulation across the chromosomes, indicating that they can detect and utilize many different chromatin environments.