Inhibition of RNA polymerase II allows controlled mobilisation of retrotransposons for plant breeding
© The Author(s). 2017
Received: 13 February 2017
Accepted: 27 June 2017
Published: 7 July 2017
Retrotransposons play a central role in plant evolution and could be a powerful endogenous source of genetic and epigenetic variability for crop breeding. To ensure genome integrity several silencing mechanisms have evolved to repress retrotransposon mobility. Even though retrotransposons fully depend on transcriptional activity of the host RNA polymerase II (Pol II) for their mobility, it was so far unclear whether Pol II is directly involved in repressing their activity.
Here we show that plants defective in Pol II activity lose DNA methylation at repeat sequences and produce more extrachromosomal retrotransposon DNA upon stress in Arabidopsis and rice. We demonstrate that combined inhibition of both DNA methylation and Pol II activity leads to a strong stress-dependent mobilization of the heat responsive ONSEN retrotransposon in Arabidopsis seedlings. The progenies of these treated plants contain up to 75 new ONSEN insertions in their genome which are stably inherited over three generations of selfing. Repeated application of heat stress in progeny plants containing increased numbers of ONSEN copies does not result in increased activation of this transposon compared to control lines. Progenies with additional ONSEN copies show a broad panel of environment-dependent phenotypic diversity.
We demonstrate that Pol II acts at the root of transposon silencing. This is important because it suggests that Pol II can regulate the speed of plant evolution by fine-tuning the amplitude of transposon mobility. Our findings show that it is now possible to study induced transposon bursts in plants and unlock their use to induce epigenetic and genetic diversity for crop breeding.
KeywordsEpigenetics DNA methylation Genome integrity Evolution Oryza sativa Arabidopsis thaliana
Like retroviruses, long terminal repeat (LTR) retrotransposons (class I elements), which represent the most abundant class of transposable elements (TEs) in eukaryotes, transpose via a copy and paste mechanism. This process requires the conversion of a full length RNA polymerase II (Pol II) transcript into extrachromosomal complementary DNA (ecDNA) by reverse transcription . In their life cycle LTR retrotransposons can produce extrachromosomal circular DNA (eccDNA), which is an indicator for their ongoing activity . In plants, TEs are increasingly seen as a source of genetic and epigenetic variability and thus important drivers of evolution [3–6]. However, plants have evolved several regulatory pathways to retain control over the activity of these potentially harmful mobile genetic elements. Cytosine methylation (mC) plays a central role in TE silencing in plants . In addition, plants have evolved two Pol II-related RNA polymerases, Pol IV and Pol V, that are essential to provide specific silencing signals leading to RNA-directed DNA methylation (RdDM) at TEs , thereby limiting their mobility [9–11]. More recently, various additional non-canonical Pol IV-independent RdDM pathways have been described . Notably it was found that Pol II itself also plays an important role in RdDM [13, 14] by feeding template RNAs into downstream factors such as RNA-DEPENDENT RNA POLYMERASE 6 (RDR6), resulting in dicer-dependent or -independent initiation and establishment of TE-specific DNA methylation . Beyond that, recent work suggests a new “non-canonical” branch of RdDM that specializes in targeting transcriptionally active full-length TEs . This pathway functions independently of RDRs via Pol II transcripts that are directly processed by DCL3 into small interfering RNAs (siRNAs).
To better understand the mechanisms by which the drugs enhanced the activation of ONSEN after HS at the DNA level, we assessed how they influenced DNA methylation at the genome-wide level using whole-genome bisulfite sequencing (WGBS) after CS. Overall, we found that all drug treatments affected global DNA methylation levels. While the treatment with Z affected all sequence contexts, we observed that inhibition of Pol II primarily affected cytosine methylation in the CHG and CHH sequence contexts (where H is an A, T or G). The combined A&Z treatment had a slight additive de-methylating effect in the CHG and CHH contexts compared to A or Z alone (Fig. 2a, b). DNA methylation levels at one ONSEN locus (AT1TE12295) is depicted in Fig 2c. Treatment with A led to a slight decrease in DNA methylation, which was more apparent in Z- and A&Z-treated plants. We then checked by northern blot whether the degree of reduction in DNA methylation would coincide with increased ONSEN transcript levels directly after HS. We found that treatment with Z alone resulted in the highest ONSEN transcript level after HS (Fig. 2d). Considering the data obtained on ONSEN ecDNA (Fig. 1c), we concluded that a substantial proportion of these Z-induced transcripts were not suitable templates for ONSEN ecDNA synthesis.
In Drosophila, it has been shown that Pol II-mediated antisense transcription results in the production of TE-derived siRNAs in a Dicer-2-dependent manner . In support of this in Arabidopsis, a recent publication pointed out the importance of DCL3 in regulating ONSEN in the ddm1 background . To elucidate whether the effect of Pol II inhibition was also dicer-dependent, we grew both rdr6 and dcl2/3/4 triple mutant plants on A, applied HS and measured ONSEN ecDNA levels. Strikingly, we found that A still enhanced ecDNA accumulation in rdr6 plants, whereas inhibition of Pol II had no additional effect in the dcl2/3/4 triple mutant (Fig. 2e).
TE insertions can interrupt genes or alter their expression by recruiting epigenetic marks or by stress-dependent readout transcription from the 3′ LTR into flanking regions . To test this, we grew the S2 generation of the selected hc lines under long- and short-day conditions. Interestingly, we observed that many hc lines showed clear and homogenous phenotypes in response to the different growth conditions (plant size, chlorophyll content and flowering time; Fig. 3c, d).
To demonstrate that ONSEN insertions could directly influence such developmental phenotypes, we closely investigated hc line 3, which produced white seeds (Fig. 4a). Using a candidate gene approach, we found that an ONSEN insertion in transparent testa 6 (TT6, AT3G51240; Fig. 4b) was responsible for the recessive white seed phenotype [23, 24]. This was confirmed by segregation analysis of the F2 generation of a cross between WT and hc line 3 (Fig. 4a) followed by genotyping (Fig. 4c).
In this study, we show the importance of Pol II in the repression of TE mobility in plants. By choosing the well-characterized heat inducible ONSEN retrotransposon, we were able to specifically address the role of Pol II in silencing transcriptionally active endogenous TEs in WT plants. Recent studies propose Pol II as the primary source for the production of TE-silencing signals that can then feed into the RNA silencing and DNA methylation pathways . Our data strongly support these findings at two levels. First, we found that inhibition of Pol II activity reduced the degree of DNA methylation at ONSEN, demonstrating its distinct role in this process, and that Pol II also contributes to reinforcing silencing at the genome-wide level, primarily in the CHH but also in the CHG context. Second, our finding that DCL enzymes are sufficient to process the silencing signal produced by Pol II suggest that Pol II acts at very early steps in the TE silencing pathway by providing substrates to these enzymes. The observation that inhibition of Pol II in the rdr6 background still further enhanced ONSEN accumulation after HS supports the notion that Pol II plays a central role in the previously proposed expression-dependent RdDM pathway .
Using mobilome sequencing we confirmed previous findings  that this approach is a powerful diagnostic tool to detect mobile retrotransposons: we detected highest levels of eccDNA of ONSEN in HS and drug-treated Arabidopsis seedlings and found new insertions in successive generations of these plants. Using the same approach on rice we were able to detect production of Houba eccDNA after drug treatments, suggesting that the progeny will then contain novel Houba insertions. This is still to be confirmed and may be hampered by the already very high Houba copy number present in the genome .
Our findings may indicate that Pol II is primarily involved in silencing young, recently active retrotransposons and perhaps to a lesser extent other tightly silenced TEs. Indeed, there are indications of very recent natural transposition events for ONSEN  and Houba  in the Arabidopsis and rice genomes, respectively. For instance, the annual temperature range has and may still contribute to contrasting ONSEN mobilization events in different Arabidopsis accessions . Houba is the most abundant TE of the copia family in rice and has been active in the last 500,000 years .
Overall, our findings lead to the question of when plants lower their guard: under what conditions could Pol II be less effective in silencing TEs? Certain stresses that affect the cell cycle have been reported to lead to the inactivation of Pol II [31, 32]; this would provide a window of opportunity for TEs to be mobilized. Therefore, combined stresses that affect the cell cycle and activate TEs may lead to actual TE bursts under natural growth conditions. Interestingly, it has been reported that retrotransposon-derived short interspersed element (SINE) transcripts can inhibit Pol II activity . This strongly suggests the presence of an ongoing arms race between retrotransposons and Pol II. Considering that almost all organisms analyzed so far have TEs  and RNA polymerases  and the reliance of TEs on host RNA polymerases, it may—from an evolutionary point of view—not come as a surprise that Pol II also has a function as an important regulator of retrotransposon activity. Strikingly, it has been shown in both Saccharomyces cerevisiae and Drosophila melanogaster that Pol II-dependent intra-element antisense transcription plays an important role in TE silencing [21, 35]. In addition, we observed a discrepancy in ONSEN transcript accumulation and measured ecDNA after HS in seedlings that were treated with zebularine only. This substantiates the notion that both the quantity and quality of transcripts affect regulation, reverse transcription and successful integration of retrotransposons. This is well in line with previous observations demonstrating that different TE-derived transcripts have distinct functions in the regulation of TE activity . As a next step it will be of great interest to investigate if Pol II-dependent antisense transcription of TEs and subsequent dicer-dependent processing may be the key to solve “the chicken and the egg problem” of de novo silencing functional retrotransposons in eukaryotes.
Finally, our findings will allow future studies on the potential beneficial role TEs play in adaptation to stresses. Indeed, two recent studies point out the adaptive potential of retrotraonsposon and, more specifically, ONSEN copy number variation in natural accessions  and RdDM mutant backgrounds of Arabidopsis . Upon mobilization, the heat-response elements in the LTRs of ONSEN  can create new gene regulatory networks responding to heat stress . Therefore, it will now be of great interest to test if the ONSEN hc lines obtained in this study are better adapted to heat stress. This will allow us to test if retrotransposon-induced genetic and epigenetic changes more rapidly create beneficial alleles than would occur by random mutagenesis. Furthermore, the observation that HS did not lead to a stronger activation of ONSEN in hc lines compared to WT plants suggests that genome stability is not compromised in these lines. This result can be explained by at least two possible mechanisms: (i) the occurrence of insertions of inverted duplications of ONSEN, such as has been observed for the Mu killer locus in maize —such insertions will lead to the production of double-stranded RNA feeding into gene silencing and thereby limit the activity of that TE; and (ii) balancing of TE activity and integrated copy number as has been described for EVADE in Arabidopsis . In this case, when a certain TE copy number threshold is reached robust transcriptional gene silencing takes over, thereby limiting TE mobility and ensuring genome stability. The stability of new TE insertions is an important aspect in light of the future use of TEs in crop breeding and trait stability.
TEs are important contributors to genome evolution. The ability to mobilize them in plants and possibly in other eukaryotes in a controlled manner with straightforward drug application, as shown here, opens the possibility to study their importance in inducing genetic and epigenetic changes resulting from external stimuli. Because the induced transposition of ONSEN can efficiently produce developmental changes in Arabidopsis, it will be very interesting to test if specific stress-induced TE activation can be used for directed crop breeding for better stress tolerance in the near future.
Prior to germination, Arabidopsis seeds were stratified for 2 days at 4 °C. Before and during stress treatments plants were grown under controlled conditions in a Sanyo MLR-350 growth chamber on solid ½ MS medium (1% sucrose, 0.5% Phytagel (Sigma), pH 5.8) under long day conditions (16 h light) at 24 °C (day) and 22 °C (night) (Arabidopsis) and 12 h at 28 °C (day) and 27 °C (night) (O. sativa).
To analyze successive generations, seedlings were transferred to soil and grown under long day conditions (16 h light) at 24 °C (day) and 22 °C (night) (Arabidopsis) in a Sanyo MLR-350 growth chamber until seed maturity.
For phenotyping, Arabidopsis plants were grown under long day conditions (16 h light) at 24 °C (day) and 22 °C (night) and short day conditions (10 h light) at 21 °C (day) and 18 °C (night).
Stress and chemical treatments
Surface sterilized seeds of Arabidopsis and O. sativa were germinated and grown on solid ½ MS medium that was supplemented with sterile filtered zebularine (Sigma; stock, 5 mg/ml in DMSO), α-amanitin (Sigma; stock, 1 mg/ml in water) or a combination of both chemicals. Control stresses (6 °C for 24 h followed by control conditions for 24 h, CS) and heat stresses (6 °C for 24 h followed by 37 °C for 24 h, HS) of Arabidopsis seedlings were conducted as described previously .
For qPCR and prior to digestions, total DNA from Arabidopsis plants was extracted with the DNeasy Plant Mini Kit (Qiagen) following the manufacturer’s recommendations. For the qPCRs to measure the ONSEN copy number following HS and chemical treatments the aerial parts of at least ten Arabidopsis plants per replicate were pooled prior to DNA extraction. To track ONSEN copy numbers in the S1–3 generations of controls (only HS) and hc lines (HS + A&Z treatment) DNA from true leaves was extracted. For the estimation of the ONSEN transposition frequency, total DNA of pools consisting of at least eight seedlings of the progeny of HS + A&Z-treated plants was isolated. The DNA concentration was measured with a Qubit Fluorometer (Thermo Fisher Scientific). The copy numbers of ONSEN were determined with qPCRs on total DNA using a TaqMan master mix (Life Technologies) in a final volume of 10 μl in the Light-Cycler 480 (Roche). ACTIN2 (AT3G18780) was used to normalize DNA levels. Primer sequences are given in Additional file 1: Table S1.
For the mobilome-seq analysis total DNA from the pooled aerial parts of three 10-day-old O. sativa seedlings was extracted as previously reported . Genomic DNA (5 μg) for each sample was purified using a Geneclean kit (MPBio, USA) according to the manufacturer’s instructions. ecDNA was isolated from the GeneClean product using PlasmidSafe DNase (Epicentre, USA) according to the manufacturer’s instructions, except that the 37 °C incubation was performed for 17 h. DNA samples were precipitated by adding 0.1 volume of 3 M sodium acetate (pH 5.2), 2.5 volumes of ethanol and 1 μl of glycogen (Fisher, USA) and incubating overnight at −20 °C. The precipitated circular DNA was amplified by random rolling circle amplification using the Illustra TempliPhi kit (GE Healthcare, USA) according to the manufacturer’s instructions except that the incubation was performed for 65 h at 28 °C. The DNA concentration was determined using the DNA PicoGreen kit (Invitrogen, USA) using a LightCycler480 (Roche, USA). One nanogram of amplified ecDNA from each sample was used to prepare the libraries using the Nextera XT library kit (Illumina, USA) according to the manufacturer’s instructions. DNA quality and concentration were determined using a high sensitivity DNA Bioanalyzer chip (Agilent Technologies, USA). Samples were pooled and loaded onto a MiSeq platform (Illumina, USA) and 2 × 250-nucleotide paired-end sequencing was performed. Quality control of FASTQ files was done using the FastQC tool (version 0.10.1). To remove any read originating from organelle circular genomes, reads were mapped against the mitochondria and chloroplast genomes using the program Bowtie2 version 2.2.2 71 with --sensitive local mapping. Unmapped reads were mapped against the reference genome IRGSP1.0 (http://rgp.dna.affrc.go.jp/E/IRGSP/Build5/build5.html) using the following parameters: --sensitive local, -k 1. DNA from both mitochondria and chloroplast genomes integrated in nuclear genomes was masked (1,697,400 bp). The TE-containing regions cover 194,224,800 bp in O. sativa. Finally, the bam alignment files were normalized and compared using deeptools  and visualized with the Integrative Genomics Viewer (IGV) software (https://www.broadinstitute.org/igv/). Data from the mobilome analysis were submitted to GEO (accession number GSE90484).
The presence of circular Houba copies was tested by an inverse PCR on 7 ng of the rolling-circle amplified template that was also used for sequencing. A PCR specific to chloroplast DNA served as a loading control. PCR products were separated on a 1% agarose gel that was stained with a Midori Green Nucleic Acid Staining Solution (Nippon Genetics Europe). Primer sequences are given in Additional file 1: Table S1.
The integration of additional copies of ONSEN into the genome of heat stressed and treated plants was ascertained by a simplified transposon display based on the GenomeWalker Universal kit (Clontech Laboratories), as previously described  with the following modifications: 300 ng of total DNA from adult plants in the S2 generation of heat stressed and A&Z-treated plants was extracted with a DNeasy Plant Mini Kit (QIAGEN) and digested with blunt cutter restriction enzyme DraI (NEB). After purification with a High Pure PCR Product Purification Kit (Roche) digested DNA was ligated to the annealed GenWalkAdapters 1&2. The PCR was performed with the adaptor-specific primer AP1 and the ONSEN-specific primer Copia78 3′ LTR. The PCR products were separated on a 2% agarose gel that was stained with Midori Green. For primer sequence information, see Additional file 1: Table S1.
Cloning, sequencing and genotyping of new insertions
To identify the genomic region of new ONSEN insertions, the PCR product of the transposon display was purified using a High Pure PCR Product Purification Kit (Roche), ligated into a pGEM-T vector (Promega) and transformed into Escherichia coli. After a blue white selection, positive clones were used for the insert amplification and sequencing (StarSEQ). The obtained sequences were analyzed with Geneious 8.2.1 and blasted against the Arabidopsis reference genome. The standard genotyping PCRs to prove novel ONSEN insertions were performed with combinations of the ONSEN-specific primer Copia78 3′ LTR and primers listed in Additional file 1: Table S1.
RNA analysis and northern blotting
Total RNA from the aerial part of at least ten Arabidopsis seedlings was isolated using the TRI Reagent (Sigma) according to the manufacturer’s recommendations. RNA concentration was measured (Qubit RNA HS Assay Kit, Thermo Fisher) and 15 μg of RNA was separated on a denaturing 1.5% agarose gel, blotted on a Hybond-N+ (GE Healthcare) membrane and hybridized with 25 ng of a gel-purified and P32-labelled probe (Megaprime DNA Labelling System, GE Healthcare) specific to the full length ONSEN transcript (see Additional file 1: Table S1 for primer sequences). Northern blots were repeated in three independent experiments with the same results.
Whole-genome DNA methylation analysis
Whole-genome bisulfite sequencing library preparation and DNA conversion were performed as previously reported . Bisulphite read mapping and methylation value extraction were done on the Arabidopsis TAIR10 genome sequence using BSMAP v2.89 . Following mapping of the reads the fold coverages of the genome for CS, CS + A, CS + Z and CS + A&Z were 13.4, 13.2, 18.4 and 16.3, respectively. Data from the bisulphite sequencing analysis have been submitted to GEO (accession number GSE99396).
Statistical analyses were performed with SigmaPlot (v. 11.0). Depending on the normality of the data, either an H-test or a one-way ANOVA was performed. The Student-Newman-Keuls method was used for multiple comparisons.
We wish to thank Emilija Hristova for her support at the beginning of this project. We thank Christel Llauro for technical support on the production of the mobilomes and Todd Blevins for providing the dcl2/3/4 triple mutant line. We thank the IMAC platform from the Structure Fédérative de Recherche ‘Qualité et Santé du Végétal’ (SFR QUASAV) for their technical support (Illumina sequencing).
This work was supported by grants provided by the European Commission (PITN-GA-2013-608422–IDP BRIDGES to MT and ERC grant 725701 BUNGEE to EB) and the region of Pays de la Loire (ConnecTalent EPICENTER project awarded to EB).
Availability of data and materials
The mobilome sequencing data and whole-genome DNA methylation analysis data are available in the GEO (accession numbers GSE90484 and GSE99396, respectively).
MT and EB conceived the study. MT, SL, SB and MM performed experiments. ND performed methylome analyses. MT and EB wrote the paper with contributions from SL and MM. All authors read and approved the final manuscript.
Ethics approval and consent to participate
MT and EB declare that a patent application based on the presented discoveries has been submitted to the European Patent Office (PCT/EP2016/079276). EB is CEO of epibreed Ltd, a company that has an exclusive use license for this patent.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Schulman AH. Retrotransposon replication in plants. Curr Opin Virol. 2013;3:604–14.View ArticlePubMedGoogle Scholar
- Lanciano S, Carpentier M-C, Llauro C, Jobet E, Robakowska-Hyzorek D, Lasserre E, et al. Sequencing the extrachromosomal circular mobilome reveals retrotransposon activity in plants. PLoS Genet. 2017;13:e1006630.View ArticlePubMedPubMed CentralGoogle Scholar
- Paszkowski J. Controlled activation of retrotransposition for plant breeding. Curr Opin Biotechnol. 2015;32C:200–6.View ArticleGoogle Scholar
- Huang CRL, Burns KH, Boeke JD. Active transposition in genomes. Annu Rev Genet. 2012;46:651–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Belyayev A. Bursts of transposable elements as an evolutionary driving force. J Evol Biol. 2014;27:2573–84.View ArticlePubMedGoogle Scholar
- Lisch D. How important are transposons for plant evolution? Nat Rev Genet. 2013;14:49–61.View ArticlePubMedGoogle Scholar
- Miura A, Yonebayashi S, Watanabe K, Toyama T, Shimada H, Kakutani T. Mobilization of transposons by a mutation abolishing full DNA methylation in Arabidopsis. Nature. 2001;411:212–4.View ArticlePubMedGoogle Scholar
- Matzke MA, Mosher RA. RNA-directed DNA methylation: an epigenetic pathway of increasing complexity. Nat Rev Genet. 2014;15:394–408.View ArticlePubMedGoogle Scholar
- Tsukahara S, Kobayashi A, Kawabe A, Mathieu O, Miura A, Kakutani T. Bursts of retrotransposition reproduced in Arabidopsis. Nature. 2009;461:423–U125.View ArticlePubMedGoogle Scholar
- Mirouze M, Reinders J, Bucher E, Nishimura T, Schneeberger K, Ossowski S, et al. Selective epigenetic control of retrotransposition in Arabidopsis. Nature. 2009;461:427–30.View ArticlePubMedGoogle Scholar
- Ito H, Gaubert H, Bucher E, Mirouze M, Vaillant I, Paszkowski J. An siRNA pathway prevents transgenerational retrotransposition in plants subjected to stress. Nature. 2011;472:115–9.View ArticlePubMedGoogle Scholar
- Matzke MA, Kanno T, Matzke AJM. RNA-directed DNA methylation: the evolution of a complex epigenetic pathway in flowering plants. Annu Rev Plant Biol. 2015;66:243–67.View ArticlePubMedGoogle Scholar
- Gao Z, Liu H-L, Daxinger L, Pontes O, He X, Qian W, et al. An RNA polymerase II- and AGO4-associated protein acts in RNA-directed DNA methylation. Nature. 2010;465:106–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Zheng B, Wang Z, Li S, Yu B, Liu J-Y, Chen X. Intergenic transcription by RNA Polymerase II coordinates Pol IV and Pol V in siRNA-directed transcriptional gene silencing in Arabidopsis. Genes Dev. 2009;23:2850–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Cuerda-Gil D, Slotkin RK. Non-canonical RNA-directed DNA methylation. Nature Plants. 2016;2:16163.View ArticlePubMedGoogle Scholar
- Panda K, Ji L, Neumann DA, Daron J, Schmitz RJ, Slotkin RK. Full-length autonomous transposable elements are preferentially targeted by expression-dependent forms of RNA-directed DNA methylation. Genome Biol. 2016;17:1–19.View ArticleGoogle Scholar
- Lindell TJ, Weinberg F, Morris PW, Roeder RG, Rutter WJ. Specific inhibition of nuclear RNA polymerase II by alpha-Amanitin. Science. 1970;170:447–9.View ArticlePubMedGoogle Scholar
- Haag JR, Ream TS, Marasco M, Nicora CD, Norbeck AD, Pasa-Tolić L, et al. In vitro transcription activities of Pol IV, Pol V, and RDR2 reveal coupling of Pol IV and RDR2 for dsRNA synthesis in plant RNA silencing. Mol Cell. 2012;48:811–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Baubec T, Pecinka A, Rozhon W, Mittelsten SO. Effective, homogeneous and transient interference with cytosine methylation in plant genomic DNA by zebularine. Plant J. 2009;57:542–54.View ArticlePubMedGoogle Scholar
- Flavell AJ, Ish-Horowicz D. Extrachromosomal circular copies of the eukaryotic transposable element Copia in cultured Drosophila cells. Nature. 1981;292:591–5.View ArticlePubMedGoogle Scholar
- Russo J, Harrington AW, Steiniger M. Antisense transcription of retrotransposons in Drosophila: an origin of endogenous small interfering RNA precursors. Genetics. 2016;202:107–21.View ArticlePubMedGoogle Scholar
- Matsunaga W, Ohama N, Tanabe N, Masuta Y, Masuda S, Mitani N, et al. A small RNA mediated regulation of a stress-activated retrotransposon and the tissue specific transposition during the reproductive period in Arabidopsis. Front Plant Sci. 2015;6:48.View ArticlePubMedPubMed CentralGoogle Scholar
- Rosso MG, Li Y, Strizhov N, Reiss B, Dekker K, Weisshaar B. An Arabidopsis thaliana T-DNA mutagenized population (GABI-Kat) for flanking sequence tag-based reverse genetics. Plant Mol Biol. 2003;53:247–59.View ArticlePubMedGoogle Scholar
- Appelhagen I, Thiedig K, Nordholt N, Schmidt N, Huep G, Sagasser M, et al. Update on transparent testa mutants from Arabidopsis thaliana: characterisation of new alleles from an isogenic collection. Planta. 2014;240:955–70.View ArticlePubMedGoogle Scholar
- Kawahara Y, la Bastide de M, Hamilton JP, Kanamori H, Mccombie WR, Ouyang S. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice. 2013;6:4–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Panaud O, Vitte C, Hivert J, Muziak S, Talag J, Brar D, et al. Characterization of transposable elements in the genome of rice (Oryza sativa L.) using Representational Difference Analysis (RDA). Mol Genet Genomics. 2002;268:113–21.View ArticlePubMedGoogle Scholar
- Vitte C, Panaud O. LTR retrotransposons and flowering plant genome size: emergence of the increase/decrease model. Cytogenet Genome Res. 2005;110:91–107.View ArticlePubMedGoogle Scholar
- Quadrana L, Silveira AB, Mayhew GF, LeBlanc C, Martienssen RA, Jeddeloh JA, et al. The Arabidopsis thaliana mobilome and its impact at the species level. elife. 2016;5:e15716.View ArticlePubMedPubMed CentralGoogle Scholar
- Vitte C, Panaud O, Quesneville H. LTR retrotransposons in rice (Oryza sativa, L.): recent burst amplifications followed by rapid DNA loss. BMC Genomics. 2007;8:218–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Wicker T, Keller B. Genome-wide comparative analysis of copia retrotransposons in Triticeae, rice, and Arabidopsis reveals conserved ancient evolutionary lineages and distinct dynamics of individual copia families. Genome Res. 2007;17:1072–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Oelgeschläger T. Regulation of RNA polymerase II activity by CTD phosphorylation and cell cycle control. J Cell Physiol. 2001;190:160–9.View ArticleGoogle Scholar
- Palancade B, Bensaude O. Investigating RNA polymerase II carboxyl-terminal domain (CTD) phosphorylation. Eur J Biochem. 2003;270:3859–70.View ArticlePubMedGoogle Scholar
- Pai DA, Kaplan CD, Kweon HK, Murakami K, Andrews PC, Engelke DR. RNAs nonspecifically inhibit RNA polymerase II by preventing binding to the DNA template. RNA. 2014;20:644–55.View ArticlePubMedPubMed CentralGoogle Scholar
- Lazcano A, Fastag J, Gariglio P, Ramirez C, Oro J. On the early evolution of RNA-polymerase. J Mol Evol. 1988;27:365–76.View ArticlePubMedGoogle Scholar
- Berretta J, Pinskaya M, Morillon A. A cryptic unstable transcript mediates transcriptional trans-silencing of the Ty1 retrotransposon in S. cerevisiae. Genes Dev. 2008;22:615–26.View ArticlePubMedPubMed CentralGoogle Scholar
- Chang W, Jääskeläinen M, Li S-P, Schulman AH. BARE retrotransposons are translated and replicated via distinct RNA pools. PLoS One. 2013;8:e72270–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Ito H, Kim J-M, Matsunaga W, Saze H, Matsui A, Endo TA, et al. A Stress-activated transposon in Arabidopsis induces transgenerational abscisic acid insensitivity. Sci Rep. 2016;6:23181.View ArticlePubMedPubMed CentralGoogle Scholar
- Pietzenuk B, Markus C, Gaubert H, Bagwan N, Merotto A, Bucher E, et al. Recurrent evolution of heat-responsiveness in Brassicaceae COPIA elements. Genome Biol. 2016;17:209.View ArticlePubMedPubMed CentralGoogle Scholar
- Slotkin R, Freeling M, Lisch D. Heritable transposon silencing initiated by a naturally occurring transposon inverted duplication. Nat Genet. 2005;37:641–4.View ArticlePubMedGoogle Scholar
- Marí-Ordóñez A, Marchais A, Etcheverry M, Martin A. Reconstructing de novo silencing of an active plant retrotransposon. Nat Genet. 2013;45:1029–39.View ArticlePubMedGoogle Scholar
- Herr AJ, Jensen MB, Dalmay T, Baulcombe DC. RNA polymerase IV directs silencing of endogenous DNA. Science. 2005;308:118–20.View ArticlePubMedGoogle Scholar
- Peragine A, Yoshikawa M, Wu G, Albrecht H, Poethig R. SGS3 and SGS2/SDE1/RDR6 are required for juvenile development and the production of trans-acting siRNAs in Arabidopsis. Genes Dev. 2004;18:2368–79.View ArticlePubMedPubMed CentralGoogle Scholar
- Blevins T, Pontes O, Pikaard CS, Meins F. Heterochromatic siRNAs and DDM1 independently silence aberrant 5S rDNA transcripts in Arabidopsis. PLoS One. 2009;4:e5932–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Mette MF, van der Winden J, Matzke MA, Matzke AJ. Production of aberrant promoter transcripts contributes to methylation and silencing of unlinked homologous promoters in trans. EMBO J. 1999;18:241–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Ramírez F, Dündar F, Diehl S, Grüning BA, Manke T. deepTools: a flexible platform for exploring deep-sequencing data. Nucleic Acids Res. 2014;42:W187–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Becker C, Hagmann J, Müller J, Koenig D, Stegle O, Borgwardt K, et al. Spontaneous epigenetic variation in the Arabidopsis thaliana methylome. Nature. 2011;480:245–9.View ArticlePubMedGoogle Scholar
- Xi Y, Li W. BSMAP: whole genome bisulfite sequence MAPping program. BMC Bioinformatics. 2009;10:232.View ArticlePubMedPubMed CentralGoogle Scholar