Heteroplasmic Variants Research Articles

Characterization of the Mitochondrial Genetic Landscape in Abdominal Aortic Aneurysm.

Background Abdominal aortic aneurysm (AAA) is a vascular disease with a mortality rate of >80% if ruptured. Mitochondrial dysfunction has been previously implicated in AAA pathogenesis. In this study, we aimed to characterize the mitochondrial genetic landscape in AAA. Methods and Results Whole mitochondrial genome sequencing and bioinformatics analysis were performed in comorbidity matched 48 cases without AAA and 48 cases with AAA, objectively diagnosed, and selected from a cohort of 65-year-old men recruited for a screening program. We identified differential mutational landscapes in men with and without AAA, with errors in mitochondrial DNA replication or repair as potential sources. Heteroplasmic insertions and overall heteroplasmy of structural rearrangements were significantly elevated in AAA cases. Three heteroplasmic variants were associated with risk factors of AAA: leukocyte concentration, plasma glucose, and cholesterol levels, respectively. Interestingly, mutations were more prevalent in regulatory part of the mitochondria, the displacement loop region, in AAA as compared with controls (P value <0.05), especially in the conserved and critical mitochondrial extended termination-associated sequenceregion. Moreover, we report a novel 24 bp mitochondrial DNA duplication present exclusively in cases with AAA (4%) and 75% of the unmatched AAA biopsies. Finally, the haplogroup cluster JTU was overrepresented in AAA and significantly associated with a positive family history of AAA (odds ratio, 2.9 [95% CI, 1.1-8.1]). Conclusions This is the first study investigating the mitochondrial genome in AAA, where important genetic alterations and haplogroups associated with AAA and clinical risk factors were identified. Our findings have the potential to fill in gaps in the missing genetic information on AAA.

Quantitative haplotype-resolved analysis of mitochondrial DNA heteroplasmy in Human single oocytes, blastoids, and pluripotent stem cells.

Maternal mitochondria are the sole source of mtDNA for every cell of the offspring. Heteroplasmic mtDNA mutations inherited from the oocyte are a common cause of metabolic diseases and associated with late-onset diseases. However, the origin and dynamics of mtDNA heteroplasmy remain unclear. We used our individual Mitochondrial Genome sequencing (iMiGseq) technology to study mtDNA heterogeneity, quantitate single nucleotide variants (SNVs) and large structural variants (SVs), track heteroplasmy dynamics, and analyze genetic linkage between variants at the individual mtDNA molecule level in single oocytes and human blastoids. Our study presented the first single-mtDNA analysis of the comprehensive heteroplasmy landscape in single human oocytes. Unappreciated levels of rare heteroplasmic variants well below the detection limit of conventional methods were identified in healthy human oocytes, of which many are reported to be deleterious and associated with mitochondrial disease and cancer. Quantitative genetic linkage analysis revealed dramatic shifts of variant frequency and clonal expansions of large SVs during oogenesis in single-donor oocytes. iMiGseq of a single human blastoid suggested stable heteroplasmy levels during early lineage differentiation of naïve pluripotent stem cells. Therefore, our data provided new insights of mtDNA genetics and laid a foundation for understanding mtDNA heteroplasmy at early stages of life.

Single-cell individual full-length mtDNA sequencing by iMiGseq uncovers unexpected heteroplasmy shifts in mtDNA editing.

The ontogeny and dynamics of mtDNA heteroplasmy remain unclear due to limitations of current mtDNA sequencing methods. We developed individual Mitochondrial Genome sequencing (iMiGseq) of full-length mtDNA for ultra-sensitive variant detection, complete haplotyping, and unbiased evaluation of heteroplasmy levels, all at the individual mtDNA molecule level. iMiGseq uncovered unappreciated levels of heteroplasmic variants in single cells well below the conventional NGS detection limit and provided accurate quantitation of heteroplasmy level. iMiGseq resolved the complete haplotype of individual mtDNA in single oocytes and revealed genetic linkage of de novo mutations. iMiGseq detected sequential acquisition of detrimental mutations, including large deletions, in defective mtDNA in NARP/Leigh syndrome patient-derived induced pluripotent stem cells. iMiGseq identified unintended heteroplasmy shifts in mitoTALEN editing, while showing no appreciable level of unintended mutations in DdCBE-mediated mtDNA base editing. Therefore, iMiGseq could not only help elucidate the mitochondrial etiology of diseases, but also evaluate the safety of various mtDNA editing strategies.

