Novel non‐neutral mitochondrial DNA mutations found in childhood acute lymphoblastic leukemia

Mitochondria produce adenosine triphosphate (ATP) for energy requirements via the mitochondrial oxidative phosphorylation (OXPHOS) system. One of the hallmarks of cancer is the energy shift toward glycolysis. Low OXPHOS activity and increased glycolysis are associated with aggressive types of cancer. Mitochondria have their own genome (mitochondrial DNA [mtDNA]) encoding for 13 essential subunits of the OXPHOS enzyme complexes. We studied mtDNA in childhood acute lymphoblastic leukemia (ALL) to detect potential pathogenic mutations in OXPHOS complexes. The whole mtDNA from blood and bone marrow samples at diagnosis and follow‐up from 36 ALL patients were analyzed. Novel or previously described pathogenic mtDNA mutations were identified in 8 out of 36 patients. Six out of these 8 patients had died from ALL. Five out of 36 patients had an identified poor prognosis genetic marker, and 4 of these patients had mtDNA mutations. Missense or nonsense mtDNA mutations were detected in the genes encoding subunits of OXPHOS complexes, as follows: MT‐ND1, MT‐ND2, MT‐ND4L and MT‐ND6 of complex I; MT‐CO3 of complex IV; and MT‐ATP6 and MT‐ATP8 of complex V. We discovered mtDNA mutations in childhood ALL supporting the hypothesis that non‐neutral variants in mtDNA affecting the OXPHOS function may be related to leukemic clones.     

Deciphering exome sequencing data: Bringing mitochondrial DNA variants to light

The expanding use of exome sequencing (ES) in diagnosis generates a huge amount of data, including untargeted mitochondrial DNA (mtDNA) sequences. We developed a strategy to deeply study ES data, focusing on the mtDNA genome on a large unspecific cohort to increase diagnostic yield. A targeted bioinformatics pipeline assembled mitochondrial genome from ES data to detect pathogenic mtDNA variants in parallel with the “in‐house” nuclear exome pipeline. mtDNA data coming from off‐target sequences (indirect sequencing) were extracted from the BAM files in 928 individuals with developmental and/or neurological anomalies. The mtDNA variants were filtered out based on database information, cohort frequencies, haplogroups and protein consequences. Two homoplasmic pathogenic variants (m.9035T>C and m.11778G>A) were identified in 2 out of 928 unrelated individuals (0.2%): the m.9035T>C (MT‐ATP6) variant in a female with ataxia and the m.11778G>A (MT‐ND4) variant in a male with a complex mosaic disorder and a severe ophthalmological phenotype, uncovering undiagnosed Leber's hereditary optic neuropathy (LHON). Seven secondary findings were also found, predisposing to deafness or LHON, in 7 out of 928 individuals (0.75%). This study demonstrates the usefulness of including a targeted strategy in ES pipeline to detect mtDNA variants, improving results in diagnosis and research, without resampling patients and performing targeted mtDNA strategies.     

Capture of Somatic mtDNA Point Mutations with Severe Effects on Oxidative Phosphorylation in Synaptosome Cybrid Clones from Human Brain

Mitochondrial DNA (mtDNA) is replicated throughout life in postmitotic cells, resulting in higher levels of somatic mutation than in nuclear genes. However, controversy remains as to the importance of low‐level mtDNA somatic mutants in cancerous and normal human tissues. To capture somatic mtDNA mutations for functional analysis, we generated synaptosome cybrids from synaptic endings isolated from fresh hippocampus and cortex brain biopsies. We analyzed the whole mtDNA genome from 120 cybrid clones derived from four individual donors by chemical cleavage of mismatch and Sanger sequencing, scanning around two million base pairs. Seventeen different somatic point mutations were identified, including eight coding region mutations, four of which result in frameshifts. Examination of one cybrid clone with a novel m.2949_2953delCTATT mutation in MT‐RNR2 (which encodes mitochondrial 16S rRNA) revealed a severe disruption of mtDNA‐encoded protein translation. We also performed functional studies on a homoplasmic nonsense mutation in MT‐ND1, previously reported in oncocytomas, and show that both ATP generation and the stability of oxidative phosphorylation complex I are disrupted. As the mtDNA remains locked against direct genetic manipulation, we demonstrate that the synaptosome cybrid approach can capture biologically relevant mtDNA mutants in vitro to study effects on mitochondrial respiratory chain function.     

