Integrated analysis of microRNA-target interactions with clinical outcomes for cancers

Background

Clinical statement alone is not enough to predict the progression of disease. Instead, the gene expression profiles have been widely used to forecast clinical outcomes. Many genes related to survival have been identified, and recently miRNA expression signatures predicting patient survival have been also investigated for several cancers. However, miRNAs and their target genes associated with clinical outcomes have remained largely unexplored.

Methods

Here, we demonstrate a survival analysis based on the regulatory relationships of miRNAs and their target genes. The patient survivals for the two major cancers, ovarian cancer and glioblastoma multiforme (GBM), are investigated through the integrated analysis of miRNA-mRNA interaction pairs.

Results

We found that there is a larger survival difference between two patient groups with an inversely correlated expression profile of miRNA and mRNA. It supports the idea that signatures of miRNAs and their targets related to cancer progression can be detected via this approach.

Conclusions

This integrated analysis can help to discover coordinated expression signatures of miRNAs and their target mRNAs that can be employed for therapeutics in human cancers.

     

SAGE detects microRNA precursors

Background

MicroRNAs (miRNAs) have been shown to play important roles in regulating gene expression. Since miRNAs are often evolutionarily conserved and their precursors can be folded into stem-loop hairpins, many miRNAs have been predicted. Yet experimental confirmation is difficult since miRNA expression is often specific to particular tissues and developmental stages.

Results

Analysis of 29 human and 230 mouse longSAGE libraries revealed the expression of 22 known and 10 predicted mammalian miRNAs. Most were detected in embryonic tissues. Four SAGE tags detected in human embryonic stem cells specifically match a cluster of four human miRNAs (mir-302a, b, c&d) known to be expressed in embryonic stem cells. LongSAGE data also suggest the existence of a mouse homolog of human and rat mir-493.

Conclusion

The observation that some orphan longSAGE tags uniquely match miRNA precursors provides information about the expression of some known and predicted miRNAs.

     

A novel two-layer SVM model in miRNA Drosha processing site detection

Background

MicroRNAs (miRNAs) are a large class of non-coding RNAs with important functions wide spread in animals, plants and viruses. Studies showed that an RNase III family member called Drosha recognizes most miRNAs, initiates their processing and determines the mature miRNAs. The Drosha processing sites identification will shed some light on both miRNA identification and understanding the mechanism of Drosha processing.

Methods

We developed a computational method for Drosha processing site predicting, named as DroshaPSP, which employs a two-layer mathematical model to integrate structure feature in the first layer and sequence features in the second layer. The performance of DroshaPSP was estimated by 5-fold cross-validation and measured by ACC (accuracy), Sn (sensitivity), Sp (specificity), P (precision) and MCC (Matthews correlation coefficient).

Results

The results of testing DroshaPSP on the miRNA data of Drosophila melanogaster indicated that the Sn, Sp, and MCC thereof reach to 0.86, 0.99 and 0.86 respectively.

Conclusions

We found the Shannon entropy, a chemical kinetics feature, is a significant feature in telling the true sites among the nearby sites and improving the performance.

     

Microrna expression signatures predict patient progression and disease outcome in pediatric embryonal central nervous system neoplasms

Background

Although, substantial experimental evidence related to diagnosis and treatment of pediatric central nervous system (CNS) neoplasms have been demonstrated, the understanding of the etiology and pathogenesis of the disease remains scarce. Recent microRNA (miRNA)-based research reveals the involvement of miRNAs in various aspects of CNS development and proposes that they might compose key molecules underlying oncogenesis. The current study evaluated miRNA differential expression detected between pediatric embryonal brain tumors and normal controls to characterize candidate biomarkers related to diagnosis, prognosis and therapy.

Methods

Overall, 19 embryonal brain tumors; 15 Medulloblastomas (MBs) and 4 Atypical Teratoid/Rabdoid Tumors (AT/RTs) were studied. As controls, 13 samples were used; The First-Choice Human Brain Reference RNA and 12 samples from deceased children who underwent autopsy and were not present with any brain malignancy. RNA extraction was carried out using the Trizol method, whilst miRNA extraction was performed with the mirVANA miRNA isolation kit. The experimental approach included miRNA microarrays covering 1211 miRNAs. Quantitative Real-Time Polymerase Chain Reaction was performed to validate the expression profiles of miR-34a and miR-601 in all 32 samples initially screened with miRNA microarrays and in an additional independent cohort of 30 patients (21MBs and 9 AT/RTs). Moreover, meta-analyses was performed in total 27 embryonal tumor samples; 19 MBs, 8 ATRTs and 121 control samples. Twelve germinomas were also used as an independent validation cohort. All deregulated miRNAs were correlated to patients’ clinical characteristics and pathological measures.

