Maternal age effect and severe germ-line bottleneck in the inheritance of human mitochondrial DNA

Strong within-host selection in a maternally inherited obligate symbiont: Buchnera and aphids

Numerous animal lineages have maternally inherited symbionts that are required for host reproduction and growth. Endosymbionts also pose a risk to their hosts because of the mutational decay of their genomes through genetic drift or to selfish mutations that favor symbiont fitness over host fitness. One model for heritable endosymbiosis is the association of aphids with their obligate bacterial symbiont, Buchnera We experimentally established heteroplasmic pea aphid matrilines containing pairs of closely related Buchnera haplotypes and used deep sequencing of diagnostic markers to measure haplotype frequencies in successive host generations. These frequencies were used to estimate the effective population size of Buchnera within hosts (i.e., the transmission bottleneck size) and the extent of within-host selection. The within-host effective population size was in the range of 10 to 20, indicating a strong potential for genetic drift and fixation of deleterious mutations. Remarkably, closely related haplotypes were subject to strong within-host selection, with selection coefficients as high as 0.5 per aphid generation. In one case, the direction of selection depended on the thermal environment and went in the same direction as between-host selection. In another, a new mutant haplotype had a strong within-host advantage under both environments but had no discernible effect on host-level fitness under laboratory conditions. Thus, within-host selection can be strong, resulting in a rapid fixation of mutations with little impact on host-level fitness. Together, these results show that within-host selection can drive evolution of an obligate symbiont, accelerating sequence evolution.

The need for high-quality oocyte mitochondria at extreme ploidy dictates mammalian germline development

Selection against deleterious mitochondrial mutations is facilitated by germline processes, lowering the risk of genetic diseases. How selection works is disputed: experimental data are conflicting and previous modeling work has not clarified the issues; here, we develop computational and evolutionary models that compare the outcome of selection at the level of individuals, cells and mitochondria. Using realistic de novo mutation rates and germline development parameters from mouse and humans, the evolutionary model predicts the observed prevalence of mitochondrial mutations and diseases in human populations. We show the importance of organelle-level selection, seen in the selective pooling of mitochondria into the Balbiani body, in achieving high-quality mitochondria at extreme ploidy in mature oocytes. Alternative mechanisms debated in the literature, bottlenecks and follicular atresia, are unlikely to account for the clinical data, because neither process effectively eliminates mitochondrial mutations under realistic conditions. Our findings explain the major features of female germline architecture, notably the longstanding paradox of over-proliferation of primordial germ cells followed by massive loss. The near-universality of these processes across animal taxa makes sense in light of the need to maintain mitochondrial quality at extreme ploidy in mature oocytes, in the absence of sex and recombination.

mtDNA Heteroplasmy: Origin, Detection, Significance, and Evolutionary Consequences

Mitochondrial DNA (mtDNA) is predominately uniparentally transmitted. This results in organisms with a single type of mtDNA (homoplasmy), but two or more mtDNA haplotypes have been observed in low frequency in several species (heteroplasmy). In this review, we aim to highlight several aspects of heteroplasmy regarding its origin and its significance on mtDNA function and evolution, which has been progressively recognized in the last several years. Heteroplasmic organisms commonly occur through somatic mutations during an individual’s lifetime. They also occur due to leakage of paternal mtDNA, which rarely happens during fertilization. Alternatively, heteroplasmy can be potentially inherited maternally if an egg is already heteroplasmic. Recent advances in sequencing techniques have increased the ability to detect and quantify heteroplasmy and have revealed that mitochondrial DNA copies in the nucleus (NUMTs) can imitate true heteroplasmy. Heteroplasmy can have significant evolutionary consequences on the survival of mtDNA from the accumulation of deleterious mutations and for its coevolution with the nuclear genome. Particularly in humans, heteroplasmy plays an important role in the emergence of mitochondrial diseases and determines the success of the mitochondrial replacement therapy, a recent method that has been developed to cure mitochondrial diseases.

Maternal age at birth and daughter's fecundability

STUDY QUESTION

Do daughters of older mothers have lower fecundability?

SUMMARY ANSWER

In this cohort study of North American pregnancy planners, there was virtually no association between maternal age ≥35 years and daughters' fecundability.

WHAT IS KNOWN ALREADY

Despite suggestive evidence that daughters of older mothers may have lower fertility, only three retrospective studies have examined the association between maternal age and daughter's fecundability.

STUDY DESIGN, SIZE, DURATION

Prospective cohort study of 6689 pregnancy planners enrolled between March 2016 and January 2020.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Pregnancy Study Online (PRESTO) is an ongoing pre-conception cohort study of pregnancy planners (age, 21-45 years) from the USA and Canada. We estimated fecundability ratios (FR) for maternal age at the participant's birth using multivariable proportional probabilities regression models.

MAIN RESULTS AND THE ROLE OF CHANCE

Daughters of mothers ≥30 years were less likely to have previous pregnancies (or pregnancy attempts) or risk factors for infertility, although they were more likely to report that their mother had experienced problems conceiving. The proportion of participants with prior unplanned pregnancies, a birth before age 21, ≥3 cycles of attempt at study entry or no follow-up was greater among daughters of mothers <25 years. Compared with maternal age 25-29 years, FRs (95% CI) for maternal age <20, 20-24, 30-34, and ≥35 were 0.72 (0.61, 0.84), 0.92 (0.85, 1.00), 1.08 (1.00, 1.17), and 1.00 (0.89, 1.12), respectively.

LIMITATIONS, REASONS FOR CAUTION

Although the examined covariates did not meaningfully affect the associations, we had limited information on the participants' mother. Differences by maternal age in reproductive history, infertility risk factors and loss to follow-up suggest that selection bias may partly explain our results.

WIDER IMPLICATIONS OF THE FINDINGS

Our finding that maternal age 35 years or older was not associated with daughter's fecundability is reassuring, considering the trend towards delayed childbirth. However, having been born to a young mother may be a marker of low fecundability among pregnancy planners.

STUDY FUNDING/COMPETING INTEREST(S)

PRESTO was funded by NICHD Grants (R21-HD072326 and R01-HD086742) and has received in-kind donations from Swiss Precision Diagnostics, FertilityFriend.com, Kindara.com, and Sandstone Diagnostics. Dr Wise is a fibroid consultant for AbbVie, Inc.

TRIAL REGISTRATION NUMBER

n/a.

Placental mitochondrial DNA mutational load and perinatal outcomes: Findings from a multi-ethnic pregnancy cohort

Mitochondria fuel placental activity, with mitochondrial dysfunction implicated in several perinatal complications. We investigated placental mtDNA mutational load using NextGen sequencing in relation to birthweight and gestational length among 358 mother-newborn pairs. We found that higher heteroplasmy, especially in the hypervariable displacement loop region, was associated with shorter gestational length. Results were similar among male and female pregnancies, but stronger in magnitude among females. With regard to growth, we observed that higher mutational load was associated with lower birthweight-for-gestational age (BWGA) among females, but higher BWGA among males. These findings support potential sex-differential fetal biological strategies for coping with increased heteroplasmies.

Mitochondrial function in development and disease

