Structural and evolutionary analysis of the copia-like elements in the Arabidopsis thaliana genome

DNA Demethylation in Response to Heat Stress in Arabidopsis thaliana

Environmental stress is one of the most important factors affecting plant growth and development. Recent studies have shown that epigenetic mechanisms, such as DNA methylation, play a key role in adapting plants to stress conditions. Here, we analyzed the dynamics of changes in the level of DNA methylation in Arabidopsis thaliana (L.) Heynh. (Brassicaceae) under the influence of heat stress. For this purpose, whole-genome sequencing of sodium bisulfite-treated DNA was performed. The analysis was performed at seven time points, taking into account the control conditions, heat stress, and recovery to control conditions after the stress treatment was discontinued. In our study we observed decrease in the level of DNA methylation under the influence of heat stress, especially after returning to control conditions. Analysis of the gene ontology enrichment and regulatory pathways showed that genes characterized by differential DNA methylation are mainly associated with stress response, including heat stress. These are the genes encoding heat shock proteins and genes associated with translation regulation. A decrease in the level of DNA methylation in such specific sites suggests that under the influence of heat stress we observe active demethylation phenomenon rather than passive demethylation, which is not locus specific.

Genomic localization of AtRE1 and AtRE2, copia-type retrotransposons, in natural variants of Arabidopsis thaliana

Retrotransposons are ubiquitous components of plant genomes. They affect genome organization, and can also affect the expression patterns of neighboring genes. Retrotransposons are therefore important elements for changing genomic information. To understand the evolution of the Arabidopsis genome, we examined the distribution of certain retrotransposons, AtRE1s and AtRE2s, in the genomes of 12 natural variants (accessions) of Arabidopsis thaliana. AtRE1 and AtRE2 are copia-type retrotransposons that are potentially active. Their copy numbers are low, and they are absent from the genomes of some accessions. We detected four loci with AtRE1s inserted in six accessions, and one locus with an insertion of a solo-LTR-like sequence derived from AtRE1 in two accessions. Seven loci with AtRE2s inserted were detected on eight accessions. These loci were distributed in euchromatic regions of chromosomes 1, 2, 3, and 4. The AtRE1 and AtRE2 sequences at some loci identified in this study have not been recorded in the database of the 1001 Genome project. The sequences of AtRE1s and those of AtRE2s in different accessions and at different loci were highly conserved. There was a complete or almost complete conservation of sequences of both long terminal repeats in each AtRE1 and in each AtRE2. These results suggest that AtRE1 and AtRE2 appeared quite recently in the Arabidopsis genome. Furthermore, sequence comparisons of AtRE1 and AtRE2 loci among accessions revealed the possibility that large deletions containing entire sequences of AtRE1 and AtRE2 have occurred in some accessions.

Genetics of cryptic speciation within an Arctic mustard, Draba nivalis

Crossing experiments indicate that hybrid sterility barriers frequently have developed within diploid, circumpolar plant species of the genus Draba. To gain insight into the rapid evolution of postzygotic reproductive isolation in this system, we augmented the linkage map of one of these species, D. nivalis, and searched for quantitative trait loci (QTLs) associated with reproductive isolation. The map adds 63 new dominant markers to a previously published dataset of 31 co-dominant microsatellites. These markers include 52 amplified fragment length polymorphisms (AFLPs) and 11 sequence-specific amplified polymorphisms (SSAPs) based on retrotransposon sequence. 22 markers displaying transmission ratio distortion were further included in the map. We resolved eight linkage groups with a total map length of 894 cM. Significant genotype-trait associations, or quantitative trait loci (QTL), were detected for reproductive phenotypes including pollen fertility (4 QTLs), seed set (3 QTLs), flowering time (3 QTLs) and number of flowers (4 QTLs). Observed patterns of inheritance were consistent with the influence of both nuclear-nuclear interactions and chromosomal changes on these traits. All seed set QTLs and one pollen fertility QTL displayed underdominant effects suggestive of the involvement of chromosomal rearrangements in hybrid sterility. Interestingly, D. nivalis is predominantly self-fertilizing, which may facilitate the establishment of underdominant loci and contribute to reproductive isolation.

