- Research article
- Open access
- Published:
RNA sequencing reveals a depletion of collagen targeting microRNAs in Dupuytren’s disease
BMC Medical Genomics volume 8, Article number: 59 (2015)
Abstract
Background
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.
Methods
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.
Results
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.
Discussion
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.
Conclusion
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.
Background
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 [6]. 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.
Methods
Tissue collection
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).
RNA extraction
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 [14]. Trimmed microRNA sequences greater than 17 nucleotides in length were then aligned to the reference genome and miRBase reference sequences using Bowtie [15]. Known microRNA expression and novel microRNA prediction and quantification were performed with miRDeep2 [16], using the CAP-miRSeq analysis pipeline [17]. 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.
RT-qPCR validation
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. [22]. 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).
Results
RNA-seq evaluation
In this investigation we applied high throughput next generation RNA-seq to quantify the expression of all known human microRNAs (N = 2252) to understand their roles in the pathogenesis of Dupuytren’s disease. Initial sequencing results showed that more than 99 % of the sequencing reads averaged across all Dupuytren’s and control samples were represented by the top 100 most abundantly expressed microRNAs (Fig. 1). This distribution is consistent across both the diseased and control tissues, and shows that a relatively small subset of all known human microRNAs account for the majority of microRNAs expressed in palmar fascia biopsies. These abundant microRNAs may therefore play a more important functional role in regulating cellular processes connected to Dupuytren’s disease pathogenesis than microRNAs with low levels of expression.
Unsupervised hierarchical clustering was performed using the Pearson correlation method to provide an unbiased assessment of the ability of microRNAs to differentiate between diseased and non-diseased palmar fascia biopsies (Fig. 2). The clustering dendrogram showed independent grouping of the control and diseased palmar fascia biopsies. There was only a single adjacent fascia biopsy that clustered with the Dupuytren’s specimens suggesting that this specimen may have in fact been diseased tissue that appeared grossly normal intra-operatively.
To identify specific microRNAs with a functional role in the pathogenesis of Dupuytren’s disease we made a fold change comparison of microRNA profiles in Dupuytren’s biopsies against control samples (excluding adjacent fascia). In the comparison of Dupuytren’s fascia and control fascia we identified 74 microRNAs with a statistically significant 2-fold enrichment in Dupuytren’s tissue, and 32 microRNAs with a statistically significant enrichment in control fascia (Additional file 2: Table S2). The diseased Dupuytren’s fascia show enrichment (fold change > 2) in a greater variety of different microRNAs compared with control fascia (Fig. 3a), and an even larger increase in the overall number of microRNA sequencing reads (Fig. 3b). To determine if the differentially expressed sequencing reads were being generated by a large number of different microRNAs, or a focused group of microRNAs we examined the distribution of differentially expressed microRNA sequencing reads. The majority of the differentially expressed microRNA reads were found to be concentrated amongst a relatively small group of differentially expressed microRNAs (Fig. 4). Notably the microRNAs upregulated in the control fascia accounting for the greatest differential in read count are heavily enriched in previously validated anti-fibrotic extracellular matrix targeting microRNAs (Table 1), including let-7 [23–25], miR-29a-3p [26], miR-26b-5p, miR-30d-5p [27, 28], miR-27a-3p, miR-27b-3p [29, 30], miR-10a-5p [31], miR-26a-5p [32–35], miR-101-3p [36–39], and miR-10b-5p [40], as well as anti-proliferative microRNAs including, miR-126-3p [41–47], miR-99a-5p [48–54], miR-125a-5p [55–59], and miR-139-5p [60–62]. MicroRNAs enriched in Dupuytren’s fascia target appear to target a diverse array of functional pathways. miR-21-5p is an established oncomiR that promotes cellular proliferation in cancer [63–65]. miR-210 is a known hypoxia inducible, fibrosis promoting microRNA that may be associated with decreased tissue perfusion in Dupuytren’s cords [66, 67]. These findings suggest that a relatively small group of microRNAs with large differences in expression act to selectively target molecular pathways regulating tissue fibrosis.
To determine if adjacent fascia exhibits a disease-like molecular phenotype, we compared the adjacent fascia biopsies with patient matched diseased Dupuytren’s fascia. A direct comparison of the three matched adjacent fascia and Dupuytren’s fascia biopsies revealed miR-181b-5p as the only microRNA with a statistically significant enrichment greater than 1.4 fold in the Dupuytren’s fascia. The lack of statistically significant microRNAs is attributable largely to the variability in the adjacent fascia biopsy specimens (Dupuytren vs control phenotype); therefore we compared upregulated microRNAs in each paired sample using a Venn diagram method (Fig. 5). This analysis showed an enrichment of anti-fibrotic extracellular matrix targeting microRNAs in two of the adjacent fascia specimens relative to their matched Dupuytren’s biopsy [23–29, 68, 69]. However the single adjacent fascia specimen that clustered with the diseased Dupuytren’s fascia during unsupervised hierarchical clustering (specimen “A”) did not show enrichment in collagen/extracellular matrix targeting microRNAs. These findings support its preferential clustering with the Dupuytren’s fascia and indicate that its molecular phenotype more closely resembles that of diseased fascia rather than control fascia despite its grossly normal appearance. A comparison of adjacent and control fascial biopsies also did not reveal any microRNAs with a statistically significant enrichment in the adjacent fascia, again highlighting the variability of the adjacent fascia biopsies. There were however seven microRNAs (miR-335-3p, miR-128, miR-224-5p, miR-28-5p, miR-191-5p, miR-423-3p, and miR-181a-5p) with a statistically significant upregulation in the controls. Two of these seven microRNAs, miR-335-3p, and miR-181a-5p are also validated extracellular matrix targeting microRNAs [30, 31, 70–72]. Their enrichment in external control specimens supports their previously described role as anti-fibrotic microRNAs, and is a reflection of a partial disease phenotype in the adjacent fascia specimens. These findings highlight the potential utility of microRNAs as biomarkers to confirm diagnosis and define optimal surgical margins to ensure all diseased fascia is removed at the time of surgery. It also implicates microRNAs as potential therapeutic targets/agents to prevent as well as treat active fibrotic disease.
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 [73]. 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.
Analysis with ComiR identified 6685 genes that were computationally predicted to be inhibited by microRNAs enriched in control fascia (p < 0.05). In the Dupuytren’s fascia there were 4969 genes computationally predicted to be targeted by enriched microRNAs. Functional annotation clustering using David 6.7 for the top 3000 most heavily targeted genes in control and diseased fascia identified candidate molecular pathways regulated by microRNAs (Fig. 6). ComiR analysis confirmed preferential targeting of extracellular matrix synthesis and cellular proliferation in the control palmar fascia. Notable microRNA inhibited pathways in control fascia included wound healing (enrichment score 1.99), extracellular matrix synthesis (enrichment score 1.48), cytokine production (enrichment score 1.48), mitochondrion (enrichment score 1.39), muscle cell differentiation (enrichment score 1.37), cell cycle (enrichment score 1.3), epithelial cell differentiation and keratinization (enrichment score 1.29), and collagen and laminin production (enrichment score 1.16). In diseased Dupuytren’s fascia, microRNAs are most notably predicted to inhibit pathways linked to lipid-binding/lipoproteins (enrichment score 2.16), cadherins (enrichment score 1.97), fatty acid biosynthesis (enrichment score 1.61), and metalloenzyme inhibitors (enrichment score 1.16).
RT-qPCR validation of microRNA gene targets
Since extracellular matrix synthesis is predicted to be inhibited by microRNAs enriched in control palmar fascia, we examined the expression of collagens, the primary constituents of extracellular matrix in palmar fascia. We performed RT-qPCR to evaluate the expression of selected collagens linked to fibrosis in diseased and non-diseased tissues (Fig. 7). We then examined which of these collagens were preferentially targeted for microRNA inhibition in control palmar fascia based on ComiR analysis (Table 2). RT-qPCR of fibrosis related collagens shows that collagens predicted to be targeted by microRNA inhibition in controls are enriched in Dupuytren’s fascia, while collagens with less microRNA targeting (e.g. COL15A1 and COL18A1) do not show statistically significant differences in expression between diseased and non-diseased tissue. These results show that the expression of collagen targeting microRNAs inversely correlates with the expression of fibrosis related collagens in palmar fascia.
Previous studies have implicated alterations in Wnt and TGFβ signaling in the pathogenesis of Dupuytren’s disease [7, 74, 75]. To determine if microRNAs regulate either of these pathways, we evaluated the expression of mRNAs associated with these two signaling pathways using RT-qPCR (Fig. 8). TGFBR2 and WNT5A were the only two genes to show statistically significant differences in expression between the Dupuytren’s and control palmar fascia. TGFBR2 paradoxically showed increased expression in the control fascia, enrichment inconsistent with its known fibrosis promoting affects. WNT5A was the only gene to show a statistically significant enrichment in the Dupuytren’s fascia. However microRNA target prediction indicates that WNT5A is preferentially inhibited by microRNAs enriched in Dupuytren’s fascia, indicating that microRNA regulation is unlikely to be driving the increased expression of WNT5A in diseased palmar fascia. These findings suggest that microRNAs, which promote mRNA degradation, do not play a major role in post-transcriptional regulation of WNT and TGFβ signaling pathways in the terminal stages of Dupuytren’s disease.
Discussion
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 [78]. 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 [77]. 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 [76]. 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 [26], miR-26b-5p, miR-30d-5p [28, 29], miR-27b-3p [30, 31], miR-10a-5p [33], miR-26a-5p [37–40], miR-101-3p [41–44], miR-27a-3p and miR-10b-5p [45]. 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.
Conclusions
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.
References
Armstrong JR, Hurren JS, Logan AM. Dermofasciectomy in the management of Dupuytren’s disease. J Bone Joint Surg Br. 2000;82:90–4.
Coert JH, Barret JP, Meek MF. Results of Partial Fasciectomy for Dupuytren Disease in 261 Consecutive Patients. Ann Plast Surg. 2006;57:13–7.
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.
Mavrogenis AF, Spyridonos SG, Ignatiadis IA, Antonopoulos D, Papagelopoulos PJ. Partial fasciectomy for Dupuytren’s contractures. J Surg Orthop Adv. 2009;18:106–10.
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.
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.
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.
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.
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.
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].
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].
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.
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.
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.
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.
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.
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.
Coronnello C, Benos PV. ComiR: Combinatorial microRNA target prediction tool. Nucleic Acids Res. 2013;41(Web Server issue):W159–64.
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.
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.
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.
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.
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.
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.
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.
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].
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.
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.
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.
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.
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].
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.
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.
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.
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].
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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: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.
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.
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.
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.
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.
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.
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.
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.
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.
Ameres SL, Martinez J, Schroeder R. Molecular basis for target RNA recognition and cleavage by human RISC. Cell. 2007;130:101–12.
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.
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.
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.
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.
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.
Acknowledgements
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).
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare no competing interests related to this publication.
Authors’ contributions
SR, DA, AJvW, and SK participated in the design of the study. SR, DA, MR and SK recruited patients into the study and collected research specimens. JT helped classify and select specimens for research analysis. SR and AD prepared samples for microRNA sequencing. EC, AD, AvW and CP carried out RT-qPCR analysis of tissue specimens. JE performed the primary bioinformatics analysis of sequencing data including sequence read alignment. PK and MKJ provided expertise on experimental strategies for analyzing tissue samples. SR, AvW and SK drafted the manuscript. All authors read and approved the final manuscript.
Authors' information
Not applicable.
Additional files
Additional file 1: Table S1.
Primer sequences used for RT-qPCR evaluation of gene expression in this investigation. (XLSX 10 kb)
Additional file 2: Table S2.
MicroRNAs showing a 2-fold, or greater statistically significant differences in expression between Dupuytren’s and control palmar fascia. (XLSX 13 kb)
Additional file 3: Table S3.
ComiR scores for microRNAs showing the greatest difference in expression based on total read count. (XLSX 4350 kb)
Rights and permissions
Open Access This 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.
About this article
Cite this article
Riester, S.M., Arsoy, D., Camilleri, E.T. et al. RNA sequencing reveals a depletion of collagen targeting microRNAs in Dupuytren’s disease. BMC Med Genomics 8, 59 (2015). https://doi.org/10.1186/s12920-015-0135-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12920-015-0135-8