mtDNA Population Research Articles

Mitochondrial Control Region Database of Hungarian Fallow Deer (Dama dama) Populations for Forensic Use.

The evidential value of an mtDNA match between biological remains and their potential donor is determined by the random match probability of the haplotype. This probability is based on the haplotype's population frequency estimate. Consequently, implementing a population study representative of the population relevant to a forensic case is vital to correctly evaluating the evidence. The emerging number of poaching cases and the limited availability of such data emphasizes the need for an improved fallow deer mtDNA population databank for forensic purposes, including targeting the entire mitochondrial control region. By sequencing a 945-base-pair-long segment of the mitochondrial control region in 138 animals from five populations in Hungary, we found four different haplotypes, including one which had not yet been described. Our results, supplemented with data already available from previous research, do not support the possibility of determining the population of origin, although some patterns of geographical separation can be distinguished. Estimates of molecular diversity indicate similarly low mtDNA diversity (Hd = 0.565 and π = 0.002) compared to data from other countries. The calculated random match probability of 0.547 shows a high probability of coincidence and, therefore, a limited capacity for exclusion. Our results indicate that despite the overall low genetic diversity of mtDNA within the Hungarian fallow deer samples, a pattern of differentiation among the regions is present, which can have relevance from a forensic point of view.

Simulating the effects of kinship and postmarital residence patterns on mitochondrial DNA diversity in mortuary contexts.

We explore the observable outcome in mtDNA diversity of different kinship systems and associated postmarital residence patterns in the archeological record, using simulations at the intrapopulation level. Four kinship systems were simulated from a set of variable fertility and mortality scenarios. Initial conditions consisted of six clusters of variable size and a random number of assigned haplotypes, with individuals migrating between groups and reproducing for 15 generations. Each 15-generation span was simulated 500 times to obtain representative intrasite mtDNA diversity distributions for each kinship system. Additional simulations were devised to consider the effect of migration and different rates of adherence to kinship norms. Matrilineal kinship generates low male and female haplotype diversities that are statistically indistinguishable from each other, while female diversity in bilateral kinship with matrilocality is significantly lower than that observed in male diversity. Furthermore, mtDNA diversity generated by patrilineal kinship is very high. The effect of noncompliance with kinship rules is low; migration has a considerable impact on diversity, eventually obscuring the effect of kinship practices. The results of the simulations can be applied to ancient mtDNA data from archeological contexts, as exemplified with data from two studies. On a broader scale, the kinship system followed by the sampled population, can lead to either over- or underestimation of mtDNA population diversity. The results of the simulations can be used in the design of inferential frameworks to discern kinship scenarios in the archeological record, based on mtDNA and other types of evidence.

Mitochondrial network structure controls cell-to-cell mtDNA variability generated by cell divisions.

Mitochondria are highly dynamic organelles, containing vital populations of mitochondrial DNA (mtDNA) distributed throughout the cell. Mitochondria form diverse physical structures in different cells, from cell-wide reticulated networks to fragmented individual organelles. These physical structures are known to influence the genetic makeup of mtDNA populations between cell divisions, but their influence on the inheritance of mtDNA at divisions remains less understood. Here, we use statistical and computational models of mtDNA content inside and outside the reticulated network to quantify how mitochondrial network structure can control the variances of inherited mtDNA copy number and mutant load. We assess the use of moment-based approximations to describe heteroplasmy variance and identify several cases where such an approach has shortcomings. We show that biased inclusion of one mtDNA type in the network can substantially increase heteroplasmy variance (acting as a genetic bottleneck), and controlled distribution of network mass and mtDNA through the cell can conversely reduce heteroplasmy variance below a binomial inheritance picture. Network structure also allows the generation of heteroplasmy variance while controlling copy number inheritance to sub-binomial levels, reconciling several observations from the experimental literature. Overall, different network structures and mtDNA arrangements within them can control the variances of key variables to suit a palette of different inheritance priorities.

Mitochondrial genome of the jack silverside, Atherinopsis californiensis (Atherinopsidae, Atheriniformes), a nearshore fish of the California Current Ecosystem

