Secondary structure of the human mitochondrial genome affects formation of deletions

BACKGROUND

Aging in postmitotic tissues is associated with clonal expansion of somatic mitochondrial deletions, the origin of which is not well understood. Such deletions are often flanked by direct nucleotide repeats, but this alone does not fully explain their distribution. Here, we hypothesized that the close proximity of direct repeats on single-stranded mitochondrial DNA (mtDNA) might play a role in the formation of deletions.

RESULTS

By analyzing human mtDNA deletions in the major arc of mtDNA, which is single-stranded during replication and is characterized by a high number of deletions, we found a non-uniform distribution with a "hot spot" where one deletion breakpoint occurred within the region of 6-9 kb and another within 13-16 kb of the mtDNA. This distribution was not explained by the presence of direct repeats, suggesting that other factors, such as the spatial proximity of these two regions, can be the cause. In silico analyses revealed that the single-stranded major arc may be organized as a large-scale hairpin-like loop with a center close to 11 kb and contacting regions between 6-9 kb and 13-16 kb, which would explain the high deletion activity in this contact zone. The direct repeats located within the contact zone, such as the well-known common repeat with a first arm at 8470-8482 bp (base pair) and a second arm at 13,447-13,459 bp, are three times more likely to cause deletions compared to direct repeats located outside of the contact zone. A comparison of age- and disease-associated deletions demonstrated that the contact zone plays a crucial role in explaining the age-associated deletions, emphasizing its importance in the rate of healthy aging.

CONCLUSIONS

Overall, we provide topological insights into the mechanism of age-associated deletion formation in human mtDNA, which could be used to predict somatic deletion burden and maximum lifespan in different human haplogroups and mammalian species.

Genomic mutation landscape of skin cancers from DNA repair-deficient xeroderma pigmentosum patients

Xeroderma pigmentosum (XP) is a genetic disorder caused by mutations in genes of the Nucleotide Excision Repair (NER) pathway (groups A-G) or in Translesion Synthesis DNA polymerase η (V). XP is associated with an increased skin cancer risk, reaching, for some groups, several thousand-fold compared to the general population. Here, we analyze 38 skin cancer genomes from five XP groups. We find that the activity of NER determines heterogeneity of the mutation rates across skin cancer genomes and that transcription-coupled NER extends beyond the gene boundaries reducing the intergenic mutation rate. Mutational profile in XP-V tumors and experiments with POLH knockout cell line reveal the role of polymerase η in the error-free bypass of (i) rare TpG and TpA DNA lesions, (ii) 3' nucleotides in pyrimidine dimers, and (iii) TpT photodimers. Our study unravels the genetic basis of skin cancer risk in XP and provides insights into the mechanisms reducing UV-induced mutagenesis in the general population.

Integrative pan-cancer genomic and transcriptomic analyses of refractory metastatic cancer

Metastatic relapse after treatment is the leading cause of cancer mortality, and known resistance mechanisms are missing for most treatments administered to patients. To bridge this gap, we analyze a pan-cancer cohort (META-PRISM) of 1,031 refractory metastatic tumors profiled via whole-exome and transcriptome sequencing. META-PRISM tumors, particularly prostate, bladder, and pancreatic types, displayed the most transformed genomes compared to primary untreated tumors. Standard-of-care resistance biomarkers were identified only in lung and colon cancers - 9.6% of META-PRISM tumors, indicating that too few resistance mechanisms have received clinical validation. In contrast, we verified the enrichment of multiple investigational and hypothetical resistance mechanisms in treated compared to non-treated patients, thereby confirming their putative role in treatment resistance. Additionally, we demonstrated that molecular markers improve six-month survival prediction, particularly in patients with advanced breast cancer. Our analysis establishes the utility of META-PRISM cohort for investigating resistance mechanisms and performing predictive analyses in cancer.

Abstract PR009: Integrative pan-cancer genomic and transcriptomic analyses of refractory metastatic cancer

