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.

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.

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.

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.

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.

Genome-wide association study in 176,678 Europeans reveals genetic loci for tanning response to sun exposure.

The skin's tendency to sunburn rather than tan is a major risk factor for skin cancer. Here we report a large genome-wide association study of ease of skin tanning in 176,678 subjects of European ancestry. We identify significant association with tanning ability at 20 loci. We confirm previously identified associations at six of these loci, and report 14 novel loci, of which ten have never been associated with pigmentation-related phenotypes. Our results also suggest that variants at the AHR/AGR3 locus, previously associated with cutaneous malignant melanoma the underlying mechanism of which is poorly understood, might act on disease risk through modulation of tanning ability.

Changes in quality of life shortly after routine cataract surgery.

OBJECTIVE: To perform a pilot study to explore use of a generic patient-reported outcome measure to assess patient-perceived improvements in quality of life within 2-4 weeks of routine cataract surgery and to explore differences after first or second eye surgery. Secondary analysis explored effects of gender and ethnicity. DESIGN: Prospective observational study. PARTICIPANTS: Consecutive patients attending a weekly nurse-led postoperative clinic. METHODS: The Glasgow Benefit Inventory (GBI), a validated, postinterventional questionnaire (not specific to one particular medical or surgical intervention), was administered. Mean scores were calculated. Scores were compared when patients were grouped by first or second eye, and by gender or ethnicity (unpaired t test). Scores potentially range from +100 (maximum benefit) to -100 (maximum detriment). RESULTS: The GBI was administered 113 times to 109 patients (4 patients were seen after both first and second eye surgery). Mean overall score was +22.8 (median +19.4; SD 19.7; 95% CI +19.2 to +26.4). Mean (SD) subscores were +30.5 (25.3), +17.8 (26.7), and -3.1 (19.9) for general, social support, and physical health subdomains, respectively. Total benefit scores were not significantly different for first or second eye surgery, or across gender (p > 0.3). Scores for patients of African (including African Caribbean) ethnicity were significantly higher than those obtained from European patients (p = 0.002). CONCLUSIONS: Patients reported significant improvements in quality of life even a few weeks after cataract surgery, as assessed by the GBI. Second eye surgery appeared to confer similar benefit to first eye surgery. The significant difference in scores between ethnic groups invites further investigation.

Integrative pathway genomics of lung function and airflow obstruction.

Chronic respiratory disorders are important contributors to the global burden of disease. Genome-wide association studies (GWASs) of lung function measures have identified several trait-associated loci, but explain only a modest portion of the phenotypic variability. We postulated that integrating pathway-based methods with GWASs of pulmonary function and airflow obstruction would identify a broader repertoire of genes and processes influencing these traits. We performed two independent GWASs of lung function and applied gene set enrichment analysis to one of the studies and validated the results using the second GWAS. We identified 131 significantly enriched gene sets associated with lung function and clustered them into larger biological modules involved in diverse processes including development, immunity, cell signaling, proliferation and arachidonic acid. We found that enrichment of gene sets was not driven by GWAS-significant variants or loci, but instead by those with less stringent association P-values. Next, we applied pathway enrichment analysis to a meta-analyzed GWAS of airflow obstruction. We identified several biologic modules that functionally overlapped with those associated with pulmonary function. However, differences were also noted, including enrichment of extracellular matrix (ECM) processes specifically in the airflow obstruction study. Network analysis of the ECM module implicated a candidate gene, matrix metalloproteinase 10 (MMP10), as a putative disease target. We used a knockout mouse model to functionally validate MMP10's role in influencing lung's susceptibility to cigarette smoke-induced emphysema. By integrating pathway analysis with population-based genomics, we unraveled biologic processes underlying pulmonary function traits and identified a candidate gene for obstructive lung disease.

Molecular mechanisms underlying variations in lung function: a systems genetics analysis.