Results

In several cases, there was a positive correlation between individual miRNA expression levels and laboratory or clinical characteristics. Based on that, miR-601 could serve as a putative tumor suppressor gene, whilst miR-34a as an oncogene. In general, miR-34a demonstrated oncogenic roles in all pediatric embryonal CNS neoplasms studied.

Conclusions

Deeper understanding of the aberrant miRNA expression in pediatric embryonal brain tumors might aid in the development of tumor-specific miRNA signatures, which could potentially afford promising biomarkers related to diagnosis, prognosis and patient targeted therapy.

     

Sex differences in the expression of lupus-associated miRNAs in splenocytes from lupus-prone NZB/W<Subscript>F1</Subscript> mice

Background

A majority of autoimmune diseases, including systemic lupus erythematosus (SLE), occur predominantly in females. Recent studies have identified specific dysregulated microRNAs (miRNAs) in both human and murine lupus, implying an important contribution of these miRNAs to lupus pathogenesis. However, to date, there is no study that examined sex differences in miRNA expression in immune cells as a plausible basis for sex differences in autoimmune disease. This study addresses this aspect in NZB/W_F1 mice, a classical murine lupus model with marked female bias, and further investigates estrogen regulation of lupus-associated miRNAs.

Methods

The Taqman miRNA assay system was used to quantify the miRNA expression in splenocytes from male and female NZB/W_F1 mice at 17–18, 23, and 30 weeks (wks) of age. To evaluate potential estrogen's effect on lupus-associated miRNAs, 6-wk-old NZB/W_F1 male mice were orchidectomized and surgically implanted with empty (placebo) or estrogen implants for 4 and 26 wks, respectively. To assess the lupus status in the NZB/W_F1 mice, serum anti-dsDNA autoantibody levels, proteinuria, and renal histological changes were determined.

Results

The sex differences in the expression of lupus-associated miRNAs, including the miR-182-96-183 cluster, miR-155, miR-31, miR-148a, miR-127, and miR-379, were markedly evident after the onset of lupus, especially at 30 wks of age when female NZB/W_F1 mice manifested moderate to severe lupus when compared to their male counterparts. Our limited data also suggested that estrogen treatment increased the expression of aforementioned lupus-associated miRNAs, with the exception of miR-155, in orchidectomized male NZB/W_F1 mice to a similar level in age-matched intact female NZB/W_F1 mice. It is noteworthy that orchiectomy, itself, did not affect the expression of lupus-associated miRNAs.

Conclusion

To our knowledge, this is the first study that demonstrated sex differences in the expression of lupus-associated miRNAs in splenocytes, especially in the context of autoimmunity. The increased expression of lupus-associated miRNA in female NZB/W_F1 mice and conceivably in estrogen-treated orchidectomized male NZB/W_F1 mice was associated with lupus manifestation. The notable increase of lupus-associated miRNAs in diseased, female NZB/W_F1 mice may be a result of both lupus manifestation and the female gender.

     

Profiling microRNAs in individuals at risk of progression to rheumatoid arthritis

Background

Individuals at risk of rheumatoid arthritis (RA) demonstrate systemic autoimmunity in the form of anti-citrullinated peptide antibodies (ACPA). MicroRNAs (miRNAs) are implicated in established RA. This study aimed to (1) compare miRNA expression between healthy individuals and those at risk of and those that develop RA, (2) evaluate the change in expression of miRNA from “at-risk” to early RA and (3) explore whether these miRNAs could inform a signature predictive of progression from “at-risk” to RA.

Methods

We performed global profiling of 754 miRNAs per patient on a matched serum sample cohort of 12 anti-cyclic citrullinated peptide (CCP)?+?“at-risk” individuals that progressed to RA. Each individual had a serum sample from baseline and at time of detection of synovitis, forming the matched element. Healthy controls were also studied. miRNAs with a fold difference/fold change of four in expression level met our primary criterion for selection as candidate miRNAs. Validation of the miRNAs of interest was conducted using custom miRNA array cards on matched samples (baseline and follow up) in 24 CCP+ individuals; 12 RA progressors and 12 RA non-progressors.