The structure, organization and radiation of Sadhu non-long terminal repeat retroelements in Arabidopsis species

Background

Sadhu elements are non-autonomous retroposons first recognized in Arabidopsis thaliana. There is a wide degree of divergence among different elements, suggesting that these sequences are ancient in origin. Here we report the results of several lines of investigation into the genomic organization and evolutionary history of this element family.

Results

We present a classification scheme for Sadhu elements in A. thaliana, describing derivative elements related to the full-length elements we reported previously. We characterized Sadhu5 elements in a set of A. thaliana strains in order to trace the history of radiation in this subfamily. Sequences surrounding the target sites of different Sadhu insertions are consistent with mobilization by LINE retroelements. Finally, we identified Sadhu elements grouping into distinct subfamilies in two related species, Arabidopsis arenosa and Arabidopsis lyrata.

Conclusions

Our analyses suggest that the Sadhu retroelement family has undergone target primed reverse transcription-driven retrotransposition during the divergence of different A. thaliana strains. In addition, Sadhu elements can be found at moderate copy number in three distinct Arabidopsis species, indicating that the evolutionary history of these sequences can be traced back at least several millions of years.

The marker SCK13(603) associated with resistance to ascochyta blight in chickpea is located in a region of a putative retrotransposon

The sequence characterized amplified region (SCAR) marker SCK13(603), associated with ascochyta blight resistance in a chickpea recombinant inbred line (RIL) population, was used as anchored sequence for genome walking. The PCRs performed in the walking steps to walk in the same direction produced eight bands in 5' direction and five bands in 3' direction with a length ranking from 530 to 2,871 bp. The assembly of the bands sequences along with the sequence of SCK13(603) resulted in 7,815 bp contig. Blastn analyses showed stretches of DNA sequence mainly distributed from the nucleotides 1,500 to 4,500 significantly similar to Medicago truncatula genomic DNA. Three open reading frames (ORFs) were identified and blastp analysis of predicted amino acids sequences revealed that ORF1, ORF2 and ORF3 had significant similarity to a CCHC zinc finger protein, to an integrase, and to a precursor of the glucoamylase s1/s2, respectively, from M. truncatula. The high homology of the putative proteins derived from ORF1 and ORF2 with retrotransposon proteins and the prediction of the existence of conserved domains usually present in retrotransposon proteins indicate that the marker SCK13(603) is located in a region of a putative retrotransposon. The information generated in this study has contributed to increase the knowledge of this important region for blight resistance in chickpea.

Novel clades of chromodomain-containing Gypsy LTR retrotransposons from mosses (Bryophyta)

Retrotransposons are the major component of plant genomes. Chromodomain-containing Gypsy long terminal repeat (LTR) retrotransposons are widely distributed in eukaryotes. Four distinct clades of chromodomain-containing Gypsy retroelements are known from the vascular plants: Reina, CRM, Galadriel and Tekay. At the same time, almost nothing is known about the repertoire of LTR retrotransposons in bryophyte genomes. We have combined a search of chromodomain-containing Gypsy retroelements in Physcomitrella genomic sequences and an experimental investigation of diverse moss species. The computer-based mining of the chromodomain-containing LTR retrotransposons allowed us to describe four different elements from Physcomitrella. Four novel clades were identified that are evolutionarily distinct from the chromodomain-containing Gypsy LTR retrotransposons of other plants.

AtCopeg1, the unique gene originated from AtCopia95 retrotransposon family, is sensitive to external hormones and abiotic stresses

Retrotransposons, the important component of eukaryotic genome, are seeds of evolution and play great role in creating new genes. The compact Arabidopsis genome harbors over 200 Copia-like retrotransposons, but mostly silent. Here we isolated an expressed gene AtCopeg1 (Copia evolved gene 1), which shows higher than 90% identity to AtCopia95_I, the consensus sequence encoding AtCopia95 polyprotein. AtCopeg1 is the unique gene evolved from AtCopia95 family. It is an intron-containing gene with two alternative 3' ends. The transcript accumulation of AtCopeg1 is tissue-specific, also significantly affected by external hormones and abiotic stresses. The presence of regulatory elements in its promoter region (originating from AtCopia95_I and AtCopia95 long terminal repeat), is adequate for conferring its essential expression feature. Thus, AtCopeg1 is a versatile functional gene involved in many developmental and adaptive processes probably including the signaling crosstalk of hormone and nutrient stress. Our work highlighted the role of transposable elements in creating new functional genes, and will incite the enthusiasm for isolation and functional characterization of plant genes evolved from those previously considered as selfish and junk DNA.

