Stochastic survival of the densest and mitochondrial DNA clonal expansion in aging.

The expansion of mitochondrial DNA molecules with deletions has been associated with aging, particularly in skeletal muscle fibers; its mechanism has remained unclear for three decades. Previous accounts have assigned a replicative advantage (RA) to mitochondrial DNA containing deletion mutations, but there is also evidence that cells can selectively remove defective mitochondrial DNA. Here we present a spatial model that, without an RA, but instead through a combination of enhanced density for mutants and noise, produces a wave of expanding mutations with speeds consistent with experimental data. A standard model based on RA yields waves that are too fast. We provide a formula that predicts that wave speed drops with copy number, consonant with experimental data. Crucially, our model yields traveling waves of mutants even if mutants are preferentially eliminated. Additionally, we predict that mutant loads observed in single-cell experiments can be produced by de novo mutation rates that are drastically lower than previously thought for neutral models. Given this exemplar of how spatial structure (multiple linked mtDNA populations), noise, and density affect muscle cell aging, we introduce the mechanism of stochastic survival of the densest (SSD), an alternative to RA, that may underpin other evolutionary phenomena.

Mitochondrial DNA heteroplasmy is modulated during oocyte development propagating mutation transmission.

Heteroplasmic mitochondrial DNA (mtDNA) mutations are a common cause of inherited disease, but a few recurrent mutations account for the vast majority of new families. The reasons for this are not known. We studied heteroplasmic mice transmitting m.5024C>T corresponding to a human pathogenic mutation. Analyzing 1167 mother-pup pairs, we show that m.5024C>T is preferentially transmitted from low to higher levels but does not reach homoplasmy. Single-cell analysis of the developing mouse oocytes showed the preferential increase in mutant over wild-type mtDNA in the absence of cell division. A similar inheritance pattern is seen in human pedigrees transmitting several pathogenic mtDNA mutations. In m.5024C>T mice, this can be explained by the preferential propagation of mtDNA during oocyte maturation, counterbalanced by purifying selection against high heteroplasmy levels. This could explain how a disadvantageous mutation in a carrier increases to levels that cause disease but fails to fixate, causing multigenerational heteroplasmic mtDNA disorders.

Oxygen tension modulates the mitochondrial genetic bottleneck and influences the segregation of a heteroplasmic mtDNA variant in vitro.

Most humans carry a mixed population of mitochondrial DNA (mtDNA heteroplasmy) affecting ~1-2% of molecules, but rapid percentage shifts occur over one generation leading to severe mitochondrial diseases. A decrease in the amount of mtDNA within the developing female germ line appears to play a role, but other sub-cellular mechanisms have been implicated. Establishing an in vitro model of early mammalian germ cell development from embryonic stem cells, here we show that the reduction of mtDNA content is modulated by oxygen and reaches a nadir immediately before germ cell specification. The observed genetic bottleneck was accompanied by a decrease in mtDNA replicating foci and the segregation of heteroplasmy, which were both abolished at higher oxygen levels. Thus, differences in oxygen tension occurring during early development likely modulate the amount of mtDNA, facilitating mtDNA segregation and contributing to tissue-specific mutation loads.

Mitochondrial Network State Scales mtDNA Genetic Dynamics

Mitochondrial DNA (mtDNA) mutations cause severe congenital diseases but may also be associated with healthy aging. mtDNA is stochastically replicated and degraded, and exists within organelles which undergo dynamic fusion and fission. The role of the resulting mitochondrial networks in the time evolution of the cellular proportion of mutated mtDNA molecules (heteroplasmy), and cell-to-cell variability in heteroplasmy (heteroplasmy variance), remains incompletely understood. Heteroplasmy variance is particularly important since it modulates the number of pathological cells in a tissue. Here, we provide the first wide-reaching theoretical framework which bridges mitochondrial network and genetic states. We show that, under a range of conditions, the (genetic) rate of increase in heteroplasmy variance and de novo mutation are proportionally modulated by the (physical) fraction of unfused mitochondria, independently of the absolute fission-fusion rate. In the context of selective fusion, we show that intermediate fusion:fission ratios are optimal for the clearance of mtDNA mutants. Our findings imply that modulating network state, mitophagy rate, and copy number to slow down heteroplasmy dynamics when mean heteroplasmy is low could have therapeutic advantages for mitochondrial disease and healthy aging.

