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.

Genome-Wide Analysis of Parent-of-Origin Allelic Expression in Endosperms of Brassicaceae Species, Brassica rapa.

Uniparental gene expression, observed in both animals and plants, is termed genomic imprinting. Genomic imprinting is a well-known epigenetic phenomenon regulated through epigenetic modifications such as DNA methylation and histone modifications. Recent genome-wide studies of endosperm transcription have revealed the rapid change of imprinted genes between species, suggesting the flexibility of this phenomenon. Although the functional significance and evolutionary trends of imprinted genes are still obscure, it can be clarified by inter-species comparisons. In this study, we analyzed the pattern of genomic imprinting in Brassica rapa, a species related to Arabidopsis thaliana. Compared with the ancient karyotype of A. thaliana and B. rapa, B. rapa has a triplicated genome. Many imprinted genes, beyond the estimated number previously reported in other species, were observed. Several imprinted genes have been conserved among species in Brassicaceae. We also observed rapid molecular evolution of imprinted genes compared to non-imprinted genes in B. rapa. Especially, imprinted gene overlapping between species showed more rapid molecular evolution and preferential expression in endosperms. It may imply that a small number of imprinted genes have retained functional roles among diverged species and have been the target of natural selection.

Characterization of a heat-activated retrotransposon in Vigna angularis.

In plants, several transposable elements are conserved across species. We found a homolog of ONSEN, which is a heat-activated retrotransposon originally isolated from Arabidopsis thaliana, in Vigna. The ONSEN-like elements (VaONS) were detected in all the analyzed Japanese accessions of Vigna angularis (adzuki bean) by Southern blot analysis. However, VaONS sequences were observed to be polymorphic in the different accessions. Interestingly, extrachromosomal DNA (ecDNA) was detected in some accessions of adzuki bean, indicating the conserved heat-activation of VaONS. Furthermore, we successfully induced retrotransposition of VaONS in adzuki plant regenerated through callus. Findings of our study should provide a new tool for molecular breeding of adzuki bean.

Exploring the History of Chloroplast Capture in Arabis Using Whole Chloroplast Genome Sequencing.

Chloroplast capture occurs when the chloroplast of one plant species is introgressed into another plant species. The phylogenies of nuclear and chloroplast markers from East Asian Arabis species are incongruent, which indicates hybrid origin and shows chloroplast capture. In the present study, the complete chloroplast genomes of A. hirsuta, A. nipponica, and A. flagellosa were sequenced in order to analyze their divergence and their relationships. The chloroplast genomes of A. nipponica and A. flagellosa were similar, which indicates chloroplast replacement. If hybridization causing chloroplast capture occurred once, divergence between recipient species would be lower than between donor species. However, the chloroplast genomes of species with possible hybrid origins, A. nipponica and A. stelleri, differ at similar levels to possible maternal donor species A. flagellosa, which suggests that multiple hybridization events have occurred in their respective histories. The mitochondrial genomes exhibited similar patterns, while A. nipponica and A. flagellosa were more similar to each other than to A. hirsuta. This suggests that the two organellar genomes were co-transferred during the hybridization history of the East Asian Arabis species.

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.

Bursts of retrotransposition reproduced in Arabidopsis.

Retrotransposons, which proliferate by reverse transcription of RNA intermediates, comprise a major portion of plant genomes. Plants often change the genome size and organization during evolution by rapid proliferation and deletion of long terminal repeat (LTR) retrotransposons. Precise transposon sequences throughout the Arabidopsis thaliana genome and the trans-acting mutations affecting epigenetic states make it an ideal model organism with which to study transposon dynamics. Here we report the mobilization of various families of endogenous A. thaliana LTR retrotransposons identified through genetic and genomic approaches with high-resolution genomic tiling arrays and mutants in the chromatin-remodelling gene DDM1 (DECREASE IN DNA METHYLATION 1). Using multiple lines of self-pollinated ddm1 mutant, we detected an increase in copy number, and verified this for various retrotransposons in a gypsy family (ATGP3) and copia families (ATCOPIA13, ATCOPIA21, ATCOPIA93), and also for a DNA transposon of a Mutator family, VANDAL21. A burst of retrotransposition occurred stochastically and independently for each element, suggesting an additional autocatalytic process. Furthermore, comparison of the identified LTR retrotransposons in related Arabidopsis species revealed that a lineage-specific burst of retrotransposition of these elements did indeed occur in natural Arabidopsis populations. The recent burst of retrotransposition in natural population is targeted to centromeric repeats, which is presumably less harmful than insertion into genes. The ddm1-induced retrotransposon proliferations and genome rearrangements mimic the transposon-mediated genome dynamics during evolution and provide experimental systems with which to investigate the controlling molecular factors directly.

