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.

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.

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.

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.

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 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.

Genome-wide association study identifies nine novel loci for 2D:4D finger ratio, a putative retrospective biomarker of testosterone exposure in utero.

The ratio of the length of the index finger to that of the ring finger (2D:4D) is sexually dimorphic and is commonly used as a non-invasive biomarker of prenatal androgen exposure. Most association studies of 2D:4D ratio with a diverse range of sex-specific traits have typically involved small sample sizes and have been difficult to replicate, raising questions around the utility and precise meaning of the measure. In the largest genome-wide association meta-analysis of 2D:4D ratio to date (N = 15 661, with replication N = 75 821), we identified 11 loci (9 novel) explaining 3.8% of the variance in mean 2D:4D ratio. We also found weak evidence for association (ß = 0.06; P = 0.02) between 2D:4D ratio and sensitivity to testosterone [length of the CAG microsatellite repeat in the androgen receptor (AR) gene] in females only. Furthermore, genetic variants associated with (adult) testosterone levels and/or sex hormone-binding globulin were not associated with 2D:4D ratio in our sample. Although we were unable to find strong evidence from our genetic study to support the hypothesis that 2D:4D ratio is a direct biomarker of prenatal exposure to androgens in healthy individuals, our findings do not explicitly exclude this possibility, and pathways involving testosterone may become apparent as the size of the discovery sample increases further. Our findings provide new insight into the underlying biology shaping 2D:4D variation in the general population.

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.

Comparison of Associations with Different Macular Inner Retinal Thickness Parameters in a Large Cohort: The UK Biobank.

PURPOSE: To describe and compare associations with macular retinal nerve fiber layer (mRNFL), ganglion cell complex (GCC), and ganglion cell-inner plexiform layer (GCIPL) thicknesses in a large cohort. DESIGN: Cross-sectional study. PARTICIPANTS: We included 42 044 participants in the UK Biobank. The mean age was 56 years. METHODS: Spectral-domain OCT macular images were segmented and analyzed. Corneal-compensated intraocular pressure (IOPcc) was measured with the Ocular Response Analyzer (Reichert, Corp., Buffalo, NY). Multivariable linear regression was used to examine associations with mean mRNFL, GCC, and GCIPL thicknesses. Factors examined were age, sex, ethnicity, height, body mass index (BMI), smoking status, alcohol intake, Townsend deprivation index, education level, diabetes status, spherical equivalent, and IOPcc. MAIN OUTCOME MEASURES: Thicknesses of mRNFL, GCC, and GCIPL. RESULTS: We identified several novel independent associations with thinner inner retinal thickness. Thinner inner retina was associated with alcohol intake (most significant for GCIPL: -0.46 µm for daily or almost daily intake compared with special occasion only or never [95% confidence interval (CI), 0.61-0.30]; P = 1.1×10-8), greater social deprivation (most significant for GCIPL: -0.28 µm for most deprived quartile compared with least deprived quartile [95% CI, -0.42 to -0.14]; P = 6.6×10-5), lower educational attainment (most significant for mRNFL: -0.36 µm for less than O level compared with degree level [95% CI, -0.45 to 0.26]; P = 2.3×10-14), and nonwhite ethnicity (most significant for mRNFL comparing blacks with whites: -1.65 µm [95% CI, -1.86 to -1.43]; P = 2.4×10-50). Corneal-compensated intraocular pressure was associated most significantly with GCIPL (-0.04 µm/mmHg [95% CI, -0.05 to -0.03]; P = 4.0×10-10) and was not associated significantly with mRNFL (0.00 µm/mmHg [95% CI, -0.01 to 0.01]; P = 0.77). The variables examined explained a greater proportion of the variance of GCIPL (11%) than GCC (6%) or mRNFL (7%). CONCLUSIONS: The novel associations we identified may be important to consider when using inner retinal parameters as a diagnostic tool. Associations generally were strongest with GCIPL, particularly for IOP. This suggests that GCIPL may be the superior inner retinal biomarker for macular pathophysiologic processes and especially for glaucoma.

