- Open Access
Phased secondary small interfering RNAs in Panaxnotoginseng
© The Author(s) 2018
Published: 19 January 2018
Recent results demonstrated that either non-coding or coding genes generate phased secondary small interfering RNAs (phasiRNAs) guided by specific miRNAs. Till now, there is no studies for phasiRNAs in Panax notoginseng (Burk.) F.H. Chen (P. notoginseng), an important traditional Chinese herbal medicinal plant species.
Here we performed a genome-wide discovery of phasiRNAs and its host PHAS loci in P. notoginseng by analyzing small RNA sequencing profiles. Degradome sequencing profile was used to identify the trigger miRNAs of these phasiRNAs and potential targets of phasiRNAs. We also used RLM 5’-RACE to validate some of the identified phasiRNA targets.
After analyzing 24 small RNA sequencing profiles of P. notoginseng, 204 and 90 PHAS loci that encoded 21 and 24 nucleotide (nt) phasiRNAs, respectively, were identified. Furthermore, we found that phasiRNAs produced from some pentatricopeptide repeat-contain (PPR) genes target another layer of PPR genes as validated by both the degradome sequencing profile and RLM 5’-RACE analysis. We also found that miR171 with 21 nt triggers the generations of 21 nt phasiRNAs from its conserved targets.
We validated that some phasiRNAs generated from PPRs and TASL genes are functional by targeting other PPRs in trans. These results provide the first set of PHAS loci and phasiRNAs in P. notoginseng, and enhance our understanding of PHAS in plants.
Plant small RNAs, with 21 to 24 nucleotides (nt), play crucial roles in a variety of biological processes, including development, stress responses, defense and epigenetic modifications [1, 2]. Based on the origin and biogenesis, small RNA in plants can be divided into two main categories: microRNA (miRNA) and small interfering RNA (siRNA); siRNA can be further divided into three classes: heterochromatic siRNA (hc-siRNAs), natural antisense transcript siRNA (NAT-siRNA) and phased, secondary siRNA (phasiRNA) [3, 4]. The common feature of these small RNA is that they are generated by different DICER-LIKE (DCL) member from their long precursors, and then are incorporated into an effector complex to guide target recognition and subsequent silencing events in a sequence-specific manner [4, 5].
MicroRNAs are derived from their precursors, which are non-coding RNAs and can form hairpin-like secondary structures. The miRNA precursors are generated by RNA polymerase II (Pol II); RNase III enzyme DCL1 is responsible for the biogenesis of the mature miRNAs via processing of miRNA precursors. Most miRNAs negatively regulate their target genes through homolog based mRNA cleavage or translation inhibition at post-transcriptional level [2, 6], however some miRNAs may activate their targets through different mechanisms [7–9]. Small interfering RNAs (siRNAs) are characterized by their biogenesis depending on different RNA-dependent RNA polymerase (RDR) member. hc-siRNAs are 24 nucleotides in length, generated by DCL3 activity . hc-siRNAs play a crucial roles in directing DNA methylation and histone modification in an AGO4-dependent pathway, and regulate target gene expression in transcription level [11, 12]. NAT-siRNAs derives from dsRNA precursors formed by bidirectional transcription of two partially overlapping genes [13–15]. PhasiRNA is a kind of secondary siRNA. Phase simple means these siRNAs are generated precisely in phased pattern initiated at a specific nucleotide. Its biogenesis requires an initiative cleavage on the phasiRNA precursor transcript (PHAS) by a specific miRNA in either “one-hit" or “two-hit" manner, then one of the cleaved product is made double stranded by RDR6 and then catalyzed by DCL4 into 21 nt siRNA in a phased pattern [16–20]. Some of the phasiRNAs may further target their parental genes in cis or other genes in trans .
PhasiRNAs can be generated from either long noncoding RNAs or coding genes. Arabidopsis miR173 (TAS1 and TAS2), miR390 (TAS3), and miR828 (TAS4) can function as guides on non-coding primary transcripts to initiate tasiRNA1-2, tasiRNA3 and tasiRNA4 processing, respectively. These tasiRNAs can further target pentatricopeptide repeat (PPR) family member, auxin response factors (ARFs) and MYB transcription factor in trans manner [17, 19, 22, 23]. Among them, tasiRNA3 is much conversed in land plants [16, 17]. Recently, 21- and 24-nt phasiRNAs derived from long noncoding RNAs have been reported in rice and maize male reproductive organs, they are trigged by miR2118 and miR2275, and cleaved by DCL4 and DCL5 (also known as DCL3b), respectively [20, 24]. A noncoding PHAS locus, triggered by miR4392, was found accumulating preferentially in soybean anthers .
PhasiRNAs, triggered by miRNAs, are also generated from protein-coding loci in many plants, and first found in Arabidopsis . It is worthy to point out that a significant number of PPR, nucleotide binding site-leucine-rich repeat (NB-LRR) and MYB family members are PHAS loci [21, 26–29]. PPR family is one of the largest families in Arabidopsis, containing at least 448 PPR-related genes. Many PPR-P clade transcripts have been identified as PHAS loci triggered by miR161.1, miR161.2, miR400 and tasiRNAs produced by miR173-TAS1/2 [3, 26]. PPR PHAS loci have been found in 9 eudicots, and triggered by miR7122, miR1509 and fve-PPRtri1/2 . MYB transcription factor are target by miR159/miR828/miR858 in apple, peach and lotus, MYB-derived phasiRNAs could target a variety of genes with different function [27, 29, 30]. NB-LRR is the largest family targeted by small RNA. Arikit et al.  found 500 PHAS loci in soybean, amazingly, 208 loci are NB-LRR genes; Zhai et al.  found 114 PHAS loci in Medicago, among them, 79 loci are NB-LRR PHAS loci and 74 NB-LRR loci are targeted by miR1507, miR2109 and miR2118. NB-LRR PHAS loci are found in many plant species, indicating a conserved role of phasiRNA in regulation of NB-LRR [29, 31]. NB-LRR can be targeted by both miRNAs and phasiRNAs . Besides these three gene families, AP2-like gene was found as PHAS loci triggered by miR156 and miR172 in Arabidopsis, DCL2 and Suppressor of Gene Silencing 3 (SGS3) were designed as PHAS loci trigged by miR1515 and miR2118, respectively, in Medicago and soybean .
To investigate phasiRNAs in Panax notoginseng, a precious traditional Chinese herb, we analyzed 24 small RNA profiles of P. notoginseng, and identified 204 PHAS loci generating 21 nt phasiRNAs, and 90 PHAS loci generating 24 nt phasiRNAs. Via integrating miRNAs and parallel analysis of RNA end (PARE) data, we found that some PPR genes generate phasiRNAs, moreover, we found these phasiRNAs could target other PPR loci in trans manner. Using RLM 5’-RACE technique, we further validated the cleavage on PPR transcripts induced by phasiRNAs. We present direct evidence that phasiRNAs originated from PPR genes are functional to regulate other PPR genes in trans. Our results also suggest that miR171, with 21 nt, could trigger phasiRNA generation in wide type plants, although existing study found that miR171 with additional nucleotide at 3’ end could trigger phasiRNAs in a mutant.
Materials and small RNA sequencing profiles
We sequenced 7 small RNA sequencing profiles of P. notoginseng grown in Wenshan County, Yunnan, China (Additional file 1: Table S1). Two of these 7 samples was mixed RNA sample of the leaf and root. Other 5 samples were RNA samples extracted from roots of 5 different plants. The tissue samples collected were stored in liquid nitrogen immediately. The samples were kept at -80 °C until RNAs were extracted. Total RNAs were extracted from samples using the TRIzol reagent (Invitrogen, Thermo Fisher Scientific Inc., USA) according to the manufacturer’s protocol. The integrities of the RNAs were checked using an ultraviolet spectrophotometer (Hoefer, MA, USA), based on the ratio of the optical density at 260 nm to that at 280 nm (OD260/280) and were also assessed by electrophoresis in a denaturing formaldehyde agarose gel, based on visual comparison with the strength of 18S and 28S ribosomal RNAs. The RNA samples with OD260/280 between 1.8 and 2.0 were checked for the total quantities. Only samples with at least 20 μg were chosen for preparation of sRNA sequencing libraries. The small RNAs of the samples were isolated from total RNAs and were sequenced using Illumina HiSeq 2000 sequencer. The 7 obtained small RNA sequencing profiles were deposited into the NCBI GEO database and are accessible with series ID GSE98118.
Because each plant sample may only express some parts of the whole gene set, it is very helpful to identify PHAS loci as completely as possible by analyzing sRNA profiles from multiple samples. Therefore, the 17 small RNA sequencing profiles downloaded from the NCBI SRA database with accession ID SRP082250 and the 7 profiles obtained in this study were combined to get unique sRNA sequences and normalized frequencies in each of the profiles.
Computational steps for identifying PHAS loci
where k was the number of phased unique 21 nt (or 24 nt) sRNAs in the window, n was the total number of unique 21 nt (or 24 nt) sRNAs in the window, and m was the number of phases. In this study, m was set to 10.
We set the window sizes of PHAS loci according to the sizes of reported PHAS loci and the distances between the complementary sites of miRNA triggers to the phasiRNAs in other species. As reported previously [16, 17, 29], the functional tasiRNAs derived from TAS3 are at the seventh or eighth 21 nt segments from the 3’ complementary site of miR390. TAS4 derived functional and conserved tasiRNA, TAS4-siR81(-), is also very close the miR828 complementary site [23, 29]. For other PHAS loci and TAS1/2 in Arabidopsis, it has not been reported that these PHAS loci encode conserved phasiRNAs. But the phasiRNAs from these PHAS and TAS1/2 were very close to complementary sites of their miRNA triggers [17, 21]. Therefore, we set the window sizes of PHAS loci as ten 21 or 24 nt segments.
where m was the number of phases in the window, k was the number of unique phased siRNAs in the window, U i was the number of non-phased reads at the ith phase from the position, and P i was the number of phased reads at the ith phase from the position. m was 10 in this study.
If a window has a P-value less than 0.05, it was extended 100bp at both 5’- and 3’-ends. Then we merged the overlapped windows. We next calculated the false discovery rates with the P-values of the merged windows with the method in . We selected the merged windows with a maximal phase scores of 5 or more and multiple test corrected P-values less than 0.05 as putative PHAS loci. We named the predicted PHAS loci with its transcript ID and a increasing serial number. We reported the phased siRNAs of the predicted PHAS loci as phasiRNAs. We named the phasiRNAs of a PHAS loci by adding siR and a serial number to the name of the PHAS loci.
The degradome of a P. notoginseng root sample in our previous study was retrieved from NCBI SRA database using the accession number SRP087606. Briefly, the root of a P. notoginseng plants grown in a shading greenhouse at Kunming University of Science and Technology (24° 51 ′ 0 ′′ N, 102° 52 ′ 2 ′′ E, altitude 1835 m), Kunming, Yunnan, China, was frozen in liquid nitrogen immediately after harvesting. During the experiment, the daily average temperature was 25 °C, the daily maximum difference in temperature was 10 °C and humidity was 60-80%. The total RNA from the root of a P. notoginseng plant was extracted using the TRIzol reagent (Invitrogen, Thermo Fisher Scientific Inc., USA) according to the manufacturer’s protocol. The integrity of the RNA was checked with an ultraviolet spectrophotometry (Hoefer, MA, USA) and 2100 BioAnalyzer (Agilent Technologies, Santa Clara, CA, USA). The degradome of polyadenylated transcripts was sequenced using Illumina HiSeq 4000 sequencer. The obtained degradome sequencing profile was filtered to remove low quality reads that have low scored nucleotides (< 20). Then, the 3’ adapters in the remaining reads were removed.
Identifying miRNA complementary sites on PHAS and phasiRNA targets
The degradome sequencing profile were downloaded from the NCBI SRA database with the accession number SRP087606. The raw degradome sequencing profile was filtered to remove low quality reads that have low scored nucleotides (< 20). Then, the 3’ adapters in the remaining reads were removed. The unique sequences were obtained and the frequencies of the unique sequences were calculated. The SeqTar algorithm  was used to predict miRNA complementary sites on the original transcripts of PHAS loci. For conserved miRNAs, the targets that have at least one valid read, i.e., read started at the 9th to 11th positions of a miRNA binding site (as defined in ), or targets that have less than 4 mismatches were used for further analysis.
We identified the targets of phasiRNAs with the SeqTar algorithm . Only targets that have at least one valid read and less than 4 mismatches were used for further analysis.
RLM 5’-RACE validations of some phasiRNA targets
Modified 5’-RACE assay was performed using the GeneRacer Kit (Invitrogen) to validate the predicted targets. Briefly, 5 μg total RNA was ligated with GeneRacer T M RNA Oligo (5’-CGACUGGAGCACGAGGACACUGACAUGGACUGAAGGAGUAGAAA-3’) and reverse transcription was performed using SuperScript T M III Reverse Transcriptase and oligo dT(18) primer. The resulting cDNA was used as template for PCR amplification with GeneRacer 5’ primer (5’-CGACTGGAGCACGAGGACACTGA-3’) and a gene-specific reverse primer. A second nested PCR was performed using GeneRacer 5’ nested primer (5’-GGACACTGACATGGACTGAAGGAGTA-3’) and a gene specific nested primer. The amplified products were run on a 2% agarose gel, bands with expected size were purified using QIAquick® Gel Extraction Kit (QIAGEN) and ligated onto T Vector pMD 19 (Simple)(Takara), following transformation and clone PCR, plasmids were isolated and subjected to Sanger sequencing. Gene specific primers used are listed in Additional file 1: Table S2.
Results and discussion
Samples and small RNA profiles of P. notoginseng
We collected 5 root samples and 2 whole plant samples of P. notoginseng grown in Wenshan County, Yunnan, China. The small RNAs of the samples were isolated from total RNAs and were sequenced using Illumina HiSeq 2000 sequencer (see Methods). The obtained sRNA profiles were combined with 17 root sRNA profiles from NCBI SRA database obtained with accession number SRP082250 to get unique sRNA sequences and their frequencies in each of the profiles.
Twenty-one nt PHAS loci in P. notoginseng
The combined sRNA profiles of 23 samples were aligned to a self-assembled transcriptome database of P. notoginseng using SOAP2 . Then, we used a self-written program to predict PHAS loci and phasiRNAs from the alignment result of SOAP2 based on methods proposed previously [33, 34] (see Materials and methods for details). Then we merged overlapped PHAS loci. As listed in Additional file 1: Table S3, we totally identified 204 PHAS loci encoding 21 nt phasiRNAs with a combined criterion of a phase score ≥ 5 and a multiple-test corrected P<0.05. These PHAS loci produced at least 2505 twenty-one nt phasiRNAs that were sequenced in the sRNA profiles used in the study (Additional file 1: Table S4).
miR171l, with 21 nt, triggers two PHAS loci from its conserved targets in the scarecrow-like transcription factor family (Additional file 1: Table S3). One of these two PHAS loci, P21_TR101805_c1_g7_i2_1, is shown in Fig. 2d. Together with the degradome profile (Fig. 2h), it is shown that miR171l triggers the generation of phasiRNAs from the center of its complementary site on P21_TR101805_c1_g7_i2_1.
More over, miR482 also triggers phasiRNAs by targeting an unknown gene and a hypothetical protein (Additional file 2: Figures S1i and S1j, respectively).
The miRNA triggers that are reported to target PHAS loci
 and this study
