- Research article
- Open Access
Characterization of cadmium-responsive MicroRNAs and their target genes in maize (Zea mays) roots
BMC Molecular Biologyvolume 20, Article number: 14 (2019)
Current research has shown that microRNAs (miRNAs) play vital roles in plant response to stress caused by heavy metals such as aluminum, arsenic, cadmium (Cd), and mercury. Cd has become one of the most hazardous pollutants in the environment. Maize can be a potential model to study phytoremediation of Cd-contaminated soil owing to its large biomass production. However, little is known about miRNAs as a response to Cd stress in maize.
To investigate the role of miRNAs in response to Cd stress, roots of seedlings of the inbred maize lines B73 and Mo17 were collected and treated with 200 mg/L CdCl2·2.5 H2O over different exposure times. Enzyme activities of superoxide dismutase and peroxidase were measured to confirm Cd stress. The expression of six candidate miRNAs and their targets were validated using quantitative real-time PCR (qRT-PCR) technology. In addition, the expression of Zma-miR171b was assessed using in situ hybridization.
Our results showed that miRNAs and their respective target genes were differentially expressed in maize seedling roots exposed to Cd stress. This research produced new insights into the molecular mechanism of miRNAs responsive to Cd stress in plants and sheds light on the latent roles of miRNAs in plants exposed to heavy metal stresses.
Contamination of soils with heavy metals has become a prevalent topic in ecology owing to their widespread release from industry, agriculture, and other human activities . Cadmium (Cd) adversely affects various physiological and biochemical processes such as nitrogen metabolism, photosynthesis, mineral utilization, carbohydrate metabolism, and water uptake, thereby inhibiting plant growth and development [2, 3]. Cd pollution has substantially increased over the least four decades worldwide, which is a matter of concern for human health . Therefore, research on the molecular mechanisms underlying Cd pollution in plants is urgently required.
Recently, several studies investigated the molecular mechanisms association with Cd tolerance and toxicity in plants [2, 3]. Exceptional microRNAs (miRNA) expression was identified in tissues and cells of stressed plants, and the involvement of miRNAs in Cd-stress responses was suggested [5,6,7]. In rapeseed under Cd stress, the miRNA bna-miR393 was specifically expressed in leaves, bna-miR156a and bna-miR167a/c were expressed in roots and leaves, respectively, and bna-miR164b and bna-miR394a/b/c were expressed in all tissues . Similarly, differential expression of numerous non-conserved and conserved miRNAs was reported to be induced by Cd stress in Brassica napus roots . In the same study, 84 conserved and non-conserved miRNAs families with differential expression were identified in shoots and roots under Cd stress . In rice, 19 Cd-responsive miRNAs and their targets including transcription factors and stress-related proteins were identified in seedlings using a microarray assay . In addition, to assess Cd-responsive miRNAs and their associated genes in radish (Raphanus sativus), small-RNA libraries were constructed from plants under Cd stress and from controls. A total of 23 miRNA families with notable differential expression (15 previously known and 8 novel miRNA families) was identified in Cd-stressed plants, and their targets were then predicted using a degradome analysis . Numerous potential heavy Cd-responsive miRNA candidates were identified in rice, B. napus, radish, and many other plants [5,6,7]. However, information on the molecular mechanism of Cd stress resistance in maize is limited. Therefore, we performed a characterization of Cd stress-mediated miRNAs and their target genes in roots of B73 and Mo17 maize under Cd stress and in negative controls. We identified target genes associated with Cd-regulated candidate miRNAs in maize seedling roots. The candidate miRNAs that respond to Cd stress will be useful for further improving Cd tolerance in maize.
Expression pattern of conserved miRNAs in plant under different heavy metal stress
Recent studies showed that miRNAs play vital roles in regulating the expression of stress response genes following stress caused by various heavy metals. In this study, we assessed the expression patterns of miRNAs in plant exposed to different heavy metals, including aluminum (Al), arsenic (As), cadmium (Cd), and mercury (Hg) [6, 8,9,10,11]. Significant differential expression of miRNAs due to heavy metal stress was observed, such as down-regulation of common miRNAs (miR159, miR160, miR166, miR319, and miR396) and up-regulation of miRNA miR393. Other miRNAs also showed differential expression in response to heavy metal stress: miR156, miR159, miR160, miR166, miR167, miR171, miR319, miR390, miR394, miR396, and miR408 were down-regulated, and miR393 was up-regulated under Cd stress (Fig. 1). A conservation analysis of those miRNAs showed that miR156/miR157, miR159/319, miR166, and miR170/miR171 were higher conserved in all organisms, followed by miR160, miR396, miR167, and miR390. The miRNAs miR398 and miR394 were less conserved in all organisms (Additional file 1: Table S4). Based on the results of significant differential expression and conservation analyses of these miRNAs under Cd stress, miR156, miR166, miR167, miR171, and miR393 were selected as candidates for differentially expressed miRNAs responsive to Cd stress. In addition, we analyzed the expression profiles of these miRNAs and their target genes in the inbred maize lines B73 and Mo17 using qRT-PCR.
Alterations in the development of maize leaves and root system in response to Cd stress
To confirm the appropriate treatment time for validation of candidate Cd-regulated miRNAs, we analyzed the enzymatic activities of superoxide dismutase (SOD) and peroxidase (POD) in Cd-treated maize roots and leaves after treatment times of 6–96 h. We found that SOD activity increased after 12, 24, and 48 h of the Cd-treatment, and particularly so after 12 h, in both Mo17 and B73. POD activity showed the opposite pattern in both phenotypes. After 48 h, SOD and POD activity gradually decreased to a constant level. SOD and POD activity was higher in roots than in leaves, in both phenotypes (Additional file 2: Figure S1). These results suggest that certain Cd-tolerant mechanisms were activated to maintain the normal physiological and metabolic activities shortly after the Cd treatment. Subsequently, RNA was isolated from the B73 and Mo17 treatments and control groups.
Cd-regulated miRNA expression levels
To assess the functions of these candidate miRNAs selected according to altered expression levels of conserved miRNAs in plants under Cd stress, we analyzed their expression profiles using qRT-PCR. Zma-miR156b, Zma-miR156k, Zma-miR166d, Zma-miR167f, Zma-miR171b, and Zma-miR393b were expressed in maize roots of both maize lines (Fig. 2). Zma-miR156b showed a trend of down-regulation in roots after the Cd treatment in B73, compared with the control, whereas no significant changes were observed in Mo17 (Fig. 2a). Zma-miR156k showed no significant alterations in B73 and Mo17 (Fig. 2b). Zma-miR166d was more down-regulated in B73 than in Mo17 roots under Cd stress (Fig. 2c), and Zma-miR167f was more down-regulated in Mo17 than in B73 (Fig. 2d). Zma-miR171b was down-regulated in both maize lines (Fig. 2e), in contrast to miR393b (Fig. 2f).
qRT-PCR analysis of target genes related to candidate microRNAs in response to Cd stress
The results suggested that the activities of miRNA-mRNA complexes could be altered by stress. To identify target genes of candidate miRNAs regulated by Cd stress, the target genes were predicted using WMD3 and PsRNA target website using the MGC database (ZmB73_v5b, Table. S5). Most of the target genes association with the candidate miRNAs belonged to transcription factor families, e.g., Zma-miRNA156 targeting the SBP transcription factor family, Zma-miR166 acting on the basic-leucine zipper (bZIP) transcription factors, Zma-miR167 regulating auxin response factors (ARFs), GRAS transcription factors were predicted targets of Zma-miR171, and F-box proteins were predicted targets of Zma-miR393. The relative expression level of target genes seemed negatively correlated with that of the respective miRNAs (Fig. 3). The target genes of Zma-miR156b (SBP-box), Zma-miR166d (bZIP), and Zma-miR171b (GRAS) were significantly up-regulated during the dynamic development of roots in B73, compared to in Mo17 (Fig. 3a, c, e). In contrast, the target gene of Zma-miR393b (F-box-like protein) was observed to be down-regulated in B73 than in Mo17 samples (Fig. 3f). However, the target genes of Zma-miR156k (SBP-box) and Zma-miR167f (ARF) showed no significant up-regulation in either maize line (Fig. 3b, d). In conclusion, candidate miRNAs such as Zma-miR156b, Zma-miR166d, Zma-miR171b, and Zma-miR393b may negatively regulate their respective target genes in maize under Cd stress, which likely involves transcriptional and post-transcriptional regulatory networks.
Expression of Zma-miR171b identified by in situ hybridization
Based on the results obtained from the qRT-PCR of candidate miRNAs and their respective targets genes, the expression of Zma-miR171b was selected and further analyzed using in situ hybridization (ISH). Black-blue positive hybrid signals of Zma-miR171b were detected throughout at 0, 4, 12, and 24 h of treatment in maize roots of B73 and Mo17. Expression of Zma-miR171b was down-regulated following Cd stress at different stages in B73 roots, compared to that in the control, and reached a minimum at 48 h; similar trends were observed in Mo17, but expression levels of Zma-miR171b were consistently higher in B73 than in Mo17 (Fig. 4). Zma-miR171b was down-regulated in both phenotypes, thus B73 likely accumulated GRAS gene products earlier in response to Cd stress than Mo17 (Figs. 2e, 3e).
Gene ontology annotation and analysis of target genes
Gene ontology (GO) analyses were performed to gain understanding of the function of target genes. Target genes that were specifically enriched in the regulation of biological processes showed transcription regulation activity and binding activity; however, those target genes were decreased in metabolism processes, and their catalytic activities were also reduced (Additional file 2: Figure S2).
miRNAs are one of the major types of endogenous, small, non-protein-coding single-stranded RNA, which act as regulators of plant growth and development at the translation and post-transcription levels by silencing gene expression [12, 13]. Previous studies demonstrated that miRNAs acts as regulators in plants exposed to abiotic stressors . Moreover, many miRNAs were found to be important for plant responses to heavy metal stress using small RNA profiling [15, 16]. Recent studies identified numerous heavy metals such as Hg, Al, As, and Cd to regulate miRNAs in plants [6, 8,9,10,11]. However, we found that heavy metals regulated miRNAs varied from different plants and metal species. As shown in Fig. 1, some miRNA families (miR396, miR390, miR171, miR166, miR162, and miR159) were down-regulated, whereas miR395 and miR393 were up-regulated by most heavy metals. The expression of miR396, miR390, and miR159 was consistently down-regulated, and miR393 was up-regulated under exposure to four heavy metals. Ding et al. demonstrated that miR390, miR171, miR168, miR166, and miR156 responded to Cd stress in rice , and Zhou et al. found that expression of miR408, miR398, miR397, and miR396 was associated with Cd exposure in Brassica . In addition, miR393 and miR171 were found to be affected by Cd stress through negatively regulating their target genes in Oryza sativa , Medicago truncatula , and B. napus . In the current study, we demonstrated that Zma-miR156b, Zma-miR166d, Zma-miR171b, and Zma-miR393b were affected by Cd exposure in maize roots, using qRT-PCR-based analyses.
Recently, identification of target genes of metal-regulated miRNAs showed that most of these miRNAs encoded key regulatory proteins or TFs that were associated with diverse metabolic pathways [20,21,22]. In our study, several candidate miRNAs such as Zma-miR156b, Zma-miR166d, Zma-miR167f, Zma-miR171b, and Zma-miR393b were found to target the genes of TF families (SBP-box, bZIP, ARFs, GRAS), and F-box-like protein, most of which have been confirmed to play vital roles in the response to Cd stress (Additional file 1: Table S5, Fig. 4). Previous studies confirmed that SBP-box transcription factors were involved in a broad range of developmental and stress response processes targeted by miR156 in Arabidopsis [23, 24]. In addition, the SBP-box gene family was characterized as metal-containing transcription factors regulating copper homeostasis in Arabidopsis . In the present study, two miRNAs belonged to the miRNA156 family (Zma-miR156b and Zma-miR156k) targeting SBP-box transcription factors. It may be inferred that miR156 may be an important regulator in Cd2+ homeostasis in maize roots by targeting SBP-box transcription factors. Moreover, for the bZIP family, bZIP62, ThbZIP1, and BjCdR15 were reported to be associated with Cd stress. Of those, the gene encoding bZIP62 was induced by Cd exposure in soybean roots, whereas ThbZIP1 showed increased expression in Tamarix hispida [26, 27]. In the current study, we found that bZIP transcription factors were regulated by Zma-miRNA166d. We suggest that Zma-miR166 acted in Cd2+ homeostasis in maize roots by targeting bZIP transcription factors. Furthermore, previous studies showed that some miRNAs responded to heavy metal stress through regulating phytohormone biosynthesis and signaling responses. For instance, As stress in B. juncea led to an increased miR167 abundance . The auxin response factors ARF6 and ARF8 targeted by miR167 were found positively correlated with adventitious rooting and development of anthers and ovules via auxin homeostasis in Arabidopsis . Moreover, ARF proteins may activate or repress the transcription of primary auxin response genes through binding to auxin-responsive elements . F-box proteins, including transport inhibitor response 1 (TIR1) and auxin-signaling F-box proteins (AFB1, AFB2, and AFB3) can perceive auxins in plants [30,31,32,33]. In addition, TIR1 encoded an E3 ubiquitin ligase containing F-box subunit were targeted by miR393  which interacts with Aux/IAA proteins as an auxin receptor, thereby resulting in Aux/IAA ubiquitination and subsequent degradation in auxin signal pathway . In our study, we demonstrated that Zma-miR393b and Zma-miR167f targeted auxin signaling F-box proteins and that ARFs transcription factors responded to the Cd2+ treatment.
In conclusion, this is the first study to characterize Cd-stress induced miRNAs and their target genes in maize roots. Six previously known miRNAs and their respective target genes were found to respond to Cd stress. The expression patterns of these differentially regulated miRNAs and association targets were shown to be affected by Cd stress in the inbred maize lines B73 and Mo17. Most of the target transcripts belonged to transcription factors that were predicted to be functionally involved in different pathways that respond to Cd stress. These results provide new insights for the functional characterization of miRNAs and their targets in plant responses to heavy metal stress.
Plant material and experimental design
The inbred maize line B73 was previously reported to be more tolerant to heavy metal stress than Mo17 , particularly regarding Cd stress. In this study, seeds collected from Mo17 and B73 were washed using distilled water and germinated at 26 °C in the dark for 3 days. Germinated seeds were grown in a growth chamber under the following conditions: a photoperiod of 14 h light/10 h dark and relative humidity of 70% at 26 °C. After 20 days, the seedlings were transplanted to a plastic container with a modified half-strength Hoagland nutrient medium (see Additional file 1: Table S1). Subsequently, the plants were randomly assigned to a Cd stress treatment (200 mg/L CdCl2·2.5 H2O) and a control (Cd-free).
Enzymatic activity assays of POD and SOD
Maize seedlings of the treatment group and the control were collected after exposed to 200 mg/L CdCl2·2.5H2O for 0, 4, 8, 12, 24, 48, 72, or 96 h. POD and SOD activities were determined according to the methods described by Beyer  and Kim , respectively. All measurements were performed three times using three biological replicates.
RNA isolation from maize roots
Based on the results of the POD and SOD enzymatic activity assays, total RNA was isolated from Cd-free and Cd-exposed roots of B73 and Mo17 after 0, 4, 8, 12, 24, and 48 h, using Trizol reagent (Life, US) following the manufacturer’s instructions. Small RNA was isolated from each sample using the mirVana™miRNA Isolation Kit (Ambion, US) according to the manufacturer’s instructions.
Validation of candidate miRNAs using real-time PCR
Total miRNA was isolated from B73 and Mo17 using a Universal Plant microRNA Extraction Kit (Spin-column, BioTeke, China). Reverse transcription reactions of miRNAs were performed, and PCRs were performed in triplicate, according to Jian et al. . Mitochondrial 5S RNA was used as a control, and specific primers were for real-time PCR (Additional file 1: Table S2). The absolute amount of each miRNA was calculated according to a standard curve using the 2−∆∆CT method .
ISH assays for Zma-miR171b
For ISH, miRCURY 5′-DIG and 3′-DIG-labeled locked nucleic acid (LNA)-miRNAs were utilized (Zma-miR171b) as detection probes. The sequence of the probe Zma-miR171b was 5′-DIG-UUGAGCCGUGCCAAUAUCAC-DIG-3′. Sections of 30 µm thickness were performed using LNA probes, and all steps were carried out as described by Urbanek et al. .
Target gene prediction of candidate miRNAs in the Cd2+ treatment
To predict the accuracy of target genes, the online tools WMD3 (http://wmd3.weigelworld.org/cgi-bin/webapp.cgi?) and psRNA (http://plantgrn.noble.org/psRNATarget/) were used with default settings [42, 43].
qRT-PCR confirmation of target genes
Specific qRT-PCR primers were designed for target genes using the software Primer Premier (version 5.0); the primers are shown in Additional file 1: Table S3. The amplification programs were performed according to the standard protocol of the ABI7500 system and were conducted in triplicate, according to the method of Jian et al. . The cycle threshold (Ct) of each tested gene were averaged for triplicate reactions, and the values were normalized according to the Ct of the control products of the Actin1 gene. Statistical analysis was performed using the 2−∆∆CT method.
GO annotation of putative target genes responsive to Cd stress
To annotate target genes we used singular enrichment analysis by examining the AgriGO database . Using a customized or pre-calculated background, enriched GO terms were removed after statistical testing. GO annotation comprised biological processes, cellular components, and component functions of the putative sequences.
quantitative real-time PCR
in situ hybridization
singular enrichment analysis
auxin response factor
Valko M, Morris H, Cronin M. toxicity and oxidative stress. Curr Med Chem. 2005;12:1161–208.
Besson-Bard A, Gravot A, Richaud P, et al. Nitric oxide contributes to cadmium toxicity in Arabidopsis by promoting cadmium accumulation in roots and by up-regulating genes related to iron uptake. Plant Physiol. 2009;149:1302–15.
DalCorso G, Farinati S, Furini A. Regulatory networks of cadmium stress in plants. Plant Signal Behav. 2010;5:663–7.
Ishikawa S, Noriharu A, Murakami M, Wagatsuma T. Is Brassica juncea a suitable plant for phytoremediation of cadmium in soils with moderately low cadmium contamination?–Possibility of using other plant species for Cd-phytoextraction. Soil Sci Plant Nutr. 2006;52:32–42.
Zhao SZ, Jian BS, Zhi MY. Genome-wide identification of Brassica napus microRNAs and their targets in response to cadmium. J Exp Bot. 2012;63:4597–613.
Ding Y, Qu A, Gong S, et al. Molecular identification and analysis of Cd-responsive microRNAs in rice. J Agric Food Chem. 2013;61:11668–75.
Xu L, Wang Y, Zhai L, et al. Genome-wide identification and characterization of cadmium-responsive microRNAs and their target genes in radish (Raphanus sativus L.) roots. J Exp Bot. 2013;64:4271–87.
Zhou ZS, Zeng HQ, Liu ZP, Yang ZM. Genome-wide identification of Medicago truncatula microRNAs and their targets reveals their differential regulation by heavy metal. Plant Cell Environ. 2012;35:86–99.
Valdéslópez O, Yang SS, Apariciofabre R, et al. MicroRNA expression profile in common bean (Phaseolus vulgaris) under nutrient deficiency stresses and manganese toxicity. New Phytol. 2010;187:805–18.
Chen L, Wang T, Zhao M, et al. Identification of aluminum-responsive microRNAs in Medicago truncatula by genome-wide high-throughput sequencing. Planta. 2012;235:375–86.
Srivastava S, Srivastava AK, Suprasanna P, D’Souza SF. Identification and profiling of arsenic stress-induced microRNAs in Brassica juncea. J Exp Bot. 2013;64:303–15.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.
Mallory AC, Vaucheret H. Functions of microRNAs and related small RNAs in plants. Nat Genet. 2006;38:S31–6.
Sunkar R, Kapoor A, Zhu J-K. Posttranscriptional induction of two Cu/Zn superoxide dismutase genes in Arabidopsis is mediated by down regulation of miR398 and important for oxidative stress tolerance. Plant Cell. 2006;18:2051–65.
Ding Y-F, Zhu C. The role of microRNAs in copper and cadmium homeostasis. Biochem Biophys Res Commun. 2009;386:6–10.
Ding Y, Chen Z, Zhu C. Microarray-based analysis of cadmium-responsive microRNAs in rice (Oryza sativa). J Exp Bot. 2011;62:3563–73.
Huang Y, Ying H, Yunxia L. Combined toxicity of copper and cadmium to six rice genotypes (Oryza sativa L.). J Environ Sci. 2009;21:647–53.
Huang SQ, Xiang AL, Che LL, et al. A set of miRNAs from Brassica napus in response to sulphate deficiency and cadmium stress. Plant Biotechnol J. 2010;8:887–99.
Van Ex F, Jacob Y, Martienssen RA. Multiple roles for small RNAs during plant reproduction. Curr Opin Plant Biol. 2011;14:588–93.
Noman A, Fahad S, Aqeel M, et al. miRNAs: major modulators for crop growth and development under abiotic stresses. Biotechol Lett. 2017;39(5):685–700.
Noman A, Aqeel M. miRNA-based heavy metal homeostasis and plant growth. Environ Sci Pollut Res. 2017;24(11):10068–82.
Yu S, Galvão VC, Zhang Y-C, et al. Gibberellin regulates the Arabidopsis floral transition through miR156-targeted SQUAMOSA PROMOTER BINDING-LIKE transcription factors. Plant Cell. 2012;24:3320–32.
Pilon M, Cohu CM, Ravet K, et al. Essential transition metal homeostasis in plants. Curr Opin Plant Biol. 2009;12:347–57.
Kropat J, Tottey S, Birkenbihl RP, et al. A regulator of nutritional copper signaling in Chlamydomonas is an SBP domain protein that recognizes the GTAC core of copper response element. PNAS. 2005;102:18730–5.
Chmielowska-Bąk J, Lefèvre I, Lutts S, Deckert J. Short term signaling responses in roots of young soybean seedlings exposed to cadmium stress. J Plant Physiol. 2013;170:1585–94.
Xu J, Wang W, Yin H, et al. Exogenous nitric oxide improves antioxidative capacity and reduces auxin degradation in roots of Medicago truncatula seedlings under cadmium stress. Plant Soil. 2010;326:321–30.
Gutierrez L, Mongelard G, Floková K, et al. Auxin controls Arabidopsis adventitious root initiation by regulating jasmonic acid homeostasis. Plant Cell. 2012;24:2515–27.
Hagen G, Guilfoyle T. Auxin-responsive gene expression: genes, promoters and regulatory factors. Plant Mol Biol. 2002;49:373–85.
Navarro L, Dunoyer P, Jay F, et al. A plant miRNA contributes to antibacterial resistance by repressing auxin signaling. Sci Signal. 2006;312:436.
Dharmasiri N, Dharmasiri S, Estelle M. The F-box protein TIR1 is an auxin receptor. Nature. 2005;435:441–5.
Dharmasiri N, Dharmasiri S, Weijers D, et al. Plant development is regulated by a family of auxin receptor F box proteins. Dev Cell. 2005;9:109–19.
Kepinski S, Leyser O. The Arabidopsis F-box protein TIR1 is an auxin receptor. Nature. 2005;435:446–51.
Chen Z-H, Bao M-L, Sun Y-Z, et al. Regulation of auxin response by miR393-targeted transport inhibitor response protein 1 is involved in normal development in Arabidopsis. Plant Mol Biol. 2011;77:619–29.
Si-Ammour A, Windels D, Arn-Bouldoires E, et al. miR393 and secondary siRNAs regulate expression of the TIR1/AFB2 auxin receptor clade and auxin-related development of Arabidopsis leaves. Plant Physiol. 2011;157:683–91.
Baxter IR, Gustin JL, Settles AM, Hoekenga OA. Ionomic Characterization of Maize Kernels in the Intermated B73 × Mo17 Population. Crop Sci. 2013;53:208–20.
Beyer WF, Fridovich I. Assaying for superoxide dismutase activity: some large consequences of minor changes in conditions. Anal Biochem. 1987;161:559–66.
Kim YH, Yoo YJ. Peroxidase production from carrot hairy root cell culture. Enzyme Microbial Technol. 1996;18:531–5.
Gao J, Luo M, Zhang C, et al. A putative pathogen-resistant regulatory pathway between MicroRNAs and candidate target genes in maize. J Plant Biol. 2015;58(4):211–9.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3:1101–8.
Urbanek MO, Nawrocka AU, Krzyzosiak WJ. Small RNA detection by in situ hybridization methods. Int J Mol Sci. 2015;16:13259–86.
Dai X, Zhao PX. psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res. 2011;39(suppl 2):W155–9.
Ossowski S, Schwab R, Weigel D. Gene silencing in plants using artificial microRNAs and other small RNAs. Plant J. 2008;53:674–90.
Du Z, Zhou X, Ling Y, et al. agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010;38(suppl 2):W64–70.
JG and ML designed research and they contributed equally to this study; HP and WL performed research. WL and JG analyzed data and conducted qRT-PCR experiments; FC and JG wrote the paper. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analysed during this study are included in this published article [and its additional information files].
Consent for publication
Ethics approval and consent to participate
Normal Program supported this research for Science and Technology Research Program of Sichuan Municipal Education Commission (2018JY0170), Special Project of Science and Technology Innovation for Social Undertaking and People’s Livelihood Guarantee of CSTC (No. cstc2017shms-xdny80074).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.