Automated sequencing of mitochondrial DNA

The Role of Mitochondrial DNA Variation in Drug Response: A Systematic Review

Background: The triad of drug efficacy, toxicity and resistance underpins the risk-benefit balance of all therapeutics. The application of pharmacogenomics has the potential to improve the risk-benefit balance of a given therapeutic via the stratification of patient populations based on DNA variants. A growth in the understanding of the particulars of the mitochondrial genome, alongside the availability of techniques for its interrogation has resulted in a growing body of literature examining the impact of mitochondrial DNA (mtDNA) variation upon drug response. Objective: To critically evaluate and summarize the available literature, across a defined period, in a systematic fashion in order to map out the current landscape of the subject area and identify how the field may continue to advance. Methods: A systematic review of the literature published between January 2009 and December 2020 was conducted using the PubMed database with the following key inclusion criteria: reference to specific mtDNA polymorphisms or haplogroups, a core objective to examine associations between mtDNA variants and drug response, and research performed using human subjects or human in vitro models. Results: Review of the literature identified 24 articles reporting an investigation of the association between mtDNA variant(s) and drug efficacy, toxicity or resistance that met the key inclusion criteria. This included 10 articles examining mtDNA variations associated with antiretroviral therapy response, 4 articles examining mtDNA variants associated with anticancer agent response and 4 articles examining mtDNA variants associated with antimicrobial agent response. The remaining articles covered a wide breadth of medications and were therefore grouped together and referred to as "other." Conclusions: Investigation of the impact of mtDNA variation upon drug response has been sporadic to-date. Collective assessment of the associations identified in the articles was inconclusive due to heterogeneous methods and outcomes, limited racial/ethnic groups, lack of replication and inadequate statistical power. There remains a high degree of idiosyncrasy in drug response and this area has the potential to explain variation in drug response in a clinical setting, therefore further research is likely to be of clinical benefit.

Deep-Coverage MPS Analysis of Heteroplasmic Variants within the mtGenome Allows for Frequent Differentiation of Maternal Relatives

Distinguishing between maternal relatives through mitochondrial (mt) DNA sequence analysis has been a longstanding desire of the forensic community. Using a deep-coverage, massively parallel sequencing (DCMPS) approach, we studied the pattern of mtDNA heteroplasmy across the mtgenomes of 39 mother-child pairs of European decent; haplogroups H, J, K, R, T, U, and X. Both shared and differentiating heteroplasmy were observed on a frequent basis in these closely related maternal relatives, with the minor variant often presented as 2–10% of the sequencing reads. A total of 17 pairs exhibited differentiating heteroplasmy (44%), with the majority of sites (76%, 16 of 21) occurring in the coding region, further illustrating the value of conducting sequence analysis on the entire mtgenome. A number of the sites of differentiating heteroplasmy resulted in non-synonymous changes in protein sequence (5 of 21), and to changes in transfer or ribosomal RNA sequences (5 of 21), highlighting the potentially deleterious nature of these heteroplasmic states. Shared heteroplasmy was observed in 12 of the 39 mother-child pairs (31%), with no duplicate sites of either differentiating or shared heteroplasmy observed; a single nucleotide position (16093) was duplicated between the data sets. Finally, rates of heteroplasmy in blood and buccal cells were compared, as it is known that rates can vary across tissue types, with similar observations in the current study. Our data support the view that differentiating heteroplasmy across the mtgenome can be used to frequently distinguish maternal relatives, and could be of interest to both the medical genetics and forensic communities.

Development and assessment of an optimized next-generation DNA sequencing approach for the mtgenome using the Illumina MiSeq

