Mitochondrial Genome Analysis of Primary Open Angle Glaucoma Patients

Primary open angle glaucoma (POAG) is a multi-factorial optic disc neuropathy characterized by accelerating damage of the retinal ganglion cells and atrophy of the optic nerve head. The vulnerability of the optic nerve damage leading to POAG has been postulated to result from oxidative stress and mitochondrial dysfunction. In this study, we investigated the possible involvement of the mitochondrial genomic variants in 101 patients and 71 controls by direct sequencing of the entire mitochondrial genome. The number of variable positions in the mtDNA with respect to the revised Cambridge Reference Sequence (rCRS), have been designated “Segregating Sites”. The segregating sites present only in the patients or controls have been designated “Unique Segregating Sites (USS)”. The population mutation rate (θ = 4N_eμ) as estimated by Watterson’s θ (θ_w), considering only the USS, was significantly higher among the patients (p = 9.8×10⁻¹⁵) compared to controls. The difference in θ_w and the number of USS were more pronounced when restricted to the coding region (p<1.31×10⁻²¹ and p = 0.006607, respectively). Further analysis of the region revealed non-synonymous variations were significantly higher in Complex I among the patients (p = 0.0053). Similar trends were retained when USS was considered only within complex I (frequency 0.49 vs 0.31 with p<0.0001 and mutation rate p-value <1.49×10⁻⁴³) and ND5 within its gene cluster (frequency 0.47 vs 0.23 with p<0.0001 and mutation rate p-value <4.42×10⁻⁴⁷). ND5 is involved in the proton pumping mechanism. Incidentally, glaucomatous trabecular meshwork cells have been reported to be more sensitive to inhibition of complex I activity. Thus mutations in ND5, expected to inhibit complex I activity, could lead to generation of oxidative stress and favor glaucomatous condition.     

Mitochondrial DNA Variant Discovery and Evaluation in Human Cardiomyopathies through Next-Generation Sequencing

Mutations in mitochondrial DNA (mtDNA) may cause maternally-inherited cardiomyopathy and heart failure. In homoplasmy all mtDNA copies contain the mutation. In heteroplasmy there is a mixture of normal and mutant copies of mtDNA. The clinical phenotype of an affected individual depends on the type of genetic defect and the ratios of mutant and normal mtDNA in affected tissues. We aimed at determining the sensitivity of next-generation sequencing compared to Sanger sequencing for mutation detection in patients with mitochondrial cardiomyopathy. We studied 18 patients with mitochondrial cardiomyopathy and two with suspected mitochondrial disease. We “shotgun” sequenced PCR-amplified mtDNA and multiplexed using a single run on Roche''s 454 Genome Sequencer. By mapping to the reference sequence, we obtained 1,300× average coverage per case and identified high-confidence variants. By comparing these to >400 mtDNA substitution variants detected by Sanger, we found 98% concordance in variant detection. Simulation studies showed that >95% of the homoplasmic variants were detected at a minimum sequence coverage of 20× while heteroplasmic variants required >200× coverage. Several Sanger “misses” were detected by 454 sequencing. These included the novel heteroplasmic 7501T>C in tRNA serine 1 in a patient with sudden cardiac death. These results support a potential role of next-generation sequencing in the discovery of novel mtDNA variants with heteroplasmy below the level reliably detected with Sanger sequencing. We hope that this will assist in the identification of mtDNA mutations and key genetic determinants for cardiomyopathy and mitochondrial disease.     

Variants in the ASB10 Gene Are Associated with Primary Open Angle Glaucoma

Background

Recently nonsynonymous coding variants in the ankyrin repeats and suppressor of cytokine signaling box-containing protein 10 (ASB10) gene were found to be associated with primary open angle glaucoma (POAG) in cohorts from Oregon and Germany, but this finding was not confirmed in an independent cohort from Iowa. The aim of the current study was to assess the role of ASB10 gene variants in Pakistani glaucoma patients.

Methods

Sanger sequencing of the coding exons and splice junctions of the ASB10 gene was performed in 30 probands of multiplex POAG families, 208 sporadic POAG patients and 151 healthy controls from Pakistan. Genotypic associations of individual variants with POAG were analyzed with the Fisher’s exact or Chi-square test.

