Skip to main content

Altered microRNAs in C3H10T1/2 cells induced by p.E95K mutant IHH signaling



Indian Hedgehog (IHH), an important cell signaling protein, plays a key regulatory role in development of cartilage and chondrogenesis. Earlier studies have shown that heterozygous missense mutations in IHH gene may cause brachydactyly type A1 (BDA1), an autosomal dominant inheritance disease characterized by apparent shortness or absence of the middle phalanges of all digits. MicroRNAs (miRNAs) have been found to be significant post-transcriptional regulators of gene expression and significantly influence the process of bone-development. Therefore, it is possible that miRNAs are involved in the mechanism underlying the development of BDA1. However, the relationship between miRNAs and the pathogenesis of BDA1 remains unclear.


In this study, we used microarray-based miRNA profiling to investigate the role of miRNAs in BDA1 by characterization of differentially expressed miRNAs in C3H10T1/2 cell line induced by wild type (WT) and p.E95K mutant (MT) IHH signaling.


Our results identified 6 differentially expressed miRNAs between WT and control (CT) group and 5 differentially expressed miRNAs between MT and CT groups. In particular, miR-135a-1-3p was found to be a significantly differentially expressed miRNA between WT and CT group. Results of dual-luciferase reporter gene experiment successfully discovered Hoxd10 was one of the target gene of miR-135a-1-3p. Additionally, our pathway analysis revealed that the targets of these miRNAs of interest were highly involved with Runx1/2, Notch and collagen-related pathways.


Taken together, our findings provided important clue for future study of the process of miRNA-regulation in IHH signaling and novel insights into the regulatory role of miRNA in pathogenesis of BDA1.


Hedgehog (HH) signaling pathway is known to be a key regulator in development of cartilage in all the animals that possess bilateral symmetry [1, 2]. In mammals, there are three homologous proteins that are actively involved in the HH signaling pathway: sonic hedgehog (SHH), desert hedgehog (DHH) and Indian Hedgehog (IHH). In particular, IHH signaling pathway is an essential developmental pathway and plays a critical regulatory role in many physiological processes [3]. Members of this pathway are mainly expressed by chondrocytes in mammal development of cartilage [4]. IHH protein interacts with patched 1 (PTCH1), a transmembrane receptor protein of its receptor. The interaction between IHH and PTCH1 results in the suppression of smoothened (SMO). The suppression of SMO allows transducing activation signal to the downstream glioma-associated oncogene homolog (Gli) transcription factors [5, 6]. IHH/Gli signaling pathway induces differentiation of osteoblasts and plays an important role in endochondral osteogenesis [7]. Any abnormalities in IHH/Gli signaling pathway leads to abnormal bone development.

MicroRNAs (miRNAs) are single-stranded non-coding RNAs (17-25nt) that are known to play a key role in various cellular pathways, including pathways of growth and proliferation, apoptosis, and developmental timing [8,9,10]. MiRNAs are significant post-transcriptional regulators of gene expression that have been founded in a great number of organisms [11, 12]. MiRNAs negatively regulate specific messager RNA (mRNA) translation by direct binding to partially complementary sequences in target mRNA. Previous studies have shown that miRNAs negatively regulate osteoblast differentiation and bone formation by targeting related signaling pathways or downstream transcription factors [13, 14]. For example, miRNA-467g has been found to inhibit new bone regeneration through suppression of HH signaling pathway by targeting Runx-2[13]. However, the mechanism of miRNA-mediated regulation of abnormalities in IHH pathway remains unclear.

Brachydactyly type A1 (BDA1) is an autosomal dominant inheritance disease caused by missense mutations of heterozygotes in IHH. It is a dysplasia/aplastic disorder characterized by the shortening or missing of the middle phalanges [15, 16]. A number of earlier studies have suggested that p.E95K in the N-terminal fragment of IHH protein, a mutation found in BDA1 individuals, affected the interaction between PTCH1 and its partners. This interaction of PTCH1 may be able to prevent the induction of the signaling [17]. The p.E95K mutation has been found to affect the differentiation of chondrocyte and severely impair the development process of cartilage, resulting in a delayed formation of endochondral bone and brachydactyly phenotype [6]. It has also been found that the p.E95K mutant IHH protein affected Gli-mediated downstream regulation through a series of signal transduction [18]. However, the modulation of miRNAs in p.E95K mutant IHH signaling remains unknown.

