HmtVar: a new resource for human mitochondrial variations and pathogenicity data

The genomic mosaic of mitochondrial dysfunction: Decoding nuclear and mitochondrial epigenetic contributions to maternally inherited diabetes and deafness pathogenesis

Aims

Maternally inherited diabetes and deafness (MIDD) is a complex disorder characterized by multiorgan clinical manifestations, including diabetes, hearing loss, and ophthalmic complications. This pilot study aimed to elucidate the intricate interplay between nuclear and mitochondrial genetics, epigenetic modifications, and their potential implications in the pathogenesis of MIDD.

Main methods

A comprehensive genomic approach was employed to analyze a Sicilian family affected by clinically characterized MIDD, negative to the only known causative m.3243 A > G variant, integrating whole-exome sequencing and whole-genome bisulfite sequencing of both nuclear and mitochondrial analyses.

Key findings

Rare and deleterious variants were identified across multiple nuclear genes involved in retinal homeostasis, mitochondrial function, and epigenetic regulation, while complementary mitochondrial DNA analysis revealed a rich tapestry of genetic diversity across genes encoding components of the electron transport chain and ATP synthesis machinery. Epigenetic analyses uncovered significant differentially methylated regions across the genome and within the mitochondrial genome, suggesting a nuanced landscape of epigenetic modulation.

Significance

The integration of genetic and epigenetic data highlighted the potential crosstalk between nuclear and mitochondrial regulation, with specific mtDNA variants influencing methylation patterns and potentially impacting the expression and regulation of mitochondrial genes. This pilot study provides valuable insights into the complex molecular mechanisms underlying MIDD, emphasizing the interplay between nucleus and mitochondrion, tracing the way for future research into targeted therapeutic interventions and personalized approaches for disease management.

Association between mitochondrial DNA genotype and sperm motility in humans

The relationship between genetic alterations in mitochondrial DNA (mtDNA) and progressive motility (PR) and rapid progressive motility (grade A) of ejaculated human spermatozoa remains unclear. In this study, we explored the association between human mtDNA genotype and sperm PR and grade A by analyzing mtDNA copy number, loci, haplogroup, rearrangement, deletions, and duplications and sperm motility parameters. Human sperm mtDNA copy number, loci and haplogroups were not associated with human sperm motility PR or A grade. However, the cumulative frequency of human sperm mtDNA rearrangements (including deletions and duplications) in participants with high PR and grade A ratio was higher than in participants with low PR and grade A ratio. Additional studies are needed to understand the relationship between mtDNA genotypes, including deletions and duplications, and human sperm motility.

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.

An in silico approach to develop potential therapies against Middle East Respiratory Syndrome Coronavirus (MERS-CoV)

A deadly respiratory disease Middle East Respiratory Syndrome (MERS) is caused by a perilous virus known as MERS-CoV, which has a severe impact on human health. Currently, there is no approved vaccine, prophylaxis, or antiviral therapeutics for preventing MERS-CoV infection. Due to its inexorable and integral role in the maturation and replication of the MERS-CoV virus, the 3C-like protease is unavoidly a viable therapeutic target. In this study, 2369 phytoconstituents were enlisted from Japanese medicinal plants, and these compounds were screened against 3C-like protease to identify feasible inhibitors. The best three compounds were identified as Kihadanin B, Robustaflavone, and 3-beta-O- (trans-p-Coumaroyl) maslinic acid, with binding energies of -9.8, -9.4, and -9.2 kcal/mol, respectively. The top three potential candidates interacted with several active site residues in the targeted protein, including Cys145, Met168, Glu169, Ala171, and Gln192. The best three compounds were assessed by in silico technique to determine their drug-likeness properties, and they exhibited the least harmful features and the greatest drug-like qualities. Various descriptors, such as solvent-accessible surface area, root-mean-square fluctuation, root-mean-square deviation, hydrogen bond, and radius of gyration, validated the stability and firmness of the protein-ligand complexes throughout the 100ns molecular dynamics simulation. Moreover, the top three compounds exhibited better binding energy along with better stability and firmness than the inhibitor (Nafamostat), which was further confirmed by the binding free energy calculation. Therefore, this computational investigation could aid in the development of efficient therapeutics for life-threatening MERS-CoV infections.

MmisAT and MmisP: an efficient and accurate suite of variant analysis toolkit for primary mitochondrial diseases

Recent advances in next-generation sequencing (NGS) technology have greatly accelerated the need for efficient annotation to accurately interpret clinically relevant genetic variants in human diseases. Therefore, it is crucial to develop appropriate analytical tools to improve the interpretation of disease variants. Given the unique genetic characteristics of mitochondria, including haplogroup, heteroplasmy, and maternal inheritance, we developed a suite of variant analysis toolkits specifically designed for primary mitochondrial diseases: the Mitochondrial Missense Variant Annotation Tool (MmisAT) and the Mitochondrial Missense Variant Pathogenicity Predictor (MmisP). MmisAT can handle protein-coding variants from both nuclear DNA and mtDNA and generate 349 annotation types across six categories. It processes 4.78 million variant data in 76 min, making it a valuable resource for clinical and research applications. Additionally, MmisP provides pathogenicity scores to predict the pathogenicity of genetic variations in mitochondrial disease. It has been validated using cross-validation and external datasets and demonstrated higher overall discriminant accuracy with a receiver operating characteristic (ROC) curve area under the curve (AUC) of 0.94, outperforming existing pathogenicity predictors. In conclusion, the MmisAT is an efficient tool that greatly facilitates the process of variant annotation, expanding the scope of variant annotation information. Furthermore, the development of MmisP provides valuable insights into the creation of disease-specific, phenotype-specific, and even gene-specific predictors of pathogenicity, further advancing our understanding of specific fields.

Effects of gamma-radiation on microbial, nutritional, and functional properties of Katimon mango peels: A combined biochemical and in silico studies

Gamma radiation has notable impacts on the flesh of mangoes. In this research, Katimon mangoes were subjected to different levels of irradiation (0.5, 1.0, 1.5, and 2.0 kGy) using a60Co irradiator. The results showed that irradiation significantly reduced the microbial population in the mango peels, with the 1.5 kGy dose showing the most significant reduction. Irradiation also delayed ripening and extended the shelf life of the mango peels. The total fat, protein, ash, moisture, and sugar content of the mango peels were all affected by irradiation. The total protein content, ash content and moisture content increased after irradiation, while the fat content remained relatively unchanged. The sugar content increased in all samples after storage, but the non-irradiated samples had higher sugar levels than the irradiated ones. The dietary fiber content of the mango peels was not significantly affected by irradiation. The vitamin C content decreased in all samples after storage. The titratable acidity and total soluble solids content of the mango peels increased after storage, but there were no significant differences between the irradiated and non-irradiated samples. Antioxidant activity and cytotoxicity assessment highlighted the antioxidant potential and reduced toxicity of irradiated samples. Additionally, the antimicrobial effectiveness of irradiated mango peels was evaluated. The most substantial inhibitory zones (measuring 16.90 ± 0.35) against Pseudomonas sp. were observed at a radiation dose of 1.5 kGy with 150 μg/disc. To identify potential antimicrobial agents, the volatile components of mangoes irradiated with 1.5 kGy were analyzed through GC-MS. Subsequently, these compounds were subjected to in silico studies against a viable protein, TgpA, of Pseudomonas sp. (PDB ID: 6G49). Based on molecular dynamic simulations and ADMET properties, (-)-Carvone (-6.2), p-Cymene (-6.1), and Acetic acid phenylmethyl ester (-6.1) were identified as promising compounds for controlling Pseudomonas sp.

Wide diagnostic and genotypic spectrum in patients with suspected mitochondrial disease

BACKGROUND

Mitochondrial Diseases (MDs) are a diverse group of neurometabolic disorders characterized by impaired mitochondrial oxidative phosphorylation and caused by pathogenic variants in more than 400 genes. The implementation of next-generation sequencing (NGS) technologies helps to increase the understanding of molecular basis and diagnostic yield of these conditions. The purpose of the study was to investigate diagnostic and genotypic spectrum in patients with suspected MD. The comprehensive analysis of mtDNA variants using Sanger sequencing was performed in the group of 83 unrelated individuals with clinically suspected mitochondrial disease. Additionally, targeted next generation sequencing or whole exome sequencing (WES) was performed for 30 patients of the study group.

