Transposon diversity in Arabidopsis thaliana

The Genomic Shock Hypothesis: Genetic and Epigenetic Alterations of Transposable Elements after Interspecific Hybridization in Plants

Transposable elements (TEs) are major components of plant genomes with the ability to change their position in the genome or to create new copies of themselves in other positions in the genome. These can cause gene disruption and large-scale genomic alterations, including inversions, deletions, and duplications. Host organisms have evolved a set of mechanisms to suppress TE activity and counter the threat that they pose to genome integrity. These includes the epigenetic silencing of TEs mediated by a process of RNA-directed DNA methylation (RdDM). In most cases, the silencing machinery is very efficient for the vast majority of TEs. However, there are specific circumstances in which TEs can evade such silencing mechanisms, for example, a variety of biotic and abiotic stresses or in vitro culture. Hybridization is also proposed as an inductor of TE proliferation. In fact, the discoverer of the transposons, Barbara McClintock, first hypothesized that interspecific hybridization provides a "genomic shock" that inhibits the TE control mechanisms leading to the mobilization of TEs. However, the studies carried out on this topic have yielded diverse results, showing in some cases a total absence of mobilization or being limited to only some TE families. Here, we review the current knowledge about the impact of interspecific hybridization on TEs in plants and the possible implications of changes in the epigenetic mechanisms.

Virus-like transposons cross the species barrier and drive the evolution of genetic incompatibilities

Horizontal gene transfer, the movement of genetic material between species, has been reported across all major eukaryotic lineages. However, the underlying mechanisms of transfer and their impact on genome evolution are still poorly understood. While studying the evolutionary origin of a selfish element in the nematode Caenorhabditis briggsae, we discovered that Mavericks, ancient virus-like transposons related to giant viruses and virophages, are one of the long-sought vectors of horizontal gene transfer. We found that Mavericks gained a novel herpesvirus-like fusogen in nematodes, leading to the widespread exchange of cargo genes between extremely divergent species, bypassing sexual and genetic barriers spanning hundreds of millions of years. Our results show how the union between viruses and transposons causes horizontal gene transfer and ultimately genetic incompatibilities in natural populations.

Regulatory mechanism of a heat-activated retrotransposon by DDR complex in Arabidopsis thaliana

The RNA-directed DNA methylation (RdDM) pathway plays an essential role in the transposon silencing mechanism; the DDR complex, consisting of DRD1, DMS3, and RDM1, is an essential component of the RdDM pathway. ONSEN, identified in Arabidopsis, is a retrotransposon activated by heat stress at 37°C; however, studies on the regulation of ONSEN are limited. In this study, we analyzed the regulation of ONSEN activity by the DDR complex in Arabidopsis. We elucidated that loss of any component of the DDR complex increased ONSEN transcript levels. Transgenerational transposition of ONSEN was observed in the DDR-complex mutants treated with heat stress for 48 h. Furthermore, the DDR complex components DRD1, DMS3, and RDM1 played independent roles in suppressing ONSEN transcription and transposition. Moreover, we found that the duration of heat stress affects ONSEN activity. Therefore, the results of this study provide new insights into the retrotransposon regulatory mechanisms of the DDR complex in the RdDM pathway.

GmRAV confers ecological adaptation through photoperiod control of flowering time and maturity in soybean

Photoperiod strictly controls vegetative and reproductive growth stages in soybean (Glycine max). A soybean GmRAV (Related to ABI3/VP1) transcription factor containing both AP2 and B3 domains was shown to be a key component of this process. We identified six polymorphisms in the GmRAV promoter that showed significant association with flowering time and maturity of soybean in one or multiple environments. Soybean varieties with minor polymorphism exhibited a longer growth period contributing to soybean adaptation to lower latitudes. The cis-acting element GT1CONSENSUS motif of the GmRAV promoter controlled the growth period, and the major allele in this motif shortened duration of late reproductive stages by reducing GmRAV expression levels. Three GmRAV-overexpressing (GmRAV-ox) transgenic lines displayed later flowering time and maturity, shorter height and fewer numbers of leaves compared with control plants, whereas transgenic inhibition of GmRAV expression resulted in earlier flowering time and maturity and increased plant height. Combining DNA affinity purification sequencing and RNA sequencing analyses revealed 154 putative target genes directly bound and transcriptionally regulated by GmRAV. Two GmRAV binding motifs [C(A/G)AACAA(G/T)A(C/T)A(G/T)] and [C(T/A)A(C)C(T/G)CTG] were identified, and acting downstream of E3E4, GmRAV repressed GmFT5a transcriptional activity through binding a CAACA motif, thereby delaying soybean growth and extending both vegetative and reproductive phases.

An Overview of Duplicated Gene Detection Methods: Why the Duplication Mechanism Has to Be Accounted for in Their Choice

Gene duplication is an important evolutionary mechanism allowing to provide new genetic material and thus opportunities to acquire new gene functions for an organism, with major implications such as speciation events. Various processes are known to allow a gene to be duplicated and different models explain how duplicated genes can be maintained in genomes. Due to their particular importance, the identification of duplicated genes is essential when studying genome evolution but it can still be a challenge due to the various fates duplicated genes can encounter. In this review, we first describe the evolutionary processes allowing the formation of duplicated genes but also describe the various bioinformatic approaches that can be used to identify them in genome sequences. Indeed, these bioinformatic approaches differ according to the underlying duplication mechanism. Hence, understanding the specificity of the duplicated genes of interest is a great asset for tool selection and should be taken into account when exploring a biological question.

Prediction of regulatory long intergenic non-coding RNAs acting in trans through base-pairing interactions

Background

Long intergenic non-coding RNAs (lincRNAs) can act as regulators of expression of protein-coding genes. Trans-natural antisense transcripts (trans-NATs) are a type of lincRNAs that contain sequence complementary to mRNA from other loci. The regulatory potential of trans-NATs has been poorly studied in eukaryotes and no example of trans-NATs regulating gene expression in plants are reported. The goal of this study was to identify lincRNAs, and particularly trans-NATs, in Arabidopsis thaliana that have a potential to regulate expression of target genes in trans at the transcriptional or translational level.

Results

We identified 1001 lincRNAs using an RNAseq dataset from total polyA⁺ and polysome-associated RNA of seedlings grown under high and low phosphate, or shoots and roots treated with different phytohormones, of which 550 were differentially regulated. Approximately 30% of lincRNAs showed conservation amongst Brassicaceae and 25% harbored transposon element (TE) sequences. Gene co-expression network analysis highlighted a group of lincRNAs associated with the response of roots to low phosphate. A total of 129 trans-NATs were predicted, of which 88 were significantly differentially expressed under at least one pairwise comparison. Five trans-NATs showed a positive correlation between their expression and target mRNA steady-state levels, and three showed a negative correlation. Expression of four trans-NATs positively correlated with a change in target mRNA polysome association. The regulatory potential of these trans-NATs did not implicate miRNA mimics nor siRNAs. We also looked for lincRNAs that could regulate gene expression in trans by Watson-Crick DNA:RNA base pairing with target protein-encoding loci. We identified 100 and 81 with a positive or negative correlation, respectively, with steady-state level of their predicted target. The regulatory potential of one such candidate lincRNA harboring a SINE TE sequence was validated in a protoplast assay on three distinct genes containing homologous TE sequence in their promoters. Construction of networks highlighted other putative lincRNAs with multiple predicted target loci for which expression was positively correlated with target gene expression.

