Mitochondrial Variants in Schizophrenia, Bipolar Disorder, and Major Depressive Disorder

Abstract

Background

Mitochondria provide most of the energy for brain cells by the process of oxidative phosphorylation. Mitochondrial abnormalities and deficiencies in oxidative phosphorylation have been reported in individuals with schizophrenia (SZ), bipolar disorder (BD), and major depressive disorder (MDD) in transcriptomic, proteomic, and metabolomic studies. Several mutations in mitochondrial DNA (mtDNA) sequence have been reported in SZ and BD patients.

Methodology/Principal Findings

Dorsolateral prefrontal cortex (DLPFC) from a cohort of 77 SZ, BD, and MDD subjects and age-matched controls (C) was studied for mtDNA sequence variations and heteroplasmy levels using Affymetrix mtDNA resequencing arrays. Heteroplasmy levels by microarray were compared to levels obtained with SNaPshot and allele specific real-time PCR. This study examined the association between brain pH and mtDNA alleles. The microarray resequencing of mtDNA was 100% concordant with conventional sequencing results for 103 mtDNA variants. The rate of synonymous base pair substitutions in the coding regions of the mtDNA genome was 22% higher (p = 0.0017) in DLPFC of individuals with SZ compared to controls. The association of brain pH and super haplogroup (U, K, UK) was significant (p = 0.004) and independent of postmortem interval time.

Conclusions

Focusing on haplogroup and individual susceptibility factors in psychiatric disorders by considering mtDNA variants may lead to innovative treatments to improve mitochondrial health and brain function.

About the Authors:

Brandi Rollins

Affiliation: Department of Psychiatry & Human Behavior, University of California Irvine, Irvine, California, United States of America

Maureen V. Martin

Affiliation: Department of Psychiatry & Human Behavior, University of California Irvine, Irvine, California, United States of America

P. Adolfo Sequeira

Affiliation: Department of Psychiatry & Human Behavior, University of California Irvine, Irvine, California, United States of America

Emily A. Moon

Affiliation: Department of Psychiatry & Human Behavior, University of California Irvine, Irvine, California, United States of America

Ling Z. Morgan

Affiliation: Department of Psychiatry & Human Behavior, University of California Irvine, Irvine, California, United States of America

Stanley J. Watson

Affiliation: Molecular and Behavioral Neuroscience Institute, University of Michigan, Ann Arbor, Michigan, United States of America

Alan Schatzberg

Affiliation: Department of Psychiatry, Stanford University, Palo Alto, California, United States of America

Huda Akil

Affiliation: Molecular and Behavioral Neuroscience Institute, University of Michigan, Ann Arbor, Michigan, United States of America

Richard M. Myers

Affiliation: Hudson Alpha Institute for Biotechnology, Huntsville, Alabama, United States of America

Edward G. Jones

Affiliation: Neuroscience Center, University of California Davis, Davis, California, United States of America

Douglas C. Wallace

Affiliation: Molecular and Mitochondrial Medicine and Genetics, University of California Irvine, Irvine, California, United States of America

William E. Bunney

Affiliation: Department of Psychiatry & Human Behavior, University of California Irvine, Irvine, California, United States of America

Marquis P. Vawter

* E-mail: [email protected]

Affiliation: Department of Psychiatry & Human Behavior, University of California Irvine, Irvine, California, United States of America

Introduction

Mitochondrial DNA (mtDNA) mutations can alter functions involving apoptosis, calcium buffering and signaling, and the brain's capacity to meet critical energy oxidative phosphorylation (OXPHOS) demands. The inheritance of mtDNA is normally matrilineal [1]. High rates of maternal offspring of SZ [2], [3], [4], [5] and BD [6] compared to paternal rates have been primarily found in family studies [7], [8] but not in the general population [9], [10], [11]. This supports the hypothesis that increased risk for these disorders might be related to mitochondrial dysfunction. Failure to find a pattern of maternal inheritance as well as evidence of paternal inheritance have also been reported in SZ and BD [12], [13].

Tissues with high energy requirements seem to accumulate mtDNA mutations more readily; thus, in general, the accumulation of mtDNA mutations is proportional to their metabolic rate [14]. Multiple mechanisms across the life span exist for transmission of mtDNA mutations [15], [16]. A cell may become heteroplasmic (increase of intra individual sequence variability in mtDNA) through a germ line or an acquired somatic mutation; thus, heteroplasmy may be inherited or acquired [17], [18]. MtDNA has been suggested as a potential “weak point” of the genome [19] due to susceptibility to somatic deletions as mtDNA is not protected by histones and lacks the repair machinery associated with nuclear DNA [20] which is one reason there is a higher rate of mutation or base pair substitutions in mtDNA compared to nuclear DNA. Heteroplasmic events in brain appear to be common in several studies [21], [22], [23], [24], [25]. Small changes in the level of heteroplasmy may be clinically relevant, and can appear in vulnerable cells and tissues in neurodegenerative disorders such as Parkinson's, Huntington's, and Alzheimer's Disease [16], [26], [27], suggesting a critical threshold of heteroplasmic mutation above which functional deficits and neuropathological damage occurs.

Evidence for dysfunction of mitochondria in psychiatric disorders

Psychiatric symptoms have been documented in subjects with mitochondrial disease [28], [29], [30]. Conversely, increased psychiatric symptoms have been associated with declines in mitochondrial functional activity [31], [32], [33], [34]. The neuroimaging literature suggests that energy metabolism abnormalities are widespread in the brains of subjects with SZ [35], [36], [37], [38], BD [39], [40], [41], [42], and MDD [43].

Common, novel, and ethnic specific defining haplogroup mitochondrial SNPs are associated with psychiatric disorders

Sequence variations in the mitochondrial genome are implicated in the pathogenesis of SZ and BD [44], [45], [46], [47], [48], [49], [50], [51], [52], [53], [54], [55], [56], [57], [58]. ND1, ND3, ND4, and ND5 SNPs have been associated with SZ or BD [51], [52], [59], [60]. An interesting SNP, the ND3 A10398G polymorphism, may be related to lithium response in BD [61] as well as increased pH in cybrid cells [62].

Ethnic specific defining haplogroup SNPs are also reported to be associated with psychiatric disorders. The sub-haplogroup F1ba of the haplogroup N, may harbor a haplotype risk factor for atypical psychosis [63] and in an Israeli sample of SZ there was an over-representation of the HV mtDNA haplogroup [64]. Several attempts have been made to assess association of certain mtDNA genetic variants with BD [44] in haplogroup T. A recent report indicated that the number of synonymous substitutions was increased in mtDNA obtained from blood samples in SZ [58], independent of haplogroup.

Based upon the studies implicating mitochondrial dysfunction in psychiatric disorders, the present study measured mtDNA variants in brain tissue, including substitutions, synonymous and non-synonymous, and rare variants, which could predispose to BD, MDD, and SZ. Secondly, we compared heteroplasmy detection by microarray using a published method [65] to conventional methods. Thirdly, this study examined the association between brain pH and mtDNA alleles. The hypothesis that mtDNA variants might influence brain pH is supported by the haplogroup defining ND3 A10398G polymorphism associated with an increased pH in cybrid cells [62]. Further, microarray studies have demonstrated differential nuclear-encoded mitochondrial gene expression [45], [66], [67] as well as mitochondrial gene expression [67] related to mean group pH differences. 31P magnetic resonance spectroscopy findings have also indicated a decrease in intracellular pH in BD [68] increasing the potential relevance of such an effect on psychiatric disorders.

Results

Microarray resequencing arrays and capillary electrophoresis sequencing of mtDNA

The Affymetrix Genetic Sequence Analysis Software (GSEQ) assigned high quality base pair sequencing calls (T,C,G,A) to 98.2% of the mtDNA genome on all 77 chips analyzed. There was a 100% concordance between microarrays resequencing the same individual repeated de novo. The call rate did not depend on the source of starting genomic DNA, as call rates from 8 freshly prepared lymphocyte samples were also 98.2% on average in the investigator's laboratory. Three individual mtDNA genomes were previously sequenced on an ABI 3130 capillary electrophoresis (CE) sequencer. All 103 calls that were different from the rCRS were 100% concordant between microarray and CE sequencing in these three individuals. The present mitochondrial microarray resequencing results are consistent with additional experimental verification from other investigators showing at least 99.99% microarray accuracy compared to sequencing [65], [69], [70], [71], [72], [73], [74], [75]. In addition, the microarray procedure detected the known heteroplasmic mutation A5793G with 45% heteroplasmy [65]. The same subject was found to be 30% heteroplasmic by PCR-RFLP analysis. Taken together, these initial results confirmed the mitochondria resequencing array is a useful tool for resequencing.

Haplogroup differences in psychiatric disorders

Upon completion of the resequencing arrays, the 77 mtDNA sequences from DLPFC were aligned with the rCRS. A phylogenetic tree (Fig. 1) was reconstructed with the MRBAYES program, which is a variant of the Markov chain Monte Carlo method [76], [77].

[Figure omitted. See PDF.]

Figure 1. Phylogenetic tree of human mtDNA coding sequence variants demonstrating regional association of haplogroups using mtDNA sequence from DLPFC.