In this study, we performed a microarray-based screening of miRNA expression in response to wild type (WT) or p.E95K mutant (MT) IHH signaling to identify the differentially expressed miRNAs. Our results will help to explain the mediation role of miRNA in the regulation in IHH signaling and provide novel insight into the molecular mechanisms underlying p.E95K IHH induced BDA1.

Materials and methods

IHH protein expression and purification

The pGEX-2T-based recombinant human IHH-N-terminus proteins (IhhN), including wild type protein (WT, amino acid residues 28-202) and a mutant protein (WT, p.E95k), were expressed by the BL21 (DE3) strain of Escherichia coli (E. coli). We obtained the two transformed strains from Bio-X Institutes, Shanghai Jiao Tong University, which had been described in earlier publication [19]. After coated plates, monoclonal colony was selected and incubated in 5ml LB liquid medium with ampicillin at 37℃ overnight. Then the medium was transferred to 400 ml 2×YT liquid medium with ampicillin and incubated at 37℃ for 4 h. When the value of OD600 reached 0.6, the strains were induced with 0.2 mM IPTG at 32℃ for 5 h. Pellets were collected by centrifugation at 8 000× g for 10 min at 4℃ and re-suspended in 200 mM NaCl and 20 mM Tris-HCl, pH 7.5 with cocktail protein inhibitors (Sigma-Aldrich), followed by sonication. After centrifugation at 25,000× g for 30 min at 4℃, the supernatant was collected. Then the proteins were purified according to the Bulk and RediPack GST Purification Modules (GE-Healthcare), followed by cleavage of the GST tag using thrombin. The purified IhhN proteins were verified by western blotting with IHH Antibody (SC-1782, Santa Cruz) and filtered by 0.22 μm filter (Millipore).

IHH activation of C3H10T1/2 cells

The induction was performed using the method that had been described in previous study [19]. C3H10T1/2 cells (obtained from ATCC) were chosen to analyze IHH signaling in the process of osteoblast differentiation of the mesenchymal cell line [20]. C3H10T1/2 cells were cultured in 12-well plates in a growth medium containing MEM/EBSS, 10 % fetal bovine serum (FBS), and 2 % penicillin/streptomycin. WT or MT IhhN (E95K) was added into the growth mediums (1x = 750 nM IhhN protein finally) when cells grew to 5 × 106 per well. After 48 h of further incubation at 37℃ (5 % CO2), the cells were harvested for further analysis. Expression of Gli1, PTCH1 and IHH, as markers of IHH signaling pathway, were detected by real time PCR. All the induction assays were performed in triplicate.

Microarray-based miRNA Expression Analysis

We used TRIZOL™ solution (Invitrogen) to extract the total RNA from cell samples according to the manufacturer’s instruction. RNAs were examined by Agilent 2100 Bioanalyzer to assure its RIN to be more than 7. MicroRNA microarray analyses were performed using 200 ng total RNA on Unrestricted Mouse miRNA Microarray v2 (Agilent Technologies). The amplification, labeling, hybridization and scanning of samples were performed by Agilent Technologies (Shanghai, China) using standard protocols. The data was extracted by Feature Extraction version 12.5 (Agilent Technologies) and analyzed by Genespring GX version 12.6 (Agilent Technologies). The raw data can be obtained from GEO (GSE74023).

Bioinformatics Analysis

Target genes regulated by differentially expressed miRNAs (DEMs) were predicted using the TargetScan database [21] (Version 7.2) and DIANA-microT database [22, 23]. The top 5 % predicted targets of each DEM in these two datasets were obtained and their overlapped genes were selected. Pathway analysis was conducted by Kyoto Encyclopedia of Genes and Genomes (KEGG) database [24,25,26] and Reactome database [27]. Rnacentral database [28] was used to explore the literatures that were relevant to miR135a-1-3p. All figures and statistics method were performed by R 4.0.

Validation target gene of mir135a by dual-luciferase report experiment