Results

In total 24 variants were identified in POAG probands and sporadic patients, including 11 novel variants and 13 known variants. 13 of the variants were nonsynonymous, 6 were synonymous, and 5 were intronic. Three nonsynonymous variants (p.Arg49Cys, p.Arg237Gly, p.Arg453Cys) identified in the probands were not segregating in the respective families. This is not surprising since glaucoma is a multifactorial disease, and multiple factors are likely to be involved in the disease manifestation in these families. However a nonsynonymous variant, p.Arg453Cys (rs3800791), was found in 6 sporadic POAG patients but not in controls, suggesting that it infers increased risk for the disease. In addition, one synonymous variant was found to be associated with sporadic POAG: p.Ala290Ala and the association of the variant with POAG remained significant after correction for multiple testing (uncorrected p-value 0.002, corrected p-value 0.047). The cumulative burden of rare, nonsynonymous variants was significantly higher in sporadic POAG patients compared to control individuals (p-value 0.000006).

Conclusions

Variants in ASB10 were found to be significantly associated with sporadic POAG in the Pakistani population. This supports previous findings that sequence variants in the ASB10 gene may act as a risk factor for glaucoma.     

Mitochondrial Genomic Analysis of Late Onset Alzheimer’s Disease Reveals Protective Haplogroups H6A1A/H6A1B: The Cache County Study on Memory in Aging

Background

Alzheimer’s disease (AD) is the most common cause of dementia and AD risk clusters within families. Part of the familial aggregation of AD is accounted for by excess maternal vs. paternal inheritance, a pattern consistent with mitochondrial inheritance. The role of specific mitochondrial DNA (mtDNA) variants and haplogroups in AD risk is uncertain.

Methodology/Principal Findings

We determined the complete mitochondrial genome sequence of 1007 participants in the Cache County Study on Memory in Aging, a population-based prospective cohort study of dementia in northern Utah. AD diagnoses were made with a multi-stage protocol that included clinical examination and review by a panel of clinical experts. We used TreeScanning, a statistically robust approach based on haplotype networks, to analyze the mtDNA sequence data. Participants with major mitochondrial haplotypes H6A1A and H6A1B showed a reduced risk of AD (p = 0.017, corrected for multiple comparisons). The protective haplotypes were defined by three variants: m.3915G>A, m.4727A>G, and m.9380G>A. These three variants characterize two different major haplogroups. Together m.4727A>G and m.9380G>A define H6A1, and it has been suggested m.3915G>A defines H6A. Additional variants differentiate H6A1A and H6A1B; however, none of these variants had a significant relationship with AD case-control status.

Conclusions/Significance

Our findings provide evidence of a reduced risk of AD for individuals with mtDNA haplotypes H6A1A and H6A1B. These findings are the results of the largest study to date with complete mtDNA genome sequence data, yet the functional significance of the associated haplotypes remains unknown and replication in others studies is necessary.     

Discovery and Functional Annotation of SIX6 Variants in Primary Open-Angle Glaucoma

Glaucoma is a leading cause of blindness worldwide. Primary open-angle glaucoma (POAG) is the most common subtype and is a complex trait with multigenic inheritance. Genome-wide association studies have previously identified a significant association between POAG and the SIX6 locus (rs10483727, odds ratio (OR) = 1.32, p = 3.87×10⁻¹¹). SIX6 plays a role in ocular development and has been associated with the morphology of the optic nerve. We sequenced the SIX6 coding and regulatory regions in 262 POAG cases and 256 controls and identified six nonsynonymous coding variants, including five rare and one common variant, Asn141His (rs33912345), which was associated significantly with POAG (OR = 1.27, p = 4.2×10⁻¹⁰) in the NEIGHBOR/GLAUGEN datasets. These variants were tested in an in vivo Danio rerio (zebrafish) complementation assay to evaluate ocular metrics such as eye size and optic nerve structure. Five variants, found primarily in POAG cases, were hypomorphic or null, while the sixth variant, found only in controls, was benign. One variant in the SIX6 enhancer increased expression of SIX6 and disrupted its regulation. Finally, to our knowledge for the first time, we have identified a clinical feature in POAG patients that appears to be dependent upon SIX6 genotype: patients who are homozygous for the SIX6 risk allele (His141) have a statistically thinner retinal nerve fiber layer than patients homozygous for the SIX6 non-risk allele (Asn141). Our results, in combination with previous SIX6 work, lead us to hypothesize that SIX6 risk variants disrupt the development of the neural retina, leading to a reduced number of retinal ganglion cells, thereby increasing the risk of glaucoma-associated vision loss.     