The jack silverside (Atherinopsis californiensis), also referred to as jacksmelt, is a neotropical silverside fish that inhabits nearshore shallow waters of the California Current Ecosystem in the Northeast Pacific Ocean, ranging from the coast of Oregon, USA, in the north to as far south as Baja California, Mexico. This fish is the sole member of its genus and is a commonly taken species when hook-and-line fishing in pelagic-neritic environments including bays, estuaries, kelp forests, and along sand beaches. Here we report the first complete mitochondrial genome of jack silverside consisting of 16,519 bp nucleotides and encoding 13 protein-coding regions, 12S and 16S rRNAs, 22 tRNAs, and an 841 bp D-loop control region. Phylogenetic analysis using all protein-coding genes of the complete mitogenome confirmed the inclusion of A. californiensis within subfamily Atherinopsinae of family Atherinopsidae, order Atheriniformes. This complete mitochondrial DNA genome will be of use for biodiversity assessments in the California Current ecosystem, while also providing a foundation for future mtDNA population genetic studies on this prominently caught species in shore- and pier-based recreational sport fishing.

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.

Independent regulation of mitochondrial DNA quantity and quality in Caenorhabditis elegans primordial germ cells.

Mitochondria harbor an independent genome, called mitochondrial DNA (mtDNA), which contains essential metabolic genes. Although mtDNA mutations occur at high frequency, they are inherited infrequently, indicating that germline mechanisms limit their accumulation. To determine how germline mtDNA is regulated, we examined the control of mtDNA quantity and quality in C. elegans primordial germ cells (PGCs). We show that PGCs combine strategies to generate a low point in mtDNA number by segregating mitochondria into lobe-like protrusions that are cannibalized by adjacent cells, and by concurrently eliminating mitochondria through autophagy, reducing overall mtDNA content twofold. As PGCs exit quiescence and divide, mtDNAs replicate to maintain a set point of ~200 mtDNAs per germline stem cell. Whereas cannibalism and autophagy eliminate mtDNAs stochastically, we show that the kinase PTEN-induced kinase 1 (PINK1), operating independently of Parkin and autophagy, preferentially reduces the fraction of mutant mtDNAs. Thus, PGCs employ parallel mechanisms to control both the quantity and quality of the founding population of germline mtDNAs.

NOT ALL CELL-FREE MITOCHONDRIAL DNA IS EQUAL IN TRAUMA PATIENTS.

ABSTRACTMitochondrial DNA (mtDNA) acts as a proinflammatory damage-associated molecular pattern that stimulates innate immune activation via Toll-like receptor 9, similarly to bacterial DNA. A number of clinical studies have measured elevated cell-free mtDNA in the plasma of trauma patients, thought to originate from tissue injury and inflammatory processes; however, the magnitude of this increase, the absolute concentration, and the association with poor outcomes varies considerably across studies. Measurements of cell-free mtDNA in healthy individuals have shown that the majority of “cell-free” mtDNA (>95%) can be centrifuged/filtered from plasma in the size range of 0.45 to 5 μm, suggesting that there are larger forms of mtDNA-containing complexes in the plasma that could be considered cell-free. Whether this is true for trauma patients (and other relevant disease states) and the clinical relevance of the larger forms of mtDNA is unknown. These findings from healthy individuals also suggest that the centrifugation speeds used to generate cell-free plasma (which are rarely consistent among studies) could result in mixed populations of cell-free mtDNA that could confound associations with outcomes. We demonstrate in this study of 25 major trauma patients that the majority of the cell-free mtDNA in trauma patient plasma (>95%) is removed after centrifugation at 16,000g. Despite the larger forms of mtDNA being predominant, they do not correlate with outcomes or expected parameters such as injury/shock severity, multiple organ failure, and markers of inflammation, whereas low-molecular-weight cell-free mtDNA correlates strongly with these variables.

Mitochondrial DNA variation in Parkinson's disease: Analysis of "out-of-place" population variants as a risk factor.

