Mutational profiling of mtDNA control region reveals tumor-specific evolutionary selection involved in mitochondrial dysfunction.

backgroundMitochondrial DNA (mtDNA) mutations alter mitochondrial function in oxidative metabolism and play an important role in tumorigenesis. A series of studies have demonstrated that the mtDNA control region (mtCTR), which is essential for mtDNA replication and transcription, represents a mutational hotspot in human tumors. However, a comprehensive pan-cancer evolutionary pattern analysis of mtCTR mutations is urgently needed.

methodsWe generated a comprehensive combined dataset containing 10026 mtDNA somatic mutations from 4664 patients, covering 20 tumor types based on public and private next-generation sequencing data.

findingsOur results demonstrated a significantly higher and much more variable mutation rate in mtCTR than in the coding region across different tumor types. Moreover, our data showed a remarkable distributional bias of tumor somatic mutations between the hypervariable segment (HVS) and non-HVS, with a significantly higher mutation density and average mutation sites in HVS. Importantly, the tumor-specific mutational pattern between mtCTR HVS and non-HVS was identified, which was classified into three evolutionary selection types (relaxed, moderate, and strict constraint types). Analysis of substitution patterns revealed that the prevalence of C

interpretationOur results suggest that somatic mutations in mtCTR may be shaped by tumor-specific selective pressure and are involved in tumorigenesis. FUNDINGS: National Natural Science Foundation of China [grants 82020108023, 81830070, 81872302], and Autonomous Project of State Key Laboratory of Cancer Biology, China [grants CBSKL2019ZZ06, CBSKL2019ZZ27].