An evaluation of polymorphisms in casein kinase 1 delta and epsilon genes in major psychiatric disorders

Disturbances of the circadian rhythm are involved in the pathophysiology of bipolar disorder (BD), schizophrenia (SCZ) and major depressive disorder (MDD). Specifically, because clock gene dysfunction is good candidate for enhancing the susceptibility to these psychiatric disorders, we selected two circadian rhythm-related genes (CSNK1D and CSNK1E) and investigated genetic associations of the genes with these three disorders.     

Mitochondrial DNA Polymerase γ Mutations and Their Implications in mtDNA Alterations in Colorectal Cancer

Mitochondrial DNA was found to be highly mutated in colorectal cancer cells. One of the key molecules involved in the maintenance of the mitochondrial genome is the nuclear‐encoded polymerase gamma. The aim of our study was to determine if there is a link between polymorphisms within the polymerase gamma gene (POLG) and somatic mutations within the mitochondrial genome in cancer cells. We investigated POLG sequence variability in 50 colorectal cancer patients whose complete mitochondrial genome sequences were determined. Relative mtDNA copy number was also determined. We identified 251 sequence variants in the POLG gene. Most of them were germline‐specific (～92%). Twenty‐one somatic changes in POLG were found in 10 colorectal cancer patients. We have found no association between the occurrence of mtDNA somatic mutations and the somatically occurring variants in POLG. MtDNA content was reduced in patients carrying somatic variants in POLG or germline nucleotide variants located in the region encoding the POLG polymerase domain, but the difference did not reach statistical significance. Our findings suggest that somatic mtDNA mutations occurring in colorectal cancer are not a consequence of somatic mutations in POLG. Nevertheless, POLG nucleotide variants may lead to a decrease in mtDNA content, and consequently result in mitochondrial dysfunction.     

False positive results of mitochondrial DNA depletion/deletion due to single nucleotide substitutions causing appearance of additional PvuII restriction sites.

Human mitochondrial diseases are usually caused by dysfunction of mitochondrial DNA (mtDNA), particularly by point mutations, deletions, or depletions. In commonly used procedures for molecular diagnostics of mitochondrial dysfunction, one of the first steps is linearization of circular mitochondrial genomes with either BamHI or PvuII restriction endonulease, which cuts human mtDNA at a unique site. Here, we describe a case of false positive results, which suggested mtDNA depletion or a large deletion in a patient's tissue sample. More detailed analysis (mtDNA sequencing) revealed that these false positive results were caused by the presence of the 12753A>G substitution in the gene coding for NADH dehydrogenase subunit 5 (ND5). This substitution results in no change in amino acid sequence of the gene product but creates an additional PvuII site. Investigating a population of 200 patients not affected by mitochondrial diseases, we found an additional case of 12753A>G, and also another substitution, 12804T>C, which also results in no change in amino acid sequence of ND5 but creates an additional PvuII site. A few cases of 12753A>G and 12804T>C substitutions were found previously in Asian, American, African, and European populations (though they were not reported to date in the MITOMAP), but those samples were used in population studies and not tested for mtDNA deletion or depletion. Therefore, we present a cautionary report indicating that these mtDNA polymorphisms exist in various human populations (and thus, they are panethnic) and may cause false positive results of standard molecular analyses, including molecular diagnostics, of human mtDNA.     

The mitochondrial DNA 4977‐bp deletion and copy number alteration in Han Chinese samples with uterine fibroids

