Selection of Reference Genes for Normalization of MicroRNA Expression by RT-qPCR in Sugarcane Buds under Cold Stress

Screening and validation of optimal miRNA reference genes in different developing stages and tissues of Lilium henryi Baker

Dynamic miRNA detection using the qRT-PCR technique requires appropriate reference genes to ensure data reliability. Previous studies have screened internal reference genes in plants during embryonic development and various stress treatment, involving relatively few tissues and organs. There is no relevant miRNA study in Lilium henryi Baker and limited research on the optimal miRNA reference genes in lilies, such as 5S, 18S, U6 and Actin. Twelve genes were selected as candidate reference genes whose expression stability was analyzed in petals at different developmental stages and other tissues using various algorithms, such as geNorm, NormFinder, BestKeeper, and Delta CT. The results revealed that the optimal combination of reference genes for Lilium henryi Baker petals at different developmental stages was osa-miR166m and osa-miR166a-3p, while that for different tissues of Lilium henryi Baker was osa-miR166g-3p and osa-miR166a-3p.Four important genes related to growth and development regulation, namely, osa-miR156a, osa-miR395b, osa-miR396a-3p, and osa-miR396a-5p, were selected for validation. The findings of the present study could contribute to future investigations onmiRNA expression and the related functions in Lilium henryi Baker while providing important references for the normalization of the miRNA expression in other varieties of lily.

Integration of mRNA and microRNA analysis reveals the molecular mechanisms underlying drought stress tolerance in maize (Zea mays L.)

Drought is among the most serious environmental issue globally, and seriously affects the development, growth, and yield of crops. Maize (Zea mays L.), an important crop and industrial raw material, is planted on a large scale worldwide and drought can lead to large-scale reductions in maize corn production; however, few studies have focused on the maize root system mechanisms underlying drought resistance. In this study, miRNA-mRNA analysis was performed to deeply analyze the molecular mechanisms involved in drought response in the maize root system under drought stress. Furthermore, preliminary investigation of the biological function of miR408a in the maize root system was also conducted. The morphological, physiological, and transcriptomic changes in the maize variety "M8186" at the seedling stage under 12% PEG 6000 drought treatment (0, 7, and 24 h) were analyzed. With prolonged drought stress, seedlings gradually withered, the root system grew significantly, and abscisic acid, brassinolide, lignin, glutathione, and trehalose content in the root system gradually increased. Furthermore, peroxidase activity increased, while gibberellic acid and jasmonic acid gradually decreased. Moreover, 32 differentially expressed miRNAs (DEMIRs), namely, 25 known miRNAs and 7 new miRNAs, and 3,765 differentially expressed mRNAs (DEMRs), were identified in maize root under drought stress by miRNA-seq and mRNA-seq analysis, respectively. Through combined miRNA-mRNA analysis, 16 miRNA-target gene pairs, comprising 9 DEMIRs and 15 DEMRs, were obtained. In addition, four metabolic pathways, namely, "plant hormone signal transduction", "phenylpropane biosynthesis", "glutathione metabolism", and "starch and sucrose metabolism", were predicted to have important roles in the response of the maize root system to drought. MiRNA and mRNA expression results were verified by real-time quantitative PCR. Finally, miR408a was selected for functional analysis and demonstrated to be a negative regulator of drought response, mainly through regulation of reactive oxygen species accumulation in the maize root system. This study helps to elaborate the regulatory response mechanisms of the maize root system under drought stress and predicts the biological functions of candidate miRNAs and mRNAs, providing strategies for subsequent mining for, and biological breeding to select for, drought-responsive genes in the maize root system.

Identification and evaluation of the novel genes for transcript normalization during female gametophyte development in sugarcane

Background

Sugarcane (Saccharum spontaneum L.), the major sugar and biofuel feedstock crop, is cultivated mainly by vegetative propagation worldwide due to the infertility of female reproductive organs resulting in the reduction of quality and output of sugar. Deciphering the gene expression profile during ovule development will improve our understanding of the complications underlying sexual reproduction in sugarcane. Optimal reference genes are essential for elucidating the expression pattern of a given gene by quantitative real-time PCR (qRT-PCR).

Method

In this study, based on transcriptome data obtained from sugarcane ovule, eighteen candidate reference genes were identified, cloned, and their expression levels were evaluated across five developmental stages ovule (AC, MMC, Meiosis, Mitosis, and Mature).