293T cells were plated onto the 48-well plates at 37 °C and cultured to reach 70 %-80 % confluence. Fragment of the 3’UTR of Hoxd10 was inserted into the psiCHECK-2 Luciferase vector (Promega, Madison, Wisconsin, USA). The designed mmu-miR-135a-1-3p mimics and NC mimics was used to co-transfect 293T cells together with the corresponding plasmids (Wide type Hoxd10) using Lipofectamine 2000 for 48 h. For assessment of action of mmu-miR-135a-1-3p on Hoxd10, cells were grouped as following: (A) NC group, cells co-transfected with WT Hoxd10 plasmid and NC mimics; (B) mmu-miR-135a-1-3p mimics, cells co‐transfected with WT Hoxd10 plasmid and mmu-miR-135a-1-3p mimics group. The luciferase report activity was measured using Dual Luciferase reporter 1000 Assay System (Promega). The data was normalized by the rellina luciferase activity of control group. Each experiment was performed in triplicate.


Differentially expressed miRNAs in activated IHH signaling pathway

Cells were divided into WT, MT and CT groups according to their treatments (activated by WT IhhN, activated by MT IhhN and unactivated respectively). There were two biological replicates within each group. 650 probe features were derived from the microarray. Among them, we selected 169 probe features that had gene information to perform further analysis. The difference between the spots and the median of the same spot across all the arrays were evaluated by the relative log expression (RLE), which was computed for every spot in the array. The majority of the spots was centered around zero, indicating that most of the miRNAs among three groups were not differentially expressed. A few miRNAs were significantly differentially expressed across the three groups (Fig. 1, Additional file 1, 2, 3 and 4). As shown in Table 1, there were 6 DEMs between WT and CT group; 5 DEMs between MT and CT groups (P-value < 0.05). In particular, miR-135a-1-3p was the only DEM between WT and MT group (Table 1).

Fig. 1
figure 1

Volcano plot for Differentially expressed miRNAs in activated IHH signaling pathway. A The differentially expressed miRNAs between WT an MT group; B The differentially expressed miRNAs between MT an CT group; C The differentially expressed miRNAs between WT an CT group; D venn chart for the three group. WT: wild type group, MT: p.E95K mutant group, CT: control group. Red dot means the differentially expressed miRNAs

Table 1 Differentially Expressed miRNAs

Predicted targets of DEMs

The top 5 % predicted targets of each DEMs using TargetScan database [21] (Version 7.2) and DIANA-microT database [22, 23] were obtained. The common predicted genes from two datasets are shown in Additional file 5. Pathway analyses were conducted based on Kyoto Encyclopedia of Genes and Genomes (KEGG) database [24,25,26] and Reactome database [27] to find the enrichment of DEM targets. The top 5 most significantly related pathways are shown in Fig. 2 A and Supplementary files (Additional files 6, 7 and 8). Runx1/2-related pathways were significantly enriched in the targets of up-regulated miRNAs (4 of the top 5 pathways, Additional file 6). Notch3 signaling pathways were significantly enriched in the targets of down-regulated miRNAs (2 of the top 5 pathways, Additional file 7). Runx1/2 and collagen-related pathways were enriched in the targets of all DEMs (Additional file 8). Notch signaling pathway was enriched in the targets of miR135-1-3p (Additional file 9). These pathways mentioned above had been reported to be closely related to chondrogenesis [29, 30].

Fig. 2
figure 2

Heatmap and the result of dual-luciferase report experiment. A Heatmap for enrichment of the differentially expressed miRNAs targets. B The dual-luciferase reporter gene experiment for miR135-1-3p and its target Hoxd10. **P < 0.01, analyzed by t test

MiR135-1-3p regulates the expression of Hoxd10 in 293T cells

In this step, we validated the relationship between miR135-1-3p and Hoxd10 to confirm if miR135-1-3p regulated the process of endochondral osteogenesis. Hoxd10 is one of the predicted target gene of miR135-1-3p. Hoxd10 has been reported to be expressed in the developing limb buds and was involved in differentiation and limb development [31]. The dual-luciferase reporter gene experiment was performed to detect if miR135-1-3p directly combined with the 3′-UTR of the Hoxd10 mRNA. Compared to the control group (Mir-NC), the luciferase activity of the Hoxd10 of the miR-135a-1-3p group was significantly lower (Fig. 2B, P = 0.0014).