Frequency and signature of somatic variants in 1461 human brain exomes.

PURPOSE: To systematically study somatic variants arising during development in the human brain across a spectrum of neurodegenerative disorders. METHODS: In this study we developed a pipeline to identify somatic variants from exome sequencing data in 1461 diseased and control human brains. Eighty-eight percent of the DNA samples were extracted from the cerebellum. Identified somatic variants were validated by targeted amplicon sequencing and/or PyroMark® Q24. RESULTS: We observed somatic coding variants present in >10% of sampled cells in at least 1% of brains. The mutational signature of the detected variants showed a predominance of C>T variants most consistent with arising from DNA mismatch repair, occurred frequently in genes that are highly expressed within the central nervous system, and with a minimum somatic mutation rate of 4.25 × 10-10 per base pair per individual. CONCLUSION: These findings provide proof-of-principle that deleterious somatic variants can affect sizeable brain regions in at least 1% of the population, and thus have the potential to contribute to the pathogenesis of common neurodegenerative diseases.

High prevalence of focal and multi-focal somatic genetic variants in the human brain.

Somatic mutations during stem cell division are responsible for several cancers. In principle, a similar process could occur during the intense cell proliferation accompanying human brain development, leading to the accumulation of regionally distributed foci of mutations. Using dual platform >5000-fold depth sequencing of 102 genes in 173 adult human brain samples, we detect and validate somatic mutations in 27 of 54 brains. Using a mathematical model of neurodevelopment and approximate Bayesian inference, we predict that macroscopic islands of pathologically mutated neurons are likely to be common in the general population. The detected mutation spectrum also includes DNMT3A and TET2 which are likely to have originated from blood cell lineages. Together, these findings establish developmental mutagenesis as a potential mechanism for neurodegenerative disorders, and provide a novel mechanism for the regional onset and focal pathology in sporadic cases.

Oligogenic genetic variation of neurodegenerative disease genes in 980 postmortem human brains.

BACKGROUND: Several studies suggest that multiple rare genetic variants in genes causing monogenic forms of neurodegenerative disorders interact synergistically to increase disease risk or reduce the age of onset, but these studies have not been validated in large sporadic case series. METHODS: We analysed 980 neuropathologically characterised human brains with Alzheimer's disease (AD), Parkinson's disease-dementia with Lewy bodies (PD-DLB), frontotemporal dementia-amyotrophic lateral sclerosis (FTD-ALS) and age-matched controls. Genetic variants were assessed using the American College of Medical Genetics criteria for pathogenicity. Individuals with two or more variants within a relevant disease gene panel were defined as 'oligogenic'. RESULTS: The majority of oligogenic variant combinations consisted of a highly penetrant allele or known risk factor in combination with another rare but likely benign allele. The presence of oligogenic variants did not influence the age of onset or disease severity. After controlling for the single known major risk allele, the frequency of oligogenic variants was no different between cases and controls. CONCLUSIONS: A priori, individuals with AD, PD-DLB and FTD-ALS are more likely to harbour a known genetic risk factor, and it is the burden of these variants in combination with rare benign alleles that is likely to be responsible for some oligogenic associations. Controlling for this bias is essential in studies investigating a potential role for oligogenic variation in neurodegenerative diseases.

Mitochondrial Heterogeneity.