Mitochondrial DNA (mtDNA), a potential source of mitochondrial dysfunction, has been implicated in Parkinson’s disease (PD). However, many previous studies investigating associations between mtDNA population variation and PD reported inconsistent or contradictory findings. Here, we investigated an alternative hypothesis to determine whether mtDNA variation could play a significant role in PD risk. Emerging evidence suggests that haplogroup-defining mtDNA variants may have pathogenic potential if they occur “out-of-place” on a different maternal lineage. We hypothesized that the mtDNA of PD cases would be enriched for out-of-place variation in genes encoding components of the oxidative phosphorylation complexes. We tested this hypothesis with a unique dataset comprising whole mitochondrial genomes of 70 African ancestry PD cases, two African ancestry control groups (n = 78 and n = 53) and a replication group of 281 European ancestry PD cases and 140 controls from the Parkinson’s Progression Markers Initiative cohort. Significantly more African ancestry PD cases had out-of-place variants than controls from the second control group (P < 0.0125), although this association was not observed in the first control group nor the replication group. As the first mtDNA study to include African ancestry PD cases and to explore out-of-place variation in a PD context, we found evidence that such variation might be significant in this context, thereby warranting further replication in larger cohorts.

Altered collective mitochondrial dynamics in the Arabidopsis msh1 mutant compromising organelle DNA maintenance.

Mitochondria form highly dynamic populations in the cells of plants (and almost all eukaryotes). The characteristics and benefits of this collective behaviour, and how it is influenced by nuclear features, remain to be fully elucidated. Here, we use a recently developed quantitative approach to reveal and analyse the physical and collective 'social' dynamics of mitochondria in an Arabidopsis msh1 mutant where the organelle DNA maintenance machinery is compromised. We use a newly created line combining the msh1 mutant with mitochondrially targeted green fluorescent protein (GFP), and characterize mitochondrial dynamics with a combination of single-cell time-lapse microscopy, computational tracking, and network analysis. The collective physical behaviour of msh1 mitochondria is altered from that of the wild type in several ways: mitochondria become less evenly spread, and networks of inter-mitochondrial encounters become more connected, with greater potential efficiency for inter-organelle exchange-reflecting a potential compensatory mechanism for the genetic challenge to the mitochondrial DNA population, supporting more inter-organelle exchange. We find that these changes are similar to those observed in friendly, where mitochondrial dynamics are altered by a physical perturbation, suggesting that this shift to higher connectivity may reflect a general response to mitochondrial challenges, where physical dynamics of mitochondria may be altered to control the genetic structure of the mtDNA population.

A world of taxonomic pain: cryptic species, inexplicable host-specificity, and host-induced morphological variation among species of Bivesicula Yamaguti, 1934 (Trematoda: Bivesiculidae) from Indo-Pacific Holocentridae, Muraenidae and Serranidae.

The taxonomy of species of Bivesicula Yamaguti, 1934 is analysed for samples from holocentrid, muraenid and serranid fishes from Japan, Ningaloo Reef (Western Australia), the Great Barrier Reef (Queensland), New Caledonia and French Polynesia. Analysis of three genetic markers (cox1 mtDNA, ITS2 and 28S rDNA) identifies three strongly supported clades of species and suggests that Bivesicula as presently recognized is not monophyletic. On the basis of combined morphological, molecular and biological data, 10 species are distinguished of which five are proposed as new. Bivesicula Clade 1 comprises seven species of which three are effectively morphologically cryptic relative to each other; all seven infect serranids and four also infect holocentrids. Bivesicula Clade 2 comprises three species of which two are effectively morphologically cryptic relative to each other; all three infect serranids and one also infects a muraenid. Bivesicula Clade 3 comprises two known species from apogonids and a pomacentrid, and forms a clade with species of Paucivitellosus Coil, Reid & Kuntz, 1965 to the exclusion of other Bivesicula species. Taxonomy in this genus is made challenging by the combination of low resolving power of ribosomal markers, the existence of regional cox1 mtDNA populations, exceptional and unpredictable host-specificity and geographical distribution, and significant host-induced morphological variation.

Mitochondrial DNA variation across 56,434 individuals in gnomAD.