Metastatic relapse after treatment is the primary cause of cancer morbidity and mortality. While genetic mechanisms of primary tumors and, to a lesser extent, metastatic cancers have been studied in large cohorts, refractory metastatic tumors are not yet sufficiently characterized. Markers of aggressiveness and resistance that molecular profiling can extract from these tumors have yet to be identified and incorporated into clinical care. In this study, we present a pan-cancer cohort of 1,031 metastatic tumors (which we refer to as META-PRISM) that are resistant to at least one systemic therapy or with no approved treatment options. We retrieved the complete clinical history of patients and performed whole-exome (n=571) and transcriptome sequencing (n=947) for this cohort. The prevalence of detected cancer biomarkers was assessed and compared to an external cohort of primary tumors. In the META-PRISM cohort, we observed an increase in (i) whole-genome duplication frequency, (ii) tumor mutational burden, (iii) germline cancer-predisposing variants, and (iv) somatic alterations in cancer genes, including KRAS, EGFR, CCND1, MYC, and TP53, as compared to the tumor type-matched primary tumors. The most extensive increase in genomic variation at metastatic stage was observed in prostate cancer. We also identified enrichment of standard-of-care resistance biomarkers in most cancer types. However, only 7.6% of tumors harbored at least one such biomarker, indicating that the current understanding of resistance mechanisms remains insufficient. Our cohort demonstrated a significantly improved 6-month survival prediction from models incorporating molecular markers over models with only clinical markers for breast cancer patients and to a lesser extent for other studied tumor types. Overall, our data establish a unique resource for investigating treatment resistance mechanisms and performing predictive analyses in cancer.

Citation Format: Yoann Pradat, Julien Viot, Konstantin Gunbin, Andrei Iurchenko, Marc Deloger, Luigi Cerbone, Guillaume Grisay, Loic Verlingue, Veronique Scott, Stefan Michiels, Antoine Hollebecque, Gerome Jules-Clement, Antoine Laine, Luc Friboulet, Laura Mezquita, Yohann Loriot, Benjamin Besse, Fabrice Andre, Paul-Henry Cournede, Daniel Gautheret, Sergey Nikolaev. Integrative pan-cancer genomic and transcriptomic analyses of refractory metastatic cancer [abstract]. In: Proceedings of the AACR Special Conference: Cancer Metastasis; 2022 Nov 14-17; Portland, OR. Philadelphia (PA): AACR; Cancer Res 2022;83(2 Suppl_2):Abstract nr PR009.

A mitochondria-specific mutational signature of aging: increased rate of A > G substitutions on the heavy strand

The mutational spectrum of the mitochondrial DNA (mtDNA) does not resemble any of the known mutational signatures of the nuclear genome and variation in mtDNA mutational spectra between different organisms is still incomprehensible. Since mitochondria are responsible for aerobic respiration, it is expected that mtDNA mutational spectrum is affected by oxidative damage. Assuming that oxidative damage increases with age, we analyse mtDNA mutagenesis of different species in regards to their generation length. Analysing, (i) dozens of thousands of somatic mtDNA mutations in samples of different ages (ii) 70053 polymorphic synonymous mtDNA substitutions reconstructed in 424 mammalian species with different generation lengths and (iii) synonymous nucleotide content of 650 complete mitochondrial genomes of mammalian species we observed that the frequency of AH > GH substitutions (H: heavy strand notation) is twice bigger in species with high versus low generation length making their mtDNA more AH poor and GH rich. Considering that AH > GH substitutions are also sensitive to the time spent single-stranded (TSSS) during asynchronous mtDNA replication we demonstrated that AH > GH substitution rate is a function of both species-specific generation length and position-specific TSSS. We propose that AH > GH is a mitochondria-specific signature of oxidative damage associated with both aging and TSSS.

Mitochondrial Pseudogenes Suggest Repeated Inter-Species Hybridization among Direct Human Ancestors