Uterine fibroids (UFs) are the most common benign neoplasms, but their pathogenesis is not completely understood. Thus far, alterations in the mitochondrial DNA (mtDNA) content and the mtDNA 4977‐bp deletion level in UFs, as well as the corresponding nontumorous tissue, have remained elusive. To test whether large mtDNA deletions and mtDNA content are involved in the pathogenesis of UFs, a total of 309 UF tissues and 28 paired adjacent myometrium from 270 UF patients were enrolled for the analysis of large mtDNA deletions and mtDNA content through the use of nested PCR and qPCR techniques, respectively. In our samples, a 4977‐bp deletion was identified: 36 out of 309 UF tissues (11.56%) and 15 out of 28 (53.57%) paired adjacent myometrium were detected to harbor the 4977‐bp deletion. In addition, a novel 4838‐bp mtDNA deletion was identified in three UF tissues, and other different sizes of deleted fragments (4910, 4926, 5135‐bp) were also found in UFs for the first time. Furthermore, older age was significantly associated with an mtDNA large deletion in the paired adjacent myometrium. We also found that increased mtDNA content and higher expression of ND1 occurred in solitary fibroids compared to adjacent myometrium. In conclusion, we identified a lower frequency of mtDNA large deletions and some novel large deletion in UFs for the first time. Furthermore, there was a general increase of mtDNA copy number during solitary UF development. Although the definite mechanism by which mtDNA was altered is supposed to be further confirmed, it will be helpful for further studies on the pathological mechanism of UFs.     

A Novel in‐Frame 18‐bp Microdeletion in MT‐CYB Causes a Multisystem Disorder with Prominent Exercise Intolerance

A novel heteroplasmic mitochondrial DNA (mtDNA) microdeletion affecting the cytochrome b gene (MT‐CYB) was identified in an Italian female patient with a multisystem disease characterized by sensorineural deafness, cataracts, retinal pigmentary dystrophy, dysphagia, postural and gait instability, and myopathy with prominent exercise intolerance. The deletion is 18‐base pair long and encompasses nucleotide positions 15,649–15,666, causing the loss of six amino acids (Ile‐Leu‐Ala‐Met‐Ile‐Pro) in the protein, but leaving the remaining of the MT‐CYB sequence in frame. The defective complex III function was cotransferred with mutant mtDNA in cybrids, thus unequivocally establishing its pathogenic role. Maternal relatives failed to show detectable levels of the deletion in blood and urinary epithelium, suggesting a de novo mutational event. This is the second report of an in‐frame intragenic deletion in MT‐CYB, which most likely occurred in early stages of embryonic development, associated with a severe multisystem disorder with prominent exercise intolerance.     

Evidence for additionally increased apoptosis in the peripheral blood mononuclear cells of major depressive patients with a high risk for suicide

Several studies have suggested a pathophysiological role of blood cell apoptosis in major depressive disorder (MDD). The aim of this study was to evaluate mRNA expression levels of Bcl‐2, Bax, and Fas in peripheral blood mononuclear cells (PBMCs) of MDD patients with a high risk for suicide relative to those without a high risk for suicide as well as healthy subjects. The mRNA expression of Bcl‐2, Bax, and Fas as well as the Bcl‐2/Bax ratio was examined in the PBMCs of 30 MDD patients with a high risk for suicide, 30 MDD patients without a high risk for suicide, and 30 healthy controls. The mRNA expression of target genes was measured using real‐time quantitative Polymerase Chain Reaction (PCR). FAS mRNA expression was significantly increased, and Bcl‐2 mRNA expression and the Bcl‐2/Bax expression ratio were significantly decreased, in the PBMCs of MDD patients with or without a high risk for suicide attempts compared to healthy controls (p < .001). However, Bax mRNA expression was significantly increased only in MDD patients with a high risk for suicide. Moreover, MDD patients with a high risk for suicide had increased Bax and FAS mRNA expression and decreased Bcl‐2 and Bcl‐2/Bax ratio when compared to patients without risk for suicide (p < .001). Our findings may support the role of both internal and external apoptotic pathways in the interplay between the immune system and depressive symptoms, especially in patients with a high risk for suicide.     

