Genome-wide association analyses identify 139 loci associated with macular thickness in the UK Biobank cohort.

The macula, located near the center of the retina in the human eye, is responsible for providing critical functions, such as central, sharp vision. Structural changes in the macula are associated with many ocular diseases, including age-related macular degeneration (AMD) and glaucoma. Although macular thickness is a highly heritable trait, there are no prior reported genome-wide association studies (GWASs) of it. Here we describe the first GWAS of macular thickness, which was measured by spectral-domain optical coherence tomography using 68 423 participants from the UK Biobank cohort. We identified 139 genetic loci associated with macular thickness at genome-wide significance (P < 5 × 10-8). The most significant loci were LINC00461 (P = 5.1 × 10-120), TSPAN10 (P = 1.2 × 10-118), RDH5 (P = 9.2 × 10-105) and SLC6A20 (P = 1.4 × 10-71). Results from gene expression demonstrated that these genes are highly expressed in the retina. Other hits included many previously reported AMD genes, such as NPLOC4 (P = 1.7 × 10-103), RAD51B (P = 9.1 × 10-14) and SLC16A8 (P = 1.7 × 10-8), further providing functional significance of the identified loci. Through cross-phenotype analysis, these genetic loci also exhibited pleiotropic effects with myopia, neurodegenerative diseases (e.g. Parkinson's disease, schizophrenia and Alzheimer's disease), cancer (e.g. breast, ovarian and lung cancers) and metabolic traits (e.g. body mass index, waist circumference and type 2 diabetes). Our findings provide the first insight into the genetic architecture of macular thickness and may further elucidate the pathogenesis of related ocular diseases, such as AMD.

Genome-wide association analyses identify new loci influencing intraocular pressure.

Elevated intraocular pressure (IOP) is a significant risk factor for glaucoma, the leading cause of irreversible blindness worldwide. While previous studies have identified numerous genetic variants associated with IOP, these loci only explain a fraction of IOP heritability. Recently established of biobank repositories have resulted in large amounts of data, enabling the identification of the remaining heritability for complex traits. Here, we describe the largest genome-wide association study of IOP to date using participants of European ancestry from the UK Biobank. We identified 671 directly genotyped variants that are significantly associated with IOP (P < 5 × 10-8). In addition to 103 novel loci, the top ranked novel IOP genes are LMX1B, NR1H3, MADD and SEPT9. We replicated these findings in an external population and examined the pleiotropic nature of these loci. These discoveries not only further our understanding of the genetic architecture of IOP, but also shed new light on the biological processes underlying glaucoma.

PleioNet: a web-based visualization tool for exploring pleiotropy across complex traits.

SUMMARY: Pleiotropy plays an important role in furthering our understanding of the shared genetic architecture of different human diseases and traits. However, exploring and visualizing pleiotropic information with currently publicly available tools is limiting and challenging. To aid researchers in constructing and digesting pleiotropic networks, we present PleioNet, a web-based visualization tool for exploring this information across human diseases and traits. This program provides an intuitive and interactive web interface that seamlessly integrates large database queries with visualizations that enable users to quickly explore complex high-dimensional pleiotropic information. PleioNet works on all modern computer and mobile web browsers, making pleiotropic information readily available to a broad range of researchers and clinicians with diverse technical backgrounds. We expect that PleioNet will be an important tool for studying the underlying pleiotropic connections among human diseases and traits. AVAILABILITY AND IMPLEMENTATION: PleioNet is hosted on Google cloud and freely available at http://www.pleionet.com/.

Polygenic Risk Score Is Associated With Intraocular Pressure and Improves Glaucoma Prediction in the UK Biobank Cohort.