Results

Our results indicated that FAB2 and MOR1 were the most stably expressed genes during sugarcane female gametophyte development. Moreover, two genes, cell cycle-related genes REC8 and CDK, were selected, and their feasibility was validated. This study provides important insights into the female gametophyte development of sugarcane and reports novel reference genes for gene expression research on sugarcane sexual reproduction.

Selection and Optimization of Reference Genes for MicroRNA Expression Normalization by qRT-PCR in Chinese Cedar (Cryptomeria fortunei) under Multiple Stresses

MicroRNA (miRNA) expression analysis is very important for investigating its functions. To date, no research on reference genes (RGs) for miRNAs in gymnosperms, including Cryptomeria fortunei, has been reported. Here, ten miRNAs (i.e., pab-miR159a, cln-miR162, cas-miR166d, pab-miR395b, ppt-miR894, cln-miR6725, novel1, novel6, novel14 and novel16) and three common RGs (U6, 5S and 18S) were selected as candidate RGs. qRT-PCR was used to analyse their expressions in C. fortunei under various experimental conditions, including multiple stresses (cold, heat, drought, salt, abscisic acid and gibberellin) and in various tissues (roots, stems, tender needles, cones and seeds). Four algorithms (delta Ct, geNorm, NormFinder and BestKeeper) were employed to assess the stability of candidate RG expression; the geometric mean and RefFinder program were used to comprehensively evaluate RG stability. According to the results, novel16, cln-miR6725, novel1 and U6 were the most stable RGs for studying C. fortunei miRNA expression. In addition, the expression of three target miRNAs (aly-miR164c-5p, aly-miR168a-5p and smo-miR396) was examined to verify that the selected RGs are suitable for miRNA expression normalisation. This study may aid further investigations of miRNA expression/function in the response of C. fortunei to abiotic stress and provides an important basis for the standardisation of miRNA expression in other gymnosperm species.

Identification and validation of miRNA reference genes in poplar under pathogen stress

Quantitative real time polymerase chain reaction (qRT-PCR) is a common method to analyze gene expression. Due to differences in RNA quantity, quality, and reverse transcription efficiency between qRT-PCR samples, reference genes are used as internal standards to normalize gene expression. However, few universal genes, especially miRNAs, have been identified as reference so far. Therefore, it is essential to identify reference genes that can be used across various experimental conditions, stress treatments, or tissues. In this study, 14 microRNAs (miRNAs) and 5.8S rRNA were assessed for expression stability in poplar trees infected with canker pathogen. Using geNorm, NormFinder and Bestkeeper reference gene analysis programs, we found that miR156g and miR156a exhibited stable expression throughout the infection process. miR156g, miR156a and 5.8S rRNA were then tested as internal standards to measure the expression of miR1447 and miR171c, and the results were compared to small RNA sequencing (RNA-seq) data. We found that when miR156a and 5.8S rRNA were used as the reference gene, the expression of miR1447 and miR171c were consistent with the small RNA-seq expression profiles. Therefore, miR156a was the most stable miRNAs examined in this study, and could be used as a reference gene in poplar under canker pathogen stress, which should enable comprehensive comparisons of miRNAs expression and avoid the bias caused by different length between detected miRNAs and traditional reference genes. The present study has expanded the miRNA reference genes available for gene expression studies in trees under biotic stress.

Selection of reference genes for normalization of microRNA expression in sugarcane stalks during its interaction with Colletotrichum falcatum

The microRNAs role in various cellular and metabolic functions is gaining more limelight in line with second-generation NGS technology. For the validation of candidate miRNA genes, the quantitative real-time PCR is the widely trusted and efficient method to follow. Sugarcane miRNAs are less explored in sugarcane defense response during their interaction with Colletotrichum falcatum inciting red rot. Further, for RT-qPCR experiments involving sugarcane miRNA expression studies, a stable internal reference gene is required. Hence, we have taken a study involving 20 candidate genes to identify stable expressing reference genes using NormFinder, geNorm, BestKeeper, and deltaCt statistical algorithms. The candidate reference genes included miRNAs and protein-coding genes. The results indicated that there is a variation in ranking among the algorithms. We found miR1862c as the stably expressed miRNA reference gene among the candidates and miR444b.2 along miR1862c formed the best reference gene pair combination, which can be used in the experiments aiming to explore sugarcane miRNAs in the defense mechanism against C. falcatum. The stable miRNA reference gene was further validated with other lesser stable reference gene candidates to assess the effect of stable reference genes during normalization. The present study evaluating the sugarcane miRNAs as reference genes for normalizing RT-qPCR expression data involving miRNAs during sugarcane × C. falcatum interaction is the first of its kind. Further, this systematic approach can be followed to assess the reference gene in various experimental conditions involving sugarcane miRNAs.