RESULTS

The overall diagnostic rate was 21.7% for the patients with suspected MD, increasing to 36.7% in the group of patients where NGS methods were applied. Mitochondrial disease was confirmed in 11 patients (13.3%), including few classical mitochondrial syndromes (MELAS, MERRF, Leigh and Kearns-Sayre syndrome) caused by pathogenic mtDNA variants (8.4%) and MDs caused by pathogenic variants in five nDNA genes. Other neuromuscular diseases caused by pathogenic variants in seven nDNA genes, were confirmed in seven patients (23.3%).

CONCLUSION

The wide spectrum of identified rare mitochondrial or neurodevelopmental diseases proves that MD suspected patients would mostly benefit from an extensive genetic profiling allowing rapid diagnostics and improving the care of these patients.

Modeling-based prediction tools for preimplantation genetic testing of mitochondrial DNA diseases: estimating symptomatic thresholds, risk, and chance of success

PURPOSE

Preimplantation genetic testing (PGT) has become a reliable tool for preventing the germline transmission of mitochondrial DNA (mtDNA) variants. However, procedures are not standardized across mtDNA variants. In this study, we aim to estimate symptomatic thresholds, risk, and chance of success for PGT for mtDNA pathogenic variant carriers.

METHODS

We performed a systematic analysis of heteroplasmy data including 455 individuals from 187 familial pedigrees with the common m.3243A>G, m.8344A>G, or m.8993T>G pathogenic variants. We applied binary logistic regression for estimating symptomatic thresholds of heteroplasmy, simplified Sewell-Wright formula and Kimura equations for predicting the risk of disease transmission, and binomial distribution for predicting minimum oocyte numbers.

RESULTS

We estimated the symptomatic thresholds of m.8993T>G and m.8344A>G as 29.86% and 16.15%, respectively. We could not determine a threshold for m.3243A>G. We established models for mothers harboring common and rare mtDNA pathogenic variants to predict the risk of disease transmission and the number of oocytes required to produce an embryo with sufficiently low variant load. In addition, we provide a table allowing the prediction of transmission risk and the minimum required oocytes for PGT patients with different variant levels.

CONCLUSION

We have established models that can determine the symptomatic thresholds of common mtDNA pathogenic variants. We also constructed universal models applicable to nearly all mtDNA pathogenic variants which can predict risk and minimum numbers for PGT patients. These models have advanced our understanding of mtDNA disease pathogenesis and will enable more effective prevention of disease transmission using PGT.

APOGEE 2: multi-layer machine-learning model for the interpretable prediction of mitochondrial missense variants

Mitochondrial dysfunction has pleiotropic effects and is frequently caused by mitochondrial DNA mutations. However, factors such as significant variability in clinical manifestations make interpreting the pathogenicity of variants in the mitochondrial genome challenging. Here, we present APOGEE 2, a mitochondrially-centered ensemble method designed to improve the accuracy of pathogenicity predictions for interpreting missense mitochondrial variants. Built on the joint consensus recommendations by the American College of Medical Genetics and Genomics/Association for Molecular Pathology, APOGEE 2 features an improved machine learning method and a curated training set for enhanced performance metrics. It offers region-wise assessments of genome fragility and mechanistic analyses of specific amino acids that cause perceptible long-range effects on protein structure. With clinical and research use in mind, APOGEE 2 scores and pathogenicity probabilities are precompiled and available in MitImpact. APOGEE 2's ability to address challenges in interpreting mitochondrial missense variants makes it an essential tool in the field of mitochondrial genetics.

Mitochondrial DNA in Human Diversity and Health: From the Golden Age to the Omics Era

Mitochondrial DNA (mtDNA) is a small fraction of our hereditary material. However, this molecule has had an overwhelming presence in scientific research for decades until the arrival of high-throughput studies. Several appealing properties justify the application of mtDNA to understand how human populations are-from a genetic perspective-and how individuals exhibit phenotypes of biomedical importance. Here, I review the basics of mitochondrial studies with a focus on the dawn of the field, analysis methods and the connection between two sides of mitochondrial genetics: anthropological and biomedical. The particularities of mtDNA, with respect to inheritance pattern, evolutionary rate and dependence on the nuclear genome, explain the challenges of associating mtDNA composition and diseases. Finally, I consider the relevance of this single locus in the context of omics research. The present work may serve as a tribute to a tool that has provided important insights into the past and present of humankind.

Assessment of Amikacin- and Capreomycin-Related Adverse Drug Reactions in Patients with Multidrug-Resistant Tuberculosis and Exploring the Role of Genetic Factors

Following the introduction of all-oral treatment regimens for patients with drug-resistant tuberculosis (TB), second-line injectable drug applications have been reduced in the last few years. However, they are still important for anti-TB therapy. This study aims to analyze the occurrence of amikacin- and capreomycin-related adverse drug reactions (ADR) in patients with multidrug-resistant tuberculosis (MDR-TB) and evaluate the role of multiple patient-, disease-, and therapy-related factors on the frequency of the observed adverse events. In addition, the possible role of genetic risk factors was studied by full-length mitochondrial DNA sequencing. Toward this aim, we retrospectively evaluated 47 patients with MDR-TB who received amikacin and/or capreomycin. In total, 16 (34.0%) patients developed ototoxicity and 13 (27.7%) developed nephrotoxicity, including 3 (6.4%) patients who experienced both adverse events. Ototoxicity development was more common in patients who received amikacin. No other factors showed a significant impact. Nephrotoxicity was likely associated with previous renal health impairment. Full mitochondrial genome sequencing did not reveal any specific ADR-associated variants, and results showed no differences in adverse event occurrence for any specific variants, mutation count, or mitochondrial haplogroup. The absence of the previously reported ototoxicity-related mtDNA variants in our patients with ototoxicity and nephrotoxicity highlighted the complex nature of the ADR occurrence.

Mitochondrial DNA variants in a cohort from Argentina with suspected Leber's hereditary optic neuropathy (LHON)

The present study investigates the spectrum and analysis of mitochondrial DNA (mtDNA) variants associated with Leber hereditary optic neuropathy (LHON) in an Argentinean cohort, analyzing 3 LHON-associated mitochondrial genes. In 32% of the cases, molecular confirmation of the diagnosis could be established, due to the identification of disease-causing variants. A total of 54 variants were observed in a cohort of 100 patients tested with direct sequencing analysis. The frequent causative mutations m.11778G>A in MT-ND4, m.3460G>A in MT-ND1, and m.14484T>C in MT-ND6 were identified in 28% of the cases of our cohort. Secondary mutations in this Argentinean LHON cohort were m.11253T>C p.Ile165Thr in MT-ND4, identified in three patients (3/100, 3%) and m.3395A>G p.Tyr30Cys in MT-ND1, in one of the patients studied (1%). This study shows, for the first time, the analysis of mtDNA variants in patients with a probable diagnosis of LHON in Argentina. Standard molecular methods are an effective first approach in order to achieve genetic diagnosis of the disease, leaving NGS tests for those patients with negative results.

Mitochondrial genome in sporadic breast cancer: A case control study and a proteomic analysis in a Sinhalese cohort from Sri Lanka

Breast cancer is the commonest malignancy in women and the majority occurs sporadically with no hereditary predisposition. However, sporadic breast cancer has been studied less intensively than the hereditary form and to date hardly any predictive biomarkers exist for the former. Furthermore, although mitochondrial DNA variants have been reported to be associated with breast cancer, findings have been inconsistent across populations. Thus we carried out a case control study on sporadic breast cancer patients and healthy controls of Sinhalese ethnicity (N = 60 matched pairs) in order to characterize coding region variants associated with the disease and to identify any potential biomarkers. Mitochondrial genome was fully sequenced in 30 pairs and selected regions were sequenced in the remaining 30 pairs. Several in-silico tools were used to assess functional significance of the variants observed. A number of variants were identified among the patients and the controls. Missense variants identified were either polymorphisms or rare variants. Their prevalence did not significantly differ between patients and the healthy controls (matched for age, body mass index and menopausal status). MT-CYB, MT-ATP6 and MT-ND2 genes showed a higher mutation rate. A higher proportion of pre-menopausal patients carried missense and pathogenic variants. Unique combinations of missense variants were seen within genes and these occurred mostly in MT-ATP6 and MT-CYB genes. Such unique combinations that occurred exclusively among the patients were common in obese patients. Mitochondrial DNA variants may have a role in breast carcinogenesis in obesity and pre-menopause. Molecular dynamic simulations suggested the mutants, G78S in MT-CO3 gene and T146A in MT-ATP6 gene are likely to be more stable than their wild type counterparts.