Genomic databases of allele frequency are extremely helpful for evaluating clinical variants of unknown significance; however, until now, databases such as the Genome Aggregation Database (gnomAD) have focused on nuclear DNA and have ignored the mitochondrial genome (mtDNA). Here, we present a pipeline to call mtDNA variants that addresses three technical challenges: (1) detecting homoplasmic and heteroplasmic variants, present, respectively, in all or a fraction of mtDNA molecules; (2) circular mtDNA genome; and (3) misalignment of nuclear sequences of mitochondrial origin (NUMTs). We observed that mtDNA copy number per cell varied across gnomAD cohorts and influenced the fraction of NUMT-derived false-positive variant calls, which can account for the majority of putative heteroplasmies. To avoid false positives, we excluded contaminated samples, cell lines, and samples prone to NUMT misalignment due to few mtDNA copies. Furthermore, we report variants with heteroplasmy ≥10%. We applied this pipeline to 56,434 whole-genome sequences in the gnomAD v3.1 database that includes individuals of European (58%), African (25%), Latino (10%), and Asian (5%) ancestry. Our gnomAD v3.1 release contains population frequencies for 10,850 unique mtDNA variants at more than half of all mtDNA bases. Importantly, we report frequencies within each nuclear ancestral population and mitochondrial haplogroup. Homoplasmic variants account for most variant calls (98%) and unique variants (85%). We observed that 1/250 individuals carry a pathogenic mtDNA variant with heteroplasmy above 10%. These mtDNA population allele frequencies are freely accessible and will aid in diagnostic interpretation and research studies.

Investigation of the correlation between mildly deleterious mtDNA Variations and the clinical progression of multiple sclerosis

BackgroundEvidence suggests that mitochondrial DNA (mtDNA) variation at a population level may influence susceptibility to, or the clinical progression of Multiple Sclerosis (MS). ObjectiveTo determine if mtDNA population variation is linked to the clinical progress of MS. MethodsUsing the complete mtDNA sequences of 217 MS patients, we applied the new ‘variant load’ model, designed as a framework by which to examine the role of mtDNA variation in the context of complex clinical disease. ResultsNo significant association was detected between mtDNA ‘variant load’and the clinical measures of progression. ConclusionOur results suggest that mtDNA population variation does not play a substantial role in the clinical progression of MS; however, modest effects and/or effects in a subgroup of patients cannot be entirely excluded. Results do not exclude the possibility of detecting an association between variation and more strictly quantified variables obtained from histopathologically-stained specimens. The results further illustrate the method's applicabilityto other disease phenotypes.

A Naturally Heteroplasmic Clam Provides Clues about the Effects of Genetic Bottleneck on Paternal mtDNA.

Mitochondrial DNA (mtDNA) is present in multiple copies within an organism. Since these copies are not identical, a single individual carries a heterogeneous population of mtDNAs, a condition known as heteroplasmy. Several factors play a role in the dynamics of the within-organism mtDNA population: among them, genetic bottlenecks, selection, and strictly maternal inheritance are known to shape the levels of heteroplasmy across mtDNAs.In Metazoa, the only evolutionarily stable exception to the strictly maternal inheritance of mitochondria is the doubly uniparental inheritance (DUI), reported in 100+ bivalve species. In DUI species, there are two highly divergent mtDNA lineages, one inherited through oocyte mitochondria (F-type) and the other through sperm mitochondria (M-type). Having both parents contributing to the mtDNA pool of the progeny makes DUI a unique system to study the dynamics of mtDNA populations. Since, in bivalves, the spermatozoon has few mitochondria (4–5), M-type mtDNA faces a tight bottleneck during embryo segregation, one of the narrowest mitochondrial bottlenecks investigated so far.Here, we analyzed the F- and M-type mtDNA variability within individuals of the DUI species Ruditapes philippinarum and investigated for the first time the effects of such a narrow bottleneck affecting mtDNA populations. As a potential consequence of this narrow bottleneck, the M-type mtDNA shows a large variability in different tissues, a condition so pronounced that it leads to genotypes from different tissues of the same individual not to cluster together. We believe that such results may help understanding the effect of low population size on mtDNA bottleneck.