Results

We report on the first study to use matched serum samples and a comprehensive miRNA array approach to identify in particular, three miRNAs (miR-22, miR-486-3p, and miR-382) associated with progression from systemic autoimmunity to RA inflammation. MiR-22 demonstrated significant fold difference between progressors and non-progressors indicating a potential biomarker role for at-risk individuals.

Conclusions

This first study using a cohort with matched serum samples provides important mechanistic insights in the transition from systemic autoimmunity to inflammatory disease for future investigation, and with further evaluation, might also serve as a predictive biomarker.

     

Estimating survival time of patients with glioblastoma multiforme and characterization of the identified microRNA signatures

Background

Though glioblastoma multiforme (GBM) is the most frequently occurring brain malignancy in adults, clinical treatment still faces challenges due to poor prognoses and tumor relapses. Recently, microRNAs (miRNAs) have been extensively used with the aim of developing accurate molecular therapies, because of their emerging role in the regulation of cancer-related genes. This work aims to identify the miRNA signatures related to survival of GBM patients for developing molecular therapies.

Results

This work proposes a support vector regression (SVR)-based estimator, called SVR-GBM, to estimate the survival time in patients with GBM using their miRNA expression profiles. SVR-GBM identified 24 out of 470 miRNAs that were significantly associated with survival of GBM patients. SVR-GBM had a mean absolute error of 0.63 years and a correlation coefficient of 0.76 between the real and predicted survival time. The 10 top-ranked miRNAs according to prediction contribution are as follows: hsa-miR-222, hsa-miR-345, hsa-miR-587, hsa-miR-526a, hsa-miR-335, hsa-miR-122, hsa-miR-24, hsa-miR-433, hsa-miR-574 and hsa-miR-320. Biological analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway on the identified miRNAs revealed their influence in GBM cancer.

Conclusion

The proposed SVR-GBM using an optimal feature selection algorithm and an optimized SVR to identify the 24 miRNA signatures associated with survival of GBM patients. These miRNA signatures are helpful to uncover the individual role of miRNAs in GBM prognosis and develop miRNA-based therapies.

     

Comprehensive annotation of microRNA expression profiles

Background

MicroRNAs (miRNAs) regulate many biological processes by post-translational gene silencing. Analysis of miRNA expression profiles is a reliable method for investigating particular biological processes due to the stability of miRNA and the development of advanced sequencing methods. However, this approach is limited by the broad specificity of miRNAs, which may target several mRNAs.

Result

In this study, we developed a method for comprehensive annotation of miRNA array or deep sequencing data for investigation of cellular biological effects. Using this method, the specific pathways and biological processes involved in Alzheimer’s disease were predicted with high correlation in four independent samples. Furthermore, this method was validated for evaluation of cadmium telluride (CdTe) nanomaterial cytotoxicity. As a result, apoptosis pathways were selected as the top pathways associated with CdTe nanoparticle exposure, which is consistent with previous studies.

Conclusions

Our findings contribute to the validation of miRNA microarray or deep sequencing results for early diagnosis of disease and evaluation of the biological safety of new materials and drugs.

     

Upregulation of microRNA 344a-3p is involved in curcumin induced apoptosis in RT4 schwannoma cells

Background

Schwannoma arising from peripheral nervous sheaths is a benign tumor.

Methods

To evaluate cell cytotoxicity, (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) tetrazolium reduction and terminal deoxynucleotidyltransferase UTP nick-end labeling (TUNEL) assays were used. A microRNA (miRNA) array was used to identify the miRNAs involved in curcumin-induced apoptosis. To examine miRNA expression, quantitative RT-PCR was used.

Results