Spectrum of germline and somatic mitochondrial DNA variants in Tuberous Sclerosis Complex

Tuberous Sclerosis Complex (TSC) is caused by loss of function variants in either TSC1 or TSC2 and is characterized by broad phenotypic heterogeneity. Currently, there is limited knowledge regarding the role of the mitochondrial genome (mtDNA) in TSC pathogenesis. In this study, we aimed to determine the prevalence and spectrum of germline and somatic mtDNA variants in TSC and identify potential disease modifiers. Analysis of mtDNA amplicon massively parallel sequencing (aMPS) data, off-target mtDNA from whole-exome sequencing (WES), and/or qPCR, revealed mtDNA alterations in 270 diverse tissues (139 TSC-associated tumors and 131 normal tissue samples) from 199 patients and six healthy individuals. Correlation of clinical features to mtDNA variants and haplogroup analysis was done in 102 buccal swabs (age: 20-71 years). No correlation was found between clinical features and either mtDNA variants or haplogroups. No pathogenic variants were identified in the buccal swab samples. Using in silico analysis, we identified three predicted pathogenic variants in tumor samples: MT-ND4 (m.11742G>A, p. Cys328Tyr, VAF: 43%, kidney angiomyolipoma), MT-CYB (m.14775T>C, p. Leu10Pro, VAF: 43%, LAM abdominal tumor) and MT-CYB (m.15555C>T, p. Pro270Leu, VAF: 7%, renal cell carcinoma). Large deletions of the mitochondrial genome were not detected. Analysis of tumors from 23 patients with corresponding normal tissue did not reveal any recurrent tumor-associated somatic variants. The mtDNA/gDNA ratio between tumors and corresponding normal tissue was also unchanged. Overall, our findings demonstrate that the mitochondrial genome is highly stable across tissues and within TSC-associated tumors.

MTSviewer: A database to visualize mitochondrial targeting sequences, cleavage sites, and mutations on protein structures