The tale of springs and streams: how different aquatic ecosystems impacted the mtDNA population structure of two riffle beetles in the Western Carpathians

The Western Carpathians are a particularly interesting part of the Carpathian Arc. According to recent molecular data upon aquatic and terrestrial taxa, this mountain area is an important biodiversity hotspot of Europe. Moreover, the W Carpathians include rich systems of karst springs inhabited by specific fauna, where molecular diversity and phylogeographic patterns are yet to be fully explored. Our study aims to compare population genetic structure and molecular diversity of two related and commonly co-occurring riffle beetles, Elmis aenea (PWJ Müller, 1806) and Limnius perrisi (Dufour, 1843) in the springs and streams of the W Carpathians using the mitochondrial DNA barcoding fragment of the cytochrome c oxidase subunit I gene (COI). The relatively stable thermal and chemical conditions of springs throughout unfavourable climatic settings make these highly specific lotic systems potentially ideal for a long-term survival of some aquatic biota. Populations of both elmid species were relatively homogeneous genetically, with a single dominant haplotype. However, we revealed that E. aenea significantly dominated in the springs, while L. perrisi preferred streams. Relative isolation of the springs and their stable conditions were reflected in significantly higher molecular diversity of the E. aenea population in comparison to L. perrisi. The results of Bayesian Skyline Plot analysis also indicated the exceptional position of springs regarding maintaining the population size of E. aenea. On the other hand, it seems that streams in the W Carpathians provide more effective dispersal channels for L. perrisi, whose population expanded much earlier compared to E. aenea. Present study points out that different demographic histories of these two closely related elmid species are manifested by their different habitat preference and molecular diversity.

Chromosomal Differentiation in Genetically Isolated Populations of the Marsh-Specialist Crocidura suaveolens (Mammalia: Soricidae).

The genus Crocidura represents a remarkable model for the study of chromosome evolution. This is the case of the lesser white-toothed shrew (Crocidura suaveolens), a representative of the Palearctic group. Although continuously distributed from Siberia to Central Europe, C. suaveolens is a rare, habitat-specialist species in the southwesternmost limit of its distributional range, in the Gulf of Cádiz (Iberian Peninsula). In this area, C. suaveolens is restricted to genetically isolated populations associated to the tidal marches of five rivers (Guadiana, Piedras, Odiel, Tinto and Guadalquivir). This particular distributional range provides a unique opportunity to investigate whether genetic differentiation and habitat specialization was accompanied by chromosomal variation. In this context, the main objective of this study was to determinate the chromosomal characteristics of the habitat-specialist C. suaveolens in Southwestern Iberia, as a way to understand the evolutionary history of this species in the Iberian Peninsula. A total of 41 individuals from six different populations across the Gulf of Cádiz were collected and cytogenetically characterized. We detected four different karyotypes, with diploid numbers (2n) ranging from 2n = 40 to 2n = 43. Two of them (2n = 41 and 2n = 43) were characterized by the presence of B-chromosomes. The analysis of karyotype distribution across lineages and populations revealed an association between mtDNA population divergence and chromosomal differentiation. C. suaveolens populations in the Gulf of Cádiz provide a rare example of true karyotypic polymorphism potentially associated to genetic isolation and habitat specialization in which to investigate the evolutionary significance of chromosomal variation in mammals and their contribution to phenotypic and ecological divergence.

Mitochondrial DNA variation in Sub-Saharan Africa: Forensic data from a mixed West African sample, Côte d’Ivoire (Ivory Coast), and Rwanda

