MtSNPscore: a combined evidence approach for assessing cumulative impact of mitochondrial variations in disease

MitoLink: A generic integrated web-based workflow system to evaluate genotype-phenotype correlations in human mitochondrial diseases: Observations from the GenomeAsia Pilot project

MitoLink is a generic, scalable and modular web-based workflow system developed to study genotype-phenotype correlations in human mitochondrial diseases. MitoLink integrates applications for assessment of genomic variation and currently houses genome-wide datasets from GenomeAsia Pilot project, gnomAD, ClinVar and DisGenNet. In this study, a reference list of nearly 3975 proteins (both nuclear and mitochondrial encoded) with mitochondrial function is reported. This protein set is mapped to disease associated variants in the GenomeAsia Pilot Project and DisGenNet and evaluated for pathogenicity as defined by ClinVar. Observations of shared genetic components in potential comorbidities are discussed from gene-disease network in Asian population, however, the platform is generic and can be applied to any population dataset. MitoLink is a unique customized workflow system that allows for systematic storage, extraction, analysis and visualization of genomic variation to understand genotype-phenotype correlations for mitochondrial diseases. Given the modularity of tool and data integration, MitoLink is a scalable system that can accommodate a diverse set of applications linked via standard data structure within the framework of Galaxy. MitoLink is built on FAIR principles and supports creation of reproducible workflows towards understanding genotype-phenotype correlations across several disease phenotypes globally.

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.

MtBrowse: An integrative genomics browser for human mitochondrial DNA

The human mitochondrion is a unique semi-autonomous organelle with a genome of its own and also requires nuclear encoded components to carry out its functions. In addition to being the powerhouse of the cell, mitochondria plays a central role in several metabolic pathways. It is therefore challenging to delineate the cause-effect relationship in context of mitochondrial dysfunction. Several studies implicate mutations in mitochondrial DNA (mtDNA) in various complex diseases. The human mitochondrial DNA (mtDNA) encodes a set of 37 genes, 13 protein coding, 22 tRNAs and two ribosomal RNAs, which are essential structural and functional components of the electron transport chain. As mentioned above, variations in these genes have been implicated in a broad spectrum of diseases and are extensively reported in literature and various databases. A large number of databases and prediction methods have been published to elucidate the role of human mitochondrial DNA in various disease phenotypes. However, there is no centralized resource to visualize this genotype-phenotype data. Towards this, we have developed MtBrowse: an integrative genomics browser for human mtDNA. As of now, MtBrowse has four categories - Gene, Disease, Reported variation and Variation prediction. These categories have 105 tracks and house data on mitochondrial reference genes, around 600 variants reported in literature with respect to various disease phenotypes and predictions for potential pathogenic variations in protein-coding genes. MtBrowse also hosts genomic variation data from over 5000 individuals on 22 disease phenotypes. MtBrowse may be accessed at http://ab-openlab.csir.res.in/cgi-bin/gb2/gbrowse.

PON-mt-tRNA: a multifactorial probability-based method for classification of mitochondrial tRNA variations

Transfer RNAs (tRNAs) are essential for encoding the transcribed genetic information from DNA into proteins. Variations in the human tRNAs are involved in diverse clinical phenotypes. Interestingly, all pathogenic variations in tRNAs are located in mitochondrial tRNAs (mt-tRNAs). Therefore, it is crucial to identify pathogenic variations in mt-tRNAs for disease diagnosis and proper treatment. We collected mt-tRNA variations using a classification based on evidence from several sources and used the data to develop a multifactorial probability-based prediction method, PON-mt-tRNA, for classification of mt-tRNA single nucleotide substitutions. We integrated a machine learning-based predictor and an evidence-based likelihood ratio for pathogenicity using evidence of segregation, biochemistry and histochemistry to predict the posterior probability of pathogenicity of variants. The accuracy and Matthews correlation coefficient (MCC) of PON-mt-tRNA are 1.00 and 0.99, respectively. In the absence of evidence from segregation, biochemistry and histochemistry, PON-mt-tRNA classifies variations based on the machine learning method with an accuracy and MCC of 0.69 and 0.39, respectively. We classified all possible single nucleotide substitutions in all human mt-tRNAs using PON-mt-tRNA. The variations in the loops are more often tolerated compared to the variations in stems. The anticodon loop contains comparatively more predicted pathogenic variations than the other loops. PON-mt-tRNA is available at http://structure.bmc.lu.se/PON-mt-tRNA/.