Enhanced sugar accumulation and regulated plant hormone signalling genes contribute to cold tolerance in hypoploid Saccharum spontaneum

BACKGROUND

Wild sugarcane Saccharum spontaneum plants vary in ploidy, which complicates the utilization of its germplasm in sugarcane breeding. Investigations on cold tolerance in relation to different ploidies in S. spontaneum may promote the exploitation of its germplasm and accelerate the improvement of sugarcane varieties.

RESULTS

A hypoploid clone 12-23 (2n = 54) and hyperploid clone 15-28 (2n = 92) of S. spontaneum were analysed under cold stress from morphological, physiological, and transcriptomic perspectives. Compared with clone 15-28, clone 12-23 plants had lower plant height, leaf length, internode length, stem diameter, and leaf width; depressed stomata and prominent bristles and papillae; and thick leaves with higher bulliform cell groups and thicker adaxial epidermis. Compared with clone 15-28, clone 12-23 showed significantly lower electrical conductivity, significantly higher water content, soluble protein content, and superoxide dismutase activity, and significantly higher soluble sugar content and peroxidase activity. Under cold stress, the number of upregulated genes and downregulated genes of clone 12-23 was higher than clone 15-28, and many stress response genes and pathways were affected and enriched to varying degrees, particularly sugar and starch metabolic pathways and plant hormone signalling pathways. Under cold stress, the activity of 6-phosphate glucose trehalose synthase, trehalose phosphate phosphatase, and brassinosteroid-signalling kinase and the content of trehalose and brassinosteroids of clone 12-23 increased.

CONCLUSIONS

Compared with hyperploid clone 15-28, hypoploid clone 12-23 maintained a more robust osmotic adjustment system through sugar accumulation and hormonal regulation, which resulted in stronger cold tolerance.

Identification of Appropriate Reference Genes for Normalizing miRNA Expression in Citrus Infected by Xanthomonas citri subsp. citri

MicroRNAs (miRNAs) are short noncoding RNA molecules that regulate gene expression at the posttranscriptional level. Reverse transcription-quantitative PCR (RT-qPCR) is one of the most common methods used for quantification of miRNA expression, and the levels of expression are normalized by comparing with reference genes. Thus, the selection of reference genes is critically important for accurate quantification. The present study was intended to identify appropriate miRNA reference genes for normalizing the level of miRNA expression in Citrus sinensis L. Osbeck and Citrus reticulata Blanco infected by Xanthomonas citri subsp. citri, which caused citrus canker disease. Five algorithms (Delta Ct, geNorm, NormFinder, BestKeeper and RefFinder) were used for screening reference genes, and two quantification approaches, poly(A) extension RT-qPCR and stem-loop RT-qPCR, were used to determine the most appropriate method for detecting expression patterns of miRNA. An overall comprehensive ranking output derived from the multi-algorithms showed that poly(A)-tailed miR162-3p/miR472 were the best reference gene combination for miRNA RT-qPCR normalization in citrus canker research. Candidate reference gene expression profiles determined by poly(A) RT-qPCR were more consistent in the two citrus species. To the best of our knowledge, this is the first systematic comparison of two miRNA quantification methods for evaluating reference genes. These results highlight the importance of rigorously assessing candidate reference genes and clarify some contradictory results in miRNA research on citrus.

Integrated mRNA and small RNA sequencing reveals microRNA regulatory network associated with internode elongation in sugarcane (Saccharum officinarum L.)

Background

Internode elongation is one of the most important traits in sugarcane because of its relation to crop productivity. Understanding the microRNA (miRNA) and mRNA expression profiles related to sugarcane internode elongation would help develop molecular improvement strategies but they are not yet well-investigated. To identify genes and miRNAs involved in internode elongation, the cDNA and small RNA libraries from the pre-elongation stage (EI), early elongation stage (EII) and rapid elongation stage (EIII) were sequenced and their expression were studied.