In this study, curcumin exerted cellular cytotoxicity against RT4 schwannoma cells, with an increase in TUNEL-positive cells. Curcumin also activated the expression of apoptotic proteins, such as polyADP ribose polymerase, caspase-3, and caspase-9. The miRNA array revealed that seven miRNAs (miRNA 350, miRNA 17-2-3p, let 7e-3p, miRNA1224, miRNA 466b-1-3p, miRNA 18a-5p, and miRNA 322-5p) were downregulated following treatment with both 10 and 20 μM curcumin in RT4 cells, while four miRNAs (miRNA122-5p, miRNA 3473, miRNA182, and miRNA344a-3p) were upregulated. Interestingly, transfection with a miRNA 344a-3p mimic downregulated the mRNA expression of Bcl2 and upregulated that of Bax, Curcumin treatment in RT 4 cells also reduced the mRNA expression of Bcl2 and enhanced expression of Bax, Overexpression of miRNA344a-3p mimic combined with curcumin treatment activated the expression of apoptotic proteins, including procaspase-9 and cleaved caspase-3 while inhibition of miRNA 344a-3p using miR344a-3p inhibitor repressed cleaved caspase-3 and -9 in curcumin treated RT-4 cells compared to control.

Conclusions

Our findings demonstrate that curcumin induces apoptosis in schwannoma cells via miRNA 344a-3p. Thus, curcumin may serve as a potent therapeutic agent for the treatment of schwannoma.

     

Androgen and estrogen receptors in placental physiology and dysfunction

Background

The placenta is recognized as an endocrine organ, largely due to its secretions of steroid hormones, including progesterone, androgens, and estrogens. Steroid hormones play an essential role in the progression of pregnancy, fetal development, and growth. Furthermore, steroids are necessary for establishment and maintenance of a normal pregnancy, preparing the endometrium for implantation, stimulating endometrial secretions, and regulating uterine blood flow, however the exact mechanism of sex steroid signaling through their receptors in placental function is unknown.

Objective

In this review, we will provide an overview of the current knowledge on sex steroid receptors in normal placental development, as well as evidence of abnormal signaling associated with placental dysfunction.

Methods

A systematic literature search was performed using the NCBI PubMed search engine, including the following key works: estrogen receptor, androgen receptor, placenta, placental development, cytotrophoblast, and differentiation.

Results

Of the over 700 articles that were returned, 125 studies focused on estrogen and androgen receptors in human placenta development and function during normal and abnormal pregnancy, as well as in rodents and ruminants placentae.

Conclusion

Receptors for both estrogens and androgens have been localized within the mammalian placenta, but surprisingly little is known about their signaling in trophoblast cell differentiation and function. An emerging picture is developing in which estrogen receptors possibly play role in cytotrophoblast proliferation and extravillous trophoblast invasion, whereas androgen receptors are involved in syncytiotrophoblast differentiation and function.

     

Bioavailability of transgenic microRNAs in genetically modified plants

Background

Transgenic expression of small RNAs is a prevalent approach in agrobiotechnology for the global enhancement of plant foods. Meanwhile, emerging studies have, on the one hand, emphasized the potential of transgenic microRNAs (miRNAs) as novel dietary therapeutics and, on the other, suggested potential food safety issues if harmful miRNAs are absorbed and bioactive. For these reasons, it is necessary to evaluate the bioavailability of transgenic miRNAs in genetically modified crops.

Results

As a pilot study, two transgenic Arabidopsis lines ectopically expressing unique miRNAs were compared and contrasted with the plant bioavailable small RNA MIR2911 for digestive stability and serum bioavailability. The expression levels of these transgenic miRNAs in Arabidopsis were found to be comparable to that of MIR2911 in fresh tissues. Assays of digestive stability in vitro and in vivo suggested the transgenic miRNAs and MIR2911 had comparable resistance to degradation. Healthy mice consuming diets rich in Arabidopsis lines expressing these miRNAs displayed MIR2911 in the bloodstream but no detectable levels of the transgenic miRNAs.

Conclusions

These preliminary results imply digestive stability and high expression levels of miRNAs in plants do not readily equate to bioavailability. This initial work suggests novel engineering strategies be employed to enhance miRNA bioavailability when attempting to use transgenic foods as a delivery platform.

     

Novel link prediction for large-scale miRNA-lncRNA interaction network in a bipartite graph

Background

Current knowledge and data on miRNA-lncRNA interactions is still limited and little effort has been made to predict target lncRNAs of miRNAs. Accumulating evidences suggest that the interaction patterns between lncRNAs and miRNAs are closely related to relative expression level, forming a titration mechanism. It could provide an effective approach for characteristic feature extraction. In addition, using the coding non-coding co-expression network and sequence data could also help to measure the similarities among miRNAs and lncRNAs. By mathematically analyzing these types of similarities, we come up with two findings that (i) lncRNAs/miRNAs tend to collaboratively interact with miRNAs/lncRNAs of similar expression profiles, and vice versa, and (ii) those miRNAs interacting with a cluster of common target genes tend to jointly target at the common lncRNAs.