LTR retrotransposons in rice (Oryza sativa, L.): recent burst amplifications followed by rapid DNA loss

Background

LTR retrotransposons are one of the main causes for plant genome size and structure evolution, along with polyploidy. The characterization of their amplification and subsequent elimination of the genomes is therefore a major goal in plant evolutionary genomics. To address the extent and timing of these forces, we performed a detailed analysis of 41 LTR retrotransposon families in rice.

Results

Using a new method to estimate the insertion date of both truncated and complete copies, we estimated these two forces more accurately than previous studies based on other methods. We show that LTR retrotransposons have undergone bursts of amplification within the past 5 My. These bursts vary both in date and copy number among families, revealing that each family has a particular amplification history. The number of solo LTR varies among families and seems to correlate with LTR size, suggesting that solo LTR formation is a family-dependent process. The deletion rate estimate leads to the prediction that the half-life of LTR retrotransposon sequences evolving neutrally is about 19 My in rice, suggesting that other processes than the formation of small deletions are prevalent in rice DNA removal.

Conclusion

Our work provides insights into the dynamics of LTR retrotransposons in the rice genome. We show that transposable element families have distinct amplification patterns, and that the turn-over of LTR retrotransposons sequences is rapid in the rice genome.

Comparative cross-species alternative splicing in plants

Alternative splicing (AS) can add significantly to genome complexity. Plants are thought to exhibit less AS than animals. An algorithm, based on expressed sequence tag (EST) pairs gapped alignment, was developed that takes advantage of the relatively small intron and exon size in plants and directly compares pairs of ESTs to search for AS. EST pairs gapped alignment was first evaluated in Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), and tomato (Solanum lycopersicum) for which annotated genome sequence is available and was shown to accurately predict splicing events. The method was then applied to 11 plant species that include 17 cultivars for which enough ESTs are available. The results show a large, 3.7-fold difference in AS rates between plant species with Arabidopsis and rice in the lower range and lettuce (Lactuca sativa) and sorghum (Sorghum bicolor) in the upper range. Hence, compared to higher animals, plants show a much greater degree of variety in their AS rates and in some plant species the rates of animal and plant AS are comparable although the distribution of AS types may differ. In eudicots but not monocots, a correlation between genome size and AS rates was detected, implying that in eudicots the mechanisms that lead to larger genomes are a driving force for the evolution of AS.

The genomic organization of retrotransposons in Brassica oleracea

We have investigated the copy numbers and genomic organization of five representative reverse transcriptase domains from retrotransposons in Brassica oleracea. Two non-homologous Pseudoviridae (Ty1/copia-like) elements, two Metaviridae (Ty3/gypsy-like) elements (one related to the Athila family) and one Retroposinae (LINE) element were hybridized to a gridded BAC library, "BoB". The results indicated that the individual LTR retrotransposons (copia and gypsy-like) were represented by between 90 and 320 copies in the haploid genome, with only evidence of a single location for the LINE. Sequence analysis of the same elements against genome survey sequence gave estimates of between 60 and 570, but no LINE was found. There was minimal evidence for clustering between any of these retroelements: only half the randomly expected number of BACs hybridized to both LTR-retrotransposon families. Fluorescent in situ hybridization showed that each of the retroelements had a characteristic genomic distribution. Our results suggest there are preferential sites and perhaps control mechanisms for the insertion or excision of different retrotransposon groups.

LTR retrotransposons and flowering plant genome size: emergence of the increase/decrease model