Association of genetic variation in CACNA1C with bipolar disorder in Han Chinese

Background

A growing body of evidence highlights the existence of shared genetic susceptibility to both major depressive disorder (MDD) and bipolar disorder (BD), suggesting some potential genetic overlap between the disorders. Genome-wide association studies have identified consistent association of single nucleotide polymorphisms of the α-1 C subunit of the L-type voltage-gated calcium channel gene (CACNA1C) with MDD and BD, suggesting CACNA1C as a promising candidate gene for susceptibility to mood disorders. In the present study, we tested the association of CACNA1C with MDD and BD in Han Chinese.

Methods

We genotyped three potentially functional polymorphisms in 635 MDD patients, 286 BD patients and 730 normal, control patients.

Results

The genotype frequencies of SNP rs1051375 showed statistically significant differences between the BD and control groups (P=0.005). At the allele level, the difference of G allele frequency of rs1051375 between BD patients and control subjects was also significant (P=0.011; OR=1.30, 95% CI: 1.06–1.58). We found that GG genotype of rs1051375 carriers had a lower age at onset than those with the AG or AA genotype, and the mean±standard deviation ages at onset of GG, AG and AA carriers were 24.04±4.22, 25.76±4.75 and 25.78±4.33 years, respectively. Neither genotype nor allele frequencies of the three polymorphisms were found to be significantly different between the MDD patients and control subjects.

Limitations

The relative small sample size in BD group should be considered a limitation of this study.

Conclusions

Our initial findings support a potential association of CACNA1C as a genetic risk factor for BD susceptibility.     

Mitochondrial genome architecture in non‐alcoholic fatty liver disease

Non‐alcoholic fatty liver disease (NAFLD) is associated with mitochondrial dysfunction, a decreased liver mitochondrial DNA (mtDNA) content, and impaired energy metabolism. To understand the clinical implications of mtDNA diversity in the biology of NAFLD, we applied deep‐coverage whole sequencing of the liver mitochondrial genomes. We used a multistage study design, including a discovery phase, a phenotype‐oriented study to assess the mutational burden in patients with steatohepatitis at different stages of liver fibrosis, and a replication study to validate findings in loci of interest. We also assessed the potential protein‐level impact of the observed mutations. To determine whether the observed changes are tissue‐specific, we compared the liver and the corresponding peripheral blood entire mitochondrial genomes. The nuclear genes POLG and POLG2 (mitochondrial DNA polymerase‐γ) were also sequenced. We observed that the liver mtDNA of patients with NAFLD harbours complex genomes with a significantly higher mutational (1.28‐fold) rate and degree of heteroplasmy than in controls. The analysis of liver mitochondrial genomes of patients with different degrees of fibrosis revealed that the disease severity is associated with an overall 1.4‐fold increase in mutation rate, including mutations in genes of the oxidative phosphorylation (OXPHOS) chain. Significant differences in gene and protein expression patterns were observed in association with the cumulative number of OXPHOS polymorphic sites. We observed a high degree of homology (～98%) between the blood and liver mitochondrial genomes. A missense POLG p.Gln1236His variant was associated with liver mtDNA copy number. In conclusion, we have demonstrated that OXPHOS genes contain the highest number of hotspot positions associated with a more severe phenotype. The variability of the mitochondrial genomes probably originates from a common germline source; hence, it may explain a fraction of the ‘missing heritability’ of NAFLD. Copyright © 2016 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.     

A tRNASer(UCN) gene in Artemia salina mitochondrial DNA: a case of mistaken identity

We have sequenced a segment of mitochondrial DNA (mtDNA) of a crustacean, the brine shrimp, Artemia salina, that includes 3 end-proximal regions of the genes for subunit 1 of the NADH dehydrogenase complex (ND1) and cytochrome b (Cyt b). From our data we conclude that in this mtDNA, as in the mtDNAs of Drosophila species, a tRNA^Ser(UCN) gene separates the ND1 and Cyt b genes. This is contrary to an earlier report that the A. salina ND1 and Cyt b genes are immediately adjacent to each other.     

