Identification and characterization of bone/cartilage-associated signatures in common fibrotic skin diseases

Background: Fibrotic skin diseases are characterized by excessive accumulation of the extracellular matrix (ECM) and activation of fibroblasts, leading to a global healthcare burden. However, effective treatments of fibrotic skin diseases remain limited, and their pathological mechanisms require further investigation. This study aims to investigate the common biomarkers and therapeutic targets in two major fibrotic skin diseases, namely, keloid and systemic sclerosis (SSc), by bioinformatics analysis.

Methods: The keloid (GSE92566) and SSc (GSE95065) datasets were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified, followed by functional enrichment analysis using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). We then constructed a protein–protein interaction (PPI) network for the identification of hub genes. We explored the possibility of further functional enrichment analysis of hub genes on the Metascape, GeneMANIA, and TissueNexus platforms. Transcription factor (TF)–hub gene and miRNA–hub gene networks were established using NetworkAnalyst. We fixed GSE90051 and GSE76855 as the external validation datasets. Student’s t-test and receiver operating characteristic (ROC) curve were used for candidate hub gene validation. Hub gene expression was assessed in vitro by quantitative real-time PCR.

Results: A total of 157 overlapping DEGs (ODEGs) were retrieved from the GSE92566 and GSE95065 datasets, and five hub genes (COL11A1, COL5A2, ASPN, COL10A1, and COMP) were identified and validated. Functional studies revealed that hub genes were predominantly enriched in bone/cartilage-related and collagen-related processes. FOXC1 and miR-335-5p were predicted to be master regulators at both transcriptional and post-transcriptional levels.

Conclusion: COL11A1, COL5A2, ASPN, COL10A1, and COMP may help understand the pathological mechanism of the major fibrotic skin diseases; moreover, FOXC1 and miR-355-5p could build a regulatory network in keloid and SSc.

Introduction

Fibrosis is characterized by fibroblast proliferation and excessive extracellular matrix deposition in the dermis, which contributes to significant morbidity and mortality (Deng et al., 2021). Fibrotic skin diseases, such as keloid and scleroderma, impose a significant physiological and psychological burden (Nietert et al., 2005; Lv et al., 2020); however, the pathogenesis of the two diseases has to be elucidated yet, and effective therapeutics are currently unavailable (Lee and Pope, 2016; Limandjaja et al., 2020). Previous studies have found similarities and interactions in fibrosis molecular mechanisms, such as the activation of the Notch signaling cascade and TGF-β pathway (Condorelli et al., 2021), shared regulators of miRNA-21 (Li et al., 2017) and miRNA-29 (Deng et al., 2017), and active association between keratinocytes and fibroblasts (Russo et al., 2020). Therefore, screening differentially expressed genes between keloid and SSc could provide an alternative route to identify the mechanisms involved in the fibrosis process.

In recent years, the validation of biomarker expression has been facilitated with the increasing use of tissue microarrays (Hassan et al., 2008). COMP (Agarwal et al., 2013), POSTN (Deng et al., 2021), KGF (Canady et al., 2013), and MMP (Kelly et al., 2010), for example, have been shown to play important roles in ECM deposition and myofibroblast differentiation. The functional role of miRNAs, such as miR-29a-3p (He et al., 2022) and miR-16-5p (Yao et al., 2020), in the progression of keloid and SSc, has been demonstrated through the integration of bioinformatics technology and clinical treatment. CircRNAs have been investigated in other fibrosing processes. Circ_0070963, circRNA_010567, and circ_0044226 have been confirmed as potential therapeutic targets for liver fibrosis (Ji et al., 2020), myocardial fibrosis (Zhou and Yu, 2017), and pulmonary fibrosis (Qi et al., 2020), respectively. However, bioinformatics analysis of fibrotic skin diseases has rarely been performed. In this study, we collected and analyzed RNA sequence data from keloid and SSc, with an aim to shed light on their pathogenesis and provide a potential treatment direction.

Materials and methods

Data collection

