International Journal of Genomics

Volume 2014 (2014), Article ID 452891, 10 pages

http://dx.doi.org/10.1155/2014/452891

Research Article

Coexpression within Integrated Mitochondrial Pathways Reveals Different Networks in Normal and Chemically Treated Transcriptomes

Cong Chen,¹ Tae Kyung Hyun,² Xiao Han,² Zhihui Feng,¹ Yuan Li,¹ Xiaolong Liu,³ and Jiankang Liu¹

¹Center for Mitochondrial Biology and Medicine, The Key Laboratory of Biomedical Information Engineering of Ministry of Education, School of Life Science and Technology and Frontier Institute of Science and Technology, Xi’an Jiaotong University, Xi’an 710049, China

²Division of Applied Life Science (Brain Korea 21-World Class University Program), Plant Molecular Biology and Biotechnology Research Center, Gyeongsang National University, Jinju 660-701, Republic of Korea

³The Liver Center of Fujian Province, Fujian Medical University, Fuzhou 350025, China

Received 22 February 2014; Revised 13 April 2014; Accepted 5 May 2014; Published 24 June 2014

Academic Editor: Graziano Pesole

Copyright © 2014 Cong Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

As energy producers, mitochondria play a pivotal role in multiple cellular processes. Although several lines of evidence suggest that differential expression of mitochondrial respiratory complexes (MRCs) has a significant impact on mitochondrial function, the role of integrated MRCs in the whole coexpression network has yet to be revealed. In this study, we construct coexpression networks based on microarray datasets from different tissues and chemical treatments to explore the role of integrated MRCs in the coexpression network and the effects of different chemicals on the mitochondrial network. By grouping MRCs as one seed target, the hypergeometric distribution allowed us to identify genes that are significantly coexpress with whole MRCs. Coexpression among 46 MRC genes (approximately 78% of MRC genes tested) was significant in the normal tissue transcriptome dataset. These MRC genes are coexpressed with genes involved in the categories “muscle system process,” “metabolic process,” and “neurodegenerative disease pathways,” whereas, in the chemically treated tissues, coexpression of these genes mostly disappeared. These results indicate that chemical stimuli alter the normal coexpression network of MRC genes. Taken together, the datasets obtained from the different coexpression networks are informative about mitochondrial biogenesis and should contribute to understanding the side effects of drugs on mitochondrial function.

1. Introduction

Mitochondria are small membrane-enclosed organelles (from 0.5 to 1.0 μM in diameter) found in most eukaryotic cells except mature red blood corpuscles [1]. Mitochondria are the powerhouses of eukaryotic cells and are involved in many cellular processes, including apoptosis; ion homeostasis; and the metabolism of glucose, lipids, and amino acids [2]. ATP, the energy currency of cell, is the final product of the respiratory chain/oxidative phosphorylation system, which consists of five protein complexes (complexes I–V) localized to the inner mitochondrial membrane [3]. Mitochondrial defects including mitochondrial DNA (mtDNA) mutations, altered expression and activity of respiratory chain subunits and glycolytic enzymes, and decreased oxidation of NADH-linked substrates have been suspected to play an important role in the development and progression of diseases, such as certain neurodegenerative diseases, diabetes, Leigh’s disease, and cancer [4–7]. These clinical phenotypes are direct evidence reflecting the essential function of mitochondria. Mitochondrial genetic disorders are caused by defects in nuclear or mtDNA that affect the expression of the mtDNA-encoded mitochondrial respiratory complexes (MRCs) and the biosynthesis of the mtDNA-encoded polypeptides [8]. Mutations in genes required for mtDNA maintenance, expression, and replication regulate genetic disorders, indicating that differential expression of MRCs and related genes has a significant impact on mitochondrial dysfunction [7, 9]. Therefore, systematic analysis of nuclear and mitochondrial gene expression in the context of well-defined disease models should provide insight into the interaction of gene regulatory networks with MRCs, improving our understating of mitochondrial disorders.

Coexpression analysis using transcriptome datasets generated by high-throughput microarray transcript profiling produces correlations that have often been considered to imply functional relationships [10, 11]. A strong correlation among transcripts for MRC components has been found by this type of coexpression analysis in plants. In the case of plant MRC genes, it has been shown that genes belonging to MRCs are clustered into the same coexpression group [10]. Similarly, several mtDNA-encoded mitochondrial genes form a small cluster with a nuclear-encoded mitochondrial gene module and the glycolysis module [12]. Coexpression analysis has also implicated several unannotated genes in cancer and mitochondrial complex I disease [13], indicating that coexpression analysis is a useful tool not only for understanding many diseases at the molecular level but also for identification of novel candidate genes involved in mitochondria-related diseases. Although several studies have demonstrated the power of coexpression analysis, few have exploited MRCs as an integrated component for analyzing the coexpression network.

To further investigate the role of integrated MRCs in the whole coexpression network, we determined the coexpression networks of normal and chemically treated human tissues by analysis of Pearson correlation factors. In the coexpression network under normal conditions, we found that the MRCs are almost fully self-connected. This self-connection indicates that whole MRCs might play a role similar to that of single genes in the coexpression network. Using the hypergeometric distribution, we considered genes with a value less than 10⁻⁸ to be coexpression “friends” with MRCs. Candidate functions for these “friends” were determined through enrichment analysis, using Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. We then explored the coexpression network between “friends” and MRCs in both normal and chemically treated tissues. The systematic coexpression network of genes interacting with whole MRCs identifies candidates that potentially participate in mitochondrial biogenesis and could serve as targets for future therapeutic interventions aimed at modulating mitochondrial function.