Molecular characterization of variants in mitochondrial DNA encoded genes using next generation sequencing analysis and mitochondrial dysfunction in women with PCOS

Emerging studies indicates mitochondrial dysfunction and involvement of mitochondrial DNA (mtDNA) variants in the pathogenesis of polycystic ovary syndrome (PCOS). Cumulative effect of mtDNA rare variants are now gaining considerable interest apart from common variants in the pathogenesis of complex diseases. Rare variants may modify the effect of polymorphism or in combination with the common variants may affect the risk of disease. With the evolution of high throughput sequencing techniques, which can be utilized to identify common as well as rare variants along with heteroplasmy levels, comprehensive characterization of the mtDNA variants is possible. Till date, few studies reported common mtDNA variants using traditional sequencing techniques but rare variants in mtDNA encoding genes remain unexplored in women with PCOS. These mtDNA variants may be responsible for mitochondrial dysfunction and may contribute in PCOS pathogenesis. In this study we determined mtDNA copy number, a biomarker of mitochondrial dysfunction and first time analysed variants in mtDNA encoded genes in women with PCOS using mitochondrial Next Generation sequencing (NGS) approach and compared allele frequency from mitochondrial 1000 genome dataset. Variant annotation and prioritization was done using highly automated pipeline, MToolBox that excludes reads mapped from nuclear mitochondrial DNA sequences (NumtS) to identify unique mtDNA reads. The present study identified significant reduction in mtDNA copy number in women with PCOS compared to non-PCOS women. A total of unique 214 prioritized common to rare variants were identified in mtDNA encoded genes, 183 variants in OXPHOS complexes, 14 variants in MT-tRNA and 17 variants in MT-rRNA genes that may be involved in mitochondrial dysfunction in PCOS. Numerous variants were heteroplasmic, pathogenic in nature and occurred in evolutionary conserved region. Heteroplasmic variants were more frequently occurred in MT-CO3 gene. Non-synonymous variants were more than synonymous variants and mainly occurred in OXPHOS complex I and IV. Few variants were found to be associated with diseases in MITOMAP database. The study provides a better understanding towards pathogenesis of PCOS from novel aspects focusing on mitochondrial genetic defects as underlying cause for contributing mitochondrial dysfunction in women with PCOS.

Mitochondrial DNA heteroplasmy distinguishes disease manifestation in PINK1/PRKN-linked Parkinson's disease.

Biallelic mutations in PINK1/PRKN cause recessive Parkinson's disease. Given the established role of PINK1/Parkin in regulating mitochondrial dynamics, we explored mitochondrial DNA integrity and inflammation as disease modifiers in carriers of mutations in these genes. Mitochondrial DNA integrity was investigated in a large collection of biallelic (n = 84) and monoallelic (n = 170) carriers of PINK1/PRKN mutations, idiopathic Parkinson's disease patients (n = 67) and controls (n = 90). In addition, we studied global gene expression and serum cytokine levels in a subset. Affected and unaffected PINK1/PRKN monoallelic mutation carriers can be distinguished by heteroplasmic mitochondrial DNA variant load (area under the curve = 0.83, CI 0.74-0.93). Biallelic PINK1/PRKN mutation carriers harbour more heteroplasmic mitochondrial DNA variants in blood (P = 0.0006, Z = 3.63) compared to monoallelic mutation carriers. This enrichment was confirmed in induced pluripotent stem cell-derived (controls, n = 3; biallelic PRKN mutation carriers, n = 4) and post-mortem (control, n = 1; biallelic PRKN mutation carrier, n = 1) midbrain neurons. Last, the heteroplasmic mitochondrial DNA variant load correlated with IL6 levels in PINK1/PRKN mutation carriers (r = 0.57, P = 0.0074). PINK1/PRKN mutations predispose individuals to mitochondrial DNA variant accumulation in a dose- and disease-dependent manner.