Conclusions

This study identified lincRNAs in Arabidopsis with potential in regulating target gene expression in trans by both RNA:RNA and RNA:DNA base pairing and highlights lincRNAs harboring TE sequences in such activity.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5946-0) contains supplementary material, which is available to authorized users.

Helitron distribution in Brassicaceae and whole Genome Helitron density as a character for distinguishing plant species

BACKGROUND

Helitron is a rolling-circle DNA transposon; it plays an important role in plant evolution. However, Helitron distribution and contribution to evolution at the family level have not been previously investigated.

RESULTS

We developed the software easy-to-annotate Helitron (EAHelitron), a Unix-like command line, and used it to identify Helitrons in a wide range of 53 plant genomes (including 13 Brassicaceae species). We determined Helitron density (abundance/Mb) and visualized and examined Helitron distribution patterns. We identified more than 104,653 Helitrons, including many new Helitrons not predicted by other software. Whole genome Helitron density is independent from genome size and shows stability at the species level. Using linear discriminant analysis, de novo genomes (next-generation sequencing) were successfully classified into Arabidopsis thaliana groups. For most Brassicaceae species, Helitron density negatively correlated with gene density, and Helitron distribution patterns were similar to those of A. thaliana. They preferentially inserted into sequence around the centromere and intergenic region. We also associated 13 Helitron polymorphism loci with flowering-time phenotypes in 18 A. thaliana ecotypes.

CONCLUSION

EAHelitron is a fast and efficient tool to identify new Helitrons. Whole genome Helitron density can be an informative character for plant classification. Helitron insertion polymorphism could be used in association analysis.

Evolution of Mutator transposable elements across eukaryotic diversity

Background

Mutator-like elements (MULEs) are a significant superfamily of DNA transposons on account of their: (i) great transpositional activity and propensity for insertion in or near gene sequences, (ii) their consequent high mutagenic capacity, and, (iii) their tendency to acquire host gene fragments. Consequently, MULEs are important genetic tools and represent a key study system for research into host-transposon interactions. Yet, while several studies have focused on the impacts of MULEs on crop and fungus genomes, their evolution remains poorly explored.

Results

We perform comprehensive bioinformatic and phylogenetic analyses to address currently available MULE diversity and reconstruct evolution for the group. For this, we mine MULEs from online databases, and combine search results with available transposase sequences retrieved from previously published studies. Our analyses uncover two entirely new MULE clades that contain elements almost entirely restricted to arthropod hosts, considerably expanding the set of MULEs known from this group, suggesting that many additional MULEs may await discovery from further arthropod genomes. In several cases, close relationships occur between MULEs recovered from distantly related host organisms, suggesting that horizontal transfer events may have played an important role in the evolution of the group. However, it is apparent that MULEs from plants remain separate from MULEs identified from other host groups. MULE structure varies considerably across phylogeny, and TIR length is shown to vary greatly both within and between MULE groups. Our phylogeny suggests that MULE diversity is clustered in well-supported groups, typically according to host taxonomy. With reference to this, we make suggestions on how MULE diversity can be partitioned to provide a robust taxonomic framework.

Conclusions

Our study represents a considerable advance in the understanding of MULE diversity, host range and evolution, and provides a taxonomic framework for the classification of further MULE elements that await discovery. Our findings also raise a number of questions relating to MULE biology, suggesting that this group will provide a rich avenue for future study.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0153-8) contains supplementary material, which is available to authorized users.

Opportunities and limitations of reduced representation bisulfite sequencing in plant ecological epigenomics

Contents Summary 738 I. Introduction 738 II. RRBS loci as genome-wide epigenetic markers 739 III. Exploiting functional annotation of RRBS loci 739 IV. Limitations of RRBS methods for nonmodel species 740 V. Maximising the impact of RRBS in plants 741 VI. Conclusions 741 Acknowledgements 741 SUMMARY: Investigating the features and implications of epigenetic mechanisms across the breadth of organisms and ecosystems is important for understanding the ecological relevance of epigenetics. Several cost-effective reduced representation bisulfite sequencing approaches (RRBS) have been recently developed and applied to different organisms that lack a well annotated reference genome. These new approaches improve the assessment of epigenetic diversity in ecological settings and may provide functional insights. We assess here the opportunities and limitations of RRBS in nonmodel plant species. Well thought out experimental designs that include complementary gene expression studies, and the improvement of genomics resources for the target group, promise to maximize the effect of future RRBS studies.

Evolution of sequence-specific anti-silencing systems in Arabidopsis

The arms race between parasitic sequences and their hosts is a major driving force for evolution of gene control systems. Since transposable elements (TEs) are potentially deleterious, eukaryotes silence them by epigenetic mechanisms such as DNA methylation. Little is known about how TEs counteract silencing to propagate during evolution. Here, we report behavior of sequence-specific anti-silencing proteins used by Arabidopsis TEs and evolution of those proteins and their target sequences. We show that VANC, a TE-encoded anti-silencing protein, induces extensive DNA methylation loss throughout TEs. Related VANC proteins have evolved to hypomethylate TEs of completely different spectra. Targets for VANC proteins often form tandem repeats, which vary considerably between related TEs. We propose that evolution of VANC proteins and their targets allow propagation of TEs while causing minimal host damage. Our findings provide insight into the evolutionary dynamics of these apparently "selfish" sequences. They also provide potential tools to edit epigenomes in a sequence-specific manner.

Comparative epigenomics in the Brassicaceae reveals two evolutionarily conserved modes of PRC2-mediated gene regulation

Background

Polycomb Repressive Complexes 2 (PRC2) are multi-protein chromatin modifiers that are evolutionarily conserved among eukaryotes and play key roles in the regulation of gene expression, notably through the trimethylation of lysine 27 of histone H3 (H3K27me3). Although PRC2-mediated gene regulation has been studied in many organisms, few studies have explored in depth the evolutionary conservation of PRC2 targets.

Results

Here, we compare the H3K27me3 epigenomic profiles for the two closely related species Arabidopsis thaliana and Arabidopsis lyrata and the more distant species Arabis alpina, three Brassicaceae that diverged from each other within the past 24 million years.

Using a robust set of gene orthologs present in the three species, we identify two classes of evolutionarily conserved PRC2 targets, which are characterized by either developmentally plastic or developmentally constrained H3K27me3 marking across species. Constrained H3K27me3 marking is associated with higher conservation of promoter sequence information content and higher nucleosome occupancy compared to plastic H3K27me3 marking. Moreover, gene orthologs with constrained H3K27me3 marking exhibit a higher degree of tissue specificity and tend to be involved in developmental functions, whereas gene orthologs with plastic H3K27me3 marking preferentially encode proteins associated with metabolism and stress responses. In addition, gene orthologs with constrained H3K27me3 marking are the predominant contributors to higher-order chromosome organization.

Conclusions

Our findings indicate that developmentally plastic and constrained H3K27me3 marking define two evolutionarily conserved modes of PRC2-mediated gene regulation that are associated with distinct selective pressures operating at multiple scales, from DNA sequence to gene function and chromosome architecture.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1333-9) contains supplementary material, which is available to authorized users.

Profuse evolutionary diversification and speciation on volcanic islands: transposon instability and amplification bursts explain the genetic paradox