Validation of potential targets of DEMs by gene expression profiles

Expressions of predicted targets were examined based on mRNA expression profiles that had been described in our previous study [19]. We successfully identified 9 target genes (Table 2). MiRNAs silence the target genes by base pairing with complementary sequences within mRNA molecules. Therefore, we considered a pair of a miRNA and its regulation target were a gene and a DEM with opposite expression states. Our result showed that miR-29-3p/Col4a5, miR-29-3p/Dpysl5, miR-712-5p/Gpr116 and miR-101a-3p/Nova1 were highly confident miRNA-gene pairs (Absolute fold change cutoff = 2). However, the function of these genes did not appear to be related to endochondral osteogenesis.

Table 2 Pairs of DEMs and their predicted target genes


Previous studies have revealed that the stability of E95K mutant IHH and wild type IHH protein were similar [17]. However, a recent study has reported a different conclusion, suggesting that the mutant IHH impaired its interactions with its partners in mouse model with digit abnormalities [6]. In addition, several studies on molecular mechanism of downstream regulation involving the mutant IHH signaling have been published. To our knowledge, this study is the first attempt to characterize the miRNA-mediated regulation in BDA1-related mutant IHH signaling and demonstrated that a number of chondrogenesis-relevant pathways, including Runx1/2-related, Notch3 signaling and collagen-related pathways, were significantly enriched in the target genes of DEMs. Our results suggested that miRNA regulation in IHH signaling might play a crucial role in chondrogenesis.

The result of differential miRNA profiling showed that Runx1/2-related pathways were enriched in the DEMs. Earlier studies have shown that Runx1 contributed to parathyroid hormone (PTH)-induced chondrogenesis and the induction of mesenchymal stem cell commitment to the early stages of chondrogenesis [29]. The knockdown of Runx1 completely blunted PTH-mediated chondrogenesis [29, 32]. IHH stimulates the expression of PTH-related peptide (PTHrP), a protein that negatively regulates terminal chondrocyte differentiation through the PTH/PTHrP receptor (PPR) [33]. Moreover, it has been reported that the duplicated copy number of PTHrP caused BDA1 phenotype [34]. Recent studies have found that the copy number variation of RUNX2 resulted in brachydactyly-liked phenotypes [35,36,37]. Taking together, our results indicated that the identified DEMs might be responsible for the pathogenesis of BDA1 by exerting influence on the regulation of Runx1/2.

Our study identified a miRNA (miR-135a-1-3p) had the significant difference in expression level between WT and MT group. The dual-luciferase reporter gene experiment verified that Hoxd10 was a direct target of miR-135a-1-3p, who may suppress the expression of Hoxd10. Hoxd10 encodes a sequence-specific transcription factor, which has been found to be expressed during the developing of limb buds. Hoxd10 is critical in limb formation and limb morphogenesis [38]. Carpenter’s team has reported that targeted disruption of Hoxd10 gene caused hindlimb-specific defects in gait and adduction in mice [39]. Hence, miR-135a-1-3p may regulate the IHH signaling by target Hoxd10.

The results of pathway analysis revealed that collagen fibrillogenesis-related pathways were enriched in the target genes of miR135-1-3p. Collagens have been recognized as candidate indicators of the cartilage phenotype in process of chondrogenesis and chondrocyte differentiation [40]. PTHrP/Ihh axis has been found to control the expression of collagen VI by proliferating chondrocytes during the development of limb [41]. Therefore, miR-135a-1-3p may be associated with the pathogenesis of BDA1 and chondrogenesis. Furthermore, it is noteworthy that the expression levels of mmu-let-7c-1-3p in both WT and MT group were lower than the one in CT group. However, the expression levels of mmu-let-7c-1-3p in WT and MT group appeared to have no difference. NOTCH3 signaling pathway was the most significant pathway enriched in target genes of mmu-let-7c-1-3p. It has been reported that Notch signaling plays an important role in the coordination of perichondrial osteoblast differentiation, regulation of the communication between chondrocytes and perichondrial osteoblasts and promotion of chondrocyte proliferation and apoptosis. Further, the inhibitory roles of Indian hedgehog (Ihh) downregulation on chondrocyte growth and differentiation [42]. Therefore, our findings indicated that mmu-let-7c-1-3p may be regulate the IHH signaling pathway to inhibit chondrocyte differentiation.


