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

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. f.

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.

Increases in computationally predicted deleterious variants of unknown significance and sperm mtDNA copy numbers may be associated with semen quality

BACKGROUND

Mitochondria are essential for sperm motility because they provide the energy required for the movement. Changes in sperm mtDNA, such as point mutations, large-scale deletions, or copy number variations, may interfere with ATP production and reduce sperm motility. However, it is not clear if changes in mtDNA are linked to semen quality.

OBJECTIVES

To explore the association between sperm mitochondrial DNA (mtDNA) changes and semen quality.

MATERIALS AND METHODS

Sixty-five oligo and/or astheno and/or terato patients (O/A/T) patients and 41 controls were recruited from couples undergoing assisted reproduction. Semen and blood samples were collected from the same individual on the day of oocyte retrieval to extract, isolate and purify mtDNA for next-generation sequencing. mtDNA copy numbers were assessed in 64 patient and 39 control sperm DNA samples using quantitative real-time PCR. The 4977 bp deletion was assessed in 20 patient and 20 control sperm DNA samples using polymerase chain reaction.

RESULTS

The mtDNA of patients was more likely to carry pathogenic variants or variants of unknown significance (VUSs) (P = 0.091) with higher heteroplasmy levels (P < 0.05) than that of controls. Interestingly, 33.85% of O/A/T patients (22 out of 65) lacked unique variants in their spermatozoa. but presented an exceptionally high mtDNA copy number (P < 0.0001). Moreover, we observed a decrease in the heteroplasmy level of common mtDNA variants shared by somatic and gamete cells (P < 0.0001) and the emergence of a very large number of de novo mtDNA variants with low-level heteroplasmy in spermatozoa.

DISCUSSION AND CONCLUSION

The increases in the number of computationally predicted deleterious VUS and mtDNA copies in spermatozoa may be associated with semen quality. Exposure to environmental mutation pressure that causes novel mtDNA variants with low-level heteroplasmy may occur during spermatogenesis. Furthermore, when a certain harmful threshold is reached, male germ cells may degrade mtDNA with mutations and replicate the correct mtDNA sequence to maintain the mitochondrial function in spermatozoa.

Heteroplasmy and Individual Mitogene Pools: Characteristics and Potential Roles in Ecological Studies

The mitochondrial genome (mitogenome or mtDNA), the extrachromosomal genome, is a multicopy circular DNA with high mutation rates due to replication and repair errors. A mitochondrion, cell, tissue, organ, or an individual body may hold multiple variants, both inherited and developed over a lifetime, which make up individual mitogene pools. This phenomenon is also called mtDNA heteroplasmy. MtDNA variants influence cellular and tissular functions and are consequently subjected to selection. Although it has long been recognized that only inheritable germline heteroplasmies have evolutionary significance, non-inheritable somatic heteroplasmies have been overlooked since they directly affect individual fitness and thus indirectly affect the fate of heritable germline variants. This review focuses on the characteristics, dynamics, and functions of mtDNA heteroplasmy and proposes the concept of individual mitogene pools to discuss individual genetic diversity from multiple angles. We provide a unique perspective on the relationship between individual genetic diversity and heritable genetic diversity and guide how the individual mitogene pool with novel genetic markers can be applied to ecological research.

Heteroplasmic mitochondrial DNA variants in cardiovascular diseases

Mitochondria are implicated in the pathogenesis of cardiovascular diseases (CVDs) but the reasons for this are not well understood. Maternally-inherited population variants of mitochondrial DNA (mtDNA) which affect all mtDNA molecules (homoplasmic) are associated with cardiometabolic traits and the risk of developing cardiovascular disease. However, it is not known whether mtDNA mutations only affecting a proportion of mtDNA molecules (heteroplasmic) also play a role. To address this question, we performed a high-depth (~1000-fold) mtDNA sequencing of blood DNA in 1,399 individuals with hypertension (HTN), 1,946 with ischemic heart disease (IHD), 2,146 with ischemic stroke (IS), and 723 healthy controls. We show that the per individual burden of heteroplasmic single nucleotide variants (mtSNVs) increases with age. The age-effect was stronger for low-level heteroplasmies (heteroplasmic fraction, HF, 5-10%), likely reflecting acquired somatic events based on trinucleotide mutational signatures. After correcting for age and other confounders, intermediate heteroplasmies (HF 10-95%) were more common in hypertension, particularly involving non-synonymous variants altering the amino acid sequence of essential respiratory chain proteins. These findings raise the possibility that heteroplasmic mtSNVs play a role in the pathophysiology of hypertension.