This study provides 398 novel complete mitochondrial control region sequences that augment the still underrepresented data from Africa by three datasets: a mixed West African sample set deriving from 12 countries (n = 145) and datasets from Côte d’Ivoire (Ivory Coast) (n = 100) as well as Rwanda (n = 153). The analysis of mtDNA variation and genetic comparisons with published data revealed low random match probabilities in all three datasets and typical West African and East African diversity, respectively. Genetic parameters indicate that the presented mixed West African dataset may serve as first forensic mtDNA control region database for West Africa in general. In addition, a strategy for responsible forensic application of precious mtDNA population samples potentially containing close maternal relatives is outlined. The datasets will be uploaded to the forensic mtDNA database EMPOP (https://empop.online) upon publication.

Population Dynamics of Mitochondria in Cells: A Minimal Mathematical Model

Mitochondria are cellular organelles that play an important role in cellular energy production, signaling, differentiation, and apoptosis. All organisms are believed to have low levels of variation in mitochondrial DNA (mtDNA), and repeated mitotic segregation and clonal expansion can enable a mitochondrion to eventually dominate the mtDNA pool. Alterations in mtDNA are connected to a range of human diseases. Therefore, understanding how changes in mtDNA accumulate over time and are correlated to changes in mitochondrial function can have a profound impact on our understanding of fundamental cell biophysics and the origins of some human diseases. Motivated by this, we develop and study a mathematical model to determine which cellular parameters have the largest impact on mtDNA population dynamics. The model consists of coupled differential equations to describe populations of healthy and dysfunctional mitochondria subject to mitochondrial fission, fusion, autophagy, and varying levels of cellular ATP. We study the time evolution of each population under specific selection biases and obtain a heat map in the parameter space of the ratio of the rates of fusion and autophagy of the healthy and dysfunctional populations. Our results may provide insights into how different mitochondrial populations survive and evolve under different selection pressures and with time.

Evolving mtDNA populations within cells.

Mitochondrial DNA (mtDNA) encodes vital respiratory machinery. Populations of mtDNA molecules exist in most eukaryotic cells, subject to replication, degradation, mutation, and other population processes. These processes affect the genetic makeup of cellular mtDNA populations, changing cell-to-cell distributions, means, and variances of mutant mtDNA load over time. As mtDNA mutant load has nonlinear effects on cell functionality, and cell functionality has nonlinear effects on tissue performance, these statistics of cellular mtDNA populations play vital roles in health, disease, and inheritance. This mini review will describe some of the better-known ways in which these populations change over time in different organisms, highlighting the importance of quantitatively understanding both mutant load mean and variance. Due to length constraints, we cannot attempt to be comprehensive but hope to provide useful links to some of the many excellent studies on these topics.

Energetic costs of cellular and therapeutic control of stochastic mitochondrial DNA populations.

The dynamics of the cellular proportion of mutant mtDNA molecules is crucial for mitochondrial diseases. Cellular populations of mitochondria are under homeostatic control, but the details of the control mechanisms involved remain elusive. Here, we use stochastic modelling to derive general results for the impact of cellular control on mtDNA populations, the cost to the cell of different mtDNA states, and the optimisation of therapeutic control of mtDNA populations. This formalism yields a wealth of biological results, including that an increasing mtDNA variance can increase the energetic cost of maintaining a tissue, that intermediate levels of heteroplasmy can be more detrimental than homoplasmy even for a dysfunctional mutant, that heteroplasmy distribution (not mean alone) is crucial for the success of gene therapies, and that long-term rather than short intense gene therapies are more likely to beneficially impact mtDNA populations.

MtDNA population variation in Myalgic encephalomyelitis/Chronic fatigue syndrome in two populations: a study of mildly deleterious variants

Myalgic Encephalomyelitis (ME), also known as Chronic Fatigue Syndrome (CFS) is a debilitating condition. There is growing interest in a possible etiologic or pathogenic role of mitochondrial dysfunction and mitochondrial DNA (mtDNA) variation in ME/CFS. Supporting such a link, fatigue is common and often severe in patients with mitochondrial disease. We investigate the role of mtDNA variation in ME/CFS. No proven pathogenic mtDNA mutations were found. We then investigated population variation. Two cohorts were analysed, one from the UK (n = 89 moderately affected; 29 severely affected) and the other from South Africa (n = 143 moderately affected). For both cohorts, ME/CFS patients had an excess of individuals without a mildly deleterious population variant. The differences in population variation might reflect a mechanism important to the pathophysiology of ME/CFS.