This is the first study attempted to reveal the role of miRNA regulation involved in the pathogenesis of BDA1. Six miRNAs were found to have significant difference in expression level between induced and control group. We discovered that miR-135a-1-3p was a significantly differentially expressed miRNA. Hoxd10 was identified as a target of miR-135a-1-3p, who is involved in chondrogenesis. Pathway analysis revealed that Runx1/2, Notch and collagen-related pathways were enriched in the top-confident predicted targets of these six miRNAs, which were closely related to chondrogenesis. Our results indicated the important role of miRNA-regulation in the pathogenesis of BDA1 and chondrogenesis process.

Availability of data and materials

All requests for raw and analyzed data will be reviewed by the corresponding author. The detailed information of the primer is shown in the supplementary data.



Indian Hedgehog


brachydactyly type A1




wild type


p.E95K mutant




sonic hedgehog


desert hedgehog




differentially expressed miRNAs


IHH-N-terminus proteins


fetal bovine serum


Kyoto Encyclopedia of Genes and Genomes


Minimum Information for Publication of Quantitative Real-Time PCR Experiment




relative log expression


parathyroid hormone


PTH-related peptide


  1. Maes C: Signaling pathways effecting crosstalk between cartilage and adjacent tissues: Seminars in cell and developmental biology: The biology and pathology of cartilage. Semin Cell Dev Biol 2017, 62:16–33.

    Article  CAS  Google Scholar 

  2. Tang J, Su N, Zhou S, Xie Y, Huang J, Wen X, Wang Z, Wang Q, Xu W, Du X et al: Fibroblast Growth Factor Receptor 3 Inhibits Osteoarthritis Progression in the Knee Joints of Adult Mice. Arthritis Rheumatol 2016, 68(10):2432–2443.

    Article  CAS  Google Scholar 

  3. Zhou J, Wei X, Wei L: Indian Hedgehog, a critical modulator in osteoarthritis, could be a potential therapeutic target for attenuating cartilage degeneration disease. Connect Tissue Res 2014, 55(4):257–261.

    Article  CAS  Google Scholar 

  4. St-Jacques B, Hammerschmidt M, McMahon AP: Indian hedgehog signaling regulates proliferation and differentiation of chondrocytes and is essential for bone formation. Genes Dev 1999, 13(16):2072–2086.

    Article  CAS  Google Scholar 

  5. Hooper JE, Scott MP: Communicating with Hedgehogs. Nat Rev Mol Cell Biol 2005, 6(4):306–317.

    Article  CAS  Google Scholar 

  6. Gao B, Hu J, Stricker S, Cheung M, Ma G, Law KF, Witte F, Briscoe J, Mundlos S, He L et al: A mutation in Ihh that causes digit abnormalities alters its signalling capacity and range. Nature 2009, 458(7242):1196–1200.

    Article  CAS  Google Scholar 

  7. Yang J, Andre P, Ye L, Yang YZ: The Hedgehog signalling pathway in bone formation. Int J Oral Sci 2015, 7(2):73–79.

    Article  CAS  Google Scholar 

  8. Lin Z, Murtaza I, Wang K, Jiao J, Gao J, Li PF: miR-23a functions downstream of NFATc3 to regulate cardiac hypertrophy. Proc Natl Acad Sci U S A 2009, 106(29):12103–12108.

    Article  Google Scholar 

  9. Wang Y, Liang Y, Lu Q: MicroRNA epigenetic alterations: predicting biomarkers and therapeutic targets in human diseases. Clin Genet 2008, 74(4):307–315.

    Article  CAS  Google Scholar 

  10. Yan LX, Huang XF, Shao Q, Huang MY, Deng L, Wu QL, Zeng YX, Shao JY: MicroRNA miR-21 overexpression in human breast cancer is associated with advanced clinical stage, lymph node metastasis and patient poor prognosis. RNA 2008, 14(11):2348–2360.

    Article  CAS  Google Scholar 

  11. Dong J, Cui X, Jiang Z, Sun J: MicroRNA-23a modulates tumor necrosis factor-alpha-induced osteoblasts apoptosis by directly targeting Fas. J Cell Biochem 2013, 114(12):2738–2745.

    Article  CAS  Google Scholar 

  12. Hassan MQ, Maeda Y, Taipaleenmaki H, Zhang W, Jafferji M, Gordon JA, Li Z, Croce CM, van Wijnen AJ, Stein JL et al: miR-218 directs a Wnt signaling circuit to promote differentiation of osteoblasts and osteomimicry of metastatic cancer cells. J Biol Chem 2012, 287(50):42084–42092.

    Article  CAS  Google Scholar 

  13. Kureel J, John AA, Dixit M, Singh D: MicroRNA-467 g inhibits new bone regeneration by targeting Ihh/Runx-2 signaling. Int J Biochem Cell Biol 2017, 85:35–43.

    Article  CAS  Google Scholar 

  14. Bae Y, Yang T, Zeng HC, Campeau PM, Chen Y, Bertin T, Dawson BC, Munivez E, Tao J, Lee BH: miRNA-34c regulates Notch signaling during bone development. Hum Mol Genet 2012, 21(13):2991–3000.

    Article  CAS  Google Scholar 

  15. Gao B, He L: Answering a century old riddle: brachydactyly type A1. Cell Res 2004, 14(3):179–187.

    Article  CAS  Google Scholar 

  16. Byrnes AM, Racacho L, Grimsey A, Hudgins L, Kwan AC, Sangalli M, Kidd A, Yaron Y, Lau YL, Nikkel SM et al: Brachydactyly A-1 mutations restricted to the central region of the N-terminal active fragment of Indian Hedgehog. Eur J Hum Genet 2009, 17(9):1112–1120.

    Article  CAS  Google Scholar 

  17. Ma G, Yu J, Xiao Y, Chan D, Gao B, Hu J, He Y, Guo S, Zhou J, Zhang L et al: Indian hedgehog mutations causing brachydactyly type A1 impair Hedgehog signal transduction at multiple levels. Cell Res 2011, 21(9):1343–1357.

    Article  CAS  Google Scholar 

  18. Guo S, Zhou J, Gao B, Hu J, Wang H, Meng J, Zhao X, Ma G, Lin C, Xiao Y et al: Missense mutations in IHH impair Indian Hedgehog signaling in C3H10T1/2 cells: Implications for brachydactyly type A1, and new targets for Hedgehog signaling. Cell Mol Biol Lett 2010, 15(1):153–176.

    Article  CAS  Google Scholar 

  19. Shen L, Ma G, Shi Y, Ruan Y, Yang X, Wu X, Xiong Y, Wan C, Yang C, Cai L et al: p.E95K mutation in Indian hedgehog causing brachydactyly type A1 impairs IHH/Gli1 downstream transcriptional regulation. BMC Genet 2019, 20(1):10.

  20. Pathi S, Pagan-Westphal S, Baker DP, Garber EA, Rayhorn P, Bumcrot D, Tabin CJ, Blake Pepinsky R, Williams KP: Comparative biological responses to human Sonic, Indian, and Desert hedgehog. Mech Dev 2001, 106(1-2):107–117.

    Article  CAS  Google Scholar 

  21. Agarwal V, Bell GW, Nam JW, Bartel DP: Predicting effective microRNA target sites in mammalian mRNAs. Elife 2015, 4.

  22. Reczko M, Maragkakis M, Alexiou P, Grosse I, Hatzigeorgiou AG: Functional microRNA targets in protein coding sequences. Bioinformatics 2012, 28(6):771–776.

    Article  CAS  Google Scholar 

  23. Paraskevopoulou MD, Georgakilas G, Kostoulas N, Vlachos IS, Vergoulis T, Reczko M, Filippidis C, Dalamagas T, Hatzigeorgiou AG: DIANA-microT web server v5.0: service integration into miRNA functional analysis workflows. Nucleic Acids Res 2013, 41(Web Server issue):W169-173.

    Google Scholar 

  24. Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K: KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res 2017, 45(D1):D353-D361.

    Article  CAS  Google Scholar 

  25. Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res 2000, 28(1):27–30.

    Article  CAS  Google Scholar 

  26. Kanehisa M, Sato Y, Furumichi M, Morishima K, Tanabe M: New approach for understanding genome variations in KEGG. Nucleic Acids Res 2019, 47(D1):D590-D595.

    Article  CAS  Google Scholar 

  27. Fabregat A, Jupe S, Matthews L, Sidiropoulos K, Gillespie M, Garapati P, Haw R, Jassal B, Korninger F, May B et al: The Reactome Pathway Knowledgebase. Nucleic Acids Res 2018, 46(D1):D649-D655.

    Article  CAS  Google Scholar 

  28. The RC, Petrov AI, Kay SJE, Kalvari I, Howe KL, Gray KA, Bruford EA, Kersey PJ, Cochrane G, Finn RD et al: RNAcentral: a comprehensive database of non-coding RNA sequences. Nucleic Acids Res 2017, 45(D1):D128-D134.

    Google Scholar 

  29. Wang Y, Belflower RM, Dong YF, Schwarz EM, O’Keefe RJ, Drissi H: Runx1/AML1/Cbfa2 mediates onset of mesenchymal cell differentiation toward chondrogenesis. J Bone Miner Res 2005, 20(9):1624–1636.

    Article  CAS  Google Scholar 

  30. Karlsson C, Lindahl A: Notch signaling in chondrogenesis. Int Rev Cell Mol Biol 2009, 275:65–88.

    Article  CAS  Google Scholar 

  31. Raines AM, Magella B, Adam M, Potter SS: Key pathways regulated by HoxA9,10,11/HoxD9,10,11 during limb development. BMC Dev Biol 2015, 15:28.

    Article  CAS  Google Scholar 

  32. Duan S, Dong X, Hai J, Jiang J, Wang W, Yang J, Zhang W, Chen C: MicroRNA-135a-3p is downregulated and serves as a tumour suppressor in ovarian cancer by targeting CCR2. Biomed Pharmacother 2018, 107:712–720.

    Article  CAS  Google Scholar 

  33. Kobayashi T, Chung UI, Schipani E, Starbuck M, Karsenty G, Katagiri T, Goad DL, Lanske B, Kronenberg HM: PTHrP and Indian hedgehog control differentiation of growth plate chondrocytes at multiple steps. Development 2002, 129(12):2977–2986.

    Article  CAS  Google Scholar 

  34. Flottmann R, Sowinska-Seidler A, Lavie J, Chateil JF, Lacombe D, Mundlos S, Horn D, Spielmann M: Duplication of PTHLH causes osteochondroplasia with a combined brachydactyly type E/A1 phenotype with disturbed bone maturation and rhizomelia. Eur J Hum Genet 2016, 24(8):1132–1136.

    Article  CAS  Google Scholar 

  35. Al-Yassin A, Calder AD, Harrison M, Lester T, Lord H, Oldridge M, Watkins S, Keen R, Wakeling EL: A three-generation family with metaphyseal dysplasia, maxillary hypoplasia and brachydactyly (MDMHB) due to intragenic RUNX2 duplication. Eur J Hum Genet 2018, 26(9):1288–1293.

    Article  CAS  Google Scholar 

  36. Avela K, Hirvinen H, Ben Amor M, Rauch F: Metaphyseal dysplasia with maxillary hypoplasia and brachydactyly in a Finnish woman: first confirmation of a duplication in RUNX2 as pathogenic variant. Eur J Med Genet 2014, 57(11-12):617–620.

    Article  Google Scholar 

  37. Moffatt P, Ben Amor M, Glorieux FH, Roschger P, Klaushofer K, Schwartzentruber JA, Paterson AD, Hu P, Marshall C, Consortium FC et al: Metaphyseal dysplasia with maxillary hypoplasia and brachydactyly is caused by a duplication in RUNX2. Am J Hum Genet 2013, 92(2):252–258.

    Article  CAS  Google Scholar 

  38. Gurnett CA, Keppel C, Bick J, Bowcock AM, Dobbs MB: Absence of HOXD10 mutations in idiopathic clubfoot and sporadic vertical talus. Clin Orthop Relat Res 2007, 462:27–31.

    Article  Google Scholar 

  39. Carpenter EM, Goddard JM, Davis AP, Nguyen TP, Capecchi MR: Targeted disruption of Hoxd-10 affects mouse hindlimb development. Development 1997, 124(22):4505–4514.

    Article  CAS  Google Scholar 

  40. Vickers SM, Squitieri LS, Spector M: Effects of cross-linking type II collagen-GAG scaffolds on chondrogenesis in vitro: dynamic pore reduction promotes cartilage formation. Tissue Eng 2006, 12(5):1345–1355.

    Article  CAS  Google Scholar 

  41. Onyekwelu I, Goldring MB, Hidaka C: Chondrogenesis, joint formation, and articular cartilage regeneration. J Cell Biochem 2009, 107(3):383–392.

    Article  CAS  Google Scholar 

  42. Deng A, Zhang H, Hu M, Liu S, Wang Y, Gao Q, Guo C: The inhibitory roles of Ihh downregulation on chondrocyte growth and differentiation. Exp Ther Med 2018, 15(1):789–794.

    PubMed  CAS  Google Scholar 