The development of molecular tools to detect and report mitochondrial DNA (mtDNA) heteroplasmy will increase the discrimination potential of the testing method when applied to forensic cases. The inherent limitations of the current state-of-the-art, Sanger-based sequencing, including constrictions in speed, throughput, and resolution, have hindered progress in this area. With the advent of next-generation sequencing (NGS) approaches, it is now possible to clearly identify heteroplasmic variants, and at a much lower level than previously possible. However, in order to bring these approaches into forensic laboratories and subsequently as accepted scientific information in a court of law, validated methods will be required to produce and analyze NGS data. We report here on the development of an optimized approach to NGS analysis for the mtDNA genome (mtgenome) using the Illumina MiSeq instrument. This optimized protocol allows for the production of more than 5 gigabases of mtDNA sequence per run, sufficient for detection and reliable reporting of minor heteroplasmic variants down to approximately 0.5-1.0% when multiplexing twelve samples. Depending on sample throughput needs, sequence coverage rates can be set at various levels, but were optimized here for at least 5000 reads. In addition, analysis parameters are provided for a commercially available software package that identify the highest quality sequencing reads and effectively filter out sequencing-based noise. With this method it will be possible to measure the rates of low-level heteroplasmy across the mtgenome, evaluate the transmission of heteroplasmy between the generations of maternal lineages, and assess the drift of variant sequences between different tissue types within an individual.

Controlling for contamination in re-sequencing studies with a reproducible web-based phylogenetic approach

Polymorphism discovery is a routine application of next-generation sequencing technology where multiple samples are sent to a service provider for library preparation, subsequent sequencing, and bioinformatic analyses. The decreasing cost and advances in multiplexing approaches have made it possible to analyze hundreds of samples at a reasonable cost. However, because of the manual steps involved in the initial processing of samples and handling of sequencing equipment, cross-contamination remains a significant challenge. It is especially problematic in cases where polymorphism frequencies do not adhere to diploid expectation, for example, heterogeneous tumor samples, organellar genomes, as well as during bacterial and viral sequencing. In these instances, low levels of contamination may be readily mistaken for polymorphisms, leading to false results. Here we describe practical steps designed to reliably detect contamination and uncover its origin, and also provide new, Galaxy-based, readily accessible computational tools and workflows for quality control. All results described in this report can be reproduced interactively on the web as described at http://usegalaxy.org/contamination.

The mitochondrial T1095C mutation increases gentamicin-mediated apoptosis

We have previously reported a heteroplasmic mtDNA mutation (T1095C) in the 12SrRNA gene of an Italian family with features of maternally-inherited parkinsonism, antibiotic-mediated deafness and peripheral neuropathy. In the present study, we demonstrate that a transmitochondrial cybrid line derived from the proband of this family shows selective depletion of mitochondrial glutathione and decreases in the activity of complex II/III. Moreover, when exposed to an aminoglycoside antibiotic these cells responded with a ten-fold increase in the number of apoptotic cells compared to controls. These results support a pathogenic role for the T1095C mutation and indicate that the mutation increases the risk for aminoglycoside-induced toxicity.

Schizophrenia: maternal inheritance and heteroplasmy of mtDNA mutations

Role of mitochondrial pathology in schizophrenia has not been fully clarified. We searched for distinctive variants in mtDNA extracted from the gray matter of postmortem brains and from peripheral blood samples. We screened mtDNA region containing 5 genes encoding subunits of cytochrome c oxidase and ATPases. Polymorphisms not already reported in databases are recorded as unregistered rare variants. Four unregistered, non-synonymous rare variants were detected in 4 schizophrenic samples. Seven registered non-synonymous variants were not previously detected in non-psychotic Japanese samples registered in the mtSNP database. These variants may contribute to disease pathophysiology. In one family, compound mutations showed co-segregation with schizophrenia. MtDNA mutations could confer a risk for schizophrenia in the Japanese population, although further analyses are needed.

Dynamics of mitochondrial heteroplasmy in three families investigated via a repeatable re-sequencing study

Background

Originally believed to be a rare phenomenon, heteroplasmy - the presence of more than one mitochondrial DNA (mtDNA) variant within a cell, tissue, or individual - is emerging as an important component of eukaryotic genetic diversity. Heteroplasmies can be used as genetic markers in applications ranging from forensics to cancer diagnostics. Yet the frequency of heteroplasmic alleles may vary from generation to generation due to the bottleneck occurring during oogenesis. Therefore, to understand the alterations in allele frequencies at heteroplasmic sites, it is of critical importance to investigate the dynamics of maternal mtDNA transmission.