Mitochondrial dysfunction is implicated in a wide array of human diseases ranging from neurodegenerative disorders to cardiovascular defects. The coordinated localization and import of proteins into mitochondria are essential processes that ensure mitochondrial homeostasis. The localization and import of most mitochondrial proteins are driven by N-terminal mitochondrial targeting sequences (MTS's), which interact with import machinery and are removed by the mitochondrial processing peptidase (MPP). The recent discovery of internal MTS's-those which are distributed throughout a protein and act as import regulators or secondary MPP cleavage sites-has expanded the role of both MTS's and MPP beyond conventional N-terminal regulatory pathways. Still, the global mutational landscape of MTS's remains poorly characterized, both from genetic and structural perspectives. To this end, we have integrated a variety of tools into one harmonized R/Shiny database called MTSviewer (https://neurobioinfo.github.io/MTSvieweR/), which combines MTS predictions, cleavage sites, genetic variants, pathogenicity predictions, and N-terminomics data with structural visualization using AlphaFold models of human and yeast mitochondrial proteomes. Using MTSviewer, we profiled all MTS-containing proteins across human and yeast mitochondrial proteomes and provide multiple case studies to highlight the utility of this database.

Mitochondrial DNA copy number changes, heteroplasmy, and mutations in plasma-derived exosomes and brain tissue of glioblastoma patients

Glioblastoma is the most common malignant tumor of the central nervous system (CNS) in adults. Glioblastoma cells show increased glucose consumption associated with poor prognosis. Since mitochondria play a crucial role in energy metabolism, mutations and copy number changes of mitochondrial DNA may serve as biomarkers. As the brain is difficult to access, analysis of mitochondria directly from the brain tissue represents a challenge. Exosome analysis is an alternative (still poorly explored) approach to investigate molecular changes in CNS tumors. We analyzed brain tissue DNA and plasma-derived exosomal DNA (exoDNA) of 44 glioblastoma patients and 40 control individuals. Quantitative real-time PCR was performed to determine mtDNA copy numbers and the Kruskal-Wallis and Mann-Whitney U test were used for statistical analysis of data. Subsequently, sequencing libraries were prepared and sequenced on the MiSeq platform to identify mtDNA point mutations. Tissue mtDNA copy number was different among controls and patients in multiple comparisons. A similar tendency was detected in exosomes. Based on NGS analysis, several mtDNA point mutations showed slightly different frequencies between cases and controls, but the clinical relevance of these observations is difficult to assess and likely less than that of overall mtDNA copy number changes. Allele frequencies of variants were used to determine the level of heteroplasmy (found to be higher in exo-mtDNA of control individuals). Despite the suggested potential, the use of such biomarkers for the screening and/or diagnosis of glioblastomas is still limited, thus further studies are needed.

MicroRNA and Metabolic Profiling of a Primary Ovarian Neuroendocrine Carcinoma Pulmonary-Type Reveals a High Degree of Similarity with Small Cell Lung Cancer

Small cell neuroendocrine carcinoma is most frequently found in the lung (SCLC), but it has been also reported, albeit with a very low incidence, in the ovary. Here, we analyze a case of primary small cell carcinoma of the ovary of pulmonary type (SCCOPT), a rare and aggressive tumor with poor prognosis, whose biology and molecular features have not yet been thoroughly investigated. The patient affected by SCCOPT had a residual tumor following chemotherapy which displayed pronounced similarity with neuroendocrine tumors and lung cancer in terms of its microRNA expression profile and mTOR-downstream activation. By analyzing the metabolic markers of the neoplastic lesion, we established a likely glycolytic signature. In conclusion, this in-depth characterization of SCCOPT could be useful for future diagnoses, possibly aided by microRNA profiling, allowing clinicians to adopt the most appropriate therapeutic strategy.

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

Background

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

Results

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

Conclusions

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

MitoVisualize: a resource for analysis of variants in human mitochondrial RNAs and DNA

SUMMARY

We present MitoVisualize, a new tool for analysis of the human mitochondrial DNA (mtDNA). MitoVisualize enables visualization of: (i) the position and effect of variants in mitochondrial transfer RNA and ribosomal RNA secondary structures alongside curated variant annotations, (ii) data across RNA structures, such as to show all positions with disease-associated variants or with post-transcriptional modifications and (iii) the position of a base, gene or region in the circular mtDNA map, such as to show the location of a large deletion. All visualizations can be easily downloaded as figures for reuse. MitoVisualize can be useful for anyone interested in exploring mtDNA variation, though is designed to facilitate mtDNA variant interpretation in particular.

AVAILABILITY AND IMPLEMENTATION

MitoVisualize can be accessed via https://www.mitovisualize.org/. The source code is available at https://github.com/leklab/mito_visualize/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Variant pathogenic prediction by locus variability: the importance of the current picture of evolution

Accurate detection of pathogenic single nucleotide variants (SNVs) is a key challenge in whole exome and whole genome sequencing studies. To date, several in silico tools have been developed to predict deleterious variants from this type of data. However, these tools have limited power to detect new pathogenic variants, especially in non-coding regions. In this study, we evaluate the use of a new metric, the Shannon Entropy of Locus Variability (SELV), calculated as the Shannon entropy of the variant frequencies reported in genome-wide population studies at a given locus, as a new predictor of potentially pathogenic variants in non-coding nuclear and mitochondrial DNA and also in coding regions with a selective pressure other than that imposed by the genetic code, e.g splice-sites. For benchmarking, SELV was compared to predictors of pathogenicity in different genomic contexts. In nuclear non-coding DNA, SELV outperformed CDTS (AUC_SELV = 0.97 in ROC curve and PR-AUC_SELV = 0.96 in Precision-recall curve). For non-coding mitochondrial variants (AUC_SELV = 0.98 in ROC curve and PR-AUC_SELV = 1.00 in Precision-recall curve) SELV outperformed HmtVar. Moreover, SELV was compared against two state-of-the-art ensemble predictors of pathogenicity in splice-sites, ada-score, and rf-score, matching their overall performance both in ROC (AUC_SELV = 0.95) and Precision-recall curves (PR-AUC = 0.97), with the advantage that SELV can be easily calculated for every position in the genome, as opposite to ada-score and rf-score. Therefore, we suggest that the information about the observed genetic variability in a locus reported from large scale population studies could improve the prioritization of SNVs in splice-sites and in non-coding regions.

Best practices for the interpretation and reporting of clinical whole genome sequencing

Whole genome sequencing (WGS) shows promise as a first-tier diagnostic test for patients with rare genetic disorders. However, standards addressing the definition and deployment practice of a best-in-class test are lacking. To address these gaps, the Medical Genome Initiative, a consortium of leading health care and research organizations in the US and Canada, was formed to expand access to high quality clinical WGS by convening experts and publishing best practices. Here, we present best practice recommendations for the interpretation and reporting of clinical diagnostic WGS, including discussion of challenges and emerging approaches that will be critical to harness the full potential of this comprehensive test.

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.

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.

Lung Tumor Growth Promotion by Tobacco-Specific Nitrosamines Involves the β2-Adrenergic Receptors-Dependent Stimulation of Mitochondrial REDOX Signaling

Aims: Lung cancer is the leading cause of cancer death worldwide, and tobacco smoking is a recognized major risk factor for lung tumor development. We analyzed the effect of tobacco-specific nitrosamines (TSNAs) on human lung adenocarcinoma metabolic reprogramming, an emergent hallmark of carcinogenesis. Results: A series of in vitro and in vivo bioenergetic, proteomic, metabolomic, and tumor biology studies were performed to analyze changes in lung cancer cell metabolism and the consequences for hallmarks of cancer, including tumor growth, cancer cell invasion, and redox signaling. The findings revealed that nicotine-derived nitrosamine ketone (NNK) stimulates mitochondrial function and promotes lung tumor growth in vivo. These malignant properties were acquired from the induction of mitochondrial biogenesis induced by the upregulation and activation of the beta-2 adrenergic receptors (β2-AR)-cholinergic receptor nicotinic alpha 7 subunit (CHRNAα7)-dependent nitrosamine canonical signaling pathway. The observed NNK metabolic effects were mediated by TFAM overexpression and revealed a key role for mitochondrial reactive oxygen species and Annexin A1 in tumor growth promotion. Conversely, ectopic expression of the mitochondrial antioxidant enzyme manganese superoxide dismutase rescued the reprogramming and malignant metabolic effects of exposure to NNK and overexpression of TFAM, underlining the link between NNK and mitochondrial redox signaling in lung cancer. Innovation: Our findings describe the metabolic changes caused by NNK in a mechanistic framework for understanding how cigarette smoking causes lung cancer. Conclusion: Mitochondria play a role in the promotion of lung cancer induced by tobacco-specific nitrosamines. Antioxid. Redox Signal. 36, 525-549.

Papillary thyroid carcinoma tall cell variant shares accumulation of mitochondria, mitochondrial DNA mutations, and loss of oxidative phosphorylation complex I integrity with oncocytic tumors

Papillary thyroid carcinoma tall cell variant (PTC-TCV), a form of PTC regarded as an aggressive subtype, shares several morphologic features with oncocytic tumors. Pathogenic homoplasmic/highly heteroplasmic somatic mitochondrial DNA (mtDNA) mutations, usually affecting oxidative phosphorylation (OXPHOS) complex I subunits, are hallmarks of oncocytic cells. To clarify the relationship between PTC-TCV and oncocytic thyroid tumors, 17 PTC-TCV and 16 PTC non-TCV controls (cPTC) were subjected to: (1) transmission electron microscopy (TEM) to assess mitochondria accumulation, (2) next-generation sequencing to analyze mtDNA and nuclear genes frequently mutated in thyroid carcinoma, and (3) immunohistochemistry (IHC) for prohibitin and complex I subunit NDUFS4 to evaluate OXPHOS integrity. TEM showed replacement of cytoplasm by mitochondria in PTC-TCV but not in cPTC cells. All 17 PTC-TCV had at least one mtDNA mutation, totaling 21 mutations; 3/16 cPTC (19%) had mtDNA mutations (p < 0.001). PTC-TCV mtDNA mutations were homoplasmic/highly heteroplasmic, 16/21 (76%) mapping within mtDNA-encoded complex I subunits. MtDNA mutations in cPTC were homoplasmic in 2 cases and at low heteroplasmy in the third case, 2/3 mapping to mtDNA-encoded complex I subunits; 2/3 cPTC with mtDNA mutations had small tall cell subpopulations. PTC-TCV showed strong prohibitin expression and cPTC low-level expression, consistent with massive and limited mitochondrial content, respectively. All 17 PTC-TCV showed NDUFS4 loss (partial or complete) and 3 of 16 cPTC (19%) had (partial) NDUFS4 loss, consistent with lack of complex I integrity in PTC-TCV (p < 0.001). IHC loss of NDUFS4 was associated with mtDNA mutations (p < 0.001). Four BRAF V600E mutated PTCs had loss of NDUSF4 expression limited to neoplastic cell subpopulations with tall cell features, indicating that PTCs first acquire BRAF V600E and then mtDNA mutations. Similar to oncocytic thyroid tumors, PTC-TCV is characterized by mtDNA mutations, massive accumulation of mitochondria, and loss of OXPHOS integrity. IHC loss of NDUFS-4 can be used as a surrogate marker for OXPHOS disruption and to reliably diagnose PTC-TCV.

The frequency of the known mitochondrial variants associated with drug-induced toxicity in a Korean population

BACKGROUND

Few studies have annotated the whole mitochondrial DNA (mtDNA) genome associated with drug responses in Asian populations. This study aimed to characterize mtDNA genetic profiles, especially the distribution and frequency of well-known genetic biomarkers associated with diseases and drug-induced toxicity in a Korean population.

METHOD

Whole mitochondrial genome was sequenced for 118 Korean subjects by using a next-generation sequencing approach. The bioinformatic pipeline was constructed for variant calling, haplogroup classification and annotation of mitochondrial mutation.

RESULTS

A total of 681 variants was identified among all subjects. The MT-TRNP gene and displacement loop showed the highest numbers of variants (113 and 74 variants, respectively). The m.16189T > C allele, which is known to reduce the mtDNA copy number in human cells was detected in 25.4% of subjects. The variants (m.2706A > G, m.3010A > G, and m.1095T > C), which are associated with drug-induced toxicity, were observed with the frequency of 99.15%, 30.51%, and 0.08%, respectively. The m.2150T > A, a genotype associated with highly disruptive effects on mitochondrial ribosomes, was identified in five subjects. The D and M groups were the most dominant groups with the frequency of 34.74% and 16.1%, respectively.

CONCLUSIONS

Our finding was consistent with Korean Genome Project and well reflected the unique profile of mitochondrial haplogroup distribution. It was the first study to annotate the whole mitochondrial genome with drug-induced toxicity to predict the ADRs event in clinical implementation for Korean subjects. This approach could be extended for further study for validation of the potential ethnic-specific mitochondrial genetic biomarkers in the Korean population.

Extensive analysis of mitochondrial DNA quantity and sequence variation in human cumulus cells and assisted reproduction outcomes

STUDY QUESTION

Are relative mitochondrial DNA (mtDNA) content and mitochondrial genome (mtGenome) variants in human cumulus cells (CCs) associated with oocyte reproductive potential and assisted reproductive technology (ART) outcomes?

SUMMARY ANSWER

Neither the CC mtDNA quantity nor the presence of specific mtDNA genetic variants was associated with ART outcomes, although associations with patient body mass index (BMI) were detected, and the total number of oocytes retrieved differed between major mitochondrial haplogroups.

WHAT IS KNOWN ALREADY

CCs fulfil a vital role in the support of oocyte developmental competence. As with other cell types, appropriate cellular function is likely to rely upon adequate energy production, which in turn depends on the quantity and genetic competence of the mitochondria. mtDNA mutations can be inherited or they can accumulate in somatic cells over time, potentially contributing to aging. Such mutations may be homoplasmic (affecting all mtDNA in a cell) or they may display varying levels of heteroplasmy (affecting a proportion of the mtDNA). Currently, little is known concerning variation in CC mitochondrial genetics and how this might influence the reproductive potential of the associated oocyte.

STUDY DESIGN, SIZE, DURATION

This was a prospective observational study involving human CCs collected with 541 oocytes from 177 IVF patients. mtDNA quantity was measured in all the samples with a validated quantitative PCR method and the entire mtGenome was sequenced in a subset of 138 samples using a high-depth massively parallel sequencing approach. Associations between relative mtDNA quantity and mtGenome variants in CCs and patient age, BMI (kg/m2), infertility diagnosis and ART outcomes were investigated.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Massively parallel sequencing permitted not only the accurate detection of mutations but also the precise quantification of levels of mutations in cases of heteroplasmy. Sequence variants in the mtDNA were evaluated using Mitomaster and HmtVar to predict their potential impact.

MAIN RESULTS AND THE ROLE OF CHANCE

The relative mtDNA CC content was significantly associated with BMI. No significant associations were observed between CC mtDNA quantity and patient age, female infertility diagnosis or any ART outcome variable. mtGenome sequencing revealed 4181 genetic variants with respect to a reference genome. The COXI locus contained the least number of coding sequence variants, whereas ATPase8 had the most. The number of variants predicted to affect the ATP production differed significantly between mitochondrial macrohaplogroups. The total number of retrieved oocytes was different between the H-V and J-T as well as the U-K and J-T macrohaplogroups. There was a non-significant increase in mtDNA levels in CCs with heteroplasmic mitochondrial mutations.

LARGE SCALE DATA

N/A.

LIMITATIONS, REASONS FOR CAUTION

Although a large number of samples were analysed in this study, it was not possible to analyse all the CCs from every patient. Also, the results obtained with respect to specific clinical outcomes and macrohaplogroups should be interpreted with caution due to the smaller sample sizes when subdividing the dataset.

WIDER IMPLICATIONS OF THE FINDINGS

These findings suggest that the analysis of mtDNA in CCs is unlikely to provide an advantage in terms of improved embryo selection during assisted reproduction cycles. Nonetheless, our data raise interesting biological questions, particularly regarding the interplay of metabolism and BMI and the association of mtDNA haplogroup with oocyte yield in ovarian stimulation cycles.

STUDY FUNDING/COMPETING INTEREST(S)

This study was funded by National Institutes of Health grant 5R01HD092550-02. D.J.N. and C.R. co-hold patent US20150346100A1 and D.J.N. holds US20170039415A1, both for metabolic imaging methods. D.W. receives support from the NIHR Oxford Biomedical Research Centre. The remaining authors have no conflicts of interest to declare.

From Forensics to Clinical Research: Expanding the Variant Calling Pipeline for the Precision ID mtDNA Whole Genome Panel

Despite a multitude of methods for the sample preparation, sequencing, and data analysis of mitochondrial DNA (mtDNA), the demand for innovation remains, particularly in comparison with nuclear DNA (nDNA) research. The Applied Biosystems™ Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific, USA) is an innovative library preparation kit suitable for degraded samples and low DNA input. However, its bioinformatic processing occurs in the enterprise Ion Torrent Suite™ Software (TSS), yielding BAM files aligned to an unorthodox version of the revised Cambridge Reference Sequence (rCRS), with a heteroplasmy threshold level of 10%. Here, we present an alternative customizable pipeline, the PrecisionCallerPipeline (PCP), for processing samples with the correct rCRS output after Ion Torrent sequencing with the Precision ID library kit. Using 18 samples (3 original samples and 15 mixtures) derived from the 1000 Genomes Project, we achieved overall improved performance metrics in comparison with the proprietary TSS, with optimal performance at a 2.5% heteroplasmy threshold. We further validated our findings with 50 samples from an ongoing independent cohort of stroke patients, with PCP finding 98.31% of TSS's variants (TSS found 57.92% of PCP's variants), with a significant correlation between the variant levels of variants found with both pipelines.