Mitochondria are organelles with vital functions in almost all eukaryotic cells. Often described as the cellular 'powerhouses' due to their essential role in aerobic oxidative phosphorylation, mitochondria perform many other essential functions beyond energy production. As signaling organelles, mitochondria communicate with the nucleus and other organelles to help maintain cellular homeostasis, allow cellular adaptation to diverse stresses, and help steer cell fate decisions during development. Mitochondria have taken center stage in the research of normal and pathological processes, including normal tissue homeostasis and metabolism, neurodegeneration, immunity and infectious diseases. The central role that mitochondria assume within cells is evidenced by the broad impact of mitochondrial diseases, caused by defects in either mitochondrial or nuclear genes encoding for mitochondrial proteins, on different organ systems. In this Review, we will provide the reader with a foundation of the mitochondrial 'hardware', the mitochondrion itself, with its specific dynamics, quality control mechanisms and cross-organelle communication, including its roles as a driver of an innate immune response, all with a focus on development, disease and aging. We will further discuss how mitochondrial DNA is inherited, how its mutation affects cell and organismal fitness, and current therapeutic approaches for mitochondrial diseases in both model organisms and humans.

Mitochondrial DNA segregation and replication restrict the transmission of detrimental mutation

Although mitochondrial DNA (mtDNA) is prone to accumulate mutations and lacks conventional DNA repair mechanisms, deleterious mutations are exceedingly rare. How the transmission of detrimental mtDNA mutations is restricted through the maternal lineage is debated. Here, we demonstrate that mitochondrial fission, together with the lack of mtDNA replication, segregate mtDNA into individual organelles in the Drosophila early germarium. After mtDNA segregation, mtDNA transcription begins, which activates respiration. Mitochondria harboring wild-type genomes have functional electron transport chains and propagate more vigorously than mitochondria containing deleterious mutations in hetreoplasmic cells. Therefore, mtDNA expression acts as a stress test for the integrity of mitochondrial genomes and sets the stage for replication competition. Our observations support selective inheritance at the organelle level through a series of developmentally orchestrated mitochondrial processes. We also show that the Balbiani body has a minor role in mtDNA selective inheritance by supplying healthy mitochondria to the pole plasm. These two mechanisms may act synergistically to secure the transmission of functional mtDNA through Drosophila oogenesis.

Interpreting NUMTs in forensic genetics: Seeing the forest for the trees

Nuclear mitochondrial DNA (mtDNA) segments (NUMTs) were discovered shortly after sequencing the first human mitochondrial genome. They have earlier been considered to represent archaic elements of ancient insertion events, but modern sequencing technologies and growing databases of mtDNA and NUMT sequences confirm that they are abundant and some of them phylogenetically young. Here, we build upon mtDNA/NUMT review articles published in the mid 2010 s and focus on the distinction of NUMTs and other artefacts that can be observed in aligned sequence reads, such as mixtures (contamination), point heteroplasmy, sequencing error and cytosine deamination. We show practical examples of the effect of the mtDNA enrichment method on the representation of NUMTs in the mapped sequence data and discuss methods to bioinformatically filter NUMTs from mtDNA reads.

Associations Between Maternal Lifetime Stress and Placental Mitochondrial DNA Mutations in an Urban Multiethnic Cohort

BACKGROUND

Disrupted placental functioning due to stress can have lifelong implications. Cumulative stress and trauma are likely to have lasting impacts on maternal physiological functioning and offspring development, resulting in increased risk for later-life complex disorders for which racial disparities exist.

METHODS

This study examined the association between maternal lifetime stress and placental mitochondrial DNA mutational load in an urban multiethnic cohort. Maternal lifetime exposure to stressful events was assessed using the validated Life Stressor Checklist-Revised. Whole mitochondrial DNA sequencing was performed and mutations were determined for 365 placenta samples with complete exposure and covariate data. Multivariable regression was used to model maternal lifetime stress in relation to placental mitochondrial DNA mutational load. Racial/ethnic differences were examined by cross-product terms and contrast statements. Gene-wise analyses were conducted.

RESULTS

We identified 13,189 heteroplasmies (Phred score > 10,000, minor allele frequency < 0.5, number of mutant reads > 1). Women experiencing increased psychosocial stress over their lifetime exhibited a higher number of total placental mitochondrial mutations (β = .23, 95% confidence interval = .03 to .42) and heteroplasmic mutations (β = .18, 95% confidence interval = .05 to .31) but not homoplasmic mutations (β = -.008, 95% confidence interval = -.03 to .01); the strongest associations were observed among Black women and genes coding for NADH dehydrogenase and cytochrome c oxidase subunits.

CONCLUSIONS

Cumulative maternal lifetime stress is associated with a greater mitochondrial mutational load, particularly among Black women. The impact of racial/ethnic differences in mutational load on placental function directly affecting offspring development and/or leading to chronic disease disparities warrants further investigation.

Mitochondrial 13513G>A Mutation With Low Mutant Load Presenting as Isolated Leber's Hereditary Optic Neuropathy Assessed by Next Generation Sequencing

Objective: Mitochondrial 13513G>A mutation presenting as isolated Leber's hereditary optic neuropathy (LHON) without any extraocular pathology has not been reported in literature. We herein evaluate the clinical characteristics and heteroplasmy of m.13513G>A mutation manifesting as isolated LHON. Methods: Seven members of a Chinese family were enrolled in this study. All subjects underwent detailed systemic and ophthalmic examinations. Mitochondrial DNA in their blood was assessed by targeted PCR amplifications, next generation sequencing (NGS), and pyrosequencing. One hundred of blood samples from ethnic-matched healthy volunteers were tested by NGS and pyrosequencing as normal controls. Results: Isolated LHON without any other ocular or extraocular pathology was identified in a 16 year old patient in this family. Heteroplasmic m.13513G>A mutation was detected by NGS of the full mtDNA genome in the patient with mutant load of 33.56%, and of 26% 3 months and 3 years after the onset of LHON, respectively. No m.13513G>A mutation was detected in all his relatives by NGS. Pyrosequencing revealed the mutant load of m.13513G>A mutation of the LHON patient, his mother, father and sister were 22.4, 1.9, 0, and 0%, respectively. None of 100 healthy control subjects was detected to harbor m.13513G>A mutation either by NGS or by pyrosequencing of the full mt DNA genome. Conclusions: We first report m.13513G>A mutation with low mutant load presenting as isolated LHON. NGS of the full mitochondrial DNA genome is highly recommended for LHON suspects when targeted PCR amplification for main primary point mutations of LHON was negative.

Human molecular evolutionary rate, time dependency and transient polymorphism effects viewed through ancient and modern mitochondrial DNA genomes

Human evolutionary genetics gives a chronological framework to interpret the human history. It is based on the molecular clock hypothesis that suppose a straightforward relationship between the mutation rate and the substitution rate with independence of other factors as demography dynamics. Analyzing ancient and modern human complete mitochondrial genomes we show here that, along the time, the substitution rate can be significantly slower or faster than the average germline mutation rate confirming a time dependence effect mainly attributable to changes in the effective population size of the human populations, with an exponential growth in recent times. We also detect that transient polymorphisms play a slowdown role in the evolutionary rate deduced from haplogroup intraspecific trees. Finally, we propose the use of the most divergent lineages within haplogroups as a practical approach to correct these molecular clock mismatches.

Detection of high heteroplasmy in complete loggerhead and hawksbill sea turtles mitochondrial genomes using RNAseq

