Smoking, Corneal Biomechanics, and Glaucoma: Results From Two Large Population-Based Cohorts.

Purpose: Smoking may influence measured IOP through an effect on corneal biomechanics, but it is unclear whether this factor translates into an increased risk for glaucoma. This study aimed to examine the association of cigarette smoking with corneal biomechanical properties and glaucoma-related traits, and to probe potential causal effects using Mendelian randomization (MR). Methods: Cross-sectional analyses within the UK Biobank (UKB) and Canadian Longitudinal Study on Aging (CLSA) cohorts. Multivariable linear and logistic regression models were used to assess associations of smoking (status, intensity, and duration) with corneal hysteresis (CH), corneal resistance factor, IOP, inner retinal thicknesses, and glaucoma. Two-sample MR analyses were performed. Results: Overall, 68,738 UKB (mean age, 56.7 years; 54.7% women) and 22 845 CLSA (mean age, 62.7 years; 49.1% women) participants were included. Compared with nonsmokers, smokers had a higher CH (UKB, +0.48 mm Hg; CLSA, +0.57 mm Hg; P < 0.001) and corneal resistance factor (UKB, +0.47 mm Hg; CLSA, +0.60 mm Hg; P < 0.001) with evidence of a dose-response effect in both studies. Differential associations with Goldmann-correlated IOP (UKB, +0.25 mm Hg; CLSA, +0.36 mm Hg; P < 0.001) and corneal-compensated IOP (UKB, -0.28 mm Hg; CLSA, -0.32 mm Hg; P ≤ 0.001) were observed. Smoking was not associated with inner retinal thicknesses or glaucoma status in either study. MR provided evidence for a causal effect of smoking on corneal biomechanics, especially higher CH. Conclusions: Cigarette smoking seems to increase corneal biomechanical resistance to deformation, but there was little evidence to support a relationship with glaucoma. This outcome may result in an artefactual association with measured IOP and could account for discordant results with glaucoma in previous epidemiological 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.

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.

Association of Behavioral and Clinical Risk Factors With Cataract: A Two-Sample Mendelian Randomization Study.

Purpose: To investigate the association of genetically determined primary open-angle glaucoma (POAG), myopic refractive error (RE), type 2 diabetes (T2D), blood pressure (BP), body mass index (BMI), cigarette smoking, and alcohol consumption with the risk of age-related cataract. Methods: To assess potential causal effects of clinical or behavioral factors on cataract risk, we conducted two-sample Mendelian randomization analyses. Genetic instruments, based on common genetic variants associated with risk factors at genome-wide significance (P < 5 × 10-8), were derived from published genome-wide association studies (GWAS). For age-related cataract, we used GWAS summary statistics from our previous GWAS conducted in the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort (28,092 cataract cases and 50,487 controls; all non-Hispanic whites) or in the UK Biobank (31,852 cataract cases and 428,084 controls; all European-descent individuals). We used the inverse-variance weighted (IVW) method as our primary source of Mendelian randomization estimates and conducted common sensitivity analyses. Results: We found that genetically determined POAG and mean spherical equivalent RE were significantly associated with cataract risk (IVW model: odds ratio [OR] = 1.04; 95% confidence interval [CI], 1.01-1.08; P = 0.018; per diopter more hyperopic: OR = 0.92; 95% CI, 0.89-0.93; P = 6.51 × 10-13, respectively). In contrast, genetically determined T2D, BP, BMI, cigarette smoking, or alcohol consumption were not associated with cataract risk (P > 0.05). Conclusions: Our results provide evidence that genetic risks for POAG and myopia may be causal risk factors for age-related cataract. These results are consistent with previous observational studies reporting associations of myopia with cataract risk. This information may support population cataract risk stratification and screening strategies.

Genetic Associations Between Smoking- and Glaucoma-Related Traits.

