Microarray Analysis of Hub Genes and Pathways in Damaged Cartilage Tissues of Knee

doi:10.21203/rs.3.rs-112849/v1

Download PDF

Research article

Microarray Analysis of Hub Genes and Pathways in Damaged Cartilage Tissues of Knee

https://doi.org/10.21203/rs.3.rs-112849/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 17 Sep, 2021

Read the published version in Medicine →

Version 1

posted

You are reading this latest preprint version

Background: The aim of this study was to identify genes and functional pathways associated with damaged cartilage tissues of knee using microarray analysis.

Methods: The gene expression profile GSE129147 including including 10 knee cartilage tissues from damaged side and 10 knee non-weight-bearing healthy cartilage was downloaded and bioinformatics analysis was made.

Results: A total of 182 DEGs including 123 up-regulated and 59 down-regulated genes were identified from the GSE129147 dataset. GO and pathway enrichment analysis confirmed that extracellular matrix organization, collagen catabolic process, antigen processing and presentation of peptide or polysaccharide antigen, and endocytic vesicle membrane were strongly associated with cartilage injury. Furthermore, 10 hub DEGs with a higher connectivity degree in PPI network were found such as POSTN, FBN1, LOX, IGFBP3, C3AR1, MMP2, ITGAM, CDKN2A, COL1A1, COL5A1.

Conclusion: These hub genes and pathways provide a new perspective for revealing the potential pathological mechanisms and therapy strategy of cartilage injury.

Orthopedic Surgery

Cartilage injury

Bioinformatics analysis

Differentially expressed genes

Protein-protein interaction network

Cartilage is an avascular, non-innervated and alymphatic tissue with limited intrinsic repair potential. Articular cartilage lesions of the knee are commonly encountered in all ages. A retrospective study of 31,516 knee arthroscopies found 63% of their patients had cartilage injuries [1]. Another retrospective study of 25,124 knee arthroscopies observed 60% of patients who suffered from knee pain had articular cartilage lesions [2]. Due to its limited capacity to spontaneously repair largely, focal chondral defects can extend and progress into degeneration of the whole joint [3, 4]. So, focal chondral lesions of the knee are regarded as a potential risk factor for later osteoarthritis [5].

Knee osteoarthritis is the result of imbalance between anabolic and catabolic metabolism of articular cartilage matrix; the main risk factor is aging, but the degeneration can accelerated by mechanical factors, such as chondral lesions [6]. Surgical treatments for chondral lesions including chondrocyte implantation, osteochondral transplantation, and microfracture, provide not fully satisfying outcomes [6–9]. In this regard, understanding the molecular mechanisms and biological pathways that regulate chondrocyte metabolism and steady of cartilage matrix and gene therapy that enables the transfer of suitable genes into cartilage tissue is necessary for the treatment of cartilage lesions. Mori et al. identified Fibroblast Growth factor-18 maintain the homeostasis of adult rats articular cartilage by decreasing glycosaminoglycan release and depletion from the cartilage, and enhancing proliferation of articular chondrocytes [10]. Delivery of genes encoding for TGF-β, IGF-1, and BMP-2 by adenovirus vector to chondrocytes has been demonstrated to enhance proteoglycan synthesis by the cells in culture [11, 12].

Whole-genome microarray analysis is an efficient and high-throughput molecular technique for analyzing general genetic alterations in various tissues [13–15]. For example, Aşık et al. compared damaged and undamaged non-weight bearing cartilage and showed differentially-expressed genes (DEGs) and multiple biological processes related to these DEGs [16]. And, genome microarray combined with bioinformatic analysis can help us to better understand the underlying mechanisms that involved in development of disease.

The aim of this study was to identify DEGs and hub genes in damaged and healthy cartilage samples using bioinformatics methods and then to investigate the potential biological processes and pathways of these changes.

Affymetrix Microarray Data