Sea turtle populations around the world face rapid decline due to the effect of anthropogenic and environmental factors. Among the affected populations are those of hawksbill turtles (Eretmochelys imbricata) and loggerhead turtles (Caretta caretta), which is why a greater effort is currently being made in their monitoring and tracing. The intragenic degree of heteroplasmic mutations, commonly associated with diseases of variable symptoms, has not been analyzed in these species. In this study, heteroplasmy in the complete mitogenome (mtDNA) of three loggerhead turtles and one hawksbill turtle was identified from data obtained by RNAseq. Individuals Cc3, Ei1, Cc1 and Cc2 presented 0.3, 1.7, 1.8 and 7.1% of heteroplasmic mutations in all their mtDNA, respectively. The protein-coding genes that presented the highest percentage of heteroplasmy were ND4 and ND5 in individual Cc2 with 16 and 38.6%, respectively. Of the tRNA genes, only tRNA^Tyr was heteroplasmic in the four individuals with 5.63% (Cc1), 25.35% (Ei1 and Cc2) and 49.3% (Cc3). In this study, we identified the critical sites of heteroplasmy in each individual and the genetic variability of their mitogenomes. The data obtained represents the baseline for future projects that evaluate the population status of these species.

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.

Contribution of Mitochondrial DNA Heteroplasmy to the Congenital Cardiac and Palatal Phenotypic Variability in Maternally Transmitted 22q11.2 Deletion Syndrome

Congenital heart disease (CHD) and palatal anomalies (PA), are among the most common characteristics of 22q11.2 deletion syndrome (22q11.2DS), but they show incomplete penetrance, suggesting the presence of additional factors. The 22q11.2 deleted region contains nuclear encoded mitochondrial genes, and since mitochondrial function is critical during development, we hypothesized that changes in the mitochondrial DNA (mtDNA) could be involved in the intrafamilial variability of CHD and PA in cases of maternally inherited 22q11.2DS. To investigate this, we studied the transmission of heteroplasmic mtDNA alleles in seventeen phenotypically concordant and discordant mother-offspring 22q11.2DS pairs. We sequenced their mtDNA and identified 26 heteroplasmic variants at >1% frequency, representing 18 transmissions. The median allele frequency change between a mother and her child was twice as much, with a wider distribution range, in PA discordant pairs, p-value = 0.039 (permutation test, 11 concordant vs. 7 discordant variants), but not in CHD discordant pairs, p-value = 0.441 (9 vs. 9). Only the variant m.9507T>C was considered to be pathogenic, but it was unrelated to the structural phenotypes. Our study is novel, yet our results are not consistent with mtDNA variation contributing to PA or CHD in 22q11.2DS. Larger cohorts and additional factors should be considered moving forward.

Pressure-Driven Mitochondrial Transfer Pipeline Generates Mammalian Cells of Desired Genetic Combinations and Fates

Generating mammalian cells with desired mitochondrial DNA (mtDNA) sequences is enabling for studies of mitochondria, disease modeling, and potential regenerative therapies. MitoPunch, a high-throughput mitochondrial transfer device, produces cells with specific mtDNA-nuclear DNA (nDNA) combinations by transferring isolated mitochondria from mouse or human cells into primary or immortal mtDNA-deficient (ρ0) cells. Stable isolated mitochondrial recipient (SIMR) cells isolated in restrictive media permanently retain donor mtDNA and reacquire respiration. However, SIMR fibroblasts maintain a ρ0-like cell metabolome and transcriptome despite growth in restrictive media. We reprogrammed non-immortal SIMR fibroblasts into induced pluripotent stem cells (iPSCs) with subsequent differentiation into diverse functional cell types, including mesenchymal stem cells (MSCs), adipocytes, osteoblasts, and chondrocytes. Remarkably, after reprogramming and differentiation, SIMR fibroblasts molecularly and phenotypically resemble unmanipulated control fibroblasts carried through the same protocol. Thus, our MitoPunch "pipeline" enables the production of SIMR cells with unique mtDNA-nDNA combinations for additional studies and applications in multiple cell types.

Lighting up single-nucleotide variation in situ in single cells and tissues

The ability to 'see' genetic information directly in single cells can provide invaluable insights into complex biological systems. In this review, we discuss recent advances of in situ imaging technologies for visualizing the subtlest sequence alteration, single-nucleotide variation (SNV), at single-cell level. The mechanism of recently developed methods for SNV discrimination are summarized in detail. With recent developments, single-cell SNV imaging methods have opened a new door for studying the heterogenous and stochastic genetic information in individual cells. Furthermore, SNV imaging can be used on morphologically preserved tissue, which can provide information on histological context for gene expression profiling in basic research and genetic diagnosis. Moreover, the ability to visualize SNVs in situ can be further developed into in situ sequencing technology. We expect this review to inspire more research work into in situ SNV imaging technologies for investigating cellular phenotypes and gene regulation at single-nucleotide resolution, and developing new clinical and biomedical applications.

STAMP: a multiplex sequencing method for simultaneous evaluation of mitochondrial DNA heteroplasmies and content

Human mitochondrial genome (mtDNA) variations, such as mtDNA heteroplasmies (the co-existence of mutated and wild-type mtDNA), have received increasing attention in recent years for their clinical relevance to numerous diseases. But large-scale population studies of mtDNA heteroplasmies have been lagging due to the lack of a labor- and cost-effective method. Here, we present a novel human mtDNA sequencing method called STAMP (sequencing by targeted amplification of multiplex probes) for measuring mtDNA heteroplasmies and content in a streamlined workflow. We show that STAMP has high-mapping rates to mtDNA, deep coverage of unique reads and high tolerance to sequencing and polymerase chain reaction errors when applied to human samples. STAMP also has high sensitivity and low false positive rates in identifying artificial mtDNA variants at fractions as low as 0.5% in genomic DNA samples. We further extend STAMP, by including nuclear DNA-targeting probes, to enable assessment of relative mtDNA content in the same assay. The high cost-effectiveness of STAMP, along with the flexibility of using it for measuring various aspects of mtDNA variations, will accelerate the research of mtDNA heteroplasmies and content in large population cohorts, and in the context of human diseases and aging.

Mitonuclear Interactions in the Maintenance of Mitochondrial Integrity

In eukaryotic cells, mitochondria originated in an α-proteobacterial endosymbiont. Although these organelles harbor their own genome, the large majority of genes, originally encoded in the endosymbiont, were either lost or transferred to the nucleus. As a consequence, mitochondria have become semi-autonomous and most of their processes require the import of nuclear-encoded components to be functional. Therefore, the mitochondrial-specific translation has evolved to be coordinated by mitonuclear interactions to respond to the energetic demands of the cell, acquiring unique and mosaic features. However, mitochondrial-DNA-encoded genes are essential for the assembly of the respiratory chain complexes. Impaired mitochondrial function due to oxidative damage and mutations has been associated with numerous human pathologies, the aging process, and cancer. In this review, we highlight the unique features of mitochondrial protein synthesis and provide a comprehensive insight into the mitonuclear crosstalk and its co-evolution, as well as the vulnerabilities of the animal mitochondrial genome.

Stable retention of chloramphenicol-resistant mtDNA to rescue metabolically impaired cells

