- Research
- Open access
- Published:
Identification and functional analysis of a three-miRNA ceRNA network in hypertrophic scars
Journal of Translational Medicine volume 19, Article number: 451 (2021)
Abstract
Background
Hypertrophic scar (HTS) is a fibrotic disorder of skins and may have repercussions on the appearance as well as functions of patients. Recent studies related have shown that competitive endogenous RNA (ceRNA) networks centering around miRNAs may play an influential role in HTS formation. This study aimed to construct and validate a three-miRNA (miR-422a, miR-2116-3p, and miR-3187-3p) ceRNA network, and explore its potential functions.
Methods
Quantitative real‑time PCR (qRT‑PCR) was used to compare expression levels of miRNAs, lncRNAs, and genes between HTS and normal skin. Target lncRNAs and genes of each miRNA were predicted using starBase as well as TargetScan database to construct a distinct ceRNA network; overlapping target lncRNAs and genes of the three miRNAs were utilized to develop a three-miRNA ceRNA network. For every network, protein–protein interaction (PPI) network analysis was performed to identify its hub genes. For each network and its hub genes, Gene Oncology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were conducted to explore their possible functions.
Results
MiR-422a, miR-2116-3p, and miR-3187-3p were all downregulated in HTS tissues and fibroblasts. MiR-422a-based ceRNA network consisted of 101 lncRNAs with 133 genes; miR-2116-3p-centered ceRNA network comprised 85 lncRNAs and 978 genes; miR-3187-3p-derived ceRNA network encompassed 84 lncRNAs as well as 1128 genes. The three-miRNA ceRNA network included 2 lncRNAs with 9 genes, where MAPK1, FOSL2, ABI2, KPNA6, CBL, lncRNA-KCNQ1OT1, and lncRNA-EBLN3P were upregulated. According to GO and KEGG analysis, these networks were consistently related to ubiquitination. Three ubiquitination-related genes (CBL, SMURF2, and USP4) were upregulated and negatively correlated with the expression levels of the three miRNAs in HTS tissues.
Conclusions
This study identified a three-miRNA ceRNA network, which might take part in HTS formation and correlate with ubiquitination.
Background
Hypertrophic scar (HTS) is a fibro-proliferative disease, charactered by excessive proliferation of fibroblasts, activation of myofibroblasts, and deposition of extracellular matrix (ECM). This disorder has the potential to cause disfiguration and dysfunction in facial areas and joints, which may take a heavy toll on the physical and psychological condition of patients. Despite increasing strategies towards this challenge, in an effort to seek more satisfactory management, the pathogenesis underlying HTS formation remains to be elucidated [1].
MicroRNAs (miRNAs) can be competitively sponged by mRNAs as well as long non-coding RNAs (lncRNAs), and these can form competitive endogenous RNA (ceRNA) networks, which have been reported to regulate multiple biological processes to engage in the pathogenesis of many diseases including HTS [2, 3]. For example, miR-422a has been shown to have the potential to inhibit cell proliferation in glioblastoma [4] and fibrotic genes in liver fibrosis [5]; miR-2116-3p has proven to hinder the proliferation of breast cancer cells [6]; miR-3187-3p has been reported to impede fibroblast proliferation [7]. Furthermore, in a preliminary microarray test, miR-422a, miR-2116-3p, and miR-3187-3p have been discovered to be downregulated in HTS tissues (fold change > 2 and P < 0.05) [8]. These studies have suggested that the three miRNAs may correlate with HTS formation and contribute positively to alleviating HTS. However, their expression patterns remain validated in tissues as well as fibroblasts, and their roles in HTS pathogenesis need to be answered.
A single lncRNA or mRNA tends to bear multiple binding sites for miRNAs, and can simultaneously sponge several miRNAs to constitute a sophisticated ceRNA network. This one-to-many pattern plays a crucial role in a range of fibrotic disorders including HTS. In identical tissues from patients suffering from HTS, dysregulated expression profiles of lncRNAs, miRNAs, and mRNAs have been observed to coexist and correlate [9]. Moreover, miR-23b-3p and miR-27b-3p have been confirmed to promote atrial fibrosis by targeting two distinct sites of transforming growth factor (TGF)-β1 receptor 3 (TGFBR3) [10]. In addition, miR-760 and miR-215-3p have been demonstrated to be competitively sponged by lncRNA NR_003923 and IL22RA1, promoting TGF-β-induced proliferation, migration, and autophagy of fibroblasts [11]. Nevertheless, it needs to be further investigated whether miR-422a, miR-2116-3p, and miR-3187-3p can be concurrently sponged by ceRNAs, contributing to the pathogenesis of fibrosis of many organs including skins.
In this study, we aimed to determine both expression profiles and possible functions of a three-miRNA (miR-422a, miR-2116-3p, and miR-3187-3p) ceRNA network in HTS. Therefore, we detected expression levels of these three miRNAs in HTS tissue samples and HTS-derived fibroblasts. In order to explore their roles in the pathogenesis of HTS, we predicted both target lncRNAs as well as genes and constructed ceRNA networks for each of the three miRNAs. These networks were subject to protein–protein interaction (PPI) network analysis to detect hub genes. GO and KEGG enrichment analysis was utilized to explore potential functions of both the networks and hub genes. Furthermore, we performed integrated bioinformatics analysis of miR-422a, miR-2116-3p, and miR-3187-3p to establish a comprehensive ceRNA network (Fig. 1), whose roles and expression patterns were analyzed through functional enrichment analysis and quantitative real-time PCR (qRT-PCR).
Methods
Tissue collection and cell culture
HTS tissues and normal skin samples were collected from 10 patients, who were informed about the aim as well as procedure of the study, signed written informed consent, and underwent HTS excision. Primary fibroblasts were isolated from HTS as well as normal skin, and cultured as previously described [12]. Briefly, tissues were minced and then placed in collagenase type IV (0.2 mg/ml, Roche) solution at 37 °C for 6 h. After filtering and centrifuging, fibroblasts were cultured in DMEM (Gibco, USA) supplemented with 10% fetal bovine serum (Gibco, USA) as well as 1% penicillin/streptomycin (Gibco, USA) in 5% CO2 at 37 °C. Fibroblasts between passage 3 and passage 9 were utilized for experiments. The study protocol was approved by the Medical Ethics Committee of The First Affiliated Hospital of Zhengzhou University and Shanghai Ninth People’s Hospital, Shanghai Jiaotong University School of Medicine.
Quantitative real-time PCR (qRT-PCR)
Total RNA and miRNA were extracted form tissues and fibroblasts through TRIzol Reagent (Invitrigen, USA). 2 ng of total RNA and miRNA was subject to cDNA synthesis. Expression levels of protein-coding genes and miRNAs were normalized to those of GAPDH and U6, respectively. Primer sequences used in this research were listed as follows: GAPDH: 5’-GGGAAGGTGAAGGTCGGAGT-3’ (forward) and 5’-GGGGTCATTGATGGCAACA-3’ (reverse); KCNQ1OT1: 5’-GGCTTCTCGTTCAGGTGACA-3' (forward) and 5'-TGGTACATGGGCAGGGGATA-3' (reverse); EBLN3P: 5'-AATTTCGGGAGCTGTGACCC-3' (forward) and 5'-CAAAAGACTGCTGCCCAGTG-3' (reverse); ABI2: 5'-ACATGAAGATGGGTGGGCTG-3' (forward) and 5'-CTCCCACTGCTGCTGTAAGT-3' (reverse); FOSL2: 5'-AAACACCCTGTTTCCTCTCCG-3' (forward) and 5'-ATCTACCCGGAATTTCTGCTG-3' (reverse); MAPK1: 5'-CGAGTGACGAGCCCATCG-3' (forward) and 5'-GATGTCTGAGCACGTCCAGT-3' (reverse); CBL: 5'-ACGTGAAGAAGAGCTCTGGG-3' (forward) and 5'-ACAACCGCACCACCTTGTC-3' (reverse); KPNA6: 5'-TACTGAAAGCTGCCGCTGAA-3' (forward) and 5'-GCTCGCCATGGTCTCCATC-3' (reverse); SMURF2: 5'-CAGAGGACAACGCAACAAGG-3' (forward) and 5'-GTTGCTAAGATCCCTGGGCA-3' (reverse); USP4: 5'-GGTACGGCTGTGTAGAAGGC-3' (forward) and 5'-TCTCGATGGTTGCAATGGTG-3' (reverse); U6: 5’-CGCTAGCACATATCGGCTA-3’ (forward) and 5’-TTCTGCGACGAATTTGTCAT-3’ (reverse); miR-422a: 5’-GCATACCGCTATGCCTAATGGTG-3’ (forward) and 5’-GTGCAGGAGGTTCCGGT-3’ (reverse); miR-2116-3p: 5’-AATCCTATGCCAAGAACTCCC-3’ (forward) and 5’- CTCTACAGCTATATTGCCAGCCA-3’ (reverse); miR-3187-3p: 5’-GCCGAGTTGGCCATGGGGCTG-3’ (forward) and 5’- CTCAACTGGTGTCGTGGA-3’ (reverse).
Construction of ceRNA networks
For each of miR-422a, miR-2116-3p, and miR-3187-3p, its target lncRNAs were predicted using starBase (http://starbase.sysu.edu.cn/); its target genes were forecasted by overlapping results from starBase as well as TargetScan (http://www.targetscan.org). Target lncRNAs and genes of every miRNA were utilized to develop its ceRNA network; overlapping target lncRNAs and genes of the three miRNAs were used to build up a three-miRNA ceRNA network. These networks were constructed through Cytoscape (version 3.7.0, http://www.cytoscape.org/).
PPI network analysis
For target genes of every miRNA, a single PPI network was established using the Search Tool for the Retrieval of Interacting Genes (STRING) database (https:// string-db.org) and visualized through Cytoscape. The Cytoscape plugin, cytoHubba, was utilized to analyze each PPI network in an effort to identify its hub genes, whose ranks were calculated according to the built-in maximal clique centrality (MCC) algorithm. For genes in the three-miRNA ceRNA network, the PPI network was developed using interacting proteins, which were predicted to interact with the mutual genes by NetworkAnalyst (http://www.networkanalyst.ca/).
GO and KEGG enrichment analysis
Using DAVID database (https://david.ncifcrf.gov/summary.jsp), Gene Ontology (GO) as well as Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis was used to explore functions of genes in the ceRNA networks, in order to predict potential roles of these networks in HTS pathogenesis. The results of GO annotation include terms about the biological process (BP), cellular component (CC), and molecular function (MF); outcomes of KEGG analysis encompass pathways involved in a variety of physiological and pathological processes. Adjusted P value < 0.05 was set as the cut-off value for GO and KEGG enrichment analysis.
Statistical analysis
Statistical analysis was carried out using GraphPad Prism (version 9.0.1). Expression levels of miRNAs, lncRNAs, and genes between HTS and control groups were presented as mean ± SD. Student’s t-test and Mann–Whitney test were performed for data with equal and unequal variances, respectively. P < 0.05 was considered to reach statistical significance.
Results
miR-422a, miR-2116-3p, and miR-3187-3p are downregulated in HTS tissues and fibroblasts
To compare expression levels of miR-422a, miR-2116-3p, and miR-3187-3p between HTS and normal skin, 10 HTS tissues as well as normal skin samples were subject to qRT-PCR. The results revealed that each expression of these three miRNAs was significantly lower in HTS tissues as opposed to normal skin samples (Fig. 2a–c). To further confirm the expression pattern, qRT-PCR was employed to dissect mRNA levels of the three miRNAs in fibroblasts derived from HTS and normal skin tissues. The outcomes were consistent with those in tissues and indicated that all of the three miRNAs were downregulated in HTS fibroblasts (Fig. 2d–f). Taken together, these findings suggest that miR-422a, miR-2116-3p, and miR-3187-3p may have a close relationship with HTS formation, and play a protective role in the process.
Construction of ceRNA networks
Given that in the pathogenesis of most diseases including HTS, ceRNA networks are the major functional pattern of miRNAs, both upstream lncRNAs and downstream genes were predicted for each of miR-422a, miR-2116-3p, and miR-3187-3p to construct ceRNA networks. As for miR-422a, a total of 101 target lncRNAs were obtained using the starBase database; 133 target genes were acquired after overlapping data from starBase and TargetScan (Fig. 3a). Using these lncRNAs and mRNAs, a ceRNA network centering around miR-422a was constructed (Fig. 3b). In identical ways, by overlapping results from starBase and TargetScan databases (Fig. 3c), the ceRNA network for miR-2116-3p was established with 85 lncRNAs and 978 target genes (Fig. 3d). Similarly, concerning miR-3187-3p, 84 target lncRNAs and 1128 genes were predicted through starBase and intersection of outcomes in starBase as well as TargetScan, respectively (Fig. 3e). The lncRNAs and target genes were then utilized to set up a ceRNA network of miR-3187-3p (Fig. 3f).
Establishment of PPI networks and identification of hub genes
In light of the fact that miRNAs tend to regulate HTS development through interacting with mRNAs of protein-coding genes, target genes in each of the three ceRNA networks were utilized to build up a distinct PPI network in order to find hub genes. The PPI network of miR-422a encompassed 87 edges and 133 nodes, and analysis of this network detected 10 hub genes (Fig. 4a). In terms of miR-2116-3p, its PPI network consisted of 3926 edges and 975 nodes, 20 of which were shown to be hub genes (Fig. 4b). Concerning miR-3187-3p, the PPI network involved 5047 edges and 1125 nodes, the analysis of which identified 20 hub genes (Fig. 4c).
GO and KEGG enrichment analysis
To gain insights into potential roles of miR-422a, miR-2116-3p, and miR-3187-3p, GO and KEGG enrichment analyses were conducted for both target and hub genes in the ceRNA networks. For the miR-422a-centered network, target and hub genes were both closely related with ubiquitin binding, PI3K/Akt/mTOR signaling pathway, and signaling pathways regulating pluripotency of stem cells (Fig. 5a, b). As to the network centering miR-2116-3p, target and hub genes were consistently enriched in proteasome-mediated ubiquitin-dependent protein catabolic process, ubiquitin-protein transferase activity, and ubiquitin-mediated proteolysis (Fig. 5c, d). Regarding miR-3187-3p, analysis of its network demonstrated that functions of target genes were similar to those of hub genes, focusing on proteasomal protein catabolic process as well as ubiquitin-protein ligase activity (Fig. 5e, f). Taken together, these findings highlight that the three miRNAs may be associated with ubiquitination during HTS formation.
Integrated analysis of a three-miRNA ceRNA network
In an effort to develop a comprehensive ceRNA network based on miR-422a, miR-2116-3p, and miR-3187-3p, their target lncRNAs and genes were intersected to reap 2 mutual lncRNAs (EBLN3P and KCNQ1OT1) (Fig. 6a) and 9 mutual genes (IPO9, ABI2, ZKSCAN1, FOSL2, GFPT1, MAPK1, CBL, CSNK1G1, KPNA6) (Fig. 6b), respectively. The resultant three-miRNA ceRNA network was demonstrated as the diagram (Fig. 6c). To explore its potential roles in HTS pathogenesis, the 9 mutual genes were subject to NetworkAnalyst (http://www.networkanalyst.ca/) to determine their interacting proteins, which were used to build up a PPI network (Fig. 6d) and analyze related GO terms as well as KEGG pathways (Fig. 6e). The enrichment analysis result implied that these genes interacting with the network may be significantly connected to MAPK signaling pathway and ubiquitin-protein ligase binding, which also existed in the outcomes of the ceRNA RNA networks of miR-422a, miR-2116-3p, and miR-3187-3p. Overall, these findings suggest that the three miRNAs may synergistically engage in HTS development through regulating the same genes and form a ceRNA network associated with ubiquitination.
Validation of the three-miRNA ceRNA network
To further confirm our hypothesis, expression levels of some lncRNAs and genes in the three-miRNA network were detected in HTS tissues and normal skin samples. The qRT-PCR results showed that 2 lncRNAs (KCNQ1OT1 and EBLN3P) and 4 protein-coding genes (MAPK1, FOSL2, ABI2, and KPNA6) were significantly overexpressed in HTS (Fig. 7). Furthermore, three ubiquitination-related genes (CBL, SMURF2, and USP4) were validated to be upregulated and negatively related to the expression levels of miR-422a, miR-2116-3p, and miR-3187-3p (Fig. 8). Collectively, these data supported the notion that the three-miRNA ceRNA network may participate in HTS formation and correlate with ubiquitination.
Discussion
Hypertrophic scars are a kind of skin fibrotic disorder and may have repercussions on the quality of life for patients, but the underlying mechanism remains poorly understood [13]. During the process of HTS formation, miRNAs have been observed to play a regulatory role [14]. A preliminary microarray test on HTS tissues has identified top 40 up and downregulated miRNAs, among which miR-422a, miR-2116-3p, and miR-3187-3p displayed relatively great variations between HTS and control samples. Moreover, according to the prediction of starBase and TargetScan databases, the three miRNAs may target molecules that have been confirmed to play an important role in hypertrophic scarring. For example, miR-422a could target MAPK1 which is a key component of the MAPK signaling pathway, which is actively involved in HTS formation; miR-2116-3p might bond with either COL1A1 whose upregulation causes excessive collagen deposition in HTS, or SERPINE1 that is also known as PAI-1 and has proven a fibrotic gene related to HTS [15]; miR-3187-3p was predicted to inhibit IGF1R whose deficiency has been shown to be significantly associated with the fibrosis of skin [16]. These data imply that miR-422a, miR-2116-3p, and miR-3187-3p may be downregulated and take part in biological processes during HTS formation. This study utilized qRT-PCR to confirm the expression pattern, where all of these three miRNAs were significantly lower in both tissues and fibroblasts of HTS than those of normal skin. These results suggest that miR-422a, miR-2116-3p, and miR-3187-3p are downregulated and might perform a role during HTS pathogenesis.
CeRNA networks, where lncRNAs and protein-coding genes competitively sponge miRNAs, have been shown to engage in multiple diseases including fibrosis. For example, miR-422a was shown to target MEF2D and be sponged by lncRNA D63785, which promoted chemoresistance in gastric cancer [17]; a LINC01433/miR-2116-3p/MYC positive feedback loop was identified to facilitate breast cancer cell proliferation, migration, and epithelial to mesenchymal transition (EMT) [6]; lncRNA H19 targeting miR-3187-3p/GAB1 axis acted as a ceRNA to promote fibroblast proliferation, migration, and invasion in HTS [7]. However, the potential roles of the miR-422a and miR-2116-3p in HTS formation remain unclear. In this study, for each of miR-422a, miR-2116-3p, and miR-3187-3p, the target lncRNAs and genes were predicted to construct a ceRNA network. In addition, the target genes of each miRNA were utilized to establish a distinct PPI network, the analysis of which identified its hub genes. Both target genes and hub genes of every miRNA were subject to GO and KEGG enrichment analysis, in order to explore possible functions of these networks.
The present results showed that miR-422a was significantly related to PI3K/Akt/mTOR signaling pathway and signaling pathways regulating pluripotency of stem cells, supported by previous studies. For instance, miR-422a was confirmed to inhibit glioma cell proliferation, invasion, and migration by targeting PIK3CA, which is a critical member of PI3K/Akt signal pathway [18]; miR-422a was proven to directly target AKT1 to inhibit PI3K/Akt pathway and colorectal cancer cell proliferation [19]. In addition, miR-422a was shown to weaken breast cancer stem cell proliferation in vitro and in vivo [20]. Although these results were obtained from cancer research, HTS development is closely highly associated with PI3K/Akt/mTOR signaling pathway [21] and can also be influenced through manipulating and stem cells [22]. Therefore, miR-422a may regulate HTS formation by impacting PI3K/Akt/mTOR signaling pathway or stem cells.
In this study, each ceRNA network for miR-422a, miR-2116-3p, and miR-3187-3p was predicted to be significantly related to ubiquitination, which was also one of the significant terms in the enrichment analysis results of the integrated ceRNA network centering around the three miRNAs. By means of qRT-PCR, we verified the negative correlation between the three miRNAs and three ubiquitination-related genes (CBL, SMURF2, and USP4), which are key components of the ubiquitin–proteasome system and have been shown to engage in the ubiquitination and degradation of some pivotal genes implicated in HTS [23, 24]. Consistent with these findings, miR-422a was reported to directly target SKP2 [25], which is an E3 ubiquitin ligase and was confirmed to promote fibroblast proliferation in pulmonary fibrosis through via ubiquitin–proteasome pathway [26]. In addition, miR-2116-3p was shown to interact with and inhibit PIWIL1 [27], and silencing of PIWIL1 was reported to upregulate FBXW7 [28], which is a component of E3 ubiquitin-protein ligase complex and was observed to induce expression of fibrotic genes including fibronectin during pulmonary fibrosis [29]. Also, miR-3187-3p was confirmed to target ASXL1 [30], a histone H2A deubiquitinase [31]. Collectively, these studies raise the possibility that miR-422a, miR-2116-3p, and miR-3187-3p may affect the pathogenesis of HTS by influencing the ubiquitin–proteasome system, which plays a pivotal role in multiple skin fibrosis diseases [32]. However, more direct evidence for complicated mechanisms between the three miRNAs and ubiquitination needs to be found in future research, by means of cell transfection experiments.
In a ceRNA network, a single lncRNA has the potential to simultaneously sponge more than one miRNA. However, this one-to-many pattern has not been investigated in research about the pathogenesis of HTS, and the most of research studying the three miRNAs has only focused on a single miRNA [4]. In this research, a three-miRNA ceRNA network was constructed using the overlapping lncRNAs and target genes of miR-422a, miR-2116-3p, and miR-3187-3p, thereby exploring potential functions that the three miRNAs may perform through regulating the same genes. According to the GO and KEGG functional enrichment analysis, in addition to ubiquitination, this integrated network was predicted to be significantly closed to MAPK signaling pathway, which plays a critical role in HTS pathogenesis. Interestingly, among the three miRNAs, only miR-422a was reported to impact MAPK signaling pathway [19], whose relationship with the other two miRNA is worthy of scholarly attention.
To validate the three-miRNA network in HTS, qRT-PCR was used to determine expression levels of some lncRNAs as well as genes in the network. The results revealed that MAPK1, FOSL2, ABI2, KPNA6, CBL, lncRNA-KCNQ1OT1, and lncRNA-EBLN3P were overexpressed in HTS tissues. MAPK1 was shown to be targeted by miR-422a [19] and a key component of MAPK signaling pathway, which is actively involved in HTS formation. FOSL2, also known as FRA-2, is an activator protein 1 (AP-1) transcription factor and was reported to promote pulmonary fibrosis [33], cardiac fibrosis [34], and skin fibrosis [35], and be regulated by lncRNA UCA1 in ceRNA networks during tumorigenesis [36, 37]. LncRNA KCNQ1OT1 silencing in vivo was confirmed to alleviate cardiac fibrosis by inhibiting the TGF-β1 signaling pathway [38]. These findings suggest that MAPK1, FOSL2, and lncRNA KCNQ1OT1 have the potential to take part in HTS pathogenesis through interacting with miR-422a, miR-2116-3p, and miR-3187-3p. However, there is a shortage of research on the roles of ABI2, KPNA6, CBL, and lncRNA EBLN3P in fibrotic disorders such as HTS; relationships between these genes and the three miRNAs still need further investigation by experiments.
Conclusions
In conclusion, this study showed that miR-422a, miR-2116-3p, and miR-3187-3p are downregulated in HTS tissues and fibroblasts. Furthermore, a distinct ceRNA network for each of the three miRNAs and a three-miRNA (miR-422a, miR-2116-3p, and miR-3187-3p) ceRNA network were proposed and analyzed using bioinformatics analysis, which implied that these networks were significantly related to ubiquitination. Moreover, the expression of the three-miRNA ceRNA network and its correlation with ubiquitination-related genes were validated by qRT-PCR. Overall, the current research suggests that miR-422a, miR-2116-3p, and miR-3187-3p may play a role in HTS formation and correlate with ubiquitination.
Availability of data and materials
The microarray data of miR-422a, miR-2116-3p, and miR-3187-3p in hypertrophic scars are openly available at https://doi.org/10.3892/mmr.2017.6104.
Change history
18 January 2022
A Correction to this paper has been published: https://doi.org/10.1186/s12967-022-03241-w
Abbreviations
- ceRNA:
-
Competitive endogenous RNA
- HTS:
-
Hypertrophic scar
- (qRT‑PCR):
-
Quantitative real‑time PCR
- miRNA:
-
MicroRNA
- lncRNA:
-
Long non-coding RNA
- PPI:
-
Protein–protein network
- GO:
-
Gene Oncology
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- STRING:
-
Search Tool for the Retrieval of Interacting Genes
- MCC:
-
Maximal clique centrality
- EBLN3P:
-
Endogenous bornavirus-like nucleoprotein
- KCNQ1OT1:
-
Potassium voltage-gated channel subfamily Q member 1 opposite strand 1
- IPO9:
-
Importin-9
- ABI2:
-
Abl interactor 2
- ZKSCAN1:
-
Zinc finger protein with KRAB and SCAN domains 1
- FOSL2:
-
Fos-related antigen 2
- GFPT1:
-
Glutamine–fructose-6-phosphate aminotransferase 1
- MAPK1:
-
Mitogen-activated protein kinase 1
- CBL:
-
E3 ubiquitin-protein ligase casitas B lymphoma
- CSNK1G1:
-
Casein kinase I isoform gamma-1
- KPNA6:
-
Importin subunit alpha-7
- EMT:
-
Epithelial to mesenchymal transition
- MEF2D:
-
Myocyte enhancer factor-2D
- PIK3CA:
-
Phosphatidylinositol 4,5-bisphosphate 3-kinase catalytic subunit alpha isoform
- AKT1:
-
RAC-alpha serine/threonine-protein kinase
- SKP2:
-
S-phase kinase-associated protein 2
- PIWIL1:
-
Piwi-like protein 1
- FBXW7:
-
F-box/WD repeat-containing protein 7
- ASXL1:
-
Polycomb group protein
References
Finnerty CC, Jeschke MG, Branski LK, Barret JP, Dziewulski P, Herndon DN. Hypertrophic scarring: the greatest unmet challenge after burn injury. Lancet. 2016;388:1427–36.
Chen L, Li J, Li Q, Yan H, Zhou B, Gao Y, et al. Non-coding RNAs: the new insight on hypertrophic scar. J Cell Biochem. 2017;118:1965–8.
Li L, Han W, Chen Y, Chen Y. MiR-3613-3p inhibits hypertrophic scar formation by down-regulating arginine and glutamate-rich 1. Mol Cell Biochem. 2021;476:1025–36.
Wang R, Zhang S, Chen X, Li N, Li J, Jia R, et al. CircNT5E Acts as a Sponge of miR-422a to Promote Glioblastoma Tumorigenesis. Cancer Res. 2018;78:4812–25.
Chen W, Zhao W, Yang A, Xu A, Wang H, Cong M, et al. Integrated analysis of microRNA and gene expression profiles reveals a functional regulatory module associated with liver fibrosis. Gene. 2017;636:87–95.
Wu M, Wu W, Ding J, Yang J. LINC01433/miR-2116-3p/MYC feedback loop promotes cell proliferation, migration, and the epithelial-mesenchymal transition in breast cancer. Cancer Biother Radiopharm. 2019;34:388–97.
Xiao M, Zou X, Li B, Zhang B. Long non-coding RNA H19 promotes the proliferation, migration and invasion while inhibits apoptosis of hypertrophic scarring fibroblasts by targeting miR-3187-3p/GAB1 axis. Burns. 2021;47:654–64.
Guo L, Xu K, Yan H, Feng H, Wang T, Chai L, et al. MicroRNA expression signature and the therapeutic effect of the microRNA-21 antagomir in hypertrophic scarring. Mol Med Rep. 2017;15:1211–21.
Tu L, Huang Q, Fu S, Liu D. Aberrantly expressed long noncoding RNAs in hypertrophic scar fibroblasts in vitro: a microarray study. Int J Mol Med. 2018;41:1917–30.
Yang Z, Xiao Z, Guo H, Fang X, Liang J, Zhu J, et al. Novel role of the clustered miR-23b-3p and miR-27b-3p in enhanced expression of fibrosis-associated genes by targeting TGFBR3 in atrial fibroblasts. J Cell Mol Med. 2019;23:3246–56.
Zhao Y, Zhang F, Pan Z, Luo H, Liu K, Duan X. LncRNA NR_003923 promotes cell proliferation, migration, fibrosis, and autophagy via the miR-760/miR-215-3p/IL22RA1 axis in human Tenon’s capsule fibroblasts. Cell Death Dis. 2019;10:594.
Xu X, Khoong YM, Gu S, Huang X, Ren J-Y, Gu Y-H, et al. Investigating the potential of LSKL peptide as a novel hypertrophic scar treatment. Biomed Pharmacother. 2020;124:109824.
Amini-Nik S, Yousuf Y, Jeschke MG. Scar management in burn injuries using drug delivery and molecular signaling: current treatments and future directions. Adv Drug Deliv Rev. 2018;123:135–54.
Zhu Z, Hou Q, Li M, Fu X. Molecular mechanism of myofibroblast formation and strategies for clinical drugs treatments in hypertrophic scars. J Cell Physiol. 2020;235:4109–19.
Fang F, Huang R-L, Zheng Y, Liu M, Huo R. Bone marrow derived mesenchymal stem cells inhibit the proliferative and profibrotic phenotype of hypertrophic scar fibroblasts and keloid fibroblasts through paracrine signaling. J Dermatol Sci. 2016;83:95–105.
Liu B, Lin L, Yu S, Xia R, Zheng L. Long non-coding RNA H19 acts as a microRNA-194 sponge to inhibit theapoptosis and promote the proliferation of hypertrophic scar fibroblasts. Can J Physiol Pharmacol. 2021. https://doi.org/10.1139/cjpp-2021-0351.
Zhou Z, Lin Z, He Y, Pang X, Wang Y, Ponnusamy M, et al. The long noncoding RNA D63785 regulates chemotherapy sensitivity in human gastric cancer by targeting miR-422a. Mol Ther Nucleic Acids. 2018;12:405–19.
Liang H, Wang R, Jin Y, Li J, Zhang S. MiR-422a acts as a tumor suppressor in glioblastoma by targeting PIK3CA. Am J Cancer Res. 2016;6:1695–707.
Wei W-T, Nian X-X, Wang S-Y, Jiao H-L, Wang Y-X, Xiao Z-Y, et al. miR-422a inhibits cell proliferation in colorectal cancer by targeting AKT1 and MAPK1. Cancer Cell Int. 2017;17:91.
Zou Y, Chen Y, Yao S, Deng G, Liu D, Yuan X, et al. MiR-422a weakened breast cancer stem cells properties by targeting PLP2. Cancer Biol Ther. 2018;19:436–44.
Shi J, Wang H, Guan H, Shi S, Li Y, Wu X, et al. IL10 inhibits starvation-induced autophagy in hypertrophic scar fibroblasts via cross talk between the IL10-IL10R-STAT3 and IL10-AKT-mTOR pathways. Cell Death Dis. 2016;7:e2133.
Sierra-Sánchez Á, Montero-Vilchez T, Quiñones-Vico MI, Sanchez-Diaz M, Arias-Santiago S. Current Advanced Therapies Based on Human Mesenchymal Stem Cells for Skin Diseases. Frontiers in cell and developmental biology. Front Cell Dev Biol; 2021;9. https://pubmed.ncbi.nlm.nih.gov/33768095/
Zhang Z, Liu C, Chen B, Tang W, Liu Z, Cao W, et al. Smad7 down-regulation via ubiquitin degradation mediated by Smurf2 in fibroblasts of hypertrophic scars in burned patients. Burns. 2021;47:1333–41.
Zhang J, Na S, Pan S, Long S, Xin Y, Jiang Q, et al. Inhibition of USP4 attenuates pathological scarring by downregulation of the TGF-β/Smad signaling pathway. Mol Med Rep. 2019;20:1429–35.
Du S, Wang S, Zhang F, Lv Y. SKP2, positively regulated by circ_ODC1/miR-422a axis, promotes the proliferation of retinoblastoma. J Cell Biochem. 2020;121:322–31.
Inui N, Sakai S, Kitagawa M. Molecular Pathogenesis of pulmonary fibrosis, with focus on pathways related to TGF-β and the ubiquitin-proteasome pathway. Int J Mol Sci. 2021;22:6107.
Jiang H, Liu H, Jiang B. Long non-coding RNA FALEC promotes colorectal cancer progression via regulating miR-2116-3p-targeted PIWIL1. Cancer Biol Ther. 2020;21:1025–32.
Huang H, Yu X, Han X, Hao J, Zhao J, Bebek G, et al. Piwil1 regulates glioma stem cell maintenance and glioblastoma progression. Cell Rep. 2021;34:108522.
Chen X, Li H-D, Bu F-T, Li X-F, Chen Y, Zhu S, et al. Circular RNA circFBXW4 suppresses hepatic fibrosis via targeting the miR-18b-3p/FBXW7 axis. Theranostics. 2020;10:4851–70.
Yang G, Zhang T, Ye J, Yang J, Chen C, Cai S, et al. Circ-ITGA7 sponges miR-3187-3p to upregulate ASXL1, suppressing colorectal cancer proliferation. Cancer Manag Res. 2019;11:6499–509.
Scheuermann JC, de Ayala Alonso AG, Oktaba K, Ly-Hartig N, McGinty RK, Fraterman S, et al. Histone H2A deubiquitinase activity of the Polycomb repressive complex PR-DUB. Nature. 2010;465:243–7.
Shen W, Zhang Z, Ma J, Lu D, Lyu L. The ubiquitin proteasome system and skin fibrosis. Mol Diagn Ther. 2021;25:29–40.
Ucero AC, Bakiri L, Roediger B, Suzuki M, Jimenez M, Mandal P, et al. Fra-2-expressing macrophages promote lung fibrosis in mice. J Clin Invest. 2019;129:3293–309.
Seidenberg J, Stellato M, Hukara A, Ludewig B, Klingel K, Distler O, et al. The AP-1 transcription factor fosl-2 regulates autophagy in cardiac fibroblasts during myocardial fibrogenesis. Int J Mol Sci. 2021;22:1861.
Huang J, Maier C, Zhang Y, Soare A, Dees C, Beyer C, et al. Nintedanib inhibits macrophage activation and ameliorates vascular and fibrotic manifestations in the Fra2 mouse model of systemic sclerosis. Ann Rheum Dis. 2017;76:1941–8.
Chen X, Wang Z, Tong F, Dong X, Wu G, Zhang R. lncRNA UCA1 promotes gefitinib resistance as a ceRNA to target FOSL2 by sponging miR-143 in non-small cell lung cancer. Mol Ther Nucleic Acids. 2020;19:643–53.
Li Z, Niu H, Qin Q, Yang S, Wang Q, Yu C, et al. lncRNA UCA1 mediates resistance to cisplatin by regulating the miR-143/FOSL2-signaling pathway in ovarian cancer. Mol Ther Nucleic Acids. 2019;17:92–101.
Yang F, Qin Y, Lv J, Wang Y, Che H, Chen X, et al. Silencing long non-coding RNA Kcnq1ot1 alleviates pyroptosis and fibrosis in diabetic cardiomyopathy. Cell Death Dis. 2018;9:1000.
Acknowledgements
Zewei Zhang especially thanks Wenhui Liu, Luosha Gu, Qingyan Du, and Hang Zhou for their guidance on bioinformatics, writing, and conceptualization.
Funding
This work was supported by the National Natural Science Foundation of China (Grant Number: 81772086, 82072177) and Key Scientific Research Projects of Colleges and Universities in Henan Province (Grant Number: 20A320033).
Author information
Authors and Affiliations
Contributions
ZWZ, XH, and JHY designed the research, analyzed the data, and contributed equally to this work. ZWZ performed the experiments and wrote the paper. SCG, YXZ, and YHL performed the experimental validation. YMK, SQW, and SYL assisted in the study design and revised the manuscript. TZ and GSL supervised the project and revised the paper. All authorship read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
This study was approved by the Medical Ethics Committee of The First Affiliated Hospital of Zhengzhou University and Shanghai Ninth People’s Hospital, Shanghai Jiaotong University School of Medicine. Written informed consent was acquired from all patients or their guardians.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This article has been updated to correct the funding.
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 http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Zhang, Z., Huang, X., Yang, J. et al. Identification and functional analysis of a three-miRNA ceRNA network in hypertrophic scars. J Transl Med 19, 451 (2021). https://doi.org/10.1186/s12967-021-03091-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12967-021-03091-y