PURPOSE: Elevated intraocular pressure (IOP) is an important risk factor for glaucoma. We constructed polygenic risk scores (PRSs) for IOP using the UK Biobank (UKB) data set to test whether the PRSs are associated with IOP and whether using them improves glaucoma prediction. METHODS: We conducted this study using 435,678 European participants from the UKB. We constructed weighted and unweighted PRSs using single nucleotide polymorphisms (SNPs) derived from the UKB data and previously reported IOP SNPs. We examined the associations of the PRSs with IOP and primary open-angle glaucoma (POAG) using linear and logistic regression, respectively. To quantify the discriminatory ability of the PRSs on POAG, we used the area under the receiver operating characteristic curve (AUC). RESULTS: The weighted PRS was significantly associated with IOP (P ∼ 10-200), after adjusting for age and sex. The PRS explained an additional 4% of variance in IOP. The weighted PRS was also significantly associated with POAG (P = 1.8 × 10-77). Subjects in the top quintile of the IOP PRS were 6.34 (95% confidence interval [CI]: 4.82-8.33; P = 2.1 × 10-57) times more likely to have POAG, compared to those in the bottom category. The weighted PRS improved the discriminatory power for POAG (AUC increased by 5%, P = 6.2 × 10-22) when added to the other covariates. The unweighted PRS exhibited similar results. CONCLUSIONS: We determined that IOP PRSs are significantly associated with IOP and improve the prediction of POAG. TRANSLATIONAL RELEVANCE: PRSs help reduce the burden of glaucoma by early detection of genetically susceptible individuals.

The Effect of Micro- and Nanoscale Surface Topographies on Silk on Human Corneal Limbal Epithelial Cell Differentiation.

We previously reported that micro- and nano-scale topographic pitch created on silk films mimic features of the corneal basement membrane by providing biophysical cues to direct corneal epithelial cell adherence and migration. However, the effect of these topographical features on corneal limbal epithelial cell differentiation has not been explored. We hypothesize in the current study that various topographical pitch created on silk may affect corneal epithelial stem cell differentiation and alter the expression of genes involved in cell differentiation and self-renewal. We patterned silk films with different topographic pitch via soft lithography and observed human corneal limbal epithelial cell behavior. Colony forming assay demonstrated increased colony forming efficiency on patterned silk films. Cells cultured on nanoscale patterned silk films also expressed lower levels of putative keratocyte differentiation markers and higher levels of putative limbal stem cell markers. RNA-Seq analysis further implicated the involvement of pathways related to stem cell differentiation and self-renewal, including Notch, ERK/MAPK and Wnt/ß-catenin signaling. We conclude that patterned silk film substrates can be used as scaffolds and provide biophysical cues to corneal limbal stem cells that may maintain corneal epithelial stem cells at a less differentiated state.

Differential Expression of Coding and Long Noncoding RNAs in Keratoconus-Affected Corneas.

Purpose: Keratoconus (KC) is the most common corneal ectasia. We aimed to determine the differential expression of coding and long noncoding RNAs (lncRNAs) in human corneas affected with KC. Methods: From the corneas of 10 KC patients and 8 non-KC healthy controls, 200 ng total RNA was used to prepare sequencing libraries with the SMARTer Stranded RNA-Seq kit after ribosomal RNA depletion, followed by paired-end 50-bp sequencing with Illumina Sequencer. Differential analysis was done using TopHat/Cufflinks with a gene file from Ensembl and a lncRNA file from NONCODE. Pathway analysis was performed using WebGestalt. Using the expression level of differentially expressed coding and noncoding RNAs in each sample, we correlated their expression levels in KC and controls separately and identified significantly different correlations in KC against controls followed by visualization using Cytoscape. Results: Using |fold change| ≥ 2 and a false discovery rate ≤ 0.05, we identified 436 coding RNAs and 584 lncRNAs with differential expression in the KC-affected corneas. Pathway analysis indicated the enrichment of genes involved in extracellular matrix, protein binding, glycosaminoglycan binding, and cell migration. Our correlation analysis identified 296 pairs of significant KC-specific correlations containing 117 coding genes enriched in functions related to cell migration/motility, extracellular space, cytokine response, and cell adhesion. Our study highlighted the potential roles of several genes (CTGF, SFRP1, AQP5, lnc-WNT4-2:1, and lnc-ALDH3A2-2:1) and pathways (TGF-ß, WNT signaling, and PI3K/AKT pathways) in KC pathogenesis. Conclusions: Our RNA-Seq-based differential expression and correlation analyses have identified many potential KC contributing coding and noncoding RNAs.