Phylogenomic analyses based on the plastid genome and concatenated nrDNA sequence data reveal cytonuclear discordance in genus Atractylodes (Asteraceae: Carduoideae).

Atractylodes species are widely distributed across East Asia and are cultivated as medicinal herbs in China, Japan, and Korea. Their unclear morphological characteristics and low levels of genetic divergence obscure the taxonomic relationships among these species. In this study, 24 plant samples were collected representing five species of Atractylodes located in China; of these, 23 belonged to members of the A. lancea complex. High-throughput sequencing was used to obtain the concatenated nrDNA sequences (18S-ITS1-5.8S-ITS2-28S) and plastid genomes. The concatenated nrDNA sequence lengths for all the Atractylodes species were 5,849 bp, and the GC content was 55%. The lengths of the whole plastid genome sequences ranged from 152,138 bp (A. chinensis) to 153,268 bp (A. lancea), while their insertion/deletion sites were mainly distributed in the intergenic regions. Furthermore, 33, 34, 36, 31, and 32 tandem repeat sequences, as well as 30, 30, 29, 30, and 30 SSR loci, were detected in A. chinensis, A. koreana, A. lancea, A. japonica, and A. macrocephala, respectively. In addition to these findings, a considerable number of heteroplasmic variations were detected in the plastid genomes, implying a complicated phylogenetic history for Atractylodes. The results of the phylogenetic analysis involving concatenated nrDNA sequences showed that A. lancea and A. japonica formed two separate clades, with A. chinensis and A. koreana constituting their sister clade, while A. lancea, A. koreana, A. chinensis, and A. japonica were found based on plastid datasets to represent a mixed clade on the phylogenetic tree. Phylogenetic network analysis suggested that A. lancea may have hybridized with the common ancestor of A. chinensis and A. japonica, while ABBA-BABA tests of SNPs in the plastid genomes showed that A. chinensis was more closely related to A. japonica than to A. lancea. This study reveals the extensive discordance and complexity of the relationships across the members of the A. lancea complex (A. lancea, A. chinensis, A. koreana, and A. japonica) according to cytonuclear genomic data; this may be caused by interspecific hybridization or gene introgression.

Mitochondrial genomic integrity and the nuclear epigenome in health and disease.

Bidirectional crosstalk between the nuclear and mitochondrial genomes is essential for proper cell functioning. Mitochondrial DNA copy number (mtDNA-CN) and heteroplasmy influence mitochondrial function, which can influence the nuclear genome and contribute to health and disease. Evidence shows that mtDNA-CN and heteroplasmic variation are associated with aging, complex disease, and all-cause mortality. Further, the nuclear epigenome may mediate the effects of mtDNA variation on disease. In this way, mitochondria act as an environmental biosensor translating vital information about the state of the cell to the nuclear genome. Cellular communication between mtDNA variation and the nuclear epigenome can be achieved by modification of metabolites and intermediates of the citric acid cycle and oxidative phosphorylation. These essential molecules (e.g. ATP, acetyl-CoA, ɑ-ketoglutarate and S-adenosylmethionine) act as substrates and cofactors for enzymes involved in epigenetic modifications. The role of mitochondria as an environmental biosensor is emerging as a critical modifier of disease states. Uncovering the mechanisms of these dynamics in disease processes is expected to lead to earlier and improved treatment for a variety of diseases. However, the influence of mtDNA-CN and heteroplasmy variation on mitochondrially-derived epigenome-modifying metabolites and intermediates is poorly understood. This perspective will focus on the relationship between mtDNA-CN, heteroplasmy, and epigenome modifying cofactors and substrates, and the influence of their dynamics on the nuclear epigenome in health and disease.