Methods

In this work, we developed a novel group preference Bayesian collaborative filtering model called GBCF for picking up a top-k probability ranking list for an individual miRNA or lncRNA based on the known miRNA-lncRNA interaction network.

Results

To evaluate the effectiveness of GBCF, leave-one-out and k-fold cross validations as well as a series of comparison experiments were carried out. GBCF achieved the values of area under ROC curve of 0.9193, 0.8354+/??0.0079, 0.8615+/??0.0078, and 0.8928+/??0.0082 based on leave-one-out, 2-fold, 5-fold, and 10-fold cross validations respectively, demonstrating its reliability and robustness.

Conclusions

GBCF could be used to select potential lncRNA targets of specific miRNAs and offer great insights for further researches on ceRNA regulation network.

     

Deregulation of neuronal miRNAs induced by amyloid-β or TAU pathology

Background

Despite diverging levels of amyloid-β (Aβ) and TAU pathology, different mouse models, as well as sporadic AD patients show predictable patterns of episodic memory loss. MicroRNA (miRNA) deregulation is well established in AD brain but it is unclear whether Aβ or TAU pathology drives those alterations and whether miRNA changes contribute to cognitive decline.

Methods

miRNAseq was performed on cognitively intact (4 months) and impaired (10 months) male APPtg (APP^swe/PS1^L166P) and TAUtg (THY-Tau22) mice and their wild-type littermates (APPwt and TAUwt). We analyzed the hippocampi of 12 mice per experimental group (n =?96 in total), and employed a 2-way linear model to extract differentially expressed miRNAs. Results were confirmed by qPCR in a separate cohort of 4 M and 10 M APPtg and APPwt mice (n =?7–9 per group) and in human sporadic AD and non-demented control brain. Fluorescent in situ hybridization identified their cellular expression. Functional annotation of predicted targets was performed using GO enrichment. Behavior of wild-type mice was assessed after intracerebroventricular infusion of miRNA mimics.

Results

Six miRNAs (miR-10a-5p, miR-142a-5p, miR-146a-5p, miR-155-5p, miR-211-5p, miR-455-5p) are commonly upregulated between APPtg and TAUtg mice, and four of these (miR-142a-5p, miR-146a-5p, miR-155-5p and miR-455-5p) are altered in AD patients. All 6 miRNAs are strongly enriched in neurons. Upregulating these miRNAs in wild-type mice is however not causing AD-related cognitive disturbances.

Conclusion

Diverging AD-related neuropathologies induce common disturbances in the expression of neuronal miRNAs. 4 of these miRNAs are also upregulated in AD patients. Therefore these 4 miRNAs (miR-142a-5p, miR-146a-5p, miR-155-5p and miR-455-5p) appear part of a core pathological process in AD patients and APPtg and TAUtg mice. They are however not causing cognitive disturbances in wild-type mice. As some of these miRNA target AD relevant proteins, they may be, in contrast, part of a protective response in AD.

     

FMSM: a novel computational model for predicting potential miRNA biomarkers for various human diseases

Background

MicroRNA (miRNA) plays a key role in regulation mechanism of human biological processes, including the development of disease and disorder. It is necessary to identify potential miRNA biomarkers for various human diseases. Computational prediction model is expected to accelerate the process of identification.

Results

Considering the limitations of previously proposed models, we present a novel computational model called FMSM. It infers latent miRNA biomarkers involved in the mechanism of various diseases based on the known miRNA-disease association network, miRNA expression similarity, disease semantic similarity and Gaussian interaction profile kernel similarity. FMSM achieves reliable prediction performance in 5-fold and leave-one-out cross validations with area under ROC curve (AUC) values of 0.9629+/??0.0127 and 0.9433, respectively, which outperforms the state-of-the-art competitors and classical algorithms. In addition, 19 of top 25 predicted miRNAs have been validated to have associations with Colonic Neoplasms in case study.

Conclusions

A factored miRNA similarity based model and miRNA expression similarity substantially contribute to the well-performing prediction. The list of the predicted most latent miRNA biomarkers of various human diseases is publicized. It is anticipated that FMSM could serve as a useful tool guiding the future experimental validation for those promising miRNA biomarker candidates.

     

Grouping miRNAs of similar functions via weighted information content of gene ontology

