Small RNA-Mediated De Novo Silencing of Ac/Ds Transposons Is Initiated by Alternative Transposition in Maize.

Although transposable elements (TEs) comprise a major fraction of many higher eukaryotic genomes, most TEs are silenced by host defense mechanisms. The means by which otherwise active TEs are recognized and silenced remains poorly understood. Here we analyzed two independent cases of spontaneous silencing of the active maize Ac/Ds transposon system. This silencing is initiated by alternative transposition, a type of aberrant transposition event that engages the termini of two nearby separate TEs. Alternative transposition during DNA replication can generate Composite Insertions that contain inverted duplications of the transposon sequences. We show that the inverted duplications of two Composite Insertions are transcribed to produce double-stranded RNAs that trigger the production of two distinct classes of small interfering RNAs: a 24-nt class complementary to the TE terminal inverted repeats and noncoding subterminal regions, and a 21- to 22-nt class corresponding to the TE transcribed regions. Plants containing these small interfering RNA-generating Composite Insertions exhibit decreased levels of Ac transcript and heritable repression of Ac/Ds transposition. Further, we demonstrate that Composite Insertions can heritably silence otherwise active elements in trans This study documents the first case of transposon silencing induced by alternative transposition and may represent a general initiating mechanism for silencing of DNA transposons.

Introduction

Transposable elements (TEs) are often silenced by their hosts, but how TEs are initially recognized for silencing remains unclear. Here we describe two independent loci that induce de novo heritable silencing of maize Ac/Ds transposons. Plants containing these loci produce double-stranded RNA and Ac-homologous small interfering RNAs, and exhibit decreased levels of Ac transcript and heritable repression of Ac/Ds transposition. We show that these loci comprise inverted duplications of TE sequences generated by alternative transposition coupled with DNA rereplication. This study documents the first case of transposon silencing induced by alternative transposition and may represent a general initiating mechanism for TE silencing.

Transposable elements (TEs) comprise a large proportion of eukaryotic genomes, including those of important crop plants such as rice (>35%, International Rice Genome Sequencing Project, Matsumoto ), sorghum (62%) (Paterson ), and maize (85%) (Schnable ). Over time, multiple sporadic TE proliferations altered the number and distribution of TE sequences, enhancing genome diversity between species and even among varieties of the same species (Tikhonov ; Sanmiguel and Vitte 2009). These observations indicate that TEs have been, and continue to be, a major natural force driving genome evolution.

Although uncontrolled TEs can be deleterious to the host (Kidwell 1985), the majority of maize transposons are highly methylated (Regulski ; West ) and transcriptionally silenced (Anderson ), which effectively preserves genome integrity. Small RNA-mediated silencing is an efficient means to initiate and maintain repression of both class 1 (RNA) and class 2 (DNA) TEs. For example, the DNA polymerase IV–RNA-dependent RNA polymerase 2 (Pol IV-RDR2) pathway generates 24-nt small interfering RNAs (siRNAs) to mediate RNA-directed DNA methylation (RdDM) at the homologous DNA targets. Because it relies on transcripts from methylated templates, the Pol IV-RDR2 pathway likely serves to reinforce or maintain silencing of previously silenced TEs [reviewed by Law and Jacobsen (2010), Haag and Pikaard (2011), Castel and Martienssen (2013), and Matzke and Mosher (2014)]. Although the mechanisms of maintenance of TE silencing have been extensively studied, the mechanisms responsible for de novo silencing remain relatively obscure. The first evidence for a locus that can silence an active TE came from analysis of Mu killer (Muk), a locus that can heritably silence MuDR transposons in trans in maize. Muk comprises an inverted duplication of the MuDR 5′ terminal inverted repeat (TIR) and a portion of the mudrA gene, which is required for element excision (Lisch 2002). This long inverted repeat is transcribed into double-stranded RNAs (dsRNAs) which initiate siRNA-mediated silencing of intact MuDR elements (Slotkin , 2005; Li ). Recent studies of de novo silencing of active LTR retroelements in Arabidopsis demonstrated the involvement of RdDM (Marí-Ordóñez ; McCue ; Duan ; Panda ). Active TEs are transcribed by Pol II, followed by second-strand synthesis by RDR6. The resulting dsRNA is then processed into 21- to 22-nt siRNAs (Nuthikattu ); these 21- to 22-nt siRNAs can then induce post-transcriptional silencing as well as transcriptional gene silencing, which is associated with DNA methylation (Pontier ; McCue ; Matzke and Mosher 2014). Once RdDM is initiated, silencing can then be maintained via the classic RdDM pathway, which involves Pol IV-mediated transcription from previously methylated sequences (Matzke and Mosher 2014).

Maize Ac/Ds TEs were the first transposons discovered and characterized (McClintock 1948, 1949, 1950, 1951). As members of the class 2 hAT transposon superfamily, Ac/Ds elements are less numerous than class 1 retroelements, which are often highly amplified (Sanmiguel and Vitte 2009). However, Ac/Ds elements can strongly affect gene expression because of their preferential insertions into genes (Vollbrecht ) and the induction of genome rearrangements via alternative transpositions (Zhang and Peterson 1999). Unlike standard transposition reactions, which act on the 5′ and 3′ termini of a single element, alternative transposition acts on the termini of two separate, usually nearby elements (Gray 2000). During alternative transposition, the termini of two nearby, separate TEs can interact with the transposase and insert into linked or unlinked sites, resulting in various chromosomal rearrangements, including duplications, deletions, inversions, and translocations (Zhang and Peterson 2004; Huang and Dooner 2008; Zhang ; Wang ). Moreover, alternative transposition during DNA replication can generate novel structures termed Composite Insertions, which are composed of transposon termini surrounding host genome sequences that were copied from the transposon donor site. Because of their complex and heterogeneous structures, Composite Insertions represent an ongoing source of diverse Ac sequence configurations in the genome (Zhang )

In a previous study of Ac/Ds alternative transposition, Zhang and Peterson (1999) identified an allele termed p1-ww-id1 that induced significant repression of Ac/Ds transposition. Because this stock contains only one copy of Ac, the observed Ac/Ds repression could not be explained by the classic Ac negative dosage effect, in which increased Ac copy number results in reduced frequency and developmental delay of Ac/Ds transposition (McClintock 1949; Brink Nilan 1952). Subsequently, we isolated a second allele, p1-ww-id4, that was independently derived by alternative transposition and that exhibits a similar repression of Ac/Ds transposition.