Whole‐Exome Sequencing Identifies a Variant of the Mitochondrial MT‐ND1 Gene Associated with Epileptic Encephalopathy: West Syndrome Evolving to Lennox–Gastaut Syndrome

We describe a West syndrome (WS) patient with unidentified etiology that evolved to Lennox–Gastaut syndrome. The mitochondrial respiratory chain of the patient showed a simple complex I deficiency in fibroblasts. Whole‐exome sequencing (WES) uncovered two heterozygous mutations in NDUFV2 gene that were reassigned to a pseudogene. With the WES data, it was possible to obtain whole mitochondrial DNA sequencing and to identify a heteroplasmic variant in the MT‐ND1 (MTND1) gene (m.3946G>A, p.E214K). The expression of the gene in patient fibroblasts was not affected but the protein level was significantly reduced, suggesting that protein stability was affected by this mutation. The lower protein level also affected assembly of complex I and supercomplexes (I/III₂/IV and I/III₂), leading to complex I deficiency. While ATP levels at steady state under stress conditions were not affected, the amount of ROS produced by complex I was significantly increased.     

Memory mood congruency phenomenon in bipolar I disorder and major depression disorder patients

The objective of the present study was to evaluate memory performance in tasks with and without affective content (to confirm the mood congruency phenomenon) in acutely admitted patients with bipolar I disorder (BD) and major depression disorder (MDD) and in healthy participants. Seventy-eight participants (24 BD, 29 MDD, and 25 healthy controls) were evaluated. Three word lists were used as the memory task with affective content (positive, negative and indifferent). Psychiatric symptoms were also evaluated with rating scales (Young Mania Rating Scale for mania and Hamilton Depression Rating Scale for depression). Patients were selected during the first week of hospitalization. BD patients showed higher scores in the word span with positive tone than MDD patients and healthy controls (P = 0.002). No other difference was observed for tests with affective tone. MDD patients presented significantly lower scores in the Mini-Mental State Exam, logical memory test, visual recognition span, and digit span, while BD patients presented lower scores in the visual recognition test and digit span. Mood congruency effect was found for word span with positive tone among BD patients but no similar effect was observed among MDD patients for negative items. MDD patients presented more memory impairment than BD patients, but BD patients also showed memory impairment.     

No association of genetic variants in BDNF with major depression: A meta‐ and gene‐based analysis

Major depressive disorder (MDD) is a complex psychiatric condition with strong genetic predisposition. The association of MDD with genetic polymorphisms, such as Val66Met (rs6265), in the brain derived neurotrophic factor (BDNF), have been reported in many studies and the results were conflicting. In this study, we performed a systematic literature search and conducted random‐effects meta‐analysis to evaluate genetic variants in BDNF with MDD. A gene‐based analysis was also conducted to investigate the cumulative effects of genetic polymorphisms in BDNF. A total of 28 studies from 26 published articles were included in our analysis. Meta‐analysis yielded an estimated odds ratio (OR) of 0.96 (95% CI: 0.89–1.05; P = 0.402) for Val66Met (rs6265), 0.83 (95% CI: 0.67–1.04; P = 0.103) for 11757C/G, 1.16 (95% CI: 0.74–1.82; P = 0.527) for 270T/C, 1.03 (95% CI: 0.18–5.75; P = 0.974) for 712A/G and 0.98 (95% CI: 0.85–1.14; P = 0.831) for rs988748. The gene‐based analysis indicated that BDNF is not associated with MDD (P > 0.21). Our updated meta‐ and novel gene‐based analyses provide no evidence of the association of BDNF with major depression. © 2012 Wiley Periodicals, Inc.     

Analysis of Heteroplasmic Variants in the Cardiac Mitochondrial Genome of Individuals with Down Syndrome