Explicit and Implicit Emotional Expression in Bulimia Nervosa in the Acute State and after Recovery

Expression of emotional state is considered to be a core facet of an individual''s emotional competence. Emotional processing in BN has not been often studied and has not been considered from a broad perspective. This study aimed at examining the implicit and explicit emotional expression in BN patients, in the acute state and after recovery. Sixty-three female participants were included: 22 BN, 22 recovered BN (R-BN), and 19 healthy controls (HC). The clinical cases were drawn from consecutive admissions and diagnosed according to DSM-IV-TR diagnostic criteria. Self reported (explicit) emotional expression was measured with State-Trait Anger Expression Inventory-2, State-Trait Anxiety Inventory, and Symptom Check List-90 items-Revised. Emotional facial expression (implicit) was recorded by means of an integrated camera (by detecting Facial Feature Tracking), during a 20 minutes therapeutic video game. In the acute illness explicit emotional expression [anxiety (p<0.001) and anger (p<0.05)] was increased. In the recovered group this was decreased to an intermediate level between the acute illness and healthy controls [anxiety (p<0.001) and anger (p<0.05)]. In the implicit measurement of emotional expression patients with acute BN expressed more joy (p<0.001) and less anger (p<0.001) than both healthy controls and those in the recovered group. These findings suggest that there are differences in the implicit and explicit emotional processing in BN, which is significantly reduced after recovery, suggesting an improvement in emotional regulation.     

Stable heteroplasmy at the single-cell level is facilitated by intercellular exchange of mtDNA

Eukaryotic cells carry two genomes, nuclear (nDNA) and mitochondrial (mtDNA), which are ostensibly decoupled in their replication, segregation and inheritance. It is increasingly appreciated that heteroplasmy, the occurrence of multiple mtDNA haplotypes in a cell, plays an important biological role, but its features are not well understood. Accurately determining the diversity of mtDNA has been difficult, due to the relatively small amount of mtDNA in each cell (<1% of the total DNA), the intercellular variability of mtDNA content and mtDNA pseudogenes (Numts) in nDNA. To understand the nature of heteroplasmy, we developed Mseek, a novel technique to purify and sequence mtDNA. Mseek yields high purity (>90%) mtDNA and its ability to detect rare variants is limited only by sequencing depth, providing unprecedented sensitivity and specificity. Using Mseek, we confirmed the ubiquity of heteroplasmy by analyzing mtDNA from a diverse set of cell lines and human samples. Applying Mseek to colonies derived from single cells, we find heteroplasmy is stably maintained in individual daughter cells over multiple cell divisions. We hypothesized that the stability of heteroplasmy could be facilitated by intercellular exchange of mtDNA. We explicitly demonstrate this exchange by co-culturing cell lines with distinct mtDNA haplotypes. Our results shed new light on the maintenance of heteroplasmy and provide a novel platform to investigate features of heteroplasmy in normal and diseased states.     

Abnormal Mitochondrial L-Arginine Transport Contributes to the Pathogenesis of Heart Failure and Rexoygenation Injury

Background

Impaired mitochondrial function is fundamental feature of heart failure (HF) and myocardial ischemia. In addition to the effects of heightened oxidative stress, altered nitric oxide (NO) metabolism, generated by a mitochondrial NO synthase, has also been proposed to impact upon mitochondrial function. However, the mechanism responsible for arginine transport into mitochondria and the effect of HF on such a process is unknown. We therefore aimed to characterize mitochondrial L-arginine transport and to investigate the hypothesis that impaired mitochondrial L-arginine transport plays a key role in the pathogenesis of heart failure and myocardial injury.

Methods and Results