Results

Based on the sequencing results, 499,495,518 reads and 80,745 unigenes were identified from stem internodes of sugarcane. The comparisons of EI vs. EII, EI vs. EIII, and EII vs. EIII identified 493, 5035 and 3041 differentially expressed genes, respectively. Further analysis revealed that the differentially expressed genes were enriched in the GO terms oxidoreductase activity and tetrapyrrole binding. KEGG pathway annotation showed significant enrichment in “zeatin biosynthesis”, “nitrogen metabolism” and “plant hormone signal transduction”, which might be participating in internode elongation. miRNA identification showed 241 known miRNAs and 245 novel candidate miRNAs. By pairwise comparison, 11, 42 and 26 differentially expressed miRNAs were identified from EI and EII, EI and EIII, and EII and EIII comparisons, respectively. The target prediction revealed that the genes involved in “zeatin biosynthesis”, “nitrogen metabolism” and “plant hormone signal transduction” pathways are targets of the miRNAs. We found that the known miRNAs miR2592-y, miR1520-x, miR390-x, miR5658-x, miR6169-x and miR8154-x were likely regulators of genes with internode elongation in sugarcane.

Conclusions

The results of this study provided a global view of mRNA and miRNA regulation during sugarcane internode elongation. A genetic network of miRNA-mRNA was identified with miRNA-mediated gene expression as a mechanism in sugarcane internode elongation. Such evidence will be valuable for further investigations of the molecular regulatory mechanisms underpinning sugarcane growth and development.

RNA-seq-based selection of reference genes for RT-qPCR analysis of pitaya

Reverse‐transcription quantitative real‐time PCR (RT‐qPCR) is a primary tool for measuring gene expression levels, and selection of appropriate reference genes is crucial for accurate and reproducible results of gene expression under various experimental conditions. However, no systematic evaluation of reference genes in pitaya (Hylocereus undatus Britt.) has been performed. Here, we examined the expression of five candidate reference genes, namely elongation factor 1‐alpha (HuEF1‐α), 18S ribosomal RNA (Hu18S rRNA), ubiquitin (HuUBQ), actin (HuACT), and ubiquitin‐conjugating enzyme (HuUQT), under different conditions in pitaya. The expression stabilities of these five genes were evaluated using two computation programs: geNorm and NormFinder. The results were further validated by normalizing the expression of the phosphoglycerate kinase (HuPGK) and ethylene‐responsive transcription factor (HuERF) genes. Our results indicate that combined use of HuUBQ and HuUQT is the most stable reference under all of the experimental conditions examined. HuEF1‐α, HuUBQ, and HuUQT are the top three most stable reference genes under salt stress, drought stress, and heat stress, and across different cultivars. HuEF1‐α, HuACT, and HuUQT exhibited the most stable expression patterns across different tissues. Our results will allow researchers to select the most appropriate reference genes for gene expression studies of pitaya under different conditions.

OsmiR528 Enhances Cold Stress Tolerance by Repressing Expression of Stress Response-related Transcription Factor Genes in Plant Cells

Background:

MicroRNAs participate in many molecular mechanisms and signaling trans-duction pathways that are associated with plant stress tolerance by repressing expression of their target genes. However, how microRNAs enhance tolerance to low temperature stress in plant cells remains elusive.

Objective:

In this investigation, we demonstrated that overexpression of the rice microRNA528 (Os-miR528) increases cell viability, growth rate, antioxidants content, ascorbate peroxidase (APOX) activi-ty, and superoxide dismutase (SOD) activity and decreases ion leakage rate and thiobarbituric acid reac-tive substances (TBARS) under low temperature stress in Arabidopsis (Arabidopsis thaliana), pine (Pi-nus elliottii), and rice (Oryza sativa).

Methods:

To investigate the potential mechanism of OsmiR528 in increasing cold stress tolerance, we examined expression of stress-associated MYB transcription factors OsGAMYB-like1, OsMYBS3, OsMYB4, OsMYB3R-2, OsMYB5, OsMYB59, OsMYB30, OsMYB1R, and OsMYB20 in rice cells by qRT-PCR.

Results:

Our experiments demonstrated that OsmiR528 decreases expression of transcription factor OsMYB30 by targeting a F-box domain containing protein gene (Os06g06050), which is a positive regulator of OsMYB30. In OsmiR528 transgenic rice, reduced OsMYB30 expression results in in-creased expression of BMY genes OsBMY2, OsBMY6, and OsBMY10. The transcript levels of the OsBMY2, OsBMY6, and OsBMY10 were elevated by OsMYB30 knockdown, but decreased by Os-MYB30 overexpression in OsmiR528 transgenic cell lines, suggesting that OsmiR528 increases low temperature tolerance by modulating expression of stress response-related transcription factor.

Conclusion:

Our experiments provide novel information in increasing our understanding in molecular mechanisms of microRNAs-associated low temperature tolerance and are valuable in plant molecular breeding from monocotyledonous, dicotyledonous, and gymnosperm plants.

Selection and validation of reference genes for quantitative expression analysis of miRNAs and mRNAs in Poplar

BACKGROUND

Quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) is a rapid and sensitive approach to identify miRNA and protein-coding gene expression in plants. However, because of the specially designated reverse transcription and shorter PCR products, very few reference genes have been identified for the quantitative analysis of miRNA expression in plants, and different internal reference genes are needed to normalize the expression of miRNAs and mRNA genes respectively. Therefore, it is particularly important to select the suitable common reference genes for normalization of quantitative PCR of miRNA and mRNA.

RESULTS

In this study, a modified reverse transcription PCR protocol was adopted for selecting and validating universal internal reference genes of mRNAs and miRNAs. Eight commonly used reference genes, four stably expressed novel genes in Populus tremula, three small noncoding RNAs and three conserved miRNAs were selected as candidate genes, and the stability of their expression was examined across a set of 38 tissue samples from four developmental stages of poplar clone 84K (Populus alba × Populus glandulosa). The expression stability of these candidate genes was evaluated systematically by four algorithms: geNorm, NormFinder, Bestkeeper and DeltaCt. The results showed that Eukaryotic initiation factor 4A III (EIF4A) and U6-2 were suitable for samples of the callus stage; U6-1 and U6-2 were best for the seedling stage; Protein phosphatase 2A-2 (PP2A-2) and U6-1 were best for the plant stage; and Protein phosphatase 2A-2 (PP2A-2) and Oligouridylate binding protein 1B (UBP) were the best reference genes in the adventitious root (AR) regeneration stage.

CONCLUSIONS

The purpose of this study was to identify the most appropriate reference genes for qRT-PCR of miRNAs and mRNAs in different tissues at several developmental stages in poplar. U6-1, EIF4A and PP2A-2 were the three most appropriate reference genes for qRT-PCR normalization of miRNAs and mRNAs during the plant regeneration process, and PP2A-2 and UBP represent the best reference genes in the AR regeneration stage of poplar. This work will benefit future studies of expression and function analysis of miRNAs and their target genes in poplar.

Identifying optimal reference genes for the normalization of microRNA expression in cucumber under viral stress

Cucumber green mottle mosaic virus (CGMMV) is an economically important pathogen and causes significant reduction of both yield and quality of cucumber (Cucumis sativus). Currently, there were no satisfied strategies for controlling the disease. A better understanding of microRNA (miRNA) expression related to the regulation of plant-virus interactions and virus resistance would be of great assistance when developing control strategies for CGMMV. However, accurate expression analysis is highly dependent on robust and reliable reference gene used as an internal control for normalization of miRNA expression. Most commonly used reference genes involved in CGMMV-infected cucumber are not universally expressed depending on tissue types and stages of plant development. It is therefore crucial to identify suitable reference genes in investigating the role of miRNA expression. In this study, seven reference genes, including Actin, Tubulin, EF-1α, 18S rRNA, Ubiquitin, GAPDH and Cyclophilin, were evaluated for the most accurate results in analyses using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Gene expression was assayed on cucumber leaves, stems and roots that were collected at different days post inoculation with CGMMV. The expression data were analyzed using algorithms including delta-Ct, geNorm, NormFinder, and BestKeeper as well as the comparative tool RefFinder. The reference genes were subsequently validated using miR159. The results showed that EF-1α and GAPDH were the most reliable reference genes for normalizing miRNA expression in leaf, root and stem samples, while Ubiquitin and EF-1α were the most suitable combination overall.