Evolution and control of imprinted FWA genes in the genus Arabidopsis.

A central question in genomic imprinting is how a specific sequence is recognized as the target for epigenetic marking. In both mammals and plants, imprinted genes are often associated with tandem repeats and transposon-related sequences, but the role of these elements in epigenetic gene silencing remains elusive. FWA is an imprinted gene in Arabidopsis thaliana expressed specifically in the female gametophyte and endosperm. Tissue-specific and imprinted expression of FWA depends on DNA methylation in the FWA promoter, which is comprised of two direct repeats containing a sequence related to a SINE retroelement. Methylation of this element causes epigenetic silencing, but it is not known whether the methylation is targeted to the SINE-related sequence itself or the direct repeat structure is also necessary. Here we show that the repeat structure in the FWA promoter is highly diverse in species within the genus Arabidopsis. Four independent tandem repeat formation events were found in three closely related species. Another related species, A. halleri, did not have a tandem repeat in the FWA promoter. Unexpectedly, even in this species, FWA expression was imprinted and the FWA promoter was methylated. In addition, our expression analysis of FWA gene in vegetative tissues revealed high frequency of intra-specific variation in the expression level. In conclusion, we show that the tandem repeat structure is dispensable for the epigenetic silencing of the FWA gene. Rather, SINE-related sequence is sufficient for imprinting, vegetative silencing, and targeting of DNA methylation. Frequent independent tandem repeat formation events in the FWA promoter led us to propose that they may be a consequence, rather than cause, of the epigenetic control. The possible significance of epigenetic variation in reproductive strategies during evolution is also discussed.

High diversity due to balancing selection in the promoter region of the Medea gene in Arabidopsis lyrata.

Molecular imprinting is the differential expression and/or silencing of alleles according to their parent of origin [1, 2]. Conflicts between parents, or parents and offspring, should cause "arms races," with accelerated evolution of the genes involved in imprinting. This should be detectable in the evolution of imprinting genes' protein sequences and in the promoter regions of imprinted genes. Previous studies, however, found no evidence of more amino acid substitutions in imprinting genes [1, 3]. We have analyzed sequence diversity of the Arabidopsis lyrata Medea (MEA) gene and divergence from the A. thaliana sequence, including the first study of the promoter region. In A. thaliana, MEA is imprinted, with paternal alleles silenced in endosperm cells [4, 5], and also functions in the imprinting machinery [4, 6]; MEA protein binding at the MEA promoter region indicates self-regulated imprinting [7-9]. We find the same paternal MEA allele silencing in A. lyrata endosperm but no evidence for adaptive evolution in the coding region, whereas the 5' flanking region displays high diversity, with distinct haplotypes, suggesting balancing selection in the promoter region.

Author Correction: Genome defense against integrated organellar DNA fragments from plastids into plant nuclear genomes through DNA methylation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Epigenetic regulation of spurious transcription initiation in Arabidopsis.

In plants, epigenetic regulation is critical for silencing transposons and maintaining proper gene expression. However, its impact on the genome-wide transcription initiation landscape remains elusive. By conducting a genome-wide analysis of transcription start sites (TSSs) using cap analysis of gene expression (CAGE) sequencing, we show that thousands of TSSs are exclusively activated in various epigenetic mutants of Arabidopsis thaliana and referred to as cryptic TSSs. Many have not been identified in previous studies, of which up to 65% are contributed by transposons. They possess similar genetic features to regular TSSs and their activation is strongly associated with the ectopic recruitment of RNAPII machinery. The activation of cryptic TSSs significantly alters transcription of nearby TSSs, including those of genes important for development and stress responses. Our study, therefore, sheds light on the role of epigenetic regulation in maintaining proper gene functions in plants by suppressing transcription from cryptic TSSs.

High DNA sequence diversity in pericentromeric genes of the plant Arabidopsis lyrata.

Differences in neutral diversity at different loci are predicted to arise due to differences in mutation rates and from the "hitchhiking" effects of natural selection. Consistent with hitchhiking models, Drosophila melanogaster chromosome regions with very low recombination have unusually low nucleotide diversity. We compared levels of diversity from five pericentromeric regions with regions of normal recombination in Arabidopsis lyrata, an outcrossing close relative of the highly selfing A. thaliana. In contrast with the accepted theoretical prediction, and the pattern in Drosophila, we found generally high diversity in pericentromeric genes, which is consistent with the observation in A. thaliana. Our data rule out balancing selection in the pericentromeric regions, suggesting that hitchhiking is more strongly reducing diversity in the chromosome arms than the pericentromere regions.

Divergence of RNA editing among Arabidopsis species.