Results

Here we sequenced, at high coverage, mtDNA from blood and buccal tissues of nine individuals from three families with a total of six maternal transmission events. Using simulations and re-sequencing of clonal DNA, we devised a set of criteria for detecting polymorphic sites in heterogeneous genetic samples that is resistant to the noise originating from massively parallel sequencing technologies. Application of these criteria to nine human mtDNA samples revealed four heteroplasmic sites.

Conclusions

Our results suggest that the incidence of heteroplasmy may be lower than estimated in some other recent re-sequencing studies, and that mtDNA allelic frequencies differ significantly both between tissues of the same individual and between a mother and her offspring. We designed our study in such a way that the complete analysis described here can be repeated by anyone either at our site or directly on the Amazon Cloud. Our computational pipeline can be easily modified to accommodate other applications, such as viral re-sequencing.

Macrophages that survive hyperoxia exposure have higher superoxide dismutase activities in their mitochondria

Prolonged exposure to hyperoxia, which is routinely used in patients with severe respiratory failure, leads to the generation of excessive reactive oxygen species, resulting in lung injury. In the present study, we focused on macrophages and their survival, superoxide dismutase (SOD) activity in mitochondria (Mn-SOD activity), and mitochondrial DNA (mtDNA) mutation after exposure to hyperoxia. Macrophages were cultured under two different conditions: normoxia and intermittent hyperoxia. The number of cells exposed to intermittent hyperoxia for 3 weeks significantly decreased, compared with the number of cells exposed to normoxia. The Mn-SOD activity of the cells that survived intermittent hyperoxia exposure was significantly higher than that of the cells exposed to normoxia. Direct sequencing and a PCR-RFLP assay did not provide any evidence of mutation in the cells that survived intermittent hyperoxia exposure. In conclusion, an increase in the antioxidative activity of mitochondria is important for the survival of macrophages exposed to hyperoxia, and the increased activity level possibly enhances protective effects against mtDNA mutations in surviving cells.

Distilling artificial recombinants from large sets of complete mtDNA genomes

Background

Large-scale genome sequencing poses enormous problems to the logistics of laboratory work and data handling. When numerous fragments of different genomes are PCR amplified and sequenced in a laboratory, there is a high immanent risk of sample confusion. For genetic markers, such as mitochondrial DNA (mtDNA), which are free of natural recombination, single instances of sample mix-up involving different branches of the mtDNA phylogeny would give rise to reticulate patterns and should therefore be detectable.

Methodology/Principal Findings

We have developed a strategy for comparing new complete mtDNA genomes, one by one, to a current skeleton of the worldwide mtDNA phylogeny. The mutations distinguishing the reference sequence from a putative recombinant sequence can then be allocated to two or more different branches of this phylogenetic skeleton. Thus, one would search for two (or three) near-matches in the total mtDNA database that together best explain the variation seen in the recombinants. The evolutionary pathway from the mtDNA tree connecting this pair together with the recombinant then generate a grid-like median network, from which one can read off the exchanged segments.

Conclusions

We have applied this procedure to a large collection of complete human mtDNA sequences, where several recombinants could be distilled by our method. All these recombinant sequences were subsequently corrected by de novo experiments – fully concordant with the predictions from our data-analytical approach.

Identification of mitochondrial DNA polymorphisms that alter mitochondrial matrix pH and intracellular calcium dynamics