iPSC culture expansion selects against putatively actionable mutations in the mitochondrial genome

Therapeutic application of induced pluripotent stem cell (iPSC) derivatives requires comprehensive assessment of the integrity of their nuclear and mitochondrial DNA (mtDNA) to exclude oncogenic potential and functional deficits. It is unknown, to which extent mtDNA variants originate from their parental cells or from de novo mutagenesis, and whether dynamics in heteroplasmy levels are caused by inter- and intracellular selection or genetic drift. Sequencing of mtDNA of 26 iPSC clones did not reveal evidence for de novo mutagenesis, or for any selection processes during reprogramming or differentiation. Culture expansion, however, selected against putatively actionable mtDNA mutations. Altogether, our findings point toward a scenario in which intracellular selection of mtDNA variants during culture expansion shapes the mutational landscape of the mitochondrial genome. Our results suggest that intercellular selection and genetic drift exert minor impact and that the bottleneck effect in context of the mtDNA genetic pool might have been overestimated.

•

Expansion culture selects against putatively actionable mtDNA mutations in iPSCs

•

Intracellular selection on mtDNA molecules shapes the mutational landscape

•

Random genetic drift and intercellular selection exert minor impact

•

Selection acts during culture expansion but not during reprogramming or differentiation

Hearing Impairment with Monoallelic GJB2 Variants: A GJB2 Cause or Non-GJB2 Cause?

Recessive variants in GJB2 are the most common genetic cause of sensorineural hearing impairment. However, in many patients, only one variant in the GJB2 coding region is identified using conventional sequencing strategy (eg, Sanger sequencing), resulting in nonconfirmative diagnosis. Conceivably, there might be other unidentified pathogenic variants in the noncoding region of GJB2 or other deafness-causing genes in these patients. To address this, a next-generation sequencing-based diagnostic panel targeting the entire GJB2 gene and the coding regions of 158 other known deafness-causing genes was designed and applied to 95 patients with nonsyndromic sensorineural hearing impairment (including 81 Han Taiwanese and 14 Mongolian patients) in whom only a single GJB2 variant had been detected using conventional Sanger sequencing. The panel confirmed the genetic diagnosis in 24 patients (25.3%). Twenty-two of them had causative variants in several deafness-causing genes other than GJB2, including MYO15A, MYO7A, TECTA, POU4F3, KCNQ4, SLC26A4, OTOF, MT-RNR1, MITF, WFS1, and USH2A. The other two patients had causative variants in GJB2, including a Taiwanese patient with a mosaic maternal uniparental disomy c.235delC variant (approximately 69% mosaicism) and a Mongolian patient with compound heterozygous c.35dupG and c.35delG variants, which occurred at the same site. This study demonstrates the utility of next-generation sequencing in clarifying the genetic diagnosis of hearing-impaired patients with nonconfirmative GJB2 genotypes on conventional genetic examinations.

An atlas of mitochondrial DNA genotype-phenotype associations in the UK Biobank

Mitochondrial DNA (mtDNA) variation in common diseases has been underexplored, partly due to a lack of genotype calling and quality-control procedures. Developing an at-scale workflow for mtDNA variant analyses, we show correlations between nuclear and mitochondrial genomic structures within subpopulations of Great Britain and establish a UK Biobank reference atlas of mtDNA-phenotype associations. A total of 260 mtDNA-phenotype associations were new (P < 1 × 10-5), including rs2853822 /m.8655 C>T (MT-ATP6) with type 2 diabetes, rs878966690 /m.13117 A>G (MT-ND5) with multiple sclerosis, 6 mtDNA associations with adult height, 24 mtDNA associations with 2 liver biomarkers and 16 mtDNA associations with parameters of renal function. Rare-variant gene-based tests implicated complex I genes modulating mean corpuscular volume and mean corpuscular hemoglobin. Seven traits had both rare and common mtDNA associations, where rare variants tended to have larger effects than common variants. Our work illustrates the value of studying mtDNA variants in common complex diseases and lays foundations for future large-scale mtDNA association studies.

