![Open Access](/resources/images/iconopenaccess.png)
Integrative meta-analysis of differentially expressed genes in osteoarthritis using microarray technology
- Authors:
- Published online on: May 15, 2015 https://doi.org/10.3892/mmr.2015.3790
- Pages: 3439-3445
-
Copyright: © Wang et al. This is an open access article distributed under the terms of Creative Commons Attribution License [CC BY_NC 3.0].
Abstract
Introduction
Osteoarthritis (OA) is the most prevalent joint disease and is characterized by an abnormal remodeling of joint tissues, which is predominantly driven by inflammatory mediators within the affected joint (1,2). The pathological changes of OA primarily take place in the articular cartilage, and include cartilage degeneration, matrix degradation and synovial inflammation (3–5). Clinically, features of OA include pain, stiffness, limitation of motion, swelling and deformity (6). Synovial inflammation is hypothesized to be the primary underlying etiology in OA (3). However, the biological mechanisms associated with OA remain to be elucidated.
Microarray, a high-throughput genomics technology, has been developed in order to improve the understanding of complex interactions and networks in disease development (7). Thousands of genes on a genome-wide scale have been measured using microarray technology (8). The successful identification of gene expression signatures that may provide insights into OA pathogenesis and differentiate the diseased state from a healthy state, requires an adequate sample size and heterogeneous datasets (9). Although numerous microarray studies have generated lists of differentially expressed (DE) genes, there are inconsistencies among the results of such studies, due to the limitation of the small sample sizes involved (10).
To overcome these difficulties, meta-analysis has previously been applied to publicly-available genome-wide expression datasets from studies on a number of diseases (11–13). The use of meta-analysis may improve reliability and generalizability, and permit a more precise estimation of gene expression (11). Meta-analyses provide enhanced statistical power, thereby obtaining more robust and reliable gene signatures (7,14–17). Recently, integrative meta-analysis of expression data (INMEX), a new user-friendly microarray meta-analysis tool, has been developed to support meta-analysis of multiple gene expression datasets (18).
In order to overcome the limitations of individual studies, resolve inconsistencies in results, and reduce false-positive or false-negative associations due to random errors, a microarray meta-analysis was performed in the present study. The objective was to identify differentially expressed (DE) genes and biological processes associated with gene expression signature in OA.
Materials and methods
Identification of eligible gene expression datasets of OA
A search of microarray datasets was performed that examined DE genes between OA and healthy controls. The NCBI Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/) (19) was used to identify suitable microarray datasets. The keyword 'osteoarthritis' was used for this search. Studies were included in the analysis if they: Were based on gene expression profiling of blood or synovial membrane samples; contained sufficient data to perform a meta-analysis; and included patients diagnosed with OA, based on OA classification criteria (20). The following information was extracted from each of the studies that were selected: GEO accession; sample type; platform; numbers of patients and healthy controls; and gene expression data.
Meta-analysis of microarray datasets
All available OA microarray datasets that met the inclusion criteria were downloaded from the NCBI GEO database. Data tables containing gene expression or relative gene expression values were constructed, with genes/probes in the rows, and samples in the columns. The datasets were uploaded to INMEX (http://www.inmex.ca/INMEX) (18), and the data was subsequently annotated by converting different gene or probe ID to Entrez IDs. For each probe-set, intensity values were log-transformed and/or normalized to zero mean and unit variance, which is the normalization method for high density oligonucleotide array data, as reported by Bolstad et al (21). When all datasets had been uploaded, processed and annotated, a data integrity check was performed prior to the meta-analysis stage.
The random effects model presumes that different studies present substantial diversity, and evaluates between study variance as well as within study sampling error (22,23). The random effects model is used when the between-study heterogeneity is significant (23). The INMEX program was used to conduct statistical analysis (18).
Functional analysis
The functional analysis of INMEX generates new hypotheses by exploiting characteristics of the DE gene lists identified in meta-analysis. A heat map created by 'Pattern extractor' produced gene expression profiles across different datasets/conditions.
In order to examine the functions of the genes in the gene list, gene ontology (GO) enrichment analysis was performed, which detected enriched functional attributes based on gene-associated GO terms, using the hypergeometric test (http://www.geneontology.org/) (24). Functional analysis was performed using the INMEX program (18).
Results
Studies included in the meta-analysis
Three GEO data sets, which met the inclusion criteria, were identified (Table I) (4,25,26). These datasets consisted of two synovial membrane datasets and one blood dataset, and included a total of 137 patients with OA and 52 controls. Selected details of the individual studies are summarized in Table I.
Meta-analysis of gene expression in OA
A random effects model of effect size (ES) measures was used to integrate gene expression patterns. The present study incorporated between-study heterogeneities across studies, because the estimated Q-value was not in an approximately chi-squared distribution. DE genes with P<0.05 were selected. In the current analysis, 1 'gained' gene and 13402 'lost' genes were identified (Fig. 1). Gained genes are DE genes that were only identified in the meta-analysis (26). The single gained gene exhibited relatively weak but consistent expression profiles across the three different datasets. The large sample size obtained by consisting of the datasets made it possible declare this a DE gene with increased certainty. Lost genes are genes which were identified as DE genes in any of the individual analyses, but not in the meta-analysis. These genes also presented conflicting changes in expression profiles, or large variations across the different studies.
Identification of differentially expressed genes in OA
A total of 85 genes were identified, which were consistently differentially expressed in OA. Among these 85 DE genes, 30 were upregulated and 55 were downregulated. A list of the top 20 upregulated and downregulated genes is shown in Table II The upregulated gene with the lowest P-value (P=5.36E-07) was S-phase kinase-associated protein 2, E3 ubiquitin protein ligase (SKP2). The downregulated gene with the lowest P-value (P=4.42E-09) was Proline rich 5 like (PRR5L).
Identification of differentially expressed genes in the synovial membrane of patients with OA
An additional meta-analysis was performed on the results from the 2 synovial membrane samples, following exclusion of the third study, which used peripheral blood. A list of the top 20 upregulated and downregulated genes is shown in Table III. The upregulated genes with the lowest P-values (both P=0.004003) were JAZF zinc finger 1 (JAZF1) and Guanine nucleotide binding protein (G-protein), β polypeptide 4 (GNB4), which are involved in coupling membrane receptors to effector proteins, such as ion channels and enzymes (27). The downregulated gene with the largest ES (ES=2.0472; P=0.004003) was multiple inositol-polyphosphate phosphatase 1 (MINPP1). A number of the downregulated genes were related to inflammatory factors (Table III).
Functional analysis
GO analysis of the DE genes was performed in order to identify the biological processes associated with changes in gene expression in OA. The analysis identified 210 significant enrichments of the DE genes, which were categorized to 10 GO terms (Fig. 2). The two enrichments with the lowest P-values were in the GO category of 'Immune response', with a P-value of 0.000129438, and 'Immune effector process', with a P-value of 0.000288619. Other significant GO categories included 'Regulation of humoral immune response' (P=0.000308832), 'Regulation of immune response' (P=0.00055514) and 'Positive regulation of immune system process' (P=0.00059351; Table IV).
Discussion
A number of genes are differentially expressed genes between patients with OA and healthy controls, and it is necessary to identify the genes that may enhance understanding of the molecular and cellular processes, which are involved in the pathogenesis of OA. Although a large quantity of data may be produced using microarray studies, the small sample size of these studies is a significant obstacle to the identification of DE genes. A meta-analysis of multiple microarray datasets increases the sample size, rendering the identification of DE genes more reliable.
In the present study, a meta-analysis was performed using three publicly available GEO datasets in order to identify common biological mechanisms involved in the pathogenesis of OA. The analysis identified 85 genes that were consistently differentially expressed in OA (30 upregulated and 55 downregulated). The upregulated gene with the largest ES was SKP2, which is known to be involved in the inhibition of cell growth and the promotion of apoptosis. Kitagawa concluded that SKP2 controls the p300–p53 signaling pathway in cancer cells (28). Furthermore, this gene encodes a member of the F-box protein family, which is characterized by a ~40 amino acid motif, the F-box (29). The downregulated gene with the lowest P-value was PRR5L, which suppresses mTOR complex 2 (mTORC2)-mediated hydrophobic motif phosphorylation of protein kinase C, but not that of protein kinase B (30). In addition, the PRR5L protein may function to modulate the activity of mTORC2 in a substrate-dependent manner (30). Actinin α 1 (ACTN1), an upregulated gene, encodes an actin-binding protein, which exerts multiple effects in a variety of cell types. ACTN1 may protect osteoclasts from tumor necrosis factor-α (TNF-α); induce apoptosis through increasing the expression of the anti-apoptotic protein, Bcl-2; activate survival signals; and promote Akt phosphorylation and NF-κB activation (31). Although it is currently unclear exactly how these genes contribute to OA, they may be useful as potential biomarkers to facilitate early diagnosis or to monitor the efficacy of treatment in this disease. A number of these genes provide insights into the molecular mechanisms underlying the pathophysiology of OA.
Although osteoarthritis (OA) is understood to be a degradative articular cartilage disease, there is increasing data demonstrating the involvement of the immune system. In recent epidemiological studies involving a large number of patients with OA, an inflammatory synovium has been shown to be involved in increased damage to the cartilage (32) and pain (33). Immune cells, such as T cells, B cells and macrophages, have been identified in the synovial tissue of patients with OA (34–36). Furthermore, immunoglobulins and immune complexes against cartilage components have been detected in the plasma, synovium and cartilage of patients with OA (37), and it has been shown that the synovium is involved in complement activation in OA (38). In the present study, 210 significantly enriched GO terms associated with the DE genes were identified using a meta-analysis. The three enriched terms with the lowest P-values were 'Immune response', 'Immune effector process' and 'Regulation of humoral immune response', which were all involved in the immune system. The identified GO terms may be grouped into a smaller number of categories: 'Response to stimulus', 'Signal transduction', 'Regulation of response to stimulus', 'Immune system process', 'Immune response', 'Regulation of apoptotic process', 'Regulation of programmed cell death', 'Regulation of cell death', 'Apoptotic process' and others. Although it is difficult to identify all the significant functional categories that are expressed differentially in OA, the GO categories identified here, merit further investigation in subsequent studies.
There were certain limitations to the present study, which ought to be considered. Firstly, heterogeneity and confounding factors may have distorted the analysis. Clinical samples may have been heterogeneous with respect to clinical activity, severity or gender. Secondly, there are differences in gene expression between tissues, such as blood and synovial membrane, that were not considered. Although an additional subgroup analysis of the synovial membrane samples was performed, this only included two studies. By contrast, the initial meta-analysis integrated the results obtained from different tissues, which should have enabled detection of the genes that may have been missed in an analysis of two studies only.
In conclusion, the meta-analysis of microarray studies that was performed in the present study, provided an overview of differential gene expression in OA; identifying 85 differential expressed genes (30 upregulated and 55 downregulated genes). Future studies to validate these genes as markers for the diagnosis and response to biological therapy for OA may provide further insight into their involvement in the development and progression of OA.
Acknowledgments
This study was supported by the Key Scientific and Technological Innovation Special Projects of Shaanxi '13115' of China (grant no. 2009ZDKG-79) and the National Natural Science Foundation of China (grant no. 81472924).
References
Aigner T, Rose J, Martin J and Buckwalter J: Aging theories of primary osteoarthritis: From epidemiology to molecular biology. Rejuvenation Res. 7:134–145. 2004. View Article : Google Scholar : PubMed/NCBI | |
Loeser RF, Goldring SR, Scanzello CR and Goldring MB: Osteoarthritis: A disease of the joint as an organ. Arthritis Rheum. 64:1697–1707. 2012. View Article : Google Scholar : PubMed/NCBI | |
Benito MJ, Veale DJ, FitzGerald O, van den Berg WB and Bresnihan B: Synovial tissue inflammation in early and late osteoarthritis. Ann Rheum Dis. 64:1263–1267. 2005. View Article : Google Scholar : PubMed/NCBI | |
Ma CH, Lv Q, Cao Y, et al: Genes relevant with osteoarthritis by comparison gene expression profiles of synovial membrane of osteoarthritis patients at different stages. Eur Rev Med Pharmacol Sci. 18:431–439. 2014.PubMed/NCBI | |
Duan C, Guo X, Zhang XD, et al: Comparative analysis of gene expression profiles between primary knee osteoarthritis and an osteoarthritis endemic to Northwestern China, Kashin-Beck disease. Arthritis Rheum. 62:771–780. 2010. View Article : Google Scholar : PubMed/NCBI | |
Sellam J and Berenbaum F: The role of synovitis in pathophysiology and clinical symptoms of osteoarthritis. Nat Rev Rheumatol. 6:625–635. 2010. View Article : Google Scholar : PubMed/NCBI | |
Choi SJ, Rho YH, Ji JD, Song GG and Lee YH: Genome scan meta-analysis of rheumatoid arthritis. Rheumatology (Oxford). 45:166–170. 2006. View Article : Google Scholar | |
Golub TR, Slonim DK, Tamayo P, et al: Molecular classification of cancer: Class discovery and class prediction by gene expression monitoring. Science. 286:531–537. 1999. View Article : Google Scholar : PubMed/NCBI | |
Ramasamy A, Mondry A, Holmes CC and Altman DG: Key issues in conducting a meta-analysis of gene expression microarray datasets. PLoS Med. 5:e1842008. View Article : Google Scholar : PubMed/NCBI | |
Siddiqui AS, Delaney AD, Schnerch A, Griffith OL, Jones SJ and Marra MA: Sequence biases in large scale gene expression profiling data. Nucleic Acids Res. 34:e832006. View Article : Google Scholar : PubMed/NCBI | |
Griffith OL, Melck A, Jones SJ and Wiseman SM: Meta-analysis and meta-review of thyroid cancer gene expression profiling studies identifies important diagnostic biomarkers. J Clin Oncol. 24:5043–5051. 2006. View Article : Google Scholar : PubMed/NCBI | |
Rung J and Brazma A: Reuse of public genome-wide gene expression data. Nat Rev Genet. 14:89–99. 2013. View Article : Google Scholar | |
Cahan P, Rovegno F, Mooney D, Newman JC, St Laurent G III and McCaffrey TA: Meta-analysis of microarray results: Challenges, opportunities, and recommendations for standardization. Gene. 401:12–18. 2007. View Article : Google Scholar : PubMed/NCBI | |
Lee YH and Nath SK: Systemic lupus erythematosus susceptibility loci defined by genome scan meta-analysis. Hum Genet. 118:434–443. 2005. View Article : Google Scholar : PubMed/NCBI | |
Lee YH, Witte T, Momot T, et al: The mannose-binding lectin gene polymorphisms and systemic lupus erythematosus: Two case-control studies and a meta-analysis. Arthritis Rheum. 52:3966–3974. 2005. View Article : Google Scholar : PubMed/NCBI | |
Lee YH, Rho YH, Choi SJ, Ji JD and Song GG: Ankylosing spondylitis susceptibility loci defined by genome-search meta-analysis. J Hum Genet. 50:453–459. 2005. View Article : Google Scholar : PubMed/NCBI | |
Lee YH, Rho YH, Choi SJ, Ji JD and Song GG: Meta-analysis of genome-wide linkage studies for bone mineral density. J Hum Genet. 51:480–486. 2006. View Article : Google Scholar : PubMed/NCBI | |
Xia J, Fjell CD, Mayer ML, Pena OM, Wishart DS and Hancock RE: INMEX - a web-based tool for integrative meta-analysis of expression data. Nucleic Acids Res. 41:W63–W70. 2013. View Article : Google Scholar : PubMed/NCBI | |
Barrett T, Troup DB, Wilhite SE, et al: NCBI GEO: Archive for functional genomics data sets - 10 years on. Nucleic Acids Res. 39:D1005–D1010. 2011. View Article : Google Scholar | |
Altman R, Asch E, Bloch D, et al: Development of criteria for the classification and reporting of osteoarthritis. Classification of osteoarthritis of the knee. Diagnostic and Therapeutic Criteria Committee of the American Rheumatism Association. Arthritis Rheum. 29:1039–1049. 1986. View Article : Google Scholar : PubMed/NCBI | |
Bolstad BM, Irizarry RA, Astrand M and Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 19:185–193. 2003. View Article : Google Scholar : PubMed/NCBI | |
Choi JK, Yu U, Kim S and Yoo OJ: Combining multiple microarray studies and modeling interstudy variation. Bioinformatics. 19(Suppl 1): i84–i90. 2003. View Article : Google Scholar : PubMed/NCBI | |
DerSimonian R and Laird N: Meta-analysis in clinical trials. Control Clin Trials. 7:177–188. 1986. View Article : Google Scholar : PubMed/NCBI | |
Falcon S and Gentleman R: Using GOstats to test gene lists for GO term association. Bioinformatics. 23:257–258. 2007. View Article : Google Scholar | |
Ramos YF, Bos SD, Lakenberg N, et al: Genes expressed in blood link osteoarthritis with apoptotic pathways. Ann Rheum Dis. 73:1844–1853. 2013. View Article : Google Scholar : PubMed/NCBI | |
Lambert C, Dubuc JE, Montell E, et al: Gene expression pattern of cells from inflamed and normal areas of osteoarthritis synovial membrane. Arthritis Rheumatol. 66:960–968. 2014. View Article : Google Scholar : PubMed/NCBI | |
Ruiz-Velasco V, Ikeda SR and Puhl HL: Cloning, tissue distribution, and functional expression of the human G protein beta 4-subunit. Physiol Genomics. 8:41–50. 2002. View Article : Google Scholar : PubMed/NCBI | |
Kitagawa M, Lee SH and McCormick F: Skp2 suppresses p53-dependent apoptosis by inhibiting p300. Mol Cell. 29:217–231. 2008. View Article : Google Scholar : PubMed/NCBI | |
Bai C, Sen P, Hofmann K, et al: SKP1 connects cell cycle regulators to the ubiquitin proteolysis machinery through a novel motif, the F-box. Cell. 86:263–274. 1996. View Article : Google Scholar : PubMed/NCBI | |
Gan X, Wang J, Wang C, et al: PRR5L degradation promotes mTORC2-mediated PKC-δ phosphorylation and cell migration downstream of Gα12. Nat Cell Biol. 14:686–696. 2012. View Article : Google Scholar : PubMed/NCBI | |
Triplett JW and Pavalko FM: Disruption of alpha-actinin-integrin interactions at focal adhesions renders osteoblasts susceptible to apoptosis. Am J Physiol Cell Physiol. 291:C909–C921. 2006. View Article : Google Scholar : PubMed/NCBI | |
Ayral X, Pickering EH, Woodworth TG, Mackillop N and Dougados M: Synovitis: A potential predictive factor of structural progression of medial tibiofemoral knee osteoarthritis - results of a 1 year longitudinal arthroscopic study in 422 patients. Osteoarthritis Cartilage. 13:361–367. 2005. View Article : Google Scholar : PubMed/NCBI | |
Hill CL, Hunter DJ, Niu J, et al: Synovitis detected on magnetic resonance imaging and its relation to pain and cartilage loss in knee osteoarthritis. Ann Rheum Dis. 66:1599–1603. 2007. View Article : Google Scholar : PubMed/NCBI | |
Revell PA, Mayston V, Lalor P and Mapp P: The synovial membrane in osteoarthritis: A histological study including the characterisation of the cellular infiltrate present in inflammatory osteoarthritis using monoclonal antibodies. Ann Rheum Dis. 47:300–307. 1988. View Article : Google Scholar : PubMed/NCBI | |
Sakkas LI, Scanzello C, Johanson N, et al: T cells and T-cell cytokine transcripts in the synovial membrane in patients with osteoarthritis. Clin Diagn Lab Immunol. 5:430–437. 1998.PubMed/NCBI | |
Nakamura H, Yoshino S, Kato T, Tsuruha J and Nishioka K: T-cell mediated inflammatory pathway in osteoarthritis. Osteoarthritis Cartilage. 7:401–402. 1999. View Article : Google Scholar : PubMed/NCBI | |
Jasin HE: Immune mediated cartilage destruction. Scand J Rheumatol Suppl. 76:111–116. 1988. View Article : Google Scholar : PubMed/NCBI | |
Wang Q, Rozelle AL, Lepus CM, et al: Identification of a central role for complement in osteoarthritis. Nat Med. 17:1674–1679. 2011. View Article : Google Scholar : PubMed/NCBI |