Global reach, higher impact
Materials and methods
Int J Biol Sci 2015; 11(7):781-793. doi:10.7150/ijbs.11619
Identification and Characterization of Novel Maize Mirnas Involved in Different Genetic Background
Key Laboratory of Crop Biology of Anhui Province, Anhui, Agricultural University, Hefei 230036, China
This is an open access article distributed under the terms of the Creative Commons Attribution (CC BY-NC) License. See http://ivyspring.com/terms for full terms and conditions.
How to cite this article:
Sheng L, Chai W, Gong X, Zhou L, Cai R, Li X, Zhao Y, Jiang H, Cheng B. Identification and Characterization of Novel Maize Mirnas Involved in Different Genetic Background. Int J Biol Sci 2015; 11(7):781-793. doi:10.7150/ijbs.11619. Available from http://www.ijbs.com/v11p0781.htm
MicroRNAs (miRNAs) are a class of small, non-coding regulatory RNAs that regulate gene expression by guiding target mRNA cleavage or translational inhibition in plants and animals. At present there is relatively little information regarding the role of miRNAs in the response to drought stress in maize. In this study, two small RNA libraries were sequenced, and a total of 11,973,711 and 14,326,010 raw sequences were generated from growing leaves of drought-tolerant and drought-sensitive maize seedlings, respectively. Further analysis identified 192 mature miRNAs, which include 124 known maize (zma) miRNAs and 68 potential novel miRNA candidates. Additionally, 167 target genes (259 transcripts) of known and novel miRNAs were predicted to be differentially expressed between two maize inbred lines. Of these, three novel miRNAs were up-regulated and two were down-regulated under drought stress. The expression of these five miRNAs and nine target genes was confirmed using quantitative reverse transcription PCR. The expression of three of the miRNAs and their putative target genes exhibited an inverse correlation, and expression analysis suggested that all five may play important roles in maize leaves. Finally, GO annotations of the target genes indicated a potential role in photosynthesis, may therefore contribute to the drought stress response. This study describes the identification and characterization of novel miRNAs that are the differentially expressed in drought-tolerant and drought-sensitive inbred maize lines. This provides the foundation for further investigation into the mechanism of miRNA function in response to drought stress in maize.
Keywords: Maize, Drought stress, MicroRNA, Target genes, High-throughput sequencing, qRT-PCR
MicroRNAs (miRNAs) are small, endogenous RNAs that regulate gene expression in plants and animals [1-5] at the post-transcriptional level by translational repression or target degradation and gene silencing. Plant miRNAs function in diverse processes associated with growth and development such as developmental transitions, organ polarity, auxin signaling, boundary formation or organ separation, leaf and stem growth, floral organ identity and reproductive development . Several miRNAs are involved in the regulation of root development in Arabidopsis, including miR167 that targets auxin response factors (ARFs) in order to regulate the emergence of shoot-borne roots [7-9], and miR164 that directly cleaves transcription factor NAC1 which inhibits lateral root development. Other miRNAs such as miR156 and miR172 affect flowering time when over-expressed in Arabidopsis [9-12], and miR172 exhibits a similar temporal expression pattern in maize, where it targets Glossy15, a gene required for the expression of juvenile epidermal traits . Moreover, miR395 regulates sulfur assimilation and translocation by adjusting the mRNA levels of ATP sulfurylase and a low affinity sulfur transporter [14, 15], while miR399 regulates phosphate homeostasis in Arabidopsis by suppressing the ubiquitin-conjugating E2 enzyme, PHO2 [16, 17]. In addition, accumulating data demonstrate that miRNAs play major roles in response to biotic and abiotic stress [18, 19, 20], such as low nitrate condition [21, 22], low phosphorus stress , salt stress  and so on.
Bioinformatics approaches have been successful for the discovery of conserved miRNAs in numerous plant species including cotton , maize , caster , sorghum , and soybean . At the time of writing, there are 18,226 mature miRNAs deposited in miRBase . As an important model system for basic biological research, maize has contributed significantly to our understanding of plant development and evolution, and this knowledge has been employed to elucidate the regulatory functions of miRNA genes . The total number of miRNAs in an individual organism has been estimated to represent approximately 1% of the total number of coding genes [31, 32]. To date, only 150 genes within 26 miRNA families have been identified in maize , and many more remain to be discovered.
Maize is one of the most important food crops worldwide, and is also used for animal feed, silage and in industrial products. Drought stress is often the limiting factor for plant growth and agricultural productivity, and this is especially true for maize. Plant microRNAs have been shown to play important functions in plant growth and development. However, role of miRNAs in the response to drought stress is poorly understood in maize. In recent years, high-throughput sequencing technologies and bioinformatics approaches have together identified many novel miRNAs [33, 34, 47]. In the present work, drought-tolerant (Hz4) and drought-sensitive (3189) inbred maize lines  were investigated using Illumina sequencing to identify novel drought-responsive miRNAs. Expression of the identified miRNAs was confirmed using qRT-PCR and target genes were elucidated. Probing of two small RNA libraries from young leaves led to the identification of 124 known maize (zma) miRNAs, and a further 68 novel maize miRNAs were discovered. Of these, five novel miRNAs that were differentially expressed in the two inbred lines were characterized, and the results suggest they may play an important role in response to drought stress.
Materials and methods
Plant materials and drought stress treatment
Drought-tolerant (Hz4) and drought-sensitive (3189) maize inbred lines were germinated in a greenhouse and grown under standard conditions until seedlings had developed three leaves. Seedlings were then subjected to drought conditions. Specifically, the soil water content was reduced from 90% to 70%, and approximately 5 g of leaves was collected from each seedling and immediately frozen in liquid nitrogen.
Construction of small RNA libraries for high-throughput sequencing
To construct small RNA libraries, total RNA was extracted from leaves using Trizol reagent (Invitrogen, USA) according to the manufacturer's protocol as modified by LC Sciences (Houston, USA). The small RNA fraction (10-40 nt) was isolated by polyacrylamide gel electrophoresis and ligated with proprietary adaptors (Illumina). Short RNAs were then reverse-transcribed into cDNA by RT-PCR. Small RNA libraries were constructed for each maize line and sequenced using a Genome Analyzer GA-IIX (Illumina) following the manufacturer's instructions. Small RNA sequencing data was processed with the in-house program ACGT101-miR (LC Sciences, USA). Program parameters were set as previously described  with modifications to adapt the program for plant miRNA prediction.
Sequencing data processing
To begin with, we compared the filtered small RNA sequences with known miRNAs in miRBase 20.0 [http://www.mirbase.org  and a 0-2 nucleotide mismatch was considered. Meanwhile, the filtered small RNA sequences were compared with ZmB73_4a.53 (ZmB73_AGPv1)  using BLASTn and no mismatches were permitted. Then ESTs that contained small RNA sequences were selected and their secondary structures were predicted using MFOLD (http://mfold.rna.albany.edu/?q=mfold/RNA-Folding-Form) for identifying potential novel miRNAs. At last, small RNAs with corresponding star (miRNA*) sequences were counted to be novel miRNAs according to the criteria for miRNA definition .
Prediction and validation of miRNA target genes
The target genes of miRNAs were predicted using the method of target prediction proposed by Allen et al. (2005)  and Schwab et al. (2005) . The plant small RNA target analysis server psRNATarget  was employed for prediction through searching against the maize full-length cDNA sequences. In addition, the rules of target prediction were referred to universal criteria previously suggested [40, 41] and the parameters of psRNATarget server were the default.
Validations of 18 novel miRNAs which low abundance and more than 24 nt in length were carried out by stem-loop reverse transcription-PCR (RT-PCR). The stem-loop RT-PCR was using the Applied Bio systems 7300 Real-Time PCR System (Applied Bio systems, Foster City, CA) and all reactions were run in triplicate. 18S rRNA was used as the internal control for stem-loop RT-PCR. In addition, quantitative real-time RT-PCR (qRT-PCR) was performed to validate the target genes of five novel miRNAs that were differentially expressed between the two inbred lines using three biological replicates. Total RNA from the leaves of the two maize inbred lines was extracted using Trizol reagent (Invitrogen, USA) according to the manufacturer's protocol as modified by LC Sciences (Houston, USA). Samples were collected at the same time as those for miRNAs sequencing. The RNase-free DNase (Promega, USA) was used to remove contaminating DNA from total RNA and first strand cDNA was synthesized according to the manufacturer's instructions (Takara, Japan). The PCR system contained 2 μL cDNA, 12.5 μL Platinum SYBR Green qPCR SuperMix-UDG (Takara, Japan), 2 μL 10 μM primer mixture and 8.5 μL RNase-free water. Sample cycle threshold (Ct) values were standardized for each template based on a GAPDH control reaction and the comparative Ct method (2-ΔΔCt) was used to determine the relative transcript abundance of each gene . Gene accession and primer sequences are listed in Additional File 3: Table S4.
Gene ontology (GO) analysis
Target genes were subjected to Gene Ontology (GO) analysis to uncover the miRNA-gene regulatory network on the basis of biological process and molecular function. GO annotations were extracted using AgriGO (http://bioinfo.cau.edu.cn/agriGO/) and included alignments of individual transcripts representing diverse functions, processes or components. A graphical representation of the GO ontology distribution of the biological processes, cellular component and molecular function categories can be found in Online Resource 8.
Overview of sequencing results
Two small RNA libraries constructed from the leaves of drought-tolerant and drought-sensitive maize seedlings were sequenced using high-throughput methods, and a total of 11,973,711 and 14,326,010 raw sequences were generated from the drought-tolerant and drought-sensitive libraries, respectively. After removing low quality reads and clipping adaptor sequences, 8,983,222 and 7,948,375 clean reads were remained with lengths ranging from 17 to 44 nt (Table 1). The reads of 34 nt in length were the most abundant among total reads of each small RNA library (Figure 1A). In addition, a length of 33 nt was the most abundant among unique reads of SL_HZ4 and a length of 24 nt was the most abundant among unique reads of SL_3189 (Figure 1B). To identify known miRNAs, small RNA sequences were used with BLASTn to search the miRBase database (release 20.0, June 2013). A total of 124 known maize (zma) miRNAs were identified, of which 76 were detected in Hz4 and 124 were detected in 3189. These miRNAs belong to different miRNA families such as the miR156, miR166, and miR167 families (Figure 2). The miR156 family was the most represented, indicating a higher level of expression in maize seedlings, as was previously reported . Uridine (U) was the most common nucleotide at the 5′ end (>65%), whereas this nucleotide was comparatively rare at positions 2-4 (Figure 3).
Summary of small RNA sequencing.
Overview of reads from raw data to cleaned sequences.
3ADT&length filter: reads removed due to 3ADT not found and length with <17 nt and >25 nt were removed (for plants); length with<16 and >29 were remove(for animals)
Junk reads:Junk: >=2N, >=7A, >=8C, >=6G, >=7T, >=10Dimer, >=6Trimer, or >=5Tetramer
Rfam:Collection of many common non-coding RNA families except micro RNA; http://rfam.janelia.org
Repeats:Prototypic sequences representing repetitive DNA from different eukaryotic species; http://www.girinst.org/repbase.
Notes:There is overlap in mapping of reads with mRNA, rRNA, tRNA, snRNA, snoRNA and repeats.
Length distribution of small RNAs from sequencing of the two inbred maize libraries. (A) Size distribution of total sequences. (B) Size distribution of unique sequences.(Click on the image to enlarge.)
Abundance of conserved miRNA families in the two inbred maize libraries.(Click on the image to enlarge.)
Relative nucleotide bias at each position of the known maize miRNAs.(Click on the image to enlarge.)
Identification of novel miRNAs in the maize genome
A computational prediction method was employed to identify novel miRNAs in the maize genome. After excluding sRNAs that matched known miRNAs, rRNAs, tRNAs, snRNAs, and snoRNAs, as well as those overlapping protein-coding genes, remaining sRNAs exhibiting a perfect match to maize genomic sequences were used to predict potential novel miRNAs. The selection of novel miRNAs was based on the ability of flanking sequences to fold back in a hairpin structure . Predicted secondary structures of the flanking sequences were generated using Mfold. These putative miRNA precursors were then used to extract miRNA*s, which are considered strong evidence for DICER-LIKE-1 (DCL1)-derived products . A total of 68 regions satisfied these criteria and were designated as novel miRNA gene candidates (Table 2). The lengths of novel miRNAs varied from 20 to 37 nt, with 24 nt the most abundant length (18 out of 68), which is typical of DCL digestion products. Of these, 24 were found to be expressed in both maize inbred lines, with 10 weakly expressed, as indicated by a normalized sequencing read of less than 1 in both samples. In addition, the sequence frequency of most novel miRNA candidates was extremely low. It has been reported that miRNAs may be expressed at very low levels and may be present only in particular cell types and/or under particular circumstances . None of the novel miRNAs were previously reported in other species and are therefore likely to be maize-specific. Then 18 novel miRNAs which low abundance and more than 24 nt in length were selected for validation by stem-loop real-time PCR (Additional File 1: Table S1). According to our RT-PCR analysis, 4 novel miRNAs (PC-3p-201205, PC-3p-359801, PC-5p-864621, zma-MIR166b-p5) had no expression quantity in both the leaves of two maize inbred lines. In addition, 5 novel miRNAs (PC-3p-201205, PC-5p-1164336, zma-MIR164e-p5, zma-MIR166b-p5, zma-MIR393a-p3) had no expression quantity in the leaves of Hz4 and 1 novel miRNA (PC-5p-864621) had no expression quantity in the leaves of 3189, which was consistent with the results of the high-throughput sequencing (Additional File 1: Table S1). Expression levels of a few members—PC-3p-129630, PC-3p-420698, PC-3p-552502, PC-3p-564523, PC-3p-793235, zma-MIR164e-p3, zma-MIR166n-p3, zma-MIR169i-p3 and zma-MIR393c-p3—were slightly different (Figure 4A). This difference may be due to variations in Hz4 and 3189 sampling times, or differences in sensitivity and specificity of the two technologies. Expression patterns of the remaining novel miRNAs (PC-3p-1072330, PC-5p-1134476, PC-5p-1164336, zma-MIR164e-p5, zma-MIR393a-p3) obtained by stem-loop real-time RT-PCR were in accordance with the sequencing data (Figure 4B).
Validation of 18 novel miRNAs which low abundance and more than 24 nt in length using stem-loop real-time PCR. (A) different expression levels of miRNAs compared with the sequencing data. (B) similar expression levels of miRNAs compared with the sequencing data. Ordinates indicate relative expression levels.(Click on the image to enlarge.)
Novel maize miRNAs identified in this study.
Clustering of miRNAs
Clusters of miRNAs are often present in the genome where multiple miRNAs are aligned in the same orientation and transcribed as a polycistronic structure to facilitate cooperative action . A total of 69 of the miRNAs identified in this research were located within 22 miRNA clusters, with chromosomes 1, 2, 3, 4, 5, 6, 7, 8 and 10 containing 5, 1, 5, 1, 4, 2, 1, 2 and 1 clusters, respectively (Table 3). The largest miRNA cluster contained 7 miRNAs, all of which belonged to the same miRNA family (miRNA395). Some miRNAs in the same cluster may be transcribed by the same primary miRNA, although their expression levels can differ, as was observed with PC-3p-1105743/PC-3p-897448, PC-3p-552502/PC-5p-62968 and PC-5p-1164336/PC-3p-104764. This is presumably due to these miRNAs undergoing different transcriptional regulation during maturation.
Differentially expressed miRNAs
To identify miRNAs involved in the response to drought stress, differential expression in the two libraries was estimated from the read counts of the high-throughput sequencing. In the two libraries, miRNAs exhibiting a log2 (Hz4/3189) fold change higher than 1 were designated up-regulated, while a log2 (Hz4/3189) fold change less than -1 was designated as down-regulated. In total, 29 miRNAs were differentially expressed between the two inbred lines, with 18 down-regulated and 11 up-regulated (Table 4). Among these differentially expressed miRNAs, five novel miRNAs were selected for expression analysis using qRT-PCR. These were PC-3p-190, PC-3p-104764, PC-3p-129630, PC-3p-552502 and PC-5p-139812. Of these, two were up-regulated and three were down-regulated. Secondary structures were predicted to be of a typical hairpin shape, which is different from siRNAs and piRNAs (Figure 5). The expression patterns of these five miRNAs, obtained by qRT-PCR experiments, were consistent with the results of the high-throughput sequencing reads (Figure 6). However the fold changes obtained from the qRT-PCR data were much lower than those estimated from the high-throughput sequencing data, presumably due to differences in sensitivity and specificity between the experimental approaches .
Identification and confirmation of miRNA targets using qRT-PCR
In most cases, the identified miRNAs were predicted to cleave two or more different targets. With a specific value of 0.85 and a signal-to-noise value of 6.28, a total of 167 genes were predicted to be potential targets for 23 of the 29 differentially expressed miRNAs. Among these, nine genes were predicted to be targets of the five novel miRNAs (Table 5). Interestingly, PC-3p-104764, PC-3p-129630, PC-3p-552502 and PC-5p-139812 shared a common target (GRMZM2G448344), and PC-3p-104764 and PC-3p-129630 shared three common targets (GRMZM2G360821, GRMZM2G385635, GRMZM2G308907), suggesting differnet miRNAs may function in combination in gene regulatory networks.
To confirm whether the miRNAs regulated their potential targets under drought stress conditions, the expression patterns of putative target genes were analyzed using qRT-PCR. It is known that miRNAs regulate gene expression at the post-transcriptional level by translational repression or target degradation and gene silencing. As expected, an inverse relation was observed between the expression of three miRNAs (PC-3p-190, PC-3p-552502 and PC-5p-139812) and their putative target genes, suggesting miRNA-mediated regulation of their potential targets (Figure 6A). The other two novel miRNAs (PC-3p-104764 and PC-3p-129630) did not exhibit this inverse relationship with their putative target genes (Figure 6B). On account of the regulatory mechanism between miRNAs and their target genes is so complicated that one target gene may be correspond to multiple miRNAs. Therefore, we speculated that these two novel miRNAs had little effect on their target genes and the expression of the target genes might be affected by other regulators.
GO functional and transcriptomics analysis of target genes
The results of GO analysis demonstrated that most of the target genes were involved in cellular or metabolic processes such as response to stimuli or biological regulation among other categories (Table 5). Of the nine target genes characterized, seven could be classified into 21 biological processes including carbon fixation, photosynthesis and proteolysis, while the remaining two were of unknown function. For example, PC-3p-190 targets a gene involved in RNA-binding and PC-3p-129630 targets a gene involved in carbon fixation. PC-3p-552502 targets a gene involved in the light reaction of photosynthesis and the target genes of PC-3p-104764 are ATP-dependent peptidase and serine-type endopeptidase. Interestingly, the common target (GRMZM2G448344) of four of the five novel miRNAs, and the common targets (GRMZM2G360821 and GRMZM2G308907) of PC-3p-104764 and PC-3p-129630 are all involved in photosynthesis.
To gain more insight into the expression patterns of the target genes of the five novel miRNAs, a comprehensive expression analysis was performed using the publicly available maize transcriptome data. A total of 18 tissues from five organs have been analyzed using RNA sequencing (RNA-Seq) . Distinct expression profiles were identified for nine target genes from the transcriptome data (Additional File 2: Table S2). Transcriptome data were imported into R and Bioconductor (http://www.bioconductor.org/) for expression analysis to generate a heat map (Figure 7). Based on the heat map, all nine target genes were appeared to be expressed only in the leaves. GRMZM2G360821 and GRMZM2G308907 exhibited the highest expression levels, suggesting important roles for these target genes. These results indicated that these genes may participate in photosynthesis, which was consistent with the GO analysis.
Genome location clusters of pre-miRNAs.
miRNAs differentially expressed in the two maize inbred lines.
Gene Ontology (GO) analysis of potential targets of the five novel miRNAs.
Secondary structure of five novel miRNA precursors. Mature miRNA sequences are shown in yellow.(Click on the image to enlarge.)
Quantitative real-time RT-PCR analysis of five novel miRNAs and their target genes. Expression levels of miRNAs were normalized against 18S rRNA. Fold changes in expression level were estimated using the 2-ΔΔCT method. Data are reported as mean ± SE for three independent experiments. (A) The inverse relationship between three novel miRNAs (PC-3p-190, PC-3p-552502 and PC-5p-139812) and their putative target genes. (B) The uniform relationship between two novel miRNAs (PC-3p-104764 and PC-3p-129630) and their putative target genes.(Click on the image to enlarge.)
(Click on the image to enlarge.)
Identification of maize miRNAs by high-throughput sequencing
Identification of miRNAs in model plants using high-throughput sequencing or miRNA arrays has been previously reported. High-throughput sequencing has proven particularly successful for identifying plant miRNAs involved abiotic stress responses [49, 50, 51]. In this study libraries constructed from the leaves of two maize inbred lines were used to identify novel drought-associated miRNAs and their mechanisms of action were subsequently investigated.
Of the miRNAs identified using high-throughput sequencing, 85.48% of those already known were expressed at low levels (less than 10 raw reads; Table 1). This suggests that high-throughput sequencing is a powerful strategy for identifying poorly expressed miRNAs in plants. The miR159, miR167 and miR396 families, which are believed to target MYB transcription factors, auxin responsive factors and GRF transcription factors, respectively, were abundantly represented in both libraries. By comparing the expression levels of all members of an miRNA family, dominant members could be found, such as zma-MIR396d in the miR396 family, zma-MIR167h in the miR167 family and zma-MIR169a in the miR169 family. These dominant members may perform key regulatory roles in response to abiotic stress. Some family members exhibited comparable expression levels, such as zma-MIR166b/c/e/f/g/h/I in the miR166 family and zma-MIR168a and zma-MIR168b in the miR168 family. This indicates that several members of a family may have a synergistic effect in the regulation network.
Potential targets of differentially expressed miRNAs
Earlier research in Arabidopsis demonstrated that many of the potential targets of candidate miRNAs were transcription factors  and the majority of miRNAs were involved in multiple diverse biological processes . In this study, 29 differentially expressed miRNAs were identified using high-throughput sequencing, and potential targets of 23 of these were successfully predicted (Additional File 3: Table S3). No targets were identified for the remaining 6 miRNAs, and this may be due to differences in the spatial or temporal expression of a miRNA and its target that could result in insufficient degradation of the target . As anticipated, many of the potential targets were involved in transcription regulation, including MYB, CRP, PPA, HRT, GRF, SR, C2H2, SLS, SMH, RAV, GLH and Phytocyanin-type transcription factors. These transcription factor families are conserved in other plant species [53-56] and are known to regulate plant development.
Potential targets of drought-associated miRNAs
This study identified five novel miRNAs that were differentially expressed in drought-resistant and drought-sensitive maize inbred lines. GO analysis was performed to infer whether further investigation into these miRNAs would be appropriate, and many of the potential target genes play important roles in diverse biological processes, but especially in photosynthesis. Transcriptomics analysis suggested that the nine target genes were expressed only in leaves, which are the primary location for photosynthesis. The chemical constituents of leaves are known to change in response to drought stress, therefore we propose that these target genes are involved in drought stress. The results presented here will guide further functional analysis of the role of these novel drought-associated miRNAs and their target genes in the response to drought stress in maize.
Expression profiles of nine target genes. A heat map was generated by hierarchical clustering using a dedicated heat map package . Expression data were normalized and hierarchically clustered with average linkage. The color scale in the top right corner represents the relative gene expression level, where red, yellow and blue indicate high, medium and low levels of gene expression, respectively.(Click on the image to enlarge.)
Differentially expressed miRNAs were detected in the leaves of two maize inbred lines, and 124 known and 68 novel miRNAs were identified using a high-throughput sequencing approach. The sRNA-seq data has been deposited into NCBI GEO database and the GEO accession number is GSE61700. The novel drought-associated miRNAs identified and their predicted target genes may play important roles in the regulation of stress responses in this important food plant. These findings have enriched the maize miRNA databases and provide the foundation for further analysis of their roles in maize growth and development.
Supplementary MaterialAdditional File 1
Supplementary Table S1.Additional File 2
Supplementary Table S2.Additional File 3
Supplementary Tables S3-S4.
This work was supported by grants from the National Basic Research Program of China (201303001), National Key Technology Support Program (2012BAD20B00), National Natural Science Foundation of China (31371285, 31201217) and Natural Science Foundation of Anhui Province (1508085QC64).
The authors have declared that no competing interest exists.
1. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281-297
2. Jones-Rhoades MW, Bartel DP, Bartel B. MicroRNAs and their regulatory roles in plants. Ann Rev Plant Biol. 2006;57:19-53
3. Voinnet O. Origin, biogenesis, and activity of plant microRNAs. Cell. 2009;136:669-687
4. Wu L, Zhang Q, Zhou H, Ni F, Wu X, Qi Y. Rice microRNA effector complexes and targets. Plant Cell. 2009;21:3421-3435
5. Chuck G, Candela H, Hake S. Big impacts by small RNAs in plant development. Plant Biology. 2009;12:81-86
6. Sun G. MicroRNAs and their diverse functions in plants. Plant Mol Biol. 2012;80:17-36
7. Yang JH, Han SJ, Yoon EK, Lee WS. Evidence of an auxin signal pathway, microRNA167-ARF8-GH3, and its response to exogenous auxin in cultured rice cells. Nucleic Acids Research. 2006;34:1892-1899
8. Gutierrez L, Bussell JD, Pacurar DI, Schwambach J, Pacurar M, Bellini C. Phenotypic plasticity of adventitious rooting in Arabidopsis is controlled by complex regulation of AUXIN RESPONSE FACTOR transcripts and microRNA abundance. Plant Cell. 2009;21:3119-3132
9. Wang JW, Czech B, Weigel D. miR156-regulated SPL transcription factors define an endogenous flowering pathway in Arabidopsis thaliana. Cell. 2009;138:738-749
10. Chen X. A microRNA as a translational repressor of APETALA2 in Arabidopsis flower development. Science. 2004;303:2022-2025
11. Wu G, Poethig RS. Temporal regulation of shoot development in Arabidopsis thaliana by MIR156 and its target SPL3. Development. 2006;133:3539-3547
12. Wang JW, Schwab R, Czech B, Mica E, Weigel D. Dual effects of miR156-targeted SPL genes and CYP78A5/KLUH on plastochron length and organ size in Arabidopsis thaliana. Plant Cell. 2008;20:1231-1243
13. Lauter N, Kampani A, Carlson S, Goebe M, Moose SP. microRNA172 down-regulates glossy15 to promote vegetative phase change in maize. Proceedings of the National Academy of Sciences of the United States of America. 2005;102:9412-9417
14. Jones-Rhoades MW, Bartel DP. Computational identification of plant microRNAs and their targets, including a stress induced miRNA. Mol Cell. 2004;14:787-799
15. Kawashima CG, Yoshimoto N, Maruyama-Nakashita A, Tsuchiya YN, Saito K, Takahashi H, Dalmay T. Sulphur starvation induces the expression of microRNA-395 and one of its target genes but in different cell types. Plant J. 2009;57:313-321
16. Fujii H, Chiou TJ, Lin SI, Aung K, Zhu JK. A miRNA involved in phosphate-starvation response in Arabidopsis. Curr Biol. 2005;15:2038-2043
17. Chiou TJ, Aung K, Lin SI, Wu CC, Chiang SF, Su CL. Regulation of phosphate homeostasis by microRNA in Arabidopsis. Plant Cell. 2006;18:412-421
18. Chen X. MicroRNA biogenesis and function in plants. FEBS Lett. 2005;579:5923-5931
19. Qin Y, Duan Z, Xia X, Yin W. Expression profiles of precursor and mature microRNAs under dehydration and high salinity shock in Populus euphratica. Plant Cell Rep. 2011;30:1893-1907
20. Zhang B, Wang Q, Pan X. MicroRNAs and their regulatory roles in animals and plants. J Cell Physiol. 2007;210:279-289
21. Xu Z, Zhong S, Li X, Li W, Rothstein SJ, Zhang S, Bi Y. and Xie C: Genome-Wide Identification of MicroRNAs in Response to Low Nitrate Availability in Maize Leaves and Roots. PLoS ONE. 2011;6(11):e28009
22. Zhao Y, Xu Z, Mo Q, Zou C, Li W, Xu Y, Xie C. Combined small RNA and degradome sequencing reveals novel miRNAs and their targets in response to low nitrate availability in maize. Annals of Botany. 2013;112:633-642
23. Zhang Z, Lin H, Shen Y, Gao J, Xiang K, Liu L, Yuan G, Lan H, Zhou S, Gao S, Rong T, Pan G. Cloning and characterization of miRNAs from maize seedling roots under low phosphorus stress. Mol Biol Rep. 2012;39:8137-8146
24. Zhu J, Li W, Yang W, Qi L. Identification of microRNAs in Caragana intermedia by highthroughput sequencing and expression analysis of 12 microRNAs and their targets under salt stress. Plant Cell Rep. 2013
25. Ruan MB, Zhao YT, Meng ZH, Wang XJ, Yang WC. Conserved miRNA analysis in Gossypium hirsutum through small RNA sequencing. Genomics. 2009;94:263-268
26. Zhang LF, Chia JM, Kumari S, Stein JC, Liu ZJ, Narechania A, Maher CA, Guill K, McMullen MD, Ware D. A genomewide characterization of microRNA genes in Maize. PLoS Genet. 2009;5:11
27. Zeng CY, Wang WQ, Zheng Y, Chen X, Bo WP, Song S, Zhang WX, Peng M. Conservation and divergence of microRNAs and their functions in Euphorbiaceous plants. Nucleic Acids Res. 2010;38:981-995
28. Calvino M, Bruggmann R, Messing J. Characterization of the small RNA component of the transcriptome from grain and sweet sorghum stems. BMC Genomics. 2011;12:356
29. Turner M, Yu O, Subramanian S. Genome organization and characteristics of soybean microRNAs. BMC Genomics. 2012;13:169
30. Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ. miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008;36:D154-D158
31. Sunkar R, Zhu JK. Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell. 2004;16:2001-2019
32. Grad Y, Aach J, Hayes GD, Reinhart BJ, Church GM, Ruvkun G, Kim J. Computational and experimental identification of C. elegans microRNAs. Mol Cell. 2003;11:1253-1263
33. Zhang Z, Wei L, Zou X, Tao Y, Liu Z, Zheng Y. Submergence-responsive MicroRNAs are Potentially Involved in the Regulation of Morphological and Metabolic Adaptations in Maize Root Cells. Annals of Botany. 2008;102(4):509-519
34. Liu H, Qin C, Chen Z, Zuo T, Yang X, Zhou H, Xu M, Cao S, Shen Y, Lin H, He X, Zhang Y, Li L, Ding H, Lübberstedt T, Zhang Z, Pan G. Identification of miRNAs and their target genes in developing maize ears by combined small RNA and degradome sequencing. BMC Genomics. 2014;15:25
35. Zhang HY, Cui LN, Dong ST, Gao RQ, Sun AQ. Selection for Drought Resistance of 37 Maize Inbred Lines. SHANDONG AGRICULTURAL SCIENCES. 2010 (1)
36. Li M, Xia Y, Gu Y, Zhang K, Lang Q, Chen L, Guan J, Luo Z, Chen H, Li Y, Li Q, Li X, Jiang AA, Shuai S, Wang J, Zhu Q, Zhou X, Gao X. MicroRNAome of porcine pre- and postnatal development. PLoS ONE. 2010;5:e11541
37. Kozomara A, Griffiths-Jones S. miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011;39:D152-D157
38. Schnable PS, Ware D, Fulton RS, Stein JC, Wei F. et al. The B73 maize genome: complexity, diversity and dynamics. Science. 2009;326:1112-1115
39. Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, Griffiths-Jones S, Jacobsen SE, Mallory AC, Martienssen RA, Poethig RS, Qi Y, Vaucheret H, Voinnet O, Watanabe Y, Weigel D, Zhu JK. Criteria for annotation of plant MicroRNAs. Plant Cell. 2008;20:3186-3190
40. Allen E, Xie Z, Gustafson AM, Carrington JC. microRNA directed phasing during trans-acting siRNA biogenesis in plants. Cell. 2005;121:207-221
41. Schwab R, Palatnik JF, Riester M, Schommer C, Schmid M, Weigel D. Specific effects of microRNA on the plant transcriptome. Dev Cell. 2005;8:517-527
42. Dai X, Zhao PX. psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res. 2011;39:W155-W159
43. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-DDCT Method. Methods. 2001;25:402-408
44. Ambros V, Bartel B, Bartel DP, Burge CB, Carrington JC, Chen X, Dreyfuss G, Eddy SR, Griffiths-Jones S, Marshall M, Matzke M, Ruvkun G, Tuschl T. A uniform system for microRNA annotation. RNA. 2003;9(3):277-279
45. Chen X, Li QB, Wang J, Guo X, Jiang XR, Ren ZJ, Weng CY, Sun GX, Wang XQ, Liu YP, Ma LJ, Chen JY, Wang J, Zen K, Zhang JF, Zhang CY. Identification and characterization of novel am-phioxus microRNAs by Solexa sequencing. Genome Biology. 2009;10(7):R78
46. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281-297
47. Zeng QY, Yang CY, Ma QB, Li XP, Dong WW, Nian H. Identification of wild soybean miRNAs and their target genes responsive to aluminum stress. BMC Plant Biology. 2012;12:182
48. Sekhon RS, Lin H, Childs KL, Hansey CN, Robin Buell C, De Leon N, Kaeppler SM. Genome-wide atlas of transcription during maize development. Plant J. 2011;66:553-563
49. Navarro L, Dunoyer P, Jay F, Arnold B, Dharmasiri N, Estelle M, Voinnet O, Jones JDG. A plant miRNA contributes to antibacterial resistance by repressing auxin signaling. Science. 2006;312:436-439
50. Kurihara Y, Kaminuma E, Matsui A, Kawashima M, Tanaka M, Morosawa T, Ishida J, Mochizuki Y, Shinozaki K, Toyoda T, Seki M. Transcriptome Analyses Revealed Diverse Expression Changes in ago1 and hyl1 Arabidopsis Mutants. Plant and Cell Phydiology. 2009;50:1715-1720
51. Chen L, Wang T, Zhao M, Tian Q, Zhang WH. Identification of aluminumresponsive microRNAs in Medicago truncatula by genome-wide highthroughput sequencing. Planta. 2012;235(2):375-386
52. Sunkar R, Zhu JK. Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell. 2004;16:2001-2019
53. Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP. Prediction of plant microRNA targets. Cell. 2002;110(4):513-520
54. Jones-Rhoades MW, Bartel DP. Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004;14(6):787-799
55. Sunkar R, Girke T, Jain PK, Zhu JK. Cloning and characterization of microRNAs from rice. Plant Cell. 2005;17(5):1397-1411
56. Buhtz A, Springer F, Chappell L, Baulcombe DC, Kehr J. Identification and characterization of small RNAs from the phloem of Brassica napus. Plant J. 2008;53(5):739-749
57. Sekhon RS, Briskine R, Hirsch CN, Myers CL, Springer NM, Buell CR, de Leon N, Kaeppler SM. Maize gene atlas developed by RNA sequencing and comparative evaluation of transcriptomes based on RNA sequencing and microarrays. PLoS One. 2013:8 (4)
Corresponding authors: Email: hyjiangedu.cn or beijiuchengedu.cn