Whole exome sequencing implicates eye development, the unfolded protein response and plasma membrane homeostasis in primary open-angle glaucoma

Purpose To identify biological processes associated with POAG and its subtypes, high-tension (HTG) and normal-tension glaucoma (NTG), by analyzing rare potentially damaging genetic variants. Methods A total of 122 and 65 unrelated HTG and NTG participants, respectively, with early onset advanced POAG, 103 non-glaucoma controls and 993 unscreened ethnicity-matched controls were included in this study. Study participants without myocilin disease-causing variants and non-glaucoma controls were subjected to whole exome sequencing on an Illumina HiSeq2000. Exomes of participants were sequenced on an Illumina HiSeq2000. Qualifying variants were rare in the general population (MAF < 0.001) and potentially functionally damaging (nonsense, frameshift, splice or predicted pathogenic using SIFT or Polyphen2 software). Genes showing enrichment of qualifying variants in cases were selected for pathway and network analysis using InnateDB. Results POAG cases showed enrichment of rare variants in camera-type eye development genes (p = 1.40×10–7, corrected p = 3.28×10–4). Implicated eye development genes were related to neuronal or retinal development. HTG cases were significantly enriched for key regulators in the unfolded protein response (UPR) (p = 7.72×10–5, corrected p = 0.013). The UPR is known to be involved in myocilin-related glaucoma; our results suggest the UPR has a role in non-myocilin causes of HTG. NTG cases showed enrichment in ion channel transport processes (p = 1.05×10–4, corrected p = 0.027) including calcium, chloride and phospholipid transporters involved in plasma membrane homeostasis. Network analysis also revealed enrichment of the MHC Class I antigen presentation pathway in HTG, and the EGFR1 and cell-cycle pathways in both HTG and NTG. Conclusion This study suggests that mutations in eye development genes are enriched in POAG. HTG can result from aberrant responses to protein misfolding which may be amenable to molecular chaperone therapy. NTG is associated with impaired plasma membrane homeostasis increasing susceptibility to apoptosis.

[1]  Deborah L. Johnson,et al.  Epidermal Growth Factor Receptor 1 (EGFR1) and Its Variant EGFRvIII Regulate TATA-Binding Protein Expression through Distinct Pathways , 2008, Molecular and Cellular Biology.

[2]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[3]  Binbin Wang,et al.  A novel mutation in GJA8 causing congenital cataract–microcornea syndrome in a Chinese pedigree , 2010, Molecular vision.

[4]  A. Izzotti,et al.  The Dysfunction of the Trabecular Meshwork During Glaucoma Course , 2015, Journal of cellular physiology.

[5]  R. Ritch,et al.  Adult-Onset Primary Open-Angle Glaucoma Caused by Mutations in Optineurin , 2002, Science.

[6]  B. Morquette,et al.  The molecular basis of retinal ganglion cell death in glaucoma , 2012, Progress in Retinal and Eye Research.

[7]  S. Henikoff,et al.  Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm , 2009, Nature Protocols.

[8]  J. Roth,et al.  Sodium 4-phenylbutyrate acts as a chemical chaperone on misfolded myocilin to rescue cells from endoplasmic reticulum stress and apoptosis. , 2007, Investigative Ophthalmology and Visual Science.

[9]  Andrew J. Hill,et al.  Analysis of protein-coding genetic variation in 60,706 humans , 2015, bioRxiv.

[10]  Christian Darabos,et al.  Genome-Wide Genetic Interaction Analysis of Glaucoma Using Expert Knowledge Derived from Human Phenotype Networks , 2014, Pacific Symposium on Biocomputing.

[11]  Janey L. Wiggs,et al.  Common variants near ABCA1, AFAP1 and GMDS confer risk of primary open-angle glaucoma , 2014, Nature Genetics.

[12]  T. Mackay,et al.  Overexpression of Myocilin in the Drosophila Eye Activates the Unfolded Protein Response: Implications for Glaucoma , 2009, PloS one.

[13]  Ching-Chyuan Su,et al.  Mutations in the CYP1B1 gene may contribute to juvenile-onset open-angle glaucoma , 2012, Eye.

[14]  Adriana I. Iglesias,et al.  Meta‐analysis of Genome‐Wide Association Studies Identifies Novel Loci Associated With Optic Disc Morphology , 2015, Genetic epidemiology.