Download references


The authors gratefully acknowledge the support and advice by Lianyi Yang at Shanghai Jiao Tong University in Shanghai.

Suggested reviewers

Jingsong Ma,, Westlake University.

Ren Cai,, Shanghai Ninth People’s Hospital.

Jinghang zhu,, Anhui Medical University.

Wei Zhang,, Department of Clinical Pharmacology, Xiangya Hospital, Central South University.

Hailiang Huang,, Analytic and Translational Genetics Unit, Department of Medicine, Massachusetts General Hospital and Harvard.


This work was supported by grants from the Interdisciplinary Program of Shanghai Jiao Tong University (project number YG2019QNA51), the 863 Program (2012AA02A515, 2012AA021802), the National Nature Science Foundation of China (81773818, 81273596, 30900799, 81671326), National key research and development program (2017YFC0909303, 2016YFC0905000, 2016YFC0905002, 2016YFC1200200, 2016YFC0906400), The 4th Three-year Action Plan for Public Health of Shanghai (The Project No.: 15GWZK0101), Shanghai Pujiang Program (17PJD020), Shanghai Key Laboratory of Psychotic Disorders (13dz2260500), Guangdong Province Science and Technology Program (grant no: 2020A0505100062), Clinical Innovation Research Program of Guangzhou Regenerative Medicine and Health Guangdong Laboratory(2018GZR0201002).

