It is well documented that transposable elements (TEs) can regulate the expression of neighbouring genes. However, their ability to act in trans and influence ectopic loci has been reported rarely. We searched in rice transcriptomes for tissue-specific expression of TEs and found them to be regulated developmentally. They often shared sequence homology with co-expressed genes and contained potential microRNA-binding sites, which suggested possible contributions to gene regulation. In fact, we have identified a retrotransposon that is highly transcribed in roots and whose spliced transcript constitutes a target mimic for miR171. miR171 destabilizes mRNAs encoding the root-specific family of SCARECROW-Like transcription factors. We demonstrate that retrotransposon-derived transcripts act as decoys for miR171, triggering its degradation and thus results in the root-specific accumulation of SCARECROW-Like mRNAs. Such transposon-mediated post-transcriptional control of miR171 levels is conserved in diverse rice species.https://doi.org/10.7554/eLife.30038.001
An organism’s genome contains all of the DNA the individual needs to survive and grow. Transposons are pieces of DNA that can move around the genome. They make up almost half of human DNA and over 85% of the DNA of major crop plants like maize, barley and wheat.
When transposons move they can cause harmful changes in the regions where they insert into the DNA and so cells have mechanisms in place to tightly control the activities of the transposons. However, some transposons cause changes to DNA that are beneficial to the organism. Thus, the relationship between transposons and their host organisms is an example of a delicate but mostly peaceful coexistence. Although the cellular mechanisms controlling transposons are quite well known, the extent to which the transposons affect the ability of organisms to survive, develop and reproduce is poorly understood.
A family of proteins known as the SCARECROW-like transcription factors are important for the roots of plants to develop properly. In other organs such as the leaves or flowers these proteins can cause developmental defects, so the plants carefully control where the proteins are made. Thus, plants normally produce a molecule called miR171 in leaves and flowers, but not in roots, that inhibits protein production by binding to and destabilising the RNA molecules that act as templates to make these proteins.
Cho and Paszkowski have now identified a transposon that produces an RNA molecule with similarities to the RNA templates used to make the SCARECROW-like transcription factors. The experiments show that this transposon RNA is found in very high amounts in roots and mimics the transcription factor RNA so well that miR171 binds to it. This inactivates miR171 in roots to allow the SCARECROW-like transcription factors to be produced.
These findings reveal a new mechanism by which transposons may regulate how plants develop and provide possible new approaches for boosting the growth of rice and other crop plants. Similar regulatory interactions between transposons and their host DNA may also be present in animals and other organisms.https://doi.org/10.7554/eLife.30038.002
Transposable elements (TEs) constitute a large fraction of eukaryotic genomes. Given their mutagenic potential and largely unknown functions, they were often considered as genomic parasites that are silenced by host epigenetic mechanisms (Fultz et al., 2015; Girard and Hannon, 2008). However, there is increasing evidence that TEs contribute to various chromosomal functions, to the evolution of genomes by increasing genetic variation, and to the direct regulation of genes (Lisch, 2013). Several studies have revealed that TEs in plants endow genes with both coding and regulatory sequences (Lisch, 2013). For example, the Arabidopsis transcription factors FHY3 and FAR1, involved in light signalling, are derived from the transposase of the Mutator-like DNA transposon (Hudson et al., 2003). The domestication of hAT and Mutator-like transposases contributed to the evolution of the DAYSLEEPER and MUSTANG gene families, respectively. DAYSLEEPER was shown to play a critical role in plant development (Bundock and Hooykaas, 2005; Cowan et al., 2005; Knip et al., 2013; Knip et al., 2012). More recently, a protein derived from the transposase of the Pif/Harbinger transposon family was shown to be an inhibitor of POLYCOMB REPRESSIVE COMPLEX 2 (Liang et al., 2015).
TEs residing outside protein-coding regions of genes can influence their expression by interfering with promoters, providing enhancers, or altering RNA processing and/or epigenetic regulation. For example, TEs residing in introns or UTRs may alter the availability of splicing sites and/or splicing efficiencies. They can also shift polyadenylation signals or supply binding sites for miRNA and RNA-binding proteins (Feschotte, 2008).
In contrast to the numerous examples of local influence on gene regulation in cis, examples of TEs mediating the regulation of distant genes are rare. For example, the Arabidopsis ddm1 mutant, which is impaired in epigenetic suppression of transposon-derived transcription, accumulates 21-nt small RNAs derived from Athila retrotransposons. These small RNAs impair the levels and the translation of the mRNA of OLIGOURIDYLATE BINDING PROTEIN1 (UBP1) activated by abiotic stress (McCue et al., 2012). Interestingly, although more than 20 genes in Arabidopsis have putative binding sites for transposon-derived small RNAs that would allow regulation analogous to that of UBP1, such a network of interactions has not been documented so far (McCue et al., 2013).
TE transcripts often contain features resembling micro RNA (miRNA) genes (Li et al., 2011) or sequences that are potential targets of miRNAs (Creasey et al., 2014). miRNAs are a class of small non-coding RNAs that, directed by their sequences, selectively repress gene expression by translation inhibition or cleavage of mRNA (Rogers and Chen, 2013). Interestingly, miRNAs can interact both with their target mRNAs and also with other RNAs containing similar binding sites. Such ‘rival’ RNAs, which were seen as competing endogenous RNAs (ceRNAs) (Kartha and Subramanian, 2014; Salmena et al., 2011; Tay et al., 2014), can be derived from pseudogenes or long non-coding RNAs and also from protein-coding mRNAs. They may also appear in the form of circular RNAs. Although studies of their transgenic overexpression support activity as miRNA sponges, a possible biological role has not been demonstrated so far by their loss of function (Thomson and Dinger, 2016).
Although, transposon-derived transcripts were not considered previously to be a source of ceRNAs, we decided to search rice transcriptomes for indications of this activity. Approximately 35% of the rice genome consists of TEs, which is significantly higher than in Arabidopsis (14%) (International Rice Genome Sequencing Project, 2005; Arabidopsis Genome Initiative, 2000). More important and similar to maize (Erhard et al., 2009; Hollick et al., 2005; Parkinson et al., 2007), rice mutants impaired in epigenetic silencing of transposons, such as DICER-LIKE 4 (DCL4) or RNA-DEPENDENT RNA POLYMERASE 6 (RDR6), show severe developmental abnormalities (Liu et al., 2007; Song et al., 2012), whilst the corresponding mutants of Arabidopsis have no apparent morphological aberrations (Xie et al., 2005). This dissimilarity suggests that restrained TE-derived transcription is important for rice development (Liu et al., 2007; Song et al., 2012; Wei et al., 2014).
Here we have specifically investigated TE-derived transcripts as potential regulators of rice development. We found that numerous TEs display patterns of transcriptional activity that are associated with particular plant tissues. Remarkably, a significant proportion of TE-derived transcripts correlate with the mRNA levels of genes transcribed in the same tissues and the two classes of transcripts often share patches of homology. Notably, the sequences of these patches appear to be significantly enriched for miRNA-binding sites. Therefore, we investigated whether some of the transposon-derived transcripts act as ceRNAs. Experiments to test this hypothesis led to the identification of a novel domesticated retrotransposon that is highly expressed in rice roots and that acts as a ceRNA post-transcriptionally controlling the level of miR171. This particular ceRNA is also a target mimic of miR171, which potentially enhances its sponging activity towards miR171. Tissue-specific adjustment of miR171 levels is essential to the proper development of roots and this appears to be regulated by a retrotransposon-derived ceRNA. We demonstrated that mutations in its miR171-binding site result in an abnormal root system.
To examine tissue-specific abundance of TE-derived transcripts in rice, we accessed publicly available RNA sequencing (RNA-seq) datasets for various tissues of rice (Figure 1A). We considered only the datasets of Japonica rice, cultivar Nipponbare and applied the same data-processing pipeline to raw sequencing results generated in different laboratories (details in the Materials and methods). This way we achieved consistent results and samples representing particular tissues were clustered together (Figure 1A). Such combined dataset yielded 2961 transcribed TEs (filtered for maximal RPKM (Reads Per Kilobase per Million reads)>1). Remarkably, the TEs were transcribed in most rice tissues and their transcriptomes exhibit clear tissue specificity (Figure 1A and Figure 1—figure supplement 1A). The rice expression patterns differ from those of Arabidopsis, where TEs are activated in a non-selective way and only in seed endosperm and the vegetative cells of pollen grains (Figure 1—figure supplement 1B) (Slotkin et al., 2009). Thus, in rice, the two-dimensional correlation matrix of TE transcriptomes showed distinct TE groups reflecting their RNA abundance in various tissues and at different developmental stages (Figure 1—figure supplement 1A). In contrast, Arabidopsis TEs exhibit more uniform expression patterns (Figure 1—figure supplement 1B).
We detected TE-derived transcripts in rice tissues that do not contribute to the germ line (e.g. endosperm, leaves and roots). These transposon activities, even when resulting in insertions, would not be transmitted to the next generation. In the case of such apparently unproductive TE activity, reactivated TEs may possibly be regulatory or, as in Arabidopsis, may be the RNA substrates of small RNAs involved in TEs silencing (Creasey et al., 2014), or may simply reflect insignificant transcriptional noise.
For regulatory activity that influences gene expression, we assumed that TE transcripts would share homologies with the transcripts of co-expressed genes. Indeed, multiple alignments revealed homology patches in approximately 64% of co-transcribed TEs, while silent TEs matched only by 41% (Figure 1—figure supplement 1C). Moreover, sequence comparison of expressed TEs and corresponding genes showed a strong bias towards sense direction, by which the aligned regions of TEs were much longer than for TEs matching gene transcript in the antisense direction (Figure 1—figure supplement 1D, left panel). Intriguingly, orientation bias was observed when the alignments were sought against processed mRNAs of genes but not when introns were included (Figure 1—figure supplement 1D).
Next, we sought sequence features differentially enriched in the co-expressed TEs and found enrichment for miRNA-binding sites (Figure 1B). This raised the possibility that some TE transcripts interfere with miRNA-mediated gene regulation, possibly by competing for miRNA binding. As two different RNAs interacting with the same miRNA would need to be co-expressed in the same tissue, we compared tissue-specific expression of TEs and matching genes. We found 763 and 400 TE-gene pairs in sense and antisense orientation, respectively, including 282 sense and 111 antisense pairs with correlation coefficient above 0.5 (Figure 1C). Such correlated expression patterns between mRNAs of genes and TE-derived transcripts in sense orientation was most evident in roots (Figure 1—figure supplement 1E). Collectively, the results of our examination of tissue-specific transcriptomes are consistent with the hypothesis that TEs regulate gene expression by miRNA sequestration.
To test this hypothesis, we rigorously re-analysed 61 root-specific rice transcriptomes and selected a particular TE, which we named MIKKI (‘decoy’ in Korean), for further investigation (Figure 2A).
First, we validated RNA-seq results of root-specific transcription of MIKKI by RT-qPCR (Figure 2B). To distinguish the spliced transcript from precursor mRNA (pre-mRNA), we designed primers across exon junction or within the intron, respectively (Figure 2B, left and right panel). The RT-qPCR results confirmed that the mature MIKKI transcript is highly abundant in roots, present at low levels in leaves, and almost absent in panicles (Figure 2B, left panel). A similar expression pattern was observed for unspliced RNA but at much lower levels (Figure 2B, right panel). These results are consistent with tissue-specific regulation of MIKKI at the transcriptional level.
MIKKI is a TE-derived locus which includes Osr29 Long Terminal Repeat (LTR) retrotransposon. Based on sequence divergence between the two LTRs, an Osr29 element transposed about 3.7 million years ago (mya, Figure 2C and Figure 2—figure supplement 1A). We also found sequences of three further retrotransposons, BAJIE, Osr30 and Osr34, inserted subsequently into Osr29 (Figure 2C). Advanced degeneracy prevented estimate of the insertion times of Osr30- and Osr34-related sequences; however, the generation time of the solo LTR derived from the BAJIE family was estimated to be approximately 1.2 mya (Figure 2C and Figure 2—figure supplement 1B). The MIKKI gene product was predicted to encode just a partial reverse-transcriptase (RTase) protein and no other protein domains were found (Figure 2—figure supplement 1C). Given that several amino acid residues essential for catalytic activity of RTase are mutated in MIKKI’s RTase (Figure 2—figure supplement 1D), it seems unlikely that MIKKI’s RTase domain would be active. Thus, we concluded that MIKKI is not expected to have regulatory role at the protein level. Most important, the mature transcript of such a rearranged Osr29 (MIKKI) was found to contain an imperfect binding site for miR171, generated by a splicing event joining BAJIE solo LTR sequences to specific sequences of Osr29 (Figure 2C,E and Figure 2—figure supplement 1E). miR171 is one of the miRNAs conserved across the plant kingdom and previous studies revealed that Arabidopsis miR171 (ath-miR171) is abundant in flowers but sparse in roots (Figure 2—figure supplement 2C–E; Llave et al., 2002). Rice miR171 (osa-miR171) displays a similar expression pattern (Figure 2D, left panel and Figure 2—figure supplement 2B). Thus, miR171 levels seem to be similar in particular tissues of these two distant species, highest in reproductive organs and lowest in roots.
It is well documented that ath-miR171 targets mRNAs encoding SCARECROW-Like (SCL) transcription factors for cleavage and, thus, SCL transcript levels display patterns opposite to miR171 (Llave et al., 2002). The same SCL transcript distribution was observed in rice (Figure 2D, right panel), implying the regulation of SCL transcript stability also by osa-miR171. Moreover, the sequence identity of rice and Arabidopsis SCL mRNAs across the miR171-binding region is also consistent with the evolutionary conservation of miR171-mediated cleavage of SCL transcripts (Figure 2—figure supplement 2F). Indeed, analyses of the RNA degradome in rice panicles (Wu et al., 2009) revealed specific cleavage of OsSCL21 mRNAs at the osa-miR171 binding region (Figure 2E, left panel). We also examined whether the MIKKI transcript is also targeted by osa-miR171 but found no signals indicative of site-directed cleavage of MIKKI transcripts at the putative miR171-binding site (Figure 2E, right panel). Importantly, there are two mismatches in the miR171-binding region of MIKKI at positions 11th and 14th. Conservation of nucleotides at these sites is known to be essential for target RNA cleavage (Jeong et al., 2013; Liu et al., 2014; Llave et al., 2002). It is, therefore, possible that the mismatches around the cleavage sites in MIKKI transcripts attenuate the cleavage activity of osa-miR171. Altogether, these data are consistent with the possibility that MIKKI is a target mimic of osa-miR171 in rice roots.
To examine the target mimicry of the MIKKI transcript towards miR171, we overexpressed MIKKI in both rice and Arabidopsis (Figure 3—figure supplement 1A and B, top panel), and applied RT-qPCR analyses. miR171 levels were downregulated in independent transgenic lines generated from both plant species (Figure 3A, top panel and Figure 3—figure supplement 1B, middle panel) and the transcript levels of target genes were markedly upregulated (Figure 3A, bottom panel and Figure 3—figure supplement 1B, bottom panel). Previous studies revealed abnormalities in floral organs of Arabidopsis plants in which ath-miR171 levels were decreased by overexpression of artificial target mimics (Ivashuta et al., 2011; Todesco et al., 2010). Consistent with these observations, plants ectopically overproducing MIKKI transcripts also displayed severe defects in reproductive organs and low fertility (Figure 3B,C and Figure 3—figure supplement 1C).
To address directly the developmental role of the MIKKI retrotransposon, we generated the MIKKI mutants mikki-1 and mikki-2 using CRISPR-Cas9 (Miao et al., 2013). To ensure the targeting specificity of guide RNAs (gRNA), we designed them to target the unique junction region between Osr29 and BAJIE (Figure 3—figure supplement 2A). Transgenic plants were examined by sequencing for mutations in this region and two independent alleles were found (Figure 3—figure supplement 2A and B). The mikki-1 allele had a 2 bp deletion at the splice donor site that resulted in retention of the intron. Intron retention disrupted the miR171-binding site and generated multiple premature stop codons (Figure 3—figure supplement 2A). This transcript is likely recognized by a nonsense-mediated mRNA decay pathway and rapidly turned over (Shoemaker and Green, 2012). Indeed, the RT-qPCR analyses revealed thousand-fold reduction of MIKKI transcripts in the mikki-1 mutant (Figure 3—figure supplement 2C). The mikki-2 allele has an 8 bp deletion in the region containing the miR171-binding site (Figure 3—figure supplement 2A). This deletion did not alter RNA levels but was predicted to lose target recognition by osa-miR171 (Figure 3—figure supplement 2A and C).
Next, we performed RT-qPCR on the wildtype (wt) and the mutants. The levels of osa-miR171 were high in both mikki-1 and mikki-2 (Figure 3D, top panel and Figure 3—figure supplement 2E). This correlated with a decrease in RNA levels of OsSCL21 targeted by osa-miR171 (Figure 3D, bottom panel). In Arabidopsis, mutation of SCLs leads to defects in root development (Wang et al., 2010). From a Korean rice seed bank we obtained two independent mutant alleles of OsSCL21 that showed the highest transcript levels among OsSCLs targeted by miR171 (Figure 3—figure supplement 3). Similar to Arabidopsis, the roots of both osscl21 mutants were shorter than wt (Figure 3—figure supplement 2D). Subsequently, we examined the development of mikki-1 and mikki-2 roots. Root lengths were affected in both mutants, resembling mutants in the OsSCL21 gene (Figure 3E and Figure 3—figure supplement 2D). Histological analyses were also performed to observe the cellular consequences of MIKKI mutation. Both mutants showed reduced cell elongation above meristematic region, while the cell widths were similar to wt (Figure 3F and G). These data are consistent with the hypothesis that MIKKI negatively regulates osa-miR171 levels in rice roots, acting through a ceRNA containing target mimic site for osa-miR171.
Next, we asked whether post-transcriptional regulation by a ceRNA with target mimicry is the major regulatory mechanism governing tissue-specific levels of osa-miR171. For this, we determined the levels of the primary transcript of osa-miR171 (pri-osa-miR171) in MIKKI overexpression and mutant plants (Figure 4A,B and Figure 4—figure supplement 1A). The abundance of pri-osa-miR171 was similar in different rice tissues and was not affected by the alteration of MIKKI transcript levels or mutation, implying that mature osa-miR171 is regulated post-transcriptionally by the activity of MIKKI.
In contrast to rice, the tissue-specific distribution of primary transcripts of miR171 in Arabidopsis was the same as the mature miRNA, which is consistent with transcriptional regulation and thus an entirely different regulatory mechanism (Figure 4C).
MIKKI is present and has a conserved structure in AA-genome Oryza species (Figure 4—figure supplement 2A and B), suggesting strong selective advantage of this particular transposon. MIKKI is present in the genomes of Oryza sativa ssp. indica, O. rufipogon, O. nivara, O. barthii, and O. glaberrima (Figure 4—figure supplement 2A). Furthermore, insertion of the BAJIE-derived solo LTR and the resulting intron with a miR171-binding site at the splice junction are perfectly conserved (Figure 4—figure supplement 2B), implying that the formation of a splicing-dependent miR171 binding site retained in these related species. We examined MIKKI splicing in five of these species using available RNA-seq data (Zhai et al., 2013; Zhang et al., 2016; Zhang et al., 2014) and detected identical splicing patterns of the critical intron 4 (Figure 4—figure supplement 3A). Moreover, we found that the MIKKI homolog of Indica rice displays a developmental expression pattern similar to Japonica rice (Figure 4—figure supplement 3B).
We also examined tissue-specific levels of primary transcripts and mature miR171 in monocotyledonous Brachypodium (Figure 4—figure supplement 1B and C). As in rice, the primary transcripts of miR171 were high in all tissues examined, suggesting analogous post-transcriptional control of miR171 levels. However, so far we have not identified an MIKKI-related element in the genome of Brachypodium.
Since the transcription of TEs is usually controlled by epigenetic mechanisms, we examined DNA methylation and selected histone modifications associated with MIKKI in different rice tissues (Figure 5). In roots, where MIKKI is actively transcribed, its upstream region was enriched in lysine 4 tri-methylation of histone H3 (H3K4me3) and depleted of suppressive H3K9me2 (Figure 5B). DNA methylation levels were also lower than in panicles (Figure 5C). Analysis on the public RNA-seq data generated from rice OsDCL3a RNAi knock-down lines also showed derepression of MIKKI (Figure 5—figure supplement 1). These epigenetic signatures were well correlated with tissue-specific transcription of MIKKI.
In summary, we propose a model by which MIKKI influences rice root development via the regulation of osa-miR171 levels by tissue-specific expression of a ceRNA encoding target mimicry of miR171 (Figure 6).
In Arabidopsis, transposable elements are mostly silenced by epigenetic mechanisms preventing their transcription during development of the sporophyte. In contrast, in plant species such as maize or rice, transposon-derived transcripts are detected during specific developmental transitions or in various organs (Li et al., 2010; Tamaki et al., 2015). Since most of the tissues examined do not contribute to the formation of gametophytes and thus to transgenerational inheritance, the benefits of transposon transcription remain unclear. The prevalent view is that their transcripts are a source of mobile small RNAs that, if transported into germline progenitor cells, would contribute to silencing of transposons there and thus prevent their transgenerational accumulation (Calarco et al., 2012; Creasey et al., 2014; Slotkin et al., 2009). An alternative explanation for the developmental regulation of transposon-derived transcription, however, is that their transcripts have particular functions in a specific tissue or organ. To examine this latter possibility, we systematically analysed tissue-specific transcriptomes of rice transposons by re-analysing available raw RNA sequencing data. These analyses uncovered a surprisingly high fraction of TE-derived transcripts in specific tissues or developmental stages of rice plants. We also observed that transposon-derived transcripts are enriched in miRNA binding sites. It has been proposed that Arabidopsis miRNAs trigger generation of transposon-derived small RNAs that contribute subsequently to transposon silencing. These small RNAs may later spread to other cell types (Creasey et al., 2014; Martínez et al., 2016). This scenario is consistent with the transposon defence hypothesis described above and similar mechanisms could certainly operate in rice.
Alternatively, transposon-derived RNAs containing binding sites for miRNAs could also act as ceRNAs and there are experimental examples in Arabidopsis supporting this possibility. The first and the physiologically important example of an Arabidopsis ceRNA was non-coding RNA INDUCED BY PHOSPHATE STARVATION 1 (IPS1, Franco-Zorrilla et al., 2007). The IPS1 locus is transcriptionally activated upon phosphate starvation and encodes RNA that binds to miR399. miR399 also targets the mRNA PHO2 gene that encodes an E2 ubiquitin-conjugating-like enzyme affecting the phosphate content of shoots. Importantly, the miR399-binding region in IPS1 RNA has a 3-nt bulge in the cleavage site and, similar to MIKKI, it is resistant to cleavage. This unproductive binding of miR399 to IPS1 RNA triggers miRNA degradation (Yan et al., 2012), thus reducing its level. The mechanism of such a miRNA decoy was termed ‘target mimicry’. It has been suggested that the small RNA-specific nucleases in Arabidopsis reduce miRNA levels when target mimics are overexpressed (Ramachandran and Chen, 2008; Yan et al., 2012). However, the mechanism of miRNA degradation has not been fully elucidated and it is not known how miRNA is recognized for degradation when associated with RNA encoding target mimics but not when associated with mRNAs encoding its bona fide targets.
Subsequently, considerable efforts have been made to identify miRNA target mimics in genomes and transcriptomes of plants (Fan et al., 2015; Meng et al., 2012) and mammals (Clark et al., 2014; Helwak and Tollervey, 2014; Imig et al., 2015). In addition, further examples of their biological activities have received experimental support (Franco-Zorrilla et al., 2007; Wang et al., 2015); H.-J. Wu et al., 2013). Notably, experimental support for the activity of ceRNAs, including those from plants containing ‘target mimicry’ sites, is based on transgenic overexpression of micro RNAs or transcripts containing their target mimics. Unfortunately, these assays significantly alter the natural stoichiometry of such regulatory systems. As a consequence, the role of ceRNAs in nature has been queried, given that their relatively low abundance may be insufficient to significantly alter the levels of very dynamically regulated miRNAs (Thomson and Dinger, 2016). In addition, mathematical modelling of miRNA target competition supports the notion that target mimics and miRNAs must be at particular levels for maximal effect of target mimicry (Bosia et al., 2013; Figliuzzi et al., 2013; Yip et al., 2014; Yuan et al., 2015). Since MIKKI transcripts in rice roots are in ample excess over OsSCL mRNAs (Figure 2B and D), the proportions of the components of the MIKKI-miR171-OsSCLs module appear to be naturally sufficient for highly effective regulatory activity. Moreover, we have directly examined the biological role of MIKKI mRNA as a ceRNA containing target mimics of osa-miR171 by using site-directed mutation of its miRNA binding site or mutations in the splicing site that result in its depletion. These experiments, circumventing artificial changes in the stoichiometry of the interacting components, revealed a regulatory function of MIKKI transcripts in the proper development of rice roots.
The tissue-specific posttranscriptional control of miR171 in rice contrasts with restriction of miR171 availability in Arabidopsis roots, which seems to be determined by transcription of miR171. Activation of miR171 is thought to be by AtSCL proteins binding directly to promoters (Xue et al., 2014). Given that the promoter sequences of miRNA-encoding genes are generally less conserved than miRNA-coding regions (Zhu et al., 2015), corresponding genes in different plant species may be the subject of different transcriptional regulation. The levels of miR171 in different plant tissues are decisive for plant development and posttranscriptional control, implemented by a domesticated retrotransposon, reinforces differential organ distribution of mature miR171. This mechanism seems to be highly conserved among distantly related rice species, suggesting a strong selective advantage. The fact that osa-miR171 and MIKKI transcript levels were not affected by OsSCL21 mutation (Figure 3—figure supplement 3D and E) further supports the hypothesis of independent and unique post-transcriptional controlling mechanism emerged. Interestingly, miR171 in the model grass Brachypodium also seems to be regulated posttranscriptionally but a potential ceRNA has not yet been identified.
It has been stated frequently that transposition bursts make a large contribution to host plant genome structure and function (Lisch, 2013). However, examples of transposon-mediated control of plant development through the regulation of distantly located genes have not been reported. Although TEs are an ample source of miRNAs and can potentially act as potent target mimic, they have not been examined for such function. Our study provides the first example of the TE-derived target mimic and thus a novel mechanism for TEs acting as trans-acting regulators of genes. However, for an effective target mimic several conditions should be met, including high transcript levels, good binding affinity and advantageous stoichiometry to miRNA target genes, therefore it is currently still difficult to predict how general this regulatory mechanism will be. We have discovered many more putative TE target mimics in the rice genome (listed in Supplementary file 1). However, genetic interference with these hypothetical regulatory loops is difficult due to multicopy components and thus genetic redundancy, as is the case for most transposons. Potentially, the increasingly efficient site-specific alteration of genomes by CRISPR-Cas9 or a population genetics approach using natural variation may improve functional accessibility to the transposon-derived fractions of plant genomes and reveal the extent of their regulatory input.
Husks of rice seeds were removed and the seeds were surface sterilized in 20% bleach and germinated in ½-strength Murashige and Skoog media. The 2-week-old seedlings were transferred to soil and grown to maturity in a greenhouse. Root and leaf samples were harvested from 2-week-old seedlings and panicles collected immediately after heading. The rice strains used in this study were Oryza sativa ssp. japonica cv. Nipponbare, O. sativa ssp. japonica cv. Hwayoung and O. sativa ssp. indica cv. IR64. The mutant lines of osscl21 were identified from a T-DNA tagged population established at Kyunghee University, Korea and genotyped for the selection of homozygotes. (http://cbi.khu.ac.kr/RISD_DB.html/).
For rice transcriptome analysis, raw FASTQ files of the following RNA sequencing datasets were downloaded: GSE16631, DRA000385, SRP008821, DRA002310, SRP028376 and GSE50778. The adapter-trimmed clean reads were mapped to the reference genome of MSU7 using TOPHAT 2.0, with most of the options set to the default but with some optimization (-g 300). Cufflinks was used to call the RPKM. All the downstream analyses and plotting, for example, heat-maps of transcript levels and correlation matrix, box/violin plots, and statistical analysis were performed in R studio. The transcriptome data of different Oryza species were obtained from PRJNA264484, PRJNA264480, PRJNA264485, SRP070627, GSE41797 and analysed as described above. Arabidopsis transcriptome data were obtained from PRJNA314076 and analysed as for the rice transcriptome using the TAIR10 reference genome.
For the degradome and small RNA analysis, raw FASTQ files from GSE18251 and GSE16350 were downloaded and mapped uniquely to the MSU7 reference genome using BOWTIE2. The resulting BAM file was visualized by IGV. Arabidopsis small RNA-seq data were obtained from GSE28591 and analysed as for the degradome. HTSeq was used to calculate the read counts for each miRNA.
Total RNA was extracted using the RNeasy Plant mini kit (Qiagen, Hilden, Germany) following the manufacturer’s recommendations. Reverse transcription reactions were performed using the VILO RT kit (Invitrogen, California, USA) with a random hexamer for priming. Real-time quantitative PCR was carried out using the Roche Light-Cycler (Roche, Basel, Switzerland) in a volume of 10 µl and analysed by the ΔΔCt method. All data in this study are the average of three biological replicates performed in technical triplicate ± standard deviation and normalized against eEF1α. An Ncode miRNA first-strand cDNA synthesis kit (Invitrogen, California, USA) was used for miRNA quantification; normalization was against miR166, which is expressed constitutively in rice (Figure 2—figure supplement 2B). Sequences of primers used are listed in Supplementary file 2.
A total of 15 µg of RNA was separated on 15% urea-TBE gels (Thermo Fisher Scientific, Massachusetts, USA), transferred to Hybond N + nylon membranes (Amersham Biosciences, Buckinghamshire, UK) and fixed chemically using EDC (Sigma-Aldrich, Missouri, USA). The membranes were prehybridized for 1 hr and hybridized for at least 16 hr in DIG Easy Hyb buffer (Roche, Basel, Switzerland) at 37°C. Membranes were washed twice with 2X SSC (saline sodium citrate), 0.1% SDS. Immunological detection of DIG-labeled probe was performed using DIG wash and block buffer set (Roche, Basel, Switzerland) and DIG luminescent detection kit (Roche, Basel, Switzerland). Luminescent signal was detected with Amersham Imager 600 (Amersham Biosciences, Buckinghamshire, UK).
Roots from 4-day-old rice plants were fixed in FAA (formaldehyde, acetic acid, ethanol) solution overnight in cold room, wax embedded using Leica ASP300 tissue processor (Leica Biosystems, Wetzlar, Germany) and sectioned by 4 µm using microtome (Leica Biosystems, Wetzlar, Germany). After dewaxing and ethanol washing, samples were stained with Calcofluor White (Sigma-Aldrich, Missouri, USA) to visualize cell wall. Images were taken with a Zeiss LSM 700 confocal microscope (Leica Biosystems, Wetzlar, Germany).
Leaf and root samples were collected from Oryza sativa ssp. japonica cv. Nipponbare plants grown for 2 weeks under short-day conditions (10 hr light/14 hr dark). Panicles were harvested immediately after heading from plants grown in the greenhouse. Samples were crosslinked with 1% formaldehyde, flash-frozen, and ground in liquid nitrogen. Chromatin was fragmented by sonication and immunoprecipitated using the following antibodies: H3K4me3 (ab8580; abcam, Cambridge, UK) and H3K9me2 (ab1220; abcam, Cambridge, UK). The immunoprecipitated DNA was quantified by qPCR and normalized against levels of input and the reference genes indicated. All the results are presented as means ± standard deviation (sd) of three biological replicates performed in technical triplicate. Sequences of primers used are listed in Supplementary file 2.
Genomic DNA was extracted from rice tissues using the DNeasy plant mini kit (Qiagen, Hilden, Germany). The Epitect bisulfite kit (Qiagen, Hilden, Germany) was used for bisulfite conversion of unmethylated cytosines. The primer design and data analysis used kismet (Gruntman et al., 2008). At least 15 clones from each sample were analysed. Sequences of primers used are listed in Supplementary file 2.
A cDNA fragment of MIKKI from the start to the stop codon was amplified, cloned into the pUN1901 and pGPTVII binary vectors, and transformed into rice and Arabidopsis, respectively (Walter et al., 2004; Wang et al., 2004). For rice transformation, embryo-derived 2-week-old calli were immersed in agrobacterium-containing media. Transgenic rice plants were obtained after antibiotic selection and differentiation of plantlets. The detailed procedure was as described previously (Nishimura et al., 2006). Arabidopsis transformation was by the floral dip method as described previously (Clough and Bent, 1998).
The oligonucleotide of the designed guide RNA was inserted into the pOs-sgRNA entry vector and shuttled to the pH-Ubi-cas9-7 destination vector by the LR recombination reaction. The resulting binary vector was transformed into rice as described above (Miao et al., 2013). To detect mutation by CRISPR-Cas9, the region containing the targeted region from genomic DNA was amplified, cloned into the pGEM T-easy vector (Promega, Wisconsin, USA) and sequenced. Selected mutant lines were cultured to the next generation to segregate away the T-DNA and individuals homozygous for mutant allele were selected. Sequences of primers used are listed in Supplementary file 2.
Genome sequences of the selected Oryza species were obtained from Ensembl Plants (http://plants.ensembl.org/). Local BLAST analysis was performed manually using the MIKKI genomic sequence, followed by multiple sequence alignment in ClustalW2 and visualization by FigTree v.1.4.2 and boxshade v.3.21.
LTR retrotransposon age was estimated as described previously (Ma and Bennetzen, 2004). Briefly, for Osr29, the divergence was calculated from sequence degeneracy of two LTRs. Age of insertion was computed using the equation: T = D/2 t, where T is the time since insertion, D is the divergence and t is the substitution rate of 1.3 × 10−8 per site per year as proposed previously (Ma and Bennetzen, 2004). For BAJIE solo LTR, the sequence was compared to the consensus of BAJIE LTR sequences. We assumed that the consensus BAJIE LTR sequence represents the youngest copy and used the equation of T = D/t.
MIKKI, LOC_Os06g02304; OsSCL8, LOC_Os02g44360; OsSCL21, LOC_Os04g46860; OsSCL29, LOC_Os06g01620.
MUSTANG is a novel family of domesticated transposase genes found in diverse angiospermsMolecular Biology and Evolution 22:2084–2089.https://doi.org/10.1093/molbev/msi202
Transposable elements and the evolution of regulatory networksNature Reviews Genetics 9:397–405.https://doi.org/10.1038/nrg2337
Target mimicry provides a new mechanism for regulation of microRNA activityNature Genetics 39:1033–1037.https://doi.org/10.1038/ng2079
Silencing of active transposable elements in plantsCurrent Opinion in Plant Biology 27:67–76.https://doi.org/10.1016/j.pbi.2015.05.027
Conserved themes in small-RNA-mediated transposon controlTrends in Cell Biology 18:136–148.https://doi.org/10.1016/j.tcb.2008.01.004
miR-CLIP capture of a miRNA targetome uncovers a lincRNA H19-miR-106a interactionNature Chemical Biology 11:107–114.https://doi.org/10.1038/nchembio.1713
How important are transposons for plant evolution?Nature Reviews Genetics 14:49–61.https://doi.org/10.1038/nrg3374
Oryza sativa dicer-like4 reveals a key role for small interfering RNA silencing in plant developmentThe Plant Cell Online 19:2705–2718.https://doi.org/10.1105/tpc.107.052209
Targeted mutagenesis in rice using CRISPR-Cas systemCell Research 23:1233–1236.https://doi.org/10.1038/cr.2013.123
A protocol for agrobacterium-mediated transformation in riceNature Protocols 1:2796–2802.https://doi.org/10.1038/nprot.2006.469
Twine: display and analysis of cis-regulatory modulesBioinformatics 29:1690–1692.https://doi.org/10.1093/bioinformatics/btt264
Biogenesis, turnover, and mode of action of plant microRNAsThe Plant Cell 25:2383–2399.https://doi.org/10.1105/tpc.113.113159
Endogenous microRNA sponges: evidence and controversyNature Reviews Genetics 17:272–283.https://doi.org/10.1038/nrg.2016.20
A practical vector for efficient knockdown of gene expression in rice (Oryza sativa L.)Plant Molecular Biology Reporter 22:409–417.https://doi.org/10.1007/BF02772683
Widespread long noncoding RNAs as endogenous target mimics for microRNAs in plantsPlant Physiology 161:1875–1884.https://doi.org/10.1104/pp.113.215962
Christian S. HardtkeReviewing Editor; University of Lausanne, Switzerland
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Regulation of rice root development by a retrotransposon acting as a microRNA sponge" for consideration by eLife. Your article has been reviewed by two peer reviewers, one of whom, Christian Hardtke, acted as the Reviewing and Senior Editor. The reviewers have opted to remain anonymous.
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
As you can see from the detailed reviews below, the reviewers were quite positive about your paper, but would like you to address a number of issues regarding data presentation and documentation before publication. Please see the comments below for details.
Cho and Paszkowski present evidence for microRNA sponge action of a particular transposon in rice. This is a timely subject and reveals another layer of uncommon gene regulation, and is thus interesting for a wide audience. Apparently the manuscript has gone through in depth peer review previously and been transferred to this journal, so the present version is already revised and a response to previous reviews was provided.
Overall, I found the manuscript to be of high quality, and the interpretation of the data credible. Rice is a more difficult system to work with than, say, Arabidopsis, and although I would have liked to see some additional experiments (for example, higher resolution expression analyses), I acknowledge that in the rice system, it is not always easily feasible. For instance, the genetics is in part limit because no complementation lines are provided, but this is offset by the analysis of two independent alleles in each case. A few suggestions for further improvement of the manuscript:
– Statistical tests are not always clear. Which has been used, and why? Please amend figure legends and Materials and methods. In some cases, no test is provided. For example, Figure 1—figure supplement 1C: is the difference between expressed and silenced TEs in homology patches really significant? A simple F-test should give the answer. Also, Figure 2B, etc. In some figures, error bars are missing? Or not visible because too small (e.g. Figure 4B)?
– How was the Arabidopsis expression correlation matrix produced? From published data, following the same procedure as for rice? Please clarify.
– Figure 3A: miR171 was not at all detectable in MIKKI O.E. lines? Or value too small? Maybe show with log scale in cases of strong expression differences.
– Regarding miRNA expression: maybe it would be generally more correct to talk about miRNA abundance in the different backgrounds, rather than expression, which most would associate with gene expression in this case (which, as authors show, is not affected in the different backgrounds).
– The whole aspect of differential miR171 regulation in rice versus Arabidopsis or Brachypodium sort of distracts without further investigation. SCL control of miR171 in Arabidopsis has been shown and can remain, but the Brachypodium data are preliminary and not really helpful, could be removed.
– The most important part of the paper is the demonstration that miR171 abundance is increased in the mikki mutant alleles. Could the authors provide small RNA blots to demonstrate that in mikki roots miR171 is specifically up-regulated as compared to a control miRNA? qRT-PCR speaks for it, but given this is a key result, a technically independent verification would be desirable.
This manuscript describes the role of a domesticated rice transposable element in regulating rice development through modulation of miRNA171. The findings are interesting and the data are convincing. As the manuscript already underwent rigorous peer review, I have relatively little to add:
1) A more detailed explanation of the plots in Figure 1—figure supplement 1A-B would be very helpful. It's not clear what is being correlated, and thus the distinction between TE expression in rice and Arabidopsis does not appear well-supported.
2) Please explain how TE age was estimated in the Materials and methods section (in addition to citing the reference).
3) I believe the authors are not using the term "positive selection" as intended. Positive selection means increased frequency of a new, advantageous mutation in a population. The authors appear to argue that mutations that alter MIKKI structure are disfavored – this is negative selection.https://doi.org/10.7554/eLife.30038.053
- Jungnam Cho
- Jerzy Paszkowski
- Jerzy Paszkowski
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
We thank Drs. Zhengming Wang and Weibing Yang for technical support and advice on small RNA blots and histological analyses. This work was supported by European Research Council (EVOBREED) ; Gatsby Fellowship [AT3273/GLE].
- Christian S. Hardtke, Reviewing Editor, University of Lausanne, Switzerland
© 2017, Cho et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.