Purpose: The purpose of this study was to describe the genetic relationship between smoking and glaucoma. Methods: We used summary-level genetic data for smoking initiation, smoking intensity (cigarettes per day [CPD]), intraocular pressure (IOP), vertical cup-disc ratio, and open-angle glaucoma (OAG) to estimate global genetic correlations (rg) and perform two-sample Mendelian randomization (MR) experiments that explored relations between traits. Finally, we examined associations between smoking genetic risk scores (GRS) and smoking traits with measured IOP and OAG in Rotterdam Study participants. Results: We identified weak inverse rg between smoking- and glaucoma-related traits that were insignificant after Bonferroni correction. However, MR analysis revealed that genetically predicted smoking initiation was associated with lower IOP (-0.18 mm Hg per SD, 95% confidence interval [CI] = -0.30 to -0.06, P = 0.003). Furthermore, genetically predicted smoking intensity was associated with decreased OAG risk (odds ratio [OR] = 0.74 per SD, 95% CI = 0.61 to 0.90, P = 0.002). In the Rotterdam Study, the smoking initiation GRS was associated with lower IOP (-0.09 mm Hg per SD, 95% CI = -0.17 to -0.01, P = 0.04) and lower odds of OAG (OR = 0.84 per SD, 95% CI = 0.73 to 0.98, P = 0.02) in multivariable-adjusted analyses. In contrast, neither smoking history nor CPD was associated with IOP (P ≥ 0.38) or OAG (P ≥ 0.54). Associations between the smoking intensity GRS and glaucoma traits were null (P ≥ 0.13). Conclusions: MR experiments and GRS generated from Rotterdam Study participants support an inverse relationship between smoking and glaucoma. Translational Relevance: Understanding the genetic drivers of the inverse relationship between smoking and glaucoma could yield new insights into glaucoma pathophysiology.

One-year surgical outcomes of the PreserFlo MicroShunt in glaucoma: a multicentre analysis.

BACKGROUND/AIMS: To evaluate the efficacy and safety of the PreserFlo MicroShunt glaucoma device in a multicentre cohort study. METHODS: All consecutive patients who received the microshunt with mitomycin-C (MMC) 0.4 mg/mL from May 2019 to September 2020 in three UK tertiary centres. Primary outcome at 1 year was a complete success, with failure defined as intraocular pressure (IOP) >21 mmHg or <20% reduction, IOP≤5 mmHg with any decreased vision on two consecutive visits, reoperation or loss of light perception vision. Secondary outcomes were IOP, best-corrected visual acuity, medications, complications, interventions and reoperations. We also performed subgroup analyses for severe glaucoma and assessed risk factors for failure. RESULTS: 104 eyes had 1-year follow-up. Complete and qualified success at 1 year were achieved in 51.9% (N=54) and 16.4% (N=17), respectively, and failure occurred in 31.7% (N=33). There was a significant reduction in IOP (mmHg) from preoperatively (23.4±0.8, N=104) to 12 months (14.7±0.6, N=104) (p<0.0001). Antiglaucoma medications also decreased from preoperatively (3.4±0.1, N=104) to 12 months (0.7±0.1, N=104) (p<0.0001). Multivariate analyses showed an association between higher mean deviation and failure (HR 1.055, 95% CI 1.0075 to 1.11, p=0.0227). Complications were hypotony (19.2%; N=20), choroidal detachments (10.6%; N=11), hyphaema (5.8%; N=6) and bleb leak (5.8%; N=6). Needling and 5-fluorouracil injections were performed in 12.5% (N=13) and 33.7% (N=35), respectively, and 11.5% (N=12) required revision surgery. CONCLUSION: The PreserFlo MicroShunt with MMC 0.4 mg/mL showed an overall success rate of 68.3% at 1 year, and led to significant IOP and medication reduction with a low rate of adverse effects.

Genome-Wide Association Analysis of Over 170,000 Individuals from the UK Biobank Identifies Seven Loci Associated with Dietary Approaches to Stop Hypertension (DASH) Diet.

Diet is a modifiable risk factor for common chronic diseases and mental health disorders, and its effects are under partial genetic control. To estimate the impact of diet on individual health, most epidemiological and genetic studies have focused on individual aspects of dietary intake. However, analysing individual food groups in isolation does not capture the complexity of the whole diet pattern. Dietary indices enable a holistic estimation of diet and account for the intercorrelations between food and nutrients. In this study we performed the first ever genome-wide association study (GWA) including 173,701 individuals from the UK Biobank to identify genetic variants associated with the Dietary Approaches to Stop Hypertension (DASH) diet. DASH was calculated using the 24 h-recall questionnaire collected by UK Biobank. The GWA was performed using a linear mixed model implemented in BOLT-LMM. We identified seven independent single-nucleotide polymorphisms (SNPs) associated with DASH. Significant genetic correlations were observed between DASH and several educational traits with a significant enrichment for genes involved in the AMP-dependent protein kinase (AMPK) activation that controls the appetite by regulating the signalling in the hypothalamus. The colocalization analysis implicates genes involved in body mass index (BMI)/obesity and neuroticism (ARPP21, RP11-62H7.2, MFHAS1, RHEBL1). The Mendelian randomisation analysis suggested that increased DASH score, which reflect a healthy diet style, is causal of lower glucose, and insulin levels. These findings further our knowledge of the pathways underlying the relationship between diet and health outcomes. They may have significant implications for global public health and provide future dietary recommendations for the prevention of common chronic diseases.