Here, we show that both of these alleles (1) cause de novo and heritable repression of Ac/Ds transposition and Ac mRNA accumulation; (2) contain Composite Insertions with inverted duplications of Ac sequence; (3) produce Ac-homologous dsRNA transcripts driven by either a flanking host gene promoter or the Ac promoter; and (4) accumulate 21-, 22-, and 24-nt siRNAs corresponding to the region of dsRNAs transcribed from each Composite Insertion. The siRNA profile includes two distinct classes: a 24-nt class corresponding to the TIR and subterminal region of Ac/Ds, and a 21- to 22-nt class homologous to portions of the transcribed region of Ac. These data provide the first evidence for heritable siRNA-mediated silencing of Ac/Ds activity. Our results support and extend a previous model of alternative transposition-induced DNA rereplication to generate Composite Insertions (Zhang ), and also show that the resulting Composite Insertions can trigger heritable silencing of otherwise active elements in trans. This study is the first demonstration of de novo TE silencing induced by alternative transposition, which may represent a general mechanism of self-repression of class 2 transposons.

Materials and Methods

Maize stocks and screen

The maize p1 gene regulates pigmentation in floral organs, and p1 alleles are identified by a two-letter suffix indicating their expression in kernel pericarp (the maternal tissue surrounding the seed) and cob glumes (e.g., “w” for white, “r” for red, and “v” for variegated). The progenitor p1-vv-9D9A allele (Zhang and Peterson 1999) contains an active Ac element and a fractured Ac element (fAc; a terminally deleted Ac element) inserted into p1 intron 2. The p1-vv-9D9A allele was introgressed into maize inbred B73, which has a p1-wr genotype. To screen for variants in Ac activity derived from p1-vv-9D9A, silks of plants of genotype p1-vv-9D9A/p1-wr were crossed by pollen from an Ac tester line of genotype p1-ww; rm-3::Ds. This tester is recessive for a null allele of p1 and is also homozygous for an Ds insertion allele of r1, which exhibits excisions in the aleurone of Ds in the presence of Ac transposase, resulting in colored sectors (Kermicle 1980). The mature ears were screened for multikernel sectors or individual kernels with developmentally delayed (small) purple spots on kernel aleurone due to repressed Ac activity. These candidate kernels were then planted and backcrossed to B73 for more than six generations to generate the material for study. After introgression, the final genotype of plants carrying these candidates is p1-ww-id/p1-wr; r1/r1 in B73 background.

Genomic PCR and DNA gel blot hybridizations

Total DNA was prepared from seedling shoots by using a modified cetyltrimethylammonium bromide extraction protocol (Allen ). HotMaster Taq polymerase (5PRIME) was used in the PCR reactions. PCR reactions were heated at 94° for 2 min; followed by 35 cycles of 94° for 20 sec, 60–68° annealing (depending on primers) for 30 sec, and 65° for 1 min/1 kb expected product length; followed by a final cycle at 65° for 8 min. The sequences of oligonucleotide primers are listed in Supplemental Material, Table S1. For Southern blots, total DNA extracted from seedling shoots was digested with restriction enzymes from Promega and electrophoresed through 0.8% agarose gels. Blotting and hybridization were performed according to standard protocols (Sambrook ); blots were washed in stringent conditions (0.5% SDS, 0.5×SSC at 60°).

Plant growth and tissue collection

Seeds were germinated in SB300 Universal Mix soil mix, and grown in a PGW-40 growth chamber (Percival Scientific, Perry, IA) at 25° for 15 hr in the light and 20° for 9 hr in the dark. Both the above-ground tissues (shoot) and below-ground tissues (root) were harvested from seedlings 14 days after sowing. Nine random plants were pooled for each genotype and tissue type.

RNA isolation and complementary DNA preparation

Total RNA was extracted from seedling shoots and roots by two methods: (1) using the RNeasy Plant mini kit (QIAGEN, Valencia, CA) for small-scale RNA isolation, and (2) using the PureLink Plant RNA Reagent (Life Technologies) for large-scale RNA isolation. The RNA extracted by both methods was then treated with DNase I (New England Biolabs) to remove residual genomic DNA. The RNA extracted by RNeasy Plant mini kit was then converted to complementary DNA (cDNA) by Omniscript RT kit (QIAGEN) primed with oligo-dT and was used as template in the RT-PCR and quantitative RT-PCR experiment. The RNA used in experiments involving small RNA-sequencing and detection of dsRNAs was extracted using the PureLink Plant RNA Reagent.

Quantitative RT-PCR

Quantitative RT-PCR was performed using the Stratagene Mx4000 multiplex quantitative PCR system. Total RNA and cDNA were prepared as described above. PCR was catalyzed by SsoFast EvaGreen Supermixes (Bio-Rad, Hercules, CA) with two technical repeats and three biological repeats. Ac transcript levels measured by primers ac10 and ac11 were normalized by comparison to levels of ubiquitin transcript measured by primers Ubi-f and Ubi-r in the same sample as an internal control. The relative expression level of Ac in each sample was then calculated by comparison to levels of Ac transcript in p1-vv-9D9A, which contains a single active Ac element. Standard deviations were calculated among replications. Student’s t-test was then performed to evaluate the differences among samples at significance level of 0.05.

Detection of dsRNAs

Total RNA, prepared as described above, was treated with RNases A/T1 (Thermo Fisher Scientific) with concentrations of 0, 1.5, and 15 units at 37° for 15 min. The treated RNA was precipitated and reverse transcribed into cDNA by SuperScript III Reverse Transcriptase (Life Technologies) primed with random primers (Thermo Fisher Scientific) as the template for PCR. A seminested PCR was used to amplify regions inside the expected dsRNA (primers ac12 + ac6 followed by ac4 + ac6) to detect dsRNA. Genomic DNA contamination was not detected by using primer ac13 + ac14.

Small RNA high-throughput sequencing and data analysis

Total RNA was extracted as described above and previously (Zuo ); library preparation and sequencing on Illumina platform HiSeq 2000 were performed by Beijing Genomics Institute. Adapter sequences, contamination, and low-quality reads were filtered from raw data. The small RNA sequences were mapped to Ac full-length DNA sequence (4565 bp) by Bowtie (Langmead ). Only perfectly matched short reads were included in the analysis. The mapped reads from each library were normalized to read counts per million reads.

Data availability

The small RNA high-throughput sequencing data analyzed here is deposited in the NCBI under accession number SRP062285. Maize stocks are available upon request. All data necessary for confirming the conclusions of the article were presented within the article, figures, and tables. Supplemental material available at figshare: https://doi.org/10.25386/genetics.12150459.

