Authors: Venugopal Panga aff001; Ashwin Adrian Kallor aff001; Arunima Nair aff001; Shilpa Harshan aff001; Srivatsan Raghunathan aff001 Authors place of work: Institute of Bioinformatics and Applied Biotechnology (IBAB), Bengaluru, Karnataka, India aff001; Manipal Academy of Higher Education, Manipal, Karnataka, India aff002 Published in the journal: PLoS ONE 14(11) Category: Research Article doi: https://doi.org/10.1371/journal.pone.0224632
Several studies have reported mitochondrial dysfunction in rheumatoid arthritis (RA). Many nuclear DNA (nDNA) encoded proteins translocate to mitochondria, but their participation in the dysfunction of this cell organelle during RA is quite unclear. In this study, we have carried out an integrative analysis of gene expression, protein-protein interactions (PPI) and gene ontology data. The analysis has identified potential implications of the nDNA encoded proteins in RA mitochondrial dysfunction. Firstly, by analysing six synovial microarray datasets of RA patients and healthy controls obtained from the gene expression omnibus (GEO) database, we found differentially expressed nDNA genes that encode mitochondrial proteins. We uncovered some of the roles of these genes in RA mitochondrial dysfunction using literature search and gene ontology analysis. Secondly, by employing gene co-expression from microarrays and collating reliable PPI from seven databases, we created the first mitochondrial PPI network that is specific to the RA synovial joint tissue. Further, we identified hubs of this network, and moreover, by integrating gene expression and network analysis, we found differentially expressed neighbours of the hub proteins. The results demonstrate that nDNA encoded proteins are (i) crucial for the elevation of mitochondrial reactive oxygen species (ROS) and (ii) involved in membrane potential, transport processes, metabolism and intrinsic apoptosis during RA. Additionally, we proposed a model relating to mitochondrial dysfunction and inflammation in the disease. Our analysis presents a novel perspective on the roles of nDNA encoded proteins in mitochondrial dysfunction, especially in apoptosis, oxidative stress-related processes and their relation to inflammation in RA. These findings provide a plethora of information for further research.
Gene expression – Genetic networks – inflammation – mitochondria – apoptosis – Microarrays – Gene ontologies – Protein interaction networks
Mitochondrial dysfunction prevails among numerous diseases, including RA, Sjøgren’s syndrome, neurodegenerative diseases, diabetes, cancer and obesity [1–5]. Genomic technologies and computational approaches played a vital role in our understanding of mitochondrial dysfunction in several diseases like Leigh syndrome, cardiovascular diseases, obesity and infantile-onset mitochondrial encephalopathy [6–12]. These approaches have also discerned the mechanics of calcium uniporter in mitochondrial biology and associated diseases [7, 13–15]. Further, investigations into metabolic profiling and whole-exome sequencing data point to metabolic abnormalities concerned with mitochondria and biallelic mutations leading to instability in mitoribosomal subunits in Leigh syndrome [6, 8].
Mitochondria, which are membrane-bound cell organelles, are the primary generators of adenosine triphosphate (ATP). The respiratory chain complexes, which are part of the mitochondrial oxidative phosphorylation (OxPhos), are necessary for the production of ATP. The genome of this organelle has 13 protein-coding genes, which are associated with the OxPhos pathway. It is understood that 1158 nDNA encoded proteins get translocated to this cell organelle [16], and some of them are crucial for the OxPhos pathway. However, the functional roles of many of these proteins in RA mitochondrial dysfunction are uncertain, creating a serious lacuna in our understanding of this disease. An integrative analysis of these proteins using gene expression, PPI, gene ontology and network theory offers an excellent opportunity for deducing some of their roles.
About 1% of the world’s population is affected by RA [17]. It is a chronic inflammatory disease that usually affects the small synovial joints of the hands and feet. The disease synovium gets inflamed (a condition called synovitis) and invades articular cartilage and bone, forming a layer of granulation tissue called pannus. Further, synovitis causes irreversible damage to the synovium in joints [18]. Moreover, the cells of the RA synovium (synoviocytes) secrete inflammatory cytokines and articular cartilage-degrading enzymes, such as matrix metalloproteinases (MMPs), which further aggravate the disease.
The composition of cell types in a healthy synovium is different to that of RA. The healthy synovium primarily contains two cell types, macrophage-like synoviocytes (MLS) and fibroblast-like synoviocytes (FLS) [19]. Other cell types such as leucocytes can be seen in small numbers [19]. In contrast, the RA synovium is expanded and forms pannus and contains resident MLS and FLS as well as heavily infiltrated leucocytes [20–21].
The pannus in RA, like a tumour, increases demand for energy (ATP) in the synovium. Additionally, the dysregulated synovial microvasculature results in a poor supply of oxygen to the tissue, causing hypoxia. Both the increased energy demand on mitochondrial electron transport and hypoxia could lead to an enhanced production of ROS, creating oxidative stress in synoviocytes [1]. Further, an inverse correlation between synovitis and the partial pressure of oxygen in the synovium testifies to the role of hypoxia in arthritis [22]. Moreover, hypoxia might induce proinflammatory pathways, through hypoxia-inducible factor-1α (HIF-1α), nuclear factor κB (NF-κB), Janus kinase-signal transducer and activator of transcription (JAK-STAT), activator protein 1 (AP-1) and Notch. Most notably, anti-tumour necrosis factor therapy has significantly decreased synovial hypoxia in vivo, indicating that it is a crucial event in arthritis [1, 22–25]. This elucidates that hypoxia and ROS are relevant to RA mitochondrial dysfunction.
Superoxide anion (O2·-), hydrogen peroxide (H2O2) and hydroxyl radical (·OH) are collectively called ROS [26–27]. The components of ROS can damage DNA, proteins, lipids and many other molecules. Synovial fluid (SF) and plasma samples as well as blood lymphocytes and polymorphonuclear leucocytes from RA patients have significantly higher mitochondrial DNA (mtDNA) and oxidatively damaged DNA adduct, 8-hydroxyl-2′-deoxyguanosine (8-oxodG), than non-arthritic samples. Further, both the mtDNA and 8-oxodG levels in SF correlate with the presence of rheumatoid factor in RA patients [28–29]. This underlines the existence of ROS-mediated damage of mitochondria in this disease. Other oxidative stress markers, such as protein carbonyls are significantly higher in the serum of RA patients compared to healthy controls. Treatment of these patients with infliximab resulted in a significant decrease of the carbonyls [30]. Iron, a catalyst for the formation of ·OH from H2O2 via the Fenton’s reaction, is present in the diseased synovium [31]. Oxidised low-density lipoproteins and lipid peroxidation as well as the latter’s correlation with the concentration of Iron ions were observed in SF of RA patients [32–33]. Furthermore, hyaluronate-derived small oligosaccharides are present in the inflamed disease joints, revealing the activity of ROS [34]. The ROS-mediated damage of mitochondria might also result in angiogenesis and cartilage destruction, the latter of which ensues through the up-regulation of MMPs [1, 25, 35–36]. Besides, there is an inverse association between the dietary intake of antioxidants and the prevalence of RA as well as the levels of antioxidants and the disease inflammation [37–42]. Moreover, an element of ROS, O2·- reacts with nitric oxide (NO) to form peroxynitrite (ONOO-), which is a component of the reactive nitrogen species (RNS). This reactive species plays a role in the NF-κB-mediated production of inflammatory mediators, such as tumour necrosis factor (TNF), interleukin-1 beta (IL-1β) and inducible nitric oxide synthase (iNOS) [43].
To summarise, the pannus increases ATP demand and the dysregulated microvasculature creates hypoxia. Both the conditions can generate ROS in RA synovial mitochondria and the immediate targets of these free radicals are mtDNA, proteins and lipids. Supporting this phenomenon in RA, elevated levels of the damaged mtDNA, proteins and lipids were observed in SF, plasma and leucocytes of the patients. Additionally, both the hypoxia and ROS are known to induce pro-inflammatory HIF-1α, NF-κB, JAK-STAT, AP-1 and Notch pathways. As stated earlier, 1158 nDNA encoded proteins get translocated to mitochondria and several of them could be involved in the pathways concerned with the generation of ROS. So, it is of great pathophysiological relevance to elucidate the roles of these proteins in mitochondrial dysfunction and their connection to ROS-mediated damage, hypoxia, ATP synthesis and inflammation in RA.
In RA, apoptosis is required to control synovial hyperplasia. Apoptosis can occur by two different pathways, the extrinsic and the intrinsic, of which the latter could be initiated in mitochondria in response to oxidative stress. Both the pathways culminate in the activation of a cascade of proteases, called caspases. It has been shown that the extrinsic pathway is inactive in RA FLS [44]. Fas, which is a pro-apoptotic molecule and known to be involved in the extrinsic pathway, has been found to induce inflammation rather than apoptosis in RA FLS. However, this process depends on caspase-8 (CASP8) activity and FLICE-like inhibitory protein (FLIP) expression [45]. Therefore, studying the intrinsic pathway might give clues on the regulation of synovial hyperplasia. Hence it is important to understand the roles of the nDNA encoded proteins that could be implicated in this pathway.
In the current study, we followed an integrated approach that uses microarray data, PPI, gene ontology and network analysis. Six microarray gene expression datasets related to RA and healthy synovium were obtained from GEO, and they were analysed to discover differentially expressed genes (DEGs) encoding mitochondrial proteins. Further, a mitochondrion-specific PPI network has been created based on the information from seven publicly available databases. We also performed gene ontology analysis (GO) using the Search Tool for the Retrieval of Interacting Genes (STRING) database for identifying significantly enriched biological processes (BP), molecular functions (MF) and cellular components (CC). In addition to a discussion on the roles of nDNA encoded proteins in RA mitochondrial dysfunction based on available information in the literature, a model for the relation between mitochondrial dysfunction and the disease inflammation has also been framed.
The mRNA expression datasets (GSE77298, GSE7307, GSE12021, GSE55235 and GSE55457) were retrieved from GEO, which is a public National Center for Biotechnology Information (NCBI) database. Table 1 gives a detailed account of the mRNA expression datasets. All the datasets were downloaded in raw data file format for analysis.
The mitochondrial PPI network in RA synovium was created by pooling the experimentally determined interactions in human cells. They were obtained from seven publicly available resources, namely the biological general repository for interaction datasets (BioGRID), IntAct, the molecular interaction (MINT), STRING, the human protein reference database (HPRD), the database of interacting proteins (DIP) and CRG [46–52]. Among them, the first four databases have confidence scores for each interaction. The higher the score the more is the confidence for the interaction to occur. For the current study, from each of these four, we have got more reliable interactions by putting a cut-off to the confidence scores. The cut-off was decided in such a way that the interactions having a confidence score more than the median of the score distributions were selected. From DIP, the interactions with the core quality status were considered. From HPRD and CRG, which do not have confidence scores, only those interactions which have at least two publication evidences were considered. Collectively, a total of 387,242 interactions were obtained from all the seven resources. Then, to create the mitochondrial PPI network, only the interactions of those proteins that get localised to mitochondria were chosen using MitoCarta [16], which is a compendium of 1158 nDNA genes that encode mitochondrial proteins.
Furthermore, to make this interactome specific to the synovial tissue, we measured the co-expression of the interacting partners of these interactions using the gene expression data from six microarray datasets. Table 1 gives detailed information about these datasets. For each dataset, the raw intensities were normalised using the RMA algorithm. For the interacting partners of each interaction, we computed the Pearson correlation coefficient of the normalised expression values across all disease samples. Only those interactions with a Pearson correlation coefficient > 0.7 between the partners, in at least one microarray dataset, were considered co-expressed in the synovial tissues. The resulting interactions were used to create the undirected mitochondrial PPI network, using the ‘igraph’ package in R. The hubs of this network were identified using the same package in R.
The microarray experiments, considered in this study, were carried out on RA and normal synovial tissues by other workers (Table 1). The RA samples used in these studies were obtained by tissue excision upon joint replacement/synovectomy surgery from RA patients. Similarly, the control samples were obtained from either postmortem or traumatic joint injury cases. In four of the six datasets (GSE12021 (HGU133A), GSE12021 (HGU133B), GSE55235 and GSE55457), for which the information on duration and severity of the disease is available, the duration of the disease in the patients was reported to be a mean of at least 12 years. The number of American rheumatism association (ARA) (now, American college of Rheumatology) criteria for RA was reported to be a mean of at least five [53–54]. The patients, who participated in five of the six studies, are from the Netherlands and Germany. For one study (GSE7307), the demography of patients is not available.
We re-analysed all the datasets using the R/Bioconductor statistical package. The intensities were normalised using two algorithms, MAS5 and RMA, separately. The differential expression of the genes between RA and control groups was computed using the two sample independent t-test. A p-value < 0.05 and a fold-change of > 1.5 in the up or down direction were taken as the cut-off values for differential expression. Further, the following conditions were imposed for deciding a differentially expressed gene across the datasets:
For one dataset, if the gene is selected by both the normalisation methods in the same direction (up or down)
For multiple datasets, if the gene is selected by both the normalisation methods in at least one dataset or by complementary normalisation methods in at least two datasets.
If the gene is up-regulated in at least one dataset and not down-regulated in any of the remaining, we call it a consistently up-regulated gene. A similar criterion was applied for a down-regulated gene. On the other hand, if a gene shows up-regulation in some and down-regulation in the other datasets, we call it a mixed-regulated gene.
Since we are particularly interested in nuclear genes that encode mitochondrial proteins, and in order to maximise the DEGs of mitochondrial proteins, we did not correct the p-value. However, for most of the analyses performed, the genes that were selected in at least two or three datasets were considered. Further, the DEGs were used for integrative analysis and hence the false positives might be reduced.
The GO and pathway enrichment analyses were carried out using the STRING database. These analyses identify enriched GO terms and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways for a given list of genes by employing a hypergeometric test that was discussed elsewhere [55–56]. A false discovery rate (FDR) < 0.01 was considered as the cut-off for the significantly enriched GO terms and pathways.
High-confident PPI from seven public resources, namely BioGRID, IntAct, MINT, STRING, HPRD, DIP and CRG, were used to construct this network [46–52]. The first four databases provide a confidence score, which is a measure of reliability, for each interaction. From each of these databases, only the interactions with the scores above the median of the confidence score distributions were extracted (Fig 1). From DIP, only those interactions which have core quality status (a reliability parameter specific to DIP) were extracted. From HPRD and CRG, the interactions with at least two publication evidences were selected. This resulted in 387,242 reliable interactions which constitute ~40% of the total interactions in these databases. Then, we applied a filter requiring both the interacting partners to be nDNA encoded mitochondrial proteins, the list of which could be found in MitoCarta [16]. This has returned an interactome with 7023 interactions, representing 926 of the 1158 MitoCarta genes (79.96%). In order to make the interactome specific to the RA synovium, we computed the co-expressions between the interacting partners using RA synovial microarray data (Table 1). For each pair of interacting proteins, their expression levels across the RA disease samples in a given microarray dataset were used to compute the Pearson’s correlation coefficient (ρ) as a measure of co-expression. The interacting partners with a ρ > 0.7 between their gene expression values in at least one of the six microarray datasets were considered to be co-expressed in the synovium. This selection criterion, which was chosen to maximise the number of PPI with the co-expressed interacting partners, resulted in an interactome with 2708 interactions and 665 genes, representing 57.42% of MitoCarta genes. In this interactome, on an average, each protein is connected to four other proteins. None of the interacting partners were co-expressed in all the six datasets. Of the 2708 interactions, 13 have the co-expressed interacting partners in five microarray datasets; 65 in four; 184 in three; 618 in two and 1828 in one datasets. With these interactions, we have created the mitochondrial PPI network using the ‘igraph’ package in R. In order to identify the localisation of the network proteins in mitochondria, a GO analysis for the cellular component (CC) term was performed using the STRING database [55–56]. It was observed that the majority of the network proteins get translocated to the mitochondrial inner membrane and matrix (S1 Fig).
Differential expression analysis between RA and healthy human synovial tissue samples was carried out using the six microarray datasets obtained from GEO (Table 1). The requirement for differential expression in a dataset was set to be a fold-change > 1.5 (up or down regulation) and a p-value < 0.05. Of the 665 PPI network genes, 131 were found to be differentially expressed in at least one RA synovial microarray dataset (S1 Table). The criterion of a gene having a differential expression in at least one of the six datasets was decided so as to maximise the selection of mitochondrial DEGs. The whole network indicating the up and down DEGs, which can be visualised using the Cytoscape tool, is in S1 File. The 131 genes include 46 consistently up-, 73 consistently down - and 12 mixed-regulated genes (for methodological details, see ‘Methods‘ section). Another 77 of the MitoCarta members, which are not part of the network, were also differentially expressed in at least one dataset (S2 Table). Thus, it makes a total of 208 mitochondrial DEGs (83 up-, 111 down - and 14 mixed-regulated). Their differential expression across the six studies is as follows; only one of the 208 genes was a DEG in all the six studies; three (1.44%) genes in five studies; four (1.92%) in four; 15 (7.21%) in three; 60 (28.84%) in two and 125 (60.09%) in one study.
Mitochondrial DEGs that were found in at least three datasets are in Table 2 (13 up-, 6 down - and 4 mixed-regulated). Among them, the following up-regulated genes are highlighted in respect of their functions in mitochondria. Acyl-CoA thioesterase 7 (ACOT7) is an enzyme that hydrolyses long-chain fatty acids such as palmitoyl-CoA. Kynurenine 3-monooxygenase (KMO) is an enzyme that catalyses the hydroxylation of kynurenine to form 3-hydroxykynurenine. This enzyme has been reported to be involved in the generation of oxidative radicals as well as in cytokine-mediated inflammation [57]. Leucine amino peptidase 3 (LAP3) is involved in the degradation of glutathione, a scavenger of free radicals [58], indicating the likely impairment in the detoxification of ROS. Pyruvate dehydrogenase kinase 1(PDK1) inhibits pyruvate dehydrogenase activity and is known to play an integral role against hypoxia - and oxidative stress-mediated apoptosis [59]. Interferon alpha inducible protein 27 (IFI27) is known to be involved in cytokine signalling and apoptosis. Interestingly, this protein activates an apoptotic caspase, CASP8, which was also up-regulated in the microarray data [60]. Uncoupling protein 2 (UCP2) is implicated in the transfer of anions from the inner to the outer membrane and protons from the outer to the inner membrane, and it is known to control ROS [61]. Peroxiredoxin-4 (PRDX4) is an antioxidant enzyme which detoxifies H2O2 and regulates NF-κB activation [62]. Nonetheless, because of its high reactivity, this enzyme is susceptible to overoxidation and inactivation by H2O2 [63]. YME1 like 1 ATPase (YME1L1), which is an ATP-dependent metalloprotease, is known to function in the maintenance of mitochondrial morphology and accumulation of respiratory chain subunits [64–65]. Isocitrate dehydrogenase 2 (IDH2) is implicated in the production of NADPH and in the protection of cells from ROS. DnaJ heat shock protein family (Hsp40) member C15 (DNAJC15) negatively regulates respiratory chain and generation of ATP.
Similarly, the six genes that were down-regulated in at least three datasets participate in the following functions. MCCC1 encodes α subunit of 3-methylcrotonoyl-CoA carboxylase (3-MCC), which is an enzyme that is involved in the breakdown of leucine. Monoamine oxidase A (MAOA) catalyses the oxidative deamination of amines, such as serotonin, norepinephrine and dopamine, and its deficiency is known to induce aggression [66–67]. The transcription factors, specificity protein 1 (SP1), GATA binding protein 2 (GATA2) and TATA box binding protein (TBP) regulate the expression of this gene [68]. EF-hand domain-containing protein 1 (EFHD1) is a calcium ion sensor. Some of the other down-regulated genes are AKR1B10, PDK4 and C10orf2.
The mixed-regulated gene, solute carrier family 16 member 7 (SLC16A7) is involved in the transport of metabolites such as monocarboxylates and pyruvate. Similarly, adenylate kinase 4 (AK4) is implicated in the metabolism of nucleotides.
The largest genome-wide association study meta-analysis of RA cases and controls has identified 98 disease risk genes [69]. Six of them, C1QBP, SUOX, ACSL6, UNG, CYP27B1 and CASP8 are MitoCarta genes. Among these, only CASP8 was a DEG, and C1QBP and CYP27B1 are part of the created mitochondrial network. The role of these genes in RA and their involvement in mitochondrial dysfunction remain to be ascertained.
Of the six open-source microarray datasets we analysed, RA patients in one (GSE7307) were not treated with therapies, while the patients belonging to three others (GSE12021 (HGU133A), GSE12021 (HGU133B) and GSE55457) underwent different combinations of medical therapies. The information on medications is not available for two datasets (GSE55235 and GSE77298). All the details of medical therapies available for the datasets are listed in Table 3.
It is seen within a dataset that some patients have received the same combination of medical therapies whereas others received different combinations. To test if the gene expressions are influenced by these medical therapies, the samples in each microarray dataset were hierarchically clustered based on the mRNA levels of the DEGs that were differentially expressed in at least three microarray datasets (Table 2). The cluster results are shown as heatmaps with dendrograms (S2–S7 Figs). In GSE7307, GSE55235 and GSE55457, RA and control samples were clustered into separate groups (S2–S4 Figs). In GSE12021 (HGU133A) and GSE12021 (HGU133B), some RA samples were clustered into a separate group while others were clustered with control samples (S5 and S6 Figs), showing that there is a drug effect. In GSE77298, some RA samples were clustered with healthy controls but drug therapies are not available for this dataset (S7 Fig).
In order to find the effect of medical therapies on the differential expression of genes, we removed the RA samples that were clustered with healthy controls from GSE12021 (HGU133A) and GSE12021 (HGU133B) datasets and repeated the differential expression analysis for the 23 genes listed in Table 2. Surprisingly, with the same selection criteria of differential expression, all the 23 genes were retained. The heatmaps of the expression levels of the genes in these two datasets after eliminating the RA samples that clustered with healthy controls are shown in S8 and S9 Figs. We notice the complete separation of controls from RA samples in the clusters.
From the above analysis, we find that the 23 genes were differentially expressed in at least three datasets in both of the following cases.
Case 1: all RA and control samples in all the six studies
Case 2: all RA samples except those that clustered with healthy controls in GSE12021 (HGU133A) and GSE12021 (HGU133B) and all controls in all the six studiesSince the 23 genes were differentially expressed in at least three datasets in both the cases, we conclude that these genes were not affected by the therapy initiation.
In addition to the above analysis, we analysed two other microarray datasets (GSE77344 and GSE11237) where patients with diseases unrelated to RA were treated with prednisone or celecoxib. Prednisone is the prodrug form of prednisolone, while celecoxib is a COX-2 inhibitor. Prednisolone and COX-2 inhibitors are part of the therapies received by RA patients shown in Table 3. In the dataset GSE77344 [70], whole blood was collected from patients with chronic obstructive pulmonary disease who were either treated (n = 31) or not treated (n = 103) with prednisone. GSE11237 [71] contained colorectal primary adenocarcinomas surgically removed from 23 patients, 11 of whom received 400 mg celecoxib two times per day for seven days prior to surgery and 12 who did not receive the treatment. In addition to this, we also analysed GSE45867 [72] which had paired synovial tissue biopsies from 8 early RA patients naive to methotrexate or DMARDs. The samples were collected before and 12 weeks after the initiation of methotrexate therapy. Hierarchical clustering of the samples based on the expression levels of the 23 genes normalized across samples did not show any separation of treated and nontreated samples, or pre and post treatment samples (S10–S12 Figs). Differential expression analysis with the criteria used for the RA datasets (fold change > |1.5| and p value < = 0.05) revealed one gene out of the 23 was upregulated in GSE77344 (MAOA, fold change = 3.9, pvalue = 0.001), while no differential regulation was found for any of the 23 genes in GSE11237 and GSE45867. Thus we believe that the effects of these specific treatments on the candidate genes are negligible, and the differential regulation observed in the RA datasets is more likely due to the disease itself.
To further elucidate the properties of the mitochondrial PPI network, we performed network analysis. For each node, we chose to measure the network parameter ‘degree’ which is the number of edges a node can have. The probability distribution of the degree of nodes in the created mitochondrial PPI network along with power-law fit to the data is shown in Fig 2. The degree distribution of the network follows a power law P(k) ~ k-α (with the degree coefficient, α = 1.82), which is a property of scale-free networks [73]. From this network, we identify a small number of important nodes, called hubs, which are directly connected to a large number of interacting partners. Analogous to social networks, the hub proteins with a higher number of neighbours are crucial to PPI networks as their removal causes dysfunction of the system.
The immediate neighbours of all the 665 proteins in the mitochondrial PPI network were determined. The top 50 proteins in the decreasing order of the number of their immediate neighbours are listed in Table 4. The entire list of all the network proteins and the number of their immediate neighbours—including the extent of DEGs among them—could be found in S3 Table. The distributions of the proteins in terms of the total number of neighbours and the proportion of DEGs among them are shown in S13 Fig.
Each of the network proteins has at least one neighbour. Among them, 167 have at least 10 neighbours. Most of the network proteins are connected to one or a few DEGs. The scatter plot between the number of neighbours and the number of DEGs among the neighbours for individual proteins is shown in Fig 3. It would be interesting to look at the hubs with a high number of neighbours containing higher number of DEGs among them. For example, the upper right-side rectangle of the figure has the hubs connected to at least 27 neighbours having a minimum of seven DEGs among them. The hubs which have fallen into this rectangle are given in Table 5. They could be considered crucial for mitochondrial functions in the RA diseased synovium because of a high number of DEGs among the neighbours.
We chose four representative hubs from the top right-side rectangle to draw their subnetworks with immediate connecting proteins (points marked red in the scatter plot, Fig 3). One of these hubs, UQCR10, which is a subunit of the respiratory chain complex III, is connected to 50 proteins, including two up - and six down-regulated DEGs (Table 5). The DEGs include the up-regulated complex 1 subunit NDUFB7 and complex III subunit UQCR11; the down-regulated complex IV subunit COX7A1, complex I subunits, NDUFA4, NDUFB4, NDUFB6 and NDUFB9, and complex III subunit UQCRFS1 (S14 Fig). Another hub, NDUFV2 is linked to 49 proteins, including eight DEGs, seven of which, namely NDUFA4, NDUFB4, NDUFB6, NDUFB9, NDUFS4, PITRM1 and UQCRFS1, were down-regulated and one gene was mixed-regulated (S15 Fig).
The third hub, NDUFS6 links to 11 DEGs: four of these were up-regulated and are part of complex III (UQCR11), complex I (NDUFB7), complex IV (COX15) and complex V (ATP5E); six were down-regulated, of which ECHDC2 is implicated in the lyase activity, and the rest (NDUFA4, NDUFB4, NDUFB6, NDUFB9 and NDUFS4) are complex 1 subunits (S16 Fig).
The hub protein UQCRQ links to nine neighbours, of which three were up-regulated (UQCR11, NDUFB7 and ATP5E); and five were down-regulated (NDUFA4, NDUFB4, NDUFB6, NDUFS4 and UQCRFS1) (S17 Fig).
The GO and pathway analyses were performed for the mitochondrial DEGs that were selected in at least two datasets (83 genes), using the STRING database. The results illustrated their roles at different levels, including MF, BP and CC. A total of 14 MF, 63 BP and 16 CC GO terms, and 3 KEGG pathways were significantly enriched (FDR < 0.01). Several of these terms and pathways were shared by both the up - and down-regulated DEGs. Some of the significantly enriched BP and MF GO terms are shown in Fig 4. The detailed lists of all the GO terms and KEGG pathways could be found in S4 Table and Tables 6–8.
Among the BP terms, the processes pertaining to metabolism, mitochondrial membrane permeability, regulation of membrane potential, oxidative stress, mitochondrial transport and apoptotic processes were enriched with DEGs (Fig 4 and S4 Table). Similarly, among the MF terms, the functions related to the binding of coenzymes and cofactors were enormously enriched with DEGs (Fig 4 and Table 6). Additionally, DEGs were enriched in many mitochondrial CC terms (Table 7).
Among the enriched KEGG pathways, ‘metabolic pathways’ (KEGG pathway ID: 1100) is highly enriched with 20 DEGs. Of the others, glycine, serine, threonine and glutathione metabolisms were affected (Table 8).
From the enriched MF items, it is understood that the oxidoreductase activity is getting affected in RA. This activity is associated with the OxPhos complexes of mitochondria. There are five complexes: complex I, II, III, IV and V. Electrons from NADH and FADH2 pass through the first four complexes and eventually reduce O2 to water at complex IV. Overall, the complexes have 97 subunits, 84 of which are encoded by nDNA. The created mitochondrial PPI network contains 81 of the 84 subunits. Totally, 11 of them were DEGs in one or two datasets (Table 9). The complex1 DEGs, NDUFB4, NDUFB6, NDUFB9 and NDUFS4 were down-regulated. At least one gene each in complex III, complex IV and complex V was down or up-regulated. UQCRFS1 of complex III, COX6A1 and COX7A1 of complex IV, and ATP5G3 of complex V were down-regulated. The mitochondrial protein-coding genes of OxPhos were either missed or non-DEGs in the microarray studies. From these observations, it can be deduced that OxPhos is getting disrupted, perhaps impaired due to the decreased expression of some of the OxPhos genes in the RA synovium. This might be concerned with the escape of electrons from OxPhos, leading to the formation of ROS. Nevertheless, since 11 of the subunits are DEGs in either one or two datasets, it may be difficult to come to a concrete conclusion.
The ROS generating NADPH oxidases (NOXs), such as NOX1, NOX2 and NOX3 were not DEGs, and NOX4 was down in two and up in one microarray datasets. This might be indicating that it is the mitochondria, and not these enzymes, that could be the primary source of ROS in RA. The detoxifiers of ROS, such as catalase (CAT), glutathione peroxidase 4 (GPX4) and superoxide dismutase 1 (SOD1) were down-regulated in one dataset, specifying impaired detoxification. Further, LAP3, which degrades glutathione, was overexpressed in five datasets. This gives a clue for the accumulation of mitochondrial ROS (mtROS), which might induce apoptosis in RA synoviocytes. In support of this, the pro-apoptotic BAX was up-regulated in two datasets. The executor of apoptosis, CASP8 was up in three. However, CASP8, in the presence of FLIP, is known to induce inflammation rather than apoptosis [45]. Further, H2O2 might trigger apoptosis by causing elevation of intracellular Ca2+ levels via a pathway that includes spleen tyrosine kinase (SYK), Bruton’s tyrosine kinase (BTK), the B cell linker protein (BLNK) and phospholipase Cγ2 (PLCγ2) [74–75]. Interestingly, these four genes were up-regulated in at least three datasets. SYK, BTK, BLNK and PLCγ2 were up-regulated in six, three, five and six datasets, respectively. But, Bcl-2-like 1 (BCL2L1) protein, which controls the production of ROS by regulating membrane potential, was mixed-regulated, making it difficult to conclude on its role. Further, the anti-apoptotic protein Bcl-2 (BCL2) was mixed-regulated. Surprisingly, the pro-apoptotic Bcl-2-like 13 (BCL2L13) was down in one dataset. Collectively, the results may suggest (i) generation of mtROS, (ii) impaired ROS detoxification and (iii) induction of mitochondrion-initiated intrinsic apoptosis in the RA synovium.
Since NOXs are less expressed, mitochondria might be the primary generators of ROS in the RA synovium. A component of ROS, O2·- may interact with nitric oxide (NO), which is abundant in RA [76–81], resulting in the formation of ONOO-. The molecule ONOO- is involved in the activation of the IκB kinase (IKK) through a mechanism that depletes the–SH groups of glutathione [43]. Our analysis also indicated an increased degradation of glutathione, possibly by the enzyme LAP3. Active IKK degrades its substrate IκB, resulting in the translocation of NF-κB to the nucleus, where it induces the expression of inflammatory mediators, such as TNF, IL-1β and iNOS. The enzyme iNOS catalyses the production of NO from arginine and hence might further potentiate the production of ONOO- (Fig 5). Furthermore, damaged mitochondria can release molecules, called the damage-associated molecular patterns (DAMPs), which contribute to inflammation [82]. Taken together, mitochondrial production of ROS and their involvement in NF-κB activation, increased glutathione degradation, and DAMPs released from damaged mitochondria suggest that this cell organelle is associated with the induction and maintenance of inflammation in the RA synovium.
Many researchers had created PPI networks for RA at the level of a cell. For instance, in order to identify key molecules, earlier we had created a PPI network for cytokine signalling in RA [83]. Similarly, few studies had identified highly connected regions, ego networks and key genes in PPI networks in RA and other diseases [84–86]. Another study used PPI for determining the efficacy of leflunomide and ligustrazine drugs in the treatment of RA [87]. In contrast, in this study, we created, for the first time, a PPI network that is specific to RA synovial mitochondria and identified hub proteins. This network was created using reliable interactions from seven databases and gene expression data from six open-source microarray datasets. The network has 665 genes, of which 131 are DEGs.
The GO analysis of DEGs has given enriched processes, functions, cell components and KEGG pathways. ‘Oxidoreductase activity’ is the highest enriched molecular function. In general, several metabolic mechanisms, including OxPhos and pathways related to nucleotide, amino acid and glutathione seem to be affected. The analysis also identified the enriched oxidative stress, electron carrier activity, membrane permeability and apoptosis-related GO terms with DEGs.
The analysis of six microarray datasets gave 208 mitochondrial DEGs. Among them, the up-regulation of IFI27, which is involved in cytokine-mediated apoptosis, is in agreement with other RA studies [88]. The enzyme PDK1 was up-regulated. It is involved in hypoxia - and oxidative stress-mediated apoptosis, and is known to induce Akt pathway in human mast cells—which are abundantly seen in the RA synovium [89]. This enzyme also induces cell invasion and secretion of IL-1β and IL-6 in a ribosomal S6 kinase (RSK2)-dependent TNF pathway in FLS [90]. AKR1B10 was down-regulated in three datasets. There is evidence that hypoxia induces the down-regulation of this gene in RA and healthy FLS [91]. EFHD1, which is a calcium-binding protein and known to promote cell death, was down-regulated [92]. The enzyme ACOT7, which increases the concentration of free fatty acids (FFA) by hydrolysing the acyl-CoA thioester of long-chain fatty acids, such as palmitoyl-CoA, was up-regulated in five datasets. This enzyme could be implicated in the remodelling of membrane phospholipids [93]. UCP2 uncouples OxPhos pathway and has been identified a candidate risk gene for RA by a whole genome association study [94]. This protein, which gets activated by H2O2 and O2,-, was up-regulated in the microarray data. The uncoupling of OxPhos by UCP2 leads to the dissipation of mitochondrial membrane potential gradient. Moreover, FFAs produced by ACOT7 are likely to play a role in the uncoupling by increasing the production of ROS [95].
A high ratio of kynurenine/tryptophan, which was observed in sera of RA patients, is needed for the kynurenine pathway for its role in anti-inflammation [96–97]. In our analysis, since KMO, which catalyses the hydroxylation of kynurenine to 3-hydroxykynurenine, was up-regulated in five datasets; it may deplete the concentration of kynurenine thereby impairing its activity in controlling inflammation. Further, this might enhance the generation of free radicals [57]. The analysis also identified the disruption of OxPhos։ 11 subunits of OxPhos complexes were DEGs, of which seven were exclusively down-regulated, clearly indicating the negative regulation of OxPhos. This may further be involved in ROS production. Collectively, these results suggest that free radicals, negative regulation of OxPhos, and metabolic processes are highly active in RA synovial mitochondria.
Increased demand for ATP production on the respiratory machinery, and NOX enzymes are usually involved in ROS generation. But, the down-regulation of NOXs provides a strong support of ROS production by mitochondria rather than by the enzymes in the RA synovium. The observed down-regulation of the detoxifiers of ROS, such as CAT, GPX4 and SOD1 is consistent with other RA studies [98–100]. Since there is no up-regulation of any of these enzymes in any of the six microarray datasets, the detoxification of free radicals might be impaired. Further, LAP3, which degrades glutathione and plays a crucial role in cartilage and bone erosion, was up-regulated [101]. Apart from this, CASP8, an apoptotic caspase and a risk gene for RA, was up-regulated. But, in RA FLS, it is known to induce the activation of pro-inflammatory NF-κB and AP-1 transcription factors rather than apoptosis [45]. However, apoptosis may happen by a ROS-mediated increase in intracellular Ca2+ levels, preferably through the SYK, BTK, BLNK and PLCγ2-mediated pathway. But it may nevertheless be insufficient to limit synovial hyperplasia. Thus, our results bring together enhanced oxidative stress and intrinsic apoptosis, with the up-regulation of processes involved in the generation of free radicals and with their impaired detoxification. These can be envisaged as a potential therapeutic strategy for RA. With regard to the network analysis, we show that UQCR10, MRPL4 and NDUFV2 are the three top hubs of the PPI network. UQCR10 belongs to the complex III, which is the middle segment of OxPhos, and is connected to 50 neighbours, of which nine were DEGs. This suggests a key role for this gene in the OxPhos pathway in RA synoviocytes. MRPL4, a part of the large 39S subunit of the mitochondrial ribosome, forms 49 PPI with neighbours. The complex I protein NDUFV2 is connected to 49 neighbours, including eight DEGs. Another complex I protein NDUFS6 has the maximum number of neighbour DEGs in the network. Since these proteins are the top hubs and have neighbour DEGs, it would be interesting to elucidate their roles in RA.
In this study, using publicly available microarray data, we discussed the roles of nDNA encoded proteins in RA mitochondrial dysfunction. Although the literature provides support to some of the inferences we made in the study, they need to be validated using experiments on cell lines or laboratory animals or in clinical studies. The expression levels of mtDNA encoded genes that were not probed by the microarray platforms could not be assessed in this study. Further, as the composition of cell types is different in both the healthy and RA synovial tissues, the changes in gene expression between them may be a manifestation of the respective cell types present in them.
In conclusion, our study maximises the use of PPI and microarray data for studying mitochondrial dysfunction in RA. Identifying a set of nDNA encoded mitochondrial proteins implicated in the dysregulation of pathways and processes associated with this organelle in the RA synovium was the idea behind this study. Analysing microarray data for identifying DEGs and understanding their likely functions in synovial mitochondria is highly informative in this context. Using DEGs, we identified the processes pertaining to the generation of free radicals and their impaired detoxification. The study also reports the possible occurrence of the mitochondrion-mediated intrinsic apoptotic pathway in RA. We also, in particular, highlighted the roles of DEGs in the remodelling of membrane lipids, uncoupling electron transport and ATP synthesis, and amino acid and nucleotide metabolism in RA. We also proposed a model that links mitochondrial dysfunction to inflammation in RA by collating information from the literature. These insights suggest several new routes for research into the role of mitochondria in RA. Particularly, oxidative stress and intrinsic apoptotic pathways may become attractive candidates for new therapeutic interventions. However, our strategy herein was to develop a proof-of-principle method for studying mitochondrial dysfunction by integrating gene expression, PPI, gene ontology and network analysis. Even though literature search has provided the possible implications for the study findings in RA mitochondrial dysfunction, their additional validation in experimental settings is needed.
S1 File [xgmml] File for visualising the network using the Cytoscape tool.
S1 Fig [tif] A Venn diagram showing the translocation of the mitochondrial network genes to different parts of mitochondria.
S2 Fig [tif] Hierarchical clustering of RA and control samples based on the gene expression of selected DEGs () in GSE7307.
S3 Fig [tif] Hierarchical clustering of RA and control samples based on the gene expression of selected DEGs () in GSE55235.
S4 Fig [tif] Hierarchical clustering of RA and control samples based on the gene expression of selected DEGs () in GSE55457.
S5 Fig [tif] Hierarchical clustering of RA and control samples based on the gene expression of selected DEGs () in GSE12021 (HGU133A).
S6 Fig [tif] Hierarchical clustering of RA and control samples based on the gene expression of selected DEGs () in GSE12021 (HGU133B).
S7 Fig [tif] Hierarchical clustering of RA and control samples based on the gene expression of selected DEGs () in GSE77298.
S8 Fig [tif] Hierarchical clustering of RA and control samples based on the gene expression of selected DEGs () in GSE12021 (HGU133A) after removing RA samples clustered with control samples.
S9 Fig [tif] Hierarchical clustering of RA and control samples based on the gene expression of selected DEGs () in GSE12021 (HGU133B) after removing RA samples clustered with control samples.
S10 Fig [tiff] Hierarchical clustering of prednisone treated and control samples (whole blood) based on the gene expression of selected DEGs () normalized across samples in GSE77344.
S11 Fig [tiff] Hierarchical clustering of celecoxib treated and control samples (colorectal primary adenocarcinomas)based on the gene expression of selected DEGs () normalized across samples in GSE11237.
S12 Fig [tiff] Hierarchical clustering of pre and post methotrexate treatment samples (early RA synovial biopsy) based on the gene expression of selected DEGs () normalized across samples in GSE45867.
S13 Fig [tif] The distribution of the number of all neighbours and DEG neighbours of all the proteins of the mitochondrial PPI network.
S14 Fig [tif] The subnetwork of UQCR10.
S15 Fig [tif] The subnetwork of NDUFV2.
S16 Fig [tif] The subnetwork of NDUFS6.
S17 Fig [tif] The subnetwork of UQCRQ.
S1 Table [xls] The number of microarray datasets in which the mitochondrial PPI network genes were differentially expressed.
S2 Table [xlsx] The number of microarray datasets in which, the MitoCarta genes, which are not part of the network, were differentially expressed.
S3 Table [xlsx] Number of neighbours for all mitochondrial PPI network proteins.
S4 Table [xls] The significantly enriched biological processes (BP) of the RA synovial mitochondrial DEGs.
1. Fearon U, Canavan M, Biniecka M, Veale DJ. Hypoxia, mitochondrial dysfunction and synovial invasiveness in rheumatoid arthritis. Nat Rev Rheumatol. 2016; 12 : 385–97. doi: 10.1038/nrrheum.2016.69 27225300
2. Medina-Gomez G. Mitochondrial and endocrine function of adipose tissue. Best Pract Res Clin Endocrinol Metab. 2012; 26 : 791–804. doi: 10.1016/j.beem.2012.06.002 23168280
3. Pagano G, Castello G, Pallardó FV. Sjøgren’s syndrome-associated oxidative stress and mitochondrial dysfunction: prospects for chemoprevention trials. Free Radic Res. 2013; 47 : 71–3. doi: 10.3109/10715762.2012.748904 23153390
4. Ouang J, Wu M, Huang C, Cao L, Li G. Overexpression of oxidored-nitro domain containing protein 1 inhibits human nasopharyngeal carcinoma and cervical cancer cell proliferation and induces apoptosis: Involvement of mitochondrial apoptotic pathways. Oncol Rep. 2013; 29 : 79–86. doi: 10.3892/or.2012.2101 23124592
5. Osellame LD, Blacker TS, Duchen MR. Cellular and molecular mechanisms of mitochondrial function. Best Pract Res Clin Endocrinol Metab. 2012; 26 : 711–23. doi: 10.1016/j.beem.2012.05.003 23168274
6. Legault JT, Strittmatter L, Tardif J, Sharma R, Tremblay-Vaillancourt V, Aubut C, et al. A metabolic signature of mitochondrial dysfunction revealed through a monogenic form of Leigh syndrome. Cell Rep. 2015; 13(5): 981–89. doi: 10.1016/j.celrep.2015.09.054 26565911
7. Bick AG, Wakimoto H, Kamer KJ, Sancak Y, Goldberger O, Axelsson A, et al. Cardiovascular homeostasis dependence on MICU2, a regulatory subunit of the mitochondrial calcium uniporter. Proc Natl Acad Sci. 2017; 114(43): e9096–9104. doi: 10.1073/pnas.1711303114 29073106
8. Lake NJ, Webb BD, Stroud DA, Richman TR, Ruzzenente B, Compton AG, et al. Biallelic mutations in MRPS34 lead to instability of the small mitoribosomal subunit and Leigh syndrome. Am J Hum Genet. 2017; 101(2): 239–54. doi: 10.1016/j.ajhg.2017.07.005 28777931
9. Dai N, Zhao L, Wrighting D, Krämer D, Majithia A, Wang Y, et al. IGF2BP2/IMP2-deficient mice resist obesity through enhanced translation of Ucp1 mRNA and other mRNAs encoding mitochondrial proteins. Cell Metab. 2015; 21(4): 609–21. doi: 10.1016/j.cmet.2015.03.006 25863250
10. Tucker EJ, Hershman SG, Kohrer C, Belcher-Timme CA, Patel J, Goldberger OA, et al. Mutations in MTFMT underlie a human disorder of formylation causing impaired mitochondrial translation. Cell Metab. 2011; 14(3): 428–34. doi: 10.1016/j.cmet.2011.07.010 21907147
11. Fassone E, Duncan AJ, Taanman JW, Pagnamenta AT, Sadowski MI, Holand T, et al. FOXRED1, encoding an FAD-dependent oxidoreductase complex-I-specific molecular chaperone, is mutated in infantile-onset mitochondrial encephalopathy. Hum Mol Genet. 2010; 19(24): 4837–47. doi: 10.1093/hmg/ddq414 20858599
12. Gohil VM, Nilsson R, Belcher-Timme CA, Luo B, Root DE, Mootha VK. Mitochondrial and nuclear genomic responses to loss of LRPPRC expression. J Biol Chem. 2010; 285(18): 13742–7. doi: 10.1074/jbc.M109.098400 20220140
13. Kamer KJ, Mootha VK. MICU1 and MICU2 play nonredundant roles in the regulation of the mitochondrial calcium uniporter. EMBO Rep. 2014; 15(3): 299–07. doi: 10.1002/embr.201337946 24503055
14. Sancak Y, Markhard AL, Kitami T, Kovács-Bogdán E, Kamer KJ, Udeshi ND, et al. EMRE is an essential component of the mitochondrial calcium uniporter complex. Science. 2013; 342(6164): 1379–82. doi: 10.1126/science.1242993 24231807
15. Plovanich M, Bogorad RL, Sancak Y, Kamer KJ, Strittmatter L, Li AA, et al. MICU2, a paralog of MICU1, resides within the mitochondrial uniporter complex to regulate calcium handling. Plos One. 2013; 8(2): e55785. doi: 10.1371/journal.pone.0055785 23409044
16. Calvo SE, Clauser KR, Mootha VK. MitoCarta2.0: an updated inventory of mammalian mitochondrial proteins. Nucleic Acids Res. 2016; 44: D1251–7. doi: 10.1093/nar/gkv1003 26450961
17. McInnes IB, Schett G. Cytokines in the pathogenesis of rheumatoid arthritis. Nat Rev Immunol. 2007; 7(6): 429–42. doi: 10.1038/nri2094 17525752
18. Moelants EA, Mortier A, Van Damme J, Proost P. Regulation of TNF-alpha with a focus on rheumatoid arthritis. Immunol Cell Biol. 2013; 91 : 393–401. doi: 10.1038/icb.2013.15 23628802
19. Smith MD. The normal synovium. Open Rheumatol J. 2011; 5 : 100–06. doi: 10.2174/1874312901105010100 22279508
20. Bromley M, Woolley DE. Histopathology of the rheumatoid lesion. Identification of cell types at sites of cartilage erosion. Arthritis Rheum. 1984; 27(8): 857–63. doi: 10.1002/art.1780270804 6466394
21. Tak PP, Bresnihan B. The pathogenesis and prevention of joint damage in rheumatoid arthritis: advances from synovial biopsy and tissue analysis. Arthritis Rheum. 2000; 43(12): 2619–33. doi: 10.1002/1529-0131(200012)43 : 12<2619::AID-ANR1>3.0.CO;2-V 11145019
22. Ng CT, Biniecka M, Kennedy A, McCormick J, FitzGerald O, Bresnihan B, et al. Synovial tissue hypoxia and inflammation in vivo. Ann Rheum Dis. 2010; 69(7): 1389–95. doi: 10.1136/ard.2009.119776 20439288
23. Hildeman DA, Mitchell T, Kappler J, Marrack P. T cell apoptosis and reactive oxygen species. J Clin Invest. 2003; 111 : 575–81. doi: 10.1172/JCI18007 12618509
24. Kennedy A, Ng CT, Chang TC, Biniecka M, O’Sullivan JN, Heffernan E, et al. Tumor necrosis factor blocking therapy alters joint inflammation and hypoxia. Arthritis Rheum. 2011; 63(4): 923–32. doi: 10.1002/art.30221 21225682
25. Hitchon CA, El-Gabalawy HS. Oxidation in rheumatoid arthritis. Arthritis Res Ther. 2004; 6 : 265–78. doi: 10.1186/ar1447 15535839
26. Madamanchi NR, Runge MS. Mitochondrial dysfunction in atherosclerosis. Circ Res. 2007; 100 : 460–73. doi: 10.1161/01.RES.0000258450.44413.96 17332437
27. Kokoszka JE, Coskun P, Esposito LA, Wallace DC. Increased mitochondrial oxidative stress in the Sod2 (+/-) mouse results in the age-related decline of mitochondrial function culminating in increased apoptosis. Proc Natl Acad Sci. 2001; 98(5): 2278–83. doi: 10.1073/pnas.051627098 11226230
28. Hajizadeh S, DeGroot J, TeKoppele JM, Tarkowski A, Collins LV. Extracellular mitochondrial DNA and oxidatively damaged DNA in synovial fluid of patients with rheumatoid arthritis. Arthritis Res Ther. 2003; 5: R234–40. doi: 10.1186/ar787 12932286
29. Bashir S, Harris G, Denman MA, Blake DR, Winyard PG. Oxidative DNA damage and cellular sensitivity to oxidative stress in human autoimmune diseases. Ann Rheum Dis. 1993; 52 : 659–66. doi: 10.1136/ard.52.9.659 8239761
30. Lemarechal H, Allanore Y, Chenevier-Gobeaux C, Kahan A, Ekindjian OG, Borderie D, et al. Serum protein oxidation in patients with rheumatoid arthritis and effects of infliximab therapy. Clin Chim Acta. 2006; 372 : 147–53. doi: 10.1016/j.cca.2006.04.002 16716286
31. Dabbagh AJ, Trenam CW, Morris CJ, Blake DR. Iron in joint inflammation. Ann Rheum Dis. 1993; 52 : 67–73. doi: 10.1136/ard.52.1.67 8427520
32. Dai L, Lamb DJ, Leake DS, Kus ML, Jones HW, Morris CJ, et al. Evidence for oxidized low density lipoprotein in synovial fluid from rheumatoid arthritis patients. Free Radic Res. 2000; 32 : 479–86. doi: 10.1080/10715760000300481 10798713
33. Rowley D, Gutteridge JM, Blake D, Farr M, Halliwell B. Lipid peroxidation in rheumatoid arthritis: thiobarbituric acid-reactive material and catalytic iron salts in synovial fluid from rheumatoid patients. Clin Sci (Lond). 1984; 66(6): 691–5.
34. Grootveld M, Henderson EB, Farrell A, Blake DR, Parkes HG, Haycock P. Oxidative damage to hyaluronate and glucose in synovial fluid during exercise of the inflamed rheumatoid joint. Detection of abnormal low-molecular-mass metabolites by proton-n.m.r. spectroscopy. Biochem J. 1991; 273 : 459–67. doi: 10.1042/bj2730459 1991041
35. Henrotin YE, Bruckner P, Pujol JPL. The role of reactive oxygen species in homeostasis and degradation of cartilage. Osteoarthritis and Cartilage. 2003; 11 : 747–55. 13129694
36. Ushio-Fukai M, Tang Y, Fukai T, Dikalov SI, Ma Y, Fujimoto M, et al. Novel role of gp91phox-containing NAD(P)H oxidase in vascular endothelial growth factor-induced signaling and angiogenesis. Circ Res. 2002; 91 : 1160–67. doi: 10.1161/01.res.0000046227.65158.f8 12480817
37. Cerhan JR, Saag KG, Merlino LA, Mikuls TR, Criswell LA. Antioxidant micronutrients and risk of rheumatoid arthritis in a cohort of older women. Am J Epidemol. 2003; 157 : 345–54.
38. Heliövaara M, Knekt P, Aho K, Aaran R-K, Alfthan G, Aromaa A. Serum antioxidants and risk of rheumatoid arthritis. Ann Rheum Dis. 1994; 53(1): 51–53. doi: 10.1136/ard.53.1.51 8311556
39. Hagfors L, Leanderson P, Skoldstam L, Andersson J, Johansson G. Antioxidant intake, plasma antioxidants and oxidative stress in a randomized, controlled, parallel, Mediterranean dietary intervention study on patients with rheumatoid arthritis. Nutr J. 2003; 2 : 5. doi: 10.1186/1475-2891-2-5 12952549
40. Bae SC, Kim SJ, Sung MK. Inadequate antioxidant nutrient intake and altered plasma antioxidant status of rheumatoid arthritis patients. J Am Coll Nutr. 2003; 22(4): 311–5. doi: 10.1080/07315724.2003.10719309 12897046
41. Paredes S, Girona J, Hurt-Camejo E, Vallvé JC, Olivé S, Heras M, et al. Antioxidant vitamins and lipid peroxidation in patients with rheumatoid arthritis: association with inflammatory markers. J Rheumatol. 2002; 29(11): 2271–7. 12415581
42. Mulherin DM, Thurnham DI, Situnayake RD. Glutathione reductase activity, riboflavin status, and disease activity in rheumatoid arthritis. Ann Rheum Dis. 1996; 55 : 837–914. doi: 10.1136/ard.55.11.837 8976642
43. Filippin LI, Vercelino R, Marroni NP, Xavier RM. Redox signalling and the inflammatory response in rheumatoid arthritis. Clin Exp Immunol. 2008; 152(3): 415–22. doi: 10.1111/j.1365-2249.2008.03634.x 18422737
44. Li H, Wan A. Apoptosis of rheumatoid arthritis fibroblast-like synoviocytes: possible roles of nitric oxide and the thioredoxin 1. Mediators Inflamm. 2013; 2013(3): 953462.
45. Palao G, Santiago B, Galindo M, Rullas J, Alcamí J, Ramirez JC, et al. Fas activation of a proinflammatory program in rheumatoid synoviocytes and its regulation by FLIP and caspase 8 signaling. Arthritis Rheum. 2006; 54(5): 1473–81. doi: 10.1002/art.21768 16646028
46. Chatr-Aryamontri A, Breitkreutz B-J, Oughtred R, Boucher L, Heinicke S, Chen D, et al. The BioGRID interaction database: 2015 update. Nucleic Acids Res. 2015; 43: D470–8. doi: 10.1093/nar/gku1204 25428363
47. Orchard S, Ammari M, Aranda B, Breuza L, Briganti L, Broackes-Carter F, et al. The MIntAct project—IntAct as a common curation platform for 11 molecular interaction databases. Nucleic Acids Res. 2014; 42: D358–63. doi: 10.1093/nar/gkt1115 24234451
48. Chatr-Aryamontri A, Ceol A, Palazzi LM, Nardelli G, Schneider MV, Castagnoli L, et al. MINT: the molecular interaction database. Nucleic Acids Res. 2007; 35: D572–74. doi: 10.1093/nar/gkl950 17135203
49. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015; 43(D1): D447–52.
50. Prasad TSK, Goel R, Kandasamy K, Keerthikumar S, Kumar S, Mathivanan S, et al. Human protein reference database—2009 update. Nucleic Acids Res. 2009; 37: D767–72. doi: 10.1093/nar/gkn892 18988627
51. Salwinski L, Miller CS, Smith AJ, Pettit FK, Bowie JU, Eisenberg D. The database of interacting proteins: 2004 update. Nucleic Acids Res. 2004; 32: D449–51. doi: 10.1093/nar/gkh086 14681454
52. Bossi A, Lehner B. Tissue specificity and the human protein interaction network. Mol Syst Biol. 2009; 5 : 260. doi: 10.1038/msb.2009.17 19357639
53. Huber R, Hummert C, Gausmann U, Pohlers D, Koczan D, Guthke R, et al. Identification of intra-group, inter-individual, and gene-specific variances in mRNA expression profiles in the rheumatoid arthritis synovial membrane. Arthritis Res Ther. 2008; 10(14): R98.
54. Woetzel D, Huber R, Kupfer P, Pohlers D, Pfaff M, Driesch D, et al. Identification of rheumatoid arthritis and osteoarthritis patients by transcriptome-based rule set generation. Arthritis Res Ther. 2014; 16: R84. doi: 10.1186/ar4526 24690414
55. Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, et al. STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013; 41: D808–15. doi: 10.1093/nar/gks1094 23203871
56. Rivals I, Personnaz L, Taing L, Marie-Claude P. Enrichment or depletion of a GO category within a class of genes: which test? Bioinformatics. 2007; 23(4): 401–07. doi: 10.1093/bioinformatics/btl633 17182697
57. Dantzer R, O’Connor JC, Lawson MA, Kelley KW. Inflammation-associated depression: from serotonin to kynurenine. Psychoneuroendocrinology. 2011; 36(3): 426–36. doi: 10.1016/j.psyneuen.2010.09.012 21041030
58. Kaur H, Ganguli D, Bachhawat AK. Glutathione degradation by the alternative pathway (DUG pathway) in Saccharomyces cerevisiae is initiated by (Dug2p-Dug3p)2 complex, a novel glutamine amidotransferase (GATase) enzyme acting on glutathione. J Biol Chem. 2012; 287(12): 8920–31. doi: 10.1074/jbc.M111.327411 22277648
59. Kim JW, Tchernyshyov I, Semenza GL, Dang CV. HIF-1-mediated expression of pyruvate dehydrogenase kinase: A metabolic switch required for cellular adaptation to hypoxia. Cell Metab. 2006; 3(3): 177–85. doi: 10.1016/j.cmet.2006.02.002 16517405
60. Rosebeck S, Leaman DW. Mitochondrial localization and pro-apoptotic effects of the interferon-inducible protein ISG12a. Apoptosis. 2008; 13(4): 562–72. doi: 10.1007/s10495-008-0190-0 18330707
61. Arsenijevic D, Onuma H, Pecqueur C, Raimbault S, Manning BS, Miroux B, et al. Disruption of the uncoupling protein-2 gene in mice reveals a role in immunity and reactive oxygen species production. Nat Genet. 2000; 26(4): 435–9. doi: 10.1038/82565 11101840
62. Jin DY, Chae HZ, Rhee SG, Jeang KT. Regulatory role for a novel human thioredoxin peroxidise in NF-kappaB activation. J Biol Chem. 1997; 272(49): 30952–61. doi: 10.1074/jbc.272.49.30952 9388242
63. Wang X, Wang L, Wang X, Sun F, Wang CC. Structural insights into the peroxidase activity and inactivation of human peroxiredoxin 4. Biochem J. 2012; 441(1): 113–8. doi: 10.1042/BJ20110380 21916849
64. Hartmann B, Wai T, Hu H, MacVicar T, Musante L, Fischer-Zirnsak B, et al. Homozygous YME1L1 mutation causes mitochondriopathy with optic atrophy and mitochondrial network fragmentation. Elife. 2016; 5:e16078. doi: 10.7554/eLife.16078 27495975
65. Stiburek L, Cesnekova J, Kostkova O, Fornuskova D, Vinsova K, Wenchich L, et al. YME1L controls the accumulation of respiratory chain subunits and is required for apoptotic resistance, cristae morphogenesis, and cell proliferation. Mol Biol Cell. 2012; 23(6): 1010–23. doi: 10.1091/mbc.E11-08-0674 22262461
66. Vishnivetskaya GB, Skrinskaya JA, Seif I, Popova NK. Effect of MAO A deficiency on different kinds of aggression and social investigation in mice. Aggress Behav. 2007; 33(1): 1–6. doi: 10.1002/ab.20161 17441000
67. Brunner HG, Nelen M, Breakefield XO, Ropers HH, van Oost BA. Abnormal behavior associated with a point mutation in the structural gene for monoamine oxidase A. Science. 1993; 262(5133): 578–80. doi: 10.1126/science.8211186 8211186
68. Gupta V, Khan AA, Sasi BK, Mahapatra NR. Molecular mechanism of monoamine oxidase A gene regulation under inflammation and ischemia-like conditions: key roles of the transcription factors GATA2, Sp1 and TBP. J Neurochem. 2015; 134(1): 21–38. doi: 10.1111/jnc.13099 25810277
69. Okada Y, Wu D, Trynka G, Raj T, Terao C, Ikan K, et al. Genetics of rheumatoid arthritis contributes to biology and drug discovery. Nature. 2014; 506 : 376–81. doi: 10.1038/nature12873 24390342
70. Shannon CP, Balshaw R, Chen V, Hollander Z, Toma M, McManus BM, et al. Enumerateblood–an R Package to Estimate the Cellular Composition of Whole Blood from Affymetrix Gene ST Gene Expression Profiles. BMC Genomics. 2017; 18 : 43. doi: 10.1186/s12864-016-3460-1 28061752
71. Auman JT, Church R, Lee SY, Watson MA, Fleshman JW, Mcleod HL. Celecoxib Pre-Treatment in Human Colorectal Adenocarcinoma Patients Is Associated with Gene Expression Alterations Suggestive of Diminished Cellular Proliferation. Eur J Cancer. 2008; 44(12): 1754–60. doi: 10.1016/j.ejca.2008.05.010 18653328
72. Ducreux J, Durez P, Galant C, Nzeusseu Toukap A, Van den Eynde B, Houssiau FA, et al. Global Molecular Effects of Tocilizumab Therapy in Rheumatoid Arthritis Synovium. Arthritis Rheumatol. 2014; 66(1): 15–23. doi: 10.1002/art.38202 24449571
73. Barabási A-L, Gulbahce N, Loscalzo J. Network medicine: a network-based approach to human disease. Nat Rev Genet. 2011; 12 : 56–58. doi: 10.1038/nrg2918 21164525
74. Qin S, Chock PB. Bruton’s tyrosine kinase is essential for hydrogen peroxide-induced calcium signaling. Biochemistry. 2001; 40(27): 8085–91. doi: 10.1021/bi0100788 11434777
75. Olofsson MH, Havelka AM, Brnjic S, Shoshan MC, Linder S. Charting calcium-regulated apoptosis pathways using chemical biology: role of calmodulin kinase II. BMC Chem Biol. 2008;8 : 2. doi: 10.1186/1472-6769-8-2 18673549
76. Jang D, Murrell GAC. Nitric oxide in arthritis. Free Radic Biol Med. 1998; 24(9): 1511–9. doi: 10.1016/s0891-5849(97)00459-0 9641270
77. van’t Hof RJ, Hocking L, Wright PK, Ralston SH. Nitric oxide is a mediator of apoptosis in the rheumatoid joint. Rheumatology. 2000; 39 : 1004–8. doi: 10.1093/rheumatology/39.9.1004 10986306
78. Sakurai H, Kohsaka H, Liu MF, Higashiyama H, Hirata Y, Kanno K, et al. Nitric oxide production and inducible nitric oxide synthase expression in inflammatory arthritides. J Clin Invest. 1995; 96(5): 2357–63. doi: 10.1172/JCI118292 7593623
79. Dey P, Panga V, Raghunathan S. A cytokine signaling network for the regulation of inducible nitric oxide synthase expression in rheumatoid arthritis. Plos One. 2016; 11(9): e0161306. doi: 10.1371/journal.pone.0161306 27626941
80. Bauerová K, Bezek Š. Role of reactive oxygen and nitrogen species in etiopathogenesis of rheumatoid arthritis. Gen Physiol Biophys. 1999; 18 : 15–20.
81. Phillips DC, Irundika Dias HK, Kitas GD, Griffiths HR. Aberrant reactive oxygen and nitrogen species generation in rheumatoid arthritis (RA): causes and consequences for immune function, cell survival, and therapeutic intervention. Antioxid Redox Signal. 2010; 12(6): 743–85. doi: 10.1089/ars.2009.2607 19686039
82. Nakahira K, Hisata S, Choi AMK. The roles of mitochondrial damage-associated molecular patterns in diseases. Antioxid Redox Signal. 2015; 23 : 1329–50. doi: 10.1089/ars.2015.6407 26067258
83. Panga V, Raghunathan S. A cytokine protein-protein interaction network for identifying key molecules in rheumatoid arthritis. Plos One. 2018; 13(6): e0199530. doi: 10.1371/journal.pone.0199530 29928007
84. Zhu N, Hou J, Wu Y, Li G, Liu J, Ma G, et al. Identification of key genes in rheumatoid arthritis and osteoarthritis based on bioinformatics analysis. Medicine. 2018; 97(22): e10997. doi: 10.1097/MD.0000000000010997 29851858
85. Zhou WZ, Miao LG, Yuan H. Identification of significant ego networks and pathways in rheumatoid arthritis. J Can Res Ther. 2018; 14 : 1024–8.
86. Sun Z, Wang W, Yu D, Mao Y. Differentially expressed genes between systemic sclerosis and rheumatoid arthritis. Hereditas. 2019; 156 : 17. doi: 10.1186/s41065-019-0091-y 31178673
87. Zhang C, Guan D, Jiang M, Liang C, Li L, Zhao N, et al. Efficacy of leflunomide combined with ligustrazine in the treatment of rheumatoid arthritis: prediction with network pharmacology and validation in a clinical trial. Chin Med. 2019; 14 : 26. doi: 10.1186/s13020-019-0247-8 31388350
88. Toro-Domínguez D, Carmona-Sáez P, Alarcón-Riquelme ME. Shared signatures between rheumatoid arthritis, systemic lupus erythematosus and Sjögren’s syndrome uncovered through gene expression meta-analysis. Arthritis Res Ther. 2014; 16 : 489. doi: 10.1186/s13075-014-0489-x 25466291
89. Sawamukai N, Saito K, Yamaoka K, Nakayamada S, Ra C, Tanaka Y. Leflunomide inhibits PDK1/Akt pathway and induces apoptosis of human mast cells. J Immunol. 2007; 179(10): 6479–84. doi: 10.4049/jimmunol.179.10.6479 17982036
90. Sun C, Sun Y, Jiang D, Bao G, Zhu X, Xu D, et al. PDK1 promotes the inflammatory progress of fibroblast-like synoviocytes by phosphorylating RSK2. Cell Immunol. 2017; 315 : 27–33. doi: 10.1016/j.cellimm.2016.10.007 28314444
91. Del Rey MJ, Izquierdo E, Usategui A, Gonzalo E, Blanco FJ, Acquadro F, et al. The transcriptional response of normal and rheumatoid arthritis synovial fibroblasts to hypoxia. Arthritis Rheum. 2010; 62(12): 3584–94. doi: 10.1002/art.27750 20848564
92. Tominaga M, Kurihara H, Honda S, Amakawa G, Sakai T, Tomooka Y. Molecular characterization of mitocalcin, a novel mitochondrial Ca2+-binding protein with EF-hand and coiled-coil domains. J Neurochem. 2006; 96(1): 292–304. doi: 10.1111/j.1471-4159.2005.03554.x 16336229
93. Wall VZ, Barnhart S, Kramer F, Kanter JE, Vivekanandan-Giri A, Pennathur S, et al. Inflammatory stimuli induce acyl-CoA thioesterase 7 and remodeling of phospholipids containing unsaturated long (≥C20)-acyl chains in macrophages. J Lipid Res. 2017; 58 : 1174–85. doi: 10.1194/jlr.M076489 28416579
94. Tamiya G, Shinya M, Imanishi T, Ikuta T, Makino S, Okamoto K, et al. Whole genome association study of rheumatoid arthritis using 27039 microsatellites. Hum Mol Genet. 2005; 14(16): 2305–21. doi: 10.1093/hmg/ddi234 16000323
95. Fariss MW, Chan CB, Patel M, Van Houten B, Orrenius S. Role of mitochondria in toxic oxidative stress. Mol Interv. 2005; 5(2): 94–111. doi: 10.1124/mi.5.2.7 15821158
96. Schroecksnadel K, Kaser S, Ledochowski M, Neurauter G, Mur E, Herold M, et al. Increased degradation of tryptophan in blood of patients with rheumatoid arthritis. J Rheumatol. 2003; 30 : 1935–9. 12966593
97. Kolodziej L. An exploratory study of the interplay between decreased concentration of tryptophan, accumulation of kynurenines, and inflammatory arthritis. IUBMB Life. 2012; 64(12): 983–87. doi: 10.1002/iub.1092 23124849
98. Taysi S, Polat F, Gul M, Sari RA, Bakan E. Lipid peroxidation, some extracellular antioxidants, and antioxidant enzymes in serum of patients with rheumatoid arthritis. Rheumatol Int. 2002; 21(5): 200–4. doi: 10.1007/s00296-001-0163-x 11958437
99. Marklund SL, Bjelle A, Elmqvist LG. Superoxide dismutase isoenzymes of the synovial fluid in rheumatoid arthritis and in reactive arthritides. Ann Rheum Dis. 1986; 45 : 847–51. doi: 10.1136/ard.45.10.847 3789819
100. Ozturk HS, Cimen MY, Cimen OB, Kacmaz M, Durak I. Oxidant/antioxidant status of plasma samples from patients with rheumatoid arthritis. Rheumatol Int. 1999; 19 : 35–37. doi: 10.1007/s002960050097 10651080
101. Vainio U. Leucineaminopeptidase in rheumatoid arthritis. Ann Rheum Dis. 1970; 29 : 434.
Zvýšte si kvalifikáciu online z pohodlia domova
Nemáte účet? Registrujte sa
Zadajte e-mailovú adresu, s ktorou ste vytvárali účet. Budú Vám na ňu zasielané informácie k nastaveniu nového hesla.