A replication-linked mutational gradient drives somatic mutation accumulation and influences germline polymorphisms and genome composition in mitochondrial DNA

Mutations in mitochondrial DNA (mtDNA) cause maternally inherited diseases, while somatic mutations are linked to common diseases of aging. Although mtDNA mutations impact health, the processes that give rise to them are under considerable debate. To investigate the mechanism by which de novo mutations arise, we analyzed the distribution of naturally occurring somatic mutations across the mouse and human mtDNA obtained by Duplex Sequencing. We observe distinct mutational gradients in G→A and T→C transitions delimited by the light-strand origin and the mitochondrial Control Region (mCR). The gradient increases unequally across the mtDNA with age and is lost in the absence of DNA polymerase γ proofreading activity. In addition, high-resolution analysis of the mCR shows that important regulatory elements exhibit considerable variability in mutation frequency, consistent with them being mutational ‘hot-spots’ or ‘cold-spots’. Collectively, these patterns support genome replication via a deamination prone asymmetric strand-displacement mechanism as the fundamental driver of mutagenesis in mammalian DNA. Moreover, the distribution of mtDNA single nucleotide polymorphisms in humans and the distribution of bases in the mtDNA across vertebrate species mirror this gradient, indicating that replication-linked mutations are likely the primary source of inherited polymorphisms that, over evolutionary timescales, influences genome composition during speciation.

Cell reprogramming shapes the mitochondrial DNA landscape

Individual induced pluripotent stem cells (iPSCs) show considerable phenotypic heterogeneity, but the reasons for this are not fully understood. Comprehensively analysing the mitochondrial genome (mtDNA) in 146 iPSC and fibroblast lines from 151 donors, we show that most age-related fibroblast mtDNA mutations are lost during reprogramming. However, iPSC-specific mutations are seen in 76.6% (108/141) of iPSC lines at a mutation rate of 8.62 × 10⁻⁵/base pair. The mutations observed in iPSC lines affect a higher proportion of mtDNA molecules, favouring non-synonymous protein-coding and tRNA variants, including known disease-causing mutations. Analysing 11,538 single cells shows stable heteroplasmy in sub-clones derived from the original donor during differentiation, with mtDNA variants influencing the expression of key genes involved in mitochondrial metabolism and epidermal cell differentiation. Thus, the dynamic mtDNA landscape contributes to the heterogeneity of human iPSCs and should be considered when using reprogrammed cells experimentally or as a therapy.

Here the authors describe high depth mitochondrial DNA (mtDNA) sequence analysis of 146 human induced pluripotent stem cell (hiPSC) lines as well as single cell RNA-seq (scRNAseq) of hiPSCs undergoing differentiation from 125 donors; reporting mtDNA diversity and some variants favoured after reprogramming.

Evolutionary Insights Into Two Widespread Ectomycorrhizal Fungi (Pisolithus) From Comparative Analysis of Mitochondrial Genomes

The genus Pisolithus is a group of global ectomycorrhizal fungi. The characterizations of Pisolithus mitochondrial genomes have still been unknown. In the present study, the complete mitogenomes of two Pisolithus species, Pisolithus microcarpus, and Pisolithus tinctorius, were assembled and compared with other Boletales mitogenomes. Both Pisolithus mitogenomes comprised circular DNA molecules with sizes of 43,990 and 44,054 bp, respectively. Comparative mitogenomic analysis showed that the rps3 gene differentiated greatly between Boletales species, and this gene may be subjected to strong pressure of positive selection between some Boletales species. Several plasmid-derived genes and genes with unknown functions were detected in the two Pisolithus mitogenomes, which needs further analysis. The two Pisolithus species show a high degree of collinearity, which may represent the gene arrangement of the ancestors of ectomycorrhizal Boletales species. Frequent intron loss/gain events were detected in Boletales and basidiomycetes, and intron P717 was only detected in P. tinctorius out of the eight Boletales mitogenomes tested. We reconstructed phylogeny of 79 basidiomycetes based on combined mitochondrial gene dataset, and obtained well-supported phylogenetic topologies. This study served as the first report on the mitogenomes of the family Pisolithaceae, which will promote the understanding of the evolution of Pisolithus species.

Panorama of intron dynamics and gene rearrangements in the phylum Basidiomycota as revealed by the complete mitochondrial genome of Turbinellus floccosus