[15]  Benno Schwikowski,et al.  Discovering regulatory and signalling circuits in molecular interaction networks , 2002, ISMB.

[16]  HighWire Press.  The journal of neuroscience : the official journal of the Society for Neuroscience. , 1981 .

[17]  D. Mackey,et al.  Occurrence of CYP1B1 Mutations in Juvenile Open-Angle Glaucoma With Advanced Visual Field Loss. , 2015, JAMA ophthalmology.

[18]  S. Hodgson,et al.  The human PAX6 gene is mutated in two patients with aniridia , 1992, Nature Genetics.

[19]  M. Spitz,et al.  Shifting paradigm of association studies: value of rare single-nucleotide polymorphisms. , 2008, American journal of human genetics.

[20]  J. Roth,et al.  Aggregated myocilin induces russell bodies and causes apoptosis: implications for the pathogenesis of myocilin-caused primary open-angle glaucoma. , 2007, The American journal of pathology.

[21]  Michael G. Anderson,et al.  Topical ocular sodium 4-phenylbutyrate rescues glaucoma in a myocilin mouse model of primary open-angle glaucoma. , 2012, Investigative ophthalmology & visual science.

[22]  Young H. Kwon,et al.  Copy number variations on chromosome 12 q 14 in patients with normal tension glaucoma , 2011 .

[23]  Jochen Rieck.  The pathogenesis of glaucoma in the interplay with the immune system. , 2013, Investigative ophthalmology & visual science.

[24]  V. Sheffield,et al.  Identification of a Gene That Causes Primary Open Angle Glaucoma , 1997, Science.

[25]  M. Sarfarazi,et al.  Identification of three different truncating mutations in cytochrome P4501B1 (CYP1B1) as the principal cause of primary congenital glaucoma (Buphthalmos) in families linked to the GLC3A locus on chromosome 2p21. , 1997, Human molecular genetics.

[26]  Mauricio O. Carneiro,et al.  From FastQ Data to High‐Confidence Variant Calls: The Genome Analysis Toolkit Best Practices Pipeline , 2013, Current protocols in bioinformatics.

[27]  Kang Zhang,et al.  Genes of the Unfolded Protein Response Pathway Harbor Risk Alleles for Primary Open Angle Glaucoma , 2011, PloS one.

[28]  D. Mackey,et al.  Genome-wide association studies and gene – environment interactions , 2014 .

[29]  Matthew R. Laird,et al.  Protein Protein Interaction Network Evaluation for Identifying Potential Drug Targets , 2009 .

[30]  D. Mackey,et al.  Occurrence of CYP 1 B 1 Mutations in Juvenile Open-Angle GlaucomaWith Advanced Visual Field Loss , 2015 .

[31]  David J. Calkins,et al.  Absence of Transient Receptor Potential Vanilloid-1 Accelerates Stress-Induced Axonopathy in the Optic Projection , 2014, The Journal of Neuroscience.

[32]  A. Lotery,et al.  Current concepts on primary open-angle glaucoma genetics: a contribution to disease pathophysiology and future treatment , 2012, Eye.

[33]  H. Hakonarson,et al.  ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data , 2010, Nucleic acids research.

[34]  M. A. Carbone,et al.  A molecular mechanism for glaucoma: endoplasmic reticulum stress and the unfolded protein response. , 2013, Trends in molecular medicine.

[35]  Hye-Won Shin,et al.  ATP9B, a P4-ATPase (a Putative Aminophospholipid Translocase), Localizes to the trans-Golgi Network in a CDC50 Protein-independent Manner* , 2011, The Journal of Biological Chemistry.

[36]  Catherine Cheng,et al.  Cataracts and Microphthalmia Caused by a Gja8 Mutation in Extracellular Loop 2 , 2012, PloS one.

[37]  Hiroko Ito,et al.  Activation of the TRPV1 channel attenuates N-methyl-D-aspartic acid-induced neuronal injury in the rat retina. , 2014, European journal of pharmacology.

[38]  D. Mackey,et al.  Copy number variations of TBK1 in Australian patients with primary open-angle glaucoma. , 2015, American journal of ophthalmology.

[39]  V. Bansal,et al.  Statistical analysis strategies for association studies involving rare variants , 2010, Nature Reviews Genetics.

