RPKM Data From GTEx
TPM Data
Normalization for Age, Gender and Cause of Death via linear regression
Expression Analysis Using Residuals as Expression Values
Across-Tissues Correlation Within-Tissues Analysis Using Random Sampling Correlation Analysis
Supplementary Figure 1. Workflow of the expression pattern correlation analysis. See also dataset S3 for Python scripts used during the analysis. Supplementary Figure 2. Tissue-dependent OXPHOS gene expression. A clustered binomial heat map demonstrating similarity in the expression pattern of OXPHOS genes which are not included in the 'main OXPHOS cluster' to the expression of genes within this cluster in specific tissues (red). Dissimilarities of OXPHOS genes from the 'main OXPHOS cluster' are indicated in blue. The orange and green rectangles show the two genes groups, (COX7A1 and COX6A2) and (COX7A2, NDUFS4, NDUFC2, NDUFA11 and ATPIF1), respectively. The gray rectangles show brain tissues.
. *** Spearman's Rank Correlation Coefficient Correlation Rank Spearman's
intra-OXPHOS OXPHOS-Genome
Supplementary Figure 3. nDNA-encoded OXPHOS genes show significantly higher co-expression with mtDNA-encoded OXPHOS genes. A box-plot showing the values of the correlation coefficients between pairs of OXPHOS genes (intra-OXPHOS) and between OXPHOS genes and randomly chosen protein- coding genes (OXPHOS-genome). *** p <1E-100. Supplementary Figure 4. The correlation distribution of nDNA- and mtDNA-encoded OXPHOS genes is bi-modal. A histogram comparing the expression distribution of nDNA- encoded OXPHOS genes (green) and all other genes in the genome (blue) with medial mtDNA expression in all tissues. Median TPM TPM values of the two LDH subunit-coding genes (LDHA and LDHB) across all human tissues inspected in this in inspected tissues human all study. across LDHB) and (LDHA genes subunit-coding LDH median two the the showing of bar-plot values A TPM tissues. human across subunits LDH of Expression 5. Figure Supplementary (A)
Transcription Factor Names CCNT2 CEBPB CHD2 CREB1 CTCF E2F1 E2F4 E2F6 ELF1 ELK1 FOXP2 GTF2F1 HDAC1 IRF1 MAFK MAX MAZ MXI1 MYC NFYA NFYB NRF1 PHF8 POLR2A RBBP5 REST RUNX3 SIN3A SIN3AK20 SMC3 SP1 SREBP1 TAF1 TBP ZKSCAN1 DGF ATF3
Bases upstreamto thetranslationstartsite(TSS) −1000
−900
−800
−700
−600
−500
−400
−300
−200
−100
0 ATG (B)
Gene Names NDUFA8 NDUFA9 NDUFAB1 NDUFB1 NDUFB10 NDUFB11 NDUFB2 NDUFB3 NDUFB4 NDUFB5 NDUFB6 NDUFB7 NDUFB8 NDUFB9 NDUFC1 NDUFC2 NDUFS1 NDUFS2 NDUFS3 NDUFS4 NDUFS5 NDUFS6 NDUFS7 NDUFS8 NDUFV1 NDUFV2 NDUFV3 SDHA SDHB SDHC SDHD UQCR10 UQCR11 UQCRB UQCRC1 UQCRFS1 UQCRH UQCRQ ATP5A1 ATP5B ATP5C1 ATP5D ATP5E ATP5F1 ATP5G1 ATP5G2 ATP5G3 ATP5H ATP5I ATP5J ATP5J2 ATP5L ATP5L2 ATP5O ATPIF1 COX4I1 COX4I2 COX5A COX5B COX6A1 COX6A2 COX6B1 COX6B2 COX6C COX7A1 COX7A2 COX7B COX7B2 COX7C COX8A COX8C CYC1 NDUFA1 NDUFA10 NDUFA11 NDUFA12 NDUFA13 NDUFA2 NDUFA3 NDUFA4 NDUFA5 NDUFA6 NDUFA7
Bases upstreamto thetranslationstartsite(TSS) −1000
−750
−500
−250
0 ATG Supplementary Figure 6. Graphical representation of the in vivo identification of candidate promoters in nuclear DNA-encoded OXPHOS genes. (A) A map of the upstream region of a representative OXPHOS subunit (NDUFA8) tested for TF binding in DNaseI Genomic Footprinting (DGF) sites. X-axis represent the positions in the promoter region, relative to the translation start site (TSS). DGF sites are indicated by blue lines in the upper row and below, all the listed TFs binding to the 1000 bp promoter region are shown, with their ChIP-seq peaks shown as black lines. Note that only TFs with a ChIP-seq peak overlapping a minimum of one DGF are indicated. (B) Representation of the results for all nDNA-encoded OXPHOS genes. X-axis represent the positions in the promoter region, relative to the translation start site (TSS). DGFs are shown by blue lines if they do not overlap with any identified TF-binding peak in the ChIP-seq experiments and in black and blue lines if they overlap with at least one such TF-binding peak. Supplementary Table 1. A list of OXPHOS genes analyzed in this study, following HUGO Gene Nomenclature Committee (HGNC). OXPHOS structural and assembly factors are indicated, along with their Ensemble accession numbers.
Structural genes Structural genes Structural genes Structural genes Assembly Factors Name Ensemble ID Name Ensemble ID Name Ensemble ID Name Ensemble ID Name Ensemble ID MT-ND1 ENSG00000198888 NDUFB2 ENSG00000090266 UQCR10 ENSG00000184076 ATP5F1 ENSG00000116459 COX15 ENSG00000014919 MT-ND2 ENSG00000198763 NDUFB3 ENSG00000119013 UQCR11 ENSG00000127540 ATP5G1 ENSG00000159199 COX17 ENSG00000138495 MT-ND3 ENSG00000198840 NDUFB4 ENSG00000065518 NDUFA4 ENSG00000189043 ATP5G2 ENSG00000135390 COX18 ENSG00000163626 MT-ND4 ENSG00000198886 NDUFB5 ENSG00000136521 COX4I1 ENSG00000131143 ATP5G3 ENSG00000154518 COX19 ENSG00000240230 MT-ND4L ENSG00000212907 NDUFB6 ENSG00000165264 COX4I2 ENSG00000131055 ATP5H ENSG00000167863 COX20 ENSG00000203667 MT-ND5 ENSG00000198786 NDUFB7 ENSG00000099795 COX5A ENSG00000178741 ATP5I ENSG00000169020 ECSIT ENSG00000130159 MT-ND6 ENSG00000198695 NDUFB8 ENSG00000166136 COX5B ENSG00000135940 ATP5J ENSG00000154723 FOXRED1 ENSG00000110074 NDUFS1 ENSG00000023228 NDUFB9 ENSG00000147684 COX6A1 ENSG00000111775 ATP5J2 ENSG0000024146 NDUFAF1 ENSG00000137806 NDUFS2 ENSG00000158864 NDUFB10 ENSG00000140990 COX6A2 ENSG00000156885 ATP5L ENSG00000167283 NDUFAF2 ENSG00000164182 NDUFS3 ENSG00000213619 NDUFB11 ENSG00000147123 COX6B1 ENSG00000126267 ATP5L2 ENSG00000249222 NDUFAF3 ENSG00000178057 NDUFS7 ENSG00000115286 NDUFC1 ENSG00000109390 COX6B2 ENSG00000160471 ATP5O ENSG00000241837 NDUFAF4 ENSG00000123545 NDUFS8 ENSG00000110717 NDUFC2 ENSG00000151366 COX6C ENSG00000164919 MT-ATP6 ENSG00000198899 NDUFAF5 ENSG00000101247 NDUFV1 ENSG00000167792 NDUFS4 ENSG00000164258 COX7A1 ENSG00000161281 MT-ATP8 ENSG00000228253 NDUFAF6 ENSG00000156170 NDUFV2 ENSG00000178127 NDUFS5 ENSG00000168653 COX7A2 ENSG00000112695 Assembly Factors NDUFAF7 ENSG00000003509 NDUFAB1 ENSG00000004779 NDUFS6 ENSG00000145494 COX7B ENSG00000131174 Name Ensemble ID NUBPL ENSG00000151413 NDUFA1 ENSG00000125356 NDUFV3 ENSG00000160194 COX7B2 ENSG00000170516 ACAD9 ENSG00000177646 SCO1 ENSG00000133028 NDUFA2 ENSG00000131495 SDHA ENSG00000073578 COX7C ENSG00000127184 ATPAF1 ENSG00000123472 SCO2 ENSG00000130489 NDUFA3 ENSG00000170906 SDHB ENSG00000117118 COX8A ENSG00000176340 ATPAF2 ENSG00000171953 SDHAF1 ENSG00000205138 NDUFA5 ENSG00000128609 SDHC ENSG00000143252 COX8C ENSG00000187581 BCS1L ENSG00000074582 SDHAF2 ENSG00000167985 NDUFA6 ENSG00000184983 SDHD ENSG00000204370 MT-CO1 ENSG00000198804 COA1 ENSG00000106603 SDHAF3 ENSG00000196636 NDUFA7 ENSG00000267855 CYC1 ENSG00000179091 MT-CO2 ENSG00000198712 COA3 ENSG00000183978 SDHAF4 ENSG00000154079 NDUFA8 ENSG00000119421 MT-CYB ENSG00000198727 MT-CO3 ENSG00000198938 COA4 ENSG00000181924 SURF1 ENSG00000148290 NDUFA9 ENSG00000139180 UQCRB ENSG00000156467 ATPIF1 ENSG00000130770 COA5 ENSG00000183513 TIMMDC1 ENSG00000113845 NDUFA10 ENSG00000130414 UQCRC1 ENSG00000010256 ATP5A1 ENSG00000152234 COA6 ENSG00000168275 TMEM126B ENSG00000171204 NDUFA11 ENSG00000174886 UQCRC2 ENSG00000140740 ATP5B ENSG00000110955 COA7 ENSG00000162377 UQCC1 ENSG00000101019 NDUFA12 ENSG00000184752 UQCRFS1 ENSG00000169021 ATP5C1 ENSG00000165629 COX10 ENSG00000006695 UQCC2 ENSG00000137288 NDUFA13 ENSG00000186010 UQCRH ENSG00000173660 ATP5D ENSG00000099624 COX11 ENSG00000166260 UQCC3 ENSG00000204922 NDUFB1 ENSG00000183648 UQCRQ ENSG00000164405 ATP5E ENSG00000124172 COX14 ENSG00000178449
Supplementary Table 3. Summary statistics of the expression differences between subunits within or between OXPHOS complexes.
Complex MWU p-value Median Median
intra-complex inter-complexes
CI 1024010 0.00141 0.6171 0.5988
CII 2454 0.00687 0.56 0.4773
CIII 23912 0.07942 0.655 0.635
CIV 38126 6.21E-06 0.6957 0.6379
CV 125622 0.01499 0.6348 0.6193
Total 3724116 9.96E-07 0.6249 0.6048
Supplementary Table 4. A mathematical representation of the effect of negative expression correlation between mtDNA- and nDNA-encoded OXPHOS genes on the
OXPHOS complex assembly range. M0, N0 = Initial expression values of mtDNA (M0) and nDNA (N0) genes, ΔM, ΔN = Changes in expression of mtDNA (ΔM) and nDNA (ΔN) genes.
Correlation Initial ΔM ΔN Scenario Complexes
Direction State Complexes Range
Range (SR)
(SCR)
M0>N0 + + N0≤SCR≤∞ + 0≤SR≤∞ M0>N0 - - 0≤SCR≤N0
M0>N0 + - 0≤SCR≤N0 - 0≤SR Supplementary Table 7. Statistically significant Spearman's rank correlation coefficient positively affecting TFs. * Candidate promoter regions. Transcription Factor nDNA OXPHOS genes' Spearman's p-value promoter coverage* (%) correlation-coefficient response (SCU) POLR2A 92.59 0.574848 <1E-6 PHF8 86.42 0.173793 <1E-6 MAZ 81.48 0.055696 0.000002 JUND 72.84 0.042717 0.000273 KDM5B 69.14 0.043484 <1E-6 CREB1 65.43 0.042579 <1E-6 TBL1XR1 55.56 0.032961 0.000203 CCNT2 53.09 0.035662 0.000034 CEBPB 50.62 0.036093 <1E-6 TCF12 49.38 0.031816 0.000188 BRCA1 40.74 0.033497 0.000016 ZBTB33 39.51 0.053593 <1E-6 FOS 35.80 0.036259 0.000022 BCLAF1 29.63 0.038006 0.000270 FOSL2 24.69 0.047823 0.000065