Personalized Model to Predict Keratoconus Progression From Demographic, Topographic, and Genetic Data.

PURPOSE: To generate a prognostic model to predict keratoconus progression to corneal crosslinking (CXL). DESIGN: Retrospective cohort study. METHODS: We recruited 5025 patients (9341 eyes) with early keratoconus between January 2011 and November 2020. Genetic data from 926 patients were available. We investigated both keratometry or CXL as end points for progression and used the Royston-Parmar method on the proportional hazards scale to generate a prognostic model. We calculated hazard ratios (HRs) for each significant covariate, with explained variation and discrimination, and performed internal-external cross validation by geographic regions. RESULTS: After exclusions, model fitting comprised 8701 eyes, of which 3232 underwent CXL. For early keratoconus, CXL provided a more robust prognostic model than keratometric progression. The final model explained 33% of the variation in time to event: age HR (95% CI) 0.9 (0.90-0.91), maximum anterior keratometry 1.08 (1.07-1.09), and minimum corneal thickness 0.95 (0.93-0.96) as significant covariates. Single-nucleotide polymorphisms (SNPs) associated with keratoconus (n=28) did not significantly contribute to the model. The predicted time-to-event curves closely followed the observed curves during internal-external validation. Differences in discrimination between geographic regions was low, suggesting the model maintained its predictive ability. CONCLUSIONS: A prognostic model to predict keratoconus progression could aid patient empowerment, triage, and service provision. Age at presentation is the most significant predictor of progression risk. Candidate SNPs associated with keratoconus do not contribute to progression risk.

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.

Population screening for glaucoma in UK: current recommendations and future directions.

Effective population screening for glaucoma would enable earlier diagnosis and prevention of irreversible vision loss. The UK National Screening Committee (NSC) recently published a review that examined the viability, effectiveness and appropriateness of a population-based screening programme for primary open-angle glaucoma (POAG). In our article, we summarise the results of the review and discuss some future directions that may enable effective population screening for glaucoma in the future. Two key questions were addressed by the UK NSC review; is there a valid, accurate screening test for POAG, and does evidence exist that screening reduces morbidity from POAG compared with standard care. Six new studies were identified since the previous 2015 review. The review concluded that screening for glaucoma in adults is not recommended because there is no clear evidence for a sufficiently accurate screening test or for better outcomes with screening compared to current care. The next UK NSC review is due to be conducted in 2023. One challenge for POAG screening is that the relatively low disease prevalence results in too many false-positive referrals, even with an accurate test. In the future, targeted screening of a population subset with a higher prevalence of glaucoma may be effective. Recent developments in POAG polygenic risk prediction and deep learning image analysis offer potential avenues to identifying glaucoma-enriched sub-populations. Until such time, opportunistic case finding through General Ophthalmic Services remains the primary route for identification of glaucoma in the UK and greater public awareness of the service would be of benefit.

A large multiethnic GWAS meta-analysis of cataract identifies new risk loci and sex-specific effects.

Cataract is the leading cause of blindness among the elderly worldwide and cataract surgery is one of the most common operations performed in the United States. As the genetic etiology of cataract formation remains unclear, we conducted a multiethnic genome-wide association meta-analysis, combining results from the GERA and UK Biobank cohorts, and tested for replication in the 23andMe research cohort. We report 54 genome-wide significant loci, 37 of which were novel. Sex-stratified analyses identified CASP7 as an additional novel locus specific to women. We show that genes within or near 80% of the cataract-associated loci are significantly expressed and/or enriched-expressed in the mouse lens across various spatiotemporal stages as per iSyTE analysis. Furthermore, iSyTE shows 32 candidate genes in the associated loci have altered gene expression in 9 different gene perturbation mouse models of lens defects/cataract, suggesting their relevance to lens biology. Our work provides further insight into the complex genetic architecture of cataract susceptibility.