BACKGROUND

Species-rich adaptive radiations arising from rare plant and animal colonizers are common on remote volcanic archipelagoes. However, they present a paradox. The severe genetic bottleneck of founder events and effects of inbreeding depression, coupled with the inherently stressful volcanic environment, would seem to predict reduced evolutionary potential and increased risk of extinction, rather than rapid adaptive divergence and speciation. Significantly, eukaryotic genomes harbor many families of transposable elements (TEs) that are mobilized by genome shock; these elements may be the primary drivers of genetic reorganization and speciation on volcanic islands.

PRESENTATION OF THE HYPOTHESIS

Here I propose that a central factor in the spectacular radiation and diversification of the endemic Hawaiian Drosophila and other terrestrial lineages on the Hawaiian and other oceanic islands has been repeated bursts of transposition of multiple TEs induced by the unique ecological features of volcanic habitats. Founder individuals and populations on remote volcanic islands experience significant levels of physiological and genomic stress as a consequence of both biotic and abiotic factors. This results in disruption of the usual epigenetic suppression of TEs, unleashing them to proliferate and spread, which in turn gives rise to novel genetic variation and remodels genomic regulatory circuits, facilitating rapid morphological, ecological and behavioral change, and adaptive radiation.

TESTING THE HYPOTHESIS

To obtain empirical support for the hypothesis, test organisms should be exposed to prolonged heat stress, high levels of carbon dioxide and other volcanic gases, along with inbreeding. Data from subsequent whole genome sequencing and bioinformatics screening for TE numbers and locations would then be compared with initial pre-exposure TE information for the test strains, a labor-intensive project. Several predicted outcomes arising from the hypothesis are discussed. Currently available data are consistent with the proposed concept of stress-induced TE mobilization as a trigger of evolutionary diversification and speciation on volcanic islands.

IMPLICATIONS OF THE HYPOTHESIS

The main implication is that both TEs and species should proliferate at a much higher rate on volcanic islands than elsewhere. Second, the evolvability of a lineage may correlate with the abundance and distribution of TEs in the genome. Successful colonizers of volcanic habitats with high genomic proportions of TEs may be best poised to found a speciose lineage that gives rise to a dramatic adaptive radiation. Colonizers that are depauperate in TEs are likely to be evolutionarily constrained and diversify little, if at all.

REVIEWERS

This article was reviewed by Dr. James Shapiro and Dr. Wolfgang Miller (nominated by Editorial Board member Dr. I. King Jordan).

Activation of Tag1 transposable elements in Arabidopsis dedifferentiating cells and their regulation by CHROMOMETHYLASE 3-mediated CHG methylation

Dedifferentiation, that is, the acquisition of stem cell-like state, commonly induced by stress (e.g., protoplasting), is characterized by open chromatin conformation, a chromatin state that could lead to activation of transposable elements (TEs). Here, we studied the activation of the Arabidopsis class II TE Tag1, in which two copies, situated close to each other (near genes) on chromosome 1 are found in Landsberg erecta (Ler) but not in Columbia (Col). We first transformed protoplasts with a construct in which a truncated Tag1 (ΔTag1 non-autonomous) blocks the expression of a reporter gene AtMBD5-GFP and found a relatively high ectopic excision of ΔTag1 accompanied by expression of AtMBD5-GFP in protoplasts derived from Ler compared to Col; further increase was observed in ddm1 (decrease in DNA methylation1) protoplasts (Ler background). Ectopic excision was associated with transcription of the endogenous Tag1 and changes in histone H3 methylation at the promoter region. Focusing on the endogenous Tag1 elements we found low level of excision in Ler protoplasts, which was slightly and strongly enhanced in ddm1 and cmt3 (chromomethylase3) protoplasts, respectively, concomitantly with reduction in Tag1 gene body (GB) CHG methylation and increased Tag1 transcription; strong activation of Tag1 was also observed in cmt3 leaves. Notably, in cmt3, but not in ddm1, Tag1 elements were excised out from their original sites and transposed elsewhere in the genome. Our results suggest that dedifferentiation is associated with Tag1 activation and that CMT3 rather than DDM1 plays a central role in restraining Tag1 activation via inducing GB CHG methylation.

Phylogenetic and Genomic Analyses Resolve the Origin of Important Plant Genes Derived from Transposable Elements

Once perceived as merely selfish, transposable elements (TEs) are now recognized as potent agents of adaptation. One way TEs contribute to evolution is through TE exaptation, a process whereby TEs, which persist by replicating in the genome, transform into novel host genes, which persist by conferring phenotypic benefits. Known exapted TEs (ETEs) contribute diverse and vital functions, and may facilitate punctuated equilibrium, yet little is known about this process. To better understand TE exaptation, we designed an approach to resolve the phylogenetic context and timing of exaptation events and subsequent patterns of ETE diversification. Starting with known ETEs, we search in diverse genomes for basal ETEs and closely related TEs, carefully curate the numerous candidate sequences, and infer detailed phylogenies. To distinguish TEs from ETEs, we also weigh several key genomic characteristics including repetitiveness, terminal repeats, pseudogenic features, and conserved domains. Applying this approach to the well-characterized plant ETEs MUG and FHY3, we show that each group is paraphyletic and we argue that this pattern demonstrates that each originated in not one but multiple exaptation events. These exaptations and subsequent ETE diversification occurred throughout angiosperm evolution including the crown group expansion, the angiosperm radiation, and the primitive evolution of angiosperms. In addition, we detect evidence of several putative novel ETE families. Our findings support the hypothesis that TE exaptation generates novel genes more frequently than is currently thought, often coinciding with key periods of evolution.

Helitrons, the Eukaryotic Rolling-circle Transposable Elements

Helitrons, the eukaryotic rolling-circle transposable elements, are widespread but most prevalent among plant and animal genomes. Recent studies have identified three additional coding and structural variants of Helitrons called Helentrons, Proto-Helentron, and Helitron2. Helitrons and Helentrons make up a substantial fraction of many genomes where nonautonomous elements frequently outnumber the putative autonomous partner. This includes the previously ambiguously classified DINE-1-like repeats, which are highly abundant in Drosophila and many other animal genomes. The purpose of this review is to summarize what we have learned about Helitrons in the decade since their discovery. First, we describe the history of autonomous Helitrons, and their variants. Second, we explain the common coding features and difference in structure of canonical Helitrons versus the endonuclease-encoding Helentrons. Third, we review how Helitrons and Helentrons are classified and discuss why the system used for other transposable element families is not applicable. We also touch upon how genome-wide identification of candidate Helitrons is carried out and how to validate candidate Helitrons. We then shift our focus to a model of transposition and the report of an excision event. We discuss the different proposed models for the mechanism of gene capture. Finally, we will talk about where Helitrons are found, including discussions of vertical versus horizontal transfer, the propensity of Helitrons and Helentrons to capture and shuffle genes and how they impact the genome. We will end the review with a summary of open questions concerning the biology of this intriguing group of transposable elements.

Creating Order from Chaos: Epigenome Dynamics in Plants with Complex Genomes