Long Terminal Repeat (LTR) retrotransposons are ubiquitous components of plant genomes. Because of their copy-and-paste mode of transposition, these elements tend to increase their copy number while they are active. In addition, it is now well established that the differences in genome size observed in the plant kingdom are accompanied by variations in LTR retrotransposon content, suggesting that LTR retrotransposons might be important players in the evolution of plant genome size, along with polyploidy. The recent availability of large genomic sequences for many crop species has made it possible to examine in detail how LTR retrotransposons actually drive genomic changes in plants. In the present paper, we provide a review of the recent publications that have contributed to the knowledge of plant LTR retrotransposons, as structural components of the genomes, as well as from an evolutionary genomic perspective. These studies have shown that plant genomes undergo genome size increases through bursts of retrotransposition, while there is a counteracting process that tends to eliminate the transposed copies from the genomes. This process involves recombination mechanisms that occur either between the LTRs of the elements, leading to the formation of solo-LTRs, or between direct repeats anywhere in the sequence of the element, leading to internal deletions. All these studies have led to the emergence of a new model for plant genome evolution that takes into account both genome size increases (through retrotransposition) and decreases (through solo-LTR and deletion formation). In the conclusion, we discuss this new model and present the future prospects in the study of plant genome evolution in relation to the activity of transposable elements.

Genomic neighborhoods for Arabidopsis retrotransposons: a role for targeted integration in the distribution of the Metaviridae

BACKGROUND

Retrotransposons are an abundant component of eukaryotic genomes. The high quality of the Arabidopsis thaliana genome sequence makes it possible to comprehensively characterize retroelement populations and explore factors that contribute to their genomic distribution.

RESULTS

We identified the full complement of A. thaliana long terminal repeat (LTR) retroelements using RetroMap, a software tool that iteratively searches genome sequences for reverse transcriptases and then defines retroelement insertions. Relative ages of full-length elements were estimated by assessing sequence divergence between LTRs: the Pseudoviridae were significantly younger than the Metaviridae. All retroelement insertions were mapped onto the genome sequence and their distribution was distinctly non-uniform. Although both Pseudoviridae and Metaviridae tend to cluster within pericentromeric heterochromatin, this association is significantly more pronounced for all three Metaviridae sublineages (Metavirus, Tat and Athila). Among these, Tat and Athila are strictly associated with pericentromeric heterochromatin.

CONCLUSIONS

The non-uniform genomic distribution of the Pseudoviridae and the Metaviridae can be explained by a variety of factors including target-site bias, selection against integration into euchromatin and pericentromeric accumulation of elements as a result of suppression of recombination. However, comparisons based on the age of elements and their chromosomal location indicate that integration-site specificity is likely to be the primary factor determining distribution of the Athila and Tat sublineages of the Metaviridae. We predict that, like retroelements in yeast, the Athila and Tat elements target integration to pericentromeric regions by recognizing a specific feature of pericentromeric heterochromatin.

Genome-wide comparative analysis of the transposable elements in the related species Arabidopsis thaliana and Brassica oleracea

Transposable elements (TEs) are the major component of plant genomes where they contribute significantly to the >1,000-fold genome size variation. To understand the dynamics of TE-mediated genome expansion, we have undertaken a comparative analysis of the TEs in two related organisms: the weed Arabidopsis thaliana (125 megabases) and Brassica oleracea ( approximately 600 megabases), a species with many crop plants. Comparison of the whole genome sequence of A. thaliana with a partial draft of B. oleracea has permitted an estimation of the patterns of TE amplification, diversification, and loss that has occurred in related species since their divergence from a common ancestor. Although we find that nearly all TE lineages are shared, the number of elements in each lineage is almost always greater in B. oleracea. Class 1 (retro) elements are the most abundant TE class in both species with LTR and non-LTR elements comprising the largest fraction of each genome. However, several families of class 2 (DNA) elements have amplified to very high copy number in B. oleracea where they have contributed significantly to genome expansion. Taken together, the results of this analysis indicate that amplification of both class 1 and class 2 TEs is responsible, in part, for B. oleracea genome expansion since divergence from a common ancestor with A. thaliana. In addition, the observation that B. oleracea and A. thaliana share virtually all TE lineages makes it unlikely that wholesale removal of TEs is responsible for the compact genome of A. thaliana.

The diversity of retroelements in diploid and allotetraploid Brassica species