Haplogroups, groups of related haplotypes, are derived from a founding haplotype, harboring characteristic mtDNA sequence polymorphisms. Each haplogroup is designated by a letter in this tree, the RCRS is the reference sequence; the ticks around the perimeter of the circle represent the individual mtDNA sequences from DLPFC. The faint numbers represent the probability of support for the branches which was 0.9–1.The mtDNA haplogroups have proven to be highly geographically associated [109]. The haplogroup assignments were calculated at MITOMASTER (http://mammag.web.uci.edu/twiki/bin/view/Mitomaster/WebHome) using all mtDNA sequence base pairs.

https://doi.org/10.1371/journal.pone.0004913.g001

The rCRS clustered far from the ancestral L haplogroup, and was correctly placed within the H haplogroup with an H2 subject as the nearest neighbor (Fig. 1). Individuals within each haplogroup are joined neighbors except for the node labeled M, N1B1, pre HV, which placed these three haplogroups as closest neighbors. One reason these haplogroups did not cluster separately as expected was that each of these samples also carried a high rate of synonymous substitutions (shown in the next section). The majority of 17 individuals (88%–100%) in this cluster carried mtSNPs different from rCRS at the following positions: 9540, 10652, 12705, 14783, 15043, 16311, 16325. The underlined mtSNPs were found in the M, N, or D haplogroups (/http://www.mitomap.org/mitomap-phylogeny.pdf). The 10652 mtSNP was found in one L2b individual (Mitokor 568 sequence [78]); 16325 and 16311 were not found in a specific haplogroup.

Previous haplogroup findings that were reported for BD [44] indicated a slight over-representation in the T haplogroup, thus the T haplogroup was explored further in the current dataset. The opposite finding emerged, controls were significantly over-represented in haplogroup T compared to BP+SZ (p = 0.04 permuted p-value). There was complete concordance of findings for all nine SNPs defining this haplogroup (709, 1888, 4917, 8697, 10463, 13368, 14905, 15607, 15928). This concordance illustrates confidence for haplogroup assignment using microarray resequencing.

Another exploratory difference in the distribution of psychiatric groups by haplogroup (Table 1) was seen in analysis of the combined BD+SZ distribution of subjects found more frequently in the pre HV group (Table 2) compared to all haplogroups (Chi square = 13.7, df (6), p = 0.032). To further test for association of mtDNA SNPs we analyzed the SZ, BD, MDD mtDNA SNPs dataset (484 polymorphisms) using PLINK [79].

[Figure omitted. See PDF.]

Table 1. The haplogroup determination of mtDNA derived from DLPFC showed 31 diverse haplogroups that could be reconstructed.

https://doi.org/10.1371/journal.pone.0004913.t001

[Figure omitted. See PDF.]

Table 2. The haplogroup determination of mtDNA derived from DLPFC was condensed from Table 1 and showed that BD and SZ were over-represented in Pre HV haplogroup and MDD was over-represented in the “other” haplogroups.

https://doi.org/10.1371/journal.pone.0004913.t002

Case-control differences in polymorphisms

An allelic association analysis was conducted across polymorphic sites in mtDNA sequence from the 77 subjects. Permutation analysis (n = 50,000) for each polymorphism was run using PLINK [79] to establish an empirical p-value for SZ v C, MDD v C, and BD v C comparisons. There were 484 nominal polymorphic sites, and in this first pass analysis, all SNPs were calculated for association, without regard to haplogroup specificity or minor allele frequency. Only empirical p-values<0.05 are reported for case-control comparisons (Table 3). There were two 12S rRNA rare polymorphisms found in SZ and not in C. The 12S mtDNA G750A and G1438A SNPs had frequencies of 0.143 and 0.214 in SZ, respectively. Both substitutions have been reported in 2,704 mtDNA sequences at a rate of 0.008 for G750A and 0.032 for G1438A [80]. The empirical p-values for association with G750A and G1438A were 0.037 and 0.012, respectively. The ND2 gene G4769A synonymous substitution showed a significant allelic association with SZ (empirical p value 0.037).

[Figure omitted. See PDF.]

Table 3. Mitochondrial polymorphisms were analyzed in PLINK for allelic association with 50,000 case-control permutations to obtain empirical p-values.

https://doi.org/10.1371/journal.pone.0004913.t003

For BD, there were three allelic associations significant in the D-loop control region: T114C (p = 0.032), C195T (p = 0.008), and G16300A (p = 0.012). The T114C and G16300A minor alleles are rare at 0.004 and 0.003 in a database of 2,704 mtDNA genomes [80]. Three associations for MDD are shown in Table 3.

Microarray Heteroplasmy

The present microarray data shows the ND4L T10652C variant appears within the N (N1B1, pre HV) and M (D1) haplogroups. There were 17 subjects within these haplogroups that had a heteroplasmy level greater than 50% by microarray calculation as previously described [65]. This ND4L T10652C transition is rare and was not present in a human mtDNA database of 2,704 sequences distributed among major haplogroups [80]. In addition it was not observed in 12 mtDNA samples from control blood measured by microarray resequencing, reducing the likelihood of a microarray artifact. There was a common haplotype for the putative T10652C individuals at four mtSNPs – T10652C, T14783C, G15043A, and T9540C for 15 of 16 subjects. These three mtSNPs appear early in the phylogenetic mtDNA evolutionary tree at L2b, and the roots of M, and N. A post hoc analysis of each SNP was run separately and pooled for the three psychiatric groups (SZ, BD and MDD) and compared to controls. A trend for an association for each SNP was found (T14783C p = 0.06, G15043A p = 0.07, T9540C p = 0.16). These observations prompted further investigation of the variant by sequencing, SNaPshot, and allele specific PCR. The ND4L T10652C variant was confirmed by direct sequencing of mtDNA for the highest heteroplasmic individual (81% C heteroplasmy by microarray calculation) within the D1 haplogroup who was an MDD subject and had the C allele (Fig. 2). The other 16 individuals within the N (N1B1, pre HV) and M (D1) haplogroups showed the T allele with direct sequencing. The heteroplasmy detection limit of conventional direct CE sequencing is 20%–30% [81], [82], indicating that the microarray results were inflated for these subjects. All subjects were resequenced 100 bp on either side of the ND4L T10652C for insertions or deletions. There were no insertions or deletions found suggesting that the microarray results were not due to probe mismatches. Primers used for sequencing and allelic specific PCR did not amplify genomic DNA extracted from ρ° cells depleted of mtDNA [83].

[Figure omitted. See PDF.]

Figure 2. Electropherogram results for direct sequencing of the ND4L T10652C variant.

(A) is the ‘T’ allele sequencing result from the lowest heteroplasmic individual used to make the ‘T’ allele clone standard for SNaPshot and allele specific PCR. (B) is the ‘C’ allele sequencing result from the highest heteroplasmic individual detected with microarray (MDD subject from haplogroup D1) used to make the ‘C’ allele clone standard.

https://doi.org/10.1371/journal.pone.0004913.g002

SNaPshot multiplex was tried for validation of the microarray calls at T10652C. Twenty four subjects with varying heteroplasmy levels by microarray were run and all individuals showed the T allele indicating their heteroplasmy levels, if present, were below the detection limit of this method. Although we ran a standard curve consisting of dilutions from 100% down to 0.1% for the C allele, we were only reliably able to detect 5% or greater heteroplasmy using SNaPshot.

A more sensitive method for validation of the calculated heteroplasmy from the microarray resequencing data was performed with an allelic specific real-time detection and quantification experiment for the T10652C variant (Fig. 3). The subject with the highest heteroplasmy in the initial microarray experiment (81%) showed similar levels in the confirmation experiment (88%). Because the allelic specific experiment can detect low levels of heteroplasmy (0.25%), we were able to detect several other subjects with heteroplasmy >0.25%. Two SZ subjects (from haplogroups pre HV and M) with 69% and 58% heteroplasmy in the microarray experiment had 0.47% heteroplasmy in the real-time allelic specific qPCR experiment. Despite the failed confirmation of some of the microarray results with 50% and higher heteroplasmy, the three subjects with somewhat higher levels of heteroplasmy were all psychiatric subjects.

[Figure omitted. See PDF.]

Figure 3. Scatter plot of the heteroplasmy levels of ND4L T10652C calculated from the Affymetrix mitochondrial resequencing chip by method of Coon et al. [65] on the y-axis and allelic specific PCR (x-axis truncated).

https://doi.org/10.1371/journal.pone.0004913.g003

Substitution bias in subjects with psychiatric disorders and controls

One measure of neutrality for mtDNA variants was calculated by comparing the number of non synonymous substitutions (dN) to the number of synonymous substitutions (dS) in each group. Although it is generally accepted that S changes are greater than N changes in mtDNA, we tested for differences between groups in this rate using pair wise comparisons. Values of p<0.05 were considered significant and tabulated (Table 4, column [A] shows the number of significant pair wise comparisons). Table 4 shows the number of times the null hypothesis of strict-neutrality (dN = dS) was accepted in [B] or rejected in favor of the alternative hypothesis (dN<dS, column [A]) using pair wise comparisons for all possible subject combinations of BD-C, SZ-C, MDD-C, and C-C. Occurrences that were significant are counted in column [A]. The higher rate of synonymous substitutions in schizophrenia by 22% compared to the controls was significant (p = 0.0017).

[Figure omitted. See PDF.]

Table 4. Codon-based test of purifying selection for pair wise analysis between sequences.

https://doi.org/10.1371/journal.pone.0004913.t004

A second analysis of the overall transition/transversion bias (R) for all DLPFC mtDNA sequences was conducted using the Maximum Composite Likelihood method in MEGA4 [84], [85] for each of the four groups separately. Codon positions included were 1st, 2nd, 3rd. The observed transition/transversion ratios were 2.25 fold higher in SZ compared to C (Table 5). This suggests a higher number of substitutions in DLPFC in SZ compared to controls in the 13 mtDNA transcripts. Two caveats apply to this difference in bias: the overall sample size is small (n = 77) and the total number of transversions in the entire sample was modest (n = 22) which can inflate R. Further, we used a homogeneous model to assume that all sites are variable, and may over estimate the number of substitutions based upon pairs of sequences. With these caveats, Tables 4 and 5 consistently showed the number of substitutions in SZ pair wise comparisons exceeded the rate in control comparisons and the overall R rate is higher in SZ than in C. The present R ratio (transition/transversion) data was examined in the coding region functional grouping level by complex, the rate for R in mitochondrial complex I was 10 fold higher in SZ compared to controls.

[Figure omitted. See PDF.]

Table 5. Transition/transversion bias calculation for mtDNA coding sequences shows difference between SZ and control group substitution.

https://doi.org/10.1371/journal.pone.0004913.t005

Association with postmortem pH

Postmortem pH was next analyzed for association with mtDNA SNPs. Association and permutation analysis were conducted in PLINK and the permuted p-values that were most significant for the association of pH were three mtDNA SNPs in tRNA Leu, ND4, and ND5 (association p-values 0.002–0.0004, Table 6). These same three mtDNA SNPs are ethnic-specific haplotype-defining polymorphisms for the super haplogroup cluster (U, K, UK) and all subjects with these three mtDNA SNPs were assigned to these haplogroups in this study using MITOMASTER (www.mitomaster.org). This super haplogroup (U, K, UK) had a significantly higher postmortem pH (7.00±0.18 SD) compared to the remaining subjects pH (6.8±0.18 SD) and the permuted p-value for the pH difference was 0.01 for 5,000 random group assignments for membership in the super haplogroup cluster (U, K, UK).

[Figure omitted. See PDF.]

Table 6. pH in postmortem brain showed significant association with three mtDNA SNPs.

https://doi.org/10.1371/journal.pone.0004913.t006

The present study design intentionally avoided some peri-mortem artifacts by pre-selecting subjects with rapid deaths and without recorded agonal factors in the autopsy database. The pH association did not appear spurious due to differences in assignment of cases and controls as the distribution was equally distributed for psychiatric cases and controls between U, K, UK super haplogroup and the remaining haplogroups. The pH shift towards the right for the U, K, UK super haplogroup and all other subjects (Fig. 4) showed that the lowest nine brain pH readings were all from subjects in the other super haplogroups besides U, K, or UK. Removing the lowest five pH subjects did not change the association of pH nor reduce the bootstrap significance of the difference between haplogroups. However, to exclude if the effect could be primarily related to factors like postmortem interval or agonal factors, we controlled for both effects. The PMI was not different between the U, K, UK super haplogroup and the remaining haplogroups, and including PMI in ANCOVA did not reduce the significance of the pH differences (p = 0.0007 with PMI as a covariate).

[Figure omitted. See PDF.]

Figure 4. The super haplogroup (U, K, UK) showed a shift in postmortem brain pH, and this finding was significant following permutation analysis.

The histogram shows the frequency of the pH (bars) for the super haplogroup (red) compared to all other matrilineages (blue). The accumulated percentage (right y-axis) is shown as two lines with the super haplogroup (red line) and all other matrilineages (blue line). There were no subjects with a prolonged death, or agonal factors as rated according to the Hardy et al. scale [98]. PMI differences did not account for this significant effect, as PMI was equivalent between the super haplogroup and other haplogroups.

https://doi.org/10.1371/journal.pone.0004913.g004

Discussion

The present study analyzed mtDNA sequence data collected from 17 major yet diverse haplogroups residing in Southern California. Over 1.6×106 bp of high quality mtDNA genetic data was generated. Mitochondrial array resequencing data was 99.999% identical to conventional CE sequencing [69] indicating that the mitochondrial resequencing array is a reliable method. The mitochondrial microarray, however, showed a large N call rate, which may be improved by using similar PCR conditions as the National Institute of Standards and Technology. Additionally, changing the base pair calling algorithm [69] could increase the average call rate to >99% across the entire mtDNA genome [69], [70], [86]. In addition, the mitochondrial microarray showed a high heteroplasmy level using the calculation method of Coon et al. [65] for the ND4L T10652C SNP which could not be validated entirely with SNaPshot or allele specific PCR.

The main finding relevant to schizophrenia was an increased number of synonymous base substitutions in mtDNA in DLPFC. The first method used a codon purifying selection method [84] in coding regions only, and examined pair wise concordances of sequences to establish, for each group, whether there were synonymous substitutions in codons. The rate of synonymous substitutions in schizophrenia was 22% higher compared to the control rate (p = 0.0017).

The second calculation method was a transversion/transition rate bias (R) [84] performed for each group separately. R was increased 2.25 fold in SZ compared to controls. The present R ratio (transition/transversion) data was examined on individual mitochondrial complex functional grouping levels, the rate for R in mitochondrial complex I was 10 fold higher in SZ compared to controls. This suggests an increase in the number of substitutions in DLPFC in SZ compared to controls preferentially in mtDNA transcripts in complex I. There is a previous report of age related increases of synonymous substitutions in brain [87]. It has been previously demonstrated that there is an increase in the rate of substitutions in somatic mtDNA in brain[21]. Increased synonymous substitutions have also been reported in blood mtDNA in SZ [58]. A possible mechanism could be a result of nuclear gene defects in genes such as POLG. POLG is a proof reading enzyme for mtDNA which was shown in association to Parkinson's disease [88], [89]. POLG knock down mice have increased mtDNA synonymous substitution rates [90].

Heteroplasmy results by microarray at ND4L 10652 were largely not concordant with two separate validation attempts. One MDD subject was shown to have 81% heteroplasmy with microarray, 88% heteroplasmy with allele specific PCR, the C allele with direct sequencing and was not run with SNaPshot. Only two additional subjects were confirmed to have ‘significant’ heteroplasmy levels for this rare variant, 0.47% compared to the microarray resequencing results which were 69% and 58%. These three subjects had psychiatric disorders and were within the M, N and pre HV haplogroups. Further work is required to determine why the microarray T10652C mtDNA variants could not be validated, since the putative carriers had a four-fold higher relative risk of psychiatric disorders compared to controls, and died at younger ages both in control and psychiatric comparisons.

In a study of 18 Alzheimer cases and 19 controls there was a statistically significant difference reported for the T10652 expected reference signal (p<0.05), although it was not specified if it was the C allele substitution [65]. Given the fact that this SNP was not present in 2,704 sequences distributed among major haplogroups and only reported in one L2b individual (Mitokor 568 sequence), the high point prevalence in both our study and the Alzheimer's study is suspect since both studies used Affymetrix mitochondrial resequencing arrays.

With the above caveat in mind, we nevertheless found a common haplotype of four SNPs – T10652C, T14783C, G15043A, and T9540C in 15 of 16 subjects in the M, N and pre HV haplogroups. Individuals with this haplotype cluster separately from the remainder of subjects near the L1 and L2 haplogroups suggesting an ancient origin. It is tempting to speculate that the common haplotype might be a marker for the increased substitutions in brain mtDNA which is associated with younger ages of morbidity in both psychiatric disorder and control groups. The haplotype (T14783C, G15043A, and T9540C) could be further investigated for functional mitochondrial differences. A trend for an association for each SNP was found (T14783C, G15043A, T9540C) suggesting that although these SNPs did not pass permutation corrected p values<0.05, there is a putative haplotype association with risk for psychiatric disorders that can be tested in larger sample sizes.

There was a recent finding that schizophrenia was over-represented in the HV haplogroup with increased relative risk (odds ratio) of 1.8 [64]. Our findings for haplogroup association with psychiatric disorders were significant for pre HV showing increased risk, T showing a protective factor. While based upon a small sample size these findings clearly warrant replication in a larger sample size.

Brain pH showed an association with super haplogroup (U, K, UK) which had a significantly higher pH, i.e. a 1.54 fold decrease in H+ concentration compared to the remaining haplogroups. If haplogroup specific variations in mtDNA can produce alterations of pH levels in human brain, then it could be anticipated that disease associations might also be seen in a larger study. In fact, there was a protective effect of the haplogroup UJKT [91] in Parkinson's disease, which taken together with evidence of higher brain pH in the U, UK, and K subjects might be an indicator of more adaptive brain physiological responses to hypoxia. It has been hypothesized that certain haplogroups are more tightly coupled in the following order from most to least (H>J = T>U>UK) [92]. The haplogroups that are more tightly coupled will generate more reactive oxygen and the less coupled U and UK would have looser coupling efficiency and generate the least ROS in this group. This loose coupling could explain the higher pH, due to less excess mitochondrial oxidation and decreased H+ ion gradients in the outer membrane. A mitochondrial component in diathesis to schizophrenia and mood disorders might involve more tightly coupled OXPHOS processes that occur in H haplogroups on the left of the inequality expression. This tendency towards association of psychiatric disorders with H haplogroups was observed in this study and other studies [58], [64].

Synonymous base pair changes were observed in mtDNA genes in DLPFC: ND2, ND4L, ND6, and Cytb. Three of these ‘associated’ genes ND2, ND4L, and ND6, are located in complex I of the mitochondrion. Associations with complex I mtDNA SNPs such as ND4 T12027C and ND5 C12403T have been reported in schizophrenia [59], [60] and ND1 T3644C mutation in complex I in BD [52]. A decreased expression of complex I mitochondrial transcripts in SZ was recently shown in DLFPC [93]. In addition, a decreased expression of nuclear encoded subunits of complex I was shown in SZ in prefrontal cortex and striatum [94]. This is consistent with inhibition of complex I by various antipsychotics. Haloperidol potently inhibits complex I, followed by chlorpromazine, fluphenazine and risperidone while the atypical neuroleptic, clozapine, did not inhibit complex I activity in mouse brain slices [95]. Neuroleptics such as haloperidol and flupenthixol reduce complex I activity in brain of subjects with SZ [96]. The production of superoxide anion O2−, the principle agent of oxidative stress, is reduced in brain homogenates from rats treated with flupenthixol and reduced in brain tissue and lymphocytes from schizophrenics receiving neuroleptic medication [96]. These studies support the view that a classical neuroleptic such as haloperidol inhibits mitochondrial complex I and reactive oxygen species (ROS) [97]. Complex I deactivation pathophysiology is observed in postmortem brain at both the protein and mRNA levels and the question of whether medications and interactions with mtSNPs requires additional studies. It might be that mtSNPs modulate the responses to antipsychotic treatment at the mitochondrial level as suggested by a recent mtDNA SNP pharmacogenomic study [61].

In summary, this study analyzed mitochondrial DNA variation in brain and found increased mitochondrial DNA substitutions in subjects with SZ compared to controls similar to findings in blood [58] for SZ. Haplogroup associations with risk and protective factors for SZ and BP were found that require further replication. One MDD subject carried essentially a homoplasmic mutation in DLPFC at ND4L T10652C, and two subjects with SZ showed less than 1% heteroplasmy; however, the pathological significance of the low levels of this non synonymous heteroplasmic mutation is unknown. This study replicated that low levels of heteroplasmic non-synonymous mutations are found in the brain [87]. A third finding showed that the UK macro haplogroup maintained a higher pH in response to terminal brain hypoxia. The present findings taken together could be useful in mitochondrial medicine for earlier diagnosis and treatment of brain related disorders based upon recognition of haplogroup and individual susceptibility factors.

Methods

Ethics Statement

Informed consent to allow genetic testing was obtained from the next of kin for each sample and the collected samples were approved by the University of California, Irvine Institutional Review Board.

Synonymous and non synonymous substitutions in mtDNA in all coding and control regions were identified in dorsolateral prefrontal cortex (DLPFC) by microarray resequencing of schizophrenia (SZ), bipolar disorder (BD), major depressive disorder (MDD), and control (C) subjects. The demographics for the study are shown for the mitochondrial DNA resequencing analysis (Table 7). A total of 77 subjects (12 with BD, 14 with SZ, 15 with MDD, and 36 controls) were initially studied. The majority of subjects were identified as Caucasian (71), the remaining were African - American (1), Asian (2), Other (2), and Pacific Islander (1). Age, gender, pH, PMI and DNA crossing threshold (Ct) of each subject are summarized in Table 7. The average pH for MDD subjects was slightly higher than all other groups and the mean age for SZ was lower than all other groups. We included subjects with a rapid death, and no known agonal factors related to chronic illness. Rapidity of death was scored according to the Hardy et al. scale [98]. A rapid death occurs within 1 hr from agonal onset based upon coroner's office estimates and is rated on the Hardy scale as 1 or 2. For qualitative purposes, examples of agonal factors are medical conditions at the time of death such as coma, pyrexia, skull fracture, hypoxia, dehydration, multi-organ failure, seizures and ingestion of neurotoxic substances [98], [99], [100], [101]. Each case and control was screened by an extensive review of the medical examiner's conclusions, coroner's investigation, medical records, toxicology results, and interviews of the decedent's next-of-kin.

[Figure omitted. See PDF.]

Table 7. Demographics of subjects from which DLPFC mtDNA was resequenced.

https://doi.org/10.1371/journal.pone.0004913.t007

Suicidal death occurred in 17 SZ, MDD, and BD cases. Brain samples from SZ, BD, MDD, and controls were collected previously from the Orange County Medical Examiner's Office [102]. Each brain underwent a gross neuropathological examination to rule out lesions, infarcts, and gross abnormalities [103]. The neuropathology written opinion was based upon careful examination of high resolution digital images of 28 coronal slices per whole brain. The 1 cm coronal slabs were frozen at −140°C, placed into sealed bags, and transferred to −80°C storage. A best estimate clinical diagnosis was based upon DSM-IV criteria applied to medical and hospital records, and next of kin interviews conducted by two psychologists, and were reviewed by a senior psychiatrist.

The dorsolateral prefrontal cortex (DLPFC Brodmann areas 9/46) was dissected from coronal slabs on the left side using neuroanatomical landmarks rostral to the corpus callosum from the medial frontal gyrus on a platform cooled with dry ice. The dissected prefrontal regions from each brain were weighed frozen (∼100 mg) and processed. DNA was extracted from brain tissue using a standard Trizol protocol. Briefly, the DNA was extracted from the phenol phase and precipitated with ethanol. The DNA was then washed twice with 0.1 M sodium citrate in 10% ethanol and solubilized in 8 mM NaOH. The extraction technique provided excellent recovery and quality of nucleic acid assessed by qPCR (Table 7).

Affymetrix Mitochondrial Resequencing Arrays

Genomic DNA obtained from DLPFC was used on Affymetrix GeneChip Mitochondrial Resequencing Arrays 2.0 following manufacturer's protocol (Affymetrix, Santa Clara, CA). Briefly, the mitochondrial resequencing chip has both strands of the entire human mitochondrial coding sequence (15,451 bp) arrayed; both strands of an additional 12,935 bp (84% of coding DNA) are arrayed in duplicate. We used 300 ng of genomic DNA to amplify the mitochondrial genome in three overlapping long PCR fragments and pooled equimolar amounts of amplicon using the protocol from the Affymetrix GeneChip Resequencing Reagent Kit. The reliability of mtDNA resequencing was determined by comparing mtDNA sequence on the same subjects measured twice by the resequencing microarray. The accuracy of the resequencing microarray results was compared to conventional ABI direct sequencing for three subjects.

Statistical Analysis

Frequencies of assigned allele calls were compared by using the allele chi-square test in PLINK for a case - control analysis of MT chromosome, and running permutation analysis (n = 50,000 permutations) for determining empirical p-values (URL: http://pngu.mgh.harvard.edu/purcell/plink/) [79]. Associations between mtDNA ethnic specific haplogroup variants and postmortem pH values were tested in PLINK.

Tests of substitution rate differences between groups were conducted as follows. The pair wise comparisons of substitutions between SZ, BD, MDD, and controls were calculated with Molecular Evolutionary Genetics Analysis (MEGA software version 4.0) [84] using the Nei-Gojobori method [84], [85], [103]. The probability of rejecting the null hypothesis of strict-neutrality of selection (dN = dS) in favor of the alternative hypothesis (dN<dS) was calculated for the test statistic (dS−dN). The variance of the difference was computed using the bootstrap method (10,000 replicates). Values of p<0.05 were considered significant and tabulated for each group and group frequencies were tested with Pearson's chi-square. An estimation of the transition/transversion bias for each group was also calculated in MEGA4 [84].

Haplogroup assignment and mtDNA transversions/transitions were determined with the MitoMaster database [104], [105] in the Center for Molecular and Mitochondrial Medicine and Genetics. This uses neighbor joining distances to assign haplogroups based upon specific mtDNA haplotype.

Heteroplasmy was screened by microarray calculation with the method of Coon et al. [65]. The microarray signal intensity data from all four alleles (A, C, G, and T), on both the sense and antisense strands formed the relative expected allele intensity. The expected allele was called using the revised Cambridge Reference Sequence (rCRS) [106] (Accession # REFSEQ AC_000021.2 gi∶115315570).

Direct Sequencing and Cloning of ND4L T10652C Variant

To validate the ND4L T10652 variant, PCR was performed on the 17 DLPFC samples that showed the highest levels of heteroplasmy at this locus with the resequencing arrays. These subjects were in the N (N1B1, pre HV) and M (D1) haplogroups. Primers were designed using Primer3 software (Totowa, NJ) (forward 5′-TCGCTCACACCTCATATCCTC-3′ and reverse 5′-GGCCATATGTGTTGGAGATTG-3′). PCR reactions were run on a 2% agarose gel and fragments were cut out and purified using the QIAquick Gel Extraction Kit (Qiagen, Valencia, CA). Sequencing reactions were run using the forward PCR primer and the Big Dye Terminator Kit (Applied Biosystems, Foster City, CA). Products were purified using a DyeEx 2.0 Spin Kit (Qiagen) and analyzed on an ABI Prism 3100 capillary sequencer. Sequences were read in Geneious (Biomatters, Ltd., New Zealand) for alignment around mitochondrial locus 10652. The sample showing the highest heteroplasmy level (81%) by microarray (the MDD sample in haplogroup D1) was validated by capillary sequencing and used to create the T allele clone for the standard curve for SNaPshot analysis. A second sample, showing very low heteroplasmy by microarray, was chosen for the C allele clone for the standard curve. Both clones were used for PCR with the ND4L primers and PCR products were gel purified and cloned into the pCR®2.1-TOPO® vector using the TOPO-TA Cloning® Kit (Invitrogen Life Technologies, Carlsbad, CA). Transformations used One Shot® Mach1™-T1R competent cells and were plated out on LB plates with X-gal and kanamycin. White and light blue colonies were selected and minipreps were purified using the PureLink™ Quick Plasmid Miniprep Kit (Invitrogen) according to the manufacturer's instructions. Clones were tested for correct insert sequence and orientation using the M13 forward and reverse primers.

SNaPshot Standard Curve

Established clones for the T and C alleles at ND4L 10652 were used to create the standard curve for SNaPshot analysis. The ABI Prism SNaPshot Multiplex Kit (ABI, Foster City, CA) was used according to the manufacturer's protocol with 0.10 pmols of template. Seven dilutions were made for the standard curve: 100% C allele, 83% C allele, 10% C allele, 5% C allele, 1% C allele, .1% C allele and 100% T allele. The SNaPshot single-base extension reactions were run on the 3100 Genetic Analyzer (ABI, Foster City, CA). DNA fragment data was analyzed using the Peak Scanner Software v1.0 (ABI, Foster City, CA) to obtain the ratio of the peak heights corresponding to each allele for the different levels of heteroplasmy in the standard curve. Once a calibration curve was established with varying dilutions of each clone, 24 mtDNA samples from DLPFC that showed varying levels of heteroplasmy with microarray were run using the SNaPshot Multiplex Kit.

Allele Specific PCR

An allelic specific real-time detection and quantification of mtDNA variants was implemented to more precisely confirm the percentage heteroplasmy observed in the microarray experiment. Low levels of heteroplasmy down to 0.02% were detected using a method derived from Strand et al. [107]. This involved an allele specific amplification by PCR [108], followed by the quantification of the amounts of T and C alleles using the Light Cycler 480 (Roche, Indianapolis, IN). The standard curve was constructed from the 10652 clones (see above) using absolute copy number proportions of each allele (500 to 1 million copies). Copy numbers of each allele for the 77 DLPFC samples were then calculated from the standard curve after an allelic specific PCR to determine the unknown “Cp” of each allele. The heteroplasmy was calculated as the % C over the total (C+T). High specificity and sensitivity for the mutant C allele were achieved by blocking non-specific amplification of the wild type allele with an amine blocked T specific primer.

Acknowledgments

We appreciate the assistance of Ling Shao, Ph.D. for DNA extraction, and Preston Cartagena, Psy.D, David Walsh, Psy.D, and Richard Stein, Ph.D. for their contributions to postmortem collections and clinical characterization of subjects. We acknowledge Kathleen Burke as well as Jacque Berndt and the investigators and medical examiners at the Orange County Coroner's Office for procurement of brain tissue. F. Warren Lovell, M.D, performed a neuropathological evaluation of the postmortem brains. Tissue specimens were processed and stored at the University of California, Davis and the Human Brain and Spinal Fluid Resource Center, Veteran's Medical Center, Los Angeles under the direction of Wallace W. Tourtellotte, M.D., Ph.D.

Author Contributions

Conceived and designed the experiments: PAS MV. Performed the experiments: BR PAS EAM LZM. Analyzed the data: MVM PAS DW MV. Contributed reagents/materials/analysis tools: SJW AS HA RM EJ DW WEB. Wrote the paper: BR DW MV.

Citation: Rollins B, Martin MV, Sequeira PA, Moon EA, Morgan LZ, Watson SJ, et al. (2009) Mitochondrial Variants in Schizophrenia, Bipolar Disorder, and Major Depressive Disorder. PLoS ONE4(3): e4913. https://doi.org/10.1371/journal.pone.0004913

References

1. Giles RE, Blanc H, Cann HM, Wallace DC (1980) Maternal inheritance of human mitochondrial DNA. Proc Natl Acad Sci U S A 77: 6715–6719.RE GilesH. BlancHM CannDC Wallace1980Maternal inheritance of human mitochondrial DNA.Proc Natl Acad Sci U S A7767156719

2. Gottesman II, Bertelsen A (1991) Confirming unexpressed genotypes for schizophrenia. Risks in the offspring of Fischer's Danish identical and fraternal disocordant twins. Arch Gen Psychiatry 46: 867–872.II GottesmanA. Bertelsen1991Confirming unexpressed genotypes for schizophrenia. Risks in the offspring of Fischer's Danish identical and fraternal disocordant twins.Arch Gen Psychiatry46867872

3. Goldstein JM, Faraone SV, Chen WJ, Tsuang MT (1992) Gender and the familial risk for schizophrenia. Disentangling confounding factors. Schizophr Res 7: 135–140.JM GoldsteinSV FaraoneWJ ChenMT Tsuang1992Gender and the familial risk for schizophrenia. Disentangling confounding factors.Schizophr Res7135140

4. Wolyniec PS, Pulver AE, McGrath JA, Tam D (1992) Schizophrenia: gender and familial risk. J Psychiatr Res 26: 17–27.PS WolyniecAE PulverJA McGrathD. Tam1992Schizophrenia: gender and familial risk.J Psychiatr Res261727

5. Shimizu A, Kurachi M, Yamaguchi N, Torii H, Isaki K (1987) Morbidity risk of schizophrenia to parents and siblings of schizophrenic patients. Jpn J Psychiatry Neurol 41: 65–70.A. ShimizuM. KurachiN. YamaguchiH. ToriiK. Isaki1987Morbidity risk of schizophrenia to parents and siblings of schizophrenic patients.Jpn J Psychiatry Neurol416570

6. McMahon FJ, Stine OC, Meyers DA, Simpson SG, DePaulo JR (1995) Patterns of maternal transmission in bipolar affective disorder. Am J Hum Genet 56: 1277–1286.FJ McMahonOC StineDA MeyersSG SimpsonJR DePaulo1995Patterns of maternal transmission in bipolar affective disorder.Am J Hum Genet5612771286

7. Gershon ES, Badner JA, Detera-Wadleigh SD, Ferraro TN, Berrettini WH (1996) Maternal inheritance and chromosome 18 allele sharing in unilineal bipolar illness pedigrees. Am J Med Genet 67: 202–207.ES GershonJA BadnerSD Detera-WadleighTN FerraroWH Berrettini1996Maternal inheritance and chromosome 18 allele sharing in unilineal bipolar illness pedigrees.Am J Med Genet67202207

8. Lin JP, Bale SJ (1997) Parental transmission and D18S37 allele sharing in bipolar affective disorder. Genet Epidemiol 14: 665–668.JP LinSJ Bale1997Parental transmission and D18S37 allele sharing in bipolar affective disorder.Genet Epidemiol14665668

9. Kato T, Winokur G, Coryell W, Rice J, Endicott J, et al. (1998) Failure to demonstrate parent-of-origin effect in transmission of bipolar II disorder. J Affect Disord 50: 135–141.T. KatoG. WinokurW. CoryellJ. RiceJ. Endicott1998Failure to demonstrate parent-of-origin effect in transmission of bipolar II disorder.J Affect Disord50135141

10. Grigoroiu-Serbanescu M, Martinez M, Nothen MM, Grinberg M, Sima D, et al. (2001) Different familial transmission patterns in bipolar I disorder with onset before and after age 25. Am J Med Genet 105: 765–773.M. Grigoroiu-SerbanescuM. MartinezMM NothenM. GrinbergD. Sima2001Different familial transmission patterns in bipolar I disorder with onset before and after age 25.Am J Med Genet105765773

11. Grigoroiu-Serbanescu M, Martinez M, Nothen MM, Propping P, Milea S, et al. (1998) Patterns of parental transmission and familial aggregation models in bipolar affective disorder. Am J Med Genet 81: 397–404.M. Grigoroiu-SerbanescuM. MartinezMM NothenP. ProppingS. Milea1998Patterns of parental transmission and familial aggregation models in bipolar affective disorder.Am J Med Genet81397404

12. DeLisi LE, Razi K, Stewart J, Relja M, Shields G, et al. (2000) No evidence for a parent-of-origin effect detected in the pattern of inheritance of schizophrenia. Biol Psychiatry 48: 706–709.LE DeLisiK. RaziJ. StewartM. ReljaG. Shields2000No evidence for a parent-of-origin effect detected in the pattern of inheritance of schizophrenia.Biol Psychiatry48706709

13. Kornberg JR, Brown JL, Sadovnick AD, Remick RA, Keck PE Jr, et al. (2000) Evaluating the parent-of-origin effect in bipolar affective disorder. Is a more penetrant subtype transmitted paternally? J Affect Disord 59: 183–192.JR KornbergJL BrownAD SadovnickRA RemickPE Keck Jr2000Evaluating the parent-of-origin effect in bipolar affective disorder. Is a more penetrant subtype transmitted paternally?J Affect Disord59183192

14. Shoffner JM, Wallace D (1995) Oxidative phosphorylation diseases. In: Scriver CR, Beaudet AL, Sly WS, Valle D, editors. The Metabolic and Molecular Bases of Inherited Disease, 7th edition, vol. I. New York: McGraw-Hill. pp. 1535–1609.JM ShoffnerD. Wallace1995Oxidative phosphorylation diseases.CR ScriverAL BeaudetWS SlyD. ValleThe Metabolic and Molecular Bases of Inherited Disease, 7th edition, vol. INew YorkMcGraw-Hill15351609

15. Wallace DC (2005) A mitochondrial paradigm of metabolic and degenerative diseases, aging, and cancer: a dawn for evolutionary medicine. Annu Rev Genet 39: 359–407.DC Wallace2005A mitochondrial paradigm of metabolic and degenerative diseases, aging, and cancer: a dawn for evolutionary medicine.Annu Rev Genet39359407

16. Horton TM, Graham BH, Corral-Debrinski M, Shoffner JM, Kaufman AE, et al. (1995) Marked increase in mitochondrial DNA deletion levels in the cerebral cortex of Huntington's disease patients. Neurology 45: 1879–1883.TM HortonBH GrahamM. Corral-DebrinskiJM ShoffnerAE Kaufman1995Marked increase in mitochondrial DNA deletion levels in the cerebral cortex of Huntington's disease patients.Neurology4518791883

17. Fan H, Civalier C, Booker JK, Gulley ML, Prior TW, et al. (2006) Detection of common disease-causing mutations in mitochondrial DNA (mitochondrial encephalomyopathy, lactic acidosis with stroke-like episodes MTTL1 3243 A>G and myoclonic epilepsy associated with ragged-red fibers MTTK 8344A>G) by real-time polymerase chain reaction. J Mol Diagn 8: 277–281.H. FanC. CivalierJK BookerML GulleyTW Prior2006Detection of common disease-causing mutations in mitochondrial DNA (mitochondrial encephalomyopathy, lactic acidosis with stroke-like episodes MTTL1 3243 A>G and myoclonic epilepsy associated with ragged-red fibers MTTK 8344A>G) by real-time polymerase chain reaction.J Mol Diagn8277281

18. Fan W, Waymire KG, Narula N, Li P, Rocher C, et al. (2008) A mouse model of mitochondrial disease reveals germline selection against severe mtDNA mutations. Science 319: 958–962.W. FanKG WaymireN. NarulaP. LiC. Rocher2008A mouse model of mitochondrial disease reveals germline selection against severe mtDNA mutations.Science319958962

19. Wallace DC (1992) Diseases of the mitochondrial DNA. Annual Review of Biochemistry 61: 1175–1212.DC Wallace1992Diseases of the mitochondrial DNA.Annual Review of Biochemistry6111751212

20. Kato T (2001) The other, forgotten genome: mitochondrial DNA and mental disorders. Mol Psychiatry 6: 625–633.T. Kato2001The other, forgotten genome: mitochondrial DNA and mental disorders.Mol Psychiatry6625633

21. Jazin EE, Cavelier L, Eriksson I, Oreland L, Gyllensten U (1996) Human brain contains high levels of heteroplasmy in the noncoding regions of mitochondrial DNA. Proc Natl Acad Sci U S A 93: 12382–12387.EE JazinL. CavelierI. ErikssonL. OrelandU. Gyllensten1996Human brain contains high levels of heteroplasmy in the noncoding regions of mitochondrial DNA.Proc Natl Acad Sci U S A931238212387

22. Cavelier L, Jazin E, Jalonen P, Gyllensten U (2000) MtDNA substitution rate and segregation of heteroplasmy in coding and noncoding regions. Hum Genet 107: 45–50.L. CavelierE. JazinP. JalonenU. Gyllensten2000MtDNA substitution rate and segregation of heteroplasmy in coding and noncoding regions.Hum Genet1074550

23. Jazin E, Soodyall H, Jalonen P, Lindholm E, Stoneking M, et al. (1998) Mitochondrial mutation rate revisited: hot spots and polymorphism. Nat Genet 18: 109–110.E. JazinH. SoodyallP. JalonenE. LindholmM. Stoneking1998Mitochondrial mutation rate revisited: hot spots and polymorphism.Nat Genet18109110

24. Lindholm E, Cavelier L, Howell WM, Eriksson I, Jalonen P, et al. (1997) Mitochondrial sequence variants in patients with schizophrenia. Eur J Hum Genet 5: 406–412.E. LindholmL. CavelierWM HowellI. ErikssonP. Jalonen1997Mitochondrial sequence variants in patients with schizophrenia.Eur J Hum Genet5406412

25. Cavelier L, Jazin EE, Eriksson I, Prince J, Bave U, et al. (1995) Decreased cytochrome-c oxidase activity and lack of age-related accumulation of mitochondrial DNA deletions in the brains of schizophrenics. Genomics 29: 217–224.L. CavelierEE JazinI. ErikssonJ. PrinceU. Bave1995Decreased cytochrome-c oxidase activity and lack of age-related accumulation of mitochondrial DNA deletions in the brains of schizophrenics.Genomics29217224

26. Coskun PE, Beal MF, Wallace DC (2004) Alzheimer's brains harbor somatic mtDNA control-region mutations that suppress mitochondrial transcription and replication. Proc Natl Acad Sci U S A 101: 10726–10731.PE CoskunMF BealDC Wallace2004Alzheimer's brains harbor somatic mtDNA control-region mutations that suppress mitochondrial transcription and replication.Proc Natl Acad Sci U S A1011072610731

27. Corral-Debrinski M, Horton T, Lott MT, Shoffner JM, McKee AC, et al. (1994) Marked changes in mitochondrial DNA deletion levels in Alzheimer brains. Genomics 23: 471–476.M. Corral-DebrinskiT. HortonMT LottJM ShoffnerAC McKee1994Marked changes in mitochondrial DNA deletion levels in Alzheimer brains.Genomics23471476

28. Fattal O, Budur K, Vaughan AJ, Franco K (2006) Review of the literature on major mental disorders in adult patients with mitochondrial diseases. Psychosomatics 47: 1–7.O. FattalK. BudurAJ VaughanK. Franco2006Review of the literature on major mental disorders in adult patients with mitochondrial diseases.Psychosomatics4717

29. Fattal O, Link J, Quinn K, Cohen BH, Franco K (2007) Psychiatric comorbidity in 36 adults with mitochondrial cytopathies. CNS Spectr 12: 429–438.O. FattalJ. LinkK. QuinnBH CohenK. Franco2007Psychiatric comorbidity in 36 adults with mitochondrial cytopathies.CNS Spectr12429438

30. DiMauro S, Moraes CT (1993) Mitochondrial encephalomyopathies. Archives of Neurology 50: 1197–1208.S. DiMauroCT Moraes1993Mitochondrial encephalomyopathies.Archives of Neurology5011971208

31. Gardner A, Johansson A, Wibom R, Nennesmo I, von Dobeln U, et al. (2003) Alterations of mitochondrial function and correlations with personality traits in selected major depressive disorder patients. J Affect Disord 76: 55–68.A. GardnerA. JohanssonR. WibomI. NennesmoU. von Dobeln2003Alterations of mitochondrial function and correlations with personality traits in selected major depressive disorder patients.J Affect Disord765568

32. Gardner A, Boles RG (2008) Symptoms of somatization as a rapid screening tool for mitochondrial dysfunction in depression. Biopsychosoc Med 2: 7.A. GardnerRG Boles2008Symptoms of somatization as a rapid screening tool for mitochondrial dysfunction in depression.Biopsychosoc Med27

33. Gardner A, Boles RG (2008) Mitochondrial energy depletion in depression with somatization. Psychother Psychosom 77: 127–129.A. GardnerRG Boles2008Mitochondrial energy depletion in depression with somatization.Psychother Psychosom77127129

34. Burnett BB, Gardner A, Boles RG (2005) Mitochondrial inheritance in depression, dysmotility and migraine? J Affect Disord 88: 109–116.BB BurnettA. GardnerRG Boles2005Mitochondrial inheritance in depression, dysmotility and migraine?J Affect Disord88109116

35. Steen RG, Hamer RM, Lieberman JA (2005) Measurement of brain metabolites by 1H magnetic resonance spectroscopy in patients with schizophrenia: a systematic review and meta-analysis. Neuropsychopharmacology 30: 1949–1962.RG SteenRM HamerJA Lieberman2005Measurement of brain metabolites by 1H magnetic resonance spectroscopy in patients with schizophrenia: a systematic review and meta-analysis.Neuropsychopharmacology3019491962

36. Buchsbaum MS, Hazlett EA (1998) Positron emission tomography studies of abnormal glucose metabolism in schizophrenia. Schizophr Bull 24: 343–364.MS BuchsbaumEA Hazlett1998Positron emission tomography studies of abnormal glucose metabolism in schizophrenia.Schizophr Bull24343364

37. Volz HP, Rzanny R, May S, Hegewald H, Preussler B, et al. (1997) 31P magnetic resonance spectroscopy in the dorsolateral prefrontal cortex of schizophrenics with a volume selective technique–preliminary findings. Biol Psychiatry 41: 644–648.HP VolzR. RzannyS. MayH. HegewaldB. Preussler199731P magnetic resonance spectroscopy in the dorsolateral prefrontal cortex of schizophrenics with a volume selective technique–preliminary findings.Biol Psychiatry41644648

38. Kegeles LS, Humaran TJ, Mann JJ (1998) In vivo neurochemistry of the brain in schizophrenia as revealed by magnetic resonance spectroscopy. Biol Psychiatry 44: 382–398.LS KegelesTJ HumaranJJ Mann1998In vivo neurochemistry of the brain in schizophrenia as revealed by magnetic resonance spectroscopy.Biol Psychiatry44382398

39. Yildiz-Yesiloglu A, Ankerst DP (2006) Neurochemical alterations of the brain in bipolar disorder and their implications for pathophysicology: A systematic review of the in vivo proton magnetic resonance spectroscopy findings. Prog Neuropsychopharmacol Biol Psychiatry 30: 969–995.A. Yildiz-YesilogluDP Ankerst2006Neurochemical alterations of the brain in bipolar disorder and their implications for pathophysicology: A systematic review of the in vivo proton magnetic resonance spectroscopy findings.Prog Neuropsychopharmacol Biol Psychiatry30969995

40. Strakowski SM, DelBello MP, Adler C, Cecil KM, Sax KW (2000) Neuroimaging in bipolar disorder. Bipolar Disord 2: 148–164.SM StrakowskiMP DelBelloC. AdlerKM CecilKW Sax2000Neuroimaging in bipolar disorder.Bipolar Disord2148164

41. Deicken RF, Fein G, Weiner MW (1995) Abnormal frontal lobe phosphorous metabolism in bipolar disorder. Am J Psychiatry 152: 915–918.RF DeickenG. FeinMW Weiner1995Abnormal frontal lobe phosphorous metabolism in bipolar disorder.Am J Psychiatry152915918

42. Hamakawa H, Kato T, Shioiri T, Inubushi T, Kato N (1999) Quantitative proton magnetic resonance spectroscopy of the bilateral frontal lobes in patients with bipolar disorder. Psychol Med 29: 639–644.H. HamakawaT. KatoT. ShioiriT. InubushiN. Kato1999Quantitative proton magnetic resonance spectroscopy of the bilateral frontal lobes in patients with bipolar disorder.Psychol Med29639644

43. Videbech P (2000) PET measurements of brain glucose metabolism and blood flow in major depressive disorder: a critical review. Acta Psychiatr Scand 101: 11–20.P. Videbech2000PET measurements of brain glucose metabolism and blood flow in major depressive disorder: a critical review.Acta Psychiatr Scand1011120

44. McMahon FJ, Chen YS, Patel S, Kokoszka J, Brown MD, et al. (2000) Mitochondrial DNA sequence diversity in bipolar affective disorder. Am J Psychiatry 157: 1058–1064.FJ McMahonYS ChenS. PatelJ. KokoszkaMD Brown2000Mitochondrial DNA sequence diversity in bipolar affective disorder.Am J Psychiatry15710581064

45. Iwamoto K, Bundo M, Kato T (2005) Altered expression of mitochondria-related genes in postmortem brains of patients with bipolar disorder or schizophrenia, as revealed by large-scale DNA microarray analysis. Hum Mol Genet 14: 241–253.K. IwamotoM. BundoT. Kato2005Altered expression of mitochondria-related genes in postmortem brains of patients with bipolar disorder or schizophrenia, as revealed by large-scale DNA microarray analysis.Hum Mol Genet14241253

46. Kato T, Kakiuchi C, Iwamoto K (2007) Comprehensive gene expression analysis in bipolar disorder. Can J Psychiatry 52: 763–771.T. KatoC. KakiuchiK. Iwamoto2007Comprehensive gene expression analysis in bipolar disorder.Can J Psychiatry52763771

47. Kato T (2008) Role of mitochondrial DNA in calcium signaling abnormality in bipolar disorder. Cell Calcium 44: 92–102.T. Kato2008Role of mitochondrial DNA in calcium signaling abnormality in bipolar disorder.Cell Calcium4492102

48. Munakata K, Fujii K, Nanko S, Kunugi H, Kato T (2007) Sequence and functional analyses of mtDNA in a maternally inherited family with bipolar disorder and depression. Mutat Res 617: 119–124.K. MunakataK. FujiiS. NankoH. KunugiT. Kato2007Sequence and functional analyses of mtDNA in a maternally inherited family with bipolar disorder and depression.Mutat Res617119124

49. Kato T (2006) The role of mitochondrial dysfunction in bipolar disorder. Drug News Perspect 19: 597–602.T. Kato2006The role of mitochondrial dysfunction in bipolar disorder.Drug News Perspect19597602

50. Washizuka S, Kakiuchi C, Mori K, Tajima O, Akiyama T, et al. (2005) Expression of mitochondria-related genes in lymphoblastoid cells from patients with bipolar disorder. Bipolar Disord 7: 146–152.S. WashizukaC. KakiuchiK. MoriO. TajimaT. Akiyama2005Expression of mitochondria-related genes in lymphoblastoid cells from patients with bipolar disorder.Bipolar Disord7146152

51. Munakata K, Iwamoto K, Bundo M, Kato T (2005) Mitochondrial DNA 3243A>G mutation and increased expression of LARS2 gene in the brains of patients with bipolar disorder and schizophrenia. Biol Psychiatry 57: 525–532.K. MunakataK. IwamotoM. BundoT. Kato2005Mitochondrial DNA 3243A>G mutation and increased expression of LARS2 gene in the brains of patients with bipolar disorder and schizophrenia.Biol Psychiatry57525532

52. Munakata K, Tanaka M, Mori K, Washizuka S, Yoneda M, et al. (2004) Mitochondrial DNA 3644T–>C mutation associated with bipolar disorder. Genomics 84: 1041–1050.K. MunakataM. TanakaK. MoriS. WashizukaM. Yoneda2004Mitochondrial DNA 3644T–>C mutation associated with bipolar disorder.Genomics8410411050

53. Washizuka S, Kametani M, Sasaki T, Tochigi M, Umekage T, et al. (2006) Association of mitochondrial complex I subunit gene NDUFV2 at 18p11 with schizophrenia in the Japanese population. Am J Med Genet B Neuropsychiatr Genet 141: 301–304.S. WashizukaM. KametaniT. SasakiM. TochigiT. Umekage2006Association of mitochondrial complex I subunit gene NDUFV2 at 18p11 with schizophrenia in the Japanese population.Am J Med Genet B Neuropsychiatr Genet141301304

54. Kakiuchi C, Ishiwata M, Kametani M, Nelson C, Iwamoto K, et al. (2005) Quantitative analysis of mitochondrial DNA deletions in the brains of patients with bipolar disorder and schizophrenia. Int J Neuropsychopharmacol 8: 515–522.C. KakiuchiM. IshiwataM. KametaniC. NelsonK. Iwamoto2005Quantitative analysis of mitochondrial DNA deletions in the brains of patients with bipolar disorder and schizophrenia.Int J Neuropsychopharmacol8515522

55. Kato T, Takahashi Y (1996) Deletion of leukocyte mitochondrial DNA in bipolar disorder. J Affect Disord 37: 67–73.T. KatoY. Takahashi1996Deletion of leukocyte mitochondrial DNA in bipolar disorder.J Affect Disord376773

56. Kato T, Stine OC, McMahon FJ, Crowe RR (1997) Increased levels of a mitochondrial DNA deletion in the brain of patients with bipolar disorder. Biol Psychiatry 42: 871–875.T. KatoOC StineFJ McMahonRR Crowe1997Increased levels of a mitochondrial DNA deletion in the brain of patients with bipolar disorder.Biol Psychiatry42871875

57. Washizuka S, Kakiuchi C, Mori K, Kunugi H, Tajima O, et al. (2003) Association of mitochondrial complex I subunit gene NDUFV2 at 18p11 with bipolar disorder. Am J Med Genet B Neuropsychiatr Genet 120: 72–78.S. WashizukaC. KakiuchiK. MoriH. KunugiO. Tajima2003Association of mitochondrial complex I subunit gene NDUFV2 at 18p11 with bipolar disorder.Am J Med Genet B Neuropsychiatr Genet1207278

58. Bamne MN, Talkowski ME, Moraes CT, Manuck SB, Ferrell RE, et al. (2008) Systematic Association Studies of Mitochondrial DNA Variations in Schizophrenia: Focus on the ND5 Gene. Schizophr Bull 34: 458–465.MN BamneME TalkowskiCT MoraesSB ManuckRE Ferrell2008Systematic Association Studies of Mitochondrial DNA Variations in Schizophrenia: Focus on the ND5 Gene.Schizophr Bull34458465

59. Martorell L, Segues T, Folch G, Valero J, Joven J, et al. (2006) New variants in the mitochondrial genomes of schizophrenic patients. Eur J Hum Genet 14: 520–528.L. MartorellT. SeguesG. FolchJ. ValeroJ. Joven2006New variants in the mitochondrial genomes of schizophrenic patients.Eur J Hum Genet14520528

60. Marchbanks RM, Ryan M, Day IN, Owen M, McGuffin P, et al. (2003) A mitochondrial DNA sequence variant associated with schizophrenia and oxidative stress. Schizophr Res 65: 33–38.RM MarchbanksM. RyanIN DayM. OwenP. McGuffin2003A mitochondrial DNA sequence variant associated with schizophrenia and oxidative stress.Schizophr Res653338

61. Washizuka S, Ikeda A, Kato N, Kato T (2003) Possible relationship between mitochondrial DNA polymorphisms and lithium response in bipolar disorder. Int J Neuropsychopharmacol 6: 421–424.S. WashizukaA. IkedaN. KatoT. Kato2003Possible relationship between mitochondrial DNA polymorphisms and lithium response in bipolar disorder.Int J Neuropsychopharmacol6421424

62. Kazuno AA, Munakata K, Nagai T, Shimozono S, Tanaka M, et al. (2006) Identification of mitochondrial DNA polymorphisms that alter mitochondrial matrix pH and intracellular calcium dynamics. PLoS Genet 2: e128.AA KazunoK. MunakataT. NagaiS. ShimozonoM. Tanaka2006Identification of mitochondrial DNA polymorphisms that alter mitochondrial matrix pH and intracellular calcium dynamics.PLoS Genet2e128

63. Kazuno AA, Munakata K, Mori K, Tanaka M, Nanko S, et al. (2005) Mitochondrial DNA sequence analysis of patients with ‘atypical psychosis’. Psychiatry Clin Neurosci 59: 497–503.AA KazunoK. MunakataK. MoriM. TanakaS. Nanko2005Mitochondrial DNA sequence analysis of patients with ‘atypical psychosis’.Psychiatry Clin Neurosci59497503

64. Amar S, Shamir A, Ovadia O, Blanaru M, Reshef A, et al. (2007) Mitochondrial DNA HV lineage increases the susceptibility to schizophrenia among Israeli Arabs. Schizophr Res 94: 354–358.S. AmarA. ShamirO. OvadiaM. BlanaruA. Reshef2007Mitochondrial DNA HV lineage increases the susceptibility to schizophrenia among Israeli Arabs.Schizophr Res94354358

65. Coon KD, Valla J, Szelinger S, Schneider LE, Niedzielko TL, et al. (2006) Quantitation of heteroplasmy of mtDNA sequence variants identified in a population of AD patients and controls by array-based resequencing. Mitochondrion 6: 194–210.KD CoonJ. VallaS. SzelingerLE SchneiderTL Niedzielko2006Quantitation of heteroplasmy of mtDNA sequence variants identified in a population of AD patients and controls by array-based resequencing.Mitochondrion6194210

66. Mexal S, Berger R, Adams C, Ross R, Freedman R, et al. (2006) Brain pH has a significant impact on human postmortem hippocampal gene expression profiles. Brain Res 1106: 1–11.S. MexalR. BergerC. AdamsR. RossR. Freedman2006Brain pH has a significant impact on human postmortem hippocampal gene expression profiles.Brain Res1106111

67. Vawter MP, Tomita H, Meng F, Bolstad B, Li J, et al. (2006) Mitochondrial-related gene expression changes are sensitive to agonal-pH state: implications for brain disorders. Mol Psychiatry 11: 615, 663–679.MP VawterH. TomitaF. MengB. BolstadJ. Li2006Mitochondrial-related gene expression changes are sensitive to agonal-pH state: implications for brain disorders.Mol Psychiatry11615, 663679

68. Kato T, Kato N (2000) Mitochondrial dysfunction in bipolar disorder. Bipolar Disord 2: 180–190.T. KatoN. Kato2000Mitochondrial dysfunction in bipolar disorder.Bipolar Disord2180190

69. Hartmann A, Thieme M, Nanduri LK, Stempfl T, Moehle C, et al. (2008) Validation of microarray-based resequencing of 93 worldwide mitochondrial genomes. Hum Mutat 30: 115–122.A. HartmannM. ThiemeLK NanduriT. StempflC. Moehle2008Validation of microarray-based resequencing of 93 worldwide mitochondrial genomes.Hum Mutat30115122

70. Jakupciak JP, Maggrah A, Maragh S, Maki J, Reguley B, et al. (2008) Facile whole mitochondrial genome resequencing from nipple aspirate fluid using MitoChip v2.0. BMC Cancer 8: 95.JP JakupciakA. MaggrahS. MaraghJ. MakiB. Reguley2008Facile whole mitochondrial genome resequencing from nipple aspirate fluid using MitoChip v2.0.BMC Cancer895

71. Zhou S, Kassauei K, Cutler DJ, Kennedy GC, Sidransky D, et al. (2006) An oligonucleotide microarray for high-throughput sequencing of the mitochondrial genome. J Mol Diagn 8: 476–482.S. ZhouK. KassaueiDJ CutlerGC KennedyD. Sidransky2006An oligonucleotide microarray for high-throughput sequencing of the mitochondrial genome.J Mol Diagn8476482

72. Jakupciak JP, Wang W, Markowitz ME, Ally D, Coble M, et al. (2005) Mitochondrial DNA as a cancer biomarker. J Mol Diagn 7: 258–267.JP JakupciakW. WangME MarkowitzD. AllyM. Coble2005Mitochondrial DNA as a cancer biomarker.J Mol Diagn7258267

73. Jacobs L, Gerards M, Chinnery P, Dumoulin J, de Coo I, et al. (2007) mtDNA point mutations are present at various levels of heteroplasmy in human oocytes. Mol Hum Reprod 13: 149–154.L. JacobsM. GerardsP. ChinneryJ. DumoulinI. de Coo2007mtDNA point mutations are present at various levels of heteroplasmy in human oocytes.Mol Hum Reprod13149154

74. Leveque M, Marlin S, Jonard L, Procaccio V, Reynier P, et al. (2007) Whole mitochondrial genome screening in maternally inherited non-syndromic hearing impairment using a microarray resequencing mitochondrial DNA chip. Eur J Hum Genet 15: 1145–1155.M. LevequeS. MarlinL. JonardV. ProcaccioP. Reynier2007Whole mitochondrial genome screening in maternally inherited non-syndromic hearing impairment using a microarray resequencing mitochondrial DNA chip.Eur J Hum Genet1511451155

75. Maitra A, Cohen Y, Gillespie SE, Mambo E, Fukushima N, et al. (2004) The Human MitoChip: a high-throughput sequencing microarray for mitochondrial mutation detection. Genome Res 14: 812–819.A. MaitraY. CohenSE GillespieE. MamboN. Fukushima2004The Human MitoChip: a high-throughput sequencing microarray for mitochondrial mutation detection.Genome Res14812819

76. Huelsenbeck JP, Ronquist F (2001) MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 17: 754–755.JP HuelsenbeckF. Ronquist2001MRBAYES: Bayesian inference of phylogenetic trees.Bioinformatics17754755

77. Huelsenbeck JP, Ronquist F, Nielsen R, Bollback JP (2001) Bayesian inference of phylogeny and its impact on evolutionary biology. Science 294: 2310–2314.JP HuelsenbeckF. RonquistR. NielsenJP Bollback2001Bayesian inference of phylogeny and its impact on evolutionary biology.Science29423102314

78. Howell N, Elson JL, Turnbull DM, Herrnstadt C (2004) African Haplogroup L mtDNA sequences show violations of clock-like evolution. Mol Biol Evol 21: 1843–1854.N. HowellJL ElsonDM TurnbullC. Herrnstadt2004African Haplogroup L mtDNA sequences show violations of clock-like evolution.Mol Biol Evol2118431854

79. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.S. PurcellB. NealeK. Todd-BrownL. ThomasMA Ferreira2007PLINK: a tool set for whole-genome association and population-based linkage analyses.Am J Hum Genet81559575

80. Ingman M, Gyllensten U (2006) mtDB: Human Mitochondrial Genome Database, a resource for population genetics and medical sciences. Nucleic Acids Res 34: D749–751.M. IngmanU. Gyllensten2006mtDB: Human Mitochondrial Genome Database, a resource for population genetics and medical sciences.Nucleic Acids Res34D749751

81. Bannwarth S, Procaccio V, Paquis-Flucklinger V (2005) Surveyor Nuclease: a new strategy for a rapid identification of heteroplasmic mitochondrial DNA mutations in patients with respiratory chain defects. Hum Mutat 25: 575–582.S. BannwarthV. ProcaccioV. Paquis-Flucklinger2005Surveyor Nuclease: a new strategy for a rapid identification of heteroplasmic mitochondrial DNA mutations in patients with respiratory chain defects.Hum Mutat25575582

82. Bannwarth S, Procaccio V, Paquis-Flucklinger V (2006) Rapid identification of unknown heteroplasmic mutations across the entire human mitochondrial genome with mismatch-specific Surveyor Nuclease. Nat Protoc 1: 2037–2047.S. BannwarthV. ProcaccioV. Paquis-Flucklinger2006Rapid identification of unknown heteroplasmic mutations across the entire human mitochondrial genome with mismatch-specific Surveyor Nuclease.Nat Protoc120372047

83. King MP, Attardi G (1989) Human cells lacking mtDNA: repopulation with exogenous mitochondria by complementation. Science 246: 500–503.MP KingG. Attardi1989Human cells lacking mtDNA: repopulation with exogenous mitochondria by complementation.Science246500503

84. Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 24: 1596–1599.K. TamuraJ. DudleyM. NeiS. Kumar2007MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0.Mol Biol Evol2415961599

85. Nei M, Gojobori T (1986) Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol 3: 418–426.M. NeiT. Gojobori1986Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions.Mol Biol Evol3418426

86. Maragh S, Jakupciak JP, Wagner PD, Rom WN, Sidransky D, et al. (2008) Multiple strand displacement amplification of mitochondrial DNA from clinical samples. BMC Med Genet 9: 7.S. MaraghJP JakupciakPD WagnerWN RomD. Sidransky2008Multiple strand displacement amplification of mitochondrial DNA from clinical samples.BMC Med Genet97

87. Lin MT, Simon DK, Ahn CH, Kim LM, Beal MF (2002) High aggregate burden of somatic mtDNA point mutations in aging and Alzheimer's disease brain. Hum Mol Genet 11: 133–145.MT LinDK SimonCH AhnLM KimMF Beal2002High aggregate burden of somatic mtDNA point mutations in aging and Alzheimer's disease brain.Hum Mol Genet11133145

88. Davidzon G, Greene P, Mancuso M, Klos KJ, Ahlskog JE, et al. (2006) Early-onset familial parkinsonism due to POLG mutations. Ann Neurol 59: 859–862.G. DavidzonP. GreeneM. MancusoKJ KlosJE Ahlskog2006Early-onset familial parkinsonism due to POLG mutations.Ann Neurol59859862

89. Luoma PT, Eerola J, Ahola S, Hakonen AH, Hellstrom O, et al. (2007) Mitochondrial DNA polymerase gamma variants in idiopathic sporadic Parkinson disease. Neurology 69: 1152–1159.PT LuomaJ. EerolaS. AholaAH HakonenO. Hellstrom2007Mitochondrial DNA polymerase gamma variants in idiopathic sporadic Parkinson disease.Neurology6911521159

90. Stewart JB, Freyer C, Elson JL, Wredenberg A, Cansu Z, et al. (2008) Strong purifying selection in transmission of mammalian mitochondrial DNA. PLoS Biol 6: e10.JB StewartC. FreyerJL ElsonA. WredenbergZ. Cansu2008Strong purifying selection in transmission of mammalian mitochondrial DNA.PLoS Biol6e10

91. Pyle A, Foltynie T, Tiangyou W, Lambert C, Keers SM, et al. (2005) Mitochondrial DNA haplogroup cluster UKJT reduces the risk of PD. Ann Neurol 57: 564–567.A. PyleT. FoltynieW. TiangyouC. LambertSM Keers2005Mitochondrial DNA haplogroup cluster UKJT reduces the risk of PD.Ann Neurol57564567

92. Gibson GE, Ratan RR, Beal MF (2008) Mitochondria and oxidative stress in neurodegenerative disorders. Preface. Ann N Y Acad Sci 1147: xi–xii.GE GibsonRR RatanMF Beal2008Mitochondria and oxidative stress in neurodegenerative disorders. Preface.Ann N Y Acad Sci1147xixii

93. Shao L, Martin MV, Watson SJ, Schatzberg A, Akil H, et al. (2008) Mitochondrial involvement in psychiatric disorders. Ann Med 40: 281–295.L. ShaoMV MartinSJ WatsonA. SchatzbergH. Akil2008Mitochondrial involvement in psychiatric disorders.Ann Med40281295

94. Ben-Shachar D, Karry R (2008) Neuroanatomical pattern of mitochondrial complex I pathology varies between schizophrenia, bipolar disorder and major depression. PLoS ONE 3: e3676.D. Ben-ShacharR. Karry2008Neuroanatomical pattern of mitochondrial complex I pathology varies between schizophrenia, bipolar disorder and major depression.PLoS ONE3e3676

95. Balijepalli S, Boyd MR, Ravindranath V (1999) Inhibition of mitochondrial complex I by haloperidol: the role of thiol oxidation. Neuropharmacology 38: 567–577.S. BalijepalliMR BoydV. Ravindranath1999Inhibition of mitochondrial complex I by haloperidol: the role of thiol oxidation.Neuropharmacology38567577

96. Whatley SA, Curti D, Das Gupta F, Ferrier IN, Jones S, et al. (1998) Superoxide, neuroleptics and the ubiquinone and cytochrome b5 reductases in brain and lymphocytes from normals and schizophrenic patients. Mol Psychiatry 3: 227–237.SA WhatleyD. CurtiF. Das GuptaIN FerrierS. Jones1998Superoxide, neuroleptics and the ubiquinone and cytochrome b5 reductases in brain and lymphocytes from normals and schizophrenic patients.Mol Psychiatry3227237

97. Droge W, Schipper HM (2007) Oxidative stress and aberrant signaling in aging and cognitive decline. Aging Cell 6: 361–370.W. DrogeHM Schipper2007Oxidative stress and aberrant signaling in aging and cognitive decline.Aging Cell6361370

98. Wester P, Bateman DE, Dodd PR, Edwardson JA, Hardy JA, et al. (1985) Agonal status affects the metabolic activity of nerve endings isolated from postmortem human brain. Neurochem Pathol 3: 169–180.P. WesterDE BatemanPR DoddJA EdwardsonJA Hardy1985Agonal status affects the metabolic activity of nerve endings isolated from postmortem human brain.Neurochem Pathol3169180

99. Barton AJ, Pearson RC, Najlerahim A, Harrison PJ (1993) Pre- and postmortem influences on brain RNA. J Neurochem 61: 1–11.AJ BartonRC PearsonA. NajlerahimPJ Harrison1993Pre- and postmortem influences on brain RNA.J Neurochem61111

100. Johnston NL, Cervenak J, Shore AD, Torrey EF, Yolken RH (1997) Multivariate analysis of RNA levels from postmortem human brains as measured by three different methods of RT-PCR. Stanley Neuropathology Consortium. J Neurosci Methods 77: 83–92.NL JohnstonJ. CervenakAD ShoreEF TorreyRH Yolken1997Multivariate analysis of RNA levels from postmortem human brains as measured by three different methods of RT-PCR. Stanley Neuropathology Consortium.J Neurosci Methods778392

101. Harrison PJ, Procter AW, Barton AJ, Lowe SL, Najlerahim A, et al. (1991) Terminal coma affects messenger RNA detection in post mortem human temporal cortex. Brain Res Mol Brain Res 9: 161–164.PJ HarrisonAW ProcterAJ BartonSL LoweA. Najlerahim1991Terminal coma affects messenger RNA detection in post mortem human temporal cortex.Brain Res Mol Brain Res9161164

102. Vawter MP, Evans S, Choudary P, Tomita H, Meador-Woodruff J, et al. (2004) Gender-specific gene expression in post-mortem human brain: localization to sex chromosomes. Neuropsychopharmacology 29: 373–384.MP VawterS. EvansP. ChoudaryH. TomitaJ. Meador-Woodruff2004Gender-specific gene expression in post-mortem human brain: localization to sex chromosomes.Neuropsychopharmacology29373384

103. Mishmar D, Ruiz-Pesini EE, Golik P, Macaulay V, Clark AG, et al. (2003) Natural selection shaped regional mtDNA variation in humans. Proc Natl Acad Sci U S A 100: 171–176.D. MishmarEE Ruiz-PesiniP. GolikV. MacaulayAG Clark2003Natural selection shaped regional mtDNA variation in humans.Proc Natl Acad Sci U S A100171176

104. Brandon MC, Ruiz-Pesini E, Mishmar D, Procaccio V, Lott MT, et al. (2008) MITOMASTER: a bioinformatics tool for the analysis of mitochondrial DNA sequences. Hum Mutat 30: 1–6.MC BrandonE. Ruiz-PesiniD. MishmarV. ProcaccioMT Lott2008MITOMASTER: a bioinformatics tool for the analysis of mitochondrial DNA sequences.Hum Mutat3016

105. Ruiz-Pesini E, Lott MT, Procaccio V, Poole JC, Brandon MC, et al. (2007) An enhanced MITOMAP with a global mtDNA mutational phylogeny. Nucleic Acids Res 35: D823–828.E. Ruiz-PesiniMT LottV. ProcaccioJC PooleMC Brandon2007An enhanced MITOMAP with a global mtDNA mutational phylogeny.Nucleic Acids Res35D823828

106. Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, et al. (1999) Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet 23: 147.RM AndrewsI. KubackaPF ChinneryRN LightowlersDM Turnbull1999Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA.Nat Genet23147

107. Strand H, Ingebretsen OC, Nilssen O (2008) Real-time detection and quantification of mitochondrial mutations with oligonucleotide primers containing locked nucleic acid. Clin Chim Acta 390: 126–133.H. StrandOC IngebretsenO. Nilssen2008Real-time detection and quantification of mitochondrial mutations with oligonucleotide primers containing locked nucleic acid.Clin Chim Acta390126133

108. Newton CR, Graham A, Heptinstall LE, Powell SJ, Summers C, et al. (1989) Analysis of any point mutation in DNA. The amplification refractory mutation system (ARMS). Nucleic Acids Res 17: 2503–2516.CR NewtonA. GrahamLE HeptinstallSJ PowellC. Summers1989Analysis of any point mutation in DNA. The amplification refractory mutation system (ARMS).Nucleic Acids Res1725032516

109. Ruiz-Pesini E, Mishmar D, Brandon M, Procaccio V, Wallace DC (2004) Effects of purifying and adaptive selection on regional variation in human mtDNA. Science 303: 223–226.E. Ruiz-PesiniD. MishmarM. BrandonV. ProcaccioDC Wallace2004Effects of purifying and adaptive selection on regional variation in human mtDNA.Science303223226