Flowering plants have strikingly distinct genomes, although they contain a similar suite of expressed genes. The diversity of genome structures and organization is largely due to variation in transposable elements (TEs) and whole-genome duplication (WGD) events. We review evidence that chromatin modifications and epigenetic regulation are intimately associated with TEs and likely play a role in mediating the effects of WGDs. We hypothesize that the current structure of a genome is the result of various TE bursts and WGDs and it is likely that the silencing mechanisms and the chromatin structure of a genome have been shaped by these events. This suggests that the specific mechanisms targeting chromatin modifications and epigenomic patterns may vary among different species. Many crop species have likely evolved chromatin-based mechanisms to tolerate silenced TEs near actively expressed genes. These interactions of heterochromatin and euchromatin are likely to have important roles in modulating gene expression and variability within species.

Reproductive Mode and the Evolution of Genome Size and Structure in Caenorhabditis Nematodes

The self-fertile nematode worms Caenorhabditis elegans, C. briggsae, and C. tropicalis evolved independently from outcrossing male-female ancestors and have genomes 20-40% smaller than closely related outcrossing relatives. This pattern of smaller genomes for selfing species and larger genomes for closely related outcrossing species is also seen in plants. We use comparative genomics, including the first high quality genome assembly for an outcrossing member of the genus (C. remanei) to test several hypotheses for the evolution of genome reduction under a change in mating system. Unlike plants, it does not appear that reductions in the number of repetitive elements, such as transposable elements, are an important contributor to the change in genome size. Instead, all functional genomic categories are lost in approximately equal proportions. Theory predicts that self-fertilization should equalize the effective population size, as well as the resulting effects of genetic drift, between the X chromosome and autosomes. Contrary to this, we find that the self-fertile C. briggsae and C. elegans have larger intergenic spaces and larger protein-coding genes on the X chromosome when compared to autosomes, while C. remanei actually has smaller introns on the X chromosome than either self-reproducing species. Rather than being driven by mutational biases and/or genetic drift caused by a reduction in effective population size under self reproduction, changes in genome size in this group of nematodes appear to be caused by genome-wide patterns of gene loss, most likely generated by genomic adaptation to self reproduction per se.

MuTAnT: a family of Mutator-like transposable elements targeting TA microsatellites in Medicago truncatula

Transposable elements (TEs) are mobile DNA segments, abundant and dynamic in plant genomes. Because their mobility can be potentially deleterious to the host, a variety of mechanisms evolved limiting that negative impact, one of them being preference for a specific target insertion site. Here, we describe a family of Mutator-like DNA transposons in Medicago truncatula targeting TA microsatellites. We identified 218 copies of MuTAnTs and an element carrying a complete ORF encoding a mudrA-like transposase. Most insertion sites are flanked by a variable number of TA tandem repeats, indicating that MuTAnTs are specifically targeting TA microsatellites. Other TE families flanked by TA repeats (e.g. TAFT elements in maize) were described previously, however we identified the first putative autonomous element sharing that characteristics with a related group of short non-autonomous transposons.

Control of transposable elements in Arabidopsis thaliana

Arabidopsis thaliana serves as a very good model organism to investigate the control of transposable elements (TEs) by genetic and genomic approaches. As TE movements are potentially deleterious to the hosts, hosts silence TEs by epigenetic mechanisms, such as DNA methylation. DNA methylation is controlled by DNA methyltransferases and other regulators, including histone modifiers and chromatin remodelers. RNAi machinery directs DNA methylation to euchromatic TEs, which is under developmental control. In addition to the epigenetic controls, some TEs are controlled by environmental factors. TEs often affect expression of nearby genes, providing evolutionary sources for epigenetic, developmental, and environmental gene controls, which could even be beneficial for the host.

Global DNA cytosine methylation as an evolving trait: phylogenetic signal and correlated evolution with genome size in angiosperms

DNA cytosine methylation is a widespread epigenetic mechanism in eukaryotes, and plant genomes commonly are densely methylated. Genomic methylation can be associated with functional consequences such as mutational events, genomic instability or altered gene expression, but little is known on interspecific variation in global cytosine methylation in plants. In this paper, we compare global cytosine methylation estimates obtained by HPLC and use a phylogenetically-informed analytical approach to test for significance of evolutionary signatures of this trait across 54 angiosperm species in 25 families. We evaluate whether interspecific variation in global cytosine methylation is statistically related to phylogenetic distance and also whether it is evolutionarily correlated with genome size (C-value). Global cytosine methylation varied widely between species, ranging between 5.3% (Arabidopsis) and 39.2% (Narcissus). Differences between species were related to their evolutionary trajectories, as denoted by the strong phylogenetic signal underlying interspecific variation. Global cytosine methylation and genome size were evolutionarily correlated, as revealed by the significant relationship between the corresponding phylogenetically independent contrasts. On average, a ten-fold increase in genome size entailed an increase of about 10% in global cytosine methylation. Results show that global cytosine methylation is an evolving trait in angiosperms whose evolutionary trajectory is significantly linked to changes in genome size, and suggest that the evolutionary implications of epigenetic mechanisms are likely to vary between plant lineages.

The structure and retrotransposition mechanism of LTR-retrotransposons in the asexual yeast Candida albicans

Retrotransposons constitute a major part of the genome in a number of eukaryotes. Long-terminal repeat (LTR) retrotransposons are one type of the retrotransposons. Candida albicans have 34 distinct LTR-retrotransposon families. They respectively belong to the Ty1/copia and Ty3/gypsy groups which have been extensively studied in the model yeast Saccharomyces cerevisiae. LTR-retrotransposons carry two LTRs flanking a long internal protein-coding domain, open reading frames. LTR-retrotransposons use RNA as intermediate to synthesize double-stranded DNA copies. In this article, we describe the structure feature, retrotransposition mechanism and the influence on organism diversity of LTR retrotransposons in C. albicans. We also discuss the relationship between pathogenicity and LTR retrotransposons in C. albicans.

Detection of genomic variations and DNA polymorphisms and impact on analysis of meiotic recombination and genetic mapping

DNA polymorphisms are important markers in genetic analyses and are increasingly detected by using genome resequencing. However, the presence of repetitive sequences and structural variants can lead to false positives in the identification of polymorphic alleles. Here, we describe an analysis strategy that minimizes false positives in allelic detection and present analyses of recently published resequencing data from Arabidopsis meiotic products and individual humans. Our analysis enables the accurate detection of sequencing errors, small insertions and deletions (indels), and structural variants, including large reciprocal indels and copy number variants, from comparisons between the resequenced and reference genomes. We offer an alternative interpretation of the sequencing data of meiotic products, including the number and type of recombination events, to illustrate the potential for mistakes in single-nucleotide polymorphism calling. Using these examples, we propose that the detection of DNA polymorphisms using resequencing data needs to account for nonallelic homologous sequences.

Environmental stress activation of plant long-terminal repeat retrotransposons