apple, peach, lotus
miR2118, miR482, miR1507,
legume, cotton, lotus,
tomato, P. notoginseng
and this study
miR7122, miR1509, miR173
Twenty-four nt PHAS loci in P. notoginseng
PhasiRNA targets in P. notoginseng
We used the SeqTar algorithm  to identify targets of phasiRNAs by using a degradome profile of the P. notoginseng root. By keeping targets with at least 1 valid degradome read and have less than 4 mismatches, we found 5062 targets for 1699 twenty one nt phasiRNAs (Additional file 1: Table S7). In addition to 98 cis-targets, 4964 of these 5062 targets are trans-targets (Additional file 1: Table S7).
TR102419 |c3_g16_i2 (an unknown gene), TR99791 |c3_g3_i1 (a PPR gene), TR102500 |c1_g1_i1 (a ribosomal protein), TR83223 |c0_g1_i1 (a heat shock protein), TR82797 |c0_g1_i1 (ubiquitin-60S ribosomal protein) and TR8286 |c0_g1_i1 (a ribosomal protein) are also identified as targets of different phasiRNAs (Fig. 4d to i). Interestingly, TR102419 |c3_g16_i2 and TR99791 |c3_g3_i1 produce the same phasiRNA and may form a mutual regulation loop or self-regulation loops through phasiRNAs (Fig. 4d and e).
In soybean and Medicago truncatula, miR1509 indirectly initiates the generation of phasiRNAs from PPR genes by one or two-layers of TASL-tasiRNA interactions . In comparison, our results suggest that in addition to miR1509/TASL/PPR pathway, some PPR genes not triggered by miR1509 could also generate phasiRNAs to repress other PPRs in P. notoginseng. Presumably, these PPR PHAS loci are triggered by other unknown miRNAs or siRNAs.
One hundred twenty-nine targets of sixty one 24 nt phasiRNAs were also identified and filtered with the same criteria as those for the 21 nt phasiRNAs (Additional file 1: Table S8). Thirteen and 116 of these targets are cis- and trans-targets, respectively. But base on the degradome, only a few 24 nt phasiRNAs induce clear cleavage on the targets (Additional file 1: Table S8), suggesting these 24 nt phasiRNAs might function in other means than cleaving targets in trans, as noticed previously .
The publication costs of this article were supported by two grants (No. 31460295 and 31760314) of National Natural Science Foundation of China (http://www.nsfc.gov.cn/) and a grant (No. SKLGE-1511) of the Open Research Funds of the State Key Laboratory of Genetic Engineering, Fudan University, China, to YZ; grants (No. 2013FA031 and 2014FA003) of the Science and Technology Bureau of Yunnan Province to XC; and a starting-up grant of Henan Normal University (No. 5101049470095) and Program for Innovative Research Team in Science and Technology in University of Henan Province (No. 17IRTSTHN017) to YL. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
The 7 obtained small RNA sequencing profiles were deposited into the NCBI GEO database and are accessible with series accession ID GSE98118.
About this supplement
This article has been published as part of BMC Genomics Volume 19 Supplement 1, 2018: 16th International Conference on Bioinformatics (InCoB 2017): Genomics. The full contents of the supplement are available online at https://bmcgenomics.biomedcentral.com/articles/supplements/volume-19-supplement-1.
YZ, YL, and XC conceived of and designed the research. YZ designed the computational pipeline for analyzing small RNA sequencing profiles and degradome profile. KC, LL, XZ predicted PHAS and phasiRNAs from the small RNA profiles. SR, JG, and QW predicted phasiRNA targets using the degradome profile. PL, SL and XC collected the plant samples and extracted RNAs for high throughput sequencing. KC, LL, XZ and YZ analyzed the results. YL and YY performed the RLM 5’-RACE experiments. YZ and YL wrote the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Chen X. Small RNAs in development–insights from plants. Curr Opin Genet Dev. 2012; 22(4):361–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Sunkar R, Li Y-F, Jagadeeswaran G. Functions of microRNAs in plant stress responses. Trends Plant Sci. 2012; 17(4):196–203.View ArticlePubMedGoogle Scholar
- Fei Q, Xia R, Meyers BC. Phased, secondary, small interfering RNAs in posttranscriptional regulatory networks. Plant Cell. 2013; 25(7):2400–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Axtell MJ. Classification and comparison of small RNAs from plants. Annu Rev Plant Biol. 2013; 64:137–59.View ArticlePubMedGoogle Scholar
- Axtell MJ, Snyder JA, Bartel DP. Common functions for diverse small RNAs of land plants. Plant Cell. 2007; 19(6):1750–69.View ArticlePubMedPubMed CentralGoogle Scholar
- Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004; 116(2):281–97.View ArticlePubMedGoogle Scholar
- Vasudevan S, Tong Y, Steitz JA. Switching from repression to activation: microRNAs can up-regulate translation. Science. 2007; 318(5858):1931–4.View ArticlePubMedGoogle Scholar
- Vasudevan S. Posttranscriptional upregulation by microRNAs. Wiley Interdiscip Rev: RNA. 2012; 3(3):311–30.View ArticlePubMedGoogle Scholar
- Xiao M, Li J, Li W, Wang Y, Wu F, Xi Y, Zhang L, Ding C, Luo H, Li Y, et al.MicroRNAs activate gene transcription epigenetically as an enhancer trigger. RNA Biol. 2017; 14(10):1326–34.View ArticlePubMedGoogle Scholar
- Matzke M, Kanno T, Daxinger L, Huettel B, Matzke AJ. RNA-mediated chromatin-based silencing in plants. Curr Opin Cell Biol. 2009; 21(3):367–76.View ArticlePubMedGoogle Scholar
- Zilberman D, Cao X, Johansen LK, Xie Z, Carrington JC, Jacobsen SE. Role of Arabidopsis ARGONAUTE4 in RNA-directed DNA methylation triggered by inverted repeats. Curr Biol. 2004; 14(13):1214–20.View ArticlePubMedGoogle Scholar
- Chan SW-L, Zilberman D, Xie Z, Johansen LK, Carrington JC, Jacobsen SE. RNA silencing genes control de novo DNA methylation. Science. 2004; 303(5662):1336–1336.View ArticlePubMedGoogle Scholar
- Borsani O, Zhu J, Verslues PE, Sunkar R, Zhu J-K. Endogenous siRNAs derived from a pair of natural cis-antisense transcripts regulate salt tolerance in Arabidopsis. Cell. 2005; 123(7):1279–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Han MH, Goud S, Song L, Fedoroff N. The Arabidopsis double-stranded RNA-binding protein HYL1 plays a role in microRNA-mediated gene regulation. Proc Natl Acad Sci U S A. 2004; 101(4):1093–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang X, Xia J, Lii YE, Barrera-Figueroa BE, Zhou X, Gao S, Lu L, Niu D, Chen Z, Leung C, Wong T, Zhang H, Guo J, Li Y, Liu R, Liang W, Zhu JK, author WZ, Jin H. Genome-wide analysis of plant nat-siRNAs reveals insights into their distribution, biogenesis and function. Genome Biol. 2012; 13(3):20.View ArticleGoogle Scholar
- Axtell MJ, Jan C, Rajagopalan R, Bartel DP. A two-hit trigger for siRNA biogenesis in plants. Cell. 2006; 127(3):565–77.View ArticlePubMedGoogle Scholar
- Allen E, Xie Z, Gustafson AM, Carrington JC. microRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell. 2005; 121(2):207–21.View ArticlePubMedGoogle Scholar
- Allen E, Howell MD. miRNAs in the biogenesis of trans-acting siRNAs in higher plants. Semin Cell Dev Biol. 2010; 21(8):798–804.View ArticlePubMedGoogle Scholar
- Yoshikawa M, Peragine A, Park MY, Poethig RS. A pathway for the biogenesis of trans-acting siRNAs in Arabidopsis. Genes Dev. 2005; 19(18):2164–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Song X, Li P, Zhai J, Zhou M, Ma L, Liu B, Jeong DH, Nakano M, Cao S, Liu C, et al. Roles of DCL4 and DCL3b in rice phased small RNA biogenesis. Plant J. 2012; 69(3):462–74.View ArticlePubMedGoogle Scholar
- Zhai J, Jeong D-H, De Paoli E, Park S, Rosen BD, Li Y, González AJ, Yan Z, Kitto SL, Grusak MA, et al.MicroRNAs as master regulators of the plant NB-LRR defense gene family via the production of phased, trans-acting siRNAs. Genes Dev. 2011; 25(23):2540–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Vazquez F, Vaucheret H, Rajagopalan R, Lepers C, Gasciolli V, Mallory AC, Hilbert JL, Bartel DP, Crété P. Endogenous trans-acting siRNAs regulate the accumulation of Arabidopsis mRNAs. Mol Cell. 2004; 16(1):69–79.View ArticlePubMedGoogle Scholar
- Rajagopalan R, Vaucheret H, Trejo J, Bartel DP. A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev. 2006; 20(24):3407–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhai J, Zhang H, Arikit S, Huang K, Nan G-L, Walbot V, Meyers BC. Spatiotemporally dynamic, cell-type–dependent premeiotic and meiotic phasiRNAs in maize anthers. Proc Natl Acad Sci. 2015; 112(10):3146–151.View ArticlePubMedPubMed CentralGoogle Scholar
- Arikit S, Xia R, Kakrana A, Huang K, Zhai J, Yan Z, Valdés-López O, Prince S, Musket TA, Nguyen HT, Staceyc G, Meyers BC. An atlas of soybean small RNAs identifies phased siRNAs from hundreds of coding genes. Plant Cell. 2014; 26(12):4584–601.View ArticlePubMedPubMed CentralGoogle Scholar
- Howell MD, Fahlgren N, Chapman EJ, Cumbie JS, Sullivan CM, Givan SA, Kasschau KD, Carrington JC. Genome-Wide Analysis of the RNA-DEPENDENT RNA POLYMERASE6/DICER-LIKE4 Pathway in Arabidopsis Reveals Dependency on miRNA- and tasiRNA-Directed Targeting. Plant Cell. 2007; 19(3):926–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhu H, Xia R, Zhao B, An YQ, Dardick CD, Callahan AM, Liu Z. Unique expression, processing regulation, and regulatory network of peach (Prunus persica) miRNAs. BMC Plant Biol. 2012; 12(1):149.View ArticlePubMedPubMed CentralGoogle Scholar
- Xia R, Meyers BC, Liu Z, Beers EP, Ye S, Liu Z. MicroRNA superfamilies descended from miR390 and their roles in secondary small interfering RNA biogenesis in eudicots. Plant Cell. 2013; 25(5):1555–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Zheng Y, Wang S, Sunkar R. Genome-wide discovery and analysis of phased small interfering RNAs in chinese sacred lotus. PLoS ONE. 2014; 9(12):113790.View ArticleGoogle Scholar
- Xia R, Zhu H, An Y-q, Beers EP, Liu Z. Apple miRNAs and tasiRNAs with novel regulatory networks. Genome Biol. 2012; 13(6):47.View ArticleGoogle Scholar
- Shivaprasad PV, Chen H-M, Patel K, Bond DM, Santos BA, Baulcombe DC. A microRNA superfamily regulates nucleotide binding site–leucine-rich repeats and other mRNAs. Plant Cell. 2012; 24(3):859–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Li R, Yu C, Li Y, Lam TW, Yiu S-M, Kristiansen K, Wang J. SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009; 25(15):1966–7.View ArticlePubMedGoogle Scholar
- Chen HM, Li Y-H, Wu SH. Bioinformatic prediction and experimental validation of a microRNA-directed tandem trans-acting siRNA cascade in Arabidopsis. PNAS. 2007; 104(9):3318–23.View ArticlePubMedPubMed CentralGoogle Scholar
- De Paoli E, Dorantes-Acosta A, Zhai J, Accerbi M, Jeong DH, Park S, Meyers BC, Jorgensen RA, Green PJ. Distinct extremely abundant siRNAs associated with cosuppression in petunia. RNA. 2009; 15(11):1965–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995; 57(1):289–300.Google Scholar
- Zheng Y, Li Y-F, Sunkar R, Zhang W. SeqTar: an effective method for identifying microRNA guided cleavage sites from degradome of polyadenylated transcripts in plants. Nucleic Acids Res. 2012; 40(4):28–28.View ArticleGoogle Scholar
- Zhai J, Zhao Y, Simon SA, Huang S, Petsch K, Arikit S, Pillay M, Ji L, Xie M, Cao X, Yu B, Timmermans M, Yang B, Chen X, Meyers BC. Plant microRNAs display differential 3’ truncation and tailing modifications that are ARGONAUTE1 dependent and conserved across species. Plant Cell. 2013; 25(7):2417–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Zheng Y, Li T, Xu Z, Wai CM, Chen K, Zhang X, Wang S, Ji B, Ming R, Sunkar R. Identification of microRNAs, phasiRNAs and their targets in pineapple. Trop Plant Biol. 2016; 9(3):176–86.View ArticleGoogle Scholar
- Fei Q, Li P, Teng C, Meyers BC. Secondary siRNAs from Medicago NB-LRRs modulated via miRNA–target interactions and their abundances. Plant J. 2015; 83(3):451–65.View ArticlePubMedGoogle Scholar
- Zhu QH, Fan L, Liu Y, Xu H, Llewellyn D, Wilson I. miR482 regulation of NBS-LRR defense genes during fungal pathogen infection in cotton. PloS ONE. 2013; 8(12):84390.View ArticleGoogle Scholar