Cell-to-cell heterogeneity drives a range of (patho)physiologically important phenomena, such as cell fate and chemotherapeutic resistance. The role of metabolism, and particularly of mitochondria, is increasingly being recognized as an important explanatory factor in cell-to-cell heterogeneity. Most eukaryotic cells possess a population of mitochondria, in the sense that mitochondrial DNA (mtDNA) is held in multiple copies per cell, where the sequence of each molecule can vary. Hence, intra-cellular mitochondrial heterogeneity is possible, which can induce inter-cellular mitochondrial heterogeneity, and may drive aspects of cellular noise. In this review, we discuss sources of mitochondrial heterogeneity (variations between mitochondria in the same cell, and mitochondrial variations between supposedly identical cells) from both genetic and non-genetic perspectives, and mitochondrial genotype-phenotype links. We discuss the apparent homeostasis of mtDNA copy number, the observation of pervasive intra-cellular mtDNA mutation (which is termed "microheteroplasmy"), and developments in the understanding of inter-cellular mtDNA mutation ("macroheteroplasmy"). We point to the relationship between mitochondrial supercomplexes, cristal structure, pH, and cardiolipin as a potential amplifier of the mitochondrial genotype-phenotype link. We also discuss mitochondrial membrane potential and networks as sources of mitochondrial heterogeneity, and their influence upon the mitochondrial genome. Finally, we revisit the idea of mitochondrial complementation as a means of dampening mitochondrial genotype-phenotype links in light of recent experimental developments. The diverse sources of mitochondrial heterogeneity, as well as their increasingly recognized role in contributing to cellular heterogeneity, highlights the need for future single-cell mitochondrial measurements in the context of cellular noise studies.

Mitochondrial DNA density homeostasis accounts for a threshold effect in a cybrid model of a human mitochondrial disease.

Mitochondrial dysfunction is involved in a wide array of devastating diseases, but the heterogeneity and complexity of the symptoms of these diseases challenges theoretical understanding of their causation. With the explosion of omics data, we have the unprecedented opportunity to gain deep understanding of the biochemical mechanisms of mitochondrial dysfunction. This goal raises the outstanding need to make these complex datasets interpretable. Quantitative modelling allows us to translate such datasets into intuition and suggest rational biomedical treatments. Taking an interdisciplinary approach, we use a recently published large-scale dataset and develop a descriptive and predictive mathematical model of progressive increase in mutant load of the MELAS 3243A>G mtDNA mutation. The experimentally observed behaviour is surprisingly rich, but we find that our simple, biophysically motivated model intuitively accounts for this heterogeneity and yields a wealth of biological predictions. Our findings suggest that cells attempt to maintain wild-type mtDNA density through cell volume reduction, and thus power demand reduction, until a minimum cell volume is reached. Thereafter, cells toggle from demand reduction to supply increase, up-regulating energy production pathways. Our analysis provides further evidence for the physiological significance of mtDNA density and emphasizes the need for performing single-cell volume measurements jointly with mtDNA quantification. We propose novel experiments to verify the hypotheses made here to further develop our understanding of the threshold effect and connect with rational choices for mtDNA disease therapies.

Mitochondrial heterogeneity, metabolic scaling and cell death.

Heterogeneity in mitochondrial content has been previously suggested as a major contributor to cellular noise, with multiple studies indicating its direct involvement in biomedically important cellular phenomena. A recently published dataset explored the connection between mitochondrial functionality and cell physiology, where a non-linearity between mitochondrial functionality and cell size was found. Using mathematical models, we suggest that a combination of metabolic scaling and a simple model of cell death may account for these observations. However, our findings also suggest the existence of alternative competing hypotheses, such as a non-linearity between cell death and cell size. While we find that the proposed non-linear coupling between mitochondrial functionality and cell size provides a compelling alternative to previous attempts to link mitochondrial heterogeneity and cell physiology, we emphasise the need to account for alternative causal variables, including cell cycle, size, mitochondrial density and death, in future studies of mitochondrial physiology.