2. Methods

2.1. Selection of Datasets and Construction of Coexpression Networks

To create coexpression “friends” with MRCs, we first constructed a genome-wide coexpression networks using two different microarray datasets. Expression datasets for 65 human tissues were downloaded from the COXPRESdb website (http://coxpresdb.jp/). Transcription profiles (E-MTAB-798) of human hepatocytes treated with 130 chemical compounds including drugs such as acetaminophen, aspirin, rifampicin, metformin hydrochloride, simvastatin, and tamoxifen citrate were obtained from EMBL-EBI (http://www.ebi.ac.uk/arrayexpress/). Coexpression networks were constructed as described in Azuaje [14]. The Spearman coexpression coefficient, , was calculated for each pair of genes, and all gene pairs with were defined as gene-gene associations in the network. In the coexpression network, the nodes represent genes and the edges represent the connection with coefficient ≥0.3.

2.2. Analysis of Coexpression Significance

The MRC is designated as NADH-Coenzyme Q reductase (complex I), succinate-CoQ reductase (complex II), ubiquinol-cytochrome c reductase (complex III), cytochrome c oxidase (complex IV), and ATP synthase (complex V). Complexes I, II, III, and IV play as the electron transfer complexes, whereas complex V is known as an enzyme-conserving complex [3]. Since different datasets contain different probes mapping to different gene symbols, we used gene symbols that are present in gene platform file containing 56 MRC genes (33 genes for complex I, 4 genes for complex II, 8 genes for complex III, and 23 genes for complex IV). To determine coexpression significance, 56 MRC genes were selected as an integrated target seed for further analysis. The hypergeometric distribution was used to calculate the connection between the MRC genes and other genes in the whole coexpression network. This discrete probability distribution describes the probability of successes in draws, without replacement, from a finite population of size containing samples. For example, suppose that there are (56) MRC genes among (20,000) genes in the genome. Gene X makes connections with MRCs and connections with the genome. We define genes with a value less than 10⁻⁸ to be coexpression “friends” with MRCs.

2.3. Software Tools

The R platform (http://www.r-project.org/) was used for network generation and statistical calculations. Functional annotation of “friends” was carried out using the web based tool PANTHER (http://www.pantherdb.org/). GO enrichment and KEGG analysis were performed with DAVID (http://david.abcc.ncifcrf.gov/). DAVID was also used to analyze the functional annotations of the gene sets and modules. For pathway enrichment analysis of the MCR coexpressed genes, “GO_BP,” “KEGG_PATHWAY,” and “PANTHER_PATHWAY” were selected. The values and a modified Fisher’s exact test were used to determine the enrichment of gene sets in ontology.

3. Results and Discussion

3.1. Coexpression of MRC Genes across Different Human Tissues

Two functional entities are involved in the generation of ATP by a process called oxidative phosphorylation located in the mitochondrial inner membrane [3]. The first entity is the electron transfer chain historically defined as four complexes (I, II, III, and IV), whereas the second entity is known as the system that phosphorylates ADP to produce ATP [15]. Among these complexes, complex I is the first and largest enzyme complex of the respiratory chain and is directly involved in maintaining cellular reduction-oxidation (redox; NADH/NAD⁺) homeostasis [15]. The mammalian complex I is composed of at least 45 subunits and is the main source of reactive oxygen species, which are implicated in cell signaling, disease, and aging [16]. Its deficiency is the most frequently encountered in mitochondrial disorders [17], and the large number of genes coding for complex I subunits might explain why complex I deficiency is characterized by marked clinical and genetic heterogeneity [18]. Complex II is composed of four nuclear-encoded subunits, whereas complex III is a complex of 11 subunits [19, 20]. Complex II receives electrons via FADH2 and transfers it to complex III thought Coenzyme Q 10. Then electrons are carried by cytochrome c to complex IV, which is composed of 19 subunits. This electron transport is required for the generation of the transmembrane proton gradient in inner mitochondrial membrane which is utilized by complex V to convert ADP to ATP [21]. Defect in any of MRCs leads to impaired ATP production and results in a mitochondrial disease involving abnormality of the central nervous system and eyes, renal, muscle, heat, and haematological system, as well as diverse age-related disorders including cancer and degenerative diseases [21–24]. This indicates that these complexes have a significant impact on mitochondrial function. To investigate the connectivity of expression between MRCs and other genes, we generated coexpression networks using expression datasets for 65 human tissues. Prior genome-wide expression analyses have demonstrated significant coexpression of MRC genes under various physiological conditions in several species [25]. Similarly, our large-scale analysis across different human tissues reveals a coexpression cluster (46 out of 56 MRC genes) significantly enriched in genes belonging to mitochondrial complexes I to IV (Figure 1(a)). Of these, 37 genes belong to complex I, composed of 45 subunits [26]. Mitochondrial complex I uses NADH as a cofactor for electron transfer and translocates protons across the inner mitochondrial membrane [27]. The genes NDUFB5 () and NDUFA7 (), two subunits of complex I, exhibited the lowest values (Supplementary Table S1 available online at http://dx.doi.org/10.1155/2014/452891), indicating that these genes are significantly coexpressed with other MRCs. In many organisms, the complexes I, III, and IV can associate into supercomplexes [28–30]. Among the various types of association, the I + III₂ + IV1–4 supercomplex or the respirasome is one of the most intriguing supercomplexes, because it considered the minimal unit to perform complete respiration from NADH to oxygen [29, 31]. This supercomplex has also been detected by inhibitor titration in bovine mitochondria, suggesting that the two mitochondrial electron transfer complexes specifically interact to form this supercomplex. In our coexpression network (Figure 2(a)), NDUFB5 is coexpressed with SDHC of complex II and UQCRC2 of complex III. In addition, UQCRC1 in mitochondrial complex III is significantly coexpressed with the complex I subunits NDUFA13, NDUFAF1, and NDUFS7. It has been shown that the absence in complex III results in a dramatic loss of complex I in humans, and complex I is necessary for fully assembled complex III [32–34], indicating that supercomplex formation is necessary for assembly and stability of individual components [35]. Taken together, these results suggest that complex I is tightly coexpressed with complex III compared to other complexes and that this coexpression might be required for maintaining the supercomplex.

(a)

(b)

Figure 1: Genes significantly coexpressed with MRCs in normal and chemically treated tissues. (a) Self-connections among mitochondrial respiratory complexes. (b) Coexpressed cellular genes.

(a)

(b)

Figure 2: Genes coexpressed with MRCs in (a) the tissue-based network and (b) the chemical treatment-based network.

3.2. Coexpression of MRC Genes in Chemically Treated Tissues

Drug-induced liver injury is a common side effect of certain pharmaceutical therapies. Drugs can be metabolized into electrophilic chemicals or free radicals, which have direct effects upon mitochondrial proteins. Damage to mitochondrial proteins decreases their affinity for substrates, resulting in mitochondrial dysfunction [36, 37]. Therefore, understanding drug-induced mitochondrial toxicity is critical for the development of safe drugs. To investigate the effect of chemical toxicity on the coexpression network of MRC genes in liver, transcription profiles of human hepatocytes treated with 130 chemical compounds were used as input. Stress-related stimuli induce the remodeling of coexpression networks, resulting in the large-scale alteration of cellular function, involving a shift of resources from growth and metabolism to protection and maintenance [38, 39]. As shown in Figure 1(b), coexpression of MRC genes in treated human hepatocytes was significantly lower than in nontreated human tissues. NDUFAF4, a subunit of complex I, exhibited significant coexpression with other MRCs under normal conditions (; Supplementary Table S1). However, the value between NDUFAF4 and other MRCs increased substantially () after chemical treatment, indicating a major change in the coexpression network. UQCRC2, a complex III subunit, is also tightly connected with other MRCs in normal tissue (). Again, this connection disappeared under chemical treatment (). A similar loss of coexpression was observed for 23 MRC genes, although the remaining 23 genes (19 from complex I, 2 from complex III, and 2 from complex IV) were still highly coexpressed (Figure 2(b)). NDUFA1 is coexpressed with complex I subunit NDUFA4 and complex III subunit UQCRB. One possible explanation for these changes in coexpression in response to chemical treatment is that these compounds directly or indirectly influence MRC gene expression. Indeed, differential expression of mitochondrial genes has been induced by manipulating the agonal-pH state and through drug treatment [40, 41]. Furthermore, some compounds might modulate cellular redox levels or dissipate the mitochondrial membrane gradient by facilitating anion flux across the mitochondrial inner membrane, as suggested by Toogood [42], resulting in remodeling of the coexpression network.

3.3. Coexpression of MRC and Cellular Genes

In normal tissue, 1,422 genes are coexpressed significantly with MRC genes. Of these, coexpression of 1,308 genes was observed in normal tissue but not in chemically treated tissue. The remaining 114 genes are coexpressed in both normal and treated tissues. To identify the function of these 114 genes, we analyzed their associated GO terms using the online PANTHER tool (http://www.pantherdb.org/geneListAnalysis.do). A total of 51 (44.7%) genes were assigned to “metabolic process” (Figure 3(a)), indicating that this process is closely related to mitochondrial function. “Immune system process” and “apoptosis” were represented by 4.26% and 1.42% of these genes, respectively. The 1,308 genes that are coexpressed only in normal tissues were also categorized using PANTHER (Figure 3(b)). Of these, 26.56% were assigned to “metabolic processes” and 14.16% to “cellular processes.” “Immune system process” and “apoptosis” were assigned to 4.84% and 2.50% of these genes, respectively. AIFM1 (apoptosis-inducing factor, mitochondrion-associated 1) ENDOG (endonuclease G) is involved in intrinsic (mitochondria-associated) pathway for cancer cell apoptosis, for example, tightly coexpressed with MRC genes ( and 1.7139, resp.) in normal tissue, but this coexpression disappeared in treated tissues (Supplementary Table S1). AIFM1 is known to be important for the assembly and stability of complexes I and III [43]. In addition, the mutation or inhibition of MRC is widespread in cancer and intimately connected to apoptosis resistance [44], indicating that the MRC plays as a modulator of apoptosis for the treatment of cancer [45, 46]. Coexpression among MRC genes and the mechanistic target of rapamycin (mTOR) gene also disappeared after chemical treatment (; Supplementary Table S1). Disruption of the mTOR complex by treatment with the mTOR inhibitor rapamycin reduces mitochondrial membrane potential, oxygen consumption, and ATP synthetic capacity, indicating that formation of the mTOR complex is required for overall mitochondrial activity [47]. Taken together, these findings indicate that coexpression of MRC genes with cellular genes such as AIFM1 might be required to maintain the mitochondrial complexes. Furthermore, the disruption of this coexpression by chemical treatment suggests that similar disruptions might be responsible for mitochondria-related side effects of pharmaceuticals.

(a)

(b)

(c)

Figure 3: PANTHER analysis of the functional categories of genes coexpressed with MRCs. (a) Genes coexpressed in both normal tissues and chemically treated tissues. (b) Genes coexpressed with MRCs only in the tissue-based network. (c) Genes coexpressed with MRCs only in the chemically treated tissue.

In contrast to the coexpression network in normal tissues, only 238 genes were significantly enriched in chemically treated tissues. Of these, 114 are coexpressed in normal tissues, whereas 124 are coexpressed only in treated tissues. This finding suggests that chemical treatment alters the coexpression network between MRC genes and cellular genes. Of the genes coexpressed only in treated tissues, almost half (48.91%) were associated with the term “metabolic process” (Figure 3(c)). The gene for 2,4-dienoyl CoA reductase 2 (DECR2), an auxiliary enzyme in the mitochondrial beta-oxidation of unsaturated fatty acids, is coexpressed with a value of (Supplementary Table S1). Three acyl-CoA thioesterases (ACOTs), ACOT11 (), ACOT13 (), and ACOT2 (), involved in peroxisomal lipid metabolism [48], were highly coexpressed with MRCs in normal tissue, whereas ACOT8 () is coexpressed with MRCs in treated tissues (Supplementary Table S1).

Coexpression network, which is the reconstruction of biological networks from high-throughput data, can be used to identify higher-level features of gene-gene relationships based on graph theoretic considerations such as clustering coefficient or node degree [49, 50]. However, large-scale analyses only provide clues that help in forming a hypothesis [51]. Although the differences among coexpression networks (Supplementary Table S1) should help identify and prioritize candidate genes to determine the effects of drugs on mitochondria, further study is required to determine the relationship between the coexpressed genes and specific mitochondrial functions.

3.4. Functional Enrichment Analysis of MRC “Friends”

To investigate the biological processes represented by genes significantly coexpressed with MRC genes in normal tissue compared to treated tissues, we performed GO term enrichment analyses using the functional annotation tool DAVID (http://david.abcc.ncifcrf.gov/). For normal tissue, annotations for 1,308 genes were enriched in 17 terms, including “muscle system process,” “cellular metabolic process,” and “carboxylic acid metabolic process” (Table 1). For chemically treated tissues, GO term enrichment of 124 genes coexpressed with MRC genes found only ribosome biogenesis terms, such as “translational elongation and translation” (Table 2). A sufficient supply of ATP is required to maintain the contractile function of muscle [52], suggesting the importance of mitochondria during muscle contraction. Myosin provides energy and plays a vital role in muscle contraction. Myosin genes could be divided into several classifications [53], such as myosin heavy chain 1 (MYH1), myosin light chain 3 (MYL3), and myosin binding protein C2 (MYBPC2). Mutations in myosin genes lead to hypertrophic cardiomyopathy [54]. These findings indicate that muscle system closely interacts with MRC genes for improving mitochondrial function.

Table 1: GO enrichment of genes coexpressed with MRCs only in normal tissue transcriptome.

Term
Count
Pop
value
Genes

Muscle system process
68
168

MYBPC2, TNNC2, MYBPC1, TNNC1, GNA11, MYBPC3,   ANKRD2, PGAM2, KCNJ12, TTN, DES, DYSF, CHRNA1,   MAP2K6, MB, ACTA1, MYH1 , CRYAB, MYH2, TBCE, MYH4,  

MYLK2, ACTN2, PDE4D, MYH7, MYH6, CACNG1, MYH8,  

TNNT2, TRDN, TNNT3, TNNT1, PSEN2, RYR1, SMPX, RYR2,  

KBTBD10, STBD1, CASQ2, SGCA, CLCN1, ALDOA, MYL7,  

MYL4, TCAP, MYL3, MYL1, DAG1, MYOT, MYOM2, CKMT2,  

CAMK2D, MYOM1, SCN5A, HRC, ACTC1, MSTN, TRIM63,  

TNNI3, HOMER1, CACNA1S, TNNI2, SLC6A8, CHRNB4,  

GAMT, CHRNB1, SCN4A, SNTA1

Cellular ketone metabolic process
106
567

Carboxylic acid metabolic process
99
556

Oxoacid metabolic process
99
556

Glucose metabolic process
41
153

PRKAG3, ALDOA, LDHA, PHKA1, PRKAG2, PGAM2, OGDH, PDHB, ACN9, HIBADH, PPP1R3C, PPP1R1A, GYS1,   ENO3, PDHA1, GAPDH, AGL, PDK2, WDTC1, CRYAB,  

PHKG1, PDK4, EPM2A, BAD, DLAT, PPP1R3A, PFKM,  

FBP2, PPARGC1A, PPP1CB, GPI, PPP1R2, GBE1, PYGM,  

PGM1, GPT, PGK1, DCXR, UGP2, MDH2, MDH1

Hexose metabolic process
43
192

Acetyl-CoA metabolic process
17
31

Fatty acid metabolic process
39
198

PRKAG3, ACOX1, PPARA, TYRP1, ACADSB, CYP2J2,   PTGES2, ECH1, CPT2 , PRKAG1, PRKAG2, ACOT2,  

NDUFAB1, ECHS1, HADHA, HADHB, PEX7, ACSL1,  

ACOT11, ETFDH, GNPAT, PRKAA2, HADH, CPT1B, LPL,  

PLA2G15, ACADM, ACADS, MCAT, PRKAB2, CRAT,  

PPARGC1A , C9ORF3, ACADVL, UCP3, ANKRD23, MLYCD,  

FABP3, MECR

Table 2: GO enrichment of genes coexpressed with MRCs only in chemically treated tissue.

Term
Count
Pop hits
value
Genes

Translational elongation
12
101

RPS19, RPL32, RPL14, RPS29, RPL22, FAU, RPS10, RPL38, RPL12, RPS21, UBA52, RPL29

Translation
16
331

RPL14, MRPS21, RPL38, RPL29, MRPL11, RPS19, RPL32, RPS29, EIF3H, RPL22, FAU, RPS10, MRPL48, RPL12, RPS21, UBA52

Metabolic process
73
7647

Protein metabolic process
36
2812

Primary metabolic process
66
6923

Macromolecule metabolic process
56
5710

For the GO term “fatty acid metabolic process,” 39 genes were identified as “friends” of MRCs. Of these, carnitine palmitoyltransferase 2 (CPT2) localizes to the inner leaflet of the inner mitochondrial membrane, where it oxidizes long-chain fatty acids to produce substrates for the mitochondrial fatty acid beta-oxidation pathway [55]. Most of the genes required for mitochondrial biogenesis are controlled by DNA-binding transcription factors and coregulators [56]. Peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PPARGC1A) is a transcriptional coactivator that regulates various metabolic processes including mitochondrial biogenesis and respiration [57]. Therefore, coexpression of MRC genes with functional genes such as CPT2 and PPARGC1A might be required for the function and maintenance of mitochondria.

3.5. Coexpression of MRC Genes in Neurodegenerative Disease Pathways

Several molecular, cellular, biochemical, and animal model studies have suggested that mitochondrial dysfunction closely relates to the progression of several neurodegenerative diseases [58]. Using KEGG pathway enrichment analysis, we found that 19 MRC genes (41% of MRC genes in the coexpression network), belonging to complexes I to IV, are coregulated with Parkinson’s, Alzheimer’s, and Huntington’s disease pathways (Table 3), indicating the importance of MRCs in neurodegenerative disease pathways. In Parkinson’s disease, PTEN-induced putative kinase 1 (PINK1), a mitochondrial serine/threonine-protein kinase, was found to be a coexpression “friend” with MRC genes in normal tissue () but not in treated tissues () (Supplementary Table S1). PINK1 loss-of-function causes mitochondrial dysfunction and Parkinsonism [59]. In addition, the subunit of dihydropyridine receptor (CACNA1S), calmodulin-like 6 (CALML6), presenilin protein 2 (PSEN2), BCL2-associated agonist of cell death (BAD), and lipoprotein lipase (LPL) in Alzheimer’s disease pathway are significantly coexpressed with MRC genes (Table 3), whereas clathrin light chain B (CLTB), clathrin heavy polypeptide-like 1 (CLTCL1), DNA-directed RNA polymerase II subunits (POLR2E, POLR2L), and PPARGC1A were identified as coexpression “friends” with MRC genes in Huntington’s disease pathway. These findings indicate that MRC genes are directly or indirectly linked with neurodegenerative disease pathways. Mitochondrial dysfunction causes not only neurodegenerative diseases but also hypertrophic and dilated cardiomyopathy [60]. MRC “friend” genes were enriched in the hypertrophic and dilated cardiomyopathy pathways in normal tissue (Table 3), whereas coexpression again disappeared after treatment with chemicals. KEGG pathway enrichment analysis thus suggests that the coexpression network revealed by using MRC genes as seed genes provides a possible link between mitochondria and various disease pathways.

Table 3: KEGG enrichment analysis of genes coexpressed with MRCs.

Term
Count
Pop
value
Genes

Parkinson’s disease
40
128

UQCRC2, ATP5D, NDUFB4, NDUFB6, NDUFB7, UBE2G1, CYC1, NDUFAB1, PINK1 , ATP5G2, UQCRFS1, COX5A, NDUFB1, NDUFB2, UQCR11, ATP5O, ATP5H, NDUFS1, ATP6, NDUFA5, NDUFA2, COX7A1, SLC25A4, NDUFA6, CYCS, NDUFC2, ATP5F1, COX4I1, NDUFC1, NDUFA10, VDAC2, VDAC3, VDAC1, SDHA, NDUFV3, PPID, NDUFV2, SDHD, COX6A2, ATP5A1

Hypertrophic cardiomyopathy (HCM)
32
85

PRKAG3, MYL2, TNNC1, MYL3, PRKAG1, MYBPC3, PRKAG2, CACNB1, DAG1, TPM2, TTN, TPM3, DES, DMD, ITGB6, PRKAA2, ACTC1, CACNG6, PRKAB2, MYH7, MYH6, CACNG1, TNNI3, CACNA1S , TNNT2, ATP2A2, SGCG, ITGA7, SGCD, RYR2, SGCA, SGCB

Alzheimer’s disease
41
163

UQCRC2, ATP5D, NDUFB4, NDUFB6, NDUFB7, CYC1, NDUFAB1, ATP5G2, UQCRFS1, COX5A, NDUFB1, NDUFB2, UQCR11, CALML6 , ATP5O, ATP5H, GAPDH, NDUFS1, ATP6, NDUFA5, LPL , NDUFA2, COX7A1, NDUFA6, CYCS, NDUFC2, ATP5F1, COX4I1, BAD , NDUFC1, NDUFA10, CACNA1S, SDHA, NDUFV3, ATP2A2, ATP2A1, NDUFV2, SDHD, PSEN2 , COX6A2, ATP5A1

Dilated cardiomyopathy
30
92

MYL2, TNNC1, MYL3, MYBPC3, CACNB1, DAG1, TTN, TPM2, TPM3, DES, DMD, ITGB6, PRKACA, ACTC1, CACNG6, MYH7, MYH6, CACNG1, TNNI3, CACNA1S, TNNT2, ADCY9, ATP2A2, SGCG, PLN, ITGA7, SGCD, RYR2, SGCA, SGCB

Huntington’s disease
43
180

UQCRC2, ATP5D, NDUFB4, POLR2E , NDUFB6, CLTB , NDUFB7, POLR2LL , CYC1, NDUFAB1, ATP5G2, UQCRFS1, COX5A, NDUFB1, NDUFB2, UQCR11, ATP5O, ATP5H, NDUFS1, ATP6, NDUFA5, NDUFA2, COX7A1, SLC25A4, NDUFA6, CYCS, ATP5F1, NDUFC2, COX4I1, NDUFC1, NDUFA10, VDAC2, VDAC3, PPARGC1A , VDAC1, SDHA, NDUFV3, PPID, NDUFV2, SDHD, COX6A2, ATP5A1, CLTCL1

4. Conclusion

In this study, we have compared the networks of genes coexpressed with MRC genes in normal and chemically treated tissues. We find a differential distribution of coexpression after chemical treatment. These differences might be mediated by chemical-related stimuli, suggesting that coexpression network analysis can provide helpful information for understanding side effects of drugs on mitochondrial functions.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

Cong Chen and Tae Kyung Hyun contributed equally to this paper.

Acknowledgments

This original study was partly supported by the National Natural Science Foundation of China (Key Program, no. 31200620), China Postdoctoral Science Foundation (CPSF; no. 2013M540739), Shaanxi Postdoctoral Science Foundation, and Xi’an Jiaotong University (985 and 211 Projects).

References

M. P. Yaffe, “The machinery of mitochondrial inheritance and behavior,” Science, vol. 283, no. 5407, pp. 1493–1497, 1999.

D. D. Newmeyer and S. Ferguson-Miller, “Mitochondria: releasing power for life and unleashing the machineries of death,” Cell, vol. 112, no. 4, pp. 481–490, 2003.

B. Muster, W. Kohl, I. Wittig et al., “Respiratory chain complexes in dynamic mitochondria display a patchy distribution in life cells,” PLoS ONE, vol. 5, no. 7, Article ID e11910, 2010.

J. S. Carew and P. Huang, “Mitochondrial defects in cancer,” Molecular Cancer, vol. 1, article 9, 2002.

H. Lu, V. Koshkin, E. M. Allister, A. V. Gyulkhandanyan, and M. B. Wheeler, “Molecular and metabolic evidence for mitochondrial defects associated with β-cell dysfunction in a mouse model of type 2 diabetes,” Diabetes, vol. 59, no. 2, pp. 448–459, 2010.

S. Ohsawa, Y. Sato, M. Enomoto, M. Nakamura, A. Betsumiya, and T. Igaki, “Mitochondrial defect drives non-autonomous tumour progression through Hippo signalling in Drosophila,” Nature, vol. 490, no. 7421, pp. 547–551, 2012.

J. Labbadia and R. I. Morimoto, “Huntington's disease: underlying molecular mechanisms and emerging concepts,” Trends in Biochemical Sciences, vol. 38, no. 8, pp. 378–385, 2013.

D. J. M. Fernández-Ayala, S. Chen, E. Kemppainen, K. M. C. O'Dell, and H. T. Jacobs, “Gene expression in a Drosophila model of mitochondrial disease,” PLoS ONE, vol. 5, no. 1, Article ID e8549, 2010.

V. S. Dhillon and M. Fenech, “Mutations that affect mitochondrial functions and their association with neurodegenerative diseases,” Mutation Research, vol. 759, pp. 1–13, 2014.

S. Movahedi, M. van Bel, K. S. Heyndrickx, and K. Vandepoele, “Comparative co-expression analysis in plant biology,” Plant, Cell and Environment, vol. 35, no. 10, pp. 1787–1798, 2012.

X. Han, C. Chen, T. K. Hyun, R. Kumar, and J. Y. Kim, “Metabolic module mining based on independent component analysis in Arabidopsis thaliana,” Molecules and Cells, vol. 34, no. 3, pp. 295–304, 2012.

N. J. Hudson, A. Reverter, Y. H. Wang, P. L. Greenwood, and B. P. Dalrymple, “Inferring the transcriptional landscape of bovine skeletal muscle by integrating co-expression networks,” PLoS ONE, vol. 4, no. 10, Article ID e7249, 2009.

S. van Dam, R. Cordeiro, T. Craig, J. van Dam, S. H. Wood, and J. P. de Magalhães, “GeneFriends: an online co-expression analysis tool to identify novel gene targets for aging and complex diseases,” BMC Genomics, vol. 13, no. 1, article 535, 2012.

F. Azuaje, L. Zhang, C. Jeanty, S. Puhl, S. Rodius, and D. R. Wagner, “Analysis of a gene co-expression network establishes robust association between Col5a2 and ischemic heart disease,” BMC Medical Genomics, vol. 6, no. 1, article 13, 2013.

A. Mourier and N. G. Larsson, “Tracing the trail of protons through complex I of the mitochondrial respiratory chain,” PLoS Biology, vol. 9, no. 8, Article ID e1001129, 2011.

M. P. Murphy, A. Holmgren, N. Larsson et al., “Unraveling the biological roles of reactive oxygen species,” Cell Metabolism, vol. 13, no. 4, pp. 361–366, 2011.

H. Pagniez-Mammeri, S. Loublier, A. Legrand, P. Bénit, P. Rustin, and A. Slama, “Mitochondrial complex I deficiency of nuclear origin. I. Structural genes,” Molecular Genetics and Metabolism, vol. 105, no. 2, pp. 163–172, 2012.

E. Fassone and S. Rahman, “Complex I deficiency: clinical features, biochemistry and molecular genetics,” Journal of Medical Genetics, vol. 49, no. 9, pp. 578–590, 2012.

A. Benchoua, Y. Trioulier, D. Zala et al., “Involvement of mitochondrial complex II defects in neuronal death produced by N-terminus fragment of mutated huntingtin,” Molecular Biology of the Cell, vol. 17, no. 4, pp. 1652–1663, 2006.

P. Gaignard, M. Menezes, M. Schiff et al., “Mutations in CYC1, encoding cytochrome c1 subunit of respiratory chain complex III, cause insulin-responsive hyperglycemia,” The American Journal of Human Genetics, vol. 93, no. 2, pp. 384–389, 2013.

S. Jain-Ghai, J. M. Cameron, A. Al Maawali et al., “Complex II deficiency—a case report and review of the literature,” American Journal of Medical Genetics A, vol. 161, no. 2, pp. 285–294, 2013.

D. C. Wallace, “Mitochondrial diseases in man and mouse,” Science, vol. 283, no. 5407, pp. 1482–1488, 1999.

E. J. Lesnefsky and C. L. Hoppel, “Oxidative phosphorylation and aging,” Ageing Research Reviews, vol. 5, no. 4, pp. 402–433, 2006.

F. Scaglia, J. A. Towbin, W. J. Craigen et al., “Clinical spectrum, morbidity, and mortality in 113 pediatric patients with mitochondrial disease,” Pediatrics, vol. 114, no. 4, pp. 925–931, 2004.

H. K. Lee, A. K. Hsu, J. Sajdak, J. Qin, and P. Pavlidis, “Coexpresion analysis of human genes across many microarray data sets,” Genome Research, vol. 14, no. 6, pp. 1085–1094, 2004.

C. M. Sue and E. A. Schon, “Mitochondrial respiratory chain diseases and mutations in nuclear DNA: a promising start?” Brain Pathology, vol. 10, no. 3, pp. 442–450, 2000.

A. F. Santidrian, A. Matsuno-Yagi, M. Ritland et al., “Mitochondrial complex I activity and NAD+/NADH balance regulate breast cancer progression,” The Journal of Clinical Investigation, vol. 123, no. 3, pp. 1068–1081, 2013.

R. Acín-Pérez, P. Fernández-Silva, M. L. Peleato, A. Pérez-Martos, and J. A. Enriquez, “Respiratory active mitochondrial supercomplexes,” Molecular Cell, vol. 32, no. 4, pp. 529–539, 2008.

E. Lapuente-Brun, R. Moreno-Loshuertos, R. Aciń-Pérez et al., “Supercomplex assembly determines electron flux in the mitochondrial electron transport chain,” Science, vol. 340, no. 6140, pp. 1567–1570, 2013.

C. Blanchi, M. L. Genova, G. P. Castelli, and G. Lenaz, “The mitochondrial respiratory chain is partially organized in a supercomplex assembly: kinetic evidence using flux control analysis,” The Journal of Biological Chemistry, vol. 279, no. 35, pp. 36562–36569, 2004.

N. V. Dudkina, M. Kudryashev, H. Stahlberg, and E. J. Boekema, “Interaction of complexes I, III, and IV within the bovine respirasome by single particle cryoelectron tomography,” Proceedings of the National Academy of Sciences of the United States of America, vol. 108, no. 37, pp. 15196–15200, 2011.

R. Acín-Pérez, M. P. Bayona-Bafaluy, P. Fernández-Silva et al., “Respiratory complex III is required to maintain complex I in mammalian mitochondria,” Molecular Cell, vol. 13, no. 6, pp. 805–815, 2004.

C. Ugalde, R. J. R. J. Janssen, L. P. van den Heuvel, J. A. M. Smeitink, and L. G. J. Nijtmans, “Differences in assembly or stability of complex I and other mitochondrial OXPHOS complexes in inherited complex I deficiency,” Human Molecular Genetics, vol. 13, no. 6, pp. 659–667, 2004.

E. L. Blakely, A. L. Mitchell, N. Fisher et al., “A mitochondrial cytochrome b mutation causing severe respiratory chain enzyme deficiency in humans and yeast,” FEBS Journal, vol. 272, no. 14, pp. 3583–3592, 2005.

N. V. Dudkina, R. Kouřil, K. Peters, H. Braun, and E. J. Boekema, “Structure and function of mitochondrial supercomplexes,” Biochimica et Biophysica Acta, vol. 1797, no. 6-7, pp. 664–670, 2010.

J. Neustadt and S. R. Pieczenik, “Medication-induced mitochondrial damage and disease,” Molecular Nutrition and Food Research, vol. 52, no. 7, pp. 780–788, 2008.

M. V. Tsiper, J. Sturgis, L. V. Avramova et al., “Differential mitochondrial toxicity screening and multi-parametric data analysis,” PLoS ONE, vol. 7, no. 10, Article ID e45226, 2012.

P. V. Gopalacharyulu, V. R. Velagapudi, E. Lindfors, E. Halperin, and M. Orešič, “Dynamic network topology changes in functional modules predict responses to oxidative stress in yeast,” Molecular BioSystems, vol. 5, no. 3, pp. 276–287, 2009.

S. Lehtinen, F. X. Marsellach, S. Codlin et al., “Stress induces remodelling of yeast interaction and co-expression networks,” Molecular BioSystems, vol. 9, no. 7, pp. 1697–1707, 2013.

M. P. Vawter, H. Tomita, F. Meng et al., “Mitochondrial-related gene expression changes are sensitive to agonal-pH state: implications for brain disorders,” Molecular Psychiatry, vol. 11, no. 7, pp. 663–679, 2006.

T. Lee, M. G. Manjanatha, A. Aidoo et al., “Expression analysis of hepatic mitochondria-related genes in mice exposed to acrylamide and glycidamide,” Journal of Toxicology and Environmental Health A, vol. 75, no. 6, pp. 324–339, 2012.

P. L. Toogood, “Mitochondrial drugs,” Current Opinion in Chemical Biology, vol. 12, no. 4, pp. 457–463, 2008.

L. Delavallée, L. Cabon, P. Galán-Malo, H. K. Lorenzo, and S. A. Susin, “AIF-mediated caspase-independent necroptosis: a new chance for targeted therapeutics,” IUBMB Life, vol. 63, no. 4, pp. 221–232, 2011.

A. Lemarie and S. Grimm, “Mitochondrial respiratory chain complexes: apoptosis sensors mutated in cancer?” Oncogene, vol. 30, no. 38, pp. 3985–4003, 2011.

J. Q. Kwong, M. S. Henning, A. A. Starkov, and G. Manfredi, “The mitochondrial respiratory chain is a modulator of apoptosis,” The Journal of Cell Biology, vol. 179, no. 6, pp. 1163–1177, 2007.

S. Fulda, “Exploiting mitochondrial apoptosis for the treatment of cancer,” Mitochondrion, vol. 10, no. 6, pp. 598–603, 2010.

S. M. Schieke, D. Phillips, J. P. McCoy Jr. et al., “The mammalian target of rapamycin (mTOR) pathway regulates mitochondrial oxygen consumption and oxidative capacity,” The Journal of Biological Chemistry, vol. 281, no. 37, pp. 27643–27652, 2006.

M. C. Hunt, K. Solaas, B. Frode Kase, and S. E. H. Alexson, “Characterization of an acyl-CoA thioesterase that functions as a major regulator of peroxisomal lipid metabolism,” Journal of Biological Chemistry, vol. 277, no. 2, pp. 1128–1138, 2002.

M. Mistry, J. Gillis, and P. Pavlidis, “Meta-analysis of gene coexpression networks in the post-mortem prefrontal cortex of patients with schizophrenia and unaffected controls,” BMC Neuroscience, vol. 14, article 105, 2013.

R. de Smet and K. Marchal, “Advantages and limitations of current network inference methods,” Nature Reviews Microbiology, vol. 8, no. 10, pp. 717–729, 2010.

M. Y. Hirai, “A robust omics-based approach for the identification of glucosinolate biosynthetic genes,” Phytochemistry Reviews, vol. 8, no. 1, pp. 15–23, 2009.

R. K. Dash, B. M. Bell, M. J. Kushmerick, and P. Vicini, “Estimating in vitro mitochondrial oxygen consumption during muscle contraction and recovery: a novel approach that accounts for diffusion,” Annals of Biomedical Engineering, vol. 33, no. 3, pp. 343–355, 2005.

J. R. Sellers, “Myosins: a diverse superfamily,” Biochimica et Biophysica Acta, vol. 1496, no. 1, pp. 3–22, 2000.

E. M. McNally, J. R. Golbus, and M. J. Puckelwartz, “Genetic mutations and mechanisms in dilated cardiomyopathy,” Journal of Clinical Investigation, vol. 123, no. 1, pp. 19–26, 2013.

S. Violante, L. IJlst, H. Te Brinke et al., “Carnitine palmitoyltransferase 2 and carnitine/acylcarnitine translocase are involved in the mitochondrial synthesis and export of acylcarnitines,” FASEB Journal, vol. 27, no. 5, pp. 2039–2044, 2013.

C. A. Piantadosi and H. B. Suliman, “Transcriptional control of mitochondrial biogenesis and its interface with inflammatory processes,” Biochimica et Biophysica Acta, vol. 1820, no. 4, pp. 532–541, 2012.

J. St-Pierre, S. Drori, M. Uldry et al., “Suppression of reactive oxygen species and neurodegeneration by the PGC-1 transcriptional coactivators,” Cell, vol. 127, no. 2, pp. 397–408, 2006.

P. H. Reddy, “Role of mitochondria in neurodegenerative diseases: mitochondria as a therapeutic target in Alzheimer's disease,” CNS Spectrums, vol. 14, pp. 8–13, 16–18, 2009.

D. P. Narendra, S. M. Jin, A. Tanaka et al., “PINK1 is selectively stabilized on impaired mitochondria to activate Parkin,” PLoS Biology, vol. 8, no. 1, Article ID e1000298, 2010.

D. E. Meyers, H. I. Basha, and M. K. Koenig, “Mitochondrial cardiomyopathy: pathophysiology, diagnosis, and management,” Texas Heart Institute Journal, vol. 40, pp. 385–394, 2013.