Using universal PCR primers, some 80 fragments of retroelement reverse transcriptase genes were isolated from 16 accessions of the three diploid and three derived allotetraploid species of Brassica in the triangle of U. Sequence analysis showed that the Ty1/copia and LINE-like elements were distinct, while a third clade could be sub-divided into Ty3/gypsy, Athila and virus-like branches, providing evidence that there are multiple sub-lineages within this group normally considered to be gypsy-like elements in plants. The parsimony trees showed no branches correlating with the known genome relationships for the six diploid and allotetraploid Brassica species, probably because members of the element families were present in the common ancestor of the Brassica and, unlike other repetitive sequences, there is no evidence for genome-wide homogenization, although convergent evolution or horizontal transfer cannot be ruled out. Southern hybridization suggested some sub-families were amplified in individual species. The data show that retroelement sequence data do not allow inference of phylogeny, but knowledge of evolution of such abundant sequences assists in exploitation and interpretation of data from other species including models with much smaller genomes and may provide markers.

Purcopia, a Ty1-copia truncated retroelement in the genome of Claviceps purpurea

Truncated copy of reverse transcriptase of Ty1/copia retroelement (Purcopia) was found as part of the species-specific RAPD 257(540) marker of Claviceps purpurea. A region of 94 bp with 78.9% identity to an unannotated region of the genomic clone of the rice blast fungus Pyricularia grisea (accession no. AQ162050) was found at the 5' end of the pseudogene. Comparison with database sequences revealed that Purcopia is close to the plant retroelements represented by Tto1, Ta1-3 and Bare-1, whereas the other fungal elements of the Ty1/copia type grouped with Hopscotch elements. Restriction patterns obtained by hybridization of the labeled marker to HindIII digested genomic DNA of various C. purpurea isolates contained multiple bands. The banding was individual and did not yield any species- or population-specific fragments or patterns.

Genes of the Pseudoviridae (Ty1/copia retrotransposons)

A comprehensive survey of the Pseudoviridae (Ty1/copia) retroelement family was conducted using the GenBank sequence database and completed genome sequences of several model organisms. Plant genomes were the most abundant sources of Pseudoviridae, with the Arabidopsis thaliana genome having 276 distinct elements. A reverse transcriptase amino acid sequence phylogeny indicated that the Pseudoviridae comprises highly divergent members. Coding sequences for a representative subset of elements were analyzed to identify conserved domains and differences that may underlie functional divergence. With the exception of some fungal elements (e.g., Ty1), most Pseudoviridae encode Gag and Pol on a single open reading frame. In addition to the nearly ubiquitous RNA-binding motif of nucleocapsid, three new conserved domains were identified in Gag. pol-encoded aspartic protease was similar to the retroviral enzyme and could be mapped onto the HIV-1 structure. Pol was highly conserved throughout the family. The greatest divergence among Pol sequences was seen in the C-terminus of integrase (IN). We defined a large motif (GKGY) after the IN catalytic domain that is unique to the Pseudoviridae. Additionally, the extreme C-terminus of IN is rich in simple sequence motifs. A distinct lineage of Pseudoviridae in plants have envlike genes. This lineage has undergone a large expansion of Gag characterized by an alpha-helix-rich domain containing coiled-coil motifs. In several elements, this domain is flanked on both sides by RNA-binding domains. We propose that this monophyletic lineage defines a new Pseudoviridae genus, herein referred to as the AGROVIRUS:

Genome size reduction through illegitimate recombination counteracts genome expansion in Arabidopsis

Genome size varies greatly across angiosperms. It is well documented that, in addition to polyploidization, retrotransposon amplification has been a major cause of genome expansion. The lack of evidence for counterbalancing mechanisms that curtail unlimited genome growth has made many of us wonder whether angiosperms have a "one-way ticket to genomic obesity." We have therefore investigated an angiosperm with a well-characterized and notably small genome, Arabidopsis thaliana, for evidence of genomic DNA loss. Our results indicate that illegitimate recombination is the driving force behind genome size decrease in Arabidopsis, removing at least fivefold more DNA than unequal homologous recombination. The presence of highly degraded retroelements also suggests that retrotransposon amplification has not been confined to the last 4 million years, as is indicated by the dating of intact retroelements.