Implications of mitochondrial DNA variants in pediatric B-cell acute lymphoblastic leukemia

BackgroundResearch on the role of variations in the mitochondrial genome in pathogenesis of acute lymphoblastic leukemia (ALL) has been unfolding at a rapid rate. Our laboratory has previously described higher number of copies of the mitochondrial genomes per cell in pediatric ALL patients as compared to the healthy controls. In the current study, we evaluated the pattern of mitochondrial genome variations in 20 de-novo pediatric B-ALL cases and seven controls. Quantitative real-time Polymerase Chain Reaction was used for estimation of mitochondrial genomes’ copy number in bone marrow samples of each ALL patient and peripheral blood samples of controls. The complete mitochondrial genomes of all samples were sequenced using the Illumina platform.ResultsSequencing data analysis using multiple mitochondrial genome databases revealed 325 variants in all 27 samples, out of which 221 variants were previously known while 104 were unassigned, new variants. The 325 variants consisted of 7 loss-of-function variants, 131 synonymous variants, 75 missense variants, and 112 non-coding variants. New, missense variants (n = 21) were identified in genes encoding the electron transport chain complexes with most of them encoding ND4, ND5 of complex I. Missense and loss-of-function variants were found to be deleterious by many predictor databases of pathogenicity. MuTect2 identified true somatic variants present only in tumors between patient-sibling pairs and showed overlap with missense and loss-of-function variants. Online MtDNA-server showed heteroplasmic and homoplasmic variants in mitochondrial genome.ConclusionsThe data suggest that some of these variations might have a deleterious impact on the expression of mitochondrial encoded genes with a possible functional relevance in leukemia.

A Novel MTTK Gene Variant m.8315A&gt;C as a Cause of MERRF Syndrome.

In this study, we report on a novel heteroplasmic pathogenic variant in mitochondrial DNA (mtDNA). The studied patient had myoclonus, epilepsy, muscle weakness, and hearing impairment and harbored a heteroplasmic m.8315A>C variant in the MTTK gene with a mutation load ranging from 71% to >96% in tested tissues. In muscle mitochondria, markedly decreased activities of respiratory chain complex I + III and complex IV were observed together with mildly reduced amounts of complex I and complex V (with the detection of V*- and free F1-subcomplexes) and a diminished level of complex IV holoenzyme. This pattern was previously seen in other MTTK pathogenic variants. The novel variant was not present in internal and publicly available control databases. Our report further expands the spectrum of MTTK variants associated with mitochondrial encephalopathies in adults.

Heteroplasmic mitochondrial DNA mutations in frontotemporal lobar degeneration

Frontotemporal lobar degeneration (FTLD) is a common cause of young onset dementia and is characterised by focal neuropathology. The reasons for the regional neuronal vulnerability are not known. Mitochondrial mechanisms have been implicated in the pathogenesis of FTLD, raising the possibility that frontotemporal regional mutations of mitochondrial DNA (mtDNA) are contributory causes. Here we applied dual sequencing of the entire mtDNA at high depth to identify high-fidelity single nucleotide variants (mtSNVs) and mtDNA rearrangements in post mortem brain tissue of people affected by FTLD and age-matched controls. Both mtSNVs and mtDNA rearrangements were elevated in the temporal lobe, with the greatest burden seen in FTLD. mtSNVs found in multiple brain regions also reached a higher heteroplasmy levels in the temporal lobe. The temporal lobe of people with FTLD had a higher burden of ribosomal gene variants predicted to affect intra-mitochondrial protein synthesis, and a higher proportion of missense variants in genes coding for respiratory chain subunits. In conclusion, heteroplasmic mtDNA variants predicted to affect oxidative phosphorylation are enriched in FTLD temporal lobe, and thus may contribute to the regional vulnerability in pathogenesis.