In mitochondria isolated from failing hearts (sheep rapid pacing model and mouse Mst1 transgenic model) we demonstrated a marked reduction in L-arginine uptake (p<0.05 and p<0.01 respectively) and expression of the principal L-arginine transporter, CAT-1 (p<0.001, p<0.01) compared to controls. This was accompanied by significantly lower NO production and higher 3-nitrotyrosine levels (both p<0.05). The role of mitochondrial L-arginine transport in modulating cardiac stress responses was examined in cardiomyocytes with mitochondrial specific overexpression of CAT-1 (mtCAT1) exposed to hypoxia-reoxygenation stress. mtCAT1 cardiomyocytes had significantly improved mitochondrial membrane potential, respiration and ATP turnover together with significantly decreased reactive oxygen species production and cell death following mitochondrial stress.

Conclusion

These data provide new insights into the role of L-arginine transport in mitochondrial biology and cardiovascular disease. Augmentation of mitochondrial L-arginine availability may be a novel therapeutic strategy for myocardial disorders involving mitochondrial stress such as heart failure and reperfusion injury.     

Private Mitochondrial DNA Variants in Danish Patients with Hypertrophic Cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is a genetic cardiac disease primarily caused by mutations in genes coding for sarcomeric proteins. A molecular-genetic etiology can be established in ~60% of cases. Evolutionarily conserved mitochondrial DNA (mtDNA) haplogroups are susceptibility factors for HCM. Several polymorphic mtDNA variants are associated with a variety of late-onset degenerative diseases and affect mitochondrial function. We examined the role of private, non-haplogroup associated, mitochondrial variants in the etiology of HCM. In 87 Danish HCM patients, full mtDNA sequencing revealed 446 variants. After elimination of 312 (69.9%) non-coding and synonymous variants, a further 109 (24.4%) with a global prevalence > 0.1%, three (0.7%) haplogroup associated and 19 (2.0%) variants with a low predicted in silico likelihood of pathogenicity, three variants: MT-TC: m.5772G>A, MT-TF: m.644A>G, and MT-CYB: m.15024G>A, p.C93Y remained. A detailed analysis of these variants indicated that none of them are likely to cause HCM. In conclusion, private mtDNA mutations are frequent, but they are rarely, if ever, associated with HCM.     

Pooled Sequencing of 531 Genes in Inflammatory Bowel Disease Identifies an Associated Rare Variant in BTNL2 and Implicates Other Immune Related Genes

The contribution of rare coding sequence variants to genetic susceptibility in complex disorders is an important but unresolved question. Most studies thus far have investigated a limited number of genes from regions which contain common disease associated variants. Here we investigate this in inflammatory bowel disease by sequencing the exons and proximal promoters of 531 genes selected from both genome-wide association studies and pathway analysis in pooled DNA panels from 474 cases of Crohn’s disease and 480 controls. 80 variants with evidence of association in the sequencing experiment or with potential functional significance were selected for follow up genotyping in 6,507 IBD cases and 3,064 population controls. The top 5 disease associated variants were genotyped in an extension panel of 3,662 IBD cases and 3,639 controls, and tested for association in a combined analysis of 10,147 IBD cases and 7,008 controls. A rare coding variant p.G454C in the BTNL2 gene within the major histocompatibility complex was significantly associated with increased risk for IBD (p = 9.65x10⁻¹⁰, OR = 2.3[95% CI = 1.75–3.04]), but was independent of the known common associated CD and UC variants at this locus. Rare (<1%) and low frequency (1–5%) variants in 3 additional genes showed suggestive association (p<0.005) with either an increased risk (ARIH2 c.338-6C>T) or decreased risk (IL12B p.V298F, and NICN p.H191R) of IBD. These results provide additional insights into the involvement of the inhibition of T cell activation in the development of both sub-phenotypes of inflammatory bowel disease. We suggest that although rare coding variants may make a modest overall contribution to complex disease susceptibility, they can inform our understanding of the molecular pathways that contribute to pathogenesis.     

Association of HK2 and NCK2 with Normal Tension Glaucoma in the Japanese Population