Genomic retrotransposons (RTs) are major components of most plant genomes. They spread throughout the genomes by a process termed retrotransposition, which consists of reverse transcription and reinsertion of the copied element into a new genomic location (a copy-and-paste system). Abiotic and biotic stresses activate long-terminal repeat (LTR) RTs in photosynthetic eukaryotes from algae to angiosperms. LTR RTs could represent a threat to the integrity of host genomes because of their activity and mutagenic potential by epigenetic regulation. Host genomes have developed mechanisms to control the activity of the retroelements and their mutagenic potential. Some LTR RTs escape these defense mechanisms, and maintain their ability to be activated and transpose as a result of biotic or abiotic stress stimuli. These stimuli include pathogen infection, mechanical damage, in vitro tissue culturing, heat, drought and salt stress, generation of doubled haploids, X-ray irradiation and many others. Reactivation of LTR RTs differs between different plant genomes. The expression levels of reactivated RTs are influenced by the transcriptional and post-transcriptional gene silencing mechanisms (e.g. DNA methylation, heterochromatin formation and RNA interference). Moreover, the insertion of RTs (e.g. Triticum aestivum L. Wis2-1A) into or next to coding regions of the host genome can generate changes in the expression of adjacent host genes of the host. In this paper, we review the ways that plant genomic LTR RTs are activated by environmental stimuli to affect restructuring and diversification of the host genome.

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.

Scanning of transposable elements and analyzing expression of transposase genes of sweet potato [Ipomoea batatas]

BACKGROUND

Transposable elements (TEs) are the most abundant genomic components in eukaryotes and affect the genome by their replications and movements to generate genetic plasticity. Sweet potato performs asexual reproduction generally and the TEs may be an important genetic factor for genome reorganization. Complete identification of TEs is essential for the study of genome evolution. However, the TEs of sweet potato are still poorly understood because of its complex hexaploid genome and difficulty in genome sequencing. The recent availability of the sweet potato transcriptome databases provides an opportunity for discovering and characterizing the expressed TEs.

METHODOLOGY/PRINCIPAL FINDINGS

We first established the integrated-transcriptome database by de novo assembling four published sweet potato transcriptome databases from three cultivars in China. Using sequence-similarity search and analysis, a total of 1,405 TEs including 883 retrotransposons and 522 DNA transposons were predicted and categorized. Depending on mapping sets of RNA-Seq raw short reads to the predicted TEs, we compared the quantities, classifications and expression activities of TEs inter- and intra-cultivars. Moreover, the differential expressions of TEs in seven tissues of Xushu 18 cultivar were analyzed by using Illumina digital gene expression (DGE) tag profiling. It was found that 417 TEs were expressed in one or more tissues and 107 in all seven tissues. Furthermore, the copy number of 11 transposase genes was determined to be 1-3 copies in the genome of sweet potato by Real-time PCR-based absolute quantification.

CONCLUSIONS/SIGNIFICANCE

Our result provides a new method for TE searching on species with transcriptome sequences while lacking genome information. The searching, identification and expression analysis of TEs will provide useful TE information in sweet potato, which are valuable for the further studies of TE-mediated gene mutation and optimization in asexual reproduction. It contributes to elucidating the roles of TEs in genome evolution.

Acquisition of an Archaea-like ribonuclease H domain by plant L1 retrotransposons supports modular evolution

Although a variety of non-LTR retrotransposons of the L1 superfamily have been found in plant genomes over recent decades, their diversity, distribution, and evolution have yet to be analyzed in depth. Here, we perform comprehensive comparative and evolutionary analyses of L1 retrotransposons from 29 genomes of land plants covering a wide range of taxa. We identify numerous L1 elements in these genomes and detect a striking diversity of their domain composition. We show that all known land plant L1 retrotransposons can be grouped into five major families based on their phylogenetic relationships and domain composition. Moreover, we trace the putative evolution timeline that created the current variants and reveal that evolutionary events included losses and acquisitions of diverse putative RNA-binding domains and the acquisition of an Archaea-like ribonuclease H (RNH) domain. We also show that the latter RNH domain is autonomously active in vitro and speculate that retrotransposons may play a role in the horizontal transfer of RNH between plants, Archaea, and bacteria. The acquisition of an Archaea-like RNH domain by plant L1 retrotransposons negates the hypothesis that RNH domains in non-LTR retrotransposons have a single origin and provides evidence that acquisition happened at least twice. Together, our data indicate that the evolution of the investigated retrotransposons can be mainly characterized by repeated events of domain rearrangements and identify modular evolution as a major trend in the evolution of plant L1 retrotransposons.

Site-specific gene targeting using transcription activator-like effector (TALE)-based nuclease in Brassica oleracea

Site-specific recognition modules with DNA nuclease have tremendous potential as molecular tools for genome targeting. The type III transcription activator-like effectors (TALEs) contain a DNA binding domain consisting of tandem repeats that can be engineered to bind user-defined specific DNA sequences. We demonstrated that customized TALE-based nucleases (TALENs), constructed using a method called "unit assembly", specifically target the endogenous FRIGIDA gene in Brassica oleracea L. var. capitata L. The results indicate that the TALENs bound to the target site and cleaved double-strand DNA in vitro and in vivo, whereas the effector binding elements have a 23 bp spacer. The T7 endonuclease I assay and sequencing data show that TALENs made double-strand breaks, which were repaired by a non-homologous end-joining pathway within the target sequence. These data show the feasibility of applying customized TALENs to target and modify the genome with deletions in those organisms that are still in lacking gene target methods to provide germplasms in breeding improvement.

Mobilization of a plant transposon by expression of the transposon-encoded anti-silencing factor

Transposable elements (TEs) have a major impact on genome evolution, but they are potentially deleterious, and most of them are silenced by epigenetic mechanisms, such as DNA methylation. Here, we report the characterization of a TE encoding an activity to counteract epigenetic silencing by the host. In Arabidopsis thaliana, we identified a mobile copy of the Mutator-like element (MULE) with degenerated terminal inverted repeats (TIRs). This TE, named Hiun (Hi), is silent in wild-type plants, but it transposes when DNA methylation is abolished. When a Hi transgene was introduced into the wild-type background, it induced excision of the endogenous Hi copy, suggesting that Hi is the autonomously mobile copy. In addition, the transgene induced loss of DNA methylation and transcriptional activation of the endogenous Hi. Most importantly, the trans-activation of Hi depends on a Hi-encoded protein different from the conserved transposase. Proteins related to this anti-silencing factor, which we named VANC, are widespread in the non-TIR MULEs and may have contributed to the recent success of these TEs in natural Arabidopsis populations.

New insights into nested long terminal repeat retrotransposons in Brassica species

Long terminal repeat (LTR) retrotransposons, one of the foremost types of transposons, continually change or modify gene function and reorganize the genome through bursts of dramatic proliferation. Many LTR-TEs preferentially insert within other LTR-TEs, but the cause and evolutionary significance of these nested LTR-TEs are not well understood. In this study, a total of 1.52Gb of Brassica sequence containing 2020 bacterial artificial chromosomes (BACs) was scanned, and six bacterial artificial chromosome (BAC) clones with extremely nested LTR-TEs (LTR-TEs density: 7.24/kb) were selected for further analysis. The majority of the LTR-TEs in four of the six BACs were found to be derived from the rapid proliferation of retrotransposons originating within the BAC regions, with only a few LTR-TEs originating from the proliferation and insertion of retrotransposons from outside the BAC regions approximately 5-23Mya. LTR-TEs also preferably inserted into TA-rich repeat regions. Gene prediction by Genescan identified 207 genes in the 0.84Mb of total BAC sequences. Only a few genes (3/207) could be matched to the Brassica expressed sequence tag (EST) database, indicating that most genes were inactive after retrotransposon insertion. Five of the six BACs were putatively centromeric. Hence, nested LTR-TEs in centromere regions are rapidly duplicated, repeatedly inserted, and act to suppress activity of genes and to reshuffle the structure of the centromeric sequences. Our results suggest that LTR-TEs burst and proliferate on a local scale to create nested LTR-TE regions, and that these nested LTR-TEs play a role in the formation of centromeres.