The microarray data of GSE129147 was downloaded from National Center of Biotechnology Information (NCBI) Gene Expression Omnibus database (GEO, http://www.ncbi.nlm.nih.gov/geo/), based on GPL15207 platform ([PrimeView] Affymetrix Human Gene Expression Array; Affymetrix, Inc.), including 10 knee cartilage tissues from damaged side and 10 knee non-weight-bearing healthy cartilage [16].

Identification of DEGs in damaged and healthy cartilage samples

The DEGs between damaged and healthy cartilage samples were screened by GEO2R (http://www.ncbi.nlm.nih.gov/geo/geo2r/). Genes with the adjusted p-value < 0.05 and |log FC| > 1 or |log FC| < -1 were regarded as the cutoff value for DEGs screening, which were the signature genes of damaged cartilage of knee.

GO enrichment function and pathway analysis

To determine the DEGs involved in biological processes (BPs), cellular components (CCs), and molecular functions (MFs), gene ontology (GO) enrichment analysis were used with The Database for Annotation, Visualization and Integrated Discovery (DAVID) 6.8 [17]. Pathway analysis of DEGs was also performed with DAVID online to learn the significantly altered metabolic pathways [18]. GO enrichment analysis (P < 0.05, FDR < 0.05) and pathway analysis (P < 0.05) of the DEGs was considered significant.

Construction of the Protein–protein interactions (PPI) network

Search Tool for the Retrieval of Interacting Genes (STRING, version 11.0) was used to analyze the interaction between proteins for the DEGs with default parameters [19]. Medium confidence > 0.4 was defined as the cutoff criterion. PPI network was visualized using Cytoscape (version 3.7.2, http://www.cytoscape.org/) after downloading the raw image from STRING [20]. In the network, the node represents a protein, the line represents the interaction, and the degree represents the number of interactions. And connectivity degree was analyzed by statistics in the network to obtain the important nodes, namely hub proteins [21].

Module analysis in the PPI network

Molecular Complex Detection (MCODE) was used to discover dense clique-like structures within a network in the Cytoscape platform. Modules were extracted according to the criteria: Degree Cutoff = 2, Max. Depth = 100, K-Core = 2, Node Score Cutoff = 0.2, and MCODE score ≥ 4.

DEGs identification and GO and pathway enrichment analysis

In comparing damaged to healthy cartilage samples, a total of 182 DEGs including 123 up-regulated and 59 down-regulated genes were identified from the GSE129147 dataset. After GO functional enrichment analysis, 3 significant BP terms of up-regulated DEGs including: extracellular matrix organization, collagen catabolic process, and antigen processing and presentation of peptide or polysaccharide antigen via MHC class II; 3 significant CC terms of up-regulated DEGs including: extracellular matrix, endocytic vesicle membrane, MHC class II protein complex; 1 significant MF terms of up-regulated DEGs is peptide antigen binding (Table 1). The up-regulated DEGs were enriched in 3 pathways such as Staphylococcus aureus infection, Leishmaniasis, and Phagosome. In terms of down-regulated DEGs, we did not find significant GO terms or KEGG pathways according to stated criteria.

Table 1

The enriched GO terms and KEGG pathways for the down-regulated DEGs.
Category	Term	Description	Count	Genes	P value	FDR
BP	GO:0030198	extracellular matrix organization	10	B4GALT1, COL7A1, FBN1, POSTN, COL1A1, LOX, COL8A1, LOXL1, ITGAM, COL5A1	1.38E-06	2.00E-03
BP	GO:0030574	collagen catabolic process	7	COL7A1, MMP19, COL1A1, MMP14, COL8A1, MMP2, COL5A1	1.42E-06	2.00E-03
BP	GO:0002504	antigen processing and presentation of peptide or polysaccharide antigen via MHC class II	5	HLA-DRB1, HLA-DRB3, HLA-DRB4, HLA-DRB5, HLA-DRA	1.99E-06	3.00E-03
CC	GO:0031012	extracellular matrix	13	ACTG1, HMCN1, COL7A1, MMP19, FBN1, POSTN, COL1A1, LOXL2, COL8A1, MMP14, MMP2, LOXL1, COL5A1	5.63E-08	6.95E-05
CC	GO:0030666	endocytic vesicle membrane	7	WNT5B, HLA-DRB1, HLA-DRB3, HLA-DRB4, HLA-DRB5, HLA-DRA, CD163	1.37E-06	2.00E-03
CC	GO:0042613	MHC class II protein complex	5	HLA-DRB1, HLA-DRB3, HLA-DRB4, HLA-DRB5, HLA-DRA	5.13E-06	6.00E-03
MF	GO:0042605	peptide antigen binding	6	HLA-DRB1, HLA-DRB3, HLA-DRB4, HLA-DRB5, SLC7A5, HLA-DRA	4.27E-07	5.33E-04
KEGG	hsa05150	Staphylococcus aureus infection	8	C3AR1, HLA-DRB1, HLA-DRB3, HLA-DRB4, HLA-DRB5, FCGR3A, ITGAM, HLA-DRA	1.87E-07	2.20E-04
KEGG	hsa05140	Leishmaniasis	8	IRAK1, HLA-DRB1, HLA-DRB3, HLA-DRB4, HLA-DRB5, FCGR3A, ITGAM, HLA-DRA	1.26E-06	1.49E-03
KEGG	hsa04145	Phagosome	9	ACTG1, ATP6V1A, HLA-DRB1, HLA-DRB3, HLA-DRB4, HLA-DRB5, FCGR3A, ITGAM, HLA-DRA	2.17E-05	2.55E-02

PPI network construction

The PPI network was screened, comprising 132 nodes and 215 edges (Fig. 1). The top 10 DEGs with a higher connectivity degree are all up-regulated genes (Table 2), including POSTN, FBN1, LOX, IGFBP3, C3AR1, MMP2, ITGAM, CDKN2A, COL1A1, COL5A1, which are hub proteins. Among these genes, MMP2 (degree = 20) and ITGAM (degree = 20) showed the highest node degree.

Table 2

Results of connectivity degree of top 10 DEGs in the PPI network.
Genes	Degree
MMP2	20
ITGAM	20
COL1A1	17
LOX	14
C3AR1	14
POSTN	12
FBN1	11
IGFBP3	11
CDKN2A	11
COL5A1	11

Module analysis

A total of 3 significant modules were selected with the parameter of MCODE score ≥ 4. Module A (MCODE score = 8.889), Module B (MCODE score = 5.000), Module C (MCODE score = 4.000) were shown in Fig. 2. Module A had 10 nodes involving 10 up-regulated gene (LOXL1, LOX, MMP14, LOXL2, MMP2, FBN1, COL5A1, POSTN, COL1A1, ADAM12); Module B with 5 nodes and 10 edges involving 5 up-regulated gene (CDKN3, MAD2L1, PRC1, TYMS, TK1); Module C with 4 nodes and 6 edges involving 4 up-regulated gene (FCGR3A, CD163, VSIG4, MS4A6A).

In the present study, we explored the hub genes and pathways in damaged cartilage tissues of knee by bioinformatics methods. By comparing gene expression profiles between 10 knee cartilage from damaged side and 10 knee non-weight-bearing healthy cartilage tissues, we found 182 DEGs, including 123 up-regulated and 59 down-regulated genes in damaged cartilage tissues. Subsequently, GO enrichment analysis found that DEGs were associated with various biological processes and molecular functions, such as extracellular matrix organization, collagen catabolic process, antigen processing and presentation of peptide or polysaccharide antigen via MHC class II, and endocytic vesicle membrane. In addition, a PPI network with these DEGs was constructed, and 10 hub genes, including POSTN, FBN1, LOX, IGFBP3, C3AR1, MMP2, ITGAM, CDKN2A, COL1A1, COL5A1, were identified as the hub genes with higher connectivity degree in damaged cartilage tissues.

In view of the results of GO terms enrichment analysis, we linked the up-regulated gene DEGs with extracellular matrix organization and collagen catabolic process. The structure and organization of cartilage’s extracellular matrix are the primary determinants of normal function. Most diseases involving cartilage lead to dramatic changes in the extracellular matrix. Well organized elastic networks have been found in the superficial zone of articular cartilage using immunohistochemistry or multiphoton microscopy together with histochemistry. Fibrillin-1(FBN1) as the major component of elastic networks, mainly present in the uppermost superficial zone of articular cartilage [22]. Collagens such as COL1A1 and COL5A1 and FBN1 were detectable in damaged cartilage tissues compared to undamaged controls in this study. The significant difference in expression levels of COL1A1, COL5A1, and FBN1 between damaged and undamaged sites of cartilage suggesting the attempt of chondrocytes for recovery of the damaged cartilage. Meanwhile, several matrix degrading proteases previously described to be linked to cartilage degeneration were up-regulated in damaged cartilage, such as MMP2 and POSTN. Study considered POSTN as a catabolic protein that promotes cartilage degeneration through collagen and proteoglycan degradation in OA by up-regulating MMP-13 and ADAMTS4 [23]. These differentially-expressed genes demonstrated the balance between synthetic and catabolic activities of chondrocytes is struggling after cartilage injuries.

Lysyl oxidase (LOX), is a copper-dependent amine oxidase, that catalyzes cross-linking of collagen and elastin, which plays an important role in biological functions of extracellular matrix [24]. LOXL2 was visualized in tissues from human knee and hip joints by immunofluorescence. LOXL2 is upregulated in cartilage affected by OA, may be a protective response that promotes anabolism while inhibiting specific catabolic responses in the pathophysiology of OA [25]. LOXL3, a downstream of leptin, stimulates the apoptosis, but inhibits the autophagy of chondrocytes, is a potential therapy target for osteoarthritis [26]. In this study, LOX is identified as an up-regulated gene in damaged cartilage tissues, which could promote cartilage maintenance.

IGFBP3 is one of six insulin-like growth factor binding proteins (IGFBPs) in humans. IGFBP3 modulates IGF signaling by binding to the two ligands (IGF-I and IGF-II) that primarily signal through IGF-I receptor [27]. In healthy articular cartilage from individuals without OA, IGF-I activates anabolism and inhibits catabolism of cartilage. Higher levels of IGFBP3 have been observed in OA cartilage, leading to the decrease of IGF-I’s anabolic activity and could be partially responsible for the reduced responsiveness of OA cartilage to IGF-I [28, 29]. And a meta-analysis of genome-wide association studies of hip osteoarthritis found suggestive links between IGFBP3 and hip osteoarthritis and IGFBP3 overexpression induced cartilage catabolism and osteogenic differentiation [30]. In this study, IGFBP3 is identified as an up-regulated gene in damaged cartilage tissues, which could promote cartilage degeneration.

Cells release into the extracellular environment different types of membrane vesicles called exosomes or microvesicles. Kim at al. found dendritic cells or vesicles derived from the dendritic cells expressing either secreted IL-4 or membrane-bound IL-4 are able to modulate the activity of APC and T cells in vivo through a MHC class II and partly Fas ligand/Fas-dependent mechanism, resulting in effective treatment of established collagen-induced arthritis [31]. And Kim at al. also found exosomes derived from IL-10-treated dendritic cells can suppress inflammation and collagen-induced arthritis [32]. In this study, we linked the up-regulated gene DEGs in antigen processing and presentation of peptide or polysaccharide antigen via MHC class II and endocytic vesicle membrane through GO terms enrichment analysis, which suggests extracellular vesicle and inflammation implicated in the pathogenesis of cartilage injuries.

In conclusion, our data provide a comprehensive bioinformatics analysis of DGEs and pathways, which may be involved in synthetic and catabolic metabolisms of chondrocytes after cartilage injuries. Nevertheless, additional experiments are needed to further validate the identifed genes and pathways.

DEGs: diﬀerentially-expressed genes; BPs: biological processes; CCs: cellular components; MFs: molecular functions; GO: gene ontology; PPI: protein–protein interactions; FBN1: Fibrillin-1; LOX: Lysyl oxidase; IGFBPs: insulin-like growth factor binding proteins.

Ethics approval and consent to participate

Not applicable. This paper does not involve research on humans.

Consent for publication

All authors read the final manuscript and approved for publication.

Availability of data and materials

The data was freely downloaded from the public GEO database.

Competing interests

The authors declare that they have no competing interests.

Funding

This research was supported by the Natural Science Foundation of Inner Mongolia (No. 2016MS0366; No. 2020MS08119).

Authors contributions

YYN and CXH conceived the design of the study. YZR performed and collected the data. YTS analyzed the data. Aoqierbatu and AL prepared and revised the manuscript. All authors read and approved the final content of the manuscript.

Acknowledgements

None.

Curl WW, Krome J, Gordon ES, Rushing J, Smith BP, Poehling GG: Cartilage injuries: a review of 31,516 knee arthroscopies. Arthroscopy 1997, 13(4):456-460.
Widuchowski W, Widuchowski J, Trzaska T: Articular cartilage defects: study of 25,124 knee arthroscopies. Knee 2007, 14(3):177-182.
Buckwalter JA, Mankin HJ, Grodzinsky AJ: Articular cartilage and osteoarthritis. Instr Course Lect 2005, 54:465-480.
Ossendorff R, Franke K, Erdle B, Uhl M, Südkamp NP, Salzmann GM: Clinical and radiographical ten years long-term outcome of microfracture vs. autologous chondrocyte implantation: a matched-pair analysis. Int Orthop 2019, 43(3):553-559.
Buckwalter JA, Mankin HJ: Articular cartilage: degeneration and osteoarthritis, repair, regeneration, and transplantation. Instr Course Lect 1998, 47:487-504.
Dall'Oca C, Cengarle M, Costanzo A, Giannini N, Vacchiano A, Magnan B: Current concepts in treatment of early knee osteoarthritis and osteochondral lesions; the role of biological augmentations. Acta Biomed 2017, 88(4S):5-10.
Na Y, Shi Y, Liu W, Jia Y, Kong L, Zhang T, Han C, Ren Y: Is implantation of autologous chondrocytes superior to microfracture for articular-cartilage defects of the knee? A systematic review of 5-year follow-up data. Int J Surg 2019, 68:56-62.
Tírico LEP, McCauley JC, Pulido PA, Bugbee WD: Lesion Size Does Not Predict Outcomes in Fresh Osteochondral Allograft Transplantation. Am J Sports Med 2018, 46(4):900-907.
Frehner F, Benthien JP: Microfracture: State of the Art in Cartilage Surgery? Cartilage 2018, 9(4):339-345.
Mori Y, Saito T, Chang SH, Kobayashi H, Ladel CH, Guehring H, Chung U-i, Kawaguchi H: Identification of fibroblast growth factor-18 as a molecule to protect adult articular cartilage by gene expression profiling. J Biol Chem 2014, 289(14):10192-10200.
Smith P, Shuler FD, Georgescu HI, Ghivizzani SC, Johnstone B, Niyibizi C, Robbins PD, Evans CH: Genetic enhancement of matrix synthesis by articular chondrocytes: comparison of different growth factor genes in the presence and absence of interleukin-1. Arthritis Rheum 2000, 43(5):1156-1164.
Arai Y, Kubo T, Kobayashi K, Takeshita K, Takahashi K, Ikeda T, Imanishi J, Takigawa M, Hirasawa Y: Adenovirus vector-mediated gene transduction to chondrocytes: in vitro evaluation of therapeutic efficacy of transforming growth factor-beta 1 and heat shock protein 70 gene transduction. J Rheumatol 1997, 24(9):1787-1795.
Cui S, Zhang X, Hai S, Lu H, Chen Y, Li C, Tong P, Lu F, Yuan Z: Molecular mechanisms of osteoarthritis using gene microarrays. Acta Histochem 2015, 117(1):62-68.
Li PCH: Overview of Microarray Technology. Methods Mol Biol 2016, 1368:3-4.
Zeng M, Liu J, Yang W, Zhang S, Liu F, Dong Z, Peng Y, Sun L, Xiao L: Multiple-microarray analysis for identification of hub genes involved in tubulointerstial injury in diabetic nephropathy. J Cell Physiol 2019:10.1002/jcp.28313.
Aşık MD, Gürsoy S, Akkaya M, Kozacı LD, Doğan M, Bozkurt M: Microarray analysis of cartilage: comparison between damaged and non-weight-bearing healthy cartilage. Connect Tissue Res 2019:1-9.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT et al: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 2000, 25(1):25-29.
Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res 2000, 28(1):27-30.
Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, Santos A, Doncheva NT, Roth A, Bork P et al: The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res 2017, 45(D1):D362-D368.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003, 13(11):2498-2504.
He X, Zhang J: Why do hubs tend to be essential in protein networks? PLoS Genet 2006, 2(6):e88-e88.
Yu J, Urban JP: The elastic network of articular cartilage: an immunohistochemical study of elastin fibres and microfibrils. Journal of anatomy 2010, 216(4):533-541.
Attur M, Yang Q, Shimada K, Tachida Y, Nagase H, Mignatti P, Statman L, Palmer G, Kirsch T, Beier F et al: Elevated expression of periostin in human osteoarthritic cartilage and its potential role in matrix degradation via matrix metalloproteinase-13. FASEB J 2015, 29(10):4107-4121.
Smith-Mungo LI, Kagan HM: Lysyl oxidase: properties, regulation and multiple functions in biology. Matrix biology : journal of the International Society for Matrix Biology 1998, 16(7):387-398.
Alshenibr W, Tashkandi MM, Alsaqer SF, Alkheriji Y, Wise A, Fulzele S, Mehra P, Goldring MB, Gerstenfeld LC, Bais MV: Anabolic role of lysyl oxidase like-2 in cartilage of knee and temporomandibular joints with osteoarthritis. Arthritis research & therapy 2017, 19(1):179.
Huang ZM, Du SH, Huang LG, Li JH, Xiao L, Tong P: Leptin promotes apoptosis and inhibits autophagy of chondrocytes through upregulating lysyl oxidase-like 3 during osteoarthritis pathogenesis. Osteoarthritis and cartilage 2016, 24(7):1246-1253.
Jogie-Brahim S, Feldman D, Oh Y: Unraveling insulin-like growth factor binding protein-3 actions in human disease. Endocrine reviews 2009, 30(5):417-437.
Morales TI: The quantitative and functional relation between insulin-like growth factor-I (IGF) and IGF-binding proteins during human osteoarthritis. Journal of orthopaedic research : official publication of the Orthopaedic Research Society 2008, 26(4):465-474.
Doré S, Pelletier JP, DiBattista JA, Tardif G, Brazeau P, Martel-Pelletier J: Human osteoarthritic chondrocytes possess an increased number of insulin-like growth factor 1 binding sites but are unresponsive to its stimulation. Possible role of IGF-1-binding proteins. Arthritis Rheum 1994, 37(2):253-263.
Evans DS, Cailotto F, Parimi N, Valdes AM, Castaño-Betancourt MC, Liu Y, Kaplan RC, Bidlingmaier M, Vasan RS, Teumer A et al: Genome-wide association and functional studies identify a role for IGFBP3 in hip osteoarthritis. Annals of the rheumatic diseases 2015, 74(10):1861-1867.
Kim SH, Bianco NR, Shufesky WJ, Morelli AE, Robbins PD: Effective treatment of inflammatory disease models with exosomes derived from dendritic cells genetically modified to express IL-4. Journal of immunology (Baltimore, Md : 1950) 2007, 179(4):2242-2249.
Kim SH, Lechman ER, Bianco N, Menon R, Keravala A, Nash J, Mi Z, Watkins SC, Gambotto A, Robbins PD: Exosomes derived from IL-10-treated dendritic cells can suppress inflammation and collagen-induced arthritis. Journal of immunology (Baltimore, Md : 1950) 2005, 174(10):6440-6448.

Download PDF

Journal Publication

published 17 Sep, 2021

Read the published version in Medicine →

Version 1

posted

You are reading this latest preprint version

Microarray Analysis of Hub Genes and Pathways in Damaged Cartilage Tissues of Knee

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction

Materials And Methods

Results

DEGs identification and GO and pathway enrichment analysis

PPI network construction

Module analysis

Discussion

Abbreviations

Declarations

References

Status:

Journal Publication

Version 1