Clinical Bioinformatics in Precise Diagnosis of Mitochondrial Disease

Clinical bioinformatics system is well-established for diagnosing genetic disease based on next-generation sequencing, but requires special considerations when being adapted for the next-generation sequencing-based genetic diagnosis of mitochondrial diseases. Challenges are caused by the involvement of mitochondrial DNA genome in disease etiology. Heteroplasmy and haplogroup are key factors in interpreting mitochondrial DNA variant effects. Data resources and tools for analyzing variant and sequencing data are available at MSeqDR, MitoMap, and HmtDB. Revised specifications of the American College of Medical Genetics/Association of Molecular Pathology standards and guidelines for mitochondrial DNA variant interpretation are proposed by the MSeqDr Consortium and community experts.

Mitochondrial genome variations are associated with amyotrophic lateral sclerosis in patients from mainland China

BACKGROUND

Amyotrophic lateral sclerosis (ALS) is a late-onset neurodegenerative disorder. Mitochondrial dysfunction is involved in the complex pathophysiology of ALS; however, the role of mitochondrial DNA (mtDNA) variants in ALS is poorly understood. We aimed to elucidate the role of mtDNA variants in the pathogenesis of ALS.

METHODS

The mitochondrial haplogroups of 585 ALS patients and 371 healthy controls were determined; 38 ALS patients and 42 controls underwent long-range polymerase chain reaction combined with next-generation sequencing technology to analyze whole mitochondrial genome variants.

RESULTS

A higher percentage of variants accumulated in ALS patients than in controls. Analysis of coding region variations that were further stratified by mtDNA genes revealed that nonsynonymous variants were more vulnerable in ALS patients than in controls, particularly in the ND4L, ND5, and ATP8 genes. Moreover, pathogenic nonsynonymous variants tended to over-represent in ALS patients. Unsurprisingly, nonsynonymous variants were not related to the phenotype. Haplogroup analysis did not found evidence of association between haplogroups with the risk of ALS, however, patients belonging to haplogroup Y and M7c were prone to develop later onset of ALS.

CONCLUSIONS

This is the first study to profile mtDNA variants in ALS patients from mainland China. Our results suggest that an increase in the number of nonsynonymous variants is linked to the pathogenesis of ALS. Moreover, haplogroup Y and M7c may modulate the clinical expression of ALS. Our findings provide independent, albeit limited, evidence for the role of mtDNA in the pathogenesis of ALS.

Application of Genome Sequencing from Blood to Diagnose Mitochondrial Diseases

Mitochondrial diseases can be caused by pathogenic variants in nuclear or mitochondrial DNA-encoded genes that often lead to multisystemic symptoms and can have any mode of inheritance. Using a single test, Genome Sequencing (GS) can effectively identify variants in both genomes, but it has not yet been universally used as a first-line approach to diagnosing mitochondrial diseases due to related costs and challenges in data analysis. In this article, we report three patients with mitochondrial disease molecularly diagnosed through GS performed on DNA extracted from blood to demonstrate different diagnostic advantages of this technology, including the detection of a low-level heteroplasmic pathogenic variant, an intragenic nuclear DNA deletion, and a large mtDNA deletion. Current technical improvements and cost reductions are likely to lead to an expanded routine diagnostic usage of GS and of the complementary “Omic” technologies in mitochondrial diseases.

Low frequency mitochondrial DNA heteroplasmy SNPs in blood, retina, and [RPE+choroid] of age-related macular degeneration subjects

Purpose

Mitochondrial (mt) DNA damage is associated with age-related macular degeneration (AMD) and other human aging diseases. This study was designed to quantify and characterize mtDNA low-frequency heteroplasmy single nucleotide polymorphisms (SNPs) of three different tissues isolated from AMD subjects using Next Generation Sequencing (NGS) technology.

Methods

DNA was extracted from neural retina, [RPE+choroid] and blood from three deceased age-related macular degeneration (AMD) subjects. Entire mitochondrial genomes were analyzed for low-frequency heteroplasmy SNPs using NGS technology that independently sequenced both mtDNA strands. This deep sequencing method (average sequencing depth of 30,000; range 1,000–100,000) can accurately differentiate low-frequency heteroplasmy SNPs from DNA modification artifacts. Twenty-three ‘hot-spot’ heteroplasmy mtDNA SNPs were analyzed in 222 additional blood samples.

Results

Germline homoplasmy SNPs that defined mtDNA haplogroups were consistent in the three tissues of each subject. Analyses of SNPs with <40% heteroplasmy revealed the blood had significantly greater numbers of heteroplasmy SNPs than retina alone (p≤0.05) or retina+choroid combined (p = 0.008). Twenty-three ‘hot-spot’ mtDNA heteroplasmy SNPs were present, with three being non-synonymous (amino acid change). Four ‘hot-spot’ heteroplasmy SNPs (m.1120C>T, m.1284T>C, m.1556C>T, m.7256C>T) were found in additional samples (n = 222). Five heteroplasmy SNPs (m.4104A>G, m.5320C>T, m.5471G>A, m.5474A>G, m.5498A>G) declined with age. Two heteroplasmy SNPs (m.13095T>C, m.13105A>G) increased in AMD compared to Normal samples. In the heteroplasmy SNPs, very few transversion mutations (purine to pyrimidine or vice versa, associated with oxidative damage) were found and the majority were transition changes (purine to purine or pyrimidine to pyrimidine, associated with replication errors).

Conclusion

Within an individual, the blood, retina and [RPE+choroid] contained identical homoplasmy SNPs representing inherited germline mtDNA haplogroup. NGS methodology showed significantly more mtDNA heteroplasmy SNPs in blood compared to retina and [RPE+choroid], suggesting the latter tissues have substantial protection. Significantly higher heteroplasmy levels of m.13095T>C and m.13105A>G may represent potential AMD biomarkers. Finally, high levels of transition mutations suggest that accumulation of heteroplasmic SNPs may occur through replication errors rather than oxidative damage.

Analysis of human mitochondrial genome co-occurrence networks of Asian population at varying altitudes

Networks have been established as an extremely powerful framework to understand and predict the behavior of many large-scale complex systems. We studied network motifs, the basic structural elements of networks, to describe the possible role of co-occurrence of genomic variations behind high altitude adaptation in the Asian human population. Mitochondrial DNA (mtDNA) variations have been acclaimed as one of the key players in understanding the biological mechanisms behind adaptation to extreme conditions. To explore the cumulative effects of variations in the mitochondrial genome with the variation in the altitude, we investigated human mt-DNA sequences from the NCBI database at different altitudes under the co-occurrence motifs framework. Analysis of the co-occurrence motifs using similarity clustering revealed a clear distinction between lower and higher altitude regions. In addition, the previously known high altitude markers 3394 and 7697 (which are definitive sites of haplogroup M9a1a1c1b) were found to co-occur within their own gene complexes indicating the impact of intra-genic constraint on co-evolution of nucleotides. Furthermore, an ancestral 'RSRS50' variant 10,398 was found to co-occur only at higher altitudes supporting the fact that a separate route of colonization at these altitudes might have taken place. Overall, our analysis revealed the presence of co-occurrence interactions specific to high altitude at a whole mitochondrial genome level. This study, combined with the classical haplogroups analysis is useful in understanding the role of co-occurrence of mitochondrial variations in high altitude adaptation.

Mitochondrial DNA analysis efficiently contributes to the identification of metastatic contralateral breast cancers

Purpose

In daily practice, a contralateral breast cancer (CBC) is usually considered as a new independent tumor despite the indications of several studies showing that the second neoplasia may be a metastatic spread of the primary tumor. Recognition of clonal masses in the context of multiple synchronous or metachronous tumors is crucial for correct prognosis, therapeutic choice, and patient management. Mitochondrial DNA (mtDNA) sequencing shows high informative potential in the diagnosis of synchronous neoplasms, based on the fact that somatic mtDNA mutations are non-recurrent events, whereas tumors sharing them have a common origin. We here applied this technique to reveal clonality of the CBC with respect to the first tumor.