Mitochondrial DNA (mtDNA) is highly polymorphic, and its variations in humans may contribute to individual differences in function as well as susceptibility to various diseases such as Parkinson disease, Alzheimer disease, bipolar disorder, and cancer. However, it is unclear whether and how mtDNA polymorphisms affect intracellular function, such as calcium signaling or pH regulation. Here we searched for mtDNA polymorphisms that have intracellular functional significance using transmitochondrial hybrid cells (cybrids) carrying ratiometric Pericam (RP), a fluorescent calcium indicator, targeted to the mitochondria and nucleus. By analyzing the entire mtDNA sequence in 35 cybrid lines, we found that two closely linked nonsynonymous polymorphisms, 8701A and 10398A, increased the basal fluorescence ratio of mitochondria-targeted RP. Mitochondrial matrix pH was lower in the cybrids with 8701A/10398A than it was in those with 8701G/10398G, suggesting that the difference observed by RP was mainly caused by alterations in mitochondrial calcium levels. Cytosolic calcium response to histamine also tended to be higher in the cybrids with 8701A/10398A. It has previously been reported that 10398A is associated with an increased risk of Parkinson disease, Alzheimer disease, bipolar disorder, and cancer, whereas 10398G associates with longevity. Our findings suggest that these mtDNA polymorphisms may play a role in the pathophysiology of these complex diseases by affecting mitochondrial matrix pH and intracellular calcium dynamics.

Mitochondria play important roles in energy production and regulation of intracellular calcium levels. Mitochondria have their own genetic material, mitochondrial DNA (mtDNA). In spite of its short length (16 kbp), mtDNA is highly variable among individuals and is thought to contribute to interindividual functional variability in energy-requiring activities such as intelligence and athletic performance. However, it is unclear whether mtDNA polymorphisms affect intracellular function and condition. Using transmitochondrial hybrid cells, the authors found two closely linked mtDNA polymorphisms, 10398A/G and 8701A/G, which cause alterations in mitochondrial pH and calcium concentration. Cytosolic calcium response to histamine tended to be different between transmitochondrial hybrid cells carrying these two mtDNA polymorphisms. It has been reported that the 10398A mtDNA polymorphism is a risk factor for Parkinson disease, Alzheimer disease, cancer, and bipolar disorder, whereas 10398G is associated with longevity. The present findings suggest that these mtDNA polymorphisms may play a role in the pathophysiology of these complex diseases by affecting mitochondrial matrix pH and intracellular calcium dynamics.

Phantom mutation hotspots in human mitochondrial DNA

Phantom mutations are systematic artifacts generated in the course of the sequencing process. Contra common belief these artificial mutations are nearly ubiquitous in sequencing results, albeit at frequencies that may vary dramatically. The amount of artifacts depends not only on the sort of automated sequencer and sequencing chemistry employed, but also on other lab-specific factors. An experimental study executed on four samples under various combinations of sequencing conditions revealed a number of phantom mutations occurring at the same sites of mitochondrial DNA (mtDNA) repeatedly. To confirm these and identify further hotspots for artifacts, > 5000 mtDNA electropherograms were screened for artificial patterns. Further, > 30 000 published hypervariable segment I sequences were compared at potential hotspots for phantom mutations, especially for variation at positions 16085 and 16197. Resequencing of several samples confirmed the artificial nature of these and other polymorphisms in the original publications. Single-strand sequencing, as typically executed in medical and anthropological studies, is thus highly vulnerable to this kind of artifacts. In particular, phantom mutation hotspots could easily lead to misidentification of somatic mutations and to misinterpretations in all kinds of clinical mtDNA studies.

Mitochondrial genome polymorphisms associated with type-2 diabetes or obesity

Although a strong correlation between type-2 diabetes and obesity has been found, no comparative analysis between diabetes and obesity has been performed with respect to mitochondrial DNA (mtDNA) polymorphisms. To test the hypothesis that certain mitochondrial single nucleotide polymorphisms (mtSNPs) might be associated with obesity or type-2 diabetes, we determined the entire sequences of the mitochondrial genomes from 96 patients with type-2 diabetes and those from 96 young obese adults by direct sequencing and compared the frequencies of mtSNPs between these two groups. A mtSNP, 8684C > T (T53I) in the mitochondrial ATP synthase subunit 6 gene (ATP6), was detected in 5 of the 96 patients with type-2 diabetes, whereas this substitution was not detected in any of the 96 young obese adults. Two mtSNPs, 3497C > T (A64V) in NADH dehydrogenase subunit 1 gene (ND1) and 1119T > C (472U > C) in the 12S rRNA gene, were detected in 5 of the 96 young obese adults, whereas these substitutions were not detected in any of the 96 diabetic patients. The 8684C > T transition associated with type-2 diabetes represents haplogroup M8a, and the 3497C > T and 1119T > C transitions predisposing to obesity represent haplogroup B4c. These results suggest that distinct mtSNPs contribute to susceptibility to type-2 diabetes or obesity, pointing out the necessity of large-scale case control studies.