Although family studies and genome-wide association studies have shown that genetic factors play a role in glaucoma, it has been difficult to identify the specific genetic variants involved. We tested 669 single nucleotide polymorphisms (SNPs) from the region of chromosome 2 that includes the GLC1B glaucoma locus for association with primary open-angle glaucoma (POAG) and normal tension glaucoma (NTG) in the Japanese population. We performed a two-stage case-control study. The first cohort consisted of 123 POAG cases, 121 NTG cases and 120 controls: the second cohort consisted of 187 POAG cases, 286 NTG cases, and 271 controls. Out of six SNPs showing significant association with POAG in the first round screening, seven SNPs were tested in the second round. Rs678350 in the HK2 gene coding sequence showed significant allelic (p = 0.0027 in Stage Two, 2.7XE-4 in meta-analysis) association with POAG, and significant allelic (p = 4.7XE-4 in Stage Two, 1.0XE-5 in meta-analysis) association with NTG. Although alleles in the TMEM182 gene did not show significant association with glaucoma in the second round, subjects with the A/A allele in TMEM182 rs869833 showed worse visual field mean deviation (p = 0.01). Even though rs2033008 in the NCK2 gene coding sequence did not show significant association in the first round, it had previously shown association with NTG so it was tested for association with NTG in round 2 (p = 0.0053 in Stage Two). Immunohistochemistry showed that both HK2 and NCK2 are expressed in the retinal ganglion cell layer. Once multi-testing was taken into account, only HK2 showed significant association with POAG and NTG in Stage Two. Our data also support previous reports of NCK2 association with NTG, and raise questions about what role TMEM182 might play in phenotypic variability. Our data suggest that HK2 may play an important role in NTG in the Japanese population.     

Validation of Next-Generation Sequencing of Entire Mitochondrial Genomes and the Diversity of Mitochondrial DNA Mutations in Oral Squamous Cell Carcinoma

Background

Oral squamous cell carcinoma (OSCC) is mainly caused by smoking and alcohol abuse and shows a five-year survival rate of ~50%. We aimed to explore the variation of somatic mitochondrial DNA (mtDNA) mutations in primary oral tumors, recurrences and metastases.

Methods

We performed an in-depth validation of mtDNA next-generation sequencing (NGS) on an Illumina HiSeq 2500 platform for its application to cancer tissues, with the goal to detect low-level heteroplasmies and to avoid artifacts. Therefore we genotyped the mitochondrial genome (16.6 kb) from 85 tissue samples (tumors, recurrences, resection edges, metastases and blood) collected from 28 prospectively recruited OSCC patients applying both Sanger sequencing and high-coverage NGS (~35,000 reads per base).

Results

We observed a strong correlation between Sanger sequencing and NGS in estimating the mixture ratio of heteroplasmies (r = 0.99; p<0.001). Non-synonymous heteroplasmic variants were enriched among cancerous tissues. The proportions of somatic and inherited variants in a given gene region were strongly correlated (r = 0.85; p<0.001). Half of the patients shared mutations between benign and cancerous tissue samples. Low level heteroplasmies (<10%) were more frequent in benign samples compared to tumor samples, where heteroplasmies >10% were predominant. Four out of six patients who developed a local tumor recurrence showed mutations in the recurrence that had also been observed in the primary tumor. Three out of five patients, who had tumor metastases in the lymph nodes of their necks, shared mtDNA mutations between primary tumors and lymph node metastases. The percentage of mutation heteroplasmy increased from the primary tumor to lymph node metastases.

Conclusions

We conclude that Sanger sequencing is valid for heteroplasmy quantification for heteroplasmies ≥10% and that NGS is capable of reliably detecting and quantifying heteroplasmies down to the 1%-level. The finding of shared mutations between primary tumors, recurrences and metastasis indicates a clonal origin of malignant cells in oral cancer.     

Sequence variations of mitochondrial DNA D-loop region are associated with familial nasopharyngeal carcinoma

