mtDNA eQTLs and the m1A 16S rRNA modification explain mtDNA tissue-specific gene expression pattern in humans

Expression quantitative trait loci (eQTLs) are instrumental in genome-wide identification of regulatory elements, yet were overlooked in the mitochondrial DNA (mtDNA). By analyzing 5079 RNA-seq samples from 23 tissues we identified association of ancient mtDNA SNPs (haplogroups T2, L2, J2 and V) and recurrent SNPs (mtDNA positions 263, 750, 1438 and 10398) with tissue-dependent mtDNA gene-expression. Since the recurrent SNPs independently occurred in different mtDNA genetic backgrounds, they constitute the best candidates to be causal eQTLs. Secondly, the discovery of mtDNA eQTLs in both coding and non-coding mtDNA regions, propose the identification of novel mtDNA regulatory elements. Third, we identified association between low m1A 947 MT-RNR2 (16S) rRNA modification levels and altered mtDNA gene-expression in twelve tissues. Such association disappeared in skin which was exposed to sun, as compared to sun-unexposed skin from the same individuals, thus supporting the impact of UV on mtDNA gene expression. Taken together, our findings reveal that both mtDNA SNPs and mt-rRNA modification affect mtDNA gene expression in a tissue-dependent manner.