GWAS on retinal vasculometry phenotypes.

The eye is the window through which light is transmitted and visual sensory signalling originates. It is also a window through which elements of the cardiovascular and nervous systems can be directly inspected, using ophthalmoscopy or retinal imaging. Measurements of ocular parameters may therefore offer important information on the physiology and homeostasis of these two important systems. Here we report the results of a genetic characterisation of retinal vasculature. Four genome-wide association studies performed on different aspects of retinal vasculometry phenotypes, such as arteriolar and venular tortuosity and width, found significant similarities between retinal vascular characteristics and cardiometabolic health. Our analyses identified 119 different regions of association with traits of retinal vasculature, including 89 loci associated arteriolar tortuosity, the strongest of which was rs35131825 (p = 2.00×10-108), 2 loci with arteriolar width (rs12969347, p = 3.30×10-09 and rs5442, p = 1.9E-15), 17 other loci associated with venular tortuosity and 11 novel associations with venular width. Our causal inference analyses also found that factors linked to arteriolar tortuosity cause elevated diastolic blood pressure and not vice versa.

Combined genome-wide association study of 136 quantitative ear morphology traits in multiple populations reveal 8 novel loci.

Human ear morphology, a complex anatomical structure represented by a multidimensional set of correlated and heritable phenotypes, has a poorly understood genetic architecture. In this study, we quantitatively assessed 136 ear morphology traits using deep learning analysis of digital face images in 14,921 individuals from five different cohorts in Europe, Asia, and Latin America. Through GWAS meta-analysis and C-GWASs, a recently introduced method to effectively combine GWASs of many traits, we identified 16 genetic loci involved in various ear phenotypes, eight of which have not been previously associated with human ear features. Our findings suggest that ear morphology shares genetic determinants with other surface ectoderm-derived traits such as facial variation, mono eyebrow, and male pattern baldness. Our results enhance the genetic understanding of human ear morphology and shed light on the shared genetic contributors of different surface ectoderm-derived phenotypes. Additionally, gene editing experiments in mice have demonstrated that knocking out the newly ear-associated gene (Intu) and a previously ear-associated gene (Tbx15) causes deviating mouse ear morphology.

Education interacts with genetic variants near GJD2, RBFOX1, LAMA2, KCNQ5 and LRRC4C to confer susceptibility to myopia.

Myopia most often develops during school age, with the highest incidence in countries with intensive education systems. Interactions between genetic variants and educational exposure are hypothesized to confer susceptibility to myopia, but few such interactions have been identified. Here, we aimed to identify genetic variants that interact with education level to confer susceptibility to myopia. Two groups of unrelated participants of European ancestry from UK Biobank were studied. A 'Stage-I' sample of 88,334 participants whose refractive error (avMSE) was measured by autorefraction and a 'Stage-II' sample of 252,838 participants who self-reported their age-of-onset of spectacle wear (AOSW) but who did not undergo autorefraction. Genetic variants were prioritized via a 2-step screening process in the Stage-I sample: Step 1 was a genome-wide association study for avMSE; Step 2 was a variance heterogeneity analysis for avMSE. Genotype-by-education interaction tests were performed in the Stage-II sample, with University education coded as a binary exposure. On average, participants were 58 years-old and left full-time education when they were 18 years-old; 35% reported University level education. The 2-step screening strategy in the Stage-I sample prioritized 25 genetic variants (GWAS P < 1e-04; variance heterogeneity P < 5e-05). In the Stage-II sample, 19 of the 25 (76%) genetic variants demonstrated evidence of variance heterogeneity, suggesting the majority were true positives. Five genetic variants located near GJD2, RBFOX1, LAMA2, KCNQ5 and LRRC4C had evidence of a genotype-by-education interaction in the Stage-II sample (P < 0.002) and consistent evidence of a genotype-by-education interaction in the Stage-I sample. For all 5 variants, University-level education was associated with an increased effect of the risk allele. In this cohort, additional years of education were associated with an enhanced effect of genetic variants that have roles including axon guidance and the development of neuronal synapses and neural circuits.

Electrical responses from human retinal cone pathways associate with a common genetic polymorphism implicated in myopia.