Mitochondrial DNA (mtDNA) D-loop has been identified as a frequent hot spot of mutations in various tumors. The aim here was to investigate the sequence variations of mitochondrial D-loop region in familial nasopharyngeal carcinoma (FNPC) and their possible associations with cancer risk. 29 subjects from 4 Chinese NPC families and 20 sporadic NPC as well as 122 cases of normal control were recruited. mtDNA extracted from peripheral blood was examined by PCR-based assay for D-loop sequence variations, followed by sequencing analysis. Compared with normal control, four high variations and 6 unrepoted novel polymorphisms were found. Particularly, the np16362 and 16519T to C variants show significantly higher (100%, 81.8%) and lower (0, 22.7%) frequencies in FNPC and unaffected pedigree members, respectively. The occurrence of mitochondrial microsatellite instability (mtMSI) at D310 in experimental groups was statistically significantly higher than in normal control (53.3%). Likewise, in Base Variation Rate consistent with the result, there was a statistically significant difference compared with NC (6.05%). Our results indicated that mtDNA exhibited a high rate of sequence variants in patients with NPC and pedigree members and the mtDNA np16362, np16519 variants and mtMSI at D310 are associated with an increased risk of familial nasopharyngeal carcinoma in pedigree members from families with NPC, which might be involved in the NPC carcinogenesis.     

Recent Mitochondrial DNA Mutations Increase the Risk of Developing Common Late-Onset Human Diseases

Mitochondrial DNA (mtDNA) is highly polymorphic at the population level, and specific mtDNA variants affect mitochondrial function. With emerging evidence that mitochondrial mechanisms are central to common human diseases, it is plausible that mtDNA variants contribute to the “missing heritability” of several complex traits. Given the central role of mtDNA genes in oxidative phosphorylation, the same genetic variants would be expected to alter the risk of developing several different disorders, but this has not been shown to date. Here we studied 38,638 individuals with 11 major diseases, and 17,483 healthy controls. Imputing missing variants from 7,729 complete mitochondrial genomes, we captured 40.41% of European mtDNA variation. We show that mtDNA variants modifying the risk of developing one disease also modify the risk of developing other diseases, thus providing independent replication of a disease association in different case and control cohorts. High-risk alleles were more common than protective alleles, indicating that mtDNA is not at equilibrium in the human population, and that recent mutations interact with nuclear loci to modify the risk of developing multiple common diseases.     

A novel mitochondrial DNA-like sequence in the human nuclear genome.

We describe here a nuclear mitochondrial DNA-like sequence (numtDNA) that is nearly identical in sequence to a continuous 5842 bp segment of human mitochondrial DNA (mtDNA) that spans nucleotide positions 3914 to 9755. On the basis of evolutionary divergence among modern primates, this numtDNA molecule appears to represent mtDNA from a hominid ancestor that has been translocated to the nuclear genome during the recent evolution of humans. This numtDNA sequence harbors synonymous and nonsynonymous nucleotide substitutions relative to the authentic human mtDNA sequence, including an array of substitutions that was previously found in the cytochrome c oxidase subunit 1 and 2 genes. These substitutions were previously reported to occur in human mtDNA, but subsequently contended to be present in a nuclear pseudogene sequence. We now demonstrate their exclusive association with this 5842-bp numtDNA, which we have characterized in its entirety. This numtDNA does not appear to be expressed as a mtDNA-encoded mRNA. It is present in nuclear DNA from human blood donors, in human SH-SY5Y and A431 cell lines, and in rho(0) SH-SY5Y and rho(0) A431 cell lines that were depleted of mtDNA. The existence of human numtDNA sequences with great similarities to human mtDNA renders the amplification of pure mtDNA from cellular DNA very difficult, thereby creating the potential for confounding studies of mitochondrial diseases and population genetics.     

Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

Mitochondrial DNA (mtDNA) variants are widely used in evolutionary genetics as markers for population history and to estimate divergence times among taxa. Inferences of species history are generally based on phylogenetic comparisons, which assume that molecular evolution is clock-like. Between-species comparisons have also been used to estimate the mutation rate, using sites that are thought to evolve neutrally. We directly estimated the mtDNA mutation rate by scanning the mitochondrial genome of Drosophila melanogaster lines that had undergone approximately 200 generations of spontaneous mutation accumulation (MA). We detected a total of 28 point mutations and eight insertion-deletion (indel) mutations, yielding an estimate for the single-nucleotide mutation rate of 6.2 × 10⁻⁸ per site per fly generation. Most mutations were heteroplasmic within a line, and their frequency distribution suggests that the effective number of mitochondrial genomes transmitted per female per generation is about 30. We observed repeated occurrences of some indel mutations, suggesting that indel mutational hotspots are common. Among the point mutations, there is a large excess of G→A mutations on the major strand (the sense strand for the majority of mitochondrial genes). These mutations tend to occur at nonsynonymous sites of protein-coding genes, and they are expected to be deleterious, so do not become fixed between species. The overall mtDNA mutation rate per base pair per fly generation in Drosophila is estimated to be about 10× higher than the nuclear mutation rate, but the mitochondrial major strand G→A mutation rate is about 70× higher than the nuclear rate. Silent sites are substantially more strongly biased towards A and T than nonsynonymous sites, consistent with the extreme mutation bias towards A+T. Strand-asymmetric mutation bias, coupled with selection to maintain specific nonsynonymous bases, therefore provides an explanation for the extreme base composition of the mitochondrial genome of Drosophila.     