Transposon variants and their effects on gene expression in Arabidopsis

Transposable elements (TEs) make up the majority of many plant genomes. Their transcription and transposition is controlled through siRNAs and epigenetic marks including DNA methylation. To dissect the interplay of siRNA–mediated regulation and TE evolution, and to examine how TE differences affect nearby gene expression, we investigated genome-wide differences in TEs, siRNAs, and gene expression among three Arabidopsis thaliana accessions. Both TE sequence polymorphisms and presence of linked TEs are positively correlated with intraspecific variation in gene expression. The expression of genes within 2 kb of conserved TEs is more stable than that of genes next to variant TEs harboring sequence polymorphisms. Polymorphism levels of TEs and closely linked adjacent genes are positively correlated as well. We also investigated the distribution of 24-nt-long siRNAs, which mediate TE repression. TEs targeted by uniquely mapping siRNAs are on average farther from coding genes, apparently because they more strongly suppress expression of adjacent genes. Furthermore, siRNAs, and especially uniquely mapping siRNAs, are enriched in TE regions missing in other accessions. Thus, targeting by uniquely mapping siRNAs appears to promote sequence deletions in TEs. Overall, our work indicates that siRNA–targeting of TEs may influence removal of sequences from the genome and hence evolution of gene expression in plants.

Transposable elements (TEs) are selfish DNA sequences. Together with their immobilized derivatives, they account for a large fraction of eukaryotic genomes. TEs can affect nearby gene activity, either directly by disrupting regulatory sequences or indirectly through the host mechanisms used to prevent TE proliferation. A comparison of Arabidopsis thaliana genomes reveals rapid TE degeneration. We asked what drives TE degeneration and how often TE variation affects nearby gene expression. To answer these questions, we studied the interplay between TEs, DNA sequence variation, and short interfering RNAs (siRNAs) in three A. thaliana strains. We find sequence variation in genes and adjacent TEs to be correlated, from which we conclude either that TEs insert more often near polymorphic genes or that TEs next to polymorphic genes are less efficiently purged from the genome. We also noticed that processes that cause deletions within TEs and ones that silence TEs appear to be linked, because siRNA targeting is a predictor of sequence loss in accessions. Our work provides insight into the contribution of TEs to gene expression plasticity, and it links TE silencing mechanisms to the evolution of TE variation between genomes, thereby linking TE silencing mechanisms to expression plasticity.

Control of transposon activity by a histone H3K4 demethylase in rice

Transposable elements (TEs) are ubiquitously present in plant genomes and often account for significant fractions of the nuclear DNA. For example, roughly 40% of the rice genome consists of TEs, many of which are retrotransposons, including 14% LTR- and ∼1% non-LTR retrotransposons. Despite their wide distribution and abundance, very few TEs have been found to be transpositional, indicating that TE activities may be tightly controlled by the host genome to minimize the potentially mutagenic effects associated with active transposition. Consistent with this notion, a growing body of evidence suggests that epigenetic silencing pathways such as DNA methylation, RNA interference, and H3K9me2 function collectively to repress TE activity at the transcriptional and posttranscriptional levels. It is not yet clear, however, whether the removal of histone modifications associated with active transcription is also involved in TE silencing. Here, we show that the rice protein JMJ703 is an active H3K4-specific demethylase required for TEs silencing. Impaired JMJ703 activity led to elevated levels of H3K4me3, the misregulation of numerous endogenous genes, and the transpositional reactivation of two families of non-LTR retrotransposons. Interestingly, loss of JMJ703 did not affect TEs (such as Tos17) previously found to be silenced by other epigenetic pathways. These results indicate that the removal of active histone modifications is involved in TE silencing and that different subsets of TEs may be regulated by distinct epigenetic pathways.

TIRfinder: A Web Tool for Mining Class II Transposons Carrying Terminal Inverted Repeats

Transposable elements (TEs) can be found in virtually all known genomes; plant genomes are exceptionally rich in this kind of dispersed repetitive sequences. Current knowledge on TE proliferation dynamics places them among the main forces of molecular evolution. Therefore efficient tools to analyze TE distribution in genomes are needed that would allow for comparative genomics studies and for studying TE dynamics in a genome. This was our main motivation underpinning TIRfinder construction–-an efficient tool for mining class II TEs carrying terminal inverted repeats. TIRfinder takes as an input a genomic sequence and information on structural properties of a TE family, and identifies all TEs in the genome showing the desired structural characteristics. The efficiency and small memory requirements of our approach stem from the use of suffix trees to identify all DNA segments surrounded by user-specified terminal inverse repeats (TIR) and target site duplications (TSD) which together constitute a mask. On the other hand, the flexibility of the notion of the TIR/TSD mask makes it possible to use the tool for de novo detection. The main advantages of TIRfinder are its speed, accuracy and convenience of use for biologists. A web-based interface is freely available at http:/bioputer.mimuw.edu.pl/tirfindertool/ .

A gene family derived from transposable elements during early angiosperm evolution has reproductive fitness benefits in Arabidopsis thaliana

The benefits of ever-growing numbers of sequenced eukaryotic genomes will not be fully realized until we learn to decipher vast stretches of noncoding DNA, largely composed of transposable elements. Transposable elements persist through self-replication, but some genes once encoded by transposable elements have, through a process called molecular domestication, evolved new functions that increase fitness. Although they have conferred numerous adaptations, the number of such domesticated transposable element genes remains unknown, so their evolutionary and functional impact cannot be fully assessed. Systematic searches that exploit genomic signatures of natural selection have been employed to identify potential domesticated genes, but their predictions have yet to be experimentally verified. To this end, we investigated a family of domesticated genes called MUSTANG (MUG), identified in a previous bioinformatic search of plant genomes. We show that MUG genes are functional. Mutants of Arabidopsis thaliana MUG genes yield phenotypes with severely reduced plant fitness through decreased plant size, delayed flowering, abnormal development of floral organs, and markedly reduced fertility. MUG genes are present in all flowering plants, but not in any non-flowering plant lineages, such as gymnosperms, suggesting that the molecular domestication of MUG may have been an integral part of early angiosperm evolution. This study shows that systematic searches can be successful at identifying functional genetic elements in noncoding regions and demonstrates how to combine systematic searches with reverse genetics in a fruitful way to decipher eukaryotic genomes.

The genomes of complex organisms are mostly made up not of ordinary genes but of transposable elements. Transposable elements have been called “selfish DNA” because they normally persist by copying themselves, not by helping the organism to survive or reproduce. Yet transposable elements can help organisms to evolve; for instance, transposable element genes sometimes acquire new functions that do benefit the organism. Because they are difficult to distinguish from transposable elements, little is known about these “domesticated genes.” Although studies have attempted to identify them computationally, the predictions have not been verified experimentally. Here, we examine some of the first domesticated genes to be predicted computationally, the MUSTANG family of plant genes. We show that the predictions were correct: MUSTANGs are, like ordinary genes, functional. MUSTANG mutations result in serious defects in how plants grow, flower, and reproduce. Since they are present only in flowering plants, MUSTANG probably originated when flowers first evolved, perhaps taking on a key role. This study is important both because it shows that MUSTANG is critical to plant fitness and because, in the future, a similar approach can be used to find additional domesticated genes and to better understand how transposable elements contribute to evolution.