Myopia is the commonest visual impairment. Several genetic loci confer risk, but mechanisms by which they do this are unknown. Retinal signals drive eye growth, and myopia usually results from an excessively long eye. The common variant most strongly associated with myopia is near the GJD2 gene, encoding connexin-36, which forms retinal gap junctions. Light-evoked responses of retinal neurons can be recorded noninvasively as the electroretinogram (ERG). We analyzed these responses from 186 adult twin volunteers who had been genotyped at this locus. Participants underwent detailed ERG recordings incorporating international standard stimuli as well as experimental protocols aiming to separate dark-adapted rod- and cone-driven responses. A mixed linear model was used to explore association between allelic dosage at the locus and international standard ERG parameters after adjustment for age, sex, and family structure. Significant associations were found for parameters of light-adapted, but not dark-adapted, responses. Further investigation of isolated rod- and cone-driven ERGs confirmed associations with cone-driven, but not rod-driven, a-wave amplitudes. Comparison with responses to similar experimental stimuli from a patient with a prior central retinal artery occlusion, and from two patients with selective loss of ON-bipolar cell signals, was consistent with the associated parameters being derived from signals from cone-driven OFF-bipolar cells. Analysis of single-cell transcriptome data revealed strongest GJD2 expression in cone photoreceptors; bipolar cell expression appeared strongest in OFF-bipolar cells and weakest in rod-driven ON-bipolar cells. Our findings support a potential role for altered signaling in cone-driven OFF pathways in myopia development.

A genome-wide analysis of 340 318 participants identifies four novel loci associated with the age of first spectacle wear.

Refractive errors, particularly myopia, are the most common eye conditions, often leading to serious visual impairment. The age of onset is correlated with the severity of refractive error in adulthood observed in epidemiological and genetic studies and can be used as a proxy in refractive error genetic studies. To further elucidate genetic factors that influence refractive error, we analysed self-reported age of refractive error correction data from the UK Biobank European and perform genome-wide time-to-event analyses on the age of first spectacle wear (AFSW). Genome-wide proportional hazards ratio analyses were conducted in 340 318 European subjects. We subsequently assessed the similarities and differences in the genetic architectures of refractive error correction from different causes. All-cause AFSW was genetically strongly correlated (rg = -0.68) with spherical equivalent (the measured strength of spectacle lens required to correct the refractive error) and was used as a proxy for refractive error. Time-to-event analyses found genome-wide significant associations at 44 independent genomic loci, many of which (GJD2, LAMA2, etc.) were previously associated with refractive error. We also identified six novel regions associated with AFSW, the most significant of which was on chromosome 17q (P = 3.06 × 10-09 for rs55882072), replicating in an independent dataset. We found that genes associated with AFSW were significantly enriched for expression in central nervous system tissues and were involved in neurogenesis. This work demonstrates the merits of time-to-event study design in the genetic investigation of refractive error and contributes additional knowledge on its genetic risk factors in the general population.

Whole exome sequence analysis in 51 624 participants identifies novel genes and variants associated with refractive error and myopia.

Refractive errors are associated with a range of pathological conditions, such as myopic maculopathy and glaucoma, and are highly heritable. Studies of missense and putative loss of function (pLOF) variants identified via whole exome sequencing (WES) offer the prospect of directly implicating potentially causative disease genes. We performed a genome-wide association study for refractive error in 51 624 unrelated adults, of European ancestry, aged 40-69 years from the UK and genotyped using WES. After testing 29 179 pLOF and 495 263 missense variants, 1 pLOF and 18 missense variants in 14 distinct genomic regions were taken forward for fine-mapping analysis. This yielded 19 putative causal variants of which 18 had a posterior inclusion probability >0.5. Of the 19 putative causal variants, 12 were novel discoveries. Specific variants were associated with a more myopic refractive error, while others were associated with a more hyperopic refractive error. Association with age of onset of spectacle wear (AOSW) was examined in an independent validation sample (38 100 early AOSW cases and 74 243 controls). Of 11 novel variants that could be tested, 8 (73%) showed evidence of association with AOSW status. This work identified COL4A4 and ATM as novel candidate genes associated with refractive error. In addition, novel putative causal variants were identified in the genes RASGEF1, ARMS2, BMP4, SIX6, GSDMA, GNGT2, ZNF652 and CRX. Despite these successes, the study also highlighted the limitations of community-based WES studies compared with high myopia case-control WES studies.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

Circulating metabolites modulated by diet are associated with depression.