Identification of Somatic Mitochondrial DNA Mutations, Heteroplasmy, and Increased Levels of Catenanes in Tumor Specimens Obtained from Three Endometrial Cancer Patients.

Endometrial carcinoma (EC) is the most common type of gynecologic malignant epithelial tumor, with the death rate from this disease doubling over the past 20 years. Mitochondria provide cancer cells with necessary anabolic building blocks such as amino acids, lipids, and nucleotides, and EC samples have been shown to increase mitochondrial biogenesis. In cancer, mitochondrial DNA (mtDNA) heteroplasmy studies suggest that heteroplasmic variants encode predicted pathogenic proteins. We investigated the mtDNA genotypes within peri-normal and tumor specimens obtained from three individuals diagnosed with EC. DNA extracts from peri-normal and tumor tissues were used for mtDNA-specific next-generation sequencing and analyses of mtDNA content and topoisomers. The three tumors harbor heteroplasmic somatic mutations, and at least one mutation in each carcinoma is predicted to deleteriously alter a mtDNA-encoded protein. Somatic heteroplasmy linked to two mtDNA tRNA genes was found in separate tumors, and two heteroplasmic non-coding variants were identified in a single EC tumor. While two tumors had altered mtDNA content, all three displayed increased mtDNA catenanes. Our findings support that EC cells require wild-type mtDNA, but heteroplasmic mutations may alter mitochondrial metabolism to help promote cancer cell growth and proliferation.

Neither cardiac mitochondrial DNA variation nor copy number contribute to congenital heart disease risk

The well-established manifestation of mitochondrial mutations in functional cardiac disease (e.g., mitochondrial cardiomyopathy) prompted the hypothesis that mitochondrial DNA (mtDNA) sequence and/or copy number (mtDNAcn) variation contribute to cardiac defects in congenital heart disease (CHD). MtDNAcns were calculated and rare, non-synonymous mtDNA mutations were identified in 1,837 CHD-affected proband-parent trios, 116 CHD-affected singletons, and 114 paired cardiovascular tissue/blood samples. The variant allele fraction (VAF) of heteroplasmic variants in mitochondrial RNA from 257 CHD cardiovascular tissue samples was also calculated. On average, mtDNA from blood had 0.14 rare variants and 52.9 mtDNA copies per nuclear genome per proband. No variation with parental age at proband birth or CHD-affected proband age was seen. mtDNAcns in valve/vessel tissue (320 ± 70) were lower than in atrial tissue (1,080 ± 320, p = 6.8E−21), which were lower than in ventricle tissue (1,340 ± 280, p = 1.4E−4). The frequency of rare variants in CHD-affected individual DNA was indistinguishable from the frequency in an unaffected cohort, and proband mtDNAcns did not vary from those of CHD cohort parents. In both the CHD and the comparison cohorts, mtDNAcns were significantly correlated between mother-child, father-child, and mother-father. mtDNAcns among people with European (mean = 52.0), African (53.0), and Asian haplogroups (53.5) were calculated and were significantly different for European and Asian haplogroups (p = 2.6E−3). Variant heteroplasmic fraction (HF) in blood correlated well with paired cardiovascular tissue HF (r = 0.975) and RNA VAF (r = 0.953), which suggests blood HF is a reasonable proxy for HF in heart tissue. We conclude that mtDNA mutations and mtDNAcns are unlikely to contribute significantly to CHD risk.

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 bioinformatics pipeline for estimating mitochondrial DNA copy number and heteroplasmy levels from whole genome sequencing data.