[40]  Young H. Kwon,et al.  Copy number variations on chromosome 12q14 in patients with normal tension glaucoma. , 2011, Human molecular genetics.

[41]  P. Bork,et al.  A method and server for predicting damaging missense mutations , 2010, Nature Methods.

[42]  Kenny Q. Ye,et al.  An integrated map of genetic variation from 1,092 human genomes , 2012, Nature.

[43]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[44]  Judy H. Cho,et al.  Finding the missing heritability of complex diseases , 2009, Nature.

[45]  Hye-Won Shin,et al.  ATP9B, a P4-ATPase (a putative aminophospholipid translocase), localizes to the trans-Golgi network in a CDC50 protein-independent manner. , 2015, The Journal of Biological Chemistry.

[46]  K. Abu-Amero,et al.  CYP1B1 Mutations are a Major Contributor to Juvenile-Onset Open Angle Glaucoma in Saudi Arabia , 2015, Ophthalmic genetics.

[47]  Xianjun Zhu,et al.  Phospholipid flippase ATP8A2 is required for normal visual and auditory function and photoreceptor and spiral ganglion cell survival , 2014, Journal of Cell Science.

[48]  M. Brown,et al.  Genome-wide association study identifies susceptibility loci for open angle glaucoma at TMCO1 and CDKN2B-AS1 , 2011, Nature Genetics.

[49]  David J. Lynn,et al.  Using biological networks to integrate, visualize and analyze genomics data , 2016, Genetics Selection Evolution.

[50]  C. Pang,et al.  Correction of the disease phenotype of myocilin-causing glaucoma by a natural osmolyte. , 2009, Investigative ophthalmology & visual science.

[51]  J. Morissette,et al.  Intracellular sequestration of hetero-oligomers formed by wild-type and glaucoma-causing myocilin mutants. , 2004, Investigative ophthalmology & visual science.

[52]  M. Brown,et al.  Genome-wide association study identifies susceptibility loci for open angle glaucoma at TMCO 1 and CDKN 2 BAS 1 , 2011 .

[53]  Daniel L. Koller,et al.  Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture , 2012, Nature Genetics.

[54]  Michael G. Anderson,et al.  Reduction of ER stress via a chemical chaperone prevents disease phenotypes in a mouse model of primary open angle glaucoma. , 2011, The Journal of clinical investigation.

[55]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[56]  A. Hofman,et al.  Genetic risk of primary open-angle glaucoma. Population-based familial aggregation study. , 1998, Archives of ophthalmology.

[57]  Tin Aung,et al.  Genome-wide association analyses identify multiple loci associated with central corneal thickness and keratoconus , 2013, Nature Genetics.

[58]  James Y. Zou.  Analysis of protein-coding genetic variation in 60,706 humans , 2015, Nature.

[59]  M. Coote,et al.  Australian and New Zealand Registry of Advanced Glaucoma: methodology and recruitment , 2012, Clinical & experimental ophthalmology.

[60]  J. Blangero,et al.  The path to open-angle glaucoma gene discovery: endophenotypic status of intraocular pressure, cup-to-disc ratio, and central corneal thickness. , 2010, Investigative ophthalmology & visual science.

[61]  Brittany N. Lasseigne,et al.  Exome sequencing in amyotrophic lateral sclerosis identifies risk genes and pathways , 2015, Science.

[62]  Sten Orrenius,et al.  Calcium: Regulation of cell death: the calcium–apoptosis link , 2003, Nature Reviews Molecular Cell Biology.

[63]  Xueshan Xiao,et al.  Mutation analysis of the genes associated with anterior segment dysgenesis, microcornea and microphthalmia in 257 patients with glaucoma. , 2015, International journal of molecular medicine.

[64]  D. Pascolini,et al.  Global estimates of visual impairment: 2010 , 2011, British Journal of Ophthalmology.

[65]  T. Jakobs.  Differential gene expression in glaucoma. , 2014, Cold Spring Harbor perspectives in medicine.

[66]  Eric S. Lander,et al.  A polygenic burden of rare disruptive mutations in schizophrenia , 2014, Nature.

[67]  L. Ruggiero,et al.  Diurnal, localized exposure of phosphatidylserine by rod outer segment tips in wild-type but not Itgb5−/− or Mfge8−/− mouse retina , 2012, Proceedings of the National Academy of Sciences of the United States of America.