Results

The model of sister chromatid transposition-induced DNA rereplication

We previously reported that tandem direct duplications and associated Composite Insertions can be produced by reversed end transposition, a type of alternative transposition, followed by DNA rereplication and repair (Zhang and Peterson 2004; Zhang ). The Composite Insertions are bordered by partial or full copies of the Ac transposon and may also include sequences flanking the original Ac donor site. Extending the principle of alternative transposition-induced DNA rereplication, we hypothesized that sister chromatid transposition during DNA replication can produce inverted duplications and Composite Insertions (Figure 1 and Supplemental Material Video 1). In sister chromatid transposition, the Ac transposase acts on a pair of directly oriented Ac termini that are present in the progenitor allele p1-vv-9D9A (Figure 1A). This allele contains a complete Ac element inserted 112 bp from a second element termed fAc, which contains only the 3′ half of Ac. The Ac element is known to preferentially transpose during or shortly after S phase (Greenblatt and Brink 1962; Chen ), possibly because the Ac transposase preferentially interacts with hemimethylated Ac TIRs (Ros and Kunze 2001). Following replication, the Ac 5′ and fAc 3′ termini located on sister chromatids have strand-specific hemimethylation patterns that are competent for transposition. Excision of the 5′ and 3′ termini followed by religation of the host sequences flanking these termini is expected to produce a sister chromatid fusion with a small sequence footprint at the excision site (Weil and Wessler 1993) (Figure 1B). The excised termini may then reinsert at many possible genomic sites; for example, they can reinsert at a proximal site to generate reciprocal duplication/deletion chromatids (Figure 1C) (Zhang and Peterson 1999). Insertion of the Ac termini into the target site commonly produces 8-bp target site duplications flanking the insertion site, a diagnostic feature of Ac transposition (Döring and Starlinger 1984; Peacock ; Pohlman ). When sister chromatid transposition from a replicated donor site inserts into an unreplicated target site, previously replicated sequences are joined to unreplicated DNA. As DNA synthesis continues, replication forks progress into the newly inserted Ac termini, rereplicating the Ac/fAc sequences. Replication may continue into the flanking DNA and extend for 10 kb or more (Zhang ) (Figure 1D). Ultimately, the DNA rereplication forks spontaneously abort, producing two broken ends that fuse together (McClintock 1951). In some cases, fusion of the two broken ends can generate two fAc elements in an inverted orientation, as shown in Figure 1E. In the ensuing mitosis, the two sister chromatids segregate into two daughter cells: one contains a deletion (p1-ww-def) and the other contains an inverted duplication (p1-ww-id) with a Composite Insertion. In the Composite Insertion, the two new fAcs are both derived from Ac 5′ termini and are in an inverted orientation relative to each other. We term the proximal one fAc(P), and the distal one fAc(D), to distinguish them from the original fAc(O) in the progenitor p1-vv-9D9A allele (Figure 1F). Two possible orientations of insertion are predicted by the model: (1) the 5′ terminus of Ac ligates to the proximal side of the insertion target (Figure 1, B–F), and (2) the 3′ terminus of Ac ligates to the proximal side of the insertion target (Figure 1, G–K). These two orientations are found in the p1-ww-id1 and p1-ww-id4 alleles, respectively.

Figure 1

Model of sister chromatid transposition-induced DNA rereplication. (A) Maize chromosome 1 and progenitor allele p1-vv-9D9A with DNA replication bubbles. Solid circle indicates centromere. The p1 gene contains 3 exons (black boxes) with Ac and fAc elements located in intron 2 (gray boxes with solid/open arrowheads indicate the 5′ and 3′ Ac/fAc termini, respectively). The 5′ and 3′ Ac/fAc termini involved in sister chromatid transposition are circled in red and blue, respectively. The unreplicated regions are highlighted in yellow, and the replicated regions to be rereplicated are highlighted in blue. (B and C) Sister chromatid transposition, p1-ww-id1 orientation. (B) Excision of the Ac 5′ and fAc 3′ termini results in excision footprint (marked by X) and fusion of the two sister chromatids. Curved arrows indicate insertion of the excised transposon ends into the target site (short vertical line). (C) The 5′ terminus of Ac is ligated to the centromere-proximal side (marked α) and the 3′ terminus of fAc is ligated to the distal side (marked β). This joins the unreplicated sequences at the insertion site to the previously replicated Ac/fAc sequences. The replication fork containing the sister chromatid fusion is flipped. (D) DNA replication. As DNA replication continues, Ac/fAc sequences are rereplicated (highlighted in green). (E) Completion of DNA replication. Rereplication of Ac/fAc sequences aborts, releasing two broken ends that fuse (dotted line). This produces a Composite Insertion between α and β. During mitosis, sister chromatids separate. One sister chromatid (top) has a deletion, and the other sister chromatid (bottom) carries a corresponding inverted duplication and a Composite Insertion between α and β. (F) Maize chromosome 1 with the p1-ww-id1 allele. The Composite Insertion in the α/β site contains two new fAc elements (solid red arrows) containing Ac 5′ sequences in inverted orientation: fAc(P) (proximal) and fAc(D) (distal). P1-ww-id1 also contains an inverted duplication of the segment from fAc(O) (original) to site β. (G and H) Sister chromatid transposition, p1-ww-id4 orientation. (G) Same as B, except the 3′ terminus of fAc will ligate to the centromere-proximal side (marked γ) and the 5′ terminus of Ac will ligate with the distal side (labeled as δ) of the insertion site. (H–J) Same as C–E. (K) Maize chromosome 1 containing the p1-ww-id4 allele. For animation of alternative transposition mechanism, see Supplementary Materials Video 1.

Alleles p1-ww-id1 and p1-ww-id4 contain structural hallmarks expected from the model of sister chromatid transposition-induced DNA rereplication

Figure 2 summarizes the structural features of the p1-vv-9D9A, p1-ww-id1, and p1-ww-id4 alleles (Figure 2, A–C.) As predicted by the sister chromatid transposition model, both p1-ww-id1 and p1-ww-id4 contain Composite Insertions and inverted duplications as compared to the progenitor p1-vv-9D9A. To further test the model, we examined the sequences at the junctions of the inverted duplications. As previously reported (Zhang and Peterson 1999), p1-ww-id1 contains a typical Ac excision footprint at the junction of the inverted duplication segments. For p1-ww-id4, we PCR-amplified and sequenced the excision site [using primers ac3 + p1(fAc) in Figure 2C; results shown in Figure 2D]. Similar to the footprint of p1-ww-id1 (Zhang and Peterson 1999), p1-ww-id4 shows changes at the first nucleotides flanking the Ac/fAc excision sites, as is typical for sites of Ac excision (File S1).