PpPIF-1: first isolated full-length PIF-like element from the bamboo Phyllostachys pubescens

PIF-like elements are the first-described members of a recently discovered and widespread superfamily of DNA transposons, named PIF/Harbinger. Complete and partial PIF-like elements have been isolated from hundreds of plant species. Previously, we identified 139 partial PIF-like transposases in the Bambusoideae, of which three were from the bamboo species Phyllostachys pubescens. Here we report identification and isolation of the first full-length PIF-like element (PpPIF-1) from P. pubescens; identification was made by chromosome walking, based on a modified magnetic enrichment procedure that allows efficient cloning of flanking sequences up to 3 kb in length. PpPIF-1 is 5953 bp in length, with 20-bp imperfect inverted terminal repeats and 3-bp target site duplications. This element contains two open reading frames, one encoding a putative transposase, including the complete DDE-domain typical of PIF/Harbinger elements from plants, and the other encoding a DNA-binding protein. There are seven termination codons and two frameshift mutations in the open reading frames, probably due to vertical inactivation.

Centromere-targeted de novo integrations of an LTR retrotransposon of Arabidopsis lyrata

The plant genome evolves with rapid proliferation of LTR-type retrotransposons, which is associated with their clustered accumulation in gene-poor regions, such as centromeres. Despite their major role for plant genome evolution, no mobile LTR element with targeted integration into gene-poor regions has been identified in plants. Here, we report such targeted integrations de novo. We and others have previously shown that an ATCOPIA93 family retrotransposon in Arabidopsis thaliana is mobilized when the DNA methylation machinery is compromised. Although ATCOPIA93 family elements are low copy number in the wild-type A. thaliana genome, high-copy-number related elements are found in the wild-type Arabidopsis lyrata genome, and they show centromere-specific localization. To understand the mechanisms for the clustered accumulation of the A. lyrata elements directly, we introduced one of them, named Tal1 (Transposon of Arabidopsis lyrata 1), into A. thaliana by transformation. The introduced Tal1 was retrotransposed in A. thaliana, and most of the retrotransposed copies were found in centromeric repeats of A. thaliana, suggesting targeted integration. The targeted integration is especially surprising because the centromeric repeat sequences differ considerably between A. lyrata and A. thaliana. Our results revealed unexpectedly dynamic controls for evolution of the transposon-rich heterochromatic regions.

Evolutionary Implications of Mechanistic Models of TE-Mediated Hybrid Incompatibility

New models of TE repression in plants (specifically Arabidopsis) have suggested specific mechanisms by which TE misregulation in hybrids might result in the expression of hybrid inviability. If true, these models suggest as yet undescribed consequences for (1) mechanistic connections between hybrid problems expressed at different postzygotic stages (e.g., inviability versus sterility), (2) the predicted strength, stage, and direction of isolation between diverging lineages that differ in TE activity, and (3) the association between species attributes that influence TE dynamics (e.g., mode of reproduction, geographical structure) and the rate at which they could accumulate incompatibilities. In this paper, we explore these implications and outline future empirical directions for generating data necessary to evaluate them.

The predominantly selfing plant Arabidopsis thaliana experienced a recent reduction in transposable element abundance compared to its outcrossing relative Arabidopsis lyrata

Background

Transposable elements (TEs) are major contributors to genome evolution. One factor that influences their evolutionary dynamics is whether their host reproduces through selfing or through outcrossing. According to the recombinational spreading hypothesis, for instance, TEs can spread more easily in outcrossing species through recombination, and should thus be less abundant in selfing species. We here studied the distribution and evolutionary dynamics of TE families in the predominantly selfing plant Arabidopsis thaliana and its close outcrossing relative Arabidopsis lyrata on a genome-wide scale. We characterized differences in TE abundance between them and asked which, if any, existing hypotheses about TE abundances may explain these differences.

Results

We identified 1,819 TE families representing all known classes of TEs in both species, and found three times more copies in the outcrossing A. lyrata than in the predominantly selfing A. thaliana, as well as ten times more TE families unique to A. lyrata. On average, elements in A. lyrata are younger than elements in A. thaliana. In particular, A. thaliana shows a marked decrease in element number that occurred during the most recent 10% of the time interval since A. thaliana split from A. lyrata. This most recent period in the evolution of A. thaliana started approximately 500,000 years ago, assuming a splitting time of 5 million years ago, and coincides with the time at which predominant selfing originated.

Conclusions

Our results indicate that the mating system may be important for determining TE copy number, and that selfing species are likely to have fewer TEs.

Identification of an active Mutator-like element (MULE) in rice (Oryza sativa)

Transposable elements (TEs) represent an important fraction of plant genomes and play a significant role in gene and genome evolution. Among all TE superfamilies discovered in plants, Mutator from maize (Zea mays) is the most active and mutagenic element. Mutator-like elements (MULEs) were identified in a wide range of plants. However, only few active MULEs have been reported, and the transposition mechanism of the elements is still poorly understood. In this study, an active MULE named Os3378 was discovered in rice (Oryza sativa) by a combination of computational and experimental approaches. The four newly identified Os3378 elements share more than 98% sequence identity between each other, and all of them encode transposases without any deletion derivatives, indicating their capability of autonomous transposition. Os3378 is present in the rice species with AA genome type but is absent in other non-AA genome species. A new insertion of Os3378 was identified in a rice somaclonal mutant Z418, and the element remained active in the descendants of the mutant for more than ten generations. Both germinal and somatic excision events of Os3378 were observed, and no footprint was detected after excision. Furthermore, the occurrence of somatic excision of Os3378 appeared to be associated with plant developmental stages and tissue types. Taken together, Os3378 is a unique active element in rice, which provides a valuable resource for further studying of transposition mechanism and evolution of MULEs.

Parsimonious higher-order hidden Markov models for improved array-CGH analysis with applications to Arabidopsis thaliana

Array-based comparative genomic hybridization (Array-CGH) is an important technology in molecular biology for the detection of DNA copy number polymorphisms between closely related genomes. Hidden Markov Models (HMMs) are popular tools for the analysis of Array-CGH data, but current methods are only based on first-order HMMs having constrained abilities to model spatial dependencies between measurements of closely adjacent chromosomal regions. Here, we develop parsimonious higher-order HMMs enabling the interpolation between a mixture model ignoring spatial dependencies and a higher-order HMM exhaustively modeling spatial dependencies. We apply parsimonious higher-order HMMs to the analysis of Array-CGH data of the accessions C24 and Col-0 of the model plant Arabidopsis thaliana. We compare these models against first-order HMMs and other existing methods using a reference of known deletions and sequence deviations. We find that parsimonious higher-order HMMs clearly improve the identification of these polymorphisms. Moreover, we perform a functional analysis of identified polymorphisms revealing novel details of genomic differences between C24 and Col-0. Additional model evaluations are done on widely considered Array-CGH data of human cell lines indicating that parsimonious HMMs are also well-suited for the analysis of non-plant specific data. All these results indicate that parsimonious higher-order HMMs are useful for Array-CGH analyses. An implementation of parsimonious higher-order HMMs is available as part of the open source Java library Jstacs (www.jstacs.de/index.php/PHHMM).