RNA editing altered the RNA sequence by replacing the C nucleotide to U in the organellar genomes of plants. RNA editing status sometimes differed among distant species. The pattern of conservation and variation of RNA editing status made it possible to evaluate evolutionary mechanisms impacting functional aspects of RNA editing. In this study, divergence of RNA editing in the chloroplast genome among Arabidopsis species was analyzed to determine 9 losses and 1 gain in RNA editing. All changes in A. thaliana lineage resulted from changes to the chloroplast genome sequence, whereas changes in the A. lyrata / halleri lineage were possibly due to exclusive changes in the nuclear editing factors. One loss of RNA editing in A. lyrata was caused by a deficiency in the PPR gene OTP80. The changes in RNA editing occurred approximately every two million years and were not observed at functionally important sites. These results highlight the conserved nature of RNA editing status suggesting the importance of RNA editing during evolution.

Genome defense against integrated organellar DNA fragments from plastids into plant nuclear genomes through DNA methylation.

Nuclear genomes are always faced with the modification of themselves by insertions and integrations of foreign DNAs and intrinsic parasites such as transposable elements. There is also substantial number of integrations from symbiotic organellar genomes to their host nuclear genomes. Such integration might have acted as a beneficial mutation during the evolution of symbiosis, while most of them have more or less deleterious effects on the stability of current genomes. Here we report the pattern of DNA substitution and methylation on organellar DNA fragments integrated from plastid into plant nuclear genomes. The genome analyses of 17 plants show homology-dependent DNA substitution bias. A certain number of these sequences are DNA methylated in the nuclear genome. The intensity of DNA methylation also decays according to the increase of relative evolutionary times after being integrated into nuclear genomes. The methylome data of epigenetic mutants shows that the DNA methylation of organellar DNA fragments in nuclear genomes are mainly dependent on the methylation maintenance machinery, while other mechanisms may also affect on the DNA methylation level. The DNA methylation on organellar DNA fragments may contribute to maintaining the genome stability and evolutionary dynamics of symbiotic organellar and their host's genomes.

DNA Methylation Diversification at the Integrated Organellar DNA-Like Sequence.

Plants have a lot of diversity in epigenetic modifications such as DNA methylation in their natural populations or cultivars. Although many studies observing the epigenetic diversity within and among species have been reported, the mechanisms how these variations are generated are still not clear. In addition to the de novo spontaneous epi-mutation, the intra- and inter-specific crossing can also cause a change of epigenetic modifications in their progenies. Here we report an example of diversification of DNA methylation by crossing and succeeding selfing. We traced the inheritance pattern of epigenetic modification during the crossing experiment between two natural strains Columbia (Col), and Landsberg electa (Ler) in model plant Arabidopsis thaliana to observe the inheritance of DNA methylation in two organellar DNA-like sequence regions in the nuclear genome. Because organellar DNA integration to the nuclear genome is common in flowering plants and these sequences are occasionally methylated, such DNA could be the novel source of plant genome evolution. The amplicon sequencing, using bisulfite-converted DNA and a next-generation auto-sequencer, was able to efficiently track the heredity of DNA methylation in F1 and F2 populations. One region showed hypomethylation in the F1 population and succeeding elevation of DNA methylation with large variance in the F2 population. The methylation level of Col and Ler alleles in F2 heterozygotes showed a significant positive correlation, implying the trans-chromosomal effect on DNA methylation. The results may suggest the possible mechanism causing the natural epigenetic diversity within plant populations.

The complete chloroplast genome of Sisymbrium irio.

The complete chloroplast genome of Sisymbrium irio was determined. The length of the complete chloroplast genome is 154,001 bp. The whole chloroplast genome consists of 83,891 bp long single copy (LSC) and 17,630 bp small single copy (SSC) regions, separated by a pair of 26,240 bp inverted repeat (IR) regions. The S. irio chloroplast genome encodes 112 annotated known unique genes including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The phylogenetic position of S. irio is sister to Brassiceae and Thlaspideae.

Duplication of centromeric histone H3 (HTR12) gene in Arabidopsis halleri and A. lyrata, plant species with multiple centromeric satellite sequences.

Arabidopsis halleri and lyrata have three different major centromeric satellite sequences, a unique finding for a diploid Arabidopsis species. Since centromeric histones coevolve with centromeric satellites, these proteins would be predicted to show signs of selection when new centromere satellites have recently arisen. We isolated centromeric protein genes from A. halleri and lyrata and found that one of them, HTR12 (CENP-A), is duplicated, while CENP-C is not. Phylogenetic analysis indicates that the HTR12 duplication occurred after these species diverged from A. thaliana. Genetic mapping shows that HTR12 copy B has the same genomic location as the A. thaliana gene; the other copy (A, at the other end of the same chromosome) is probably the new copy. To test for selection since the duplication, we surveyed diversity at both HTR12 loci within A. lyrata. Overall, there is no strong evidence for an "evolutionary arms race" causing multiple replacement substitutions. The A. lyrata HTR12B sequences fall into three classes of haplotypes, apparently maintained for a long time, but they all encode the same amino acid sequence. In contrast, HTR12A has low diversity, but many variants are amino acid replacements, possibly due to independent selective sweeps within populations of the species.