Figure 2

p1-ww-id1 and p1-ww-id4 contain Composite Insertions and inverted duplications. (A–C) Structure of the progenitor allele p1-vv-9D9A (A), and expected structures of p1-ww-id1 (B) and p1-ww-id4 (C) based on the model of sister chromatid transposition-induced DNA rereplication. PCR primers are labeled as arrows, and sequences homologous to hybridization probe 15 are labeled as hatched boxes. H stands for HindIII restriction site. Other symbols are as in Figure 1. (D) Gel analysis of PCR to amplify the sister chromatid transposition footprints in p1-ww-id1 and p1-ww-id4. Bands were excised from the gel and sequenced (File S1). (E) Gel analysis of PCR to amplify the sequences flanking Composite Insertions. Bands were excised from the gel and sequenced to identify the target site duplications flanking each Composite Insertion (File S2). (F) Genomic Southern blot produced by digestion with HindIII and hybridization with probe 15. (G) Gel analysis of PCRs to determine the orientation of sister chromatid transposition (see text for details).

We then isolated the sequences flanking the Composite Insertions in each allele by inverse PCR, and confirmed the insertion junctions by direct PCR using primer pairs f1(1) + ac1 for p1-ww-id1, and f4(2) + ac1 for p1-ww-id4 (Figure 2E). The sequences flanking the Composite Insertions contain perfect 8-bp target site duplications (“GCCTCGCT” in p1-ww-id1 and “GCCCGGAT” in p1-ww-id4; File S2) characteristic of Ac transposition.

The sister chromatid transposition model predicts that the p1-ww-id1 and p1-ww-id4 alleles contain inverted duplications extending from the p1 locus at 48.6 Mb of maize chromosome 1 to the Composite Insertion sites. The sequences flanking the Composite Insertions in p1-ww-id1 and p1-ww-id4 are located at positions 51.8 and 48.9 Mb, respectively, on chromosome 1, maize B73 RefGen_v4 reference genome. These results indicate that the p1-ww-id1 and p1-ww-id4 contain inverted duplications of 3.2 and 0.3 Mb, respectively.

The presence of duplicated segments in p1-ww-id1 and p1-ww-id4 is tested by Southern blot using HindIII digests and probe fragment 15 (Figure 2F). The fragments of 8 and 6 kb originate from the proximal and distal side of p1, respectively. In the progenitor allele p1-vv-9D9A (labeled as 9d), the intensities of the 8- and 6-kb bands are approximately equal, consistent with their single-copy status in p1-vv-9D9A. In the p1-ww-id1 and p1-ww-id4 alleles, the 8-kb band intensity is approximately twice that of the 6-kb band, consistent with a duplication of the 8-kb segment proximal to p1 in those alleles. Hybridization signals in lane 9d are overall weaker due to the fact that there is less DNA loaded on the gel. This experiment confirms the presence of a duplication extending to the proximal side of p1 as predicted by the sister chromatid transposition model.

Finally, we tested the orientations of the Composite Insertions in the p1-ww-id1 and p1-ww-id4 alleles. These orientations depend on whether the excised 5′ Ac and 3′ fAc(O) are joined to the proximal or distal sides of the Composite Insertion site (Figure 1, B and G), and are correlated with the presence or absence of fAc(O) adjacent to Ac in the original p1 donor site (Figure 1, F and K). Therefore, we used PCR primers p1(9d) + pp1’ to amplify from fAc(O) across Ac and into the flanking p1 gene sequence, producing a 6.6-kb product; this 6.6-kb product was observed in p1-ww-id1 and the control p1-vv-9D9A, but not p1-ww-id4 (Figure 2G). This result indicates that in p1-ww-id1, the 3′ terminus of fAc(O) is ligated to the centromere-distal side of the target site (β in Figure 1B). Conversely, primers f4(2) and pp1’ produce a 4.5-kb band from p1-ww-id4, but not p1-ww-id1 and p1-vv-9D9A (Figure 2G). This 4.5-kb band reflects the absence of fAc(O) adjacent to Ac in the p1-ww-id4 allele, indicating that in p1-ww-id4, the 5′ terminus of Ac is ligated to the centromere-distal side of target site (δ in Figure 1G). In summary, these results show that the p1-ww-id1 and p1-ww-id4 alleles contain transposon donor excision footprints, Composite Insertions in two possible orientations flanked by target site duplications, and proximal duplications, all as predicted by sister chromatid transposition with insertion into unreplicated target sites.

p1-ww-id1 and p1-ww-id4 alleles induce de novo and heritable silencing of Ac

The maize p1 gene is required for kernel pericarp (seed coat) pigmentation (Zhang and Peterson 1999). The p1-vv-9D9A allele contains an Ac element in p1 intron 2 (Figure 2A), which blocks p1 expression and results in colorless pericarp. Somatic excision of Ac can restore p1 function and thus produce red clonal sectors on the kernel pericarp (a variegated pericarp phenotype). The p1-vv-9D9A allele produces frequent kernel pericarp sectors, indicating that the Ac in p1-vv-9D9A is fully active. This is confirmed by crosses of p1-vv-9D9A to an Ac tester line of genotype r1-m3::Ds (Figure 3A). The Ac tester line contains a nonautonomous Ds element inserted in r1 (red1), a gene required for anthocyanin biosynthesis in kernel aleurone. Without an active Ac, the Ds insertion blocks r1 function, resulting in colorless aleurone; however, with an active Ac, Ds can be transposed from r1, producing purple aleurone sectors (Kermicle 1980; Lechelt ). The activity of Ac can thus be assessed by the size and frequency of purple spots on the kernel: large spots from early Ds excision events indicate a fully active Ac, while small spots from delayed Ds excisions indicate repressed Ac activity. As shown in Figure 3A, crosses of p1-vv-9D9A to the Ac tester line produce kernels with a coarsely spotted pattern typical of active Ac. In contrast, crosses of p1-ww-id1 with Ac testers produce kernels with fine spots (Figure 3B, top kernel), while p1-ww-id4 gives no spots in test crosses (Figure 3C, top kernel). These results suggest that the Ac elements in p1-ww-id1 and p1-ww-id4 are repressed. Moreover, p1-ww-id1 and p1-ww-id4 can induce trans-dominant repression of a fresh and active Ac. This can be seen in the kernels produced by crossing p1-ww-id alleles to the p1-vv-9D9A allele that contains a single active Ac. Kernels heterozygous for p1-vv-9D9A/p1-ww-id1 show very few and fine aleurone sectors (Figure 3B, second kernel), while kernels produced by crossing p1-vv-9D9A with p1-ww-id4 had no visible aleurone sectors (Figure 3C, second kernel). Note that these kernels exhibit a normal variegated pericarp phenotype because the kernel pericarp is a maternal tissue and thus is not affected by a repressive p1-ww-id allele introduced through the pollen in this cross.