Mitochondrial diseases are a heterogeneous group of disorders that can be caused by mutations in the nuclear or mitochondrial genome. Mitochondrial DNA (mtDNA) variants may exist in a state of heteroplasmy, where a percentage of DNA molecules harbor a variant, or homoplasmy, where all DNA molecules have the same variant. The relative quantity of mtDNA in a cell, or copy number (mtDNA-CN), is associated with mitochondrial function, human disease, and mortality. To facilitate accurate identification of heteroplasmy and quantify mtDNA-CN, we built a bioinformatics pipeline that takes whole genome sequencing data and outputs mitochondrial variants, and mtDNA-CN. We incorporate variant annotations to facilitate determination of variant significance. Our pipeline yields uniform coverage by remapping to a circularized chrM and by recovering reads falsely mapped to nuclear-encoded mitochondrial sequences. Notably, we construct a consensus chrM sequence for each sample and recall heteroplasmy against the sample's unique mitochondrial genome. We observe an approximately 3-fold increased association with age for heteroplasmic variants in non-homopolymer regions and, are better able to capture genetic variation in the D-loop of chrM compared to existing software. Our bioinformatics pipeline more accurately captures features of mitochondrial genetics than existing pipelines that are important in understanding how mitochondrial dysfunction contributes to disease.

Benchmarking the Effectiveness and Accuracy of Multiple Mitochondrial DNA Variant Callers: Practical Implications for Clinical Application.

Mitochondrial DNA (mtDNA) mutations contribute to human disease across a range of severity, from rare, highly penetrant mutations causal for monogenic disorders to mutations with milder contributions to phenotypes. mtDNA variation can exist in all copies of mtDNA or in a percentage of mtDNA copies and can be detected with levels as low as 1%. The large number of copies of mtDNA and the possibility of multiple alternative alleles at the same DNA nucleotide position make the task of identifying allelic variation in mtDNA very challenging. In recent years, specialized variant calling algorithms have been developed that are tailored to identify mtDNA variation from whole-genome sequencing (WGS) data. However, very few studies have systematically evaluated and compared these methods for the detection of both homoplasmy and heteroplasmy. A publicly available synthetic gold standard dataset was used to assess four mtDNA variant callers (Mutserve, mitoCaller, MitoSeek, and MToolBox), and the commonly used Genome Analysis Toolkit “best practices” pipeline, which is included in most current WGS pipelines. We also used WGS data from 126 trios and calculated the percentage of maternally inherited variants as a metric of calling accuracy, especially for homoplasmic variants. We additionally compared multiple pathogenicity prediction resources for mtDNA variants. Although the accuracy of homoplasmic variant detection was high for the majority of the callers with high concordance across callers, we found a very low concordance rate between mtDNA variant callers for heteroplasmic variants ranging from 2.8% to 3.6%, for heteroplasmy thresholds of 5% and 1%. Overall, Mutserve showed the best performance using the synthetic benchmark dataset. The analysis of mtDNA pathogenicity resources also showed low concordance in prediction results. We have shown that while homoplasmic variant calling is consistent between callers, there remains a significant discrepancy in heteroplasmic variant calling. We found that resources like population frequency databases and pathogenicity predictors are now available for variant annotation but still need refinement and improvement. With its peculiarities, the mitochondria require special considerations, and we advocate that caution needs to be taken when analyzing mtDNA data from WGS data.

Mitochondrial DNA content reduction in the most fertile spermatozoa is accompanied by increased mitochondrial DNA rearrangement.