The hypothesis that the evolution of humans involves hybridization between diverged species has been actively debated in recent years. We present the following novel evidence in support of this hypothesis: the analysis of nuclear pseudogenes of mtDNA ("NUMTs"). NUMTs are considered "mtDNA fossils" as they preserve sequences of ancient mtDNA and thus carry unique information about ancestral populations. Our comparison of a NUMT sequence shared by humans, chimpanzees, and gorillas with their mtDNAs implies that, around the time of divergence between humans and chimpanzees, our evolutionary history involved the interbreeding of individuals whose mtDNA had diverged as much as ~4.5 Myr prior. This large divergence suggests a distant interspecies hybridization. Additionally, analysis of two other NUMTs suggests that such events occur repeatedly. Our findings suggest a complex pattern of speciation in primate/human ancestors and provide one potential explanation for the mosaic nature of fossil morphology found at the emergence of the hominin lineage. A preliminary version of this manuscript was uploaded to the preprint server BioRxiv in 2017 (10.1101/134502).

Data on the time of integration of the human mitochondrial pseudogenes (NUMTs) into the nuclear genome

The data and methods presented in this article are supplementing the research article "Integration of mtDNA pseudogenes into the nuclear genome coincides with speciation of the human genus. A hypothesis", DOI: 10.1016/j.mito.2016.12.001 (Gunbin et al., 2017) [1]. Mitochondrial DNA is known to get inserted into nuclear DNA to form NUMTs, i.e. nuclear DNA pseudogenes of the mtDNA. We present here the sequences of selected NUMTs, in which time of integration can be determined with sufficient precision. We report their chromosomal positions , their position within the great ape mtDNA phylogeny, and their times of integration into the nuclear genome. The methods used to generate the data and to control their quality are also presented. The dataset is made publicly available to enable critical or extended analyzes.

Computer and Statistical Analysis of Transcription Factor Binding and Chromatin Modifications by ChIP-seq data in Embryonic Stem Cell

Advances in high throughput sequencing technology have enabled the identification of transcription factor (TF) binding sites in genome scale. TF binding studies are important for medical applications and stem cell research. Somatic cells can be reprogrammed to a pluripotent state by the combined introduction of factors such as Oct4, Sox2, c-Myc, Klf4. These reprogrammed cells share many characteristics with embryonic stem cells (ESCs) and are known as induced pluripotent stem cells (iPSCs). The signaling requirements for maintenance of human and murine embryonic stem cells (ESCs) differ considerably. Genome wide ChIP-seq TF binding maps in mouse stem cells include Oct4, Sox2, Nanog, Tbx3, Smad2 as well as group of other factors. ChIP-seq allows study of new candidate transcription factors for reprogramming. It was shown that Nr5a2 could replace Oct4 for reprogramming. Epigenetic modifications play important role in regulation of gene expression adding additional complexity to transcription network functioning. We have studied associations between different histone modification using published data together with RNA Pol II sites. We found strong associations between activation marks and TF binding sites and present it qualitatively. To meet issues of statistical analysis of genome ChIP-sequencing maps we developed computer program to filter out noise signals and find significant association between binding site affinity and number of sequence reads. The data provide new insights into the function of chromatin organization and regulation in stem cells.

Accumulation of slightly deleterious mutations in mitochondrial protein-coding genes of large versus small mammals

After the effective size of a population, N(e), declines, some slightly deleterious amino acid replacements which were initially suppressed by purifying selection become effectively neutral and can reach fixation. Here we investigate this phenomenon for a set of all 13 mitochondrial protein-coding genes from 110 mammalian species. By using body mass as a proxy for N(e), we show that large mammals (i.e., those with low N(e)) as compared with small ones (in our sample these are, on average, 369.5 kg and 275 g, respectively) have a 43% higher rate of accumulation of nonsynonymous nucleotide substitutions relative to synonymous substitutions, and an 8-40% higher rate of accumulation of radical amino acid substitutions relative to conservative substitutions, depending on the type of amino acid classification. These higher rates result in a 6% greater amino acid dissimilarity between modern species and their most recent reconstructed ancestors in large versus small mammals. Because nonsynonymous substitutions are likely to be more harmful than synonymous substitutions, and radical amino acid substitutions are likely to be more harmful than conservative ones, our results suggest that large mammals experience less efficient purifying selection than small mammals. Furthermore, because in the course of mammalian evolution body size tends to increase and, consequently, N(e) tends to decline, evolution of mammals toward large body size may involve accumulation of slightly deleterious mutations in mitochondrial protein-coding genes, which may contribute to decline or extinction of large mammals.