Background

Regulation mechanisms between miRNAs and genes are complicated. To accomplish a biological function, a miRNA may regulate multiple target genes, and similarly a target gene may be regulated by multiple miRNAs. Wet-lab knowledge of co-regulating miRNAs is limited. This work introduces a computational method to group miRNAs of similar functions to identify co-regulating miRNAsfrom a similarity matrix of miRNAs.

Results

We define a novel information content of gene ontology (GO) to measure similarity between two sets of GO graphs corresponding to the two sets of target genes of two miRNAs. This between-graph similarity is then transferred as a functional similarity between the two miRNAs. Our definition of the information content is based on the size of a GO term’s descendants, but adjusted by a weight derived from its depth level and the GO relationships at its path to the root node or to the most informative common ancestor (MICA). Further, a self-tuning technique and the eigenvalues of the normalized Laplacian matrix are applied to determine the optimal parameters for the spectral clustering of the similarity matrix of the miRNAs.

Conclusions

Experimental results demonstrate that our method has better clustering performance than the existing edge-based, node-based or hybrid methods. Our method has also demonstrated a novel usefulness for the function annotation of new miRNAs, as reported in the detailed case studies.

     

Association between microRNA-21, microRNA-150, and micro-RNA-451 expression and clinical outcome of patients with acute lymphoblastic leukemia

Background

Acute lymphoblastic leukemia (ALL) occurs owing to the defective maturation, increased proliferation, and lack of differentiation of lymphoid cells. Evaluation of the expression levels of microRNAs (miRNAs) could help in the prognosis and improve the clinical outcome of ALL patients. Given the role of miR-21, miR-150, and miR-451 as oncogenes and tumor suppressors in lymphocytes, this study explored the relation between the expression levels of these miRNAs and the clinical outcomes of ALL patients.

Methods

cDNA synthesis and RT-PCR were performed for peripheral blood samples from 41 patients with ALL, as well as for U937 and Jurkat cell lines to examine the expression of miR-451, miR-150, and miR-21 after miRNA purification. We also performed an epidemiological analysis in which Mann–Whitney and Chi-square tests were used to investigate the relationship between the expression of miRNAs and other clinical and laboratory data. Binary logistic regression models were used to estimate the odds ratio in univariate and multivariate analyses for clinical outcomes.

Results

miR-21 and miR-150 expression was found to be decreased, while miR-451 expression showed no difference compared to the control group. There was a significant relationship between miR-451 expression and hemoglobin (Hb) levels, as well as between miR-150 expression and clinical outcomes of ALL patients.

Conclusion

Increased expression of miR-451 decreased the Hb levels; reduced expression of miR-150 was associated with increased relapse rate in patients. Age, increased WBC, and decreased Hb levels were associated with increased relapse rates in ALL patients. Therefore, miR-150 could be used as a biomarker to determine the clinical outcome of ALL patients.

     

The miRNAome of ramie (<Emphasis Type="Italic">Boehmeria nivea</Emphasis> L.): identification,expression, and potential roles of novel microRNAs in regulation of cadmium stress response

Background

MicroRNAs (miRNAs) regulate numerous crucial abiotic stress processes in plants. However, information is limited on their involvement in cadmium (Cd) stress response and tolerance mechanisms in plants, including ramie (Boehmeria nivea L.) that produces a number of economic valuable as an important natural fibre crop and an ideal crop for Cd pollution remediation.

Results

Four small RNA libraries of Cd-stressed and non-stressed leaves and roots of ramie were constructed. Using small RNA-sequencing, 73 novel miRNAs were identified. Genome-wide expression analysis revealed that a set of miRNAs was differentially regulated in response to Cd stress. In silico target prediction identified 426 potential miRNA targets that include several uptake or transport factors for heavy metal ions. The reliability of small RNA sequencing and the relationship between the expression levels of miRNAs and their target genes were confirmed by quantitative PCR (q-PCR). We showed that the expression patterns of miRNAs obtained by q-PCR were consistent with those obtained from small RNA sequencing. Moreover, we demonstrated that the expression of six randomly selected target genes was inversely related to that of their corresponding miRNAs, indicating that the miRNAs regulate Cd stress response in ramie.

Conclusions

This study enriches the number of Cd-responsive miRNAs and lays a foundation for the elucidation of the miRNA-mediated regulatory mechanism in ramie during Cd stress.