Methods

We analyzed 30 sample pairs of primary breast cancers and synchronous or metachronous CBCs with detailed clinical information available and compared standard clinico-pathological criteria with mtDNA sequencing to reveal the metastatic nature of CBCs.

Results

MtDNA analysis was informative in 23% of the cases, for which it confirmed a clonal origin of the second tumor. In addition, it allowed to solve two ambiguous cases where histopathological criteria had failed to be conclusive and to suggest a clonal origin for two additional cases that had been classified as independent by pathologists.

Conclusion

Overall, the mtDNA-based classification showed a more accurate predictive power than standard histopathology in identifying cases of metastatic rather than bilateral breast cancers in our cohort, suggesting that mtDNA sequencing may be a more precise and easy-to-use method to be introduced in daily routine to support and improve histopathological diagnoses.

Possible A2E Mutagenic Effects on RPE Mitochondrial DNA from Innovative RNA-Seq Bioinformatics Pipeline

Mitochondria are subject to continuous oxidative stress stimuli that, over time, can impair their genome and lead to several pathologies, like retinal degenerations. Our main purpose was the identification of mtDNA variants that might be induced by intense oxidative stress determined by N-retinylidene-N-retinylethanolamine (A2E), together with molecular pathways involving the genes carrying them, possibly linked to retinal degeneration. We performed a variant analysis comparison between transcriptome profiles of human retinal pigment epithelial (RPE) cells exposed to A2E and untreated ones, hypothesizing that it might act as a mutagenic compound towards mtDNA. To optimize analysis, we proposed an integrated approach that foresaw the complementary use of the most recent algorithms applied to mtDNA data, characterized by a mixed output coming from several tools and databases. An increased number of variants emerged following treatment. Variants mainly occurred within mtDNA coding sequences, corresponding with either the polypeptide-encoding genes or the RNA. Time-dependent impairments foresaw the involvement of all oxidative phosphorylation complexes, suggesting a serious damage to adenosine triphosphate (ATP) biosynthesis, that can result in cell death. The obtained results could be incorporated into clinical diagnostic settings, as they are hypothesized to modulate the phenotypic expression of mtDNA pathogenic variants, drastically improving the field of precision molecular medicine.

Specifications of the ACMG/AMP standards and guidelines for mitochondrial DNA variant interpretation

Mitochondrial DNA (mtDNA) variant pathogenicity interpretation has special considerations given unique features of the mtDNA genome, including maternal inheritance, variant heteroplasmy, threshold effect, absence of splicing, and contextual effects of haplogroups. Currently, there are insufficient standardized criteria for mtDNA variant assessment, which leads to inconsistencies in clinical variant pathogenicity reporting. An international working group of mtDNA experts was assembled within the Mitochondrial Disease Sequence Data Resource Consortium and obtained Expert Panel status from ClinGen. This group reviewed the 2015 American College of Medical Genetics and Association of Molecular Pathology standards and guidelines that are widely used for clinical interpretation of DNA sequence variants and provided further specifications for additional and specific guidance related to mtDNA variant classification. These Expert Panel consensus specifications allow for consistent consideration of the unique aspects of the mtDNA genome that directly influence variant assessment, including addressing mtDNA genome composition and structure, haplogroups and phylogeny, maternal inheritance, heteroplasmy, and functional analyses unique to mtDNA, as well as specifications for utilization of mtDNA genomic databases and computational algorithms.

Ancient mitochondrial DNA pathogenic variants putatively associated with mitochondrial disease

Mitochondrial DNA variants associated with diseases are widely studied in contemporary populations, but their prevalence has not yet been investigated in ancient populations. The publicly available AmtDB database contains 1443 ancient mtDNA Eurasian genomes from different periods. The objective of this study was to use this data to establish the presence of pathogenic mtDNA variants putatively associated with mitochondrial diseases in ancient populations. The clinical significance, pathogenicity prediction and contemporary frequency of mtDNA variants were determined using online platforms. The analyzed ancient mtDNAs contain six variants designated as being "confirmed pathogenic" in modern patients. The oldest of these, m.7510T>C in the MT-TS1 gene, was found in a sample from the Neolithic period, dated 5800-5400 BCE. All six have well established clinical association, and their pathogenic effect is corroborated by very low population frequencies in contemporary populations. Analysis of the geographic location of the ancient samples, contemporary epidemiological trends and probable haplogroup association indicate diverse spatiotemporal dynamics of these variants. The dynamics in the prevalence and distribution is conceivably result of de novo mutations or human migrations and subsequent evolutionary processes. In addition, ten variants designated as possibly or likely pathogenic were found, but the clinical effect of these is not yet well established and further research is warranted. All detected mutations putatively associated with mitochondrial disease in ancient mtDNA samples are in tRNA coding genes. Most of these mutations are in a mt-tRNA type (Model 2) that is characterized by loss of D-loop/T-loop interaction. Exposing pathogenic variants in ancient human populations expands our understanding of their origin and prevalence dynamics.

A Novel Approach for the Identification of Pharmacogenetic Variants in MT-RNR1 through Next-Generation Sequencing Off-Target Data

Specific genetic variants in the mitochondrially encoded 12S ribosomal RNA gene (MT-RNR1) cause aminoglycoside-induced irreversible hearing loss. Mitochondrial DNA is usually not included in targeted sequencing experiments; however, off-target data may deliver this information. Here, we extract MT-RNR1 genetic variation, including the most relevant ototoxicity variant m.1555A>G, using the off-target reads of 473 research samples, sequenced through a capture-based, custom-targeted panel and whole exome sequencing (WES), and of 1245 diagnostic samples with clinical WES. Sanger sequencing and fluorescence-based genotyping were used for genotype validation. There was a correlation between off-target reads and mitochondrial coverage (r_customPanel = 0.39, p = 2 × 10⁻¹³ and rWES = 0.67, p = 7 × 10⁻²¹). The median read depth of MT-RNR1 m.1555 was similar to the average mitochondrial genome coverage, with saliva and blood samples giving comparable results. The genotypes from 415 samples, including three m.1555G carriers, were concordant with fluorescence-based genotyping data. In clinical WES, median MT-RNR1 coverage was 56×, with 90% of samples having ≥20 reads at m.1555 position, and one m.1494T and three m.1555G carriers were identified with no evidence for heteroplasmy. Altogether, this study shows that obtaining MT-RNR1 genotypes through off-target reads is an efficient strategy that can impulse preemptive pharmacogenetic screening of this mitochondrial gene.

Investigating Human Mitochondrial Genomes in Single Cells

Mitochondria host multiple copies of their own small circular genome that has been extensively studied to trace the evolution of the modern eukaryotic cell and discover important mutations linked to inherited diseases. Whole genome and exome sequencing have enabled the study of mtDNA in a large number of samples and experimental conditions at single nucleotide resolution, allowing the deciphering of the relationship between inherited mutations and phenotypes and the identification of acquired mtDNA mutations in classical mitochondrial diseases as well as in chronic disorders, ageing and cancer. By applying an ad hoc computational pipeline based on our MToolBox software, we reconstructed mtDNA genomes in single cells using whole genome and exome sequencing data obtained by different amplification methodologies (eWGA, DOP-PCR, MALBAC, MDA) as well as data from single cell Assay for Transposase Accessible Chromatin with high-throughput sequencing (scATAC-seq) in which mtDNA sequences are expected as a byproduct of the technology. We show that assembled mtDNAs, with the exception of those reconstructed by MALBAC and DOP-PCR methods, are quite uniform and suitable for genomic investigations, enabling the study of various biological processes related to cellular heterogeneity such as tumor evolution, neural somatic mosaicism and embryonic development.

A data and text mining pipeline to annotate human mitochondrial variants with functional and clinical information

BACKGROUND

Human mitochondrial DNA has an important role in the cellular energy production through oxidative phosphorylation. Therefore, this process may be the cause and have an effect on mitochondrial DNA mutability, functional alteration, and disease onset related to a wide range of different clinical expressions and phenotypes. Although a large part of the observed variations is fixed in a population and hence expected to be benign, the estimation of the degree of the pathogenicity of any possible human mitochondrial DNA variant is clinically pivotal.