A commonly occurring genetic variant within the NPLOC4-TSPAN10-PDE6G gene cluster is associated with the risk of strabismus.

Strabismus refers to an abnormal alignment of the eyes leading to the loss of central binocular vision. Concomitant strabismus occurs when the angle of deviation is constant in all positions of gaze and often manifests in early childhood when it is considered to be a neurodevelopmental disorder of the visual system. As such, it is inherited as a complex genetic trait, affecting 2-4% of the population. A genome-wide association study (GWAS) for self-reported strabismus (1345 cases and 65,349 controls from UK Biobank) revealed a single genome-wide significant locus on chromosome 17q25. Approximately 20 variants across the NPLOC4-TSPAN10-PDE6G gene cluster and in almost perfect linkage disequilibrium (LD) were most strongly associated (lead variant: rs75078292, OR = 1.26, p = 2.24E-08). A recessive model provided a better fit to the data than an additive model. Association with strabismus was independent of refractive error, and the degree of association with strabismus was minimally attenuated after adjustment for amblyopia. The association with strabismus was replicated in an independent cohort of clinician-diagnosed children aged 7 years old (116 cases and 5084 controls; OR = 1.85, p = 0.009). The associated variants included 2 strong candidate causal variants predicted to have functional effects: rs6420484, which substitutes tyrosine for a conserved cysteine (C177Y) in the TSPAN10 gene, and a 4-bp deletion variant, rs397693108, predicted to cause a frameshift in TSPAN10. The population-attributable risk for the locus was approximately 8.4%, indicating an important role in conferring susceptibility to strabismus.

Genome-wide association study of primary open-angle glaucoma in continental and admixed African populations.

Primary open angle glaucoma (POAG) is a complex disease with a major genetic contribution. Its prevalence varies greatly among ethnic groups, and is up to five times more frequent in black African populations compared to Europeans. So far, worldwide efforts to elucidate the genetic complexity of POAG in African populations has been limited. We conducted a genome-wide association study in 1113 POAG cases and 1826 controls from Tanzanian, South African and African American study samples. Apart from confirming evidence of association at TXNRD2 (rs16984299; OR[T] 1.20; P = 0.003), we found that a genetic risk score combining the effects of the 15 previously reported POAG loci was significantly associated with POAG in our samples (OR 1.56; 95% CI 1.26-1.93; P = 4.79 × 10-5). By genome-wide association testing we identified a novel candidate locus, rs141186647, harboring EXOC4 (OR[A] 0.48; P = 3.75 × 10-8), a gene transcribing a component of the exocyst complex involved in vesicle transport. The low frequency and high degree of genetic heterogeneity at this region hampered validation of this finding in predominantly West-African replication sets. Our results suggest that established genetic risk factors play a role in African POAG, however, they do not explain the higher disease load. The high heterogeneity within Africans remains a challenge to identify the genetic commonalities for POAG in this ethnicity, and demands studies of extremely large size.

A genome-wide association study of corneal astigmatism: The CREAM Consortium.