Figure 3

Genetic crosses indicate the in-trans and heritable repression of Ac/Ds initiated by p1-ww-id1 and p1-ww-id4. (A) Progenitor allele p1-vv-9D9A contains a single active Ac element, shown by the heavily spotted kernel aleurone from the cross of p1-vv-9D9A and Ac tester. Cross: p1-vv-9D9A x r1-m3::Ds. (B) Ac activity is repressed in p1-ww-id1, shown by the fine spots in kernel aleurone in the first cross of p1-ww-id1 and Ac tester (p1-vv-id1 x r1-m3::Ds). The second cross (p1-vv-9D9A x p1-ww-id1) introduces a fresh active Ac in p1-vv-9D9A, which is repressed in trans by p1-ww-id1 (indicated by small aleurone spots). The third cross (r1-m3::Ds x p1-vv-9D9A/p1-ww-id1) shows that silencing of Ac in p1-vv-9D9A is heritable, as small spots are maintained even after p1-vv-9D9A is segregated from p1-ww-id1. (C) Parallel crosses of p1-ww-id4 as in B. Ac activity is repressed in the line of p1-ww-id4, shown as absence of kernel aleurone spots in the first cross between p1-ww-id4 and Ac tester. The repression is trans-dominant, shown by the absence of spots in kernels produced by cross of active Ac from p1-vv-9D9A by p1-ww-id4. Ac repression is heritable, shown by the fine spotting in kernels from the third cross in which p1-vv-9D9A is segregated from p1-ww-id4.

To test whether the repression of Ac induced by p1-ww-id1 and p1-ww-id4 reflects heritable silencing, we crossed p1-vv-9D9A/p1-ww-id plants to the Ac tester line; this cross separates the Ac element in p1-vv-9D9A from each p1-ww-id allele. If the repression of Ac in p1-vv9D9A is relieved following segregation from p1-ww-id, then we would expect to see ∼50% coarsely spotted kernels (containing p1-vv-9D9A) and ∼50% weakly or nonspotted kernels (containing p1-ww-id). However, if Ac is heritably silenced, then most or all progeny kernels would again show few, small Ds excision sectors. The results show that all the kernels produced by crossing p1-vv-9D9A/p1-ww-id1 with the Ac tester show few or zero purple sectors (0% heavily spotted seeds; Figure 3B, third kernel). The corresponding ear for this cross is shown in Figure S1, and kernel count data are provided in Table S2, Cross 6. Parallel genetic tests show that Ac repression induced by p1-ww-id4 is also heritable, but is less stable than that induced by p1-ww-id1. Among 160 progeny kernels from the cross of Ac tester by p1-vv-9D9A/p1-ww-id4, 126 are weakly or nonspotted (Figure 3C, third kernel), and 34 kernels are heavily spotted. Assuming the heavily spotted kernels are derived from a reactivated Ac in the p1-vv-9D9A allele, this indicates that Ac repression is lifted in 42.5% (34/80) of kernels following segregation from p1-ww-id4 (Fisher exact test P < 0.00001, indicating a significant difference from the p1-vv-9D9A ear without the initial silencing; ear for this cross is shown in Figure S2, and kernel count data are provided in Table S2, cross 8). In other words, the silencing of Ac induced by p1-ww-id4 is heritable in only ∼57.5% of the progeny kernels that carry p1-vv-9D9A.

We also assessed the maintenance of Ac silencing in succeeding generations following segregation from the p1-ww-id1 and p1-ww-id4 alleles. Plants of genotype p1-vv-9D9A*/p1-ww (weakly spotted kernels, where * indicates Ac silenced by exposure to p1-ww-id1 or p1-ww-id4 in the prior generation) were crossed again to the Ac tester stock r1-m3::Ds. In these crosses, 50% of the kernels receive p1-ww (no Ac) and are nonspotted, while the remaining 50% of kernels receive p1-vv-9D9A* and their spotting pattern can reflect the stability of Ac silencing. Interestingly, p1-ww-id1 and p1-ww-id4 induce different levels of silencing stability. For p1-vv-9D9A* initiated by p1-ww-id1 and separated from it for two generations, all the kernels are nonspotted or weakly spotted, indicating near complete maintenance of Ac repression (Figure S1; Table S2, cross 7). In contrast, the silencing of Ac initiated by p1-ww-id4 and segregated from it for two generations is not efficiently maintained. Among 144 progeny seeds from the cross of Ac tester by p1-vv-9D9A/p1-ww, 72 seeds are expected to inherit the p1-vv-9D9A allele and thus contain Ac. In the progeny kernels we observed 52 heavily spotted seeds, meaning that Ac is reactivated in 72.2% (52/72) of the kernels (Fisher exact test P-value is 0.0116, indicating a significant difference from the p1-vv-9D9A ear without the initial silencing; Figure S2; Table S2, cross 9). Thus, the repression of Ac initiated by p1-ww-id4 is maintained in only ∼27.8% of progeny kernels after two generations.

To determine the mechanism of Ac repression, we performed quantitative RT-PCR experiment using Ac-specific primers. Compared to the progenitor allele p1-vv-9D9A, the p1-ww-id1 and p1-ww-id4 alleles express significantly decreased levels of Ac transcript (Student’s t-test P = 0.0015 and 0.0004, respectively) (Figure 4), suggesting that repression occurred via transcriptional and/or post-transcriptional gene silencing.

Figure 4

The Ac transcript level is decreased in p1-ww-id alleles. (A) Schematic structure of the full-length Ac transposable element. Solid and open triangles indicate Ac 5′ and 3′ termini, respectively. Black boxes with numbers are Ac exons 1–5; boxes with subtending lines indicate introns. Primers used in the quantitative RT PCR (ac10 and ac11) are labeled as arrows; note that primer ac11 spans Ac intron 4. (B) Quantitative RT-PCR measurement of Ac transcript levels. Bars labeled by * indicate transcript levels that are significantly different from p1-vv-9D9A by Student’s unpaired t-test with P < 0.05.

