Search

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. Author summary Vessels at the back of the eye (the “retina”) can be imaged easily. This paper reports on the largest genetic study of retinal vessel shape and size characteristics so far undertaken, to the best of our knowledge. Our study is novel in using an automated artificial intelligence imaging approach to distinguish between arteries and veins, and in demonstrating more genetic associations with vessel characteristics than any previous study (119 genetic loci in all). We also show that the tortuosity of retinal arteries is the most strongly genetically determined vessel characteristic (replicated remarkable well in a separate second large dataset). In addition, using a particular type of genetic analysis (so called “Mendelian Randomization”) we show for the first time that the tortuosity of arteries in the retina is causally related to elevated diastolic blood pressure and not the other way around. This is important as it provides unique insights into the mechanism of elevated blood pressure and hypertension, providing pointers to novel therapeutic targets for future treatment.

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.

This study investigated temporal trends in the epidemiology of primary myopia and associations with key environmental risk factors in a UK population. Data were collected at recruitment (non-cycloplegic autorefraction, year of birth, sex, ethnicity, highest educational attainment, reason and age of first wearing glasses and history of eye disease) from 107,442 UK Biobank study participants aged 40 to 69 years, born between 1939 and 1970. Myopia was defined as mean spherical equivalent (MSE) ≤-1 dioptre (D). Temporal changes in myopia frequency by birth cohort (5-year bands using date of birth) and associations with environmental factors were analysed, distinguishing both type (childhood-onset,

PurposeGenetic variants identified through population-based genome-wide studies are generally of high frequency, exerting their action in the central part of the refractive error spectrum. However, the power to identify associations with variants of lower minor allele frequency is greatly reduced, requiring considerable sample sizes. Here we aim to assess the impact of rare variants on genetic variation of refractive errors in a very large general population cohort.MethodsGenetic association analyses of non-cyclopaedic autorefraction calculated as mean spherical equivalent (SPHE) used whole-exome sequence genotypic information from 50,893 unrelated participants in the UK Biobank of European ancestry. Gene-based analyses tested for association with SPHE using an optimised SNP-set kernel association test (SKAT-O) restricted to rare variants (minor allele frequency < 1%) within protein-coding regions of the genome. All models were adjusted for age, sex and common lead variants within the same locus reported by previous genome-wide association studies. Potentially causal markers driving association at significant loci were elucidated using sensitivity analyses by sequentially dropping the most associated variants from gene-based analyses.ResultsWe found strong statistical evidence for association of SPHE with the SIX6 (p-value = 2.15 x 10-10, or Bonferroni-Corrected p = 4.41x10-06) and the CRX gene (p-value = 6.65 x 10-08, or Bonferroni-Corrected p = 0.001). The SIX6 gene codes for a transcription factor believed to be critical to the eye, retina and optic disc development and morphology, while CRX regulates photoreceptor specification and expression of over 700 genes in the retina. These novel associations suggest an important role of genes involved in eye morphogenesis in refractive error.ConclusionThe results of our study support previous research highlighting the importance of rare variants to the genetic risk of refractive error. We explain some of the origins of the genetic signals seen in GWAS but also report for the first time a completely novel association with the CRX gene.

Keratoconus is a relatively frequent disease leading to severe visual impairment. Existing therapies are imperfect and clinical management may benefit from improved understanding of mechanisms leading to this disease. We aim to investigate the replication of 11 single nucleotide polymorphisms (SNPs) with keratoconus.SNPs from loci previously found in association with keratoconus were genotyped in 165 keratoconus cases of Caucasian Czech origin (108 males and 57 females) and 193 population and gender-matched controls. They included rs1536482 (COL5A1), rs4839200 (KCND3), rs757219 and rs214884 (IMMP2L), rs1328083 and rs1328089 (DAOA), rs2721051 (FOXO1), rs4894535 (FNDC3B), rs4954218 (MAP3K19, RAB3GAP1), rs9938149 (ZNF469) and rs1324183 (MPDZ). A case-control association analysis was assessed using Fisher's exact tests.The strongest association was found for rs1324183 (allelic test OR = 1.58; 95% CI, 1.10-2.24, p = 0.01). Statistically significant values were also obtained for rs2721051 (allelic test OR = 1.72; 95% CI, 1.07-2.77, p = 0.025) and rs4954218 (allelic test OR = 1.53; 95% CI, 1.01-2.34; p = 0.047) which showed an opposite effect direction compared to previously reported one.Independent replication of association between two SNPs and keratoconus supports the association of these loci with the risks for the disease development, while the effect of rs4954218 warrants further investigation. Understanding the role of the genetic factors involved in keratoconus etiopathogenesis may facilitate development of novel therapies and an early detection.

An increasing number of genome-wide association (GWA) studies are now using the higher resolution 1000 Genomes Project reference panel (1000G) for imputation, with the expectation that 1000G imputation will lead to the discovery of additional associated loci when compared to HapMap imputation. In order to assess the improvement of 1000G over HapMap imputation in identifying associated loci, we compared the results of GWA studies of circulating fibrinogen based on the two reference panels. Using both HapMap and 1000G imputation we performed a meta-analysis of 22 studies comprising the same 91,953 individuals. We identified six additional signals using 1000G imputation, while 29 loci were associated using both HapMap and 1000G imputation. One locus identified using HapMap imputation was not significant using 1000G imputation. The genome-wide significance threshold of 5×10-8 is based on the number of independent statistical tests using HapMap imputation, and 1000G imputation may lead to further independent tests that should be corrected for. When using a stricter Bonferroni correction for the 1000G GWA study (P-value < 2.5×10-8), the number of loci significant only using HapMap imputation increased to 4 while the number of loci significant only using 1000G decreased to 5. In conclusion, 1000G imputation enabled the identification of 20% more loci than HapMap imputation, although the advantage of 1000G imputation became less clear when a stricter Bonferroni correction was used. More generally, our results provide insights that are applicable to the implementation of other dense reference panels that are under development.

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

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHRadjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHRadjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHRadjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR

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.

PURPOSE:To report the methodology and findings of a large scale investigation of burden and distribution of refractive error, from a contemporary and ethnically diverse study of health and disease in adults, in the UK. METHODS:U K Biobank, a unique contemporary resource for the study of health and disease, recruited more than half a million people aged 40-69 years. A subsample of 107,452 subjects undertook an enhanced ophthalmic examination which provided autorefraction data (a measure of refractive error). Refractive error status was categorised using the mean spherical equivalent refraction measure. Information on socio-demographic factors (age, gender, ethnicity, educational qualifications and accommodation tenure) was reported at the time of recruitment by questionnaire and face-to-face interview. RESULTS:Fifty four percent of participants aged 40-69 years had refractive error. Specifically 27% had myopia (4% high myopia), which was more common amongst younger people, those of higher socio-economic status, higher educational attainment, or of White or Chinese ethnicity. The frequency of hypermetropia increased with age (7% at 40-44 years increasing to 46% at 65-69 years), was higher in women and its severity was associated with ethnicity (moderate or high hypermetropia at least 30% less likely in non-White ethnic groups compared to White). CONCLUSIONS:Refractive error is a significant public health issue for the UK and this study provides contemporary data on adults for planning services, health economic modelling and monitoring of secular trends. Further investigation of risk factors is necessary to inform strategies for prevention. There is scope to do this through the planned longitudinal extension of the UK Biobank study.

We describe novel CHRDL1 mutations in ten families with X-linked megalocornea (MGC1). Our mutation-positive cohort enabled us to establish ultrasonography as a reliable clinical diagnostic tool to distinguish between MGC1 and primary congenital glaucoma (PCG). Megalocornea is also a feature of Neuhäuser or megalocornea-mental retardation (MMR) syndrome, a rare condition of unknown etiology. In a male patient diagnosed with MMR, we performed targeted and whole exome sequencing (WES) and identified a novel missense mutation in CHRDL1 that accounts for his MGC1 phenotype but not his non-ocular features. This finding suggests that MMR syndrome, in some cases, may be di- or multigenic. MGC1 patients have reduced central corneal thickness (CCT); however no X-linked loci have been associated with CCT, possibly because the majority of genome-wide association studies (GWAS) overlook the X-chromosome. We therefore explored whether variants on the X-chromosome are associated with CCT. We found rs149956316, in intron 6 of CHRDL1, to be the most significantly associated single nucleotide polymorphism (SNP) (p = 6.81×10(-6)) on the X-chromosome. However, this association was not replicated in a smaller subset of whole genome sequenced samples. This study highlights the importance of including X-chromosome SNP data in GWAS to identify potential loci associated with quantitative traits or disease risk.

Refractive error (RE) is a complex, multifactorial disorder characterized by a mismatch between the optical power of the eye and its axial length that causes object images to be focused off the retina. The two major subtypes of RE are myopia (nearsightedness) and hyperopia (farsightedness), which represent opposite ends of the distribution of the quantitative measure of spherical refraction. We performed a fixed effects meta-analysis of genome-wide association results of myopia and hyperopia from 9 studies of European-derived populations: AREDS, KORA, FES, OGP-Talana, MESA, RSI, RSII, RSIII and ERF. One genome-wide significant region was observed for myopia, corresponding to a previously identified myopia locus on 8q12 (p = 1.25×10(-8)), which has been reported by Kiefer et al. as significantly associated with myopia age at onset and Verhoeven et al. as significantly associated to mean spherical-equivalent (MSE) refractive error. We observed two genome-wide significant associations with hyperopia. These regions overlapped with loci on 15q14 (minimum p value = 9.11×10(-11)) and 8q12 (minimum p value 1.82×10(-11)) previously reported for MSE and myopia age at onset. We also used an intermarker linkage- disequilibrium-based method for calculating the effective number of tests in targeted regional replication analyses. We analyzed myopia (which represents the closest phenotype in our data to the one used by Kiefer et al.) and showed replication of 10 additional loci associated with myopia previously reported by Kiefer et al. This is the first replication of these loci using myopia as the trait under analysis. "Replication-level" association was also seen between hyperopia and 12 of Kiefer et al.'s published loci. For the loci that show evidence of association to both myopia and hyperopia, the estimated effect of the risk alleles were in opposite directions for the two traits. This suggests that these loci are important contributors to variation of refractive error across the distribution.

Inter-individual variation in facial shape is one of the most noticeable phenotypes in humans, and it is clearly under genetic regulation; however, almost nothing is known about the genetic basis of normal human facial morphology. We therefore conducted a genome-wide association study for facial shape phenotypes in multiple discovery and replication cohorts, considering almost ten thousand individuals of European descent from several countries. Phenotyping of facial shape features was based on landmark data obtained from three-dimensional head magnetic resonance images (MRIs) and two-dimensional portrait images. We identified five independent genetic loci associated with different facial phenotypes, suggesting the involvement of five candidate genes--PRDM16, PAX3, TP63, C5orf50, and COL17A1--in the determination of the human face. Three of them have been implicated previously in vertebrate craniofacial development and disease, and the remaining two genes potentially represent novel players in the molecular networks governing facial development. Our finding at PAX3 influencing the position of the nasion replicates a recent GWAS of facial features. In addition to the reported GWA findings, we established links between common DNA variants previously associated with NSCL/P at 2p21, 8q24, 13q31, and 17q22 and normal facial-shape variations based on a candidate gene approach. Overall our study implies that DNA variants in genes essential for craniofacial development contribute with relatively small effect size to the spectrum of normal variation in human facial morphology. This observation has important consequences for future studies aiming to identify more genes involved in the human facial morphology, as well as for potential applications of DNA prediction of facial shape such as in future forensic applications.

Intraocular pressure (IOP) is a highly heritable risk factor for primary open-angle glaucoma and is the only target for current glaucoma therapy. The genetic factors which determine IOP are largely unknown. We performed a genome-wide association study for IOP in 11,972 participants from 4 independent population-based studies in The Netherlands. We replicated our findings in 7,482 participants from 4 additional cohorts from the UK, Australia, Canada, and the Wellcome Trust Case-Control Consortium 2/Blue Mountains Eye Study. IOP was significantly associated with rs11656696, located in GAS7 at 17p13.1 (p=1.4×10(-8)), and with rs7555523, located in TMCO1 at 1q24.1 (p=1.6×10(-8)). In a meta-analysis of 4 case-control studies (total N = 1,432 glaucoma cases), both variants also showed evidence for association with glaucoma (p=2.4×10(-2) for rs11656696 and p=9.1×10(-4) for rs7555523). GAS7 and TMCO1 are highly expressed in the ciliary body and trabecular meshwork as well as in the lamina cribrosa, optic nerve, and retina. Both genes functionally interact with known glaucoma disease genes. These data suggest that we have identified two clinically relevant genes involved in IOP regulation.

Genome-wide association studies have identified numerous genetic loci for spirometic measures of pulmonary function, forced expiratory volume in one second (FEV(1)), and its ratio to forced vital capacity (FEV(1)/FVC). Given that cigarette smoking adversely affects pulmonary function, we conducted genome-wide joint meta-analyses (JMA) of single nucleotide polymorphism (SNP) and SNP-by-smoking (ever-smoking or pack-years) associations on FEV(1) and FEV(1)/FVC across 19 studies (total N = 50,047). We identified three novel loci not previously associated with pulmonary function. SNPs in or near DNER (smallest P(JMA = )5.00×10(-11)), HLA-DQB1 and HLA-DQA2 (smallest P(JMA = )4.35×10(-9)), and KCNJ2 and SOX9 (smallest P(JMA = )1.28×10(-8)) were associated with FEV(1)/FVC or FEV(1) in meta-analysis models including SNP main effects, smoking main effects, and SNP-by-smoking (ever-smoking or pack-years) interaction. The HLA region has been widely implicated for autoimmune and lung phenotypes, unlike the other novel loci, which have not been widely implicated. We evaluated DNER, KCNJ2, and SOX9 and found them to be expressed in human lung tissue. DNER and SOX9 further showed evidence of differential expression in human airway epithelium in smokers compared to non-smokers. Our findings demonstrated that joint testing of SNP and SNP-by-environment interaction identified novel loci associated with complex traits that are missed when considering only the genetic main effects.

Elevated levels of acute-phase serum amyloid A (A-SAA) cause amyloidosis and are a risk factor for atherosclerosis and its clinical complications, type 2 diabetes, as well as various malignancies. To investigate the genetic basis of A-SAA levels, we conducted the first genome-wide association study on baseline A-SAA concentrations in three population-based studies (KORA, TwinsUK, Sorbs) and one prospective case cohort study (LURIC), including a total of 4,212 participants of European descent, and identified two novel genetic susceptibility regions at 11p15.5-p13 and 1p31. The region at 11p15.5-p13 (rs4150642; p = 3.20×10(-111)) contains serum amyloid A1 (SAA1) and the adjacent general transcription factor 2 H1 (GTF2H1), Hermansky-Pudlak Syndrome 5 (HPS5), lactate dehydrogenase A (LDHA), and lactate dehydrogenase C (LDHC). This region explains 10.84% of the total variation of A-SAA levels in our data, which makes up 18.37% of the total estimated heritability. The second region encloses the leptin receptor (LEPR) gene at 1p31 (rs12753193; p = 1.22×10(-11)) and has been found to be associated with CRP and fibrinogen in previous studies. Our findings demonstrate a key role of the 11p15.5-p13 region in the regulation of baseline A-SAA levels and provide confirmative evidence of the importance of the 1p31 region for inflammatory processes and the close interplay between A-SAA, leptin, and other acute-phase proteins.

Previous studies have successfully identified genetic variants in several genes associated with human iris (eye) color; however, they all used simplified categorical trait information. Here, we quantified continuous eye color variation into hue and saturation values using high-resolution digital full-eye photographs and conducted a genome-wide association study on 5,951 Dutch Europeans from the Rotterdam Study. Three new regions, 1q42.3, 17q25.3, and 21q22.13, were highlighted meeting the criterion for genome-wide statistically significant association. The latter two loci were replicated in 2,261 individuals from the UK and in 1,282 from Australia. The LYST gene at 1q42.3 and the DSCR9 gene at 21q22.13 serve as promising functional candidates. A model for predicting quantitative eye colors explained over 50% of trait variance in the Rotterdam Study. Over all our data exemplify that fine phenotyping is a useful strategy for finding genes involved in human complex traits.

Central corneal thickness (CCT), one of the most highly heritable human traits (h(2) typically>0.9), is important for the diagnosis of glaucoma and a potential risk factor for glaucoma susceptibility. We conducted genome-wide association studies in five cohorts from Australia and the United Kingdom (total N = 5058). Three cohorts were based on individually genotyped twin collections, with the remaining two cohorts genotyped on pooled samples from singletons with extreme trait values. The pooled sample findings were validated by individual genotyping the pooled samples together with additional samples also within extreme quantiles. We describe methods for efficient combined analysis of the results from these different study designs. We have identified and replicated quantitative trait loci on chromosomes 13 and 16 for association with CCT. The locus on chromosome 13 (nearest gene FOXO1) had an overall meta-analysis p-value for all the individually genotyped samples of 4.6x10(-10). The locus on chromosome 16 was associated with CCT with p = 8.95x10(-11). The nearest gene to the associated chromosome 16 SNPs was ZNF469, a locus recently implicated in Brittle Cornea Syndrome (BCS), a very rare disorder characterized by abnormal thin corneas. Our findings suggest that in addition to rare variants in ZNF469 underlying CCT variation in BCS patients, more common variants near this gene may contribute to CCT variation in the general population.