Individuals with Down syndrome (DS, trisomy 21) exhibit a pro‐oxidative cellular environment as well as mitochondrial dysfunction. Increased oxidative stress may damage the mitochondrial DNA (mtDNA). The coexistence of mtDNA variants in a cell or tissue (i.e., heteroplasmy) may contribute to mitochondrial dysfunction. Given the evidence on mitochondrial dysfunction and the relatively high incidence of multiorganic disorders associated with DS, we hypothesized that cardiac tissue from subjects with DS may exhibit higher frequencies of mtDNA variants in comparison to cardiac tissue from donors without DS. This study documents the analysis of mtDNA variants in heart tissue samples from donors with (n = 12) and without DS (n = 33) using massively parallel sequencing. Contrary to the original hypothesis, the study's findings suggest that the cardiac mitochondrial genomes from individuals with and without DS exhibit many similarities in terms of (1) total number of mtDNA variants per sample, (2) the frequency of mtDNA variants, (3) the type of mtDNA variants, and (4) the patterns of distribution of mtDNA variants. In both groups of samples, the mtDNA control region showed significantly more heteroplasmic variants in comparison to the number of variants in protein‐ and RNA‐coding genes (P < 1.00×10^?4, ANOVA).     

Human Mitochondrial Cytochrome b Variants Studied in Yeast: Not All Are Silent Polymorphisms

Variations in mitochondrial DNA (mtDNA) cytochrome b (mt‐cyb) are frequently found within the healthy population, but also occur within a spectrum of mitochondrial and common diseases. mt‐cyb encodes the core subunit (MT‐CYB) of complex III, a central component of the oxidative phosphorylation system that drives cellular energy production and homeostasis. Despite significant efforts, most mt‐cyb variations identified are not matched with corresponding biochemical data, so their functional and pathogenic consequences in humans remain elusive. While human mtDNA is recalcitrant to genetic manipulation, it is possible to introduce human‐associated point mutations into yeast mtDNA. Using this system, we reveal direct links between human mt‐cyb variations in key catalytic domains of MT‐CYB and significant changes to complex III activity or drug sensitivity. Strikingly, m.15257G>A (p.Asp171Asn) increased the sensitivity of yeast to the antimalarial drug atovaquone, and m.14798T>C (p.Phe18Leu) enhanced the sensitivity of yeast to the antidepressant drug clomipramine. We demonstrate that while a small number of mt‐cyb variations had no functional effect, others have the capacity to alter complex III properties, suggesting they could play a wider role in human health and disease than previously thought. This compendium of new mt‐cyb‐biochemical relationships in yeast provides a resource for future investigations in humans.     

Fine Time Scaling of Purifying Selection on Human Nonsynonymous mtDNA Mutations Based on the Worldwide Population Tree and Mother–Child Pairs

A high‐resolution mtDNA phylogenetic tree allowed us to look backward in time to investigate purifying selection. Purifying selection was very strong in the last 2,500 years, continuously eliminating pathogenic mutations back until the end of the Younger Dryas (～11,000 years ago), when a large population expansion likely relaxed selection pressure. This was preceded by a phase of stable selection until another relaxation occurred in the out‐of‐Africa migration. Demography and selection are closely related: expansions led to relaxation of selection and higher pathogenicity mutations significantly decreased the growth of descendants. The only detectible positive selection was the recurrence of highly pathogenic nonsynonymous mutations (m.3394T>C‐m.3397A>G‐m.3398T>C) at interior branches of the tree, preventing the formation of a dinucleotide STR (TATATA) in the MT‐ND1 gene. At the most recent time scale in 124 mother–children transmissions, purifying selection was detectable through the loss of mtDNA variants with high predicted pathogenicity. A few haplogroup‐defining sites were also heteroplasmic, agreeing with a significant propensity in 349 positions in the phylogenetic tree to revert back to the ancestral variant. This nonrandom mutation property explains the observation of heteroplasmic mutations at some haplogroup‐defining sites in sequencing datasets, which may not indicate poor quality as has been claimed.