Composite Insertions of alleles p1-ww-id1 and p1-ww-id4 contain inverted repeats of fAc fragments

We performed genomic Southern blots to determine the structures of the Composite Insertions in p1-ww-id1 and p1-ww-id4. The results are summarized in Figure 5, and the raw data are shown in Figure S3 and Figure S4. First, the sizes of Composite Insertions are estimated by using restriction enzymes that do not cut within the Ac element. For p1-ww-id1, digestion with SalI and hybridization with probe 1D gives a band of 10–12 kb, indicating a Composite Insertions of ∼6.2–8.2 kb. Consistent with this, enzyme PstI and probe 1P give an estimated size of ∼6.7 kb. We then determined the internal Composite Insertion structure by digesting with a series of endonucleases with known sites in Ac. We observed the expected bands produced by digestion at the SspI and ScaI sites in fAc(P), and the expected bands from digestion at the SspI of fAc(O) and BglI sites of fAc(D). These results indicate that the Composite Insertion of p1-ww-id1 includes the first 1842 bp sequences from the Ac 5′ terminus of fAc(P) and the first 729 bp sequences from fAc(D). If either fAc(P) or fAc(D) contained the EcoRI site located near the center of Ac, then digestion with EcoRI and hybridization with probe 1P would yield a band of size 7.6–9.6 kb. No such band is observed, and instead we detect a larger band of ∼14 kb. This indicates that both fAc(P) and fAc(D) were truncated prior to reaching the EcoRI site at position 2486 bp from the Ac 5′ end. We conclude that the total size of Composite Insertion in p1-ww-id1 is <7.1 kb, consistent with the ∼6.7 kb size estimated by PstI digestion.

Figure 5

Summary of Southern blot results that elucidate the internal structures of Composite Insertions. Schematic structures of p1-ww-id1 and p1-ww-id4 are shown in A and B, respectively, with restriction sites (vertical lines; E: EcoRI; Sc: ScaI; P: PstI; Sp: SspI; Sl: SalI; B: BglI; N: NcoI) and probes (hatched boxes) shown. Bands observed from Southern blot (Figures S3 and S4) are indicated by braces with actual sizes labeled by each restriction enzyme.

Similar assays were performed on p1-ww-id4 (Figure 5B and Figure S4). PstI and NcoI digestions both show the total size of the Composite Insertion in p1-ww-id4 as ∼6.4–8.4 kb. The Composite Insertion contains the ScaI site in both fAc(P) and fAc(D), indicating that both fAc fragments include at least 1842 bp of Ac 5′ terminal sequences. Additionally, fAc(D) contains the EcoRI site but fAc(P) does not, indicating that fAc(P) is truncated prior to the EcoRI site at position 2486 bp from Ac 5′ terminus. Combining the results from ScaI and EcoRI digests, the size of fAc(P) could be from 1842 to 2486 bp. By subtracting the size of fAc(P) (1.8–2.5 kb) from the total size of Composite Insertion of p1-ww-id4 (6.4–8.4 kb), we infer the size of fAc(D) is 2.5–3.9 kb. Extensive attempts were made to PCR amplify and sequence the exact fusion point connecting fAc(D) and fAc(P) in both Composite Insertions; unfortunately, these attempts failed, most likely due to the long perfectly inverted repeats present in both Composite Insertions, which would give rise to self-annealed hairpins with extensive secondary structures in the template DNA.

The inverted repeats of fAc fragments are transcribed to produce dsRNAs

The inverted fAc(P) and fAc(D) fragments in both p1-ww-id1 and p1-ww-id4 retain the Ac promoter. Moreover, consistent with the characteristic preferential insertions of Ac/Ds into linked genes (Vollbrecht ), both Composite Insertions are located in the middle of annotated maize genes: the Composite Insertion in p1-ww-id1 is located at intron 1 of Zm00001d028930, while the Composite Insertion in p1-ww-id4 is located at intron 5 of Zm00001d028863 (Figure 6, A and B). The locations of Composite Insertions within genes could allow their expression by read-through transcription from the flanking genes’ promoters.

Figure 6

Composite Insertions in p1-ww-id1 and p1-ww-id4 produce dsRNA transcripts. (A) Structure of the Composite Insertion in p1-ww-id1, located in intron 1 of gene Zm00001d028930. Bent arrow indicates transcription from the host gene promoter to generate read-through transcripts of the Composite Insertion in root tissues, but not shoot. Labeled arrows indicate primers used in RT-PCR. (B) Structure of the Composite Insertion in p1-ww-id4, located in intron 5 of Zm00001d028863. The promoter of Zm00001d028863 is not active in either root or shoot tissues. Bent arrows indicate transcription from the native Ac promoters. (C) Test of genomic DNA contamination in RNA samples. Primers ac13 + ac14 are located in fAc(O) outside the dsRNA region, in Ac exons 3 and 5 (Figure 6, A and B). Bands of 744 bp and 272 bp are expected from genomic DNA and spliced cDNA, respectively. Control lanes “1” and “2” are from templates of genomic DNA and root cDNA, respectively, of p1-ww-id1; lane “w” is from a water template as a control for PCR contamination. (D) Test for chimeric transcripts by RT-PCR using primers 1tgf + ac3 in p1-ww-id1 (lanes 1–3), and 4tgr + ac3 in p1-ww-id4 (lanes 4–6). Chimeric transcript initiating from flanking gene promoter was detected only in p1-ww-id1 root. No chimeric transcript from flanking gene promoter was detected in p1-ww-id4. (E) Detection of dsRNA by RNase protection assay. Total RNA from the indicated tissues and alleles were treated with DNase1 and three quantities of RNase A/T1 (0, 1.5, and 15 units). Treated RNA samples were then analyzed by a seminested PCR (primers ac12 + ac6; followed by ac4 + ac6). Control lanes “1” and “2” are from templates of genomic DNA and root cDNA, respectively, of p1-ww-id1; lane “w” is from a water template as a control for PCR contamination.

To test for read-through transcription, we performed RT-PCR experiments using total RNA extracted from root and shoot tissues of plants containing p1-ww-id1 and p1-ww-id4. No genomic DNA contamination is detected as illustrated in Figure 6C. For p1-ww-id1, chimeric transcripts are detected in root, but not shoot (Figure 6D); this is the same pattern of expression as the flanking host gene Zm00001d028930 (Wang ) (Figure S5). The chimeric transcript includes sequences from exon 1 and intron 1 of the host gene Zm00001d028930, and sequences from the 5′ terminus of Ac (File S3), consistent with the structure expected of read-through transcripts. In contrast, gene Zm00001d028863, the host gene of the Composite Insertion of p1-ww-id4, is not expressed in either root or shoot, and we did not detect chimeric p1-ww-id4 transcripts in those tissues (Figure 6D).