Is there an association between male fertility and spermatozoa mitochondrial DNA (mtDNA) copy number and genome rearrangements? Normal spermatozoa not only have a lower mtDNA copy number but also more DNA rearrangements than spermatozoa of men with severe oligoasthenospermia (SOA). While there is a consensus that mtDNA content is decreased in the most fertile spermatozoa, the role of mtDNA sequence alteration in male infertility is unclear. High-throughput sequencing, which allows an exhaustive analysis of mtDNA rearrangements and mutations, could be helpful in this context, but has yet to be used. This is an observational study of semen samples obtained from 44 men undergoing ART at an academic infertility centre in France, from October 2018 to November 2020. The men were classified into two groups: 20 men in the SOA group and 24 men with normal semen parameters in the control group. For each patient and control, mtDNA was isolated from sperm fractions from the 40% and 90% layers of the density gradient. The average mtDNA content of each sample was assessed using digital PCR. Deep sequencing was performed using next-generation sequencing. Signal processing and base calling were performed via the embedded pre-processing pipeline, the variants were analysed using an in-house workflow and a dedicated tool, based on soft-clipping, was used to study large mtDNA rearrangements. The distribution and the type of rearrangements and variants were compared between patients with SOA and controls on one hand, and between the 40% and 90% gradient layers, on the other hand. The mtDNA content of spermatozoa in the SOA group was significantly higher than in the control group (P < 0.0001). Moreover, mtDNA content was significantly higher in spermatozoa from the 40% layer (the most fertile spermatozoa) compared to the 90% layer, both in the SOA (P = 0.02) and the control group (P < 0.0001). The frequency of large mtDNA deletions and duplications was significantly higher in the control group (P = 0.002). Most of these rearrangements are potentially related to DNA breaks and their number was reduced by the removal of the linear mtDNA from the samples. Heteroplasmic variants were found more frequently in the SOA group (P = 0.05) and in the 40% layer (P = 0.03), but none had any obvious functional consequence. Our findings are novel and significant but should be verified in larger cohorts and other types of male infertility. Our findings suggest that sperm mtDNA rearrangements are not necessarily associated with mitochondrial dysfunction and male infertility. Instead, they seem to be concomitant with the process of mtDNA content reduction in the most potentially fertile spermatozoa. Furthermore, they refute the hypothesis that, in the case of mtDNA alteration, a compensatory mechanism allows an increase in mtDNA copy number to rectify the energy deficit. The increased frequency of mtDNA rearrangements in the most fertile spermatozoa is a novel result that offers new insight into the relation between sperm quality and mtDNA. This work was supported by Angers University Hospital (grant AOI CHU Angers 2018), Angers University and the French national research centres INSERM and CNRS. There are no competing interests. N/A.

Comparison of mitochondrial DNA sequences from whole blood and lymphoblastoid cell lines

Lymphoblastoid cell lines (LCLs) provide an unlimited source of genomic DNA for genetic studies. Here, we compared mtDNA sequence variants, heteroplasmic or homplasmic, between LCL (sequenced by mitoRCA-seq method) and whole blood samples (sequenced through whole genome sequencing approach) of the same 130 participants in the Framingham Heart Study. We applied harmonization of sequence coverages and consistent quality control to mtDNA sequences. We identified 866 variation sites in the 130 LCL samples and 666 sites in the 130 blood samples. More than 94% of the identified homoplasmies were present in both LCL and blood samples while more than 70% of heteroplasmic sites were uniquely present either in LCL or in blood samples. The LCL and whole blood samples carried a similar number of homoplasmic variants (p = 0.45) per sample while the LCL carried a greater number of heteroplasmic variants than whole blood per sample (p < 2.2e−16). Furthermore, the LCL samples tended to accumulate low level heteroplasmies (heteroplasmy level in 3–25%) than their paired blood samples (p = 0.001). These results suggest that cautions should be taken in the interpretation and comparison of findings when different tissues/cell types or different sequencing technologies are applied to obtain mtDNA sequences.

Comprehensive summary of mitochondrial DNA alterations in the postmortem human brain: A systematic review.