BACKGROUND: Lung function measures reflect the physiological state of the lung, and are essential to the diagnosis of chronic obstructive pulmonary disease (COPD). The SpiroMeta-CHARGE consortium undertook the largest genome-wide association study (GWAS) so far (n=48,201) for forced expiratory volume in 1 s (FEV1) and the ratio of FEV1 to forced vital capacity (FEV1/FVC) in the general population. The lung expression quantitative trait loci (eQTLs) study mapped the genetic architecture of gene expression in lung tissue from 1111 individuals. We used a systems genetics approach to identify single nucleotide polymorphisms (SNPs) associated with lung function that act as eQTLs and change the level of expression of their target genes in lung tissue; termed eSNPs. METHODS: The SpiroMeta-CHARGE GWAS results were integrated with lung eQTLs to map eSNPs and the genes and pathways underlying the associations in lung tissue. For comparison, a similar analysis was done in peripheral blood. The lung mRNA expression levels of the eSNP-regulated genes were tested for associations with lung function measures in 727 individuals. Additional analyses identified the pleiotropic effects of eSNPs from the published GWAS catalogue, and mapped enrichment in regulatory regions from the ENCODE project. Finally, the Connectivity Map database was used to identify potential therapeutics in silico that could reverse the COPD lung tissue gene signature. FINDINGS: SNPs associated with lung function measures were more likely to be eQTLs and vice versa. The integration mapped the specific genes underlying the GWAS signals in lung tissue. The eSNP-regulated genes were enriched for developmental and inflammatory pathways; by comparison, SNPs associated with lung function that were eQTLs in blood, but not in lung, were only involved in inflammatory pathways. Lung function eSNPs were enriched for regulatory elements and were over-represented among genes showing differential expression during fetal lung development. An mRNA gene expression signature for COPD was identified in lung tissue and compared with the Connectivity Map. This in-silico drug repurposing approach suggested several compounds that reverse the COPD gene expression signature, including a nicotine receptor antagonist. These findings represent novel therapeutic pathways for COPD. INTERPRETATION: The system genetics approach identified lung tissue genes driving the variation in lung function and susceptibility to COPD. The identification of these genes and the pathways in which they are enriched is essential to understand the pathophysiology of airway obstruction and to identify novel therapeutic targets and biomarkers for COPD, including drugs that reverse the COPD gene signature in silico. FUNDING: The research reported in this article was not specifically funded by any agency. See Acknowledgments for a full list of funders of the lung eQTL study and the Spiro-Meta CHARGE GWAS.

Genetics of skin color variation in Europeans: genome-wide association studies with functional follow-up.

In the International Visible Trait Genetics (VisiGen) Consortium, we investigated the genetics of human skin color by combining a series of genome-wide association studies (GWAS) in a total of 17,262 Europeans with functional follow-up of discovered loci. Our GWAS provide the first genome-wide significant evidence for chromosome 20q11.22 harboring the ASIP gene being explicitly associated with skin color in Europeans. In addition, genomic loci at 5p13.2 (SLC45A2), 6p25.3 (IRF4), 15q13.1 (HERC2/OCA2), and 16q24.3 (MC1R) were confirmed to be involved in skin coloration in Europeans. In follow-up gene expression and regulation studies of 22 genes in 20q11.22, we highlighted two novel genes EIF2S2 and GSS, serving as competing functional candidates in this region and providing future research lines. A genetically inferred skin color score obtained from the 9 top-associated SNPs from 9 genes in 940 worldwide samples (HGDP-CEPH) showed a clear gradual pattern in Western Eurasians similar to the distribution of physical skin color, suggesting the used 9 SNPs as suitable markers for DNA prediction of skin color in Europeans and neighboring populations, relevant in future forensic and anthropological investigations.

Prevalence and risk factors of dry eye disease in a British female cohort.

BACKGROUND/AIMS: To estimate the prevalence and risk factors of dry eye disease (DED) in a female cohort in the UK. METHODS: Population-based cross-sectional association study of 3824 women from the TwinsUK cohort aged 20-87 years. A questionnaire was used to evaluate DED and several risk factors. Binary logistic regression, corrected for age, was used to examine the association between DED and risk factors. RESULTS: 9.6% of women had a DED diagnosis and concomitant use of artificial tears, and 20.8% experienced DED symptoms in the past 3 months. Risk factors that were significantly associated with DED were age, asthma, eczema, the presence of any allergy, cataract surgery, rheumatoid arthritis, osteoarthritis, migraine and stroke. The highest effect sizes were found with depression, pelvic pain, irritable bowel syndrome and chronic widespread pain syndrome (all p<0.0005). Subjects with DED symptoms scored significantly lower on self-perceived health, compared with controls (p=0.001). CONCLUSIONS: DED is common and increases with age within this cohort of female twins. We confirmed established risk factors for the first time in a British population, and found important risk factors that might relate to an underlying aetiology involving chronic pain predisposition or somatisation.