Prompted by the inverted fAc repeat structures in both p1-ww-id1 and p1-ww-id4 Composite Insertions, we searched for Ac-homologous dsRNAs using RNA protection assays. In this experiment, total RNA samples were treated with two concentrations of RNase A/T; linear RNAs can be digested by the RNase, while dsRNAs are resistant. Following RNase treatment, RNAs were reverse transcribed into cDNA, which was used as template for PCR. A seminested PCR was performed to detect sequences protected by dsRNA (primers ac12 + ac6 followed by ac4 + ac6; Figure 6, A and B). We observed bands of the expected size in p1-ww-id1 root, but not shoot. This pattern matches the expression pattern of the host gene Zm00001d028930 and of the presence of chimeric transcripts. We also detected Ac-homologous dsRNAs in both root and shoot of p1-ww-id4 (Figure 6E), suggesting that p1-ww-id4 may be transcribed from its own Ac promoter. Interestingly, both the chimeric transcript and the dsRNAs seem to retain introns (Figure 6, D and E and File S3), suggesting that the hairpin nature of these transcripts may interfere with normal splicing.

We concluded that the inverted repeats of fAc elements in the Composite Insertions of p1-ww-id1 and p1-ww-id4 are transcribed into dsRNAs, either by the flanking gene promoter as in p1-ww-id1, or by the retained Ac promoter as suggested for p1-ww-id4. The accumulation of dsRNAs can thus provide a potential template for production of siRNAs that may directly induce Ac silencing.

Small RNAs derived from Composite Insertions are detected in both p1-ww-id1 and p1-ww-id4 alleles

The presence of dsRNAs in p1-ww-id1 and p1-ww-id4 suggests a possible mechanism of Ac silencing by siRNAs. Therefore we performed high-throughput sequencing of small RNAs from root and shoot tissues of p1-ww-id1 and p1-ww-id4 plants. Progenitor allele p1-vv-9D9A and standard maize inbred B73 were used as controls. As shown in Figure 7, siRNAs of 21, 22, and 24 nt that map to Ac are dramatically enriched in both shoot and root from both p1-ww-id1 and p1-ww-id4 plants compared to the same tissues in p1-vv-9D9A and B73. In p1-ww-id1 root and shoot, the enriched small RNAs map to a 1- to 2486-nt region from the Ac 5′ terminus (Figure 8), consistent with the extent of dsRNA predicted by the Composite Insertion inverted repeat structure. In p1-ww-id4, the enriched small RNAs map to the region 1–3400 nt from the Ac 5′ terminus, again matching the expected stem-loop structure with the size range of 2486–3900 nt predicted by the Composite Insertion structure. This finding suggests that small RNAs are synthesized from both the double-strand stem and single-strand loop regions of p1-ww-id4.

Figure 7

Small RNAs mapped to Ac. Compared to inbred B73 and progenitor allele p1-vv-9D9A, small RNA levels are increased significantly in both p1-ww-id1 and p1-ww-id4. Colored bars indicate the abundance of 21-, 22-, and 24-nt small RNA classes.

Figure 8

Detection and mapping of small RNAs in p1-ww-id1 and p1-ww-id4. (A) The expected stem-loop structures in p1-ww-id1 and p1-ww-id4. The schematic structures include a ruler with corresponding nucleotide position in Ac. The loop of p1-ww-id4 is open for better presentation of the size. The transcription of Composite Insertion of p1-ww-id1 is driven by the host promoter, and the two inverted fAc fragments in p1-ww-id1 are relatively symmetrical (2100–2486 nt of each indicated by Southern blot), predicting a hairpin structure in the dsRNA contains first ∼ (2100th–2486th) nucleotides of Ac. The transcription of Composite Insertion of p1-ww-id4 is driven by the Ac promoter, and thus the transcript starts from 335th nucleotide of Ac. The inverted fAc fragments in p1-ww-id4 differ in size, with one fAc fragment of 1842–2486 nt and the other fAc of 3400–3900 nt indicated by Southern blot. This asymmetrical structure should give rise to a transcript with a stem containing 335th ∼ (1842nd–2486th) nucleotides of Ac and a loop containing (1843rd–2487th) ∼ (3400th–3900th) nucleotides of Ac in p1-ww-id4. (B) Small RNAs mapped to Ac are visualized by the Integrative Genomics Viewer (Robinson ). The y axis of the coverage track was standardized to the same scale for each sample. Colors indicate size of siRNAs: black, 21 nt; red, 22 nt; blue, 24 nt.

Interestingly, the different size-classes of siRNAs map to distinctly different regions of Ac. The 24-nt siRNAs are enriched in the 5′ TIR/sub-TIR regions of Ac (Figure S6, B–D); this region includes the native Ac promoter and the presence of these 24-nt small RNAs is consistent with a Pol IV/RDR2-based maintenance of heritable silencing (Law and Jacobsen 2010; Pikaard 2013; Matzke and Mosher 2014; Li ). In contrast, the 21- to 22-nt siRNAs predominate in the transcribed regions of Ac (Figure S6E). This size of siRNAs have been shown to be associated with the Pol II-RDR6 pathway (Pontier ; Marí-Ordóñez ; Nuthikattu ; Matzke and Mosher 2014; McCue ; West ; Duan ; Panda ) and thus may induce de novo silencing of active Ac elements. The 21- to 22-nt siRNAs are also mapped to the single-strand loop region in p1-ww-id4 (Figure 8; the expanded loop region is shown in Figure S7). This is consistent with a transitive process, by which RDR6 converts the single-stranded RNAs in the loop into dsRNAs, followed by the processing of dsRNAs into small RNAs. The ratio of 21- to 22-nt siRNA is increased in the loop region in shoot compared to root tissues (Table S3; chi-squared test gives P < 0.01 in shoot and P = 0.24 in root), suggesting tissue-specific variation in production or stability of siRNA species. Finally, we noted that 21-nt siRNAs in p1-ww-id1 are enriched in Ac intron sequences, consistent with the observation that both the chimeric Ac transcript and dsRNAs retain introns in p1-ww-id1 (Figure S8).

We observed similar profiles of siRNA in both root and shoot tissues of p1-ww-id1, even though no Composite Insertion transcripts or dsRNAs were detected in p1-ww-id1 shoots. Possibly, this might result from long-range transportation of siRNAs from root (where dsRNAs are produced) to shoot (where dsRNAs are absent), as described previously in Arabidopsis (Chitwood and Timmermans 2010; Molnar ). In any event, the fact that the siRNAs in p1-ww-id1 and p1-ww-id4 correspond well to the structure of the two Composite Insertions in each allele strongly suggests that these Composite Insertions are the source of the siRNAs.

Discussion

In this study, we characterized two naturally occurring maize alleles (p1-ww-id1 and p1-ww-id4) derived by alternative transposition of the transposable element Ac. Both alleles elicit trans-dominant and heritable repression of Ac-induced transposition, reduce accumulation of Ac mRNA, and contain inverted duplications and novel Composite Insertions derived from sister chromatid transposition-induced DNA rereplication. The inverted duplications are of 3.2 Mb in p1-ww-id1 and 0.3 Mb in p1-ww-id4; both Composite Insertions contain inverted repeats of the Ac 5′ terminal sequences, of 2.1–2.5 kb in p1-ww-id1 and 1.8–3.9 kb in p1-ww-id4. These Ac inverted repeats are transcribed either from a flanking gene promoter in the case of p1-ww-id1, or possibly from the Ac promoter itself in the case of p1-ww-id4. dsRNA transcripts are detected from both alleles, presumably due to fold back pairing of the inverted repeat transcripts.

The alternative transposition-initiated de novo and heritable silencing in the p1-ww-id alleles is consistent with mechanisms of transcriptional gene silencing and post-transcriptional gene silencing reported for class 1 elements in Arabidopsis (de novo silencing: Marí-Ordóñez ; Nuthikattu ; McCue ; Duan ; Panda ; heritable silencing: Law and Jacobsen 2010; Haag and Pikaard 2011; Castel and Martienssen 2013; Matzke and Mosher 2014). In the p1-ww-id alleles described here, the 21- to 22-nt siRNAs enriched in the transcribed region of the Composite Insertions, and particularly in the loop region of p1-ww-id4, suggests the involvement of RDR6 in the de novo silencing. The reduced levels of Ac transcript can be formally explained by mRNA degradation via post-transcriptional gene silencing, and/or transcriptional gene silencing of the Ac promoter triggered by the Pol II-RDR6 pathway. Transcriptional silencing of Ac seems most likely, as silencing is maintained even after the id1 and id4 alleles are removed by meiotic segregation (Figure 3, B and C). The production of 24-nt siRNAs corresponding to the TIR and subterminal regions of Ac signals the possible involvement of the Pol IV-RDR2 RdDM pathway in the maintenance of transcriptional gene silencing by cytosine methylation in all sequence contexts (Law and Jacobsen 2010; Zemach ; Stroud ; Li ). The siRNA-independent pathway may also participate in the maintenance of symmetrical methylation (Stroud ; Matzke and Mosher 2014).

The structure and effects of the p1-ww-id1 and p1-ww-id4 Composite Insertions are remarkably similar to the Muk allele that represses Mutator transposons in maize (Slotkin , 2005; Li ; Burgess ), and hence the p1-ww-id1 and p1-ww-id4 alleles can be considered as examples of “killers” targeting a different DNA TE superfamily. Muk and the “Ac killers” share similar features: (1) they are both initiated from naturally occurring inverted duplications of partial DNA transposon sequences; (2) the transcription of inverted duplications in Muk and Ac killer in p1-ww-id1 are both driven by nearby promoters, while p1-ww-id4 is possibly transcribed by the Ac promoter, producing dsRNAs as the precursor of 21-, 22-, and 24-nt siRNAs; and (3) both Muk and Ac killer can trigger in trans and heritable silencing of an active element. These similarities suggest a general mechanism for the heritable silencing of active DNA transposons [see Burgess ]. Active TEs are prone to rearrangements, some of which will be competent to silence otherwise active elements in trans. Thus, it is likely that if an element is active for sufficient time, it will eventually produce transposon killers. To the extent that selection favors a loss of TE activity, these killers would tend to spread within a population, eventually resulting in the heritable epigenetic silencing of homologous elements in that population. In addition to silencing resident active TEs, these killers may also serve as a reservoir of “antigens” that immunize genomes from subsequent invasion by active elements. Further, if expressed in specific cells of the male and female gametophyte, TE killers may serve as sources of small RNAs that reinforce TE silencing in germinal lineages (Martínez and Slotkin 2012).

Our discovery of siRNA-mediated de novo and heritable silencing of Ac/Ds adds another layer of complexity to the regulation of Ac/Ds activity. Previous studies of the regulation of Ac/Ds transposons reveal a complex relationship among Ac dosage, transcript abundance, transposase level, DNA methylation, and transposition frequency. Ac has a GC-rich subterminal region, and hypermethylation of this region is associated with transcriptional silencing (Kunze ; Brutnell and Dellaporta 1994; Conrad and Brutnell 2005). It is this region of the p1-ww-id alleles that we see being targeted by 24-nt small RNAs. Given that these small RNAs are also observed in B73, they likely arise from inactive Ac elements present in that genetic background. Additionally, splicing of Ac mRNA is inefficient and inaccurate, producing a large number of aberrant Ac transcripts in Arabidopsis, representing post-transcriptional regulation (Jarvis ). The aggregation of Ac transposase protein in tobacco (Kunze ), petunia, and maize (Heinlein ), may represent a type of post-translational regulation responsible for the Ac negative dosage effect first observed by McClintock (1948, 1949, 1950). Molecular analysis showed that increased dosage of Ac results in increased Ac transposase mRNA and protein levels, at least in the maize wx-m7 allele that was tested (Fußwinkel ). In this study, we find that the p1-ww-id1 and p1-ww-id4 alleles derived from alternative transposition-induced DNA rereplication initiate transgenerational silencing of active elements via trans-acting siRNAs.

These results also expand our understanding of the effect of alternative transposition on the structure of the maize genome. Previous studies have described the ability of alternative transposition to induce structural rearrangements such as deletions, inversions, translocations, duplications (Zhang and Peterson 2004; Huang and Dooner 2008; Zhang , 2013), and exon shuffling (Zhang ; Wang ). Here we show that p1-ww-id alleles derived from sister chromatid transposition contain novel inverted-repeat structures (Composite Insertions) generated by DNA rereplication. The fact that similar structures and mechanisms effect silencing in members of two unrelated transposon superfamilies (hAT and Mu) suggests that inverted repeat-induced siRNA may represent a general mechanism of spontaneous TE silencing.