In the present study, the complete mitogenome of Turbinellus floccosus was sequenced, assembled, and compared with other basidiomycete mitogenomes. The mitogenome of T. floccosus consists of a circular DNA molecule, with a size of 62,846 bp. Gene arrangement analysis indicated that large-scale gene rearrangements occurred in the levels of family and genus of basidiomycete species, and the mitogenome of T. floccosus contained a unique gene order. A significant correlation between the number of introns and the mitochondrial genome size of Basidiomycota were detected (P < 0.01). A total of 896 introns were detected in the core protein-coding genes (PCGs) of 74 basidiomycete species, and the cox1 gene was the largest host gene of basidiomycete introns. Intron position class (Pcls) P383 in the cox1 gene was the most common intron in Basidiomycota, which distributed in 40 of 74 basidiomycete species. In addition, frequent intron loss/gain events were detected in basidiomycete species. More than 50% of bases around insertion sites (- 15 bp to 15 bp) of Pcls from different species were conservative, indicating site preferences of intron insertions in Basidiomycota. Further analysis showed that 76.09% of introns tended to insert downstream to a T base in Basidiomycota. Phylogenetic analysis for 74 basidiomycetes indicated mitochondrial genes are effective molecular markers for phylogeny of basidiomycetes. The study served as the first report on the mitogenome from the family Gomphaceae, which will help to understand the intron origin and evolution in Basidiomycota. KEY POINTS: • The mitogenome of Turbinellus floccosus had a unique gene arrangement. • Intron loss/gain events were detected in the 74 basidiomycete species. • Introns tend to insert downstream of a T base in basidiomycete mitogenomes.

The mitogenomes of two saprophytic Boletales species (Coniophora) reveals intron dynamics and accumulation of plasmid-derived and non-conserved genes

The order Boletales is a group of fungi with complex life styles, which include saprophytic and ectomycorrhizal mushroom-forming fungi. In the present study, the complete mitogenomes of two saprophytic Boletales species, Coniophora olivacea, and C. puteana, were assembled and compared with mitogenomes of ectomycorrhizal Boletales. Both mitogenomes comprised circular DNA molecules with sizes of 78,350 bp and 79,655 bp, respectively. Comparative mitogenomic analysis indicated that the two saprophytic Boletales species contained more plasmid-derived (7 on average) and unknown functional genes (12 on average) than the four ectomycorrhizal Boletales species previously reported. In addition, the core protein coding genes, nad2 and rps3, were found to be subjected to positive selection pressure between some Boletales species. Frequent intron gain/loss events were detected in Boletales and Basidiomycetes, and several novel intron classes were found in two Coniophora species. A total of 33 introns were detected in C. olivacea, and most were found to have undergone contraction in the C. olivacea mitogenome. Mitochondrial genes of Coniophora species were found to have undergone large-scale gene rearrangements, and the accumulation of intra-genomic repeats in the mitogenome was considered as one of the main contributing factors. Based on combined mitochondrial gene sets, we obtained a well-supported phylogenetic tree for 76 Basidiomycetes, demonstrating the utility of mitochondrial gene analysis for inferring Basidiomycetes phylogeny. The study served as the first report on the mitogenomes of the family Coniophorineae, which will help to understand the origin and evolution patterns of Boletales species with complex lifestyles.

Mitochondrial DNA Haplotypes as Genetic Modifiers of Cancer

Mitochondria play an essential role in cellular metabolism, generation of reactive oxygen species (ROS), and the initiation of apoptosis. These properties enable mitochondria to be crucial integrators in the pathways of tumorigenesis. An open question is to what extent variation in the mitochondrial genome (mtDNA) contributes to the biological heterogeneity observed in human tumors. In this review, we summarize our current understanding of the role of mtDNA genetics in relation to human cancers.

Optimized high-fidelity 3DPCR to assess potential mitochondrial targeting by activation-induced cytidine deaminase

Activation‐induced cytidine deaminase (AID) initiates somatic hypermutation and class switch recombination of immunoglobulin genes in B cells, whereas off‐targeted AID activity contributes to oncogenic mutations and chromosomal translocations associated with B cell malignancies. Paradoxically, only a minority of AID is allowed to access the nuclear genome, but the majority of AID is retained in the cytoplasm. It is unknown whether cytoplasmic AID can access and target the mitochondrial genome [mitochondrial DNA (mtDNA)]. To address this issue, we developed high‐fidelity differential DNA denaturation PCR, which allowed the enrichment of genuine mtDNA mutations and therefore the identification of endogenous mtDNA mutation signatures in vitro. With this approach, we showed that AID targeting to mtDNA is a rare event in AID‐expressing lymphoma lines. Further biochemical and microscopic analysis revealed that a fraction of cytosol AID is associated with the outer membrane of mitochondria but unable to access the mitochondrial matrix. Together, our data suggested that the mitochondrial genome is protected from AID‐mediated mutagenesis by physical segregation of AID from accessing mtDNA within the mitochondrial matrix.