SummaryBackgroundMitochondrial DNA (mtDNA) encodes 37 genes necessary for synthesizing 13 essential subunits of the oxidative phosphorylation system. mtDNA alterations are known to cause mitochondrial disease (MitD), a clinically heterogeneous group of disorders that often present with neuropsychiatric symptoms. Understanding the nature and frequency of mtDNA alterations in health and disease could be a cornerstone in disentangling the relationship between biochemical findings and clinical symptoms of brain disorders. This systematic review aimed to summarize the mtDNA alterations in human brain tissue reported to date that have implications for further research on the pathophysiological significance of mtDNA alterations in brain functioning.MethodsWe searched the PubMed and Embase databases using distinct terms related to postmortem human brain and mtDNA up to June 10, 2021. Reports were eligible if they were empirical studies analysing mtDNA in postmortem human brains.FindingsA total of 158 of 637 studies fulfilled the inclusion criteria and were clustered into the following groups: MitD (48 entries), neurological diseases (NeuD, 55 entries), psychiatric diseases (PsyD, 15 entries), a miscellaneous group with controls and other clinical diseases (5 entries), ageing (20 entries), and technical issues (5 entries). Ten entries were ascribed to more than one group. Pathogenic single nucleotide variants (pSNVs), both homo- or heteroplasmic variants, have been widely reported in MitD, with heteroplasmy levels varying among brain regions; however, pSNVs are rarer in NeuD, PsyD and ageing. A lower mtDNA copy number (CN) in disease was described in most, but not all, of the identified studies. mtDNA deletions were identified in individuals in the four clinical categories and ageing. Notably, brain samples showed significantly more mtDNA deletions and at higher heteroplasmy percentages than blood samples, and several of the deletions present in the brain were not detected in the blood. Finally, mtDNA heteroplasmy, mtDNA CN and the deletion levels varied depending on the brain region studied.InterpretationmtDNA alterations are well known to affect human tissues, including the brain. In general, we found that studies of MitD, NeuD, PsyD, and ageing were highly variable in terms of the type of disease or ageing process investigated, number of screened individuals, studied brain regions and technology used. In NeuD and PsyD, no particular type of mtDNA alteration could be unequivocally assigned to any specific disease or diagnostic group. However, the presence of mtDNA deletions and mtDNA CN variation imply a role for mtDNA in NeuD and PsyD. Heteroplasmy levels and threshold effects, affected brain regions, and mitotic segregation patterns of mtDNA alterations may be involved in the complex inheritance of NeuD and PsyD and in the ageing process. Therefore, more information is needed regarding the type of mtDNA alteration, the affected brain regions, the heteroplasmy levels, and their relationship with clinical phenotypes and the ageing process.FundingHospital Universitari Institut Pere Mata; Institut d'Investigació Sanitària Pere Virgili; Instituto de Salud Carlos III, Ministerio de Ciencia e Innovación (PI18/00514).

Mitochondrial DNA variation across 56,434 individuals in gnomAD.

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

The minicircular and extremely heteroplasmic mitogenome of the holoparasitic plant Rhopalocnemis phalloides

The plastid and nuclear genomes of parasitic plants exhibit deeply altered architectures,1-13 whereas the few examined mitogenomes range from deeply altered to conventional.14-20 To provide further insight on mitogenome evolution in parasitic plants, we report the highly modified mitogenome of Rhopalocnemis phalloides, a holoparasite in Balanophoraceae. Its mitogenome is uniquely arranged in 21minicircular chromosomes that vary in size from 4,949 to 7,861bp, with a total length of only 130,713bp. All chromosomes share an identical 896bp conserved region, with a large stem-loop that acts as the origin of replication, flanked on each side by hypervariable and semi-conserved regions. Similar minicircular structures with shared and unique regions have been observed in parasitic animals and free-living protists,21-24 suggesting convergent structural evolution. Southern blots confirm both the minicircular structure and the replication origin of the mitochondrial chromosomes. PacBio reads provide evidence for chromosome recombination and rolling-circle replication for the R.phalloides mitogenome. Despite its small size, the mitogenome harbors a typical set of genes and introns within the unique regions of each chromosome, yet introns are the smallest among seed plants and ferns. The mitogenome also exhibits extreme heteroplasmy, predominantly involving short indels and more complex variants, many of which cause potential loss-of-function mutations for some gene copies. All heteroplasmic variants are transcribed, and functional and nonfunctional protein-coding variants are spliced and RNA edited. Our findings offer a unique perspective into how mitogenomes of parasitic plants can be deeply altered and shed light on plant mitogenome replication.