Array-based comparative genomics is a standard approach for the identification of DNA copy number polymorphisms between closely related genomes. The huge amounts of data produced by these experiments require efficient and accurate bioinformatics tools for the identification of copy number polymorphisms. Hidden Markov Models (HMMs) are frequently used for analyzing such data sets, but current models are based on first-order HMMs only having limited capabilities to model spatial dependencies between measurements of closely adjacent chromosomal regions. We develop parsimonious higher-order HMMs enabling the interpolation between a mixture model ignoring spatial dependencies and a higher-order HMM exhaustively modeling these dependencies to overcome this limitation. In an in-depth case study with Arabidopsis thaliana, we find that parsimonious higher-order HMMs clearly improve the identification of copy number polymorphisms in comparison to standard first-order HMMs and other frequently used methods. Functional analysis of identified polymorphisms revealed details of genomic differences between the accessions C24 and Col-0 of Arabidopsis thaliana. An additional study on human cell lines further indicates that parsimonious HMMs are well-suited for the analysis of Array-CGH data.

Genome evolution in diploid and tetraploid Coffea species as revealed by comparative analysis of orthologous genome segments

Sequence comparison of orthologous regions enables estimation of the divergence between genomes, analysis of their evolution and detection of particular features of the genomes, such as sequence rearrangements and transposable elements. Despite the economic importance of Coffea species, little genomic information is currently available. Coffea is a relatively young genus that includes more than one hundred diploid species and a single tetraploid species. Three Coffea orthologous regions of 470-900 kb were analyzed and compared: both subgenomes of allotetraploid Coffea arabica (contributed by the diploid species Coffea eugenioides and Coffea canephora) and the genome of diploid C. canephora. Sequence divergence was calculated on global alignments or on coding and non-coding sequences separately. A search for transposable elements detected 43 retrotransposons and 198 transposons in the sequences analyzed. Comparative insertion analysis made it possible to locate 165 TE insertions in the phylogenetic tree of the three genomes/subgenomes. In the tetraploid C. arabica, a homoeologous non-reciprocal transposition (HNRT) was detected and characterized: a 50 kb region of the C. eugenioides derived subgenome replaced the C. canephora derived counterpart. Comparative sequence analysis on three Coffea genomes/subgenomes revealed almost perfect gene synteny, low sequence divergence and a high number of shared transposable elements. Compared to the results of similar analysis in other genera (Aegilops/Triticum and Oryza), Coffea genomes/subgenomes appeared to be dramatically less diverged, which is consistent with the relatively recent radiation of the Coffea genus. Based on nucleotide substitution frequency, the HNRT was dated at 10,000-50,000 years BP, which is also the most recent estimation of the origin of C. arabica.

Roadmap for annotating transposable elements in eukaryote genomes

Current high-throughput techniques have made it feasible to sequence even the genomes of non-model organisms. However, the annotation process now represents a bottleneck to genome analysis, especially when dealing with transposable elements (TE). Combined approaches, using both de novo and knowledge-based methods to detect TEs, are likely to produce reasonably comprehensive and sensitive results. This chapter provides a roadmap for researchers involved in genome projects to address this issue. At each step of the TE annotation process, from the identification of TE families to the annotation of TE copies, we outline the tools and good practices to be used.

Miniature inverted-repeat transposable elements (MITEs) have been accumulated through amplification bursts and play important roles in gene expression and species diversity in Oryza sativa

Miniature inverted–repeat transposable elements (MITEs) are predicted to play important roles on genome evolution. We developed a BLASTN-based approach for de novo identification of MITEs and systematically analyzed MITEs in rice genome. The genome of rice cultivar Nipponbare (Oryza sativa ssp. japonica) harbors 178,533 MITE-related sequences classified into 338 families. Pairwise nucleotide diversity and phylogenetic tree analysis indicated that individual MITE families were resulted from one or multiple rounds of amplification bursts. The timing of amplification burst varied considerably between different MITE families or subfamilies. MITEs are associated with 23,623 (58.2%) genes in rice genome. At least 7,887 MITEs are transcribed and more than 3,463 were transcribed with rice genes. The MITE sequences transcribed with rice coding genes form 1,130 pairs of potential natural sense/antisense transcripts. MITEs generate 23.5% (183,837 of 781,885) of all small RNAs identified from rice. Some MITE families generated small RNAs mainly from the terminals, while other families generated small RNAs predominantly from the central region. More than half (51.8%) of the MITE-derived small RNAs were generated exclusively by MITEs located away from genes. Genome-wide analysis showed that genes associated with MITEs have significantly lower expression than genes away from MITEs. Approximately 14.8% of loci with full-length MITEs have presence/absence polymorphism between rice cultivars 93-11 (O. sativa ssp. indica) and Nipponbare. Considering that different sets of genes may be regulated by MITE-derived small RNAs in different genotypes, MITEs provide considerable diversity for O. sativa.

Isolation and characterization of seventy-nine full-length mariner-like transposase genes in the Bambusoideae subfamily

Mariner-like elements (MLEs) are the most diverse and widespread transposable elements, with members of the MLE superfamily found in fungi, plants, ciliates and animals. In a previous study, we characterized 82 MLE transposase gene fragments (average length 383 bp) in 44 bamboo species, indicating that MLEs are widespread, abundant and diverse in the Bambusoideae subfamily. In this study, we isolated 79 full-length MLE transposase genes from 63 bamboo species representing 38 genera in six subtribes mainly found in China. The transposases were highly conserved, mostly uniform in length and contained intact DNA-binding motifs and DD39D catalytic domains with few notable frameshift, indel and nonsense mutations. This suggested the MLEs are probably still mobile, not yet affected by vertical inactivation. A phylogenetic tree of the Bambusoideae subfamily established using ribosomal DNA internal transcribed spacer sequences was incongruent with a second tree based on the MLE transposase genes. This evidence, together with the presence of near-identical MLEs in distantly related species and diverse MLEs in closely related species, indicates that MLEs have evolved in a distinct manner, probably independently of speciation events in the subfamily. The evolution and diversity of MLE transposase genes in the Bambusoideae subfamily is discussed.

The catalytic domain of all eukaryotic cut-and-paste transposase superfamilies

Cut-and-paste DNA transposable elements are major components of eukaryotic genomes and are grouped into superfamilies (e.g., hAT, P) based on sequence similarity of the element-encoded transposase. The transposases from several superfamilies possess a protein domain containing an acidic amino acid triad (DDE or DDD) that catalyzes the "cut and paste" transposition reaction. However, it was unclear whether this domain was shared by the transposases from all superfamilies. Through multiple-alignment of transposase sequences from a diverse collection of previously identified and recently annotated elements from a wide range of organisms, we identified the putative DDE/D triad for all superfamilies. Furthermore, we identified additional highly conserved amino acid residues or motifs within the DDE/D domain that together form a "signature string" that is specific to each superfamily. These conserved residues or motifs were exploited as phylogenetic characters to infer evolutionary relationships among all superfamilies. The phylogenetic analysis revealed three major groups that were not previously discerned and led us to revise the classification of several currently recognized superfamilies. Taking the data together, this study suggests that all eukaryotic cut-and-paste transposable element superfamilies have a common evolutionary origin and establishes a phylogenetic framework for all future cut-and-paste transposase comparisons.