Mitochondria DNA Change and Oxidative Damage in Clinically Stable Patients with Major Depressive Disorder

BackgroundTo compare alterations of mitochondria DNA (mtDNA) copy number, single nucleotide polymorphisms (SNPs), and oxidative damage of mtDNA in clinically stable patients with major depressive disorder (MDD).MethodsPatients met DSM-IV diagnostic criteria for MDD were recruited from the psychiatric outpatient clinic at Changhua Christian Hospital, Taiwan. They were clinically stable and their medications had not changed for at least the preceding two months. Exclusion criteria were substance-induced psychotic disorder, eating disorder, anxiety disorder or illicit substance abuse. Comparison subjects did not have any major psychiatric disorder and they were medically healthy. Peripheral blood leukocytes were analyzed to compare copy number, SNPs and oxidative damage of mtDNA between the two groups.Results40 MDD patients and 70 comparison subjects were collected. The median age of the subjects was 42 years and 38 years in MDD and comparison groups, respectively. Leukocyte mtDNA copy number of MDD patients was significantly lower than that of the comparison group (p = 0.037). MDD patients had significantly higher mitochondrial oxidative damage than the comparison group (6.44 vs. 3.90, p<0.001). After generalized linear model adjusted for age, sex, smoking, family history, and psychotropic use, mtDNA copy number was still significantly lower in the MDD group (p<0.001). MtDNA oxidative damage was positively correlated with age (p<0.001) and MDD (p<0.001). Antipsychotic use was negatively associated with mtDNA copy number (p = 0.036).LimitationsThe study is cross-sectional with no longitudinal follow up. The cohort is clinically stable and generalizability of our result to other cohort should be considered.ConclusionsOur study suggests that oxidative stress and mitochondria may play a role in the pathophysiology of MDD. More large-scale studies are warranted to assess the interplay between oxidative stress, mitochondria dysfunction and MDD.     

Sequence Variation in the tRNA Genes of Human Mitochondrial DNA

Recent analyses have shown that nonsynonymous variation in human mitochondrial DNA (mtDNA) contains nonneutral variants, suggesting the presence of mildly deleterious mutations. Many of the disease-causing mutations in mtDNA occur in the genes encoding the tRNAs. Nucleotide sequence variation in these genes has not been studied in human populations, nor have the structural consequences of nucleotide substitutions in tRNA molecules been examined. We therefore determined the nucleotide sequences of the 22 tRNA genes in the mtDNA of 477 Finns and, also, obtained 435 European sequences from the MitoKor database. No differences in population polymorphism indices were found between the two data sets. We assessed selective constraints against various tRNA domains by comparing allele frequencies between these domains and the synonymous and nonsynonymous sites, respectively. All tRNA domains except the variable loop were more conserved than synonymous sites, and T stem and D stem were more conserved than the respective loops. We also analyzed the energetic consequences of the 96 polymorphisms recovered in the two data sets or in the Mitomap database. The minimum free energy (ΔG) was calculated using the free energy rules as implemented in mfold version 3.1. The ΔG’s were normally distributed among the 22 wild-type tRNA genes, whereas the 96 polymorphic tRNAs departed significantly from a normal distribution. The largest differences in ΔG between the wild-type and the polymorphic tRNAs in the Finnish population tended to be in the polymorphisms that were present at low frequencies. Allele frequency distributions and minimum free energy calculations both suggested that some polymorphisms in tRNA genes are nonneutral.Reviewing Editor: Dr. Rüdiger Cerff