A multi-parametric workflow for the prioritization of mitochondrial DNA variants of clinical interest

Assigning a pathogenic role to mitochondrial DNA (mtDNA) variants and unveiling the potential involvement of the mitochondrial genome in diseases are challenging tasks in human medicine. Assuming that rare variants are more likely to be damaging, we designed a phylogeny-based prioritization workflow to obtain a reliable pool of candidate variants for further investigations. The prioritization workflow relies on an exhaustive functional annotation through the mtDNA extraction pipeline MToolBox and includes Macro Haplogroup Consensus Sequences to filter out fixed evolutionary variants and report rare or private variants, the nucleotide variability as reported in HmtDB and the disease score based on several predictors of pathogenicity for non-synonymous variants. Cutoffs for both the disease score as well as for the nucleotide variability index were established with the aim to discriminate sequence variants contributing to defective phenotypes. The workflow was validated on mitochondrial sequences from Leber’s Hereditary Optic Neuropathy affected individuals, successfully identifying 23 variants including the majority of the known causative ones. The application of the prioritization workflow to cancer datasets allowed to trim down the number of candidate for subsequent functional analyses, unveiling among these a high percentage of somatic variants. Prioritization criteria were implemented in both standalone (http://sourceforge.net/projects/mtoolbox/) and web version (https://mseqdr.org/mtoolbox.php) of MToolBox.

Investigation of mitochondrial DNA variations among Indian Friedreich's ataxia (FRDA) patients

OBJECTIVE

The loss of function mutations (biallelic) in frataxin (FXN) has primarily been implicated in Friedreich's ataxia (FRDA), an autosomal recessive cerebellar ataxia. The protein product of FXN is a nuclear-encoded mitochondrial protein required for the biogenesis of iron- clusters (Fe-S). FRDA is characterized by neurological and non-neurological features which show variable expression in affected individuals. An inverse relationship has been demonstrated between GAA repeat size and age at onset and explains 50% variability of the age at onset. MtDNA variations and haplogroups could be one of the contributory factors to explain the remaining heterogeneity in FRDA, since mitochondrial oxidative stress is thought to be involved in the pathogenesis of FRDA.

METHODS

In our study, targeted resequencing of the D-loop and coding region of mitochondrial genes (ND1-6 and ATP) was conducted in 30 genetically confirmed FRDA patients and 62 ethnicity-matched unrelated healthy controls to identify the functionally important mtDNA variations and to trace the mitochondrial lineage of Indian FRDA patients. Cumulative mitochondrial SNP scores were computed for the identified variations in the functional region and haplogroups were determined by Haplogrep.

RESULTS

A significantly higher load of overall mitochondrial variations (with a trend toward the coding region) per individual was noted among FRDA cases rather than controls (p-value<0.03). A non-synonymous variation (p. L237M) in ND2 was over-represented among FRDA cases (p-value 0.04). This variation has a reported association with longevity and myocardial infarction. We also observed over-representation of H haplogroup (Caucasian mitochondrial haplogroup) among FRDA patients. We have not observed the influence of mitochondrial variations and haplogroup upon age at onset of FRDA.

CONCLUSIONS

Overall, our study identifies the functionally important variations and mitochondrial lineage of Indian FRDA cases and, that underscores the importance of studying the role of mitochondrial genome variations in FRDA.

Resources, challenges and way forward in rare mitochondrial diseases research

Over 300 million people are affected by about 7000 rare diseases globally. There are tremendous resource limitations and challenges in driving research and drug development for rare diseases. Hence, innovative approaches are needed to identify potential solutions. This review focuses on the resources developed over the past years for analysis of genome data towards understanding disease biology especially in the context of mitochondrial diseases, given that mitochondria are central to major cellular pathways and their dysfunction leads to a broad spectrum of diseases. Platforms for collaboration of research groups, clinicians and patients and the advantages of community collaborative efforts in addressing rare diseases are also discussed. The review also describes crowdsourcing and crowdfunding efforts in rare diseases research and how the upcoming initiatives for understanding disease biology including analyses of large number of genomes are also applicable to rare diseases.

Resources, challenges and way forward in rare mitochondrial diseases research

Over 300 million people are affected by about 7000 rare diseases globally. There are tremendous resource limitations and challenges in driving research and drug development for rare diseases. Hence, innovative approaches are needed to identify potential solutions. This review focuses on the resources developed over the past years for analysis of genome data towards understanding disease biology especially in the context of mitochondrial diseases, given that mitochondria are central to major cellular pathways and their dysfunction leads to a broad spectrum of diseases. Platforms for collaboration of research groups, clinicians and patients and the advantages of community collaborative efforts in addressing rare diseases are also discussed. The review also describes crowdsourcing and crowdfunding efforts in rare diseases research and how the upcoming initiatives for understanding disease biology including analyses of large number of genomes are also applicable to rare diseases.

Mitochondrial ATP 6 and 8 polymorphisms in irritable bowel syndrome with diarrhea

AIM: To investigate mitochondrial ATP 6 and 8 polymorphisms in the colon and ileum of patients with irritable bowel syndrome with diarrhea (IBS-D).

METHODS: Twenty-eight patients fulfilling the Rome III criteria for IBS-D and 28 healthy subjects were investigated. All study participants underwent screening colonoscopy and mucosal biopsies were obtained from the colon and/or terminal ileum. Genomic DNA was extracted from specimens based on standard protocols. Mitochondrial ATP (MT-ATP) 6 and 8 genes in specimens were polymerase chain reaction amplified and sequenced. Sequencing data were analyzed via Variant Reporter™ Software and compared with the reference sequence from Genbank (accession No. NC_012920) to indicate possible polymorphisms. The protocol was registered at www.clinicaltrials.gov as NCT01028898.

RESULTS: Twenty-five polymorphic sites of MT-ATP 6 and 8 genes were detected and 12 of them were missense mutations. A median of two polymorphic sites in MT-ATP genes was found in colon specimens of controls while a median of three polymorphic sites was noted in patients with IBS-D (Mann-Whitney test, P = 0.012). The variants of the colon and ileum specimens from the same subjects were identical in all but one case. Symptom duration in IBS was not found to be a significant factor associated with the mtDNA polymorphism (Spearman correlation, P = 0.592). The mitochondrial DNA change at 8860 was present in all cases of both groups. The frequency of the 8701 polymorphism was found to be the second most frequent; however, no statistical difference was noted between the groups (χ² test, P = 0.584).

CONCLUSION: Patients with IBS-D have a higher incidence of MT-ATP 6 and 8 polymorphisms than healthy subjects, implying that the mtDNA polymorphism may play a role in IBS-D.

MitoLSDB: a comprehensive resource to study genotype to phenotype correlations in human mitochondrial DNA variations

Human mitochondrial DNA (mtDNA) encodes a set of 37 genes which are essential structural and functional components of the electron transport chain. Variations in these genes have been implicated in a broad spectrum of diseases and are extensively reported in literature and various databases. In this study, we describe MitoLSDB, an integrated platform to catalogue disease association studies on mtDNA (http://mitolsdb.igib.res.in). The main goal of MitoLSDB is to provide a central platform for direct submissions of novel variants that can be curated by the Mitochondrial Research Community. MitoLSDB provides access to standardized and annotated data from literature and databases encompassing information from 5231 individuals, 675 populations and 27 phenotypes. This platform is developed using the Leiden Open (source) Variation Database (LOVD) software. MitoLSDB houses information on all 37 genes in each population amounting to 132397 variants, 5147 unique variants. For each variant its genomic location as per the Revised Cambridge Reference Sequence, codon and amino acid change for variations in protein-coding regions, frequency, disease/phenotype, population, reference and remarks are also listed. MitoLSDB curators have also reported errors documented in literature which includes 94 phantom mutations, 10 NUMTs, six documentation errors and one artefactual recombination. MitoLSDB is the largest repository of mtDNA variants systematically standardized and presented using the LOVD platform. We believe that this is a good starting resource to curate mtDNA variants and will facilitate direct submissions enhancing data coverage, annotation in context of pathogenesis and quality control by ensuring non-redundancy in reporting novel disease associated variants.

Human genomic disease variants: a neutral evolutionary explanation

Many perspectives on the role of evolution in human health include nonempirical assumptions concerning the adaptive evolutionary origins of human diseases. Evolutionary analyses of the increasing wealth of clinical and population genomic data have begun to challenge these presumptions. In order to systematically evaluate such claims, the time has come to build a common framework for an empirical and intellectual unification of evolution and modern medicine. We review the emerging evidence and provide a supporting conceptual framework that establishes the classical neutral theory of molecular evolution (NTME) as the basis for evaluating disease- associated genomic variations in health and medicine. For over a decade, the NTME has already explained the origins and distribution of variants implicated in diseases and has illuminated the power of evolutionary thinking in genomic medicine. We suggest that a majority of disease variants in modern populations will have neutral evolutionary origins (previously neutral), with a relatively smaller fraction exhibiting adaptive evolutionary origins (previously adaptive). This pattern is expected to hold true for common as well as rare disease variants. Ultimately, a neutral evolutionary perspective will provide medicine with an informative and actionable framework that enables objective clinical assessment beyond convenient tendencies to invoke past adaptive events in human history as a root cause of human disease.

Large scale mtDNA sequencing reveals sequence and functional conservation as major determinants of homoplasmic mtDNA variant distribution

The mitochondrial DNA (mtDNA) is highly variable, containing large numbers of pathogenic mutations and neutral polymorphisms. The spectrum of homoplasmic mtDNA variation was characterized in 730 subjects and compared with known pathogenic sites. The frequency and distribution of variants in protein coding genes were inversely correlated with conservation at the amino acid level. Analysis of tRNA secondary structures indicated a preference of variants for the loops and some acceptor stem positions. This comprehensive overview of mtDNA variants distinguishes between regions and positions which are likely not critical, mainly conserved regions with pathogenic mutations and essential regions containing no mutations at all.

Phylomedicine: an evolutionary telescope to explore and diagnose the universe of disease mutations

Modern technologies have made the sequencing of personal genomes routine. They have revealed thousands of nonsynonymous (amino acid altering) single nucleotide variants (nSNVs) of protein-coding DNA per genome. What do these variants foretell about an individual's predisposition to diseases? The experimental technologies required to carry out such evaluations at a genomic scale are not yet available. Fortunately, the process of natural selection has lent us an almost infinite set of tests in nature. During long-term evolution, new mutations and existing variations have been evaluated for their biological consequences in countless species, and outcomes are readily revealed by multispecies genome comparisons. We review studies that have investigated evolutionary characteristics and in silico functional diagnoses of nSNVs found in thousands of disease-associated genes. We conclude that the patterns of long-term evolutionary conservation and permissible sequence divergence are essential and instructive modalities for functional assessment of human genetic variations.

MitoTool: a web server for the analysis and retrieval of human mitochondrial DNA sequence variations

MitoTool, a web-based bioinformatics platform, is designed for deciphering human mitochondrial DNA (mtDNA) data in batch mode. The platform has advantages in (i) parsing diverse types of mtDNA data; (ii) automatically classifying haplogroup according to mtDNA sequences or variants; (iii) discovering possibly missing variants of the samples with claimed haplogroups status; (iv) estimating the evolutionary conservation index, protein coding effect and potential pathogenicity of certain substitutions; (v) performing statistical analysis for haplogroup distribution frequency between case and control groups. Furthermore, it offers an integrated database for retrieving five types of mitochondrion-related information. The MitoTool is freely accessed at http://www.mitotool.org.