The genomic loci of specific human tRNA genes exhibit ageing-related DNA hypermethylation.

The epigenome has been shown to deteriorate with age, potentially impacting on ageing-related disease. tRNA, while arising from only ˜46 kb (<0.002% genome), is the second most abundant cellular transcript. tRNAs also control metabolic processes known to affect ageing, through core translational and additional regulatory roles. Here, we interrogate the DNA methylation state of the genomic loci of human tRNA. We identify a genomic enrichment for age-related DNA hypermethylation at tRNA loci. Analysis in 4,350 MeDIP-seq peripheral-blood DNA methylomes (16-82 years), identifies 44 and 21 hypermethylating specific tRNAs at study-and genome-wide significance, respectively, contrasting with none hypomethylating. Validation and replication (450k array and independent targeted Bisuphite-sequencing) supported the hypermethylation of this functional unit. Tissue-specificity is a significant driver, although the strongest consistent signals, also independent of major cell-type change, occur in tRNA-iMet-CAT-1-4 and tRNA-Ser-AGA-2-6. This study presents a comprehensive evaluation of the genomic DNA methylation state of human tRNA genes and reveals a discreet hypermethylation with advancing age.

Association Between Medication-Taking and Refractive Error in a Large General Population-Based Cohort.

Purpose: Refractive errors, particularly myopia, are common and a leading cause of blindness. This study aimed to explore associations between medications and refractive error in an aging adult cohort and to determine whether childhood-onset refractive errors predict future medication use to provide novel insights into disease mechanisms. Methods: The study compared the spherical equivalent values measured in 102,318 UK Biobank participants taking the 960 most commonly used medications. The strengths of associations were evaluated against the self-reported age of spectacle wear. The causality of refractive error changes was inferred using sensitivity and Mendelian randomization analyses. Results: Anti-glaucoma drugs were associated with 1 to 2 diopters greater myopic refraction, particularly in subjects who started wearing correction in the first two decades of life, potentially due to the association of higher intraocular pressure since early years with both myopia and, later in life, glaucoma. All classes of pain-control medications, including paracetamol, opiates, non-steroidal antiinflammatory drugs, and gabapentinoids, were associated with greater hyperopia (+0.68-1.15 diopters), after correction for deprivation, education, and polypharmacy and sensitivity analyses for common diagnoses. Oral hypoglycemics (metformin, gliburonide) were associated with myopia, as was allopurinol, and participants using bronchodilators (ipratropium and salbutamol) were more hyperopic. Conclusions: This study finds for the first time, to our knowledge, that medication use is associated with refractive error in adults. The novel finding that analgesics are associated with hyperopic refraction, and the possibility that multisite chronic pain predisposes to hyperopia, deserves further research. Some drugs, such as antihyperglycemic or bronchodilators, may directly alter refractive error. Intraocular pressure appears causative for myopia.

Looking for Sunshine: Genetic Predisposition to Sun Seeking in 265,000 Individuals of European Ancestry.

Despite growing public awareness of the adverse consequences of excessive sun exposure, modifying sun-seeking behavior is challenging because it appears to be driven by addictive mechanisms. This can have effects on health because sun exposure, although beneficial, when prolonged and repeated shows a causal relationship with skin cancer risk. Using data from 2,500 United Kingdom twins, we observed sun seeking to be significantly heritable (h2 ≥ 58%). In a GWAS meta-analysis of sun-seeking behavior in 261,915 subjects of European ancestry, we identified five GWAS-significant loci previously associated with addiction, behavioral and personality traits, cognitive function, and educational attainment and enriched for CNS gene expression: MIR2113 (P = 2.08 × 10-11), FAM76B/MTMR2/CEP57 (P = 3.70 × 10-9), CADM2 (P = 9.36 × 10-9), TMEM182 (P = 1.64 × 10-8), and PLCL1/LINC01923/SATB2 (P = 3.93 × 10-8). These findings imply that the behavior concerning UV exposure is complicated by a genetic predisposition shared with neuropsychological traits. This should be taken into consideration when designing awareness campaigns and may help improve people's attitudes toward sun exposure.