Metabolome reflects the interplay of genome and exposome at molecular level and thus can provide deep insights into the pathogenesis of a complex disease like major depression. To identify metabolites associated with depression we performed a metabolome-wide association analysis in 13,596 participants from five European population-based cohorts characterized for depression, and circulating metabolites using ultra high-performance liquid chromatography/tandem accurate mass spectrometry (UHPLC/MS/MS) based Metabolon platform. We tested 806 metabolites covering a wide range of biochemical processes including those involved in lipid, amino-acid, energy, carbohydrate, xenobiotic and vitamin metabolism for their association with depression. In a conservative model adjusting for life style factors and cardiovascular and antidepressant medication use we identified 8 metabolites, including 6 novel, significantly associated with depression. In individuals with depression, increased levels of retinol (vitamin A), 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1) (lecithin) and mannitol/sorbitol and lower levels of hippurate, 4-hydroxycoumarin, 2-aminooctanoate (alpha-aminocaprylic acid), 10-undecenoate (11:1n1) (undecylenic acid), 1-linoleoyl-GPA (18:2) (lysophosphatidic acid; LPA 18:2) are observed. These metabolites are either directly food derived or are products of host and gut microbial metabolism of food-derived products. Our Mendelian randomization analysis suggests that low hippurate levels may be in the causal pathway leading towards depression. Our findings highlight putative actionable targets for depression prevention that are easily modifiable through diet interventions.

Cone-driven strong flash electroretinograms in healthy adults: Prevalence of negative waveforms.

PURPOSE: Both rod and cone-driven signals contribute to the electroretinogram (ERG) elicited by a standard strong flash in the dark. Negative ERGs usually reflect inner retinal dysfunction. However, in diseases where rod photoreceptor function is selectively lost, a negative waveform might represent the response of the dark-adapted cone system. To investigate the dark-adapted cone-driven waveform in healthy individuals, we delivered flashes on a dim blue background, designed to saturate the rods, but minimally adapt the cones. METHODS: ERGs were recorded, using conductive fibre electrodes, in adults from the TwinsUK cohort. Responses to 13 cd m-2 s white xenon flashes (similar to the standard DA 10 flash), delivered on a blue background, were analysed. Photopic and scotopic strengths of the background were 1.3 and 30 cd m-2, respectively; through a dilated pupil, this is expected to largely saturate the rods, but adapt the cones much less than the standard ISCEV background. RESULTS: Mean (SD) participant age was 62.5 (11.3) years (93% female). ERGs from 203 right and 204 left eyes were included, with mean (SD) b/a ratios of 1.22 (0.28) and 1.18 (0.28), respectively (medians, 1.19 and 1.17). Proportions with negative waveforms were 23 and 26%, respectively. Right and left eye b/a ratios were strongly correlated (correlation coefficient 0.74, p < 0.0001). We found no significant correlation of b/a ratio with age. CONCLUSIONS: Over 20% of eyes showed b/a ratios less than 1, consistent with the notion that dark-adapted cone-driven responses to standard bright flashes can have negative waveforms. The majority had ratios greater than 1. Thus, whilst selective loss of rod function can yield a negative waveform (with reduced a-wave) in some, our findings also suggest that loss of rod function can occur without necessarily yielding a negative ERG. One potential limitation is possible mild cone system adaptation by the background.

Correction: Genetic variation affects morphological retinal phenotypes extracted from UK Biobank optical coherence tomography images.

[This corrects the article DOI: 10.1371/journal.pgen.1009497.].

Genetic variation affects morphological retinal phenotypes extracted from UK Biobank optical coherence tomography images.

Optical Coherence Tomography (OCT) enables non-invasive imaging of the retina and is used to diagnose and manage ophthalmic diseases including glaucoma. We present the first large-scale genome-wide association study of inner retinal morphology using phenotypes derived from OCT images of 31,434 UK Biobank participants. We identify 46 loci associated with thickness of the retinal nerve fibre layer or ganglion cell inner plexiform layer. Only one of these loci has been associated with glaucoma, and despite its clear role as a biomarker for the disease, Mendelian randomisation does not support inner retinal thickness being on the same genetic causal pathway as glaucoma. We extracted overall retinal thickness at the fovea, representative of foveal hypoplasia, with which three of the 46 SNPs were associated. We additionally associate these three loci with visual acuity. In contrast to the Mendelian causes of severe foveal hypoplasia, our results suggest a spectrum of foveal hypoplasia, in part genetically determined, with consequences on visual function.

Adrenaline blocks key cell cycle genes and exhibits antifibrotic and vasoconstrictor effects in glaucoma surgery.