The permanent transfer of specific mtDNA sequences into mammalian cells could generate improved models of mtDNA disease and support future cell-based therapies. Previous studies documented multiple biochemical changes in recipient cells shortly after mtDNA transfer, but the long-term retention and function of transferred mtDNA remains unknown. Here, we evaluate mtDNA retention in new host cells using ‘MitoPunch’, a device that transfers isolated mitochondria into mouse and human cells. We show that newly introduced mtDNA is stably retained in mtDNA-deficient (ρ0) recipient cells following uridine-free selection, although exogenous mtDNA is lost from metabolically impaired, mtDNA-intact (ρ+) cells. We then introduced a second selective pressure by transferring chloramphenicol-resistant mitochondria into chloramphenicol-sensitive, metabolically impaired ρ+ mouse cybrid cells. Following double selection, recipient cells with mismatched nuclear (nDNA) and mitochondrial (mtDNA) genomes retained transferred mtDNA, which replaced the endogenous mutant mtDNA and improved cell respiration. However, recipient cells with matched mtDNA-nDNA failed to retain transferred mtDNA and sustained impaired respiration. Our results suggest that exogenous mtDNA retention in metabolically impaired ρ+ recipients depends on the degree of recipient mtDNA-nDNA co-evolution. Uncovering factors that stabilize exogenous mtDNA integration will improve our understanding of in vivo mitochondrial transfer and the interplay between mitochondrial and nuclear genomes.

A Forensic Genomics Approach for the Identification of Sister Marija Crucifiksa Kozulić

Sister Marija Krucifiksa Kozulić (1852–1922) was a Croatian nun who is in consideration for beatification by the Vatican, which is facilitated by the identification of her 20th-century remains. Sister Marija was buried in a tomb in Rijeka, Croatia, along with other nuns including her biological sister, Tereza Kozulić (1861–1933). When the remains were exhumed in 2011, they were found in a deteriorated state and commingled with several other sets of remains. Thus, mitochondrial genome sequencing of the long bones was performed to sort the remains by mitochondrial haplotype. Two similar but unique haplotypes belonging to haplogroup H1bu were identified, and samples from these bones were subjected to autosomal short tandem repeat (STR) and single nucleotide polymorphism (SNP) sequencing. Although only partial profiles were obtained, the data were sufficient for kinship analysis with the profile of a paternal niece of Sister Marija (Fides Kozulić). The data indicate that it is 574,195-fold more likely that the two sets of skeletal remains represent 2nd-degree relatives of Fides than sisters who are unrelated to Fides. Although it is impossible to discern which set of remains belongs to Marija and which belongs to Tereza, forensic genomics methods have enabled identification of the sisters.

GPR174 and ITM2A Gene Polymorphisms rs3827440 and rs5912838 on the X chromosome in Korean Children with Autoimmune Thyroid Disease

(1) Background: Autoimmune thyroid diseases (AITDs) are female predominant and much attention has been focused on G protein-coupled receptor 174 (GPR174) and integral membrane protein 2A (ITM2A) on the X chromosome as Grave's disease (GD) susceptible locus. (2) Methods: We genotyped four single nucleotide polymorphisms (SNPs), rs3810712, rs3810711, rs3827440, and rs5912838, of GPR174 and ITM2A in 115 Korean children with AITD (M = 25 and F = 90; GD = 74 (14.7 ± 3.6 years), HD = 41 (13.4 ± 3.2 years); GD-thyroid-associated ophthalmopathy (TAO) = 40, GD-non-TAO=34) and 204 healthy Korean individuals (M = 104 and F = 100). The data were analyzed by sex-stratified or combined. (3) Results: Three SNPs, rs3810712, rs3810711 and rs3827440, were found to be in perfect linkage disequilibrium (D' = 1, r² = 1). In AITD, HD, GD, GD-TAO, and GD-non-TAO patients, rs3827440 TT/T and rs5912838 AA/A were susceptible and rs3827440 CC/C and rs5912838 CC/C were protective genotypes. When analyzed by sex, rs3827440 TT and rs5912838 AA were susceptible and rs3827440 CC and rs5912838 CC were protective genotypes in female AITD, GD, GD-TAO, and GD-non-TAO subjects. In male AITD patients, rs3827440 T and rs5912838 A were susceptible and rs3827440 C and rs5912838 C were protective genotypes. (4) Conclusions: Polymorphisms in GPR174 and ITM2A genes on the X chromosome might be associated with AITD in Korean children.

Age-related accumulation of de novo mitochondrial mutations in mammalian oocytes and somatic tissues

Mutations create genetic variation for other evolutionary forces to operate on and cause numerous genetic diseases. Nevertheless, how de novo mutations arise remains poorly understood. Progress in the area is hindered by the fact that error rates of conventional sequencing technologies (1 in 100 or 1,000 base pairs) are several orders of magnitude higher than de novo mutation rates (1 in 10,000,000 or 100,000,000 base pairs per generation). Moreover, previous analyses of germline de novo mutations examined pedigrees (and not germ cells) and thus were likely affected by selection. Here, we applied highly accurate duplex sequencing to detect low-frequency, de novo mutations in mitochondrial DNA (mtDNA) directly from oocytes and from somatic tissues (brain and muscle) of 36 mice from two independent pedigrees. We found mtDNA mutation frequencies 2- to 3-fold higher in 10-month-old than in 1-month-old mice, demonstrating mutation accumulation during the period of only 9 mo. Mutation frequencies and patterns differed between germline and somatic tissues and among mtDNA regions, suggestive of distinct mutagenesis mechanisms. Additionally, we discovered a more pronounced genetic drift of mitochondrial genetic variants in the germline of older versus younger mice, arguing for mtDNA turnover during oocyte meiotic arrest. Our study deciphered for the first time the intricacies of germline de novo mutagenesis using duplex sequencing directly in oocytes, which provided unprecedented resolution and minimized selection effects present in pedigree studies. Moreover, our work provides important information about the origins and accumulation of mutations with aging/maturation and has implications for delayed reproduction in modern human societies. Furthermore, the duplex sequencing method we optimized for single cells opens avenues for investigating low-frequency mutations in other studies.

Counterbalancing the time-dependent effect on the human mitochondrial DNA molecular clock

BACKGROUND

The molecular clock is an important genetic tool for estimating evolutionary timescales. However, the detection of a time-dependent effect on substitution rate estimates complicates its application. It has been suggested that demographic processes could be the main cause of this confounding effect. In the present study, I propose a new algorithm for estimating the coalescent age of phylogenetically related sequences, taking into account the observed time-dependent effect on the molecular rate detected by others.

RESULTS

By applying this method to real human mitochondrial DNA trees with shallow and deep topologies, I obtained significantly older molecular ages for the main events of human evolution than were previously estimated. These ages are in close agreement with the most recent archaeological and paleontological records favoring the emergence of early anatomically modern humans in Africa 315 ± 34 thousand years ago (kya) and the presence of recent modern humans outside of Africa as early as 174 ± 48 thousand years ago. Furthermore, during the implementation process, I demonstrated that in a population with fluctuating sizes, the probability of fixation of a new neutral mutant depends on the effective population size, which is in better accordance with the fact that under the neutral theory of molecular evolution, the fate of a molecular mutation is mainly determined by random drift.

CONCLUSIONS

I suggest that the demographic history of populations has a more decisive effect than purifying selection and/or mutational saturation on the time-dependent effect observed for the substitution rate, and I propose a new method that corrects for this effect.

One Hundred Mitochondrial Genomes of Cicadas

Mitochondrial genomes can provide valuable information on the biology and evolutionary histories of their host organisms. Here, we present and characterize the complete coding regions of 107 mitochondrial genomes (mitogenomes) of cicadas (Insecta: Hemiptera: Auchenorrhyncha: Cicadoidea), representing 31 genera, 61 species, and 83 populations. We show that all cicada mitogenomes retain the organization and gene contents thought to be ancestral in insects, with some variability among cicada clades in the length of a region between the genes nad2 and cox1, which encodes 3 tRNAs. Phylogenetic analyses using these mitogenomes recapitulate a recent 5-gene classification of cicadas into families and subfamilies, but also identify a species that falls outside of the established taxonomic framework. While protein-coding genes are under strong purifying selection, tests of relative evolutionary rates reveal significant variation in evolutionary rates across taxa, highlighting the dynamic nature of mitochondrial genome evolution in cicadas. These data will serve as a useful reference for future research into the systematics, ecology, and evolution of the superfamily Cicadoidea.

Equine Aging and the Oocyte: A Potential Model for Reproductive Aging in Women

Numerous similarities in reproductive aging have been documented between the mare and woman. Aging is associated with a decline in fertility. In mares and women, oocyte transfer procedures were initially used to establish that oocyte donor age is associated with oocyte quality. Age-associated differences in oocytes include altered morphology, gene expression, and developmental potential. Reactive oxygen species and mitochondrial dysfunction are thought to be important contributors to loss of oocyte quality. In the woman, aneuploidy is a primary consideration with maternal aging. Although misalignment of chromosomes during meiosis has been observed in the mare, less is known in this area. Reproductive aging will be reviewed in the mare and compared with the woman with emphasis on factors that affect oocyte quality and developmental potential. Areas in which the mare could be used as a research model to study reproductive aging in women will be highlighted.

Evidence of Neutral Evolution of Mitochondrial DNA in Human Hepatocellular Carcinoma

Many studies have suggested that mitochondria and mitochondrial DNA (mtDNA) might be functionally associated with tumor genesis and development. Although the heterogeneity of tumors is well known, most studies were based on the analysis of a single tumor sample. The extent of mtDNA diversity in the same tumor is unclear, as is whether the diversity is influenced by selection pressure. Here, we analyzed the whole exon data from 1 nontumor sample and 23 tumor samples from different locations of one single tumor tissue from a hepatocellular carcinoma (HCC) patient. Among 18 heteroplasmic sites identified in the tumor, only 2 heteroplasmies were shared among all tumor samples. By investigating the correlations between the occurrence and frequency of heteroplasmy (Het) and sampling locations (Coordinate), relative mitochondrial copy numbers, and single-nucleotide variants in the nuclear genome, we found that the Coordinate was significantly correlated with Het, suggesting no strong purifying selection or positive selection acted on the mtDNA in HCC. By further investigating the allele frequency and proportion of nonsynonymous mutations in the tumor mtDNA, we found that mtDNA in HCC did not undergo extra selection compared with mtDNA in the adjacent nontumor tissue, and they both likely evolved under neutral selection.

A High-Resolution View of Adaptive Event Dynamics in a Plasmid

Coadaptation between bacterial hosts and plasmids frequently results in adaptive changes restricted exclusively to host genome leaving plasmids unchanged. To better understand this remarkable stability, we transformed naïve Escherichia coli cells with a plasmid carrying an antibiotic-resistance gene and forced them to adapt in a turbidostat environment. We then drew population samples at regular intervals and subjected them to duplex sequencing—a technique specifically designed for identification of low-frequency mutations. Variants at ten sites implicated in plasmid copy number control emerged almost immediately, tracked consistently across the experiment’s time points, and faded below detectable frequencies toward the end. This variation crash coincided with the emergence of mutations on the host chromosome. Mathematical modeling of trajectories for adaptive changes affecting plasmid copy number showed that such mutations cannot readily fix or even reach appreciable frequencies. We conclude that there is a strong selection against alterations of copy number even if it can provide a degree of growth advantage. This incentive is likely rooted in the complex interplay between mutated and wild-type plasmids constrained within a single cell and underscores the importance of understanding of intracellular plasmid variability.

A battle for transmission: the cooperative and selfish animal mitochondrial genomes

The mitochondrial genome is an evolutionarily persistent and cooperative component of metazoan cells that contributes to energy production and many other cellular processes. Despite sharing the same host as the nuclear genome, the multi-copy mitochondrial DNA (mtDNA) follows very different rules of replication and transmission, which translate into differences in the patterns of selection. On one hand, mtDNA is dependent on the host for its transmission, so selections would favour genomes that boost organismal fitness. On the other hand, genetic heterogeneity within an individual allows different mitochondrial genomes to compete for transmission. This intra-organismal competition could select for the best replicator, which does not necessarily give the fittest organisms, resulting in mito-nuclear conflict. In this review, we discuss the recent advances in our understanding of the mechanisms and opposing forces governing mtDNA transmission and selection in bilaterians, and what the implications of these are for mtDNA evolution and mitochondrial replacement therapy.

Deleterious mtDNA mutations are common in mature oocytes

Heritable mitochondrial DNA (mtDNA) mutations are common, yet only a few recurring pathogenic mtDNA variants account for the majority of known familial cases in humans. Purifying selection in the female germline is thought to be responsible for the elimination of most harmful mtDNA mutations during oogenesis. Here we show that deleterious mtDNA mutations are abundant in ovulated mature mouse oocytes and preimplantation embryos recovered from PolG mutator females but not in their live offspring. This implies that purifying selection acts not in the maternal germline per se, but during post-implantation development. We further show that oocyte mtDNA mutations can be captured and stably maintained in embryonic stem cells and then reintroduced into chimeras, thereby allowing examination of the effects of specific mutations on fetal and postnatal development.

Family reunion via error correction: an efficient analysis of duplex sequencing data

BACKGROUND

Duplex sequencing is the most accurate approach for identification of sequence variants present at very low frequencies. Its power comes from pooling together multiple descendants of both strands of original DNA molecules, which allows distinguishing true nucleotide substitutions from PCR amplification and sequencing artifacts. This strategy comes at a cost-sequencing the same molecule multiple times increases dynamic range but significantly diminishes coverage, making whole genome duplex sequencing prohibitively expensive. Furthermore, every duplex experiment produces a substantial proportion of singleton reads that cannot be used in the analysis and are thrown away.

RESULTS

In this paper we demonstrate that a significant fraction of these reads contains PCR or sequencing errors within duplex tags. Correction of such errors allows "reuniting" these reads with their respective families increasing the output of the method and making it more cost effective.

CONCLUSIONS

We combine an error correction strategy with a number of algorithmic improvements in a new version of the duplex analysis software, Du Novo 2.0. It is written in Python, C, AWK, and Bash. It is open source and readily available through Galaxy, Bioconda, and Github: https://github.com/galaxyproject/dunovo.

Somatic mitochondrial mutation discovery using ultra-deep sequencing of the mitochondrial genome reveals spatial tumor heterogeneity in head and neck squamous cell carcinoma

Mutations in mitochondrial DNA (mtDNA) have been linked to risk, progression, and treatment response of head and neck squamous cell carcinoma (HNSCC). Due to their clonal nature and high copy number, mitochondrial mutations could serve as powerful molecular markers for detection of cancer cells in bodily fluids, surgical margins, biopsies and lymph node (LN) metastasis, especially at sites where tumor involvement is not histologically apparent. Despite a pressing need for high-throughput, cost-effective mtDNA mutation profiling system, current methods for library preparation are still imperfect for detection of low prevalence heteroplasmic mutations. To this end, we have designed an ultra-deep amplicon-based sequencing library preparation approach that covers the entire mitochondrial genome. We sequenced mtDNA in 28 HNSCCs, matched LNs, surgical margins and bodily fluids, and applied multiregional sequencing approach on 14 primary tumors. Our results demonstrate that this quick, sensitive and cost-efficient method allows obtaining a snapshot on the mitochondrial heterogeneity, and can be used for detection of low frequency tumor-associated mtDNA mutations in LNs, sputum and serum specimens. These findings provide the foundation for using mitochondrial sequencing for risk assessment, early detection, and tumor surveillance.

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.

Characterization of background noise in MiSeq MPS data when sequencing human mitochondrial DNA from various sample sources and library preparation methods

Improved resolution of massively parallel sequencing (MPS) allows for the characterization of mitochondrial (mt) DNA heteroplasmy to levels previously unattainable with traditional sequencing approaches. An essential criterion for the reporting of heteroplasmy is the ability of the MPS method to distinguish minor sequence variants (MSVs) from system noise, or error. Therefore, an assessment of the background noise in the MPS method is desirable to identify the point at which reliable data can be reported. Substitution and sequence specific error (SSE) was evaluated for a variety of sample types and two library preparations. Substitution error rates ranged from 0.18 to 0.49 per 100 nucleotides with C positions generally having the highest rate of misincorporation. Comparison of error rates across sample types indicated a significant increase for samples with damaged DNA. The positions of error were varied across datasets (pairwise concordance 0-68%), but had greater consistency within the damaged samples (80-96%). The most commonly observed motif preceding error in forward reads was CCG, while GGT was most common in reverse reads, both consistent with previous findings. The findings illustrate that for datasets containing samples with damaged DNA, reporting thresholds for heteroplasmy may have to be modified and individual sites with error levels exceeding thresholds should be scrutinized. Collectively, the shifting error profiles observed across the various sample types and library preparation methods demonstrates the need for an assessment of error under these varying circumstances. Characterization of the applicable background noise will help to ensure that thresholds are reliably set for detection of true MSVs.

Molecular basis of Leigh syndrome: a current look

Leigh Syndrome (OMIM 256000) is a heterogeneous neurologic disorder due to damage in mitochondrial energy production that usually starts in early childhood. The first description given by Leigh pointed out neurological symptoms in children under 2 years and premature death. Following cases brought some hypothesis to explain the cause due to similarity to other neurological diseases and led to further investigation for metabolic diseases. Biochemical evaluation and specific metabolic profile suggested impairment in energy production (OXPHOS) in mitochondria. As direct approach to involved tissues is not always possible or safe, molecular analysis is a great cost-effective option and, besides biochemical results, is required to confirm the underlying cause of this syndrome face to clinical suspicion. The Next Generation Sequencing (NGS) advance represented a breakthrough in molecular biology allowing simultaneous gene analysis giving short-time results and increasing the variants underlying this syndrome, counting over 75 monogenic causes related so far. NGS provided confirmation of emerging cases and brought up diagnosis in atypical presentations as late-onset cases, which turned Leigh into a heterogeneous syndrome with variable outcomes. This review highlights clinical presentation in both classic and atypical phenotypes, the investigation pathway throughout confirmation emphasizing the underlying genetic heterogeneity and increasing number of genes assigned to this syndrome as well as available treatment.

Population-specific signatures of intra-individual mitochondrial DNA heteroplasmy and their potential evolutionary advantages

Heteroplasmy is the existence of more than one mitochondrial DNA (mtDNA) variant within a cell. The evolutionary mechanisms of heteroplasmy are not fully understood, despite being a very common phenomenon. Here we combined heteroplasmy measurements using high throughput sequencing on green turtles (Chelonia mydas) with simulations to understand how heteroplasmy modulates population diversity across generations and under different demographic scenarios. We found heteroplasmy to be widespread in all individuals analysed, with consistent signal in individuals across time and tissue. Significant shifts in haplotype composition were found from mother to offspring, signalling the effect of the cellular bottleneck during oogenesis as included in the model. Our model of mtDNA inheritance indicated that heteroplasmy favoured the increase of population diversity through time and buffered against population bottlenecks, thus indicating the importance of this phenomenon in species with reduced population sizes and frequent population bottlenecks like marine turtles. Individuals with recent haplotypes showed higher levels of heteroplasmy than the individuals with ancient haplotypes, suggesting a potential advantage of maintaining established copies when new mutations arise. We recommend using heteroplasmy through high throughput sequencing in marine turtles, as well as other wildlife populations, for diversity assessment, population genetics, and mixed stock analysis.

Pronounced somatic bottleneck in mitochondrial DNA of human hair

Heteroplasmy is the presence of variable mitochondrial DNA (mtDNA) within the same individual. The dynamics of heteroplasmy allele frequency among tissues of the human body is not well understood. Here, we measured allele frequency at heteroplasmic sites in two to eight hairs from each of 11 humans using next-generation sequencing. We observed a high variance in heteroplasmic allele frequency among separate hairs from the same individual—much higher than that for blood and cheek tissues. Our population genetic modelling estimated the somatic bottleneck during embryonic follicle development of separate hairs to be only 11.06 (95% confidence interval 0.6–34.0) mtDNA segregating units. This bottleneck is much more drastic than somatic bottlenecks for blood and cheek tissues (136 and 458 units, respectively), as well as more drastic than, or comparable to, the germline bottleneck (equal to 25–32 or 7–10 units, depending on the study). We demonstrated that hair undergoes additional genetic drift before and after the divergence of mtDNA lineages of individual hair follicles. Additionally, we showed a positive correlation between donor's age and variance in heteroplasmy allele frequency in hair. These findings have important implications for forensics and for our understanding of mtDNA dynamics in the human body.

This article is part of the theme issue ‘Linking the mitochondrial genotype to phenotype: a complex endeavour’.

Portunus pelagicus mtDNA heteroplasmy inheritance and its effect on the use of mtCR and mtCOI sequence data

Mitochondrial DNA (mtDNA) sequences, especially mitochondrial control region (mtCR) and mitochondrial cytochrome c oxidase subunit I (mtCOI), have been widely used in population and evolutionary genetic analyses of metazoan. The presence of mtDNA heteroplasmy - a mixture of mtDNA haplotypes - possibly affects these analyses. This study aimed to reveal mtDNA heteroplasmy in mtCR, mtCOI, and mtND2 (mitochondrial NADH dehydrogenase subunit 2) of Portunus pelagicus, and examine its effect on the use of mtCR and mtCOI sequences. The screening result showed that the probability of observing mtDNA heteroplasmy was approximately 8%. Across the three targeted regions, 92 heteroplasmic variants were observed from seven samples comprising three mothers and four offspring. Most inherited heteroplasmy presented transition and silence mutation. By comparing the proportion of shared variants among maternal relatives to that among non-relatives, the result suggested that most heteroplasmic variants observed in an individual are inherited. Statistical analyses carried out on the inter-generational differences suggested that random drift and purifying selection play roles in determining the offspring's heteroplasmy level. The size of the random shift varies according to the location of variants and the mothers. The phylogenetic analysis showed that the presence of mtDNA heteroplasmy in mtCR and mtCOI does not affect familial and species identification, respectively. This study firstly reported the mtDNA heteroplasmy in P. pelagicus, its inheritance pattern, and its effect on the use of mtDNA sequence data. This basic knowledge would be useful for the study based on mtDNA sequence data, especially in other invertebrates.

Bottleneck and selection in the germline and maternal age influence transmission of mitochondrial DNA in human pedigrees

Heteroplasmy-the presence of multiple mitochondrial DNA (mtDNA) haplotypes in an individual-can lead to numerous mitochondrial diseases. The presentation of such diseases depends on the frequency of the heteroplasmic variant in tissues, which, in turn, depends on the dynamics of mtDNA transmissions during germline and somatic development. Thus, understanding and predicting these dynamics between generations and within individuals is medically relevant. Here, we study patterns of heteroplasmy in 2 tissues from each of 345 humans in 96 multigenerational families, each with, at least, 2 siblings (a total of 249 mother-child transmissions). This experimental design has allowed us to estimate the timing of mtDNA mutations, drift, and selection with unprecedented precision. Our results are remarkably concordant between 2 complementary population-genetic approaches. We find evidence for a severe germline bottleneck (7-10 mtDNA segregating units) that occurs independently in different oocyte lineages from the same mother, while somatic bottlenecks are less severe. We demonstrate that divergence between mother and offspring increases with the mother's age at childbirth, likely due to continued drift of heteroplasmy frequencies in oocytes under meiotic arrest. We show that this period is also accompanied by mutation accumulation leading to more de novo mutations in children born to older mothers. We show that heteroplasmic variants at intermediate frequencies can segregate for many generations in the human population, despite the strong germline bottleneck. We show that selection acts during germline development to keep the frequency of putatively deleterious variants from rising. Our findings have important applications for clinical genetics and genetic counseling.

Quantification of mitochondrial DNA from peripheral tissues: Limitations in predicting the severity of neurometabolic disorders and proposal of a novel diagnostic test

Neurometabolic disorders stem from errors in metabolic processes yielding a neurological phenotype. A subset of those disorders encompasses mitochondrial abnormalities partially due to mitochondrial DNA (mtDNA) depletion. mtDNA depletion can be attributed to inheritance, spontaneous mutations or acquired from drug-related toxicities. In the armamentarium of diagnostic procedures, mtDNA quantification is a standard for disease classification. However, alterations in mtDNA obtained from peripheral tissues such as skin fibroblasts and blood cells do not often reflect the severity of the affected organ, in this case, the brain. The purpose of this review is to highlight the pitfalls of quantitating mtDNA from peripheral -and not limited to-tissues for diagnosing patients suffering from a variety of mtDNA depletion syndromes exhibiting neurologic abnormalities. In lieu, a qualitative test of mitochondrial substrate-level phosphorylation -even from peripheral tissues-reflecting the ability of mitochondria to rely on glutaminolysis in the presence of respiratory chain defects is proposed as a novel diagnostic assessment of mitochondrial functionality.

Varied Mechanisms and Models for the Varying Mitochondrial Bottleneck

Mitochondrial DNA (mtDNA) molecules exist in populations within cells, and may carry mutations. Different cells within an organism, and organisms within a family, may have different proportions of mutant mtDNA in these cellular populations. This diversity is often thought of as arising from a “genetic bottleneck.” This article surveys approaches to characterize and model the generation of this genetic diversity, aiming to provide an introduction to the range of concepts involved, and to highlight some recent advances in understanding. In particular, differences between the statistical “genetic bottleneck” (mutant proportion spread) and the physical mtDNA bottleneck and other cellular processes are highlighted. Particular attention is paid to the quantitative analysis of the “genetic bottleneck,” estimation of its magnitude from observed data, and inference of its underlying mechanisms. Evidence that the “genetic bottleneck” (mutant proportion spread) varies with age, between individuals and species, and across mtDNA sequences, is described. The interpretation issues that arise from sampling errors, selection, and different quantitative definitions are also discussed.

Heteroplasmy in the complete chicken mitochondrial genome

Chicken mitochondrial DNA is a circular molecule comprising ~16.8 kb. In this study, we used next-generation sequencing to investigate mitochondrial heteroplasmy in the whole chicken mitochondrial genome. Based on heteroplasmic detection thresholds at the 0.5% level, 178 cases of heteroplasmy were identified in the chicken mitochondrial genome, where 83% were due to nucleotide transitions. D-loop regionwas hot spot region for mtDNA heteroplasmy in the chicken since 130 cases of heteroplasmy were located in these regions. Heteroplasmy varied among intraindividual tissues with allele-specific, position-specific, and tissue-specific features. Skeletal muscle had the highest abundance of heteroplasmy. Cases of heteroplasmy at mt.G8682A and mt.G16121A were validated by PCR-restriction fragment length polymorphism analysis, which showed that both had low ratios of heteroplasmy occurrence in five natural breeds. Polymorphic sites were easy to distinguish. Based on NGS data for crureus tissues, mitochondrial mutation/heteroplasmy exhibited clear maternal inheritance features at the whole mitochondrial genomic level. Further investigations of the heterogeneity of the mt.A5694T and mt.T5718G transitions between generations using pyrosequencing based on pedigree information indicated that the degree of heteroplasmy and the occurrence ratio of heteroplasmy decreased greatly from the F0 to F1 generations in the mt.A5694T and mt.T5718G site. Thus, the intergenerational transmission of heteroplasmy in chicken mtDNA exhibited a rapid shift toward homoplasmy within a single generation. Our findings indicate that heteroplasmy is a widespread phenomenon in chicken mitochondrial genome, in which most sites exhibit low heteroplasmy and the allele frequency at heteroplasmic sites changes significantly during transmission events. It suggests that heteroplasmy may be under negative selection to some degree in the chicken.

Bayesian Estimation of Population Size Changes by Sampling Tajima's Trees

The large state space of gene genealogies is a major hurdle for inference methods based on Kingman's coalescent. Here, we present a new Bayesian approach for inferring past population sizes, which relies on a lower-resolution coalescent process that we refer to as "Tajima's coalescent." Tajima's coalescent has a drastically smaller state space, and hence it is a computationally more efficient model, than the standard Kingman coalescent. We provide a new algorithm for efficient and exact likelihood calculations for data without recombination, which exploits a directed acyclic graph and a correspondingly tailored Markov Chain Monte Carlo method. We compare the performance of our Bayesian Estimation of population size changes by Sampling Tajima's Trees (BESTT) with a popular implementation of coalescent-based inference in BEAST using simulated and human data. We empirically demonstrate that BESTT can accurately infer effective population sizes, and it further provides an efficient alternative to the Kingman's coalescent. The algorithms described here are implemented in the R package phylodyn, which is available for download at https://github.com/JuliaPalacios/phylodyn.

Validation of NGS for mitochondrial DNA casework at the FBI Laboratory

As a first step towards integrating next generation sequencing (NGS) technology into the FBI Laboratory's operational casework, the PowerSeq™ CRM Nested System, an NGS-based mitochondrial DNA (mtDNA) control region assay, was developmentally and internally validated. The validation studies were conducted in accordance with the Scientific Working Group on DNA Analysis Methods (SWGDAM) Validation Guidelines for Forensic DNA Analysis Methods, and the FBI's Quality Assurance Standards (QAS) for Forensic DNA Testing Laboratories. The assay was shown to be highly reproducible, with variant frequencies across intra and inter-run replicates of the same sample differing, on average, by just 0.3% for substitutions and point heteroplasmies and 1.5% for insertions and deletions. The assay was also shown to be extremely sensitive, yielding complete control region sequence data from as few as 2000 copies of mtDNA. This is a more than 20-fold increase in sensitivity when compared to the FBI Laboratory's current Sanger sequencing-based protocols and, based on mtDNA quantitation values of samples routinely encountered in mtDNA casework, suggests that the percentage of questioned samples from which full control region data can be recovered will increase from our current 20% to approximately 90% success with NGS technology. In addition, the assay requires on average only 30% of the extract volume typically required to develop control region profiles from degraded samples via Sanger sequencing. Overall, these studies establish the reliability of the PowerSeq™ CRM Nested System for accurate mtDNA control region typing and can serve as a model for laboratories seeking to validate NGS protocols for forensic mtDNA analysis.

Sharing of heteroplasmies between human liver lobes varies across the mtDNA genome

Mitochondrial DNA (mtDNA) heteroplasmy (intra-individual variation) varies among different human tissues and increases with age, suggesting that the majority of mtDNA heteroplasmies are acquired, rather than inherited. However, the extent to which heteroplasmic sites are shared across a tissue remains an open question. We therefore investigated heteroplasmy in two liver samples (one from each primary lobe) from 83 Europeans, sampled at autopsy. Minor allele frequencies (MAF) at heteroplasmic sites were significantly correlated between the two liver samples from an individual, with significantly more sharing of heteroplasmic sites in the control region than in the non-control region. We show that this increased sharing for the control region cannot be explained by recent mutations at just a few specific heteroplasmic sites or by the possible presence of 7S DNA. Moreover, we carried out simulations to show that there is significantly more sharing than would be predicted from random genetic drift from a common progenitor cell. We also observe a significant excess of non-synonymous vs. synonymous heteroplasmies in the protein-coding region, but significantly more sharing of synonymous heteroplasmies. These contrasting patterns for the control vs. the non-control region, and for non-synonymous vs. synonymous heteroplasmies, suggest that selection plays a role in heteroplasmy sharing.

New Frontiers in IVF: mtDNA and autologous germline mitochondrial energy transfer

Many infertility specialists support the existence of a relationship between the levels of mitochondrial DNA and the quality of the blastocysts. Despite the extensive use of pre-implantation genetic testing for aneuploidy, a significant percentage of euploid embryos do not implant even though the endometrium is normal. Mitochondrial DNA may be used as a new test in evaluating embryonic vitality.Ovarian aging leads to a decrease in the quantity and quality of oocytes and aged oocytes have a reduced number of mitochondria. Mitochondria are the energy factories of the cells and their lacked could leads to lower fertilization rates and poor embryonic development. Various strategies have been tested to increase the mitochondria quantity and thus improve the quality of oocytes used in in vitro fertilization. Results of ovarian rejuvenation techniques such as autologous mitochondrial transplantation have been controversial. In this review, we describe the state of the art concerning the use of mitochondrial DNA and autologous mitochondrial transplantation as new possibilities to increase success in vitro fertilization.

Regulation of nuclear epigenome by mitochondrial DNA heteroplasmy

Diseases associated with mitochondrial DNA (mtDNA) mutations are highly variable in phenotype, in large part because of differences in the percentage of normal and mutant mtDNAs (heteroplasmy) present within the cell. For example, increasing heteroplasmy levels of the mtDNA tRNA^Leu(UUR) nucleotide (nt) 3243A > G mutation result successively in diabetes, neuromuscular degenerative disease, and perinatal lethality. These phenotypes are associated with differences in mitochondrial function and nuclear DNA (nDNA) gene expression, which are recapitulated in cybrid cell lines with different percentages of m.3243G mutant mtDNAs. Using metabolic tracing, histone mass spectrometry, and NADH fluorescence lifetime imaging microscopy in these cells, we now show that increasing levels of this single mtDNA mutation cause profound changes in the nuclear epigenome. At high heteroplasmy, mitochondrially derived acetyl-CoA levels decrease causing decreased histone H4 acetylation, with glutamine-derived acetyl-CoA compensating when glucose-derived acetyl-CoA is limiting. In contrast, α-ketoglutarate levels increase at midlevel heteroplasmy and are inversely correlated with histone H3 methylation. Inhibition of mitochondrial protein synthesis induces acetylation and methylation changes, and restoration of mitochondrial function reverses these effects. mtDNA heteroplasmy also affects mitochondrial NAD⁺/NADH ratio, which correlates with nuclear histone acetylation, whereas nuclear NAD⁺/NADH ratio correlates with changes in nDNA and mtDNA transcription. Thus, mutations in the mtDNA cause distinct metabolic and epigenomic changes at different heteroplasmy levels, potentially explaining transcriptional and phenotypic variability of mitochondrial disease.

Germline selection shapes human mitochondrial DNA diversity

Approximately 2.4% of the human mitochondrial DNA (mtDNA) genome exhibits common homoplasmic genetic variation. We analyzed 12,975 whole-genome sequences to show that 45.1% of individuals from 1526 mother-offspring pairs harbor a mixed population of mtDNA (heteroplasmy), but the propensity for maternal transmission differs across the mitochondrial genome. Over one generation, we observed selection both for and against variants in specific genomic regions; known variants were more likely to be transmitted than previously unknown variants. However, new heteroplasmies were more likely to match the nuclear genetic ancestry as opposed to the ancestry of the mitochondrial genome on which the mutations occurred, validating our findings in 40,325 individuals. Thus, human mtDNA at the population level is shaped by selective forces within the female germ line under nuclear genetic control, which ensures consistency between the two independent genetic lineages.

Mutation signatures in germline mitochondrial genome provide insights into human mitochondrial evolution and disease

Variations in mitochondrial DNA (mtDNA) have been fundamental for understanding human evolution and are causative for a plethora of inherited mitochondrial diseases, but the mutation signatures of germline mtDNA and their value in understanding mitochondrial pathogenicity remain unknown. Here, we carried out a systematic analysis of mutation patterns in germline mtDNA based on 97,566 mtDNA variants from 45,494 full-length sequences and revealed a highly non-stochastic and replication-coupled mutation signature characterized by nucleotide-specific mutation pressure (G > T>A > C) and position-specific selection pressure, suggesting the existence of an intensive mutation-selection interplay in germline mtDNA. We provide evidence that this mutation-selection interplay has strongly shaped the mtDNA sequence during evolution, which not only manifests as an oriented alteration of amino acid compositions of mitochondrial encoded proteins, but also explains the long-lasting mystery of CpG depletion in mitochondrial genome. Finally, we demonstrated that these insights may be integrated to better understand the pathogenicity of disease-implicated mitochondrial variants.