METHODS

In this scenario, the establishment of standard criteria based on functional studies is required. In this context, a "data and text mining" pipeline is proposed here, developed using the programming language R, capable of extracting information regarding mitochondrial DNA functional studies and related clinical assessments from the literature, thus improving the annotation of human mitochondrial variants reported in the HmtVar database.

RESULTS

The data mining pipeline has produced a list of 1,073 Pubmed IDs (PMIDs) from which the text mining pipeline has retrieved information on 932 human mitochondrial variants regarding experimental validation and clinical features.

CONCLUSIONS

The application of the pipeline will contribute to supporting the interpretation of pathogenicity of human mitochondrial variants by facilitating diagnosis to clinicians and researchers faced with this task.

European mtDNA Variants Are Associated With Differential Responses to Cisplatin, an Anticancer Drug: Implications for Drug Resistance and Side Effects

Background: Cisplatin, a powerful antitumor agent, causes formation of DNA adducts, and activation of apoptotic pathways. Presently, cisplatin resistance develops in up to 70% of patients but the underlying molecular mechanism(s) are unclear and there are no markers to determine which patients will become resistant. Mitochondria play a significant role not only in energy metabolism but also retrograde signaling (mitochondria to nucleus) that modulates inflammation, complement, and apoptosis pathways. Maternally inherited mitochondrial (mt) DNA can be classified into haplogroups representing different ethnic populations that have diverse susceptibilities to diseases and medications. Methods: Transmitochondrial cybrids, where all cell lines possess identical nuclear genomes but either the H (Southern European) or J (Northern European) mtDNA haplogroups, were treated with cisplatin and analyzed for differential responses related to viability, oxidative stress, and expression levels of genes associated with cancer, cisplatin-induced nephrotoxicity and resistance, apoptosis and signaling pathways. Results: The cisplatin-treated-J cybrids showed greater loss of cell viability along with lower levels of reactive oxygen species and mitochondrial membrane potential compared to cisplatin-treated-H cybrids. After cisplatin treatment, J cybrids showed increased gene expression of BAX, CASP3, and CYP51A, but lower levels of SFRP1 compared to untreated-J cybrids. The cisplatin-treated-H cybrids had elevated expression of CDKN1A/P21, which has a role in cisplatin toxicity, compared to untreated-H cybrids. The cisplatin-treated H had higher transcription levels of ABCC1, DHRS2/HEP27, and EFEMP1 compared to cisplatin-treated-J cybrids. Conclusions: Cybrid cell lines that contain identical nuclei but either H mtDNA mitochondria or J mtDNA mitochondria respond differently to cisplatin treatments suggesting involvement of the retrograde signaling (from mitochondria to nucleus) in the drug-induced cell death. Varying toxicities and transcription levels of the H vs. J cybrids after cisplatin treatment support the hypothesis that mtDNA variants play a role in the expression of genes affecting resistance and side effects of cisplatin.

Synergistic Effect of Mitochondrial and Lysosomal Dysfunction in Parkinson's Disease

Crosstalk between lysosomes and mitochondria plays a central role in Parkinson’s Disease (PD). Lysosomal function may be influenced by mitochondrial quality control, dynamics and/or respiration, but whether dysfunction of endocytic or autophagic pathway is associated with mitochondrial impairment determining accumulation of defective mitochondria, is not yet understood. Here, we performed live imaging, western blotting analysis, sequencing of mitochondrial DNA (mtDNA) and senescence-associated beta-galactosidase activity assay on primary fibroblasts from a young patient affected by PD, her mother and a healthy control to analyze the occurrence of mtDNA mutations, lysosomal abundance, acidification and function, mitochondrial biogenesis activation and senescence. We showed synergistic alterations in lysosomal functions and mitochondrial biogenesis, likely associated with a mitochondrial genetic defect, with a consequent block of mitochondrial turnover and occurrence of premature cellular senescence in PARK2-PD fibroblasts, suggesting that these alterations represent potential mechanisms contributing to the loss of dopaminergic neurons.

Mitochondrial-Derived Vesicles as Candidate Biomarkers in Parkinson's Disease: Rationale, Design and Methods of the EXosomes in PArkiNson Disease (EXPAND) Study

The progressive loss of dopaminergic neurons in the nigro-striatal system is a major trait of Parkinson’s disease (PD), manifesting clinically as motor and non-motor symptoms. Mitochondrial dysfunction and oxidative stress are alleged pathogenic mechanisms underlying aggregation of misfolded α-synuclein that in turn triggers dopaminergic neurotoxicity. Peripheral processes, including inflammation, may precede and contribute to neurodegeneration. Whether mitochondrial dyshomeostasis in the central nervous system and systemic inflammation are linked to one another in PD is presently unclear. Extracellular vesicles (EVs) are delivery systems through which cells can communicate or unload noxious materials. EV trafficking also participates in mitochondrial quality control (MQC) by generating mitochondrial-derived vesicles to dispose damaged organelles. Disruption of MQC coupled with abnormal EV secretion may play a role in the pathogenesis of PD. Furthermore, due to its bacterial ancestry, circulating mitochondrial DNA can elicit an inflammatory response. Therefore, purification and characterisation of molecules packaged in, and secreted through, small EVs (sEVs)/exosomes in body fluids may provide meaningful insights into the association between mitochondrial dysfunction and systemic inflammation in PD. The EXosomes in PArkiNson Disease (EXPAND) study was designed to characterise the cargo of sEVs/exosomes isolated from the serum of PD patients and to identify candidate biomarkers for PD.

Uncovering Missing Heritability in Rare Diseases

The problem of 'missing heritability' affects both common and rare diseases hindering: discovery, diagnosis, and patient care. The 'missing heritability' concept has been mainly associated with common and complex diseases where promising modern technological advances, like genome-wide association studies (GWAS), were unable to uncover the complete genetic mechanism of the disease/trait. Although rare diseases (RDs) have low prevalence individually, collectively they are common. Furthermore, multi-level genetic and phenotypic complexity when combined with the individual rarity of these conditions poses an important challenge in the quest to identify causative genetic changes in RD patients. In recent years, high throughput sequencing has accelerated discovery and diagnosis in RDs. However, despite the several-fold increase (from ~10% using traditional to ~40% using genome-wide genetic testing) in finding genetic causes of these diseases in RD patients, as is the case in common diseases-the majority of RDs are also facing the 'missing heritability' problem. This review outlines the key role of high throughput sequencing in uncovering genetics behind RDs, with a particular focus on genome sequencing. We review current advances and challenges of sequencing technologies, bioinformatics approaches, and resources.

Bioinformatics Tools and Databases to Assess the Pathogenicity of Mitochondrial DNA Variants in the Field of Next Generation Sequencing

The development of next generation sequencing (NGS) has greatly enhanced the diagnosis of mitochondrial disorders, with a systematic analysis of the whole mitochondrial DNA (mtDNA) sequence and better detection sensitivity. However, the exponential growth of sequencing data renders complex the interpretation of the identified variants, thereby posing new challenges for the molecular diagnosis of mitochondrial diseases. Indeed, mtDNA sequencing by NGS requires specific bioinformatics tools and the adaptation of those developed for nuclear DNA, for the detection and quantification of mtDNA variants from sequence alignment to the calling steps, in order to manage the specific features of the mitochondrial genome including heteroplasmy, i.e., coexistence of mutant and wildtype mtDNA copies. The prioritization of mtDNA variants remains difficult, relying on a limited number of specific resources: population and clinical databases, and in silico tools providing a prediction of the variant pathogenicity. An evaluation of the most prominent bioinformatics tools showed that their ability to predict the pathogenicity was highly variable indicating that special efforts should be directed at developing new bioinformatics tools dedicated to the mitochondrial genome. In addition, massive parallel sequencing raised several issues related to the interpretation of very low mtDNA mutational loads, discovery of variants of unknown significance, and mutations unrelated to patient phenotype or the co-occurrence of mtDNA variants. This review provides an overview of the current strategies and bioinformatics tools for accurate annotation, prioritization and reporting of mtDNA variations from NGS data, in order to carry out accurate genetic counseling in individuals with primary mitochondrial diseases.