Adrenaline is a sympathomimetic drug used to maintain pupil dilation and to decrease the risk of bleeding. The aim of this study was to demonstrate if adrenaline could exert antifibrotic effects in glaucoma surgery. Adrenaline was tested in fibroblast-populated collagen contraction assays and there was a dose-response decrease in fibroblast contractility: matrices decreased to 47.4% (P = 0.0002) and 86.6% (P = 0.0036) with adrenaline 0.0005% and 0.01%, respectively. There was no significant decrease in cell viability even at high concentrations. Human Tenon's fibroblasts were also treated with adrenaline (0%, 0.0005%, 0.01%) for 24 h and RNA-Sequencing was performed on the Illumina NextSeq 2000. We carried out detailed gene ontology, pathway, disease and drug enrichment analyses. Adrenaline 0.01% upregulated 26 G1/S and 11 S-phase genes, and downregulated 23 G2 and 17 M-phase genes (P < 0.05). Adrenaline demonstrated similar pathway enrichment to mitosis and spindle checkpoint regulation. Adrenaline 0.05% was also injected subconjunctivally during trabeculectomy, PreserFlo Microshunt and Baerveldt 350 tube surgeries, and patients did not experience any adverse effects. Adrenaline is a safe and cheap antifibrotic drug that significantly blocks key cell cycle genes when used at high concentrations. Unless contraindicated, we recommend subconjunctival injections of adrenaline (0.05%) in all glaucoma bleb-forming surgeries.

Genome-wide association study of corneal biomechanical properties identifies over 200 loci providing insight into the genetic etiology of ocular diseases.

Corneal hysteresis and corneal resistance factor are parameters that reflect the dynamic biomechanical properties of the cornea and have been shown to be biomarkers of corneal disease. In this genome-wide association study of over 100 000 participants, we identified over 200 genetic loci, all but eight novel, significantly associated with either one or both of these traits. In addition to providing key insights into the genetic architecture underlying normal corneal function, these results identify many candidate loci in the study of corneal diseases that lead to severe visual impairment. Additionally, using Mendelian randomization, we were able to identify causal relationships between corneal biomechanics and intraocular pressure measurements, which help elucidate the relationship between corneal properties and glaucoma.

The Potential of Current Polygenic Risk Scores to Predict High Myopia and Myopic Macular Degeneration in Multiethnic Singapore Adults.

PURPOSE: To evaluate the transancestry portability of current myopia polygenic risk scores (PRSs) to predict high myopia (HM) and myopic macular degeneration (MMD) in an Asian population. DESIGN: Population-based study. PARTICIPANTS: A total of 5894 adults (2141 Chinese, 1913 Indian, and 1840 Malay) from the Singapore Epidemiology of Eye Diseases study were included in the analysis. The mean ± standard deviation age was 57.05 ± 9.31 years. A total of 361 adults had a diagnosis of HM (spherical equivalent [SE] < -5.00 diopters [D]) from refraction measurements, 240 individuals had a diagnosis of MMD graded by the International Photographic Classification and Grading System for Myopic Maculopathy criteria from fundus photographs, and 3774 individuals were control participants without myopia (SE > -0.5 D). METHODS: The PRS, derived from 687 289 HapMap3 single nucleotide polymorphisms (SNPs) from the largest genome-wide association study of myopia in Europeans to date (n = 260 974), was assessed on its ability to predict patients with HM and MMD versus control participants. MAIN OUTCOME MEASURES: The primary outcomes were the area under the receiver operating characteristic curve (AUC) to predict HM and MMD. RESULTS: The PRS had an AUC of 0.73 (95% confidence interval [CI], 0.70-0.75) for HM and 0.66 (95% CI, 0.63-0.70) for MMD versus no myopia. The inclusion of the PRS with other predictors (age, sex, educational attainment [EA], and ancestry; age-by-ancestry, sex-by-ancestry, and EA-by-ancestry interactions; and 20 genotypic principal components) increased the AUC to 0.84 (95% CI, 0.82-0.86) for HM and 0.79 (95% CI, 0.76-0.82) for MMD. Individuals with a PRS in the top 5% showed up to a 4.66 (95% CI, 3.34-6.42) times higher risk of HM developing and up to a 3.43 (95% CI, 2.27-5.05) times higher risk of MMD developing compared with the remaining 95% of individuals. CONCLUSIONS: The PRS is a good predictor for HM and facilitates the identification of high-risk children to prevent myopia progression to HM. In addition, the PRS also predicts MMD and helps to identify high-risk adults with myopia who require closer monitoring for myopia-related complications.

Alcohol, Intraocular Pressure, and Open-Angle Glaucoma: A Systematic Review and Meta-analysis.

