Advertisement
Original Research Article| Volume 17, ISSUE 4, P279-284, April 2015

Whole-mitochondrial genome sequencing in primary open-angle glaucoma using massively parallel sequencing identifies novel and known pathogenic variants

      Abstract

      Purpose

      The aim of this study was to determine whether mutations in mitochondrial DNA play a role in high-pressure primary open-angle glaucoma (OMIM 137760) by analyzing new data from massively parallel sequencing of mitochondrial DNA.

      Methods

      Glaucoma patients with high-tension primary open-angle glaucoma and ethnically matched and age-matched control subjects without glaucoma were recruited. The entire human mitochondrial genome was amplified in two overlapping fragments by long-range polymerase chain reaction and used as a template for massively parallel sequencing on an Ion Torrent Personal Genome Machine. All variants were confirmed by conventional Sanger sequencing.

      Results

      Whole-mitochondrial genome sequencing was performed in 32 patients with primary open-angle glaucoma from India (n = 16) and Ireland (n = 16). In 16 of the 32 patients with primary open-angle glaucoma (50% of cases), there were 22 mitochondrial DNA mutations consisting of 7 novel mutations and 8 previously reported disease-associated sequence variants. Eight of 22 (36.4%) of the mitochondrial DNA mutations were in complex I mitochondrial genes.

      Conclusion

      Massively parallel sequencing using the Ion Torrent Personal Genome Machine with confirmation by Sanger sequencing detected a pathogenic mitochondrial DNA mutation in 50% of the primary open-angle glaucoma cohort. Our findings support the emerging concept that mitochondrial dysfunction results in the development of glaucoma and, more specifically, that complex I defects play a significant role in primary open-angle glaucoma pathogenesis.
      Genet Med 17 4, 279–284.

      Keywords

      Introduction

      Glaucoma is a neurodegenerative disease of the eye and is one of the leading causes of visual impairment and blindness, affecting more than 60 million individuals worldwide.
      • Quigley H.A.
      • Broman A.T.
      The number of people with glaucoma worldwide in 2010 and 2020.
      Primary open-angle glaucoma (POAG; OMIM 137760), the most common form of glaucoma, is defined as a spectrum of diseases causing vision loss as a result of progressive and irreversible retinal ganglion cell (RGC) damage, optic nerve degeneration, and corresponding visual field loss.
      • Weinreb R.N.
      • Khaw P.T.
      Primary open-angle glaucoma.
      POAG is usually associated with elevated intraocular pressure (IOP), and this subtype is termed high-tension glaucoma, defined as an IOP >21 mm Hg with a clinically open anterior chamber drainage angle.
      • Weinreb R.N.
      • Khaw P.T.
      Primary open-angle glaucoma.
      POAG accounts for 9–12% of all cases of blindness in the United States.
      • Quigley H.A.
      • Vitale S.
      Models of open-angle glaucoma prevalence and incidence in the United States.
      Glaucoma is characterized by the slow, progressive degeneration of RGCs associated with structural changes at the optic nerve head in the form of pathological cupping and visual field loss.
      • Weinreb R.N.
      • Khaw P.T.
      Primary open-angle glaucoma.
      ,
      • Osborne N.N.
      • del Olmo-Aguado S.
      Maintenance of retinal ganglion cell mitochondrial functions as a neuroprotective strategy in glaucoma.
      The preferential loss of RGCs in glaucoma is also a key pathological feature seen in Leber hereditary optic neuropathy and autosomal dominant optic atrophy.
      • Yu-Wai-Man P.
      • Griffiths P.G.
      • Hudson G.
      • Chinnery P.F.
      Inherited mitochondrial optic neuropathies.
      ,
      • Carelli V.
      • Ross-Cisneros F.N.
      • Sadun A.A.
      Mitochondrial dysfunction as a cause of optic neuropathies.
      Mitochondrial abnormalities are associated with a number of optic neuropathies,
      • Yu-Wai-Man P.
      • Griffiths P.G.
      • Hudson G.
      • Chinnery P.F.
      Inherited mitochondrial optic neuropathies.
      ,
      • Carelli V.
      • Ross-Cisneros F.N.
      • Sadun A.A.
      Mitochondrial dysfunction as a cause of optic neuropathies.
      and there is evidence that mitochondrial DNA (mtDNA) defects and dysfunction play a role in POAG.
      • Lascaratos G.
      • Garway-Heath D.F.
      • Willoughby C.E.
      • Chau K.Y.
      • Schapira A.H.
      Mitochondrial dysfunction in glaucoma: understanding genetic influences.
      ,
      • He Y.
      • Ge J.
      • Tombran-Tink J.
      Mitochondrial defects and dysfunction in calcium regulation in glaucomatous trabecular meshwork cells.
      ,
      • Chrysostomou V.
      • Rezania F.
      • Trounce I.A.
      • Crowston J.G.
      Oxidative stress and mitochondrial dysfunction in glaucoma.
      ,
      • Lee S.
      • Sheck L.
      • Crowston J.G.
      Impaired complex-I-linked respiration and ATP synthesis in primary open-angle glaucoma patient lymphoblasts.
      ,
      • Noh Y.H.
      • Kim K.Y.
      • Shim M.S.
      Inhibition of oxidative stress by coenzyme Q10 increases mitochondrial mass and improves bioenergetic function in optic nerve head astrocytes.
      ,
      • McElnea E.M.
      • Quill B.
      • Docherty N.G.
      Oxidative stress, mitochondrial dysfunction and calcium overload in human lamina cribrosa cells from glaucoma donors.
      Mitochondrial genetic variants have also been linked to other neurodegenerative diseases, including Parkinson disease, Alzheimer disease, and Friedreich ataxia.
      • Chaturvedi R.K.
      • Flint Beal M.
      Mitochondrial diseases of the brain.
      The mitochondrial genome is organized as a circular, double-stranded DNA molecule in the mitochondrion and codes for only 37 genes across ~16,600 base pairs.
      • Druzhyna N.M.
      • Wilson G.L.
      • LeDoux S.P.
      Mitochondrial DNA repair in aging and disease.
      mtDNA can replicate independently of nuclear DNA, and in humans there are usually 100–10,000 separate copies of mtDNA present in each cell. Mitochondrial mutations can be inherited through the maternal line or acquired throughout life (somatic mutations). Somatic mutations accumulate more frequently in the mitochondrial genome, which lacks the protective histones and efficient DNA repair system associated with the nuclear genome.
      • Druzhyna N.M.
      • Wilson G.L.
      • LeDoux S.P.
      Mitochondrial DNA repair in aging and disease.
      ,
      • Larsen N.B.
      • Rasmussen M.
      • Rasmussen L.J.
      Nuclear and mitochondrial DNA repair: similar pathways?.
      Mitochondria may contribute to their own accelerated DNA mutational rate during periods of localized environmental stress; mitochondria both generate and act as a target for excessive reactive oxygen species from altered oxidative phosphorylation.
      • Chrysostomou V.
      • Rezania F.
      • Trounce I.A.
      • Crowston J.G.
      Oxidative stress and mitochondrial dysfunction in glaucoma.
      Several mtDNA variants can exist within a cell or tissue, a situation termed heteroplasmy. New mtDNA mutations can coexist at a low frequency with wild-type mtDNAs, making it challenging to identify all mutations in the mitochondrial genome. “Massively parallel” or “next-generation” sequencing approaches have ushered in a new era of opportunity in mtDNA mutation analysis,
      • Tang S.
      • Huang T.
      Characterization of mitochondrial DNA heteroplasmy using a parallel sequencing system.
      ,
      • Li M.
      • Schönberg A.
      • Schaefer M.
      • Schroeder R.
      • Nasidze I.
      • Stoneking M.
      Detecting heteroplasmy from high-throughput sequencing of complete human mitochondrial DNA genomes.
      enhancing understanding of mtDNA-based disorders.
      • Zaragoza M.V.
      • Fass J.
      • Diegoli M.
      • Lin D.
      • Arbustini E.
      Mitochondrial DNA variant discovery and evaluation in human cardiomyopathies through next-generation sequencing.
      Massively parallel sequencing can characterize the complete nucleotide sequence for the mitochondrial genome and provide adequate sequencing depth to identify heteroplasmic positions with enhanced sensitivity and estimate the mutant load at each nucleotide position. The aim of this study was to determine whether mutations in mtDNA play a role in POAG using massively parallel, whole-mitochondrial genome sequencing on the Ion Torrent Personal Genome Machine (PGM; Life Technologies, Paisley, UK).

      Materials and Methods

      Patients

      Patients with high-tension POAG were recruited after careful phenotyping from the departments of ophthalmology in Dublin (Mater Misericordiae University Hospital, Dublin, Ireland) and Pune (PBMA’s H. V. Desai Eye Hospital, Pune, India). Written informed consent was obtained from all subjects, and the study was performed according to the tenets of the Declaration of Helsinki and the guidelines of the Indian Council of Medical Research, New Delhi, India. All subjects underwent a standardized ophthalmic examination; the clinical parameters used for the case definition of high-tension POAG were as follows: (i) presenting IOP >21 mm Hg in at least one eye by Goldmann applanation tonometry; (ii) glaucomatous optic nerve head damage on stereoscopic optic disc examination (notching or thinning of the neuroretinal rim and/or increased cup/disc ratio in relation to the optic disc size); (iii) an open anterior chamber angle on gonioscopy; (iv) reproducible and characteristic glaucomatous visual field defect with the Humphrey 24-2 full-threshold strategy when available (Carl Zeiss Meditec, Oberkochen, Germany); and (v) absence of conditions leading to secondary glaucoma (e.g., pigmentary dispersion syndrome, pseudoexfoliation, history of steroid use, uveitis, or ocular trauma). Ethnically matched and age-matched control subjects without glaucoma (based on normal IOP) were also recruited from each study population.

      DNA Extraction and Polymerase Chain Reaction Amplification of Mitochondrial Genome

      Genomic DNA was extracted from peripheral blood leukocytes from all the subjects using a Qiagen DNA Purification Kit (Qiagen, Manchester, UK) according to the manufacturer’s instructions. The entire human mitochondrial genome was amplified in two overlapping fragments: fragment 1 (5′-AACCAAACCCCAAAGACACC-3′ and 5′-GCCAATAATGACGTGAAGTCC-3′; product size: 9,289 bp) and fragment 2 (5′-TCCCACTCCTAAACACATCC-3′ and 5′-TTTATGGGGTGATGTGAGCC-3′; product size: 7,626 bp).
      • Tang S.
      • Huang T.
      Characterization of mitochondrial DNA heteroplasmy using a parallel sequencing system.
      Polymerase chain reaction (PCR) amplification was performed using a TaKaRa LA PCR Amplification Kit (Takara Bio, Shiga, Japan) in a 25-μl volume containing 0.125 μl (5 units/μl) of Taq DNA polymerase, 2.5 μl of 10× Ex Taq buffer, 2 μl deoxynucleotide triphosphates (dNTPs) (2.5 mmol/l each), and 1 μmol/l primers. Thermal cycling was performed in a Verti thermal cycler (Life Technologies) with the following conditions: initial denaturation at 94 °C for 1 min, followed by 35 cycles at 94 °C for 30 s, 64 °C for 1 min, and 72 °C for 8 min, and a final extension at 72 °C for 5 min. The amplified PCR products were purified using AMPure beads (Agencourt AMPure XP; Beckman Coulter, Danvers, MA) and quantified using a Nanodrop-1000 spectrophotometer (Nanodrop ND-1000; Thermo Scientific, Waltham, MA).

      Massively Parallel Sequencing on the Ion Torrent PGM

      Massively parallel sequencing was performed on the Ion Torrent PGM according to the manufacturer’s protocol (Life Technologies) and included library preparation, template preparation using emulsion PCR, and PGM sequencing. Specifically, Ion Torrent sequencing libraries were prepared using the Ion Xpress Fragment Plus Library Kit (Life Technologies) according to the manufacturer’s protocol, with purification between steps performed using Agencourt AMPure XP beads (Beckman Coulter). First, equimolar amounts of mtDNA PCR fragments 1 and 2 were pooled for each sample and enzymatically sheared by incubation at 37 °C for 15 min with Ion Shear enzyme (Life Technologies). After adaptor ligation (using Ion Xpress barcodes to uniquely identify each sample), the fragments were size-selected (150–250 bp) using an E-Gel SizeSelect 2% agarose gel (Life Technologies). The library was then amplified and quantified using a Qubit Fluorometer (Life Technologies). The library size of ~200 bp was verified using a 2100 Bioanalyzer (Agilent Technologies UK, Wokingham, UK). Clonal amplification of libraries on Ion Sphere particles was achieved by emulsion PCR performed on a One Touch system using the Ion OneTouch 100-bp Template Kit according to the manufacturer’s instructions (Life Technologies). Template-positive Ion Sphere particles were isolated using the Ion OneTouch ES enrichment system. Massively parallel sequencing was performed on an Ion Torrent PGM instrument.

      PGM Data Analysis

      FASTQ output from the Ion Torrent was uploaded to Partek Flow for initial quality control and conversion to BAM files for analysis using Partek Genomics Suite version 6.6 (Partek, St. Louis, MO). The overall quality of reads was assessed using pre-alignment quality control with 10% of reads sampled per individual. Reads with less than 25 bases were removed and reads were trimmed until the quality score exceeded 20. Trimmed data were independently aligned using Bowtie version 0.12.7 (based on default settings optimized by Partek) to NCBI Reference mtDNA reference sequence (revised Cambridge Reference Sequence) NC_012920.1, GI:251831106 (hg19) with 0.2-bit reference built from FASTA sequence. mtDNA sequences for each sample were compared with the reference using single-nucleotide polymorphism (SNP) and allele percentage options.

      Sanger Sequencing of Mitochondrial Genome

      To corroborate the Ion Torrent PGM massively parallel sequencing, all patient mtDNA sequences and mtDNA from ethnically matched controls underwent conventional Sanger sequencing. Nine overlapping primer pairs were used for mtDNA amplification with 62 corresponding overlapping internal sequencing primers for mtDNA sequencing as previously described.
      • Ramos A.
      • Santos C.
      • Alvarez L.
      • Nogués R.
      • Aluja M.P.
      Human mitochondrial DNA complete amplification and sequencing: a new validated primer set that prevents nuclear DNA sequences of mitochondrial origin co-amplification.
      ,
      • Ramos A.
      • Santos C.
      • Barbena E.
      Validated primer set that prevents nuclear DNA sequences of mitochondrial origin co-amplification: a revision based on the New Human Genome Reference Sequence (GRCh37).
      Bidirectional cycle sequencing was performed using BigDye Terminator v3.1 cycle sequencing kit (Life Technologies), and amplified, targeted genetic regions electrophoresed on an ABI PRISM 3730 Genetic Analyzer (Life Technologies).

      mtDNA Variant Analysis

      Identified mtDNA sequence variants were compared with publicly accessible databases: MITOMAP (http://www.mitomap.org/MITOMAP), mtDB (http://www.mtdb.igp.uu.se/), dbSNP (Build 137) (http://www.ncbi.nlm.nih.gov/SNP/), and the 1000 Genomes Project (http://browser.1000genomes.org/index.html). Previously reported disease-associated sequence variants were identified using MITOMAP. Previously unreported or novel nucleotide changes were classified as pathogenic if they (i) were absent from ethnically matched controls (110 individuals without glaucoma), (ii) were not reported in the 24,187 human mitochondrial DNA sequences in MITOMAP, and (iii) were classified as possibly or probably damaging when considering the conservation and deleterious structural effects of amino acid substitutions on protein function using PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2/). Identified sequence variants were described according to the guidelines published by the Human Genome Variation Society (http://www.hgvs.org). Transitions were represented by purine–purine (A<->G) or pyrimidine–pyrimidine (C<->T) substitutions, and transversions were represented by purine–pyrimidine substitutions (A<->C, A<->T, G<->C, G<->T).

      Results

      Whole-mitochondrial genome sequencing was performed in 32 patients with high-tension POAG from India (n = 16) and Ireland (n = 16) on the Ion Torrent PGM. Pre-alignment quality assurance and quality control on raw reads from the PGM (FASTQ files) demonstrated an average read length of 98 bases, with average quality score of 24 and 45% guanine cytosine content across the sequence reads. FASTA-trimmed read alignments demonstrated that more than 99.98% of mtDNA bases were covered from an average read length of 93 bases, with an average depth of ×267. Unique sequence variants (n = 536) were identified in the Ion Torrent PGM sequencing data using the proportional allele model. Specifically, 334 Ion Torrent PGM sequence variants were confirmed by Sanger sequencing, whereas 57 were not confirmed by direct Sanger sequencing (likely false positives or below the heteroplasmy detection limit for Sanger sequencing), whereas the remainder were not able to be called by Sanger sequencing because of poor coverage, indels, or complex variations. Four hundred four unique SNPs were identified with log odds ratio >5.0 for the SNP-against-reference option in Partek. Eleven sequence variants were detected only by Sanger sequencing and were not identified by the algorithms used for analyzing the Ion Torrent PGM sequencing data (2.7% false-negative rate). Sequenced variants detected in ethnically matched controls (n = 110) were excluded from further analysis.
      Sixteen of the 32 POAG patients (50% of cases) carried mtDNA sequence variants. We found a total of 22 mtDNA sequence variants, consisting of seven novel mutations and eight previously reported disease-associated sequence variants (Table 1); six patients (18.75%) harbored two pathogenic mtDNA mutations. Eight of 22 (36.4%) mtDNA mutations were in complex I mitochondrial genes, and five of 22 (22.7%) were in transfer RNAs. Among the POAG patients with mutations, five of 16 (31.25%) had a mutation in the MT-ND5 gene in complex I. The eight previously reported disease-associated mtDNA sequence variants (Supplementary Table S1 online) were reported in MITOMAP. There was no evidence of a genotype–phenotype correlation between the number or location of mtDNA sequence variants and the severity of the POAG phenotype. None of the detected mtDNA sequence variants has previously been reported in POAG. Ten novel synonymous mtDNA variants were detected (Supplementary Table S2 online). In the Caucasian Irish POAG patients, nine of 16 (56.25%) had at least one pathogenic mtDNA variant, and in the Indian POAG patients, seven of 16 (43.75%) had at least one pathogenic mtDNA variant.
      Table 1Novel and known pathogenic mtDNA variants identified in Irish and Indian POAG cases using massively parallel sequencing

      Discussion

      With a prevalence of 1–3% in Europe and up to 9% in Africa, glaucoma is a major cause of visual morbidity and the second most common cause of blindness in the world, resulting in bilateral blindness in 4.5 million people.
      • Quigley H.A.
      • Broman A.T.
      The number of people with glaucoma worldwide in 2010 and 2020.
      ,
      • Leske M.C.
      Open-angle glaucoma – an epidemiologic overview.
      Visual loss with POAG is the result of progressive and irreversible RGC damage.
      • Weinreb R.N.
      • Khaw P.T.
      Primary open-angle glaucoma.
      Given the high metabolic demands of RGCs, they are heavily dependent on mitochondrial function for their function and survival.
      • Chrysostomou V.
      • Rezania F.
      • Trounce I.A.
      • Crowston J.G.
      Oxidative stress and mitochondrial dysfunction in glaucoma.
      ,
      • Lee S.
      • Van Bergen N.J.
      • Kong G.Y.
      Mitochondrial dysfunction in glaucoma and emerging bioenergetic therapies.
      Mitochondrial dysfunction and mtDNA defects have been implicated in RGC loss and the development of glaucoma.
      • Lascaratos G.
      • Garway-Heath D.F.
      • Willoughby C.E.
      • Chau K.Y.
      • Schapira A.H.
      Mitochondrial dysfunction in glaucoma: understanding genetic influences.
      ,
      • Lee S.
      • Sheck L.
      • Crowston J.G.
      Impaired complex-I-linked respiration and ATP synthesis in primary open-angle glaucoma patient lymphoblasts.
      ,
      • Lee S.
      • Van Bergen N.J.
      • Kong G.Y.
      Mitochondrial dysfunction in glaucoma and emerging bioenergetic therapies.
      ,
      • Abu-Amero K.K.
      • Morales J.
      • Bosley T.M.
      Mitochondrial abnormalities in patients with primary open-angle glaucoma.
      The expansion of massively parallel sequencing has revolutionized mitochondrial genetics.
      • Calvo S.E.
      • Compton A.G.
      • Hershman S.G.
      Molecular diagnosis of infantile mitochondrial disease with targeted next-generation sequencing.
      ,
      • Vasta V.
      • Ng S.B.
      • Turner E.H.
      • Shendure J.
      • Hahn S.H.
      Next generation sequence analysis for mitochondrial disorders.
      ,
      • Cui H.
      • Li F.
      • Chen D.
      Comprehensive next-generation sequence analyses of the entire mitochondrial genome reveal new insights into the molecular diagnosis of mitochondrial DNA disorders.
      Sanger sequencing has long been considered the gold standard for molecular diagnosis of mtDNA-based disorders. The vast majority of deleterious mtDNA mutations are heteroplasmic, and it is challenging to identify all mutations in the mitochondrial genome and simultaneously quantify the mtDNA heteroplasmy levels.
      • Tang S.
      • Huang T.
      Characterization of mitochondrial DNA heteroplasmy using a parallel sequencing system.
      Previously, the two most popular complete mitochondrial genome sequencing methods were Sanger sequencing
      • Abu-Amero K.K.
      • Morales J.
      • Bosley T.M.
      Mitochondrial abnormalities in patients with primary open-angle glaucoma.
      and the Affymetrix MitoChip,
      • van Eijsden R.G.
      • Gerards M.
      • Eijssen L.M.
      Chip-based mtDNA mutation screening enables fast and reliable genetic diagnosis of OXPHOS patients.
      but neither is sufficiently sensitive or specific to reliably detect mtDNA heteroplasmy. In comparison, massively parallel sequencing has been reported to detect mtDNA heteroplasmies of ≥5% with virtually no false-positive results.
      • Tang S.
      • Huang T.
      Characterization of mitochondrial DNA heteroplasmy using a parallel sequencing system.
      ,
      • Cui H.
      • Li F.
      • Chen D.
      Comprehensive next-generation sequence analyses of the entire mitochondrial genome reveal new insights into the molecular diagnosis of mitochondrial DNA disorders.
      ,
      • Tang S.
      • Wang J.
      • Zhang V.W.
      Transition to next generation analysis of the whole mitochondrial genome: a summary of molecular defects.
      The application of massively parallel sequencing technology to the analysis of the mitochondrial genome has significantly improved the molecular diagnosis of mtDNA disorders.
      • Calvo S.E.
      • Compton A.G.
      • Hershman S.G.
      Molecular diagnosis of infantile mitochondrial disease with targeted next-generation sequencing.
      Large mtDNA deletions have been reported in mitochondrial diseases, and the ability of massively parallel sequencing to detect mtDNA deletions has not been firmly established.
      • Kingsley C.B.
      Identification of causal sequence variants of disease in the next generation sequencing era.
      However, using an long range PCR strategy similar to that in our study, the Ion Torrent PGM was able to detect two large-scale single deletions in the mtDNA of two patients.

      Seneca S, Vancampenhout K, Van Coster R, et al. Analysis of the whole mitochondrial genome: translation of the Ion Torrent Personal Genome Machine system to the diagnostic bench? Eur J Hum Genet 2014; e-pub ahead of print 26 March 2014.

      Two other studies have used massively parallel sequencing for whole-mitochondrial sequencing in POAG. Collins et al.
      • Collins D.W.
      • Gudiseva H.V.
      • Trachtman B.T.
      Mitochondrial sequence variation in African-American primary open-angle glaucoma patients.
      used the Ion Torrent PGM to sequence the mitochondrial genome in 22 Afro-Caribbean POAG cases and 22 controls; they concluded that African-American POAG cases closely resembled controls with respect to frequency and types of mtDNA sequence variation. Using next-generation sequencing, Jeoung et al.
      • Jeoung J.W.
      • Seong M.W.
      • Park S.S.
      • Kim D.M.
      • Kim S.H.
      • Park K.H.
      Mitochondrial DNA variant discovery in normal-tension glaucoma patients by next-generation sequencing.
      identified three mtDNA SNPs (m.4883C>T ND2, m.9540T>C COX3, and m.14766C>T CYTB) with significantly different frequencies in Korean patients with normal-tension glaucoma as compared with controls. In the present study, we evaluated 32 POAG patients using massively parallel sequencing on the Ion Torrent PGM to investigate the involvement of mitochondrial genomic variants in POAG pathogenesis and detected pathogenic mtDNA mutations in 50% of patients. The mutational load was consistent between the Irish Caucasian and Indian POAG patients, with complex I mitochondrial genes being the most commonly mutated, representing 36.4% of mutations. Our study is the first to assess mtDNA defects in Caucasian subjects with POAG. This mutational load is consistent with findings from Sanger sequencing of mtDNA in Middle Eastern Arabs with POAG that identified pathogenic mtDNA mutations in 54% of patients with POAG (14/27).
      • Abu-Amero K.K.
      • Morales J.
      • Bosley T.M.
      Mitochondrial abnormalities in patients with primary open-angle glaucoma.
      Mitochondrial mutations have also been reported in other forms of glaucoma at low frequencies: primary congenital glaucoma (10–22.85%),
      • Tanwar M.
      • Dada T.
      • Sihota R.
      • Dada R.
      Mitochondrial DNA analysis in primary congenital glaucoma.
      ,
      • Kumar M.
      • Tanwar M.
      • Faiq M.A.
      Mitochondrial DNA nucleotide changes in primary congenital glaucoma patients.
      primary angle closure glaucoma (13.8%),
      • Abu-Amero K.K.
      • Morales J.
      • Osman M.N.
      • Bosley T.M.
      Nuclear and mitochondrial analysis of patients with primary angle-closure glaucoma.
      and pseudoexfoliation glaucoma (10.3%).
      • Abu-Amero K.K.
      • Bosley T.M.
      • Morales J.
      Analysis of nuclear and mitochondrial genes in patients with pseudoexfoliation glaucoma.
      All previous mtDNA sequencing studies, including our study, were performed on mtDNA isolated from peripheral blood leukocytes, and it is known that mtDNA sequence variants may exist at different levels in different tissues (tissue heteroplasmy). Without studying mtDNA defects, including functional analyses, in the target tissue within the eye, these findings cannot be extrapolated to the site of disease pathology; however, they do provide evidence that mtDNA defects are frequently present in leukocytes of POAG patients, indicating that mitochondrial dysfunction is related to disease status.
      Defects in complex I oxidative phosphorylation resulting in decreased respiration rates and adenosine triphosphate production have been reported in POAG lymphocytes, although this cohort was not screened for mtDNA defects.
      • Lee S.
      • Sheck L.
      • Crowston J.G.
      Impaired complex-I-linked respiration and ATP synthesis in primary open-angle glaucoma patient lymphoblasts.
      In our study, complex I was mutated in approximately one-third (31.25%) of the POAG patients with mtDNA mutations and, specifically, mutations in the MT-ND5 gene. Banerjee et al.
      • Banerjee D.
      • Banerjee A.
      • Mookherjee S.
      Mitochondrial genome analysis of primary open angle glaucoma patients.
      also observed a clustering of nonsynonymous mtDNA variants in MT-ND5 in a cohort of 101 patients with POAG from India evaluated by Sanger sequencing of the mitochondrial genome. Mitochondrial complex I is the largest component of the oxidative phosphorylation system and is a crucial component of the mitochondrial respiratory chain.
      • Papa S.
      • De Rasmo D.
      Complex I deficiencies in neurological disorders.
      ,
      • Roestenberg P.
      • Manjeri G.R.
      • Valsecchi F.
      • Smeitink J.A.
      • Willems P.H.
      • Koopman W.J.
      Pharmacological targeting of mitochondrial complex I deficiency: the cellular level and beyond.
      The function of complex I is essential for the maintenance of tissue reactive oxygen species levels, generation of mitochondrial membrane potential, and adenosine triphosphate production.
      • Roestenberg P.
      • Manjeri G.R.
      • Valsecchi F.
      • Smeitink J.A.
      • Willems P.H.
      • Koopman W.J.
      Pharmacological targeting of mitochondrial complex I deficiency: the cellular level and beyond.
      Complex I dysfunction leads to alterations in cell survival, apoptosis, and calcium homeostasis, and is the most frequent respiratory chain defect found in neuromuscular mitochondrial disorders.
      • Papa S.
      • De Rasmo D.
      Complex I deficiencies in neurological disorders.
      Genetic dysfunction of complex I is the most prevalent cause of Leber hereditary optic neuropathy, an inherited mitochondrial optic neuropathy.
      • Yu-Wai-Man P.
      • Griffiths P.G.
      • Chinnery P.F.
      Mitochondrial optic neuropathies - disease mechanisms and therapeutic strategies.
      Our study supports the emerging concept that mitochondrial dysfunction results in the development of glaucoma and, more specifically, that complex I and MT-ND5 defects play a significant role in POAG pathogenesis.
      • Lascaratos G.
      • Garway-Heath D.F.
      • Willoughby C.E.
      • Chau K.Y.
      • Schapira A.H.
      Mitochondrial dysfunction in glaucoma: understanding genetic influences.
      ,
      • Lee S.
      • Van Bergen N.J.
      • Kong G.Y.
      Mitochondrial dysfunction in glaucoma and emerging bioenergetic therapies.
      ,
      • Banerjee D.
      • Banerjee A.
      • Mookherjee S.
      Mitochondrial genome analysis of primary open angle glaucoma patients.
      Further studies using larger cohorts of glaucoma patients and phenotypes are required to determine the relationship between the mutational load in the mitochondrial genome and glaucoma pathogenesis. Massively parallel sequencing is a cost-effective approach with the high sensitivity, specificity, and accuracy necessary to detect a wide spectrum of mtDNA mutations while also quantifying levels of heteroplasmy.
      • Tang S.
      • Huang T.
      Characterization of mitochondrial DNA heteroplasmy using a parallel sequencing system.
      ,
      • Calvo S.E.
      • Compton A.G.
      • Hershman S.G.
      Molecular diagnosis of infantile mitochondrial disease with targeted next-generation sequencing.
      ,
      • Cui H.
      • Li F.
      • Chen D.
      Comprehensive next-generation sequence analyses of the entire mitochondrial genome reveal new insights into the molecular diagnosis of mitochondrial DNA disorders.
      ,
      • Tang S.
      • Wang J.
      • Zhang V.W.
      Transition to next generation analysis of the whole mitochondrial genome: a summary of molecular defects.
      Without functional correlates, distinguishing mildly pathogenic mtDNA mutations from benign variants is challenging and can depend on interpretation. Future studies should correlate the mtDNA mutational spectrum with functional assays to more fully explore the role mitochondrial mutations and dysfunction play in the pathogenesis of glaucoma.

      Disclosure

      The authors declare no conflict of interest.

      Acknowledgements

      This study was partially funded by the British Council for the Prevention of Blindness. C.S. was the recipient of the Sir John Wilson Prevention of Blindness Fellowship. P.S. received funding from the Department of Biotechnology, Government of India (CREST award 2010–2011). We thank the subjects who took part in this research.

      Author Contributions

      Periasamy Sundaresan and David A. Simpson: The first two authors contributed equally to this work.

      Supplementary information

      References

        • Quigley H.A.
        • Broman A.T.
        The number of people with glaucoma worldwide in 2010 and 2020.
        Br J Ophthalmol. 2006; 90: 262-267
        • Weinreb R.N.
        • Khaw P.T.
        Primary open-angle glaucoma.
        Lancet. 2004; 363: 1711-1720
        • Quigley H.A.
        • Vitale S.
        Models of open-angle glaucoma prevalence and incidence in the United States.
        Invest Ophthalmol Vis Sci. 1997; 38: 83-91
        • Osborne N.N.
        • del Olmo-Aguado S.
        Maintenance of retinal ganglion cell mitochondrial functions as a neuroprotective strategy in glaucoma.
        Curr Opin Pharmacol. 2013; 13: 16-22
        • Yu-Wai-Man P.
        • Griffiths P.G.
        • Hudson G.
        • Chinnery P.F.
        Inherited mitochondrial optic neuropathies.
        J Med Genet. 2009; 46: 145-158
        • Carelli V.
        • Ross-Cisneros F.N.
        • Sadun A.A.
        Mitochondrial dysfunction as a cause of optic neuropathies.
        Prog Retin Eye Res. 2004; 23: 53-89
        • Lascaratos G.
        • Garway-Heath D.F.
        • Willoughby C.E.
        • Chau K.Y.
        • Schapira A.H.
        Mitochondrial dysfunction in glaucoma: understanding genetic influences.
        Mitochondrion. 2012; 12: 202-212
        • He Y.
        • Ge J.
        • Tombran-Tink J.
        Mitochondrial defects and dysfunction in calcium regulation in glaucomatous trabecular meshwork cells.
        Invest Ophthalmol Vis Sci. 2008; 49: 4912-4922
        • Chrysostomou V.
        • Rezania F.
        • Trounce I.A.
        • Crowston J.G.
        Oxidative stress and mitochondrial dysfunction in glaucoma.
        Curr Opin Pharmacol. 2013; 13: 12-15
        • Lee S.
        • Sheck L.
        • Crowston J.G.
        Impaired complex-I-linked respiration and ATP synthesis in primary open-angle glaucoma patient lymphoblasts.
        Invest Ophthalmol Vis Sci. 2012; 53: 2431-2437
        • Noh Y.H.
        • Kim K.Y.
        • Shim M.S.
        Inhibition of oxidative stress by coenzyme Q10 increases mitochondrial mass and improves bioenergetic function in optic nerve head astrocytes.
        Cell Death Dis. 2013; 4: e820
        • McElnea E.M.
        • Quill B.
        • Docherty N.G.
        Oxidative stress, mitochondrial dysfunction and calcium overload in human lamina cribrosa cells from glaucoma donors.
        Mol Vis. 2011; 17: 1182-1191
        • Chaturvedi R.K.
        • Flint Beal M.
        Mitochondrial diseases of the brain.
        Free Radic Biol Med. 2013; 63: 1-29
        • Druzhyna N.M.
        • Wilson G.L.
        • LeDoux S.P.
        Mitochondrial DNA repair in aging and disease.
        Mech Ageing Dev. 2008; 129: 383-390
        • Larsen N.B.
        • Rasmussen M.
        • Rasmussen L.J.
        Nuclear and mitochondrial DNA repair: similar pathways?.
        Mitochondrion. 2005; 5: 89-108
        • Tang S.
        • Huang T.
        Characterization of mitochondrial DNA heteroplasmy using a parallel sequencing system.
        BioTechniques. 2010; 48: 287-296
        • Li M.
        • Schönberg A.
        • Schaefer M.
        • Schroeder R.
        • Nasidze I.
        • Stoneking M.
        Detecting heteroplasmy from high-throughput sequencing of complete human mitochondrial DNA genomes.
        Am J Hum Genet. 2010; 87: 237-249
        • Zaragoza M.V.
        • Fass J.
        • Diegoli M.
        • Lin D.
        • Arbustini E.
        Mitochondrial DNA variant discovery and evaluation in human cardiomyopathies through next-generation sequencing.
        PLoS ONE. 2010; 5: e12295
        • Ramos A.
        • Santos C.
        • Alvarez L.
        • Nogués R.
        • Aluja M.P.
        Human mitochondrial DNA complete amplification and sequencing: a new validated primer set that prevents nuclear DNA sequences of mitochondrial origin co-amplification.
        Electrophoresis. 2009; 30: 1587-1593
        • Ramos A.
        • Santos C.
        • Barbena E.
        Validated primer set that prevents nuclear DNA sequences of mitochondrial origin co-amplification: a revision based on the New Human Genome Reference Sequence (GRCh37).
        Electrophoresis. 2011; 32: 782-783
        • Leske M.C.
        Open-angle glaucoma – an epidemiologic overview.
        Ophthalmic Epidemiol. 2007; 14: 166-172
        • Lee S.
        • Van Bergen N.J.
        • Kong G.Y.
        Mitochondrial dysfunction in glaucoma and emerging bioenergetic therapies.
        Exp Eye Res. 2011; 93: 204-212
        • Abu-Amero K.K.
        • Morales J.
        • Bosley T.M.
        Mitochondrial abnormalities in patients with primary open-angle glaucoma.
        Invest Ophthalmol Vis Sci. 2006; 47: 2533-2541
        • Calvo S.E.
        • Compton A.G.
        • Hershman S.G.
        Molecular diagnosis of infantile mitochondrial disease with targeted next-generation sequencing.
        Sci Transl Med. 2012; 4: 118ra110
        • Vasta V.
        • Ng S.B.
        • Turner E.H.
        • Shendure J.
        • Hahn S.H.
        Next generation sequence analysis for mitochondrial disorders.
        Genome Med. 2009; 1: 100
        • Cui H.
        • Li F.
        • Chen D.
        Comprehensive next-generation sequence analyses of the entire mitochondrial genome reveal new insights into the molecular diagnosis of mitochondrial DNA disorders.
        Genet Med. 2013; 15: 388-394
        • van Eijsden R.G.
        • Gerards M.
        • Eijssen L.M.
        Chip-based mtDNA mutation screening enables fast and reliable genetic diagnosis of OXPHOS patients.
        Genet Med. 2006; 8: 620-627
        • Tang S.
        • Wang J.
        • Zhang V.W.
        Transition to next generation analysis of the whole mitochondrial genome: a summary of molecular defects.
        Hum Mutat. 2013; 34: 882-893
        • Kingsley C.B.
        Identification of causal sequence variants of disease in the next generation sequencing era.
        Methods Mol Biol. 2011; 700: 37-46
      1. Seneca S, Vancampenhout K, Van Coster R, et al. Analysis of the whole mitochondrial genome: translation of the Ion Torrent Personal Genome Machine system to the diagnostic bench? Eur J Hum Genet 2014; e-pub ahead of print 26 March 2014.

        • Collins D.W.
        • Gudiseva H.V.
        • Trachtman B.T.
        Mitochondrial sequence variation in African-American primary open-angle glaucoma patients.
        PLoS ONE. 2013; 8: e76627
        • Jeoung J.W.
        • Seong M.W.
        • Park S.S.
        • Kim D.M.
        • Kim S.H.
        • Park K.H.
        Mitochondrial DNA variant discovery in normal-tension glaucoma patients by next-generation sequencing.
        Invest Ophthalmol Vis Sci. 2014; 55: 986-992
        • Tanwar M.
        • Dada T.
        • Sihota R.
        • Dada R.
        Mitochondrial DNA analysis in primary congenital glaucoma.
        Mol Vis. 2010; 16: 518-533
        • Kumar M.
        • Tanwar M.
        • Faiq M.A.
        Mitochondrial DNA nucleotide changes in primary congenital glaucoma patients.
        Mol Vis. 2013; 19: 220-230
        • Abu-Amero K.K.
        • Morales J.
        • Osman M.N.
        • Bosley T.M.
        Nuclear and mitochondrial analysis of patients with primary angle-closure glaucoma.
        Invest Ophthalmol Vis Sci. 2007; 48: 5591-5596
        • Abu-Amero K.K.
        • Bosley T.M.
        • Morales J.
        Analysis of nuclear and mitochondrial genes in patients with pseudoexfoliation glaucoma.
        Mol Vis. 2008; 14: 29-36
        • Banerjee D.
        • Banerjee A.
        • Mookherjee S.
        Mitochondrial genome analysis of primary open angle glaucoma patients.
        PLoS ONE. 2013; 8: e70760
        • Papa S.
        • De Rasmo D.
        Complex I deficiencies in neurological disorders.
        Trends Mol Med. 2013; 19: 61-69
        • Roestenberg P.
        • Manjeri G.R.
        • Valsecchi F.
        • Smeitink J.A.
        • Willems P.H.
        • Koopman W.J.
        Pharmacological targeting of mitochondrial complex I deficiency: the cellular level and beyond.
        Mitochondrion. 2012; 12: 57-65
        • Yu-Wai-Man P.
        • Griffiths P.G.
        • Chinnery P.F.
        Mitochondrial optic neuropathies - disease mechanisms and therapeutic strategies.
        Prog Retin Eye Res. 2011; 30: 81-114