Author information

Authors and Affiliations



Wei Zhou, Lu Shen, Xiaofang Sun and Gang Ma have contributed to the conceptualization of this study. Hao Wu, Ting Wang and Yingtian Zhang have contributed to data curation. Wei Zhou, Luan Chen, Xiaofang Sun, Lu Shen, Ting Wang and Baocheng Wang have contributed to formal analysis. Hai Zhang and Lu Shen have contributed to funding acquisition. Lin He, Xiaofang Sun and Shengying Qin have contributed to the project supervision. Wei Zhou, Luan Chen, Cong Wang and Na Zhang have contributed to draft writing and reviewing. Wei Zhou, Lu Shen, Xiaofang Sun and Hai Zhang read and approved the final manuscript.

Corresponding authors

Correspondence to Xiaofang Sun, Hai Zhang or Lu Shen.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Mean signal of microarrays.

Additional file 2.

M-A plot of WT vs CT.

Additional file 3.

M-A plot of MT vs CT.

Additional file 4.

M-A plot of WT vs MT.

Additional file 5.

Pathway analysis of targets of all DEMs.

Additional file 6.

Pathway analysis of targets of all up-regulated DEMs.

Additional file 7.

Pathway analysis of targets of all down-regulated DEMs.

Additional file 8.

Pathway analysis of targets of all DEMs.

Additional file 9.

Pathway analysis of targets of miR-135a-1-3p.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhou, W., Chen, L., Wu, H. et al. Altered microRNAs in C3H10T1/2 cells induced by p.E95K mutant IHH signaling. Hereditas 158, 48 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: