Real-time assessment of mitochondrial DNA heteroplasmy dynamics at the single-cell level.

Mitochondrial DNA (mtDNA) is present in multiple copies within cells and is required for mitochondrial ATP generation. Even within individual cells, mtDNA copies can differ in their sequence, a state known as heteroplasmy. The principles underlying dynamic changes in the degree of heteroplasmy remain incompletely understood, due to the inability to monitor this phenomenon in real time. Here, we employ mtDNA-based fluorescent markers, microfluidics, and automated cell tracking, to follow mtDNA variants in live heteroplasmic yeast populations at the single-cell level. This approach, in combination with direct mtDNA tracking and data-driven mathematical modeling reveals asymmetric partitioning of mtDNA copies during cell division, as well as limited mitochondrial fusion and fission frequencies, as critical driving forces for mtDNA variant segregation. Given that our approach also facilitates assessment of segregation between intact and mutant mtDNA, we anticipate that it will be instrumental in elucidating the mechanisms underlying the purifying selection of mtDNA.

Origins of tissue and cell-type specificity in mitochondrial DNA (mtDNA) disease.

Mutations of mitochondrial (mt)DNA are a major cause of morbidity and mortality in humans, accounting for approximately two thirds of diagnosed mitochondrial disease. However, despite significant advances in technology since the discovery of the first disease-causing mtDNA mutations in 1988, the comprehensive diagnosis and treatment of mtDNA disease remains challenging. This is partly due to the highly variable clinical presentation linked to tissue-specific vulnerability that determines which organs are affected. Organ involvement can vary between different mtDNA mutations, and also between patients carrying the same disease-causing variant. The clinical features frequently overlap with other non-mitochondrial diseases, both rare and common, adding to the diagnostic challenge. Building on previous findings, recent technological advances have cast further light on the mechanisms which underpin the organ vulnerability in mtDNA diseases, but our understanding is far from complete. In this review we explore the origins, current knowledge, and future directions of research in this area.

Heteroplasmy Is Rare in Plant Mitochondria Compared with Plastids despite Similar Mutation Rates.

Plant cells harbor two membrane-bound organelles containing their own genetic material-plastids and mitochondria. Although the two organelles coexist and coevolve within the same plant cells, they differ in genome copy number, intracellular organization, and mode of segregation. How these attributes affect the time to fixation or, conversely, loss of neutral alleles is currently unresolved. Here, we show that mitochondria and plastids share the same mutation rate, yet plastid alleles remain in a heteroplasmic state significantly longer compared with mitochondrial alleles. By analyzing genetic variants across populations of the marine flowering plant Zostera marina and simulating organelle allele dynamics, we examine the determinants of allele segregation and allele fixation. Our results suggest that the bottlenecks on the cell population, e.g. during branching or seeding, and stratification of the meristematic tissue are important determinants of mitochondrial allele dynamics. Furthermore, we suggest that the prolonged plastid allele dynamics are due to a yet unknown active plastid partition mechanism. The dissimilarity between plastid and mitochondrial novel allele fixation at different levels of organization may manifest in differences in adaptation processes. Our study uncovers fundamental principles of organelle population genetics that are essential for further investigations of long-term evolution and molecular dating of divergence events.

Single-mitochondrion sequencing uncovers distinct mutational patterns and heteroplasmy landscape in mouse astrocytes and neurons.

BACKGROUND: Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting "normal" mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes, we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes. In this study, we amplified mt-genomes from 1645 single mitochondria isolated from mouse single astrocytes and neurons to (1) determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, (2) assess differences in mtDNA SNVs between neurons and astrocytes, and (3) study co-segregation of variants in the mouse mtDNA. RESULTS: (1) The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. (2) Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G > A and 9419:C > T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. (3) Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461. CONCLUSIONS: This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.

High heteroplasmy is associated with low mitochondrial copy number and selection against non-synonymous mutations in the snail Cepaea nemoralis.

Molluscan mitochondrial genomes are unusual because they show wide variation in size, radical genome rearrangements and frequently show high variation (> 10%) within species. As progress in understanding this variation has been limited, we used whole genome sequencing of a six-generation matriline of the terrestrial snail Cepaea nemoralis, as well as whole genome sequences from wild-collected C. nemoralis, the sister species C. hortensis, and multiple other snail species to explore the origins of mitochondrial DNA (mtDNA) variation. The main finding is that a high rate of SNP heteroplasmy in somatic tissue was negatively correlated with mtDNA copy number in both Cepaea species. In individuals with under ten mtDNA copies per nuclear genome, more than 10% of all positions were heteroplasmic, with evidence for transmission of this heteroplasmy through the germline. Further analyses showed evidence for purifying selection acting on non-synonymous mutations, even at low frequency of the rare allele, especially in cytochrome oxidase subunit 1 and cytochrome b. The mtDNA of some individuals of Cepaea nemoralis contained a length heteroplasmy, including up to 12 direct repeat copies of tRNA-Val, with 24 copies in another snail, Candidula rugosiuscula, and repeats of tRNA-Thr in C. hortensis. These repeats likely arise due to error prone replication but are not correlated with mitochondrial copy number in C. nemoralis. Overall, the findings provide key insights into mechanisms of replication, mutation and evolution in molluscan mtDNA, and so will inform wider studies on the biology and evolution of mtDNA across animal phyla.

Mitochondrial point heteroplasmy: insights from deep-sequencing of human replicate samples.

BACKGROUND: Human mitochondrial heteroplasmy is an extensively investigated phenomenon in the context of medical diagnostics, forensic identification and molecular evolution. However, technical limitations of high-throughput sequencing hinder reliable determination of point heteroplasmies (PHPs) with minor allele frequencies (MAFs) within the noise threshold. RESULTS: To investigate the PHP landscape at an MAF threshold down to 0.1%, we sequenced whole mitochondrial genomes at approximately 7.700x coverage, in multiple technical and biological replicates of longitudinal blood and buccal swab samples from 11 human donors (159 libraries in total). The results obtained by two independent sequencing platforms and bioinformatics pipelines indicate distinctive PHP patterns below and above the 1% MAF cut-off. We found a high inter-individual prevalence of low-level PHPs (MAF < 1%) at polymorphic positions of the mitochondrial DNA control region (CR), their tissue preference, and a tissue-specific minor allele linkage. We also established the position-dependent potential of minor allele expansion in PHPs, and short-term PHP instability in a mitotically active tissue. We demonstrate that the increase in sensitivity of PHP detection to minor allele frequencies below 1% within a robust experimental and analytical pipeline, provides new information with potential applicative value. CONCLUSIONS: Our findings reliably show different mutational loads between tissues at sub-1% allele frequencies, which may serve as an informative medical biomarker of time-dependent, tissue-specific mutational burden, or help discriminate forensically relevant tissues in a single person, close maternal relatives or unrelated individuals of similar phylogenetic background.

T cell activation contributes to purifying selection against the MELAS-associated m.3243A>G pathogenic variant in blood.

T cells have been shown to maintain a lower percentage (heteroplasmy) of the pathogenic m.3243A>G variant (MT-TL1, associated with maternally inherited diabetes and deafness [MIDD] and mitochondrial encephalomyopathy with lactic acidosis and stroke-like episodes [MELAS]). The mechanism(s) underlying this purifying selection, however, remain unknown. Here we report that purified patient memory CD4+ T cells have lower bulk m.3243A>G heteroplasmy compared to naïve CD4+ T cells. In vitro activation of naïve CD4+ m.3243A>G patient T cells results in lower bulk m.3243A>G heteroplasmy after proliferation. Finally, m.3243A>G patient T cell receptor repertoire sequencing reveals relative oligoclonality compared to controls. These data support a role for T cell activation in peripheral, purifying selection against high m.3243A>G heteroplasmy T cells at the level of the cell, in a likely cell-autonomous fashion.

Numerous insertions of mitochondrial DNA in the genome of the northern mole vole, Ellobius talpinus.

BACKGROUND: Ellobius talpinus is a subterranean rodent representing an attractive model in population ecology studies due to its highly special lifestyle and sociality. In such studies, mitochondrial DNA (mtDNA) is widely used. However, if nuclear copies of mtDNA, aka NUMTs, are present, they may co-amplify with the target mtDNA fragment, generating misleading results. The aim of this study was to determine whether NUMTs are present in E. talpinus. METHODS AND RESULTS: PCR amplification of the putative mtDNA CytB-D-loop fragment using 'universal' primers from 56 E. talpinus samples produced multiple double peaks in 90% of the sequencing chromatograms. To reveal NUMTs, molecular cloning and sequencing of PCR products of three specimens was conducted, followed by phylogenetic analysis. The pseudogene nature of three out of the seven detected haplotypes was confirmed by their basal positions in relation to other Ellobius haplotypes in the phylogenetic tree. Additionally, 'haplotype B' was basal in relation to other E. talpinus haplotypes and found present in very distant sampling sites. BLASTN search revealed 195 NUMTs in the E. talpinus nuclear genome, including fragments of all four PCR amplified pseudogenes. Although the majority of the NUMTs studied were short, the entire mtDNA had copies in the nuclear genome. The most numerous NUMTs were found for rrnL, COXI, and D-loop. CONCLUSIONS: Numerous NUMTs are present in E. talpinus and can be difficult to discriminate against mtDNA sequences. Thus, in future population or phylogenetic studies in E. talpinus, the possibility of cryptic NUMTs amplification should always be taken into account.

Mitochondrial DNA Variants at Low-Level Heteroplasmy and Decreased Copy Numbers in Chronic Kidney Disease (CKD) Tissues with Kidney Cancer.

The human mitochondrial genome (mtDNA) is a circular DNA molecule with a length of 16.6 kb, which contains a total of 37 genes. Somatic mtDNA mutations accumulate with age and environmental exposure, and some types of mtDNA variants may play a role in carcinogenesis. Recent studies observed mtDNA variants not only in kidney tumors but also in adjacent kidney tissues, and mtDNA dysfunction results in kidney injury, including chronic kidney disease (CKD). To investigate whether a relationship exists between heteroplasmic mtDNA variants and kidney function, we performed ultra-deep sequencing (30,000×) based on long-range PCR of DNA from 77 non-tumor kidney tissues of kidney cancer patients with CKD (stages G1 to G5). In total, this analysis detected 697 single-nucleotide variants (SNVs) and 504 indels as heteroplasmic (0.5% ≤ variant allele frequency (VAF) < 95%), and the total number of detected SNVs/indels did not differ between CKD stages. However, the number of deleterious low-level heteroplasmic variants (pathogenic missense, nonsense, frameshift and tRNA) significantly increased with CKD progression (p < 0.01). In addition, mtDNA copy numbers (mtDNA-CNs) decreased with CKD progression (p < 0.001). This study demonstrates that mtDNA damage, which affects mitochondrial genes, may be involved in reductions in mitochondrial mass and associated with CKD progression and kidney dysfunction.

Mitochondrial variants of complex I genes associated with leprosy clinical subtypes.

Leprosy is a chronic bacterial infection mainly caused by Mycobacterium leprae that primarily affects skin and peripheral nerves. Due to its ability to absorb carbon from the host cell, the bacillus became dependent on energy production, mainly through oxidative phosphorylation. In fact, variations in genes of Complex I of oxidative phosphorylation encoded by mtDNA have been associated with several diseases in humans, including bacterial infections, which are possible influencers in the host response to leprosy. Here, we investigated the presence of variants in the mtDNA genes encoding Complex I regarding leprosy, as well as the analysis of their pathogenicity in the studied cohort. We found an association of 74 mitochondrial variants with either of the polar forms, Pole T (Borderline Tuberculoid) or Pole L (Borderline Lepromatous and Lepromatous) of leprosy. Notably, six variants were exclusively found in both clinical poles of leprosy, including m.4158A>G and m.4248T>C in MT-ND1, m.13650C>A, m.13674T>C, m.12705C>T and m.13263A>G in MT-ND5, of which there are no previous reports in the global literature. Our observations reveal a substantial number of mutations among different groups of leprosy, highlighting a diverse range of consequences associated with mutations in genes across these groups. Furthermore, we suggest that the six specific variants exclusively identified in the case group could potentially play a crucial role in leprosy susceptibility and its clinical differentiation. These variants are believed to contribute to the instability and dysregulation of oxidative phosphorylation during the infection, further emphasizing their significance.

Association analysis of mitochondrial DNA heteroplasmic variants: methods and application.

We rigorously assessed a comprehensive association testing framework for heteroplasmy, employing both simulated and real-world data. This framework employed a variant allele fraction (VAF) threshold and harnessed multiple gene-based tests for robust identification and association testing of heteroplasmy. Our simulation studies demonstrated that gene-based tests maintained an appropriate type I error rate at α=0.001. Notably, when 5% or more heteroplasmic variants within a target region were linked to an outcome, burden-extension tests (including the adaptive burden test, variable threshold burden test, and z-score weighting burden test) outperformed the sequence kernel association test (SKAT) and the original burden test. Applying this framework, we conducted association analyses on whole-blood derived heteroplasmy in 17,507 individuals of African and European ancestries (31% of African Ancestry, mean age of 62, with 58% women) with whole genome sequencing data. We performed both cohort- and ancestry-specific association analyses, followed by meta-analysis on both pooled samples and within each ancestry group. Our results suggest that mtDNA-encoded genes/regions are likely to exhibit varying rates in somatic aging, with the notably strong associations observed between heteroplasmy in the RNR1 and RNR2 genes (p<0.001) and advance aging by the Original Burden test. In contrast, SKAT identified significant associations (p<0.001) between diabetes and the aggregated effects of heteroplasmy in several protein-coding genes. Further research is warranted to validate these findings. In summary, our proposed statistical framework represents a valuable tool for facilitating association testing of heteroplasmy with disease traits in large human populations.

Chloroplast transformation in new cultivars of tomato through particle bombardment.

A protocol has been established for genetic transformation of the chloroplasts in two new cultivars of tomato (Solanum lycopersicum L.) grown in India and Australia: Pusa Ruby and Yellow Currant. Tomato cv. Green Pineapple was also used as a control that has previously been used for establishing chloroplast transformation by other researchers. Selected tomato cultivars were finalized from ten other tested cultivars (Green Pineapple excluded) due to their high regeneration potential and better response to chloroplast transformation. This protocol was set up using a chloroplast transformation vector (pRB94) for tomatoes that is made up of a synthetic gene operon. The vector has a chimeric aadA selectable marker gene that is controlled by the rRNA operon promoter (Prrn). This makes the plant or chloroplasts resistant to spectinomycin and streptomycin. After plasmid-coated particle bombardment, leaf explants were cultured in 50 mg/L selection media. Positive explant selection from among all the dead-appearing (yellow to brown) explants was found to be the major hurdle in the study. Even though this study was able to find plastid transformants in heteroplasmic conditions, it also found important parameters and changes that could speed up the process of chloroplast transformation in tomatoes, resulting in homoplasmic plastid-transformed plants. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-024-03954-3.

Approaches and challenges in identifying, quantifying, and manipulating dynamic mitochondrial genome variations.

Mitochondria, the cellular powerhouses, possess their own unique genetic system, including replication, transcription, and translation. Studying these processes is crucial for comprehending mitochondrial disorders, energy production, and their related diseases. Over the past decades, various approaches have been applied in detecting and quantifying mitochondrial genome variations with also the purpose of manipulation of mitochondria or mitochondrial genome for therapeutics. Understanding the scope and limitations of above strategies is not only fundamental to the understanding of basic biology but also critical for exploring disease-related novel target(s), as well to develop innovative therapies. Here, this review provides an overview of different tools and techniques for accurate mitochondrial genome variations identification, quantification, and discuss novel strategies for the manipulation of mitochondria to develop innovative therapeutic interventions, through combining the insights gained from the study of mitochondrial genetics with ongoing single cell omics combined with advanced single molecular tools.

Vitamin B12 status and folic acid supplementation influence mitochondrial heteroplasmy levels in mice.

One-carbon metabolism is a complex network of metabolic reactions that are essential for cellular function including DNA synthesis. Vitamin B12 and folate are micronutrients that are utilized in this pathway and their deficiency can result in the perturbation of one-carbon metabolism and subsequent perturbations in DNA replication and repair. This effect has been well characterized in nuclear DNA but to date, mitochondrial DNA (mtDNA) has not been investigated extensively. Mitochondrial variants have been associated with several inherited and age-related disease states; therefore, the study of factors that impact heteroplasmy are important for advancing our understanding of the mitochondrial genome's impact on human health. Heteroplasmy studies require robust and efficient mitochondrial DNA enrichment to carry out in-depth mtDNA sequencing. Many of the current methods for mtDNA enrichment can introduce biases and false-positive results. Here, we use a method that overcomes these limitations and have applied it to assess mitochondrial heteroplasmy in mouse models of altered one-carbon metabolism. Vitamin B12 deficiency was found to cause increased levels of mitochondrial DNA heteroplasmy across all tissues that were investigated. Folic acid supplementation also contributed to elevated mitochondrial DNA heteroplasmy across all mouse tissues investigated. Heteroplasmy analysis of human data from the Framingham Heart Study suggested a potential sex-specific effect of folate and vitamin B12 status on mitochondrial heteroplasmy. This is a novel relationship that may have broader consequences for our understanding of one-carbon metabolism, mitochondrial-related disease and the influence of nutrients on DNA mutation rates.

AMPK Suppression Due to Obesity Drives Oocyte mtDNA Heteroplasmy via ATF5-POLG Axis.

Due to the exclusive maternal transmission, oocyte mitochondrial dysfunction reduces fertility rates, affects embryonic development, and programs offspring to metabolic diseases. However, mitochondrial DNA (mtDNA) are vulnerable to mutations during oocyte maturation, leading to mitochondrial nucleotide variations (mtSNVs) within a single oocyte, referring to mtDNA heteroplasmy. Obesity (OB) accounts for more than 40% of women at the reproductive age in the USA, but little is known about impacts of OB on mtSNVs in mature oocytes. It is found that OB reduces mtDNA content and increases mtSNVs in mature oocytes, which impairs mitochondrial energetic functions and oocyte quality. In mature oocytes, OB suppresses AMPK activity, aligned with an increased binding affinity of the ATF5-POLG protein complex to mutated mtDNA D-loop and protein-coding regions. Similarly, AMPK knockout increases the binding affinity of ATF5-POLG proteins to mutated mtDNA, leading to the replication of heteroplasmic mtDNA and impairing oocyte quality. Consistently, AMPK activation blocks the detrimental impacts of OB by preventing ATF5-POLG protein recruitment, improving oocyte maturation and mitochondrial energetics. Overall, the data uncover key features of AMPK activation in suppressing mtSNVs, and improving mitochondrial biogenesis and oocyte maturation in obese females.

Quantifying mitochondrial heteroplasmy diversity: A computational approach.

Biodiversity plays a pivotal role in sustaining ecosystem processes, encompassing diverse biological species, genetic types and the intricacies of ecosystem composition. However, the precise definition of biodiversity at the individual level remains a challenging endeavour. Hill numbers, derived from Rényi's entropy, have emerged as a popular measure of diversity, with a recent unified framework extending their application across various levels, from genetics to ecosystems. In this study, we employ a computational approach to exploring the diversity of mitochondrial heteroplasmy using real-world data. By adopting Hill numbers with q = 2, we demonstrate the feasibility of quantifying mitochondrial heteroplasmy diversity within and between individuals and populations. Furthermore, we investigate the alpha diversity of mitochondrial heteroplasmy among different species, revealing heterogeneity at multiple levels, including mitogenome components and protein-coding genes (PCGs). Our analysis explores large-scale mitochondrial heteroplasmy data in humans, examining the relationship between alpha diversity at the mitogenome components and PCGs level. Notably, we do not find a significant correlation between these two levels. Additionally, we observe significant correlations in alpha diversity between mothers and children in blood samples, exceeding the reported R2 value for allele frequency correlations. Moreover, our investigation of beta diversity and local overlay similarity demonstrates that heteroplasmy variant distributions in different tissues of children more closely resemble those of their mothers. Through systematic quantification and analysis of mitochondrial heteroplasmy diversity, this study enhances our understanding of heterogeneity at multiple levels, from individuals to populations, providing new insights into this fundamental dimension of biodiversity.

Single-Mitochondrion Sequencing Uncovers Distinct Mutational Patterns and Heteroplasmy Landscape in Mouse Astrocytes and Neurons.

Background: Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting 'normal' mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes In this study we amplified mt-genomes from 1,645 single mitochondria (mts) isolated from mouse single astrocytes and neurons to 1. determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, 2. assess differences in mtDNA SNVs between neurons and astrocytes, and 3. Study cosegregation of variants in the mouse mtDNA. Results: 1. The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. 2. Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G>A and 9419:C>T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. 3. Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461. Conclusion: This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.

Friend turned foe: selfish behavior of a spontaneously arising mitochondrial deletion in an experimentally evolved Caenorhabditis elegans population.

Selfish mitochondrial DNA (mtDNA) mutations are variants that can proliferate within cells and enjoy a replication or transmission bias without fitness benefits for the host. mtDNA deletions in Caenorhabditis elegans can reach high heteroplasmic frequencies despite significantly reducing fitness, illustrating how new mtDNA variants can give rise to genetic conflict between different levels of selection and between the nuclear and mitochondrial genomes. During a mutation accumulation experiment in C. elegans, a 1,034-bp deletion originated spontaneously and reached an 81.7% frequency within an experimental evolution line. This heteroplasmic mtDNA deletion, designated as meuDf1, eliminated portions of 2 protein-coding genes (coxIII and nd4) and tRNA-thr in entirety. mtDNA copy number in meuDf1 heteroplasmic individuals was 35% higher than in individuals with wild-type mitochondria. After backcrossing into a common genetic background, the meuDf1 mitotype was associated with reduction in several fitness traits and independent competition experiments found a 40% reduction in composite fitness. Experiments that relaxed individual selection by single individual bottlenecks demonstrated that the deletion-bearing mtDNA possessed a strong transmission bias, thereby qualifying it as a novel selfish mitotype.

Characterization and Trans-generation Dynamics of Mitogene Pool in the Silver Carp (Hypophthalmichthys molitrix).

Multi-copied mitogenome are prone to mutation during replication often resulting in heteroplasmy. The derived variants in a cell, organ or an individual animal constitute a mitogene pool. The individual mitogene pool is initiated by a small fraction of the egg mitogene pool. However, the characteristics and relationship between them has not yet been investigated. This study quantitatively analyzed the heteroplasmy landscape, genetic loads, and selection strength of the mitogene pool of egg and hatchling in the silver carp (Hypophthalmichthys molitrix) using high-throughput resequencing. The results showed heteroplasmic sites distribute across the whole mitogenome in both eggs and hatchlings. The dominant substitution was Transversion in eggs and Transition in hatching accounting for 95.23% ± 2.07% and 85.38% ± 6.94% of total HP sites, respectively. The total genetic loads were 0.293 ± 0.044 in eggs and 0.228 ± 0.022 in hatchlings (p = 0.048). The dN/dS ratio was 58.03 ± 38.98 for eggs and 9.44 ± 3.93 for hatchlings (p = 0.037). These results suggest that the mitogenomes were under strong positive selection in eggs with tolerance to variants with deleterious effects, while the selection was positive but much weaker in hatchlings showing marked quality control. Based on these findings, we proposed a trans-generation dynamics model to explain differential development mode of the two mitogene pool between oocyte maturation and ontogenesis of offspring. This study sheds light on significance of mitogene pool for persistence of populations and subsequent integration in ecological studies and conservation practices.