Centromere locations and associated chromosome rearrangements in Arabidopsis lyrata and A. thaliana.

We analyzed linkage and chromosomal positions of genes in A. lyrata ssp. petraea that are located near the centromere (CEN) regions of A. thaliana, using at least two genes from the short and long arms of each chromosome. In our map, genes from all 10 A. thaliana chromosome arms are also tightly linked in A. lyrata. Genes from the regions on the two sides of CEN5 have distant map localizations in A. lyrata (genes on the A. thaliana short-arm genes are on linkage group AL6, and long-arm genes are on AL7), but genes from the other four A. thaliana centromere regions remain closely linked in A. lyrata. The observation of complete linkage between short- and long-arm centromere genes, but not between genes in other genome regions that are separated by similar physical distances, suggests that crossing-over frequencies near the A. lyrata ssp. petraea centromere regions are low, as in A. thaliana. Thus, the centromere positions appear to be conserved between A. thaliana and A. lyrata, even though three centromeres have been lost in A. thaliana, and the core satellite sequences in the two species are very different. We can now definitively identify the three centromeres that were eliminated in the fusions that formed the A. thaliana chromosomes. However, we cannot tell whether genes were lost along with these centromeres, because such genes are absent from the A. thaliana genome, which is the sole source of markers for our mapping.

Testing for effects of recombination rate on nucleotide diversity in natural populations of Arabidopsis lyrata.

We investigated DNA sequence diversity for loci on chromosomes 1 and 2 in six natural populations of Arabidopsis lyrata and tested for the role of natural selection in structuring genomewide patterns of variability, specifically examining the effects of recombination rate on levels of silent polymorphism. In contrast with theoretical predictions from models of genetic hitchhiking, maximum-likelihood-based analyses of diversity and divergence do not suggest reduction of diversity in the region of suppressed recombination near the centromere of chromosome 1, except in a single population from Russia, in which the pericentromeric region may have undergone a local selective sweep or demographic process that reduced variability. We discuss various possibilities that might explain why nucleotide diversity in most A. lyrata populations is not related to recombination rate, including genic recombination hotspots, and low gene density in the low recombination rate region.

Analysis of nuclear mitochondrial DNAs and factors affecting patterns of integration in plant species.

Sequences homologous to organellar DNA that have been integrated into nuclear genomes are referred to as nuclear mitochondrial DNAs (NUMTs) and nuclear plastid DNAs (NUPTs). NUMTs in nine plant species were analyzed to reveal the integration patterns and possible factors involved. The cumulative lengths of NUMTs in two-thirds of species analyzed were greater than those of NUPTs observed in a previous study. The age distribution of NUMTs was similar to that of NUPTs, suggesting similar mechanisms for integration and degradation of both NUPTs and NUMTs. Nuclear genome size and the cumulative length of NUMTs showed a significant positive correlation for older but not younger NUMTs. The same correlation was also found between nuclear genome size and older NUPTs in 17 species. These results suggested that genome size is a key factor to determine the cumulative length of relatively older NUPTs/NUMTs. Although the factor(s) determining the cumulative length of younger NUPTs/NUMTs is unclear, these sequences may be more deleterious, which could explain the different manner of determining the cumulative length of younger NUPTs/NUMTs in nuclear genomes. In addition, a relationship between the cumulative length of integrated NUMTs and complexity of mitochondrial genomes (i.e., the number of repeats) was found. The results indicate that the structural complexity of both NUMTs and their original mitochondrial sequences affects integration and degradation processes.

Characterization of a heat-activated retrotransposon in natural accessions of Arabidopsis thaliana.

Natural accessions are used for studying intraspecies genetic variation in the model plant Arabidopsis thaliana in order to address fundamental questions of evolution. Transposable elements are responsible for a wide range of mutations and play significant roles in shaping a genome over evolutionary time. In the present study, we aimed to characterize ONSEN, a heat-activated long terminal repeat (LTR) retrotransposon, in natural A. thaliana accessions. Southern blot analysis demonstrated that ONSEN was present in all the studied accessions, but the copy number was diverse. Olympia-1 contained a single ONSEN copy, located in the centromere of Chromosome 3. A premature stop codon in Olympia-1 ONSEN presumably abolishes integrase activity, which in turn presumably renders the retrotransposon non-functional. Hybridization of Col-0 with Olympia-1 showed that several ONSEN copies in Col-0 were activated by heat stress and maintained their transpositional activity in the progeny.