Variants associated with HHIP expression have sex-differential effects on lung function.

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P<5x10 -8) interactions with sex on lung function, and 21 showed suggestive interactions (P<1x10 -6). The strongest signal, from rs7697189 (chr4:145436894) on forced expiratory volume in 1 second (FEV 1) (P=3.15x10 -15), was replicated (P=0.016) in SpiroMeta. The C allele increased FEV 1 more in males (untransformed FEV 1 ß=0.028 [SE 0.0022] litres) than females (ß=0.009 [SE 0.0014] litres), and this effect was not accounted for by differential effects on height, smoking or pubertal age. rs7697189 resides upstream of the hedgehog-interacting protein ( HHIP) gene and was previously associated with lung function and HHIP lung expression. We found HHIP expression was significantly different between the sexes (P=6.90x10 -6), but we could not detect sex differential effects of rs7697189 on expression. Conclusions: We identified a novel genotype-by-sex interaction at a putative enhancer region upstream of the HHIP gene. Establishing the mechanism by which HHIP SNPs have different effects on lung function in males and females will be important for our understanding of lung health and diseases in both sexes.

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.

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.

Genetic Variants Associated With Corneal Biomechanical Properties and Potentially Conferring Susceptibility to Keratoconus in a Genome-Wide Association Study.

IMPORTANCE: Keratoconus is an important cause of visual loss in young adults, but little is known about its genetic causes. Understanding the genetic determinants of corneal biomechanical factors may in turn teach us about keratoconus etiology. OBJECTIVES: To identify genetic associations with corneal biomechanical properties and to examine whether these genetic variants are associated with keratoconus. DESIGN, SETTING, AND PARTICIPANTS: A stage 1 discovery and replication genome-wide association study (GWAS) of corneal biomechanical properties was performed in 2 cross-sectional populations (6645 participants from the European Prospective Investigation into Cancer and Nutrition [EPIC]-Norfolk Eye Study and 2384 participants from the TwinsUK study). In stage 2, the association of genetic determinants identified in stage 1 with keratoconus was examined in a case-control study. A total of 752 patients with keratoconus were compared with 974 TwinsUK participants (undergoing direct sequencing) or 13 828 EPIC-Norfolk participants (undergoing genotyping and imputation) who were not part of the stage 1 analysis. Data were collected from March 1, 1993, through March 13, 2017, and analyzed from November 1, 2015, through February 1, 2018. EXPOSURES: In stage 1, allele dosage at genome-wide single-nucleotide polymorphisms (SNPs); in stage 2, allele dosage at SNPs with genome-wide significance (P < 5 × 10-8) in stage 1 and not previously reported as associated with corneal disease. MAIN OUTCOMES AND MEASURES: In stage 1, corneal hysteresis (CH) and corneal resistance factor (CRF), measured with the Ocular Response Analyzer (ORA); in stage 2, association with keratoconus compared with controls. RESULTS: Among 6645 participants in the discovery cohort (3635 women (54.7%); mean age, 69 years [range, 48-92 years]), 7 genome-wide significant loci associated with CH or CRF were identified that were independently replicated. Two further suggestive loci were identified after meta-analysis. To date, 5 of the identified loci, at ANAPC1, ADAMTS8, ADAMTS17, ABCA6, and COL6A1, have not previously been reported as associated with corneal disease. The ABCA6 locus (rs77542162) was associated with keratoconus using the TwinsUK (odds ratio [OR], 0.50; 95% CI, 0.27-0.92; P = .03) and EPIC-Norfolk controls (OR, 0.39; 95% CI, 0.22-0.70; P = .002). The other loci were associated with keratoconus using TwinsUK (OR per effect allele for ADAMTS8, 0.51 [95% CI, 0.37-0.71; P = 7.9 × 10-5]; for COL6A1, 1.65 [95% CI, 1.05-2.59; P = .03]) or EPIC-Norfolk (OR per effect allele for ANAPC1, 0.78 [95% CI, 0.68-0.89; P = 3.7 × 10-4]; for ADAMTS17, 0.82 [95% CI, 0.68-0.99; P = .04]) controls. CONCLUSIONS AND RELEVANCE: Five loci that are associated with corneal biomechanical properties and that have suggestive associations with keratoconus were reported. These findings suggest the role of type VI collagen, extracellular matrix, and connective-tissue development for corneal biomechanics and keratoconus and the role of CH and CRF as biomarkers for keratoconus.

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.