Selection of reference genes for miRNA qRT-PCR under abiotic stress in grapevine

Grapevine is among the fruit crops with high economic value, and because of the economic losses caused by abiotic stresses, the stress resistance of Vitis vinifera has become an increasingly important research area. Among the mechanisms responding to environmental stresses, the role of miRNA has received much attention recently. qRT-PCR is a powerful method for miRNA quantitation, but the accuracy of the method strongly depends on the appropriate reference genes. To determine the most suitable reference genes for grapevine miRNA qRT-PCR, 15 genes were chosen as candidate reference genes. After eliminating 6 candidate reference genes with unsatisfactory amplification efficiency, the expression stability of the remaining candidate reference genes under salinity, cold and drought was analysed using four algorithms, geNorm, NormFinder, deltaCt and Bestkeeper. The results indicated that U6 snRNA was the most suitable reference gene under salinity and cold stresses; whereas miR168 was the best for drought stress. The best reference gene sets for salinity, cold and drought stresses were miR160e + miR164a, miR160e + miR168 and ACT + UBQ + GAPDH, respectively. The selected reference genes or gene sets were verified using miR319 or miR408 as the target gene.

miRNA alteration is an important mechanism in sugarcane response to low-temperature environment

BACKGROUND

Cold is a major abiotic stress limiting the production of tropical and subtropical crops in new production areas. Sugarcane (Saccharum spp.) originates from the tropics but is cultivated primarily in the sub-tropics where it frequently encounters cold stress. Besides regulating plant growth, miRNAs play an important role in environmental adaption.

RESULTS

In this study, a total of 412 sugarcane miRNAs, including 261 known and 151 novel miRNAs, were obtained from 4 small RNA libraries through the Illumina sequencing method. Among them, 62 exhibited significant differential expression under cold stress, with 34 being upregulated and 28 being downregulated. The expression of 13 miRNAs and 12 corresponding targets was validated by RT-qPCR, with the majority being consistent with the sequencing data. GO and KEGG analysis indicated that these miRNAs were involved in stress-related biological pathways. To further investigate the involvement of these miRNAs in tolerance to abiotic stresses, sugarcane miR156 was selected for functional analysis. RT-qPCR revealed that miR156 levels increased in sugarcane during cold, salt and drought stress treatments. Nicotiana benthamiana plants transiently overexpressing miR156 exhibited better growth status, lower ROS levels, higher anthocyanin contents as well as the induction of some cold-responsive genes, suggesting its positive role in the plant cold stress response.

CONCLUSIONS

This study provides a global view of the association of miRNA expression with the sugarcane response to cold stress. The findings have enriched the present miRNA resource and have made an attempt to verify the involvement of miR156 in plant response to cold stress.

De novo characterization of microRNAs in oriental fruit moth Grapholita molesta and selection of reference genes for normalization of microRNA expression

MicroRNAs (miRNAs) are a group of endogenous non-coding small RNAs that have critical regulatory functions in almost all known biological processes at the post-transcriptional level in a variety of organisms. The oriental fruit moth Grapholita molesta is one of the most serious pests in orchards worldwide and threatens the production of Rosacea fruits. In this study, a de novo small RNA library constructed from mixed stages of G. molesta was sequenced through Illumina sequencing platform and a total of 536 mature miRNAs consisting of 291 conserved and 245 novel miRNAs were identified. Most of the conserved and novel miRNAs were detected with moderate abundance. The miRNAs in the same cluster normally showed correlated expressional profiles. A comparative analysis of the 79 conserved miRNA families within 31 arthropod species indicated that these miRNA families were more conserved among insects and within orders of closer phylogenetic relationships. The KEGG pathway analysis and network prediction of target genes indicated that the complex composed of miRNAs, clock genes and developmental regulation genes may play vital roles to regulate the developmental circadian rhythm of G. molesta. Furthermore, based on the sRNA library of G. molesta, suitable reference genes were selected and validated for study of miRNA transcriptional profile in G. molesta under two biotic and six abiotic experimental conditions. This study systematically documented the miRNA profile in G. molesta, which could lay a foundation for further understanding of the regulatory roles of miRNAs in the development and metabolism in this pest and might also suggest clues to the development of genetic-based techniques for agricultural pest control.