Purpose: To identify genes and genetic markers associated with corneal astigmatism. Methods: A meta-analysis of genome-wide association studies (GWASs) of corneal astigmatism undertaken for 14 European ancestry (n=22,250) and 8 Asian ancestry (n=9,120) cohorts was performed by the Consortium for Refractive Error and Myopia. Cases were defined as having >0.75 diopters of corneal astigmatism. Subsequent gene-based and gene-set analyses of the meta-analyzed results of European ancestry cohorts were performed using VEGAS2 and MAGMA software. Additionally, estimates of single nucleotide polymorphism (SNP)-based heritability for corneal and refractive astigmatism and the spherical equivalent were calculated for Europeans using LD score regression. Results: The meta-analysis of all cohorts identified a genome-wide significant locus near the platelet-derived growth factor receptor alpha (PDGFRA) gene: top SNP: rs7673984, odds ratio=1.12 (95% CI:1.08-1.16), p=5.55×10-9. No other genome-wide significant loci were identified in the combined analysis or European/Asian ancestry-specific analyses. Gene-based analysis identified three novel candidate genes for corneal astigmatism in Europeans-claudin-7 (CLDN7), acid phosphatase 2, lysosomal (ACP2), and TNF alpha-induced protein 8 like 3 (TNFAIP8L3). Conclusions: In addition to replicating a previously identified genome-wide significant locus for corneal astigmatism near the PDGFRA gene, gene-based analysis identified three novel candidate genes, CLDN7, ACP2, and TNFAIP8L3, that warrant further investigation to understand their role in the pathogenesis of corneal astigmatism. The much lower number of genetic variants and genes demonstrating an association with corneal astigmatism compared to published spherical equivalent GWAS analyses suggest a greater influence of rare genetic variants, non-additive genetic effects, or environmental factors in the development of astigmatism.

APLP2 Regulates Refractive Error and Myopia Development in Mice and Humans.

Myopia is the most common vision disorder and the leading cause of visual impairment worldwide. However, gene variants identified to date explain less than 10% of the variance in refractive error, leaving the majority of heritability unexplained ("missing heritability"). Previously, we reported that expression of APLP2 was strongly associated with myopia in a primate model. Here, we found that low-frequency variants near the 5'-end of APLP2 were associated with refractive error in a prospective UK birth cohort (n = 3,819 children; top SNP rs188663068, p = 5.0 × 10-4) and a CREAM consortium panel (n = 45,756 adults; top SNP rs7127037, p = 6.6 × 10-3). These variants showed evidence of differential effect on childhood longitudinal refractive error trajectories depending on time spent reading (gene x time spent reading x age interaction, p = 4.0 × 10-3). Furthermore, Aplp2 knockout mice developed high degrees of hyperopia (+11.5 ± 2.2 D, p < 1.0 × 10-4) compared to both heterozygous (-0.8 ± 2.0 D, p < 1.0 × 10-4) and wild-type (+0.3 ± 2.2 D, p < 1.0 × 10-4) littermates and exhibited a dose-dependent reduction in susceptibility to environmentally induced myopia (F(2, 33) = 191.0, p < 1.0 × 10-4). This phenotype was associated with reduced contrast sensitivity (F(12, 120) = 3.6, p = 1.5 × 10-4) and changes in the electrophysiological properties of retinal amacrine cells, which expressed Aplp2. This work identifies APLP2 as one of the "missing" myopia genes, demonstrating the importance of a low-frequency gene variant in the development of human myopia. It also demonstrates an important role for APLP2 in refractive development in mice and humans, suggesting a high level of evolutionary conservation of the signaling pathways underlying refractive eye development.

Haplotype reference consortium panel: Practical implications of imputations with large reference panels.

Recently, the Haplotype Reference Consortium (HRC) released a large imputation panel that allows more accurate imputation of genetic variants. In this study, we compared a set of directly assayed common and rare variants from an exome array to imputed genotypes, that is, 1000 genomes project (1000GP) and HRC. We showed that imputation using the HRC panel improved the concordance between assayed and imputed genotypes at common, and especially, low-frequency variants. Furthermore, we performed a genome-wide association meta-analysis of vertical cup-disc ratio, a highly heritable endophenotype of glaucoma, in four cohorts using 1000GP and HRC imputations. We compared the results of the meta-analysis using 1000GP to the meta-analysis results using HRC. Overall, we found that using HRC imputation significantly improved P values (P = 3.07 × 10-61 ), particularly for suggestive variants. Both meta-analyses were performed in the same sample size, yet we found eight genome-wide significant loci in the HRC-based meta-analysis versus seven genome-wide significant loci in the 1000GP-based meta-analysis. This study provides supporting evidence of the new avenues for gene discovery and fine mapping that the HRC imputation panel offers.