Family-Based Genome-Wide Association Study of South Indian Pedigrees Supports WNT7B as a Central Corneal Thickness Locus.

Purpose: To identify genetic risk factors contributing to central corneal thickness (CCT) in individuals from South India, a population with a high prevalence of ocular disorders. Methods: One hundred ninety-five individuals from 15 large South Indian pedigrees were genotyped using the Omni2.5 bead array. Family-based association for CCT was conducted using the score test in MERLIN. Results: Genome-wide association study (GWAS) identified strongest association for single nucleotide polymorphisms (SNPs) in the first intron of WNT7B and CCT (top SNP rs9330813; ß = -0.57, 95% confidence interval [CI]: -0.78 to -0.36; P = 1.7 × 10-7). We further investigated rs9330813 in a Latino cohort and four independent European cohorts. A meta-analysis of these data sets demonstrated statistically significant association between rs9330813 and CCT (ß = -3.94, 95% CI: -5.23 to -2.66; P = 1.7 × 10-9). WNT7B SNPs located in the same genomic region that includes rs9330813 have previously been associated with CCT in Latinos but with other ocular quantitative traits related to myopia (corneal curvature and axial length) in a Japanese population (rs10453441 and rs200329677). To evaluate the specificity of the observed WNT7B association with CCT in the South Indian families, we completed an ocular phenome-wide association study (PheWAS) for the top WNT7B SNPs using 45 ocular traits measured in these same families including corneal curvature and axial length. The ocular PheWAS results indicate that in the South Indian families WNT7B SNPs are primarily associated with CCT. Conclusions: The results indicate robust evidence for association between WNT7B SNPs and CCT in South Indian pedigrees, and suggest that WNT7B SNPs can have population-specific effects on ocular quantitative traits.

Micro- and Nanoscale Topographies on Silk Regulate Gene Expression of Human Corneal Epithelial Cells.

Purpose: Corneal basement membrane has topographical features that provide biophysical cues to direct cell adherence, migration, and proliferation. In this study, we hypothesize that varying topographic pitch created on silk films can alter epithelial cell morphology, adhesion, and the genetic expression involved in cytoskeletal dynamics-related pathways. Methods: Silicon wafers with parallel ridge widths of 2000, 1000, and 800 nm were produced and used to pattern silk films via soft lithography. Human corneal epithelial cells were cultured onto silk. After 72 hours of incubation, images were taken to study cell morphology and alignment. Cytoskeletal structures were studied by immunofluorescent staining. RNA was collected from cultured cells to perform RNA-Seq transcriptome analysis using the Illumina Hiseq 2500 sequencing system. Differentially expressed genes were identified using DNAstar Qseq then verified using quantitative real-time PCR. These genes were used to perform pathway analyses using Ingenuity Pathways Analysis. Results: Primary human corneal epithelial cell alignment to the surface pattern was the greatest on 1000-nm features. Fluorescent microscopy of f-actin staining showed cell cytoskeleton alignment either in parallel (2000 nm) or perpendicular (1000 and 800 nm) to the long feature axis. Z-stack projection of vinculin staining indicated increased focal adhesion formation localized on the cellular basal surface. RNA-seq analysis revealed differentially expressed genes involved in actin organization, integrin signaling, and focal adhesion kinase signaling (-log (P)>5). Conclusions: Patterned silk film substrates may serve as a scaffold and provide biophysical cues to corneal epithelial cells that change their gene expression, alter cellular adherence, morphology, and may offer a promising customizable material for use in ocular surface repair.