RNA sequencing reveals a depletion of collagen targeting microRNAs in Dupuytren’s disease
© Riester et al. 2015
Received: 18 June 2015
Accepted: 20 September 2015
Published: 7 October 2015
Dupuytren’s disease is an inherited disorder in which patients develop fibrotic contractures of the hand. Current treatment strategies include surgical excision or enzymatic digestion of fibrotic tissue. MicroRNAs, which are key posttranscriptional regulators of genes expression, have been shown to play an important regulatory role in disorders of fibrosis. Therefore in this investigation, we apply high throughput next generation RNA sequencing strategies to characterize microRNA expression in diseased and healthy palmar fascia to elucidate molecular mechanisms responsible for pathogenic fibrosis.
We applied high throughput RNA sequencing techniques to quantify the expression of all known human microRNAs in Dupuytren’s and control palmar fascia. MicroRNAs that were differentially expressed between diseased and healthy tissue samples were used for computational target prediction using the bioinformatics tool ComiR. Molecular pathways that were predicted to be differentially expressed based on computational analysis were validated by performing RT-qPCR on RNA extracted from diseased and non-diseased palmar fascia biopsies.
A comparison of microRNAs expressed in Dupuytren’s fascia and control fascia identified 74 microRNAs with a 2-fold enrichment in Dupuytren’s tissue, and 32 microRNAs with enrichment in control fascia. Computational target prediction for differentially expressed microRNAs indicated preferential targeting of collagens and extracellular matrix related proteins in control palmar fascia. RT-qPCR confirmed the decreased expression of microRNA targeted collagens in control palmar fascia tissues.
Control palmar fascia show decreased expression of mRNAs encoding collagens that are preferentially targeted by microRNAs enriched in non-diseased fascia. Thus alterations in microRNA regulatory networks may play an important role in driving the pathogenic fibrosis seen in Dupuytren’s disease via direct regulatory effects on extracellular matrix protein synthesis.
Dupuytren’s fascia and healthy palmar fascia can be distinguished by unique microRNA profiles, which are predicted to preferentially target collagens and other extracellular matrix proteins.
KeywordsDupuytren’s disease microRNA fibrosis RNA sequencing hand
Dupuytren’s disease is a clinically challenging disorder characterized by the formation of fibrotic bands that cause disabling contractures of the hand. If the disease is not treated, fibrosis can lead to significant functional limitations that may even necessitate amputation of the affected fingers. Current treatment strategies attempt to break up constrictive bands of fibrous tissue after collagen deposition either surgically or enzymatically with collagenase. These treatments are costly and carry a significant complication risk and are associated with a high rate of disease recurrence [1–5].
Dupuytren’s disease has a strong genetic basis and most commonly affects individuals of northern European descent . Large scale genome wide association studies have helped improve our understanding of Dupuytren’s disease, however the specific genetic abnormalities that drive disease pathogenesis have remained elusive. A variety of molecular pathways have been implicated in disease pathogenesis including alterations in Wnt signaling and mitochondrial genes [7–9].
MicroRNAs, which are small non-coding RNA molecules (20–24 nucleotides in length) that act as post transcriptional regulators of gene expression by inhibiting the translation of target mRNAs, have been shown to regulate the expression of extracellular matrix proteins in the setting of fibrosis [10–13]. Given that Dupuytren’s disease is characterized by excess collagen deposition and fibrosis, we examined the role of microRNAs as pro-fibrotic drivers of the disease process. In this investigation we applied high throughput molecular sequencing techniques to characterize all known microRNAs expressed in diseased Dupuytren’s fascia, and compared expression profiles to non-diseased palmar fascia. We also utilized differentially expressed microRNAs to identify novel pathways as well as validate mechanisms previously implicated with Dupuytren’s disease.
Dupuytren’s tissue biopsies were collected for research use from patient’s undergoing open palmar fasciectomy for the treatment of Dupuytren’s contracture. Surgical cases clinically represented end stage disease in the consolidation phase. All Dupuytren’s tissue specimens were evaluated under frozen section by trained musculoskeletal pathologists to confirm the diagnosis and to ensure representative areas of diseased tissue were selected. Samples were then snap frozen in liquid nitrogen and stored at −80 °C until use for RNA extraction. Adjacent fascia specimens were obtained from palmar fascia adjacent to diseased Dupuytren’s fascia that was deemed to be normal clinically based upon intraoperative inspection under loupe magnification. To avoid unnecessary risk to patients, adjacent fascia was only collected for research use when sufficient quantities of adjacent tissue were removed during the normal course of surgery, additional surgical procedures were not performed to acquire adjacent fascia. Control palmar fascia biopsies were obtained from healthy patients without a history of Dupuytren’s disease undergoing open carpal tunnel release. Specifically, a small section (~1 cm x 1 cm) of palmar fascia located just superficial to the transverse carpal ligament was collected. The specimen was snap frozen in liquid nitrogen, followed by storage at −80 °C until use for molecular analysis. In total, 25 tissue samples were collected and used for high resolution molecular analysis. There were 15 diseased Dupuytren’s fascia biopsies, seven external controls biopsies, and three adjacent fascia specimens. Informed consent was obtained under institutional review board approved protocols for all specimens used in this investigation (Genetic analysis of disorders of fibrosis IRB # 12–000208).
Tissue biopsies were frozen using liquid nitrogen and ground into a powder using a mortar and pestle. Crushed samples were placed into Qiazol reagent and homogenized using the TissueLyser LT (Qiagen, Hilden, Germany). MicroRNAs were extracted from research biopsies using the miRNeasy minikit (Qiagen, Hilden, Germany). Total RNA was quantified using the NanoDrop 2000 spectrophotometer (Thermo Fischer Scientific, Wilmington, Delaware).
RNA Sequencing (RNA-seq) and bioinformatics analysis
MicroRNAs were sequenced using the NEBNext Small RNA library prep kit on an Illumina HiSeq 2000. The short reads were trimmed of adapters with Cutadapt . Trimmed microRNA sequences greater than 17 nucleotides in length were then aligned to the reference genome and miRBase reference sequences using Bowtie . Known microRNA expression and novel microRNA prediction and quantification were performed with miRDeep2 , using the CAP-miRSeq analysis pipeline . Unsupervised hierarchical clustering was performed using the Pearson correlation method. ComiR, a computational tool for combinatorial microRNA target prediction was used to identify molecular pathways regulated by microRNAs that were differentially expressed between diseased and non-diseased palmar fascia [18, 19]. The Database for Annotation and Visualization and Integrated Discovery v6.7 (DAVID 6.7) [20, 21] was used to characterize functional gene clusters regulated by microRNA target genes.
The activity of pro-fibrotic pathways, including the Wnt and TGFβ signaling pathways, were assessed by measuring the expression of regulatory mRNAs using real-time quantitative polymerase chain reaction (RT-qPCR). Total RNA from 11 Dupuytren’s specimens and seven control biopsies, which were previously used for small RNA-seq, were used for cDNA synthesis. Reverse transcription and RT-qPCR reactions were performed as previously described by Dudakovic et al. . Transcript levels were normalized to AKT1, because this gene is most consistently expressed across samples within the mesenchymal lineage compared to other conventional housekeeping genes including GAPDH, HPRT and ACTB based on data we obtained for >400 different mesenchymal cell types and musculoskeletal tissues (SMR & AJvW, unpublished data). Gene expression levels were quantified using the 2-∆∆Ct method. Differences in gene expression between diseased and control samples were evaluated using a two-tailed student’s t-test. Error bars are shown as the mean ± one standard deviation statistical significance was set at p < 0.05 and is indicated by (*). Primer sequences are given in (Additional file 1: Table S1).
MicroRNAs differentially expressed between Dupuytren’s and control palmar fascia
MicroRNAs enriched in Dupuytren’s fascia
MicroRNAs enriched in external control fascia
Differential read count
Differential read count
Dupuytren’s vs external control fascia
External control vs Dupuytren’s fascia
Computational microRNA target prediction
MicroRNAs have the ability to inhibit a large number of target genes by binding to sequence specific 3’UTR regions of target mRNAs, inhibiting their translation and promoting their degradation . To confirm preferential targeting of pathways linked to extracellular matrix synthesis and cellular proliferation, and to identify novel pathways that are regulated by microRNAs in Dupuytren’s disease, we performed computational target prediction using the combinatorial miRNA target prediction tool ComiR [18, 19] (Additional file 3: Table S3). This program uses computational targets generated using miRanda, PITA, TargetScan, and mirSVR, and determines gene targets for a set of microRNAs taking into account the relative expression of each microRNA in a set of samples. We compared microRNA gene targets between Dupuytren’s fascia, and control fascia only, since adjacent fascia may represent an intermediate state between Dupuytren’s tissue and unaffected external control tissue. Abundant, differentially expressed microRNAs with an average expression level of at least 100 normalized reads per million in either the Dupuytren’s or control specimens were evaluated.
RT-qPCR validation of microRNA gene targets
Collagens targeted by microRNAs enriched in control palmar fascia - ComiR target prediction shows that the majority of collagens are preferentially inhibited by microRNAs enriched in control fascia compared with Dupuytren’s fascia. In our analysis of collagen expression using RT-qPCR, COL8A1 and COL15A1 were the only collagens that did not show a statistically significant increase in Dupuytren’s fascia. Notably these collagens also exhibited the lowest degree of microRNA targeting
ComiR score difference (Control Fascia –Dupuytren’s Fascia)
Previous studies using microarrays, which give relative abundance of mRNAs and microRNAs, have been used to study differences in microRNA expression between Dupuytren’s and non-diseased palmar fascia [76, 77]. In comparison with microarrays, RNA-seq provides a broader dynamic range for accurate quantification of differentially expressed transcripts . Therefore in this investigation we applied RNA-seq technology to quantify the expression of all known human microRNAs in Dupuytren’s and control palmar fascia biopsies. Our initial assessment of microRNA sequencing data from the Dupuytren’s and control palmar fascia biopsies using unbiased unsupervised hierarchical clustering confirmed previous findings with microarray analysis, showing that diseased and non-diseased specimens cluster based upon their microRNA profiles [76, 77].
In comparison with our study, Mosakhani et al. applied microarrays to evaluate microRNA expression in Dupuytren’s tissue samples . Our studies confirmed enrichment of microRNAs miR-10b, miR-7f, miR-101, miR-26a, miR-26b, miR-29a, and miR-30 in non-diseased palmar fascia samples. A large number of microRNAs identified as being enriched in either the Dupuytren’s or control fascia were not found to be statistically significant in our analysis. Interestingly, miR-21 the most abundant microRNA in all samples from our study, in contrary to their results, was found to be enriched in Dupuytren’s samples rather than control samples. These differences may be attributable to the fact that microarrays can become saturated with abundantly expressed transcripts, making fold change comparisons unreliable. Additional considerations include the use of fewer control samples (N = 4), and the fact that each of the control samples were taken from the transverse carpal ligament, which is deep to the true palmar fascia which gives rise to the Dupuytren’s cords. Two of the control samples were also collected from patients with acute hand trauma, which could also significantly alter microRNA profiles. Satish et al. compared transcriptomes of fibroblasts derived from Dupuytren’s fascia, palmar fascia, and the transverse carpal ligament using microarrays. They found that fibroblasts from the Dupuytren’s fascia and palmar fascia were more similar to one another than either one was to the transverse carpal ligament derived fibroblasts . Thus the comparisons used in our study comparing palmar fascia to Dupuytren’s fascia are likely to be highly informative.
The microRNAs identified in this investigation that are enriched in control fascia encompass known, as well uncharacterized, but potentially novel anti-fibrotic microRNAs. Established anti-fibrotic microRNAs identified in our analysis include let-7 [23–25], miR-29a-3p , miR-26b-5p, miR-30d-5p [28, 29], miR-27b-3p [30, 31], miR-10a-5p , miR-26a-5p [37–40], miR-101-3p [41–44], miR-27a-3p and miR-10b-5p . Additional enriched microRNAs (miR-126-3p [46–52], miR-99a-5p [53–59], miR-125a-5p [60–64], and miR-139-5p [65–67]) have been shown to affect proliferation in cancer, and may regulate the fibroproliferative activity seen in Dupuytren’s disease. Synergistic activation or inhibition of these microRNAs, which will be investigated in future studies, may collectively permit the potent attenuation or activation of fibrosis for therapeutic applications.
Previous studies have looked at single microRNA targets, and have implicated microRNAs as regulatory factors in Wnt and TGFβ signaling. In this study we used a comprehensive approach simultaneously taking into account the microRNA targets for all abundant, differentially expressed microRNAs. This analysis suggests a loss of microRNAs that target extracellular matrix synthesis in Dupuytren’s disease. MicroRNA target prediction also showed strong correlation with mRNA expression as demonstrated by the expression of extracellular matrix forming collagens. In contrast to previous studies, we did not observe strong evidence for microRNA regulation of Wnt or TGFβ signaling pathways. This finding was also exemplified by RT-qPCR expression data that did not show major difference in the expression of gene transcripts implicated in either of these signaling pathways. These findings do not rule out the possibility of post translational mechanisms (e.g. protein phosphorylation) in propagating either of these pathways. It is also important to note that the tissues examined in this investigation are from patients with advanced disease requiring surgical resection. Collagen synthesis and deposition is characteristic of late stage Dupuytren’s disease where the diseased cords are in the process of consolidating. These studies do not exclude the possibility that earlier stages of disease may still be mediated by over activation of Wnt or TGFβ signaling pathways.
Our analytical approach takes into account the fact that large changes in the expression of a small group of microRNAs may have a more dramatic change on cellular phenotype, than small changes in a large number of different microRNAs because of the lack of coordinated mRNA inhibition by differentially expressed microRNAs. This study also supports the concept that abundantly expressed microRNAs are well suited for fine tune regulation of genes translated from plentiful transcripts such as collagens and other extracellular matrix constituents. Since extracellular matrix targeting microRNAs are constitutively present in large quantities, they can act as a buffer to fine tune extracellular matrix synthesis. This is in stark contrast to alternative regulatory elements such as transcription factors that produce an all or nothing response by directly activating or suppressing mRNA transcription. Analysis of the impact of broad spectrum microRNA targeting, rather than an evaluation of their effects on individual target genes may give greater insight into the important role that microRNAs play in regulating cellular processes.
The main finding of this study is that microRNA profiles show distinct expression patterns that differentiate diseased Dupuytren’s and healthy palmar fascia. The microRNAs enriched in healthy tissue show preferential targeting of collagens and extracellular matrix proteins. This finding is strongly supported by the fact that differential collagen expression as determined by RT-qPCR, is strongly related to the degree of predicted microRNA targeting.
The microRNAs characterized in this investigation have the potential to serve as disease biomarkers that can help guide surgical management by determining optimal surgical margins during open fasciectomy. Novel RNA-therapeutics that are currently in development, also have the potential to target disease specific microRNAs and prophylactically prevent disabling fibrosis minimizing the need for invasive surgical treatments. Fibrosis related microRNAs may also play important regulatory roles in other disorders of fibrosis including scleroderma, idiopathic pulmonary fibrosis, as well as scarring and wound healing.
We thank the members of our research group for stimulating discussions, as well as the Mayo Clinic Bioinformatics Core for their assistance with high-throughput RNA sequencing and bioinformatics support. This work was supported by NIH grants R01 AR049069 (AJvW) and F32 AR066508 (AD), as well as grants from the Mayo Clinic Center for Regenerative Medicine and the Obaid Foundation. We also acknowledge research funding support from the Orthopedic Research and Education Foundation (to SMR), funding from AO North American (to SMR), and a basic science research grant from the American Society for Surgery of the Hand (to SK).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Armstrong JR, Hurren JS, Logan AM. Dermofasciectomy in the management of Dupuytren’s disease. J Bone Joint Surg Br. 2000;82:90–4.View ArticlePubMedGoogle Scholar
- Coert JH, Barret JP, Meek MF. Results of Partial Fasciectomy for Dupuytren Disease in 261 Consecutive Patients. Ann Plast Surg. 2006;57:13–7.View ArticlePubMedGoogle Scholar
- Hurst LC, Badalamente MA, Hentz VR, Hotchkiss RN, Kaplan FT, Meals RA, et al. Injectable collagenase clostridium histolyticum for Dupuytren’s contracture. N Engl J Med. 2009;361:968–79.View ArticlePubMedGoogle Scholar
- Mavrogenis AF, Spyridonos SG, Ignatiadis IA, Antonopoulos D, Papagelopoulos PJ. Partial fasciectomy for Dupuytren’s contractures. J Surg Orthop Adv. 2009;18:106–10.PubMedGoogle Scholar
- Riester SM, van Wijnen AJ, Rizzo M, Kakar S. Pathogenesis and Treatment of Dupuytren Disease. JBJS Reviews. 2014;2:e2. http://dx.doi.org/10.2106/JBJS.RVW.M.00072.View ArticleGoogle Scholar
- Hindocha S, McGrouther DA, Bayat A. Epidemiological evaluation of Dupuytren’s disease incidence and prevalence rates in relation to etiology. Hand (NY). 2009;4:256–69.View ArticleGoogle Scholar
- Dolmans GH, Werker PM, Hennies HC, Furniss D, Festen EA, Franke L, et al. Wnt Signaling in Dupuytren's Disease. N Engl J Med. 2011;365:307–17.View ArticlePubMedGoogle Scholar
- Anderson ER, Burmester JK, Caldwell MD. Evaluation of a mitochondrial DNA mutation in maternally inherited and sporadic cases of Dupuytren disease. Clin Med Res. 2012;10:122–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Hu FZ, Nystrom A, Ahmed A, Palmquist M, Dopico R, Mossberg I, et al. Mapping of an autosomal dominant gene for Dupuytren’s contracture to chromosome 16q in a Swedish family. Clin Genet. 2005;68:424–9.View ArticlePubMedGoogle Scholar
- Zeng C, Wang YL, Xie C, Sang Y, Li TJ, Zhang M, Wang R, Zhang Q, Zheng L, Zhuang SM. Identification of a novel TGF-β-miR-122-fibronectin 1/serum response factor signaling cascade and its implication in hepatic fibrogenesis. Oncotarget. 2015. [Epub ahead of print].Google Scholar
- Gras C, Ratuszny D, Hadamitzky C, Zhang H, Blasczyk R, Figueiredo C. miR-145 contributes to hypertrophic scarring of the skin by inducing myofibroblast activity. Mol Med. 2015. [Epub ahead of print].Google Scholar
- Montgomery RL, Yu G, Latimer PA, Stack C, Robinson K, Dalby CM, et al. MicroRNA mimicry blocks pulmonary fibrosis. EMBO Mol Med. 2014;6:1347–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Lian JB, Stein GS, van Wijnen AJ, Stein JL, Hassan MQ, Gaur T, et al. MicroRNA control of bone formation and homeostasis. Nat Rev Endocrinol. 2012;8:212–27.View ArticlePubMedPubMed CentralGoogle Scholar
- M. Martin. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal, North America, 17, May 2011. Available at: http://journal.embnet.org/index.php/embnetjournal/article/view/200.
- Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.View ArticlePubMedPubMed CentralGoogle Scholar
- Friedländer MR, Chen W, Adamidi C, Maaskola J, Einspanier R, Knespel S, et al. Discovering microRNAs from deep sequencing data using miRDeep. Nat Biotechnol. 2008;26:407–15.View ArticlePubMedGoogle Scholar
- Sun Z, Evans J, Bhagwate A, Middha S, Bockol M, Yan H, et al. CAP-miRSeq: a comprehensive analysis pipeline for microRNA sequencing data. BMC Genomics. 2014;15:423.View ArticlePubMedPubMed CentralGoogle Scholar
- Coronnello C, Hartmaier R, Arora A, Huleihel L, Pandit KV, Bais AS, et al. Novel modeling of combinatorial miRNA targeting identifies SNP with potential role in bone density. PLoS Comput Biol. 2012;8, e1002830.View ArticlePubMedPubMed CentralGoogle Scholar
- Coronnello C, Benos PV. ComiR: Combinatorial microRNA target prediction tool. Nucleic Acids Res. 2013;41(Web Server issue):W159–64.View ArticlePubMedPubMed CentralGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.View ArticleGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13.View ArticleGoogle Scholar
- Dudakovic A, Camilleri E, Riester SM, Lewallen EA, Kvasha S, Chen X, et al. High-resolution molecular validation of self-renewal and spontaneous differentiation in clinical-grade adipose-tissue derived human mesenchymal stem cells. J Cell Biochem. 2014;115:1816–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Park JT, Kato M, Lanting L, Castro N, Nam BY, Wang M, et al. Repression of let-7 by transforming growth factor-β1-induced Lin28 upregulates collagen expression in glomerular mesangial cells under diabetic conditions. Am J Physiol Renal Physiol. 2014;307:F1390–403.View ArticlePubMedPubMed CentralGoogle Scholar
- Makino K, Jinnin M, Hirano A, Yamane K, Eto M, Kusano T, et al. The downregulation of microRNA let-7a contributes to the excessive expression of type I collagen in systemic and localized scleroderma. J Immunol. 2013;190:3905–15.View ArticlePubMedGoogle Scholar
- Pandit KV, Corcoran D, Yousef H, Yarlagadda M, Tzouvelekis A, Gibson KF, et al. Inhibition and role of let-7d in idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 2010;182:220–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Li Z, Hassan MQ, Jafferji M, Aqeilan RI, Garzon R, Croce CM, et al. Biological functions of miR-29b contribute to positive regulation of osteoblast differentiation. J Biol Chem. 2009;284:15676–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Tu X, Zheng X, Li H, Cao Z, Chang H, Luan S, Zhu J, Chen J, Zang Y, Zhang J. MicroRNA-30 Protects Against Carbon Tetrachloride-induced Liver Fibrosis by Attenuating Transforming Growth Factor Beta Signaling in Hepatic Stellate Cells. Toxicol Sci. 2015. [Epub ahead of print].Google Scholar
- Zhou Q, Yang M, Lan H, Yu X. miR-30a negatively regulates TGF-β1-induced epithelial-mesenchymal transition and peritoneal fibrosis by targeting Snai1. Am J Pathol. 2013;183:808–19.View ArticlePubMedGoogle Scholar
- Graham JR, Williams CM, Yang Z. MicroRNA-27b targets gremlin 1 to modulate fibrotic responses in pulmonary cells. J Cell Biochem. 2014;115:1539–48.View ArticlePubMedGoogle Scholar
- Kwan P, Ding J, Tredget EE. MicroRNA 181b Regulates Decorin Production by Dermal Fibroblasts and May Be a Potential Therapy for Hypertrophic Scar. PloS One. 2015;10, e0123054.View ArticlePubMedPubMed CentralGoogle Scholar
- Li C, Zhu HY, Bai WD, Su LL, Liu JQ, Cai WX, et al. MiR-10a and miR-181c regulate collagen type I generation in hypertrophic scars by targeting PAI-1 and uPA. FEBS Lett. 2015;589:380–9.View ArticlePubMedGoogle Scholar
- Li X, Liu L, Shen Y, Wang T, Chen L, Xu D, et al. MicroRNA-26a modulates transforming growth factor beta-1-induced proliferation in human fetal lung fibroblasts. Biochem Biophys Res Commun. 2014;454:512–7 [Epub ahead of print].View ArticlePubMedGoogle Scholar
- Liang H, Gu Y, Li T, Zhang Y, Huangfu L, Hu M, et al. Integrated analyses identify the involvement of microRNA-26a in epithelial-mesenchymal transition during idiopathic pulmonary fibrosis. Cell Death Dis. 2014;5, e1238.View ArticlePubMedPubMed CentralGoogle Scholar
- Liang H, Xu C, Pan Z, Zhang Y, Xu Z, Chen Y, et al. The antifibrotic effects and mechanisms of microRNA-26a action in idiopathic pulmonary fibrosis. Mol Ther. 2014;22:1122–33.PubMedPubMed CentralGoogle Scholar
- Wei C, Kim IK, Kumar S, Jayasinghe S, Hong N, Castoldi G, et al. NF-κB mediated miR-26a regulation in cardiac fibrosis. J Cell Physiol. 2013;228:1433–42.View ArticlePubMedGoogle Scholar
- Zhao S, Zhang Y, Zheng X, Tu X, Li H, Chen J, Zang Y, Zhang J. Loss of microRNA-101 Promotes Epithelial to Mesenchymal Transition in Hepatocytes. J Cell Physiol. 2015. [Epub ahead of print].Google Scholar
- Zhao X, Wang K, Liao Y, Zeng Q, Li Y, Hu F, et al. MicroRNA-101a inhibits cardiac fibrosis induced by hypoxia via targeting TGFβRI on cardiac fibroblasts. Cell Physiol Biochem. 2015;35:213–26.View ArticlePubMedGoogle Scholar
- Tu X, Zhang H, Zhang J, Zhao S, Zheng X, Zhang Z, et al. MicroRNA-101 suppresses liver fibrosis by targeting the TGFβ signalling pathway. J Pathol. 2014;234:46–59.View ArticlePubMedGoogle Scholar
- Pan Z, Sun X, Shan H, Wang N, Wang J, Ren J, et al. MicroRNA-101 inhibited postinfarct cardiac fibrosis and improved left ventricular compliance via the FBJ osteosarcoma oncogene/transforming growth factor-β1 pathway. Circulation. 2012;126:840–50.View ArticlePubMedGoogle Scholar
- Han X, Yan S, Weijie Z, Feng W, Liuxing W, Mengquan L, et al. Critical role of miR-10b in transforming growth factor-β1-induced epithelial-mesenchymal transition in breast cancer. Cancer Gene Ther. 2014;21:60–7.View ArticlePubMedGoogle Scholar
- Zhao C, Li Y, Zhang M, Yang Y, Chang L. miR-126 inhibits cell proliferation and induces cell apoptosis of hepatocellular carcinoma cells partially by targeting Sox2. Hum Cell. 2015;28:91–9.View ArticlePubMedGoogle Scholar
- Jiang L, He A, Zhang Q, Tao C. miR-126 inhibits cell growth, invasion, and migration of osteosarcoma cells by downregulating ADAM-9. Tumour Biol. 2014;35:12645–54.View ArticlePubMedGoogle Scholar
- Yang C, Hou C, Zhang H, Wang D, Ma Y, Zhang Y, et al. miR-126 functions as a tumor suppressor in osteosarcoma by targeting Sox2. Int J Mol Sci. 2013;15:423–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu Q, Liu SL, Wang H, Shi G, Yang P, Chen XL. miR-126 Suppresses the proliferation of cervical cancer cells and alters cell sensitivity to the chemotherapeutic drug bleomycin. Asian Pac J Cancer Prev. 2014;14:6569–72.View ArticlePubMedGoogle Scholar
- Zhou Y, Feng X, Liu YL, Ye SC, Wang H, Tan WK, et al. Down-regulation of miR-126 is associated with colorectal cancer cells proliferation, migration and invasion by targeting IRS-1 via the AKT and ERK1/2 signaling pathways. PLoS One. 2013;8, e81203.View ArticlePubMedPubMed CentralGoogle Scholar
- Guo CJ, Pan Q, Xiong H, Qiao YQ, Bian ZL, Zhong W, et al. Dynamic expression of miR-126* and its effects on proliferation and contraction of hepatic stellate cells. FEBS Lett. 2013;587:3792–801.View ArticlePubMedGoogle Scholar
- Li N, Tang A, Huang S, Li Z, Li X, Shen S, et al. MiR-126 suppresses colon cancer cell proliferation and invasion via inhibiting RhoA/ROCK signaling pathway. Mol Cell Biochem. 2013;380:107–19.View ArticlePubMedGoogle Scholar
- Feng Y, Kang Y, He Y, Liu J, Liang B, Yang P, et al. microRNA-99a acts as a tumor suppressor and is down-regulated in bladder cancer. BMC Urol. 2014;14:50.View ArticlePubMedPubMed CentralGoogle Scholar
- Song Y, Dou H, Wang P, Zhao S, Wang T, Gong W, et al. A novel small-molecule compound diaporine A inhibits non-small cell lung cancer growth by regulating miR-99a/mTOR signaling. Cancer Biol Ther. 2014;15:1423–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu D, Zhou Y, Pan H, Zhou J, Fan Y, Qu P. microRNA-99a inhibiting cell proliferation, migration and invasion by targeting fibroblast growth factor receptor 3 in bladder cancer. Oncol Lett. 2014;7:1219–24. Epub 2014 Feb 12.PubMedPubMed CentralGoogle Scholar
- Zhang EB, Kong R, Yin DD, You LH, Sun M, Han L, et al. Long noncoding RNA ANRIL indicates a poor prognosis of gastric cancer and promotes tumor growth by epigenetically silencing of miR-99a/miR-449a. Oncotarget. 2014;5:2276–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang L, Chang L, Li Z, Gao Q, Cai D, Tian Y, et al. miR-99a and -99b inhibit cervical cancer cell proliferation and invasion by targeting mTOR signaling pathway. Med Oncol. 2014;31:934.View ArticlePubMedGoogle Scholar
- Jin Y, Tymen SD, Chen D, Fang ZJ, Zhao Y, Dragas D, et al. MicroRNA-99 family targets AKT/mTOR signaling pathway in dermal wound healing. PLoS One. 2013;8, e64434.View ArticlePubMedPubMed CentralGoogle Scholar
- Cui L, Zhou H, Zhao H, Zhou Y, Xu R, Xu X, et al. MicroRNA-99a induces G1-phase cell cycle arrest and suppresses tumorigenicity in renal cell carcinoma. BMC Cancer. 2012;12:546.View ArticlePubMedPubMed CentralGoogle Scholar
- Tiwari A, Shivananda S, Gopinath KS, Kumar A. MicroRNA-125a reduces proliferation and invasion of oral squamous cell carcinoma cells by targeting estrogen-related receptor α: implications for cancer therapeutics. J Biol Chem. 2014;289:32276–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Ufkin ML, Peterson S, Yang X, Driscoll H, Duarte C, Sathyanarayana P. miR-125a regulates cell cycle, proliferation, and apoptosis by targeting the ErbB pathway in acute myeloid leukemia. Leuk Res. 2014;38:402–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Ninio-Many L, Grossman H, Shomron N, Chuderland D, Shalgi R. microRNA-125a-3p reduces cell proliferation and migration by targeting Fyn. J Cell Sci. 2013;126:2867–76.View ArticlePubMedGoogle Scholar
- Bi Q, Tang S, Xia L, Du R, Fan R, Gao L, et al. Ectopic expression of MiR-125a inhibits the proliferation and metastasis of hepatocellular carcinoma by targeting MMP11 and VEGF. PLoS One. 2012;7, e40169.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang Y, Gao JS, Tang X, Tucker LD, Quesenberry P, Rigoutsos I, et al. MicroRNA 125a and its regulation of the p53 tumor suppressor gene. FEBS Lett. 2009;583:3725–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang L, Dong Y, Zhu N, Tsoi H, Zhao Z, Wu CW, et al. microRNA-139-5p exerts tumor suppressor function by targeting NOTCH1 in colorectal cancer. Mol Cancer. 2014;13:124.View ArticlePubMedPubMed CentralGoogle Scholar
- Luo HN, Wang ZH, Sheng Y, Zhang Q, Yan J, Hou J, et al. MiR-139 targets CXCR4 and inhibits the proliferation and metastasis of laryngeal squamous carcinoma cells. Med Oncol. 2014;31:789.View ArticlePubMedGoogle Scholar
- Liu R, Yang M, Meng Y, Liao J, Sheng J, Pu Y, et al. Tumor-suppressive function of miR-139-5p in esophageal squamous cell carcinoma. PLoS One. 2013;8, e77068.View ArticlePubMedPubMed CentralGoogle Scholar
- Shi C, Yang Y, Xia Y, Okugawa Y, Yang J, Liang Y, Chen H, Zhang P, Wang F, Han H, Wu W, Gao R, Gasche C, Qin H, Ma Y, Goel A. Novel evidence for an oncogenic role of microRNA-21 in colitis-associated colorectal cancer. Gut. 2015. pii:gutjnl-2014-308455. doi:https://doi.org/10.1136/gutjnl-2014-308455.
- Ma X, Conklin DJ, Li F, Dai Z, Hua X, Li Y, et al. The oncogenic microRNA miR-21 promotes regulated necrosis in mice. Nat Commun. 2015;6:7151.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang CH, Pfeffer SR, Sims M, Yue J, Wang Y, Linga VG, et al. The oncogenic microRNA-21 inhibits the tumor suppressive activity of FBXO11 to promote tumorigenesis. J Biol Chem. 2015;290:6037–46.View ArticlePubMedPubMed CentralGoogle Scholar
- Zuo J, Wen M, Lei M, Peng X, Yang X, Liu Z. MiR-210 Links Hypoxia With Cell Proliferation Regulation in Human Laryngocarcinoma Cancer. J Cell Biochem. 2015;116:1039–49.View ArticlePubMedGoogle Scholar
- Bodempudi V, Hergert P, Smith K, Xia H, Herrera J, Peterson M, et al. miR-210 promotes IPF fibroblast proliferation in response to hypoxia. Am J Physiol Lung Cell Mol Physiol. 2014;307:L283–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang C, Song X, Li Y, Han F, Gao S, Wang X, et al. Low-dose paclitaxel ameliorates pulmonary fibrosis by suppressing TGF-β1/Smad3 pathway via miR-140 upregulation. PLoS One. 2013;8, e70725.View ArticlePubMedPubMed CentralGoogle Scholar
- Ji J, Zhang J, Huang G, Qian J, Wang X, Mei S. Over-expressed microRNA-27a and 27b influence fat accumulation and cell proliferation during rat hepatic stellate cell activation. FEBS Lett. 2009;583:759–66.View ArticlePubMedGoogle Scholar
- Zheng J, Wu C, Xu Z, Xia P, Dong P, Chen B, et al. Hepatic stellate cell is activated by microRNA-181b via PTEN/Akt pathway. Mol Cell Biochem. 2015;398:1–9.View ArticlePubMedGoogle Scholar
- Brockhausen J, Tay SS, Grzelak CA, Bertolino P, Bowen DG, d’Avigdor WM, et al. miR-181a mediates TGF-β-induced hepatocyte EMT and is dysregulated in cirrhosis and hepatocellular cancer. Liver Int. 2015;35:240–53.View ArticlePubMedGoogle Scholar
- Chen C, Wu CQ, Zhang ZQ, Yao DK, Zhu L. Loss of expression of miR-335 is implicated in hepatic stellate cell migration and activation. Exp Cell Res. 2011;317:1714–25.View ArticlePubMedGoogle Scholar
- Ameres SL, Martinez J, Schroeder R. Molecular basis for target RNA recognition and cleavage by human RISC. Cell. 2007;130:101–12.View ArticlePubMedGoogle Scholar
- Krause C, Kloen P, Ten Dijke P. Elevated transforming growth factor β and mitogen-activated protein kinase pathways mediate fibrotic traits of Dupuytren’s disease fibroblasts. Fibrogenesis Tissue Repair. 2011;4:14.View ArticlePubMedPubMed CentralGoogle Scholar
- Kloen P, Jennings CL, Gebhardt MC, Springfield DS, Mankin HJ. Transforming growth factor-beta: possible roles in Dupuytren’s contracture. J Hand Surg Am. 1995;20:101–8.View ArticlePubMedGoogle Scholar
- Satish L, LaFramboise WA, Johnson S, Vi L, Njarlangattil A, Raykha C, et al. Fibroblasts from phenotypically normal palmar fascia exhibit molecular profiles highly similar to fibroblasts from active disease in Dupuytren’s Contracture. BMC Med Genomics. 2012;5:15.View ArticlePubMedPubMed CentralGoogle Scholar
- Mosakhani N, Guled M, Lahti L, Borze I, Forsman M, Pääkkönen V, et al. Unique microRNA profile in Dupuytren’s contracture supports deregulation of β-catenin pathway. Mod Pathol. 2010;23:1544–52.View ArticlePubMedGoogle Scholar
- Liu S, Lin L, Jiang P, Wang D, Xing Y. A comparison of RNA-Seq and high-density exon array for detecting differential gene expression between closely related species. Nucleic Acids Res. 2011;39:578–88.View ArticlePubMedGoogle Scholar