TOPIC: This systematic review and meta-analysis summarizes the existing evidence for the association of alcohol use with intraocular pressure (IOP) and open-angle glaucoma (OAG). CLINICAL RELEVANCE: Understanding and quantifying these associations may aid clinical guidelines or treatment strategies and shed light on disease pathogenesis. The role of alcohol, a modifiable factor, in determining IOP and OAG risk also may be of interest from an individual or public health perspective. METHODS: The study protocol was preregistered in the Open Science Framework Registries (https://osf.io/z7yeg). Eligible articles (as of May 14, 2021) from 3 databases (PubMed, Embase, Scopus) were independently screened and quality assessed by 2 reviewers. All case-control, cross-sectional, and cohort studies reporting a quantitative effect estimate and 95% confidence interval (CI) for the association between alcohol use and either IOP or OAG were included. The evidence for the associations with both IOP and OAG was qualitatively summarized. Effect estimates for the association with OAG were pooled using random effects meta-analysis. Studies not meeting formal inclusion criteria for systematic review, but with pertinent results, were also appraised and discussed. Certainty of evidence was assessed using the Grading of Recommendations Assessment, Development and Evaluation (GRADE) framework. RESULTS: Thirty-four studies were included in the systematic review. Evidence from 10 studies reporting an association with IOP suggests that habitual alcohol use is associated with higher IOP and prevalence of ocular hypertension (IOP > 21 mmHg), although absolute effect sizes were small. Eleven of 26 studies, comprising 173 058 participants, that tested for an association with OAG met inclusion criteria for meta-analysis. Pooled effect estimates indicated a positive association between any use of alcohol and OAG (1.18; 95% confidence interval [CI], 1.02-1.36; P = 0.03; I2 = 40.5%), with similar estimates for both prevalent and incident OAG. The overall GRADE certainty of evidence was very low. CONCLUSIONS: Although this meta-analysis suggests a harmful association between alcohol use and OAG, our results should be interpreted cautiously given the weakness and heterogeneity of the underlying evidence base, the small absolute effect size, and the borderline statistical significance. Nonetheless, these findings may be clinically relevant, and future research should focus on improving the quality of evidence.

Genome-Wide Association Study Identifies Two Common Loci Associated with Pigment Dispersion Syndrome/Pigmentary Glaucoma and Implicates Myopia in its Development.

PURPOSE: To identify genetic variants associated with pigment dispersion syndrome (PDS) and pigmentary glaucoma (PG) in unrelated patients and to further understand the genetic and potentially causal relationships between PDS and associated risk factors. DESIGN: A 2-stage genome-wide association meta-analysis with replication and subsequent in silico analyses including Mendelian randomization. PARTICIPANTS: A total of 574 cases with PG or PDS and 52 627 controls of European descent. METHODS: Genome-wide association analyses were performed in 4 cohorts and meta-analyzed in 3 stages: (1) a discovery meta-analysis was performed in 3 cohorts, (2) replication was performed in the fourth cohort, and (3) all 4 cohorts were meta-analyzed to increase statistical power. Two-sample Mendelian randomization was used to determine whether refractive error and intraocular pressure exert causal effects over PDS. MAIN OUTCOME MEASURES: The association of genetic variants with PDS and whether myopia exerts causal effects over PDS. RESULTS: Significant association was present at 2 novel loci for PDS/PG. These loci and follow-up analyses implicate the genes gamma secretase activator protein (GSAP) (lead single nucleotide polymorphism [SNP]: rs9641220, P = 6.0×10-10) and glutamate metabotropic receptor 5 (GRM5)/TYR (lead SNP: rs661177, P = 3.9×10-9) as important factors in disease risk. Mendelian randomization showed significant evidence that negative refractive error (myopia) exerts a direct causal effect over PDS (P = 8.86×10-7). CONCLUSIONS: Common SNPs relating to the GSAP and GRM5/TYR genes are associated risk factors for the development of PDS and PG. Although myopia is a known risk factor, this study uses genetic data to demonstrate that myopia is, in part, a cause of PDS and PG.

Genome-wide association analysis of 95 549 individuals identifies novel loci and genes influencing optic disc morphology.

Optic nerve head morphology is affected by several retinal diseases. We measured the vertical optic disc diameter (DD) of the UK Biobank (UKBB) cohort (N = 67 040) and performed the largest genome-wide association study (GWAS) of DD to date. We identified 81 loci (66 novel) for vertical DD. We then replicated the novel loci in International Glaucoma Genetic Consortium (IGGC, N = 22 504) and European Prospective Investigation into Cancer-Norfolk (N = 6005); in general the concordance in effect sizes was very high (correlation in effect size estimates 0.90): 44 of the 66 novel loci were significant at P < 0.05, with 19 remaining significant after Bonferroni correction. We identified another 26 novel loci in the meta-analysis of UKBB and IGGC data. Gene-based analyses identified an additional 57 genes. Human ocular tissue gene expression analysis showed that most of the identified genes are enriched in optic nerve head tissue. Some of the identified loci exhibited pleiotropic effects with vertical cup-to-disc ratio, intraocular pressure, glaucoma and myopia. These results can enhance our understanding of the genetics of optic disc morphology and shed light on the genetic findings for other ophthalmic disorders such as glaucoma and other optic nerve diseases.

Antisense Therapy for a Common Corneal Dystrophy Ameliorates TCF4 Repeat Expansion-Mediated Toxicity.

Fuchs endothelial corneal dystrophy (FECD) is a common disease for which corneal transplantation is the only treatment option in advanced stages, and alternative treatment strategies are urgently required. Expansion (≥50 copies) of a non-coding trinucleotide repeat in TCF4 confers >76-fold risk for FECD in our large cohort of affected individuals. An FECD subject-derived corneal endothelial cell (CEC) model was developed to probe disease mechanism and investigate therapeutic approaches. The CEC model demonstrated that the repeat expansion leads to nuclear RNA foci, with the sequestration of splicing factor proteins (MBNL1 and MBNL2) to the foci and altered mRNA processing. Antisense oligonucleotide (ASO) treatment led to a significant reduction in the incidence of nuclear foci, MBNL1 recruitment to the foci, and downstream aberrant splicing events, suggesting functional rescue. This proof-of-concept study highlights the potential of a targeted ASO therapy to treat the accessible and tractable corneal tissue affected by this repeat expansion-mediated disease.

Genetics of the human face: Identification of large-effect single gene variants.

To discover specific variants with relatively large effects on the human face, we have devised an approach to identifying facial features with high heritability. This is based on using twin data to estimate the additive genetic value of each point on a face, as provided by a 3D camera system. In addition, we have used the ethnic difference between East Asian and European faces as a further source of face genetic variation. We use principal components (PCs) analysis to provide a fine definition of the surface features of human faces around the eyes and of the profile, and chose upper and lower 10% extremes of the most heritable PCs for looking for genetic associations. Using this strategy for the analysis of 3D images of 1,832 unique volunteers from the well-characterized People of the British Isles study and 1,567 unique twin images from the TwinsUK cohort, together with genetic data for 500,000 SNPs, we have identified three specific genetic variants with notable effects on facial profiles and eyes.

Meta-analysis of genome-wide association studies identifies 8 novel loci involved in shape variation of human head hair.

Shape variation of human head hair shows striking variation within and between human populations, while its genetic basis is far from being understood. We performed a series of genome-wide association studies (GWASs) and replication studies in a total of 28 964 subjects from 9 cohorts from multiple geographic origins. A meta-analysis of three European GWASs identified 8 novel loci (1p36.23 ERRFI1/SLC45A1, 1p36.22 PEX14, 1p36.13 PADI3, 2p13.3 TGFA, 11p14.1 LGR4, 12q13.13 HOXC13, 17q21.2 KRTAP, and 20q13.33 PTK6), and confirmed 4 previously known ones (1q21.3 TCHH/TCHHL1/LCE3E, 2q35 WNT10A, 4q21.21 FRAS1, and 10p14 LINC00708/GATA3), all showing genome-wide significant association with hair shape (P < 5e-8). All except one (1p36.22 PEX14) were replicated with nominal significance in at least one of the 6 additional cohorts of European, Native American and East Asian origins. Three additional previously known genes (EDAR, OFCC1, and PRSS53) were confirmed at the nominal significance level. A multivariable regression model revealed that 14 SNPs from different genes significantly and independently contribute to hair shape variation, reaching a cross-validated AUC value of 0.66 (95% CI: 0.62-0.70) and an AUC value of 0.64 in an independent validation cohort, providing an improved accuracy compared with a previous model. Prediction outcomes of 2504 individuals from a multiethnic sample were largely consistent with general knowledge on the global distribution of hair shape variation. Our study thus delivers target genes and DNA variants for future functional studies to further evaluate the molecular basis of hair shape in humans.