The Gene Expression Omnibus ( http://www.ncbi.nlm.nih.gov/geo/) is an open international public database for next-generation functional genomic sequence datasets and high-throughput microarrays (Barrett et al., 2013). We extracted two datasets from GEO, namely, GSE92566, including four keloid and three non-lesional skin samples, and GSE95065, incorporating 18 samples of diffuse cutaneous systemic sclerosis and 15 samples of normal skin. GSE92566 and GSE95065 microarrays were obtained using an Affymetrix Human Genome U133 Plus 2.0 Array and Affymetrix Human Genome U133A 2.0 Array (HGU133A2 Hs ENTREZG 19.0.0), respectively.

Identification of DEGs

GEO2R (http://www.ncbi.nlm.nih.gov/geo/geo2r/) is an interactive tool based on limma algorithm that assists users with GEO data analysis (Davis and Meltzer, 2007). | log (fold change) | > 1 and a p-value < 0.05 were considered statistically significant parameters. Then, the volcano plot was created on a web platform (https://www.bioinformatics.com.cn). A Venn diagram was generated using Rtools 4.2 to reveal ODEGs.

Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses

GO and KEGG analyses were carried out using the Database for Annotation, Visualization and Integrated Discovery (DAVID) tool (https://david.ncifcrf.gov/), which was used for the functional annotation of DEGs and summarization of the relevant biological patterns (Huang et al., 2009). GO consisted of biological processes (BPs), cellular components (CCs), and molecular functions (MFs). A two-tailed p-value < 0.05 was set as the threshold of statistical significance.

Protein–protein interaction network construction

The protein–protein interaction network was predicted by the Search Tool for the Retrieval of Interacting Genes (STRING, version 11.0, https://string-db.org/) with a combined score > 0.4 (medium confidence). The PPI network was then analyzed in Cytoscape software (version 3.8.0), and the Molecular Complex Detection plugin (MCODE, version 2.0.0) was used to find the most significant module with the following parameters: degree cut-off = 2, node score cut-off = 0.2, k-score = 2, and max depth = 100. The top 10 genes were filtered for further analysis on the cytoHubba platform using maximal clique centrality (MCC) and edge percolated component (EPC) algorithms. Heatmaps were created to assess the two-way hierarchical clustering analysis of the hub genes in two datasets using Rtools 4.2.

Functional enrichment analysis and annotation of hub genes

To further explore the biological process (BP) of hub genes, DAVID was used to analyze them again. Metascape is intended to provide comprehensive gene list annotation and deconstruct the molecular mechanisms underlying a biological system (Zhou et al., 2019). The five hub genes were submitted to the Metascape online tool for GO cluster analysis. To perform enrichment analysis and related functional annotations of target genes, genes were divided into clusters using terms with a p-value of 0.01, a minimum count of 3, and an enrichment factor > 1.5. GeneMANIA is used to predict the relative interactive genes of hub genes and to integrate co-expression networks for these genes (Warde-Farley et al., 2010). GeneMANIA’s network consists of physical interactions, co-expression, predicted, co-localization, genetic interactions, pathway, and shared protein domains. TissueNexus is a database of 49 human tissue/cell line functional gene networks that can be used for network search, visualization, and functional analysis (Lin C. X et al., 2022). We further utilized the database to construct a network of functional annotations of the hub genes.

Construction of the regulation network of transcription factors (TFs) and miRNAs

TFs and miRNAs of the hub genes were predicted via NetworkAnalyst. The relevant TFs were obtained from the JASPAR database, and the miRNAs were acquired from miRTarBase v8.0 and TarBase v8.0 (Zheng et al., 2022). Then, the network of TFs and miRNAs was, respectively, constructed using Cytoscape 3.8.2 software.

Tissue samples

After the patients signed the informed consent, three cases of human normal tissues, keloid tissues, and SSc tissues were obtained from the Department of Dermatology, Union Hospital, Tongji Medical College (Wuhan, China). Then, they were stored at 4°C for total RNA isolation. The clinical information is shown in Table 1.

TABLE 1

TABLE 1. Clinical information of validated samples.

Validation of the performance of hub genes

The expression of the hub genes in keloid, SSc, and normal tissues was validated using Student’s t-test in GraphPad Prism 8.0.1, and their specificity and sensitivity were further evaluated by the receiver operating characteristic curve analysis.

The total RNA was isolated using TRIzol (Vazyme, China) according to the manufacturer’s instructions. Approximately 1 mg of the total RNA from each sample was reverse-transcribed into cDNA using the PrimeScriptTM RT Master Mix kit (Takara, RR036A). Real-time PCR was performed using the TB Green Premix Ex TaqTM kit (Takara, RR420A) according to the manufacturer’s instructions. The primers for human were as follows: COMP: forward 5′-CAT​CAG​GAC​TCT​CGG​GAC​AAC​T-3′ and reverse 5′-GTC​TAC​CAC​CTT​GTC​TGC​ATC​AA-3′; ASPN: forward 5′- AGT​CCC​AAC​CAA​CAT​TCC​ATT​TG-3′ and reverse 5′- TGG​GTG​AAT​CTT​CGT​TAG​CTT​GT-3′; COL5A2: forward 5′-GAA​GCC​TCC​CAG​AAC​ATC​ACT​TA-3′ and reverse 5′-CCC​ACA​TTT​CCA​TTC​CGC​TTA​GA-3′; COL10A1: forward 5′-GCA​TAA​AAG​GCC​CAC​TAC​CCA-3′ and reverse 5′- GCA​TAA​AAG​GCC​CAC​TAC​CCA-3′; COL11A1: forward 5′- AAG​ACG​GAG​ACA​AGG​GTG​AAA​TT-3′ and reverse 5′- TTC​ACC​ATC​ACC​TCC​AGC​AAT​TC-3′; and GAPDH: forward 5′-GGA​GTC​CAC​TGG​CGT​CTT​CA-3′ and reverse 5′-GTC​ATG​AGT​CCT​TCC​ACG​ATA​CC-3′. The relative mRNA expression level was standardized to GAPDH and was calculated using the 2^−ΔΔCT method. Statistically significant differences are indicated by *— *p < 0.05, **p < 0.01, and ***p < 0.001.

In addition, GSE90051 was used to validate the expression of hub genes in keloids and adjacent normal tissues, while GSE76885 was used to verify the expression pattern of hub genes in SSc and normal tissues.

Results

DEG identification

After standardization, a total of 889 DEGs (163 upregulated and 726 downregulated) with SSc in GSE95065 and 2,509 DEGs (1,020 upregulated and 1,489 downregulated) with keloids in GSE92566 were identified (Figures 1A, B). There were 157 ODEGs between GSE95065 and GSE92566, as shown in the Venn diagram (Figure 1C), including 56 consistently upregulated DEGs and 77 downregulated DEGs.

FIGURE 1

FIGURE 1. (A, B) Volcano plots of the DEGs in GSE95065 and GSE92566 datasets. Red and blue represent the upregulated and downregulated genes, respectively. (C) Venn diagram of the ODEGs in the two datasets.

GO and KEGG enrichment analyses of ODEGs

Using the p-value < 0.05 criterion, we identified the significant GO enrichment of ODEGs (Figure 2A). Among the upregulated ODEGs, BPs were highly associated with collagen fibril organization, immune response, and extracellular matrix organization. The CC was involved in the extracellular region, extracellular space, and extracellular matrix. Top MF was significantly enriched in extracellular matrix structural constituents, collagen binding, and metalloendopeptidase activity. Among the downregulated ODEGs, BP was mostly enriched in muscle contraction, positive regulation of protein kinase B signaling, and long-chain fatty acid import; CC was mostly enriched in extracellular exosome, extracellular region, and sarcolemma; and MF was mostly enriched in serine-type peptidase activity, serine-type endopeptidase activity, and protein binding. Furthermore, KEGG pathway analysis (Figure 2B) revealed that the upregulated ODEGs were involved in protein digestion and absorption and the downregulated ODEGs in hypertrophic cardiomyopathy, dilated cardiomyopathy, and arrhythmogenic right ventricular cardiomyopathy.

FIGURE 2

FIGURE 2. (A) GO analysis of the ODEGs. (B) KEGG pathway enrichment analysis of the ODEGs.

Protein–protein interaction network construction and hub gene analysis

A total of 133 consistently changed ODEGs were placed in STRING to explore the interactive networks with free proteins removed. Cytoscape was used to visualize PPI, which consisted of 67 nodes and 94 edges (Figure 3A). Based on maximal clique centrality (MCC) and edge percolated component (EPC) algorithms, the top 10 genes were identified using cytoHubba (Table 2). COLLA1, COL5A2, ASPN, COL10A1, EGF, COMP, and SYP were identified as key genes. Then, MCODE algorithm was utilized to reveal the most significant modules, as shown in Figures 3B, C. Module1 (Figure 3B) contained five upregulated DEGs (COL11A1, COL5A2, COL10A1, ASPN, and COMP); module 2 (Figure 3C) consisted of four downregulated DEGs (SGCG, SGCA, SNTB1, and DTNB). As a result of the overlap of cytoHubba and MCODE, COMP, ASPN, COL10A1, COL11A1, and COL5A2 were identified as hub genes. Then, heatmaps (Figures 3D, E) were constructed, which revealed a high degree of clustering of these hub genes.

FIGURE 3

FIGURE 3. (A) PPI networks of ODEGs constructed by Cytoscape. (B–C) Main modules obtained using MCODE. (D–E) Heatmaps of the hub genes in GSE92566 and GSE95065 datasets. Red and blue represent upregulation and downregulation of expression, respectively.

TABLE 2

TABLE 2. Top 10 genes in the network using MCC and EPC algorithms.

Enrichment analysis and functional annotation of hub genes

The five hub genes were entered into DAVID for GO analysis. BPs of these genes were most prominent in collagen fibril organization, skeletal system development, and tendon or chondrocyte development (Figure 4A).

FIGURE 4

FIGURE 4. (A) Biological processes enriched by the five hub genes. (B) Pathway and process enrichment analysis on the Metascape platform. (C) Enrichment analysis in DisGeNET. (D) Interaction network of the hub genes and their predicted genes via GeneMANIA. (E) Functional annotation of hub genes using TissueNexus.

The Metascape platform summarized the functional enrichment analysis of the hub genes. The bar graph (Figure 4B) showed that these genes mainly enriched in naba core matrisome, naba collagens, and ossification. Interestingly, DisGeNET (Figure 4C) also revealed that the hub genes were primarily enriched in bone/cartilage-related diseases such as osteoarthritis of the hip, irregular vertebral endplates, spondyloepimetaphyseal disorder, and chondrodysplasia. It could shed light on the possible pathogenesis of fibrotic skin disorders.

GeneMANIA was used to create a protein interaction network of the five hub genes and the 20 predicted genes (Figure 4D). Physical interactions (77.64%), co-expression (8.01%), predicted (5.37%), co-localization (3.63%), genetic interactions (2.87%), pathways (1.88%), and shared protein domains (0.60%) are among the networks based on GeneMANIA’s functional annotation patterns. The function of predicted genes was mostly associated with collagen formation (Table 3). It is worth noting that there were proteins related to skeletal system development (six of the predicted genes) and regulation of cartilage development (three of the predicted genes).

TABLE 3

TABLE 3. Gene functions analyzed using GeneMANIA.

TissueNexus was used to create functional gene networks. COL5A2 was discovered to be the most important player in the network (Figure 4E). Furthermore, the top three functional annotations of the hub genes were extracellular matrix organization, extracellular structure, and external encapsulating structure (Table 4). The functional annotation further demonstrated that bone/cartilage-related genes participated in the process of extracellular matrix deposition.

TABLE 4

TABLE 4. Functional annotation of hub genes using TissueNexus.

Analyses of the interactions of TF–hub genes and miRNA–hub genes

The TFs of hub genes were acquired using the JASPAR database on the NetworkAnalyst web tool. The miRNAs of the hub genes were obtained using miRTarBase v8.0 and the TarBase v8.0 via NetworkAnalyst. The regulatory network between five seeds (COMP, ASPN, COL5A2, COL10A1, and COL11A1) and their TFs was established using Cytoscape (Figure 5A). The TF–hub genes had 28 nodes and 31 edges. FOXC1 may play an important role in the network because it has the potential to regulate COL5A2, COL11A1, ASPN, and COL10A1 simultaneously. The miRNA–hub gene regulatory network had 75 nodes and 92 edges (Figure 5B). Has-miR-335-5p was most likely targeted by COL5A2, COL11A1, ASPN, and COL10A1.

FIGURE 5

FIGURE 5. (A) TF–hub gene network established using the JASPAR database on NetworkAnalyst. (B) miRNA–hub gene network acquired utilizing miRTarBase v8.0 and TarBase v8.0 via NetworkAnalyst.

Validation of hub genes

The expression of five hub genes was higher in group keloid and SSc tissues than in normal tissues (p < 0.05). Furthermore, the AUC values of the ROC curves (Figure 6A) for the five genes were all greater than 0.8, indicating high sensitivity and specificity. The hub genes (COMP, ASPN, COL10A1, COL11A1, and COL5A2) were found to be reliable biomarkers with high diagnostic accuracy (Supplementary Table S1). We validated the expression of hub genes by qRT-PCR in clinical samples (Figure 6B). Compared to normal tissues, these genes were significantly highly expressed in keloid and SSc. In addition, the hub genes were then validated using external datasets. Notably, GSE90051 described familiar changes in the hub genes, which were more abundant in keloid tissues than in normal tissues (Figure 6C). As for the dataset GSE76855, we only disclosed the expression of three hub genes (COMP, ASPN, and COL5A2). Thankfully, COMP, ASPN, and COL5A2 were all upregulated in SSc tissues compared to normal tissues, supporting our findings (Figure 6D).

FIGURE 6

FIGURE 6. Validation of hub genes. (A) ROC curves demonstrating the diagnostic value of the hub genes. (B) Expression level of hub genes validated by qPCR. (C) Hub gene validation in the GSE90051 dataset. (D) Hub gene validation in the GSE76855 dataset.

Discussion

Fibrotic skin diseases including keloid, hypertrophic scar, systemic sclerosis, and lichen sclerosus et atrophicus manifest analogous histopathology but unclear etiopathogenesis. Considering the reported relationships between keloid and SSc, screening biomarkers associated with fibrosis is a promising way to gain insights into common fibrotic skin diseases at a genomic level. To the best of our knowledge, no similar study has been reported on the common biomarkers in keloid and SSc by bioinformatics analysis, which could contribute to discovering underlying mechanisms of cutaneous fibrosis and novel therapeutic targets.

A total of 157 ODEGs obtained from keloid (GSE92566) and SSc (GSE95065) datasets were found to be related to cardiomyopathy pathways by KEGG and extracellular components/processes by GO analysis, respectively. A total of five hub genes (COMP, ASPN, COL5A2, COL10A1, and COL11A1) were identified by MCC, EPC, and MCODE algorithms relying on STRING and Cytoscape platforms. Surprisingly, the filtered hub genes were also closely associated with the excess production of bone and cartilage proteins in the extracellular matrix of fibrotic skin diseases, according to the functional enrichments.

COMP encodes cartilage oligomeric matrix protein, an extracellular matrix glycoprotein found in cartilage, tendons, ligaments, and the growth plate (Hedbom et al., 1992). Mutations in the COMP genes cause pseudoachondroplasia (PSACH) and multiple epiphyseal dysplasia (MED/EDM1) (Briggs et al., 1995). COMP has been reported to affect the secretion of type IX collagen and matrilin-3 (Hecht et al., 2005) and to bind to collagen types I and II (Södersten et al., 2005) and fibronectin (Di Cesare et al., 2002), thereby maintaining an ECM network. In clinical practice, COMP is used as a prognostic marker for joint injury, a biomarker for idiopathic pulmonary fibrosis, and a biomarker for cartilage degeneration-associated osteoarthritis and rheumatoid arthritis (Posey et al., 2018). Recently, pathogenesis of various fibrotic diseases, such as liver fibrosis (Magdaleno et al., 2016) and renal fibrosis (Kim et al., 2006), has been linked to COMP deposition. In line with previous research (Farina et al., 2006; Inui et al., 2011; Agarwal et al., 2013), our study also supported that COMP overexpression is a common feature of fibrotic skin diseases.

ASPN, a member of the small leucine-rich proteoglycan (SLRP) family of extracellular proteins, was initially discovered in human cartilage (Lorenzo et al., 2001). Copy number loss in the ASPN gene region may be involved in the pathogenesis of acetabular dysplasia (Sekimoto et al., 2017). ASPN has been identified as a novel mesenchymal stromal cell factor that affects tumor microenvironment and metastatic progression (Hughes et al., 2019), including prostate cancer (Hurley et al., 2016), breast cancer (Maris et al., 2015), colorectal cancer (Huang C. Z et al., 2022), and others. It also played a crucial biological role in the process of tissue fibrosis (Liu et al., 2021; Huang S. et al., 2022; Huang C. et al., 2022). Huang S. et al. (2022) reported that ASPN could promote lung myofibroblast differentiation by facilitating TGF-β/Smad signaling. Huang C. et al. (2022) unveiled that ASPN could attenuate fibrosis in cardiac remodeling by regulating mitochondrial bioenergetics and protecting cardiomyocytes from hypoxia–reoxygenation-mediated cell death. Liu et al. (2021) showed that ASPN participated in keloid progression by disturbing interfibroblast mechanocommunication and inhibiting matrix remodeling.

COL5A2, COL10A1, and COL11A1 are members of the collagen gene family that are highly expressed in bone/cartilage-associated diseases (Gu et al., 2014; Yang et al., 2018; Sun et al., 2020). In the early stage of osteoblast differentiation, COL5A2 may help with bone material (Tawa et al., 2021) and actin skeleton formation (Hong et al., 2010). Mutations in COL10A1, which is a specific marker of hypertrophic chondrocytes (Gu et al., 2014), result in Schmid metaphyseal chondrodysplasia (SMCD), an autosomal dominantly inherited skeletal disorder (Bateman et al., 2005). COL11A1 is a fibril-forming collagen that regulates cartilage and tendon fibrillar structure (Sun et al., 2020). Stickler syndrome type II is caused by COL11A1 gene mutations (Vogiatzi et al., 2018). Previous bioinformatics analysis has demonstrated that COL5A2, COL10A1, and COL11A1 may be used as a diagnostic or prognostic biomarker for tumors such as gastric cancer (Niu et al., 2022), breast carcinoma (Giussani et al., 2018), and non-small-cell lung cancer (Andriani et al., 2018). According to the findings of RNA sequencing, COL5A2, COL10A1, and COL11A1 have been reported to be upregulated in keloid (Bi et al., 2021; Lin P. et al., 2022), as in our study. A whole-exome sequencing study revealed that COL5A2 is relevant to the fibrotic features of diffuse cutaneous SSc (Mak et al., 2016). However, the relationship between these genes and SSc is poorly understood.

In recent years, we have considered that fibroblasts differ in their heterogeneity depending on their anatomical location (Lynch and Watt, 2018). Human skin single-cell transcriptomes indicated that fibroblasts are classified into four subpopulations: secretory-papillary, secretory-reticular, mesenchymal, and pro-inflammatory (Solé-Boldo et al., 2020). Interestingly, the upregulated genes like COMP, COL5A2, and COL11A1 were enriched in mesenchymal fibroblasts reported by single-cell RNA-seq (Deng et al., 2021). Thus, we hypothesized that mesenchymal fibroblasts secrete bone and cartilage proteins, which contribute to fibrosis.

In our study, the transcription factor (FOXC1) and miRNA (miR-335-5p) were predicted to be the key transcriptional and post-transcriptional regulators of DEGs. FOXC1 has been linked to cartilage differentiation (Xu et al., 2021), vascular development (Shu et al., 2022), tumor development, and metastasis (Zhang et al., 2020; Lin Z. et al., 2021). High miR-335-5p expression promotes bone formation and regeneration (Zhang et al., 2017) and angiogenesis (Walz et al., 2019) and regulates cancer progression by targeting MAPK10 (Gao et al., 2021) or NLRP1/7 (Lin W. et al., 2021). Furthermore, a FOXC1 and miR-335-5p network was presented in non-small-cell lung cancer (Mosharaf et al., 2022), type 2 diabetes (Rahman et al., 2020), and vascular dementia (Shu et al., 2022). The dermis of keloid and SSc is known to be made up of newly formed vessels (Trojanowska, 2010; Ogawa, 2017). The FOXC1 and miR-335-5p networks may collaborate in the keloid and SSc processes. This is the first study to look into the relationship between FOXC1 and miR-335-5p in fibrotic skin diseases.

In summary, functional enrichments showed that bone/cartilage-associated genes are overexpressed in keloid and SSc. After validation, the five hub genes can be identified as diagnostic biomarkers of fibrotic skin diseases.

Conclusion

In our present study, we screened bone/cartilage-related genes (COMP, ASPN, COL5A2, COL10A1, and COL11A1), transcription factor (FOXC1), and miRNA (miR-335-5p) in major fibrotic skin diseases including keloid and SSc. We hope to provide further group work and insights into fibrotic skin diseases and find potential therapeutic targets for clinical treatments.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding authors.

Author contributions

TW: designed the study, collected the data, and draft writing; YJ: designed the study, collected the data, and draft writing; FC: designed the study, collected the data, and draft writing; XX: collected the data; JC: collected the data; YL: collected the data; YW: collected the data; JZ: collected the data; MZ: draft writing and revising; CH: draft writing and revising.

Funding

This work was partially supported by the National Natural Science Foundation of China (Nos 81972565 and 82003324).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1121728/full#supplementary-material

References

Agarwal, P., Schulz, J-N., Blumbach, K., Andreasson, K., Heinegård, D., Paulsson, M., et al. (2013). Enhanced deposition of cartilage oligomeric matrix protein is a common feature in fibrotic skin pathologies. Matrix Biol. 32 (6), 325–331. doi:10.1016/j.matbio.2013.02.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Andriani, F., Landoni, E., Mensah, M., Facchinetti, F., Miceli, R., Tagliabue, E., et al. (2018). Diagnostic role of circulating extracellular matrix-related proteins in non-small cell lung cancer. BMC Cancer 18 (1), 899. doi:10.1186/s12885-018-4772-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrett, T., Wilhite, S. E., Ledoux, P., Evangelista, C., Kim, I. F., Tomashevsky, M., et al. (2013). NCBI geo: Archive for functional genomics data sets-update. Nucleic Acids Res. 41, D991–D995. doi:10.1093/nar/gks1193

PubMed Abstract | CrossRef Full Text | Google Scholar

Bateman, J. F., Wilson, R., Freddi, S., Lamandé, S. R., and Savarirayan, R. (2005). Mutations of COL10A1 in Schmid metaphyseal chondrodysplasia. Hum. Mutat. 25 (6), 525–534. doi:10.1002/humu.20183

PubMed Abstract | CrossRef Full Text | Google Scholar

Bi, S., Liu, R., Wu, B., He, L., and Gu, J. (2021). Bioinformatic analysis of key genes and pathways related to keloids. Biomed. Res. Int. 2021, 5897907. doi:10.1155/2021/5897907

PubMed Abstract | CrossRef Full Text | Google Scholar

Briggs, M. D., Hoffman, S. M., King, L. M., Olsen, A. S., Mohrenweiser, H., Leroy, J. G., et al. (1995). Pseudoachondroplasia and multiple epiphyseal dysplasia due to mutations in the cartilage oligomeric matrix protein gene. Nat. Genet. 10 (3), 330–336. doi:10.1038/ng0795-330

PubMed Abstract | CrossRef Full Text | Google Scholar

Canady, J., Arndt, S., Karrer, S., and Bosserhoff, A. K. (2013). Increased KGF expression promotes fibroblast activation in a double paracrine manner resulting in cutaneous fibrosis. J. Investigative Dermatology 133 (3), 647–657. doi:10.1038/jid.2012.389

CrossRef Full Text | Google Scholar

Condorelli, A. G., El Hachem, M., Zambruno, G., Nystrom, A., Candi, E., and Castiglia, D. (2021). Notch-ing up knowledge on molecular mechanisms of skin fibrosis: Focus on the multifaceted Notch signalling pathway. J. Biomed. Sci. 28 (1), 36. doi:10.1186/s12929-021-00732-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Davis, S., and Meltzer, P. S. (2007). GEOquery: A bridge between the gene expression Omnibus (GEO) and BioConductor. Bioinformatics 23 (14), 1846–1847. doi:10.1093/bioinformatics/btm254

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, C-C., Hu, Y-F., Zhu, D-H., Cheng, Q., Gu, J-J., Feng, Q-L., et al. (2021). Single-cell RNA-seq reveals fibroblast heterogeneity and increased mesenchymal fibroblasts in human fibrotic skin diseases. Nat. Commun. 12 (1), 3709. doi:10.1038/s41467-021-24110-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, Z., He, Y., Yang, X., Shi, H., Shi, A., Lu, L., et al. (2017). MicroRNA-29: A crucial player in fibrotic disease. Mol. Diagn Ther. 21 (3), 285–294. doi:10.1007/s40291-016-0253-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Di Cesare, P. E., Chen, F. S., Moergelin, M., Carlson, C. S., Leslie, M. P., Perris, R., et al. (2002). Matrix-matrix interaction of cartilage oligomeric matrix protein and fibronectin. Matrix Biol. 21 (5), 461–470. doi:10.1016/s0945-053x(02)00015-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Farina, G., Lemaire, R., Korn, J. H., and Widom, R. L. (2006). Cartilage oligomeric matrix protein is overexpressed by scleroderma dermal fibroblasts. Matrix Biol. 25 (4), 213–222. doi:10.1016/j.matbio.2006.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, Y., Wang, Y., Wang, X., Zhao, C., Wang, F., Du, J., et al. (2021). miR-335-5p suppresses gastric cancer progression by targeting MAPK10. Cancer Cell. Int. 21 (1), 71. doi:10.1186/s12935-020-01684-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Giussani, M., Landoni, E., Merlino, G., Turdo, F., Veneroni, S., Paolini, B., et al. (2018). Extracellular matrix proteins as diagnostic markers of breast carcinoma. J. Cell. Physiol. 233 (8), 6280–6290. doi:10.1002/jcp.26513

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, J., Lu, Y., Li, F., Qiao, L., Wang, Q., Li, N., et al. (2014). Identification and characterization of the novel Col10a1 regulatory mechanism during chondrocyte hypertrophic differentiation. Cell. Death Dis. 5, e1469. doi:10.1038/cddis.2014.444

PubMed Abstract | CrossRef Full Text | Google Scholar

Hassan, S., Ferrario, C., Mamo, A., and Basik, M. (2008). Tissue microarrays: Emerging standard for biomarker validation. Curr. Opin. Biotechnol. 19 (1), 19–25. doi:10.1016/j.copbio.2007.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

He, Y., Zhang, Z., Yin, B., Li, S., Wang, P., Lan, J., et al. (2022). Identifying miRNAs associated with the progression of keloid through mRNA-miRNA network analysis and validating the targets of miR-29a-3p in keloid fibroblasts. Biomed. Res. Int. 2022, 6487989. doi:10.1155/2022/6487989

PubMed Abstract | CrossRef Full Text | Google Scholar

Hecht, J. T., Hayes, E., Haynes, R., and Cole, W. G. (2005). COMP mutations, chondrocyte function and cartilage matrix. Matrix Biol. 23 (8), 525–533. doi:10.1016/j.matbio.2004.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Hedbom, E., Antonsson, P., Hjerpe, A., Aeschlimann, D., Paulsson, M., Rosa-Pimentel, E., et al. (1992). Cartilage matrix proteins. An acidic oligomeric protein (COMP) detected only in cartilage. J. Biol. Chem. 267 (9), 6132–6136. doi:10.1016/s0021-9258(18)42671-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Hong, D., Chen, H-X., Yu, H-Q., Liang, Y., Wang, C., Lian, Q-Q., et al. (2010). Morphological and proteomic analysis of early stage of osteoblast differentiation in osteoblastic progenitor cells. Exp. Cell. Res. 316 (14), 2291–2300. doi:10.1016/j.yexcr.2010.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, C. Z., Zhou, Y., Tong, Q-S., Duan, Q-J., Zhang, Q., Du, J-Z., et al. (2022a). Precision medicine-guided co-delivery of ASPN siRNA and oxaliplatin by nanoparticles to overcome chemoresistance of colorectal cancer. Biomaterials 290, 121827. doi:10.1016/j.biomaterials.2022.121827

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, S., Lai, X., Yang, L., Ye, F., Huang, C., Qiu, Y., et al. (2022b). Asporin promotes TGF-β-induced lung myofibroblast differentiation by facilitating rab11-dependent recycling of TβRI. Am. J. Respir. Cell. Mol. Biol. 66 (2), 158–170. doi:10.1165/rcmb.2021-0257OC

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, C., Sharma, A., Thakur, R., Rai, D., Katiki, M., Germano, JdF., et al. (2022c). Asporin, an extracellular matrix protein, is a beneficial regulator of cardiac remodeling. Matrix Biol. 110, 40–59. doi:10.1016/j.matbio.2022.04.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, D. W., Sherman, B. T., and Lempicki, R. A. (2009). Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 37 (1), 1–13. doi:10.1093/nar/gkn923

PubMed Abstract | CrossRef Full Text | Google Scholar

Hughes, R. M., Simons, B. W., Khan, H., Miller, R., Kugler, V., Torquato, S., et al. (2019). Asporin restricts mesenchymal stromal cell differentiation, alters the tumor microenvironment, and drives metastatic progression. Cancer Res. 79 (14), 3636–3650. doi:10.1158/0008-5472.CAN-18-2931

PubMed Abstract | CrossRef Full Text | Google Scholar

Hurley, P. J., Sundi, D., Shinder, B., Simons, B. W., Hughes, R. M., Miller, R. M., et al. (2016). Germline variants in asporin vary by race, modulate the tumor microenvironment, and are differentially associated with metastatic prostate cancer. Clin. Cancer Res. 22 (2), 448–458. doi:10.1158/1078-0432.CCR-15-0256

PubMed Abstract | CrossRef Full Text | Google Scholar

Inui, S., Shono, F., Nakajima, T., Hosokawa, K., and Itami, S. (2011). Identification and characterization of cartilage oligomeric matrix protein as a novel pathogenic factor in keloids. Am. J. Pathol. 179 (4), 1951–1960. doi:10.1016/j.ajpath.2011.06.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, D., Chen, G-F., Wang, J-C., Ji, S-H., Wu, X-W., Lu, X-J., et al. (2020). Hsa_circ_0070963 inhibits liver fibrosis via regulation of miR-223-3p and LEMD3. Aging (Albany NY) 12 (2), 1643–1655. doi:10.18632/aging.102705

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelly, B. C., Markle, L. S., Vickers, J. L., Petitt, M. S., Raimer, S. S., and McNeese, C. (2010). The imbalanced expression of matrix metalloproteinases in nephrogenic systemic fibrosis. J. Am. Acad. Dermatology 63 (3), 483–489. doi:10.1016/j.jaad.2009.09.006

CrossRef Full Text | Google Scholar

Kim, W., Moon, S-O., Lee, S. Y., Jang, K. Y., Cho, C-H., Koh, G. Y., et al. (2006). COMP-angiopoietin-1 ameliorates renal fibrosis in a unilateral ureteral obstruction model. J. Am. Soc. Nephrol. 17 (9), 2474–2483. doi:10.1681/ASN.2006020109

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J. J., and Pope, J. E. (2016). Diagnosis and management of systemic sclerosis: A practical approach. Drugs 76 (2), 203–213. doi:10.1007/s40265-015-0491-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Zhang, J., Lei, Y., Lyu, L., Zuo, R., and Chen, T. (2017). MicroRNA-21 in skin fibrosis: Potential for diagnosis and treatment. Mol. Diagn Ther. 21 (6), 633–642. doi:10.1007/s40291-017-0294-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Limandjaja, G. C., Niessen, F. B., Scheper, R. J., and Gibbs, S. (2020). The keloid disorder: Heterogeneity, histopathology, mechanisms and models. Front. Cell. Dev. Biol. 8, 360. doi:10.3389/fcell.2020.00360

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, C. X., Li, H-D., Deng, C., Guan, Y., and Wang, J. (2022a). TissueNexus: A database of human tissue functional gene networks built with a large compendium of curated RNA-seq data. Nucleic Acids Res. 50 (D1), D710–D718. doi:10.1093/nar/gkab1133

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, P., Zhang, G., Peng, R., Zhao, M., and Li, H. (2022b). Increased expression of bone/cartilage-associated genes and core transcription factors in keloids by RNA sequencing. Exp. Dermatol 31 (10), 1586–1596. doi:10.1111/exd.14630

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Z., Huang, W., He, Q., Li, D., Wang, Z., Feng, Y., et al. (2021a). FOXC1 promotes HCC proliferation and metastasis by Upregulating DNMT3B to induce DNA Hypermethylation of CTH promoter. J. Exp. Clin. Cancer Res. 40 (1), 50. doi:10.1186/s13046-021-01829-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, W., Chen, Y., Wu, B., Chen, Y., and Li, Z. (2021b). Identification of the pyroptosis-related prognostic gene signature and the associated regulation axis in lung adenocarcinoma. Cell. Death Discov. 7 (1), 161. doi:10.1038/s41420-021-00557-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L., Yu, H., Long, Y., You, Z., Ogawa, R., Du, Y., et al. (2021). Asporin inhibits collagen matrix-mediated intercellular mechanocommunications between fibroblasts during keloid progression. FASEB J. 35 (7), e21705. doi:10.1096/fj.202100111R

PubMed Abstract | CrossRef Full Text | Google Scholar

Lorenzo, P., Aspberg, A., Onnerfjord, P., Bayliss, M. T., Neame, P. J., and Heinegard, D. (2001). Identification and characterization of asporin. a novel member of the leucine-rich repeat protein family closely related to decorin and biglycan. J. Biol. Chem. 276 (15), 12201–12211. doi:10.1074/jbc.M010932200

PubMed Abstract | CrossRef Full Text | Google Scholar

Lv, W., Ren, Y., Hou, K., Hu, W., Yi, Y., Xiong, M., et al. (2020). Epigenetic modification mechanisms involved in keloid: Current status and prospect. Clin. Epigenetics 12 (1), 183. doi:10.1186/s13148-020-00981-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Lynch, M. D., and Watt, F. M. (2018). Fibroblast heterogeneity: Implications for human disease. J. Clin. Invest. 128 (1), 26–35. doi:10.1172/JCI93555

PubMed Abstract | CrossRef Full Text | Google Scholar

Magdaleno, F., Arriazu, E., Ruiz de Galarreta, M., Chen, Y., Ge, X., Conde de la Rosa, L., et al. (2016). Cartilage oligomeric matrix protein participates in the pathogenesis of liver fibrosis. J. Hepatol. 65 (5), 963–971. doi:10.1016/j.jhep.2016.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Mak, A. C. Y., Tang, P. L. F., Cleveland, C., Smith, M. H., Kari Connolly, M., Katsumoto, T. R., et al. (2016). Brief report: Whole-exome sequencing for identification of potential causal variants for diffuse cutaneous systemic sclerosis. Arthritis Rheumatol. 68 (9), 2257–2262. doi:10.1002/art.39721

PubMed Abstract | CrossRef Full Text | Google Scholar

Maris, P., Blomme, A., Palacios, A. P., Costanza, B., Bellahcène, A., Bianchi, E., et al. (2015). Asporin is a fibroblast-derived TGF-β1 inhibitor and a tumor suppressor associated with good prognosis in breast cancer. PLoS Med. 12 (9), e1001871. doi:10.1371/journal.pmed.1001871

PubMed Abstract | CrossRef Full Text | Google Scholar

Mosharaf, M. P., Reza, M. S., Gov, E., Mahumud, R. A., and Mollah, M. N. H. (2022). Disclosing potential key genes, therapeutic targets and agents for non-small cell lung cancer: Evidence from integrative bioinformatics analysis. Vaccines (Basel). 10 (5), 771. doi:10.3390/vaccines10050771

PubMed Abstract | CrossRef Full Text | Google Scholar

Nietert, P. J., Mitchell, H. C., Bolster, M. B., Curran, M. Y., Tilley, B. C., and Silver, R. M. (2005). Correlates of depression, including overall and gastrointestinal functional status, among patients with systemic sclerosis. J. Rheumatol. 32 (1), 51–57.

PubMed Abstract | Google Scholar

Niu, X., Ren, L., Hu, A., Zhang, S., and Qi, H. (2022). Identification of potential diagnostic and prognostic biomarkers for gastric cancer based on bioinformatic analysis. Front. Genet. 13, 862105. doi:10.3389/fgene.2022.862105

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogawa, R. (2017). Keloid and hypertrophic scars are the result of chronic inflammation in the reticular dermis. Int. J. Mol. Sci. 18 (3), 606. doi:10.3390/ijms18030606

PubMed Abstract | CrossRef Full Text | Google Scholar

Posey, K. L., Coustry, F., and Hecht, J. T. (2018). Cartilage oligomeric matrix protein: COMPopathies and beyond. Matrix Biol. 71-72, 161–173. doi:10.1016/j.matbio.2018.02.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, F., Li, Y., Yang, X., Wu, Y., Lin, L., and Liu, X. (2020). Hsa_circ_0044226 knockdown attenuates progression of pulmonary fibrosis by inhibiting CDC27. Aging (Albany NY) 12 (14), 14808–14818. doi:10.18632/aging.103543

PubMed Abstract | CrossRef Full Text | Google Scholar

Rahman, M. H., Peng, S., Hu, X., Chen, C., Rahman, M. R., Uddin, S., et al. (2020). A network-based bioinformatics approach to identify molecular biomarkers for type 2 diabetes that are linked to the progression of neurological diseases. Int. J. Environ. Res. Public Health 17 (3), 1035. doi:10.3390/ijerph17031035

PubMed Abstract | CrossRef Full Text | Google Scholar

Russo, B., Brembilla, N. C., and Chizzolini, C. (2020). Interplay between keratinocytes and fibroblasts: A systematic review providing a new angle for understanding skin fibrotic disorders. Front. Immunol. 11, 648. doi:10.3389/fimmu.2020.00648

PubMed Abstract | CrossRef Full Text | Google Scholar

Sekimoto, T., Ishii, M., Emi, M., Kurogi, S., Funamoto, T., Yonezawa, Y., et al. (2017). Copy number loss in the region of the ASPN gene in patients with acetabular dysplasia: ASPN CNV in acetabular dysplasia. Bone Jt. Res. 6 (7), 439–445. doi:10.1302/2046-3758.67.BJR-2016-0094.R1

CrossRef Full Text | Google Scholar

Shu, J., Wei, W., and Zhang, L. (2022). Identification of molecular signatures and candidate drugs in vascular dementia by bioinformatics analyses. Front. Mol. Neurosci. 15, 751044. doi:10.3389/fnmol.2022.751044

PubMed Abstract | CrossRef Full Text | Google Scholar

Södersten, F., Ekman, S., Eloranta, M-L., Heinegård, D., Dudhia, J., and Hultenby, K. (2005). Ultrastructural immunolocalization of cartilage oligomeric matrix protein (COMP) in relation to collagen fibrils in the equine tendon. Matrix Biol. 24 (5), 376–385. doi:10.1016/j.matbio.2005.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Solé-Boldo, L., Raddatz, G., Schütz, S., Mallm, J-P., Rippe, K., Lonsdorf, A. S., et al. (2020). Single-cell transcriptomes of the human skin reveal age-related loss of fibroblast priming. Commun. Biol. 3 (1), 188. doi:10.1038/s42003-020-0922-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, M., Luo, E. Y., Adams, S. M., Adams, T., Ye, Y., Shetye, S. S., et al. (2020). Collagen XI regulates the acquisition of collagen fibril structure, organization and functional properties in tendon. Matrix Biol. 94, 77–94. doi:10.1016/j.matbio.2020.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Tawa, G. J., Braisted, J., Gerhold, D., Grewal, G., Mazcko, C., Breen, M., et al. (2021). Transcriptomic profiling in canines and humans reveals cancer specific gene modules and biological mechanisms common to both species. PLoS Comput. Biol. 17 (9), e1009450. doi:10.1371/journal.pcbi.1009450

PubMed Abstract | CrossRef Full Text | Google Scholar

Trojanowska, M. (2010). Cellular and molecular aspects of vascular dysfunction in systemic sclerosis. Nat. Rev. Rheumatol. 6 (8), 453–460. doi:10.1038/nrrheum.2010.102

PubMed Abstract | CrossRef Full Text | Google Scholar

Vogiatzi, M. G., Li, D., Tian, L., Garifallou, J. P., Kim, C. E., Hakonarson, H., et al. (2018). A novel dominant COL11A1 mutation in a child with Stickler syndrome type II is associated with recurrent fractures. Osteoporos. Int. 29 (1), 247–251. doi:10.1007/s00198-017-4229-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Walz, J. M., Wecker, T., Zhang, P. P., Cakir, B., Gruening, B., Agostini, H., et al. (2019). Impact of angiogenic activation and inhibition on miRNA profiles of human retinal endothelial cells. Exp. Eye Res. 181, 98–104. doi:10.1016/j.exer.2019.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Warde-Farley, D., Donaldson, S. L., Comes, O., Zuberi, K., Badrawi, R., Chao, P., et al. (2010). The GeneMANIA prediction server: Biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 38, W214–W220. Web Server issue. doi:10.1093/nar/gkq537

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, P., Yu, H. V., Tseng, K-C., Flath, M., Fabian, P., Segil, N., et al. (2021). Foxc1 establishes enhancer accessibility for craniofacial cartilage differentiation. Elife 10, e63595. doi:10.7554/eLife.63595

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, F., Luo, P., Ding, H., Zhang, C., and Zhu, Z. (2018). Collagen type V a2 (COL5A2) is decreased in steroid-induced necrosis of the femoral head. Am. J. Transl. Res. 10 (8), 2469–2479.

PubMed Abstract | Google Scholar

Yao, Q., Xing, Y., Wang, Z., Liang, J., Lin, Q., Huang, M., et al. (2020). MiR-16-5p suppresses myofibroblast activation in systemic sclerosis by inhibiting NOTCH signaling. Aging (Albany NY) 13 (2), 2640–2654. doi:10.18632/aging.202308

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L., Tang, Y., Zhu, X., Tu, T., Sui, L., Han, Q., et al. (2017). Overexpression of MiR-335-5p promotes bone formation and regeneration in mice. J. Bone Min. Res. 32 (12), 2466–2475. doi:10.1002/jbmr.3230

CrossRef Full Text | Google Scholar

Zhang, Y., Liao, Y., Chen, C., Sun, W., Sun, X., Liu, Y., et al. (2020). p38-regulated FOXC1 stability is required for colorectal cancer metastasis. J. Pathol. 250 (2), 217–230. doi:10.1002/path.5362

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, Q., Wang, D., Lin, R., Lv, Q., and Wang, W. (2022). IFI44 is an immune evasion biomarker for SARS-CoV-2 and Staphylococcus aureus infection in patients with RA. Front. Immunol. 13, 1013322. doi:10.3389/fimmu.2022.1013322

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, B., and Yu, J-W. (2017). A novel identified circular RNA, circRNA_010567, promotes myocardial fibrosis via suppressing miR-141 by targeting TGF-β1. Biochem. Biophys. Res. Commun. 487 (4), 769–775. doi:10.1016/j.bbrc.2017.04.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y., Zhou, B., Pache, L., Chang, M., Khodabakhshi, A. H., Tanaseichuk, O., et al. (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 10 (1), 1523. doi:10.1038/s41467-019-09234-6

PubMed Abstract | CrossRef Full Text | Google Scholar