Mitochondrial DNA sequence analysis of patients with 'atypical psychosis'

Although classical psychopathological studies have shown the presence of an independent diagnostic category, 'atypical psychosis', most psychotic patients are currently classified into two major diagnostic categories, schizophrenia and bipolar disorder, by the Diagnostic and Statistical Manual of Mental Disorders (4th edn; DSM-IV) criteria. 'Atypical psychosis' is characterized by acute confusion without systematic delusion, emotional instability, and psychomotor excitement or stupor. Such clinical features resemble those seen in organic mental syndrome, and differential diagnosis is often difficult. Because patients with mitochondrial myopathy, encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) sometimes show organic mental disorder, 'atypical psychosis' may be caused by mutations of mitochondrial DNA (mtDNA) in some patients. In the present study whole mtDNA was sequenced for seven patients with various psychotic disorders, who could be categorized as 'atypical psychosis'. None of them had known mtDNA mutations pathogenic for mitochondrial encephalopathy. Two of seven patients belonged to a subhaplogroup F1b1a with low frequency. These results did not support the hypothesis that clinical presentation of some patients with 'atypical psychosis' is a reflection of subclinical mitochondrial encephalopathy. However, the subhaplogroup F1b1a may be a good target for association study of 'atypical psychosis'.

Mitochondrial genome variation in eastern Asia and the peopling of Japan

To construct an East Asia mitochondrial DNA (mtDNA) phylogeny, we sequenced the complete mitochondrial genomes of 672 Japanese individuals (http://www.giib.or.jp/mtsnp/index_e.html). This allowed us to perform a phylogenetic analysis with a pool of 942 Asiatic sequences. New clades and subclades emerged from the Japanese data. On the basis of this unequivocal phylogeny, we classified 4713 Asian partial mitochondrial sequences, with <10% ambiguity. Applying population and phylogeographic methods, we used these sequences to shed light on the controversial issue of the peopling of Japan. Population-based comparisons confirmed that present-day Japanese have their closest genetic affinity to northern Asian populations, especially to Koreans, which finding is congruent with the proposed Continental gene flow to Japan after the Yayoi period. This phylogeographic approach unraveled a high degree of differentiation in Paleolithic Japanese. Ancient southern and northern migrations were detected based on the existence of basic M and N lineages in Ryukyuans and Ainu. Direct connections with Tibet, parallel to those found for the Y-chromosome, were also apparent. Furthermore, the highest diversity found in Japan for some derived clades suggests that Japan could be included in an area of migratory expansion to Continental Asia. All the theories that have been proposed up to now to explain the peopling of Japan seem insufficient to accommodate fully this complex picture.

Golden mean to longevity: rareness of mitochondrial cytochrome b variants in centenarians but not in patients with Parkinson's disease

To test the hypothesis that centenarians are free from deleterious mitochondrial variations, we analyzed amino acid variations in the cytochrome b molecule of 64 Japanese centenarians. Although the frequencies of some variations, such as N260D and G251S, differed significantly between centenarians and patients with Parkinson's disease, the most striking feature of centenarian cytochrome b was the rareness of amino acid variations in contrast to the variety of amino acid replacements in patients with Parkinson's disease. These results suggest that centenarians are genetically hitting the "golden mean" (less variation from the consensus cytochrome b sequence or less mismatch with other subunits). A multiplex detection system for various deleterious variations in combination with genetic tests for longevity-associated genotypes will be necessary to predict longevity or age-related diseases.

Mitochondrial ATPase subunit 6 and cytochrome B gene polymorphisms in young obese adults

We hypothesized that the mutational strand asymmetry is more strongly exerted upon the mitochondrial cytochrome b (Cytb) gene, which is distant from the origin of the light-strand replication (Ori(L)), than upon the ATPase subunit 6 (ATP6) gene, which is close to the Ori(L). To test this hypothesis, we determined the sequences of these two genes in 96 Japanese young obese adults. The frequency of G-->A transitions was significantly higher than that of C-->T transitions in the Cytb gene, whereas the frequencies of G-->A and C-->T transitions were not significantly different in the ATP6 gene. The marked mutational strand asymmetry in the Cytb gene can be explained by the deamination of C to uracil in the long single-stranded state of the heavy strand during replication. The ratio of the nonsynonymous substitutions at the second codon positions to those at the first codon positions was significantly lower in the Cytb gene than in the ATP6 gene. The physicochemical differences between the standard and the replaced amino acid residues were significantly smaller in the Cytb gene than in ATP6 one. The present study indicates that amino acid sequences are less variable for Cytb than for ATP6 in spite of the strong mutational strand asymmetry for the Cytb gene.

Mitochondrial DNA determines the cellular response to cancer therapeutic agents

Mutations in the mitochondrial genome leading to mitochondrial dysfunction have been reported in a variety of cancers. However, the potential implication of these findings in the cellular response to cancer therapeutic agents is unclear. To examine the importance of mitochondrial DNA (mitDNA) encoded functions in cancer therapeutic response, we determined the clonogenic survival of HSL2 (Rho+, HeLa subline), and its derivative cell line lacking mitDNA (Rho0) after exposure to different anticancer agents. We found that isogenic Rho0 cells lacking mitDNA were extremely resistant to adriamycin and photodynamic therapy (PDT) induced cell death, whereas the Rho+ cell line was sensitive. However, there was no measurable difference in the responses of these cell lines to either alkylating agent or gamma-radiation. We show that the development of resistance to adriamycin was not due to changes in apoptotic cell death, cell cycle response or to the uptake of adriamycin in isogenic Rho0 cells. We also demonstrate that exposure of HeLa cells to adriamycin leads to mutations in mitDNA. These studies provide direct evidence that mitDNA plays an important role in cellular sensitivity to cancer therapeutic agents.

Mitochondrial DNA mutations and age

Apopotic cell death is reported to be prominent in the stable tissues of the failing heart, in cardiomyopathies (CM), in the sinus node of complete heart block, in B cells of diabetes mellitus, and in neurodegenerative diseases. Recently, mitochondrial (mt) control of nuclear apoptosis was demonstrated in the cell-free system. The mt bioenergetic crisis induced by exogenously added factors such as respiratory inhibitors leads to the collapse of mt transmembrane potential, to the opening of the inner membrane pore, to the release of the apoptotic protease activating factors into cytosol, and subsequently to nuclear DNA fragmentation. However, the endogenous factor for the mt bioenegertic crisis in naturally occurring cell death under the physiological conditions without vascular involvement has remained unknown. Recently devised, the total detection system for deletion demonstrates the extreme fragmentation of mtDNA in the cardiac myocytes of senescence, and mt CM harboring maternally inherited point mutations in mtDNA and on the cultured cell line with or without mtDNA disclosed that mtDNA is unexpectedly fragile to hydroxyl radial damage and hence to oxygen stress. The great majority of wild-type mtDNA fragmented into over two hundreds types of deleted mtDNA related to oxidative damage, resulting in pleioplasmic defects in the mt energy transducing system. The mtDNA fragmentation to this level is demonstrated in cardiac myocytes of normal subjects over age 80, of an mtCM patient who died at age 20 and one who died at age 19, of a recipient of heart transplantation at age 7 with severe mtCM, and in mtDNA of a cultured cell line under hyperbaric oxygen stress for two days, leading a majority of cells to apoptotic death on the third day. The extreme fragility of mtDNA could be the missing link in the apoptosis cascade that is the physiological basis of aging and geriatrics of such stable tissues as nerve and muscle.