Exome array analysis identifies CAV1/CAV2 as a susceptibility locus for intraocular pressure

Deep Learning-Based Identification of Intraocular Pressure-Associated Genes Influencing Trabecular Meshwork Cell Morphology

Purpose

Genome-wide association studies have recently uncovered many loci associated with variation in intraocular pressure (IOP). Artificial intelligence (AI) can be used to interrogate the effect of specific genetic knockouts on the morphology of trabecular meshwork cells (TMCs) and thus, IOP regulation.

Design

Experimental study.

Subjects

Primary TMCs collected from human donors.

Methods

Sixty-two genes at 55 loci associated with IOP variation were knocked out in primary TMC lines. All cells underwent high-throughput microscopy imaging after being stained with a 5-channel fluorescent cell staining protocol. A convolutional neural network was trained to distinguish between gene knockout and normal control cell images. The area under the receiver operator curve (AUC) metric was used to quantify morphological variation in gene knockouts to identify potential pathological perturbations.

Main Outcome Measures

Degree of morphological variation as measured by deep learning algorithm accuracy of differentiation from normal controls.

Results

Cells where LTBP2 or BCAS3 had been perturbed demonstrated the greatest morphological variation from normal TMCs (AUC 0.851, standard deviation [SD] 0.030; and AUC 0.845, SD 0.020, respectively). Of 7 multigene loci, 5 had statistically significant differences in AUC (P < 0.05) between genes, allowing for pathological gene prioritization. The mitochondrial channel most frequently showed the greatest degree of morphological variation (33.9% of cell lines).

Conclusions

We demonstrate a robust method for functionally interrogating genome-wide association signals using high-throughput microscopy and AI. Genetic variations inducing marked morphological variation can be readily identified, allowing for the gene-based dissection of loci associated with complex traits.

Financial Disclosures

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Molecular Mechanisms of Glaucoma Pathogenesis with Implications to Caveolin Adaptor Protein and Caveolin-Shp2 Axis

Glaucoma is a common retinal disorder characterized by progressive optic nerve damage, resulting in visual impairment and potential blindness. Elevated intraocular pressure (IOP) is a major risk factor, but some patients still experience disease progression despite IOP-lowering treatments. Genome-wide association studies have linked variations in the Caveolin1/2 (CAV-1/2) gene loci to glaucoma risk. Cav-1, a key protein in caveolae membrane invaginations, is involved in signaling pathways and its absence impairs retinal function. Recent research suggests that Cav-1 is implicated in modulating the BDNF/TrkB signaling pathway in retinal ganglion cells, which plays a critical role in retinal ganglion cell (RGC) health and protection against apoptosis. Understanding the interplay between these proteins could shed light on glaucoma pathogenesis and provide potential therapeutic targets.

10 Years of GWAS in intraocular pressure

Intraocular pressure (IOP) is the only modifiable risk factor for glaucoma, the leading cause of irreversible blindness worldwide. In this review, we summarize the findings of genome-wide association studies (GWASs) of IOP published in the past 10 years and prior to December 2022. Over 190 genetic loci and candidate genes associated with IOP have been uncovered through GWASs, although most of these studies were conducted in subjects of European and Asian ancestries. We also discuss how these common variants have been used to derive polygenic risk scores for predicting IOP and glaucoma, and to infer causal relationship with other traits and conditions through Mendelian randomization. Additionally, we summarize the findings from a recent large-scale exome-wide association study (ExWAS) that identified rare variants associated with IOP in 40 novel genes, six of which are drug targets for clinical treatment or are being evaluated in clinical trials. Finally, we discuss the need for future genetic studies of IOP to include individuals from understudied populations, including Latinos and Africans, in order to fully characterize the genetic architecture of IOP.

Calcium-Signalling in Human Glaucoma Lamina Cribrosa Myofibroblasts

Glaucoma is one of the most common causes of treatable visual impairment in the developed world, affecting approximately 64 million people worldwide, some of whom will be bilaterally blind from irreversible optic nerve damage. The optic nerve head is a key site of damage in glaucoma where there is fibrosis of the connective tissue in the lamina cribrosa (LC) extracellular matrix. As a ubiquitous second messenger, calcium (Ca²⁺) can interact with various cellular proteins to regulate multiple physiological processes and contribute to a wide range of diseases, including cancer, fibrosis, and glaucoma. Our research has shown evidence of oxidative stress, mitochondrial dysfunction, an elevated expression of Ca²⁺ entry channels, Ca²⁺-dependent pumps and exchangers, and an abnormal rise in cytosolic Ca²⁺ in human glaucomatous LC fibroblast cells. We have evidence that this increase is dependent on Ca²⁺ entry channels located in the plasma membrane, and its release is from internal stores in the endoplasmic reticulum (ER), as well as from the mitochondria. Here, we summarize some of the molecular Ca²⁺-dependent mechanisms related to this abnormal Ca²⁺-signalling in human glaucoma LC cells, with a view toward identifying potential therapeutic targets for ongoing optic neuropathy.

Caveolar and non-Caveolar Caveolin-1 in ocular homeostasis and disease

Caveolae, specialized plasma membrane invaginations present in most cell types, play important roles in multiple cellular processes including cell signaling, lipid uptake and metabolism, endocytosis and mechanotransduction. They are found in almost all cell types but most abundant in endothelial cells, adipocytes and fibroblasts. Caveolin-1 (Cav1), the signature structural protein of caveolae was the first protein associated with caveolae, and in association with Cavin1/PTRF is required for caveolae formation. Genetic ablation of either Cav1 or Cavin1/PTRF downregulates expression of the other resulting in loss of caveolae. Studies using Cav1-deficient mouse models have implicated caveolae with human diseases such as cardiomyopathies, lipodystrophies, diabetes and muscular dystrophies. While caveolins and caveolae are extensively studied in extra-ocular settings, their contributions to ocular function and disease pathogenesis are just beginning to be appreciated. Several putative caveolin/caveolae functions are relevant to the eye and Cav1 is highly expressed in retinal vascular and choroidal endothelium, Müller glia, the retinal pigment epithelium (RPE), and the Schlemm's canal endothelium and trabecular meshwork cells. Variants at the CAV1/2 gene locus are associated with risk of primary open angle glaucoma and the high risk HTRA1 variant for age-related macular degeneration is thought to exert its effect through regulation of Cav1 expression. Caveolins also play important roles in modulating retinal neuroinflammation and blood retinal barrier permeability. In this article, we describe the current state of caveolin/caveolae research in the context of ocular function and pathophysiology. Finally, we discuss new evidence showing that retinal Cav1 exists and functions outside caveolae.

A Role of Caveolae in Trabecular Meshwork Mechanosensing and Contractile Tone

Polymorphisms in the CAV1/2 gene loci impart increased risk for primary open-angle glaucoma (POAG). CAV1 encodes caveolin-1 (Cav1), which is required for biosynthesis of plasma membrane invaginations called caveolae. Cav1 knockout mice exhibit elevated intraocular pressure (IOP) and decreased outflow facility, but the mechanistic role of Cav1 in IOP homeostasis is unknown. We hypothesized that caveolae sequester/inhibit RhoA, to regulate trabecular meshwork (TM) mechanosensing and contractile tone. Using phosphorylated myosin light chain (pMLC) as a surrogate indicator for Rho/ROCK activity and contractile tone, we found that pMLC was elevated in Cav1-deficient TM cells compared to control (131 ± 10%, n = 10, p = 0.016). Elevation of pMLC levels following Cav1 knockdown occurred in cells on a soft surface (137 ± 7%, n = 24, p < 0.0001), but not on a hard surface (122 ± 17%, n = 12, p = 0.22). In Cav1-deficient TM cells where pMLC was elevated, Rho activity was also increased (123 ± 7%, n = 6, p = 0.017), suggesting activation of the Rho/ROCK pathway. Cyclic stretch reduced pMLC/MLC levels in TM cells (69 ± 7% n = 9, p = 0.002) and in Cav1-deficient TM cells, although not significantly (77 ± 11% n = 10, p = 0.059). Treatment with the Cav1 scaffolding domain mimetic, cavtratin (1 μM) caused a reduction in pMLC (70 ± 5% n = 7, p = 0.001), as did treatment with the scaffolding domain mutant cavnoxin (1 μM) (82 ± 7% n = 7, p = 0.04). Data suggest that caveolae differentially regulate RhoA signaling, and that caveolae participate in TM mechanotransduction. Cav1 regulation of these key TM functions provide evidence for underlying mechanisms linking polymorphisms in the Cav1/2 gene loci with increased POAG risk.

Variants in clock genes could be associated with lower risk of type 2 diabetes in an elderly Greek population

OBJECTIVES

Recent evidence has linked circadian rhythm dysregulation to an increased risk of metabolic disorders. This study explores a potential association between variation in genes regulating the endogenous circadian timing system (clock genes) and the risk of type 2 diabetes (T2D) in a sample of Greek elderly people.

STUDY DESIGN

Variants within and upstream or downstream of PPARA, PPARD, CLOCK/TMEM165, PER1, PER2 and PER3 genes were genotyped in 716 individuals with T2D (A) and 569 normoglycemic controls (B), and allele frequencies were compared between the groups in a case control study design.

MAIN OUTCOME MEASURES

Samples were genotyped on Illumina Human PsychArray. Permutation test analysis was implemented to determine statistical significance. To avoid the possibility of subjects with prediabetes being included in the control group, people with HbA1c <5.7% and fasting glucose <100 mg/dl comprised group C (n = 393), for whom a separate analysis was performed (secondary analysis).

RESULTS

A protective role against T2D was identified for 14 variants in the PPARA gene. The rs7291444, rs36125344, rs6008384 in PKDREJ, located upstream of PPARA, and rs2859389 in UTS2, located upstream of PER3, demonstrated a protective role against T2D in both analyses. In contrast, rs6744132, located between HES6 and PER2, was positively correlated with T2D risk. Only in the secondary analysis, rs2278637 in VAMP2, located downstream of PER1, and rs11943456 in CLOCK/TMEM165 were found to confer protection against T2D. In a recessive model analysis of all groups, PPARD variants exhibited a protective role against disease.

CONCLUSIONS

Our findings suggest a possible implication of clock genes in T2D susceptibility. Further studies are needed to clarify the mechanisms that connect circadian rhythm dysfunction and T2D pathogenesis.

Physiologic Consequences of Caveolin-1 Ablation in Conventional Outflow Endothelia

Purpose

Polymorphisms at the caveolin-1/2 locus are associated with glaucoma and IOP risk and deletion of caveolin-1 (Cav1) in mice elevates IOP and reduces outflow facility. However, the specific location/cell type responsible for Cav1-dependent regulation of IOP is unclear. We hypothesized that endothelial Cav1 in the conventional outflow (CO) pathway regulate IOP via endothelial nitric oxide synthase (eNOS) signaling.

Methods

We created a mouse with targeted deletion of Cav1 in endothelial cells (Cav1ΔEC) and evaluated IOP, outflow facility, outflow pathway distal vascular morphology, eNOS phosphorylation, and tyrosine nitration of iridocorneal angle tissues by Western blotting.

Results

Endothelial deletion of Cav1 resulted in significantly elevated IOP versus wild-type mice but not a concomitant decrease in outflow facility. Endothelial Cav1 deficiency did not alter the trabecular meshwork or Schlemm's canal morphology, suggesting that the effects observed were not due to developmental deformities. Endothelial Cav1 deletion resulted in eNOS hyperactivity, modestly increased protein nitration, and significant enlargement of the drainage vessels distal to Schlemm's canal. L-Nitro-arginine methyl ester treatment reduced outflow in Cav1ΔEC but not wild-type mice and had no effect on the size of drainage vessels. Endothelin-1 treatment decrease the outflow and drainage vessel size in both wild-type and Cav1ΔEC mice.

Conclusions

Our results suggest that hyperactive eNOS signaling in the CO pathway of both Cav1ΔEC and global Cav1 knockout mice results in chronic dilation of distal CO vessels and protein nitration, but that Cav1 expression in the trabecular meshwork is sufficient to rescue CO defects reported in global Cav1 knockout mice.

Analysis of risk allele frequencies of single nucleotide polymorphisms related to open-angle glaucoma in different ethnic groups

BACKGROUND

The prevalence of open-angle glaucoma (OAG) varies from 0.5% to 7.0% among populations of diverse ancestry, suggesting the existence of genetic differences. The purposes of this study were to provide insights into genetic causes of OAG, which can result in prevalence and phenotype differences among populations of diverse ancestry for OAG, and to compare allele frequencies of intraocular pressure (IOP) elevation-related SNPs in OAG among Koreans and other ethnic groups.

METHODS

We collected the data on a total of 135 OAG-associated single nucleotide polymorphisms (SNPs) from a genome-wide association studies (GWAS) catalog. The population-level allele frequencies of these SNPs were derived based on the 1000 Genomes Project and Korean Reference Genome Database. We used Fisher's exact test to assess whether the effect allele at a given SNP was significantly enriched or depleted.

RESULTS

European, American, and South Asian populations showed similar heatmap patterns, while African, East Asian, and Korean populations had distinct patterns. Korean population presented different profiles compared to other groups; rs1579050 (FMNL2 gene), rs2024211 (CAV2;CAV1), and rs8141433 (GNB1L;TXNRD2 gene), which are known to be associated with IOP variation, were enriched in Americans, Europeans, and Africans, and depleted in Koreans. These can be the candidates for the causative genes of differences in the prevalence of IOP variation in OAG according to ethnic groups.

CONCLUSIONS

Differences in allele frequencies associated with IOP related SNPs between Koreans and other ethnicities were observed, which may explain the high prevalence of OAG with normal IOP predominantly in Koreans and East Asians.

The vital role for nitric oxide in intraocular pressure homeostasis

Catalyzed by endothelial nitric oxide (NO) synthase (eNOS) activity, NO is a gaseous signaling molecule maintaining endothelial and cardiovascular homeostasis. Principally, NO regulates the contractility of vascular smooth muscle cells and permeability of endothelial cells in response to either biochemical or biomechanical cues. In the conventional outflow pathway of the eye, the smooth muscle-like trabecular meshwork (TM) cells and Schlemm's canal (SC) endothelium control aqueous humor outflow resistance, and therefore intraocular pressure (IOP). The mechanisms by which outflow resistance is regulated are complicated, but NO appears to be a key player as enhancement or inhibition of NO signaling dramatically affects outflow function; and polymorphisms in NOS3, the gene that encodes eNOS modifies the relation between various environmental exposures and glaucoma. Based upon a comprehensive review of past foundational studies, we present a model whereby NO controls a feedback signaling loop in the conventional outflow pathway that is sensitive to changes in IOP and its oscillations. Thus, upon IOP elevation, the outflow pathway tissues distend, and the SC lumen narrows resulting in increased SC endothelial shear stress and stretch. In response, SC cells upregulate the production of NO, relaxing neighboring TM cells and increasing permeability of SC's inner wall. These IOP-dependent changes in the outflow pathway tissues reduce the resistance to aqueous humor drainage and lower IOP, which, in turn, diminishes the biomechanical signaling on SC. Similar to cardiovascular pathogenesis, dysregulation of the eNOS/NO system leads to dysfunctional outflow regulation and ocular hypertension, eventually resulting in primary open-angle glaucoma.

Caveolin-1 Ablation Imparts Partial Protection Against Inner Retinal Injury in Experimental Glaucoma and Reduces Apoptotic Activation

Retinal ganglion cell degeneration is a characteristic feature of glaucoma, and accordingly, protection of these cells constitutes a major therapeutic objective in the disease. Here, we demonstrate the key influence of caveolin (Cav) in regulating the inner retinal homeostasis in two models of experimentally elevated intraocular pressure (IOP). Two groups of Cav-1^-/- and wild-type mice were used in the study. Animals were subjected to experimentally induced chronic and acutely elevated IOP and any changes in their retinal function were assessed by positive scotopic threshold response recordings. TUNEL and cleaved caspase-3 assays were performed to evaluate apoptotic changes in the retina while Brn3a immunostaining was used as a marker to assess and quantify ganglion cell layer (GCL) changes. H&E staining was carried out on retinal sections to evaluate histological differences in retinal laminar structure. Cav-1 ablation partially protected the inner retinal function in both chronic and acute models of elevated IOP. The protective effects of Cav-1 loss were also evident histologically by reduced loss of GCL density in both models. The phenotypic protection in Cav-1^-/- glaucoma mice paralleled with increased TrkB phosphorylation and reduced endoplasmic reticulum stress markers and apoptotic activation in the inner retinas. This study corroborated previous findings of enhanced Shp2 phosphorylation in a chronic glaucoma model and established a novel role of Cav-1 in mediating activation of this phosphatase in the inner retina in vivo. Collectively, these findings highlight the critical involvement of Cav-1 regulatory mechanisms in ganglion cells in response to increased IOP, implicating Cav-1 as a potential therapeutic target in glaucoma.

Primary open angle glaucoma genetics: The common variants and their clinical associations (Review)

Glaucoma is a group of progressive optic neuropathies that have in common characteristic optic nerve head changes, loss of retinal ganglion cells and visual field defects. Among the large family of glaucomas, primary open‑angle glaucoma (POAG) is the most common type, a complex and heterogeneous disorder with environmental and genetic factors contributing to its pathogenesis. Approximately 5% of POAG is currently attributed to single‑gene or Mendelian forms of glaucoma. Genetic linkage analysis and genome‑wide association studies have identified various genomic loci, paving the path to understanding the pathogenesis of this enigmatic, blinding disease. In this review we summarize the most common variants reported thus far and their possible clinical correlations.

Association of the CAV1-CAV2 locus with normal-tension glaucoma in Chinese and Japanese

BACKGROUND

The CAV1-CAV2 locus has been associated with primary open-angle glaucoma (POAG) and intraocular pressure. However, its association with normal-tension glaucoma (NTG) was inconclusive. Therefore, we evaluated this association in Chinese and Japanese.

METHODS

Two single-nucleotide polymorphisms (SNPs, rs4236601 and rs1052990) from previous genome-wide association studies of POAG were genotyped in a total of 2220 study subjects: a Hong Kong Chinese cohort of 537 NTG patients and 490 controls, a Shantou Chinese cohort of 102 NTG and 731 controls and an Osaka Japanese cohort of 153 NTG and 207 controls. Subgroup analysis by gender was conducted. Outcomes from different cohorts were combined using meta-analysis.

RESULTS

SNP rs4236601 was significantly associated with NTG in the two Chinese cohorts (P_meta = .0019, OR = 4.55, I² = 0). In contrast, rs4236601 was monomorphic in the Osaka cohort. The association of rs1052990 was insignificant in a meta-analysis combining Chinese and Japanese cohorts (P_meta = .81, OR = 1.05; I² = 64%), and the OR tended towards opposite directions between Chinese (OR = 1.26) and Japanese (OR = 0.69). Gender-specific effects of the SNPs were not statistically significant in the logistic regression or Breslow-day tests of ORs (P > .05), although rs4236601 was significant in males (P = .0068; OR = 10.30) but not in females (P = .14; OR = 2.65) in the meta-analysis of Chinese subjects.

CONCLUSIONS

In this study, we confirmed the association of rs4236601 at the CAV1-CAV2 locus with NTG in Chinese. SNP rs4236601 is monomorphic, and rs1052990 tends towards a different direction in the Japanese cohort. Further studies are warranted to verify the ethnic difference and gender-specific effects of this locus.

Using BXD mouse strains in vision research: A systems genetics approach

We illustrate the growing power of the BXD family of mice (recombinant inbred strains from a cross of C57BL/6J and DBA/2J mice) and companion bioinformatic tools to study complex genome-phenome relations related to glaucoma. Over the past 16 years, our group has integrated powerful murine resources and web-accessible tools to identify networks modulating visual system traits-from photoreceptors to the visual cortex. Recent studies focused on retinal ganglion cells and glaucoma risk factors, including intraocular pressure (IOP), central corneal thickness (CCT), and susceptibility of cellular stress. The BXD family was exploited to define key gene variants and then establish linkage to glaucoma in human cohorts. The power of this experimental approach to precision medicine is highlighted by recent studies that defined cadherin 11 (Cdh11) and a calcium channel (Cacna2d1) as genes modulating IOP, Pou6f2 as a genetic link between CCT and retinal ganglion cell (RGC) death, and Aldh7a1 as a gene that modulates the susceptibility of RGCs to death after elevated IOP. The role of three of these gene variants in glaucoma is discussed, along with the pathways activated in the disease process.

Caveolin-1 regulates human trabecular meshwork cell adhesion, endocytosis, and autophagy

Impaired trabecular meshwork (TM) outflow is implicated in the pathogenesis of primary open-angle glaucoma (POAG). We previously identified the association of a caveolin-1 (CAV1) variant with POAG by genome-wide association study. Here we report a study of CAV1 knockout (KO) effect on human TM cell properties. We generated human CAV1-KO TM cells by CRISPR/Cas9 technology, and we found that the CAV1-KO TM cells less adhered to the surface coating than the wildtype TM cells by 69.34% ( P < 0.05), but showed no difference in apoptosis. Higher endocytosis ability of dextran and transferrin was also observed in the CAV1-KO TM cells (4.37 and 1.89-fold respectively, P < 0.001), compared to the wildtype TM cells. Moreover, the CAV1-KO TM cells had higher expression of extracellular matrix-degrading enzyme genes ( ADMTS13 and MMP14) as well as autophagy-related genes ( ATG7 and BECN1) and protein (LC3B-II) than the wildtype TM cells. In summary, results from this study showed that the CAV1-KO TM cells have reduced adhesion with higher extracellular matrix-degrading enzyme expression, but increased endocytosis and autophagy activities, indicating that CAV1 could be involved in the regulation of adhesion, endocytosis, and autophagy in human TM cells.

Association of polymorphism rs11656696 in GAS7 with primary open-Angle Glaucoma in a Chinese Population

Background: It has been shown that genetic factors play an important role in the pathogenesis of primary open-angle glaucoma (POAG). This study was conducted to investigate the association between the polymorphism rs11656696 located in the growth arrest-specific 7 gene (GAS7) and POAG. Methods: A cohort of 799 unrelated POAG patients and 799 unrelated control subjects was enrolled in this case-control association study. The polymorphism rs11656696 was genotyped using the SNaPshot method. The genotype and allele frequencies were evaluated using the χ² tests. Results: The allele frequency distribution of rs11656696 in the GAS7 gene showed that there was significant difference between POAG cases and controls (P= .006448, OR = 0.82, 95%CI = (0.72-0.95). The minor "A" allele frequency of this polymorphism was 0.477 in the POAG cases, whereas it was 0.526 in controls, suggesting a protective effect for POAG. Significant associations were detected under the homozygous model (p = .006425, OR = 0.68, 95%CI = 0.51-0.90) and recessive model (p = .0003432, OR = 0.66, 95%CI = 0.52-0.84), indicating that subjects carrying rs11656696 AA genotype were less likely to suffer from POAG than those carrying AC/CC genotypes. Conclusion: This case-control association study showed that polymorphism rs11656696 in GAS7 is related to POAG and might be a protective factor against POAG.

The African Descent and Glaucoma Evaluation Study (ADAGES) III: Contribution of Genotype to Glaucoma Phenotype in African Americans: Study Design and Baseline Data

PURPOSE

To describe the study protocol and baseline characteristics of the African Descent and Glaucoma Evaluation Study (ADAGES) III.

DESIGN

Cross-sectional, case-control study.

PARTICIPANTS

Three thousand two hundred sixty-six glaucoma patients and control participants without glaucoma of African or European descent were recruited from 5 study centers in different regions of the United States.

METHODS

Individuals of African descent (AD) and European descent (ED) with primary open-angle glaucoma (POAG) and control participants completed a detailed demographic and medical history interview. Standardized height, weight, and blood pressure measurements were obtained. Saliva and blood samples to provide serum, plasma, DNA, and RNA were collected for standardized processing. Visual fields, stereoscopic disc photographs, and details of the ophthalmic examination were obtained and transferred to the University of California, San Diego, Data Coordinating Center for standardized processing and quality review.

MAIN OUTCOME MEASURES

Participant gender, age, race, body mass index, blood pressure, history of smoking and alcohol use in POAG patients and control participants were described. Ophthalmic measures included intraocular pressure, visual field mean deviation, central corneal thickness, glaucoma medication use, or past glaucoma surgery. Ocular conditions, including diabetic retinopathy, age-related macular degeneration, and past cataract surgery, were recorded.

RESULTS

The 3266 ADAGES III study participants in this report include 2146 AD POAG patients, 695 ED POAG patients, 198 AD control participants, and 227 ED control participants. The AD POAG patients and control participants were significantly younger (both, 67.4 years) than ED POAG patients and control participants (73.4 and 70.2 years, respectively). After adjusting for age, AD POAG patients had different phenotypic characteristics compared with ED POAG patients, including higher intraocular pressure, worse visual acuity and visual field mean deviation, and thinner corneas (all P < 0.001). Family history of glaucoma did not differ between AD and ED POAG patients.

CONCLUSIONS

With its large sample size, extensive specimen collection, and deep phenotyping of AD and ED glaucoma patients and control participants from different regions in the United States, the ADAGES III genomics study will address gaps in our knowledge of the genetics of POAG in this high-risk population.

Association analysis of exome variants and refraction, axial length, and corneal curvature in a European-American population

Refractive errors, myopia, and hyperopia are common visual disorders greatly affecting older individuals. Refraction is determined by genetic factors but only a small percentage of its variation has been explained. We performed a genetic association analysis with three ocular phenotypes: spherical equivalent (a continous measure of refraction), axial length, and corneal curvature in 1,871 European-Americans from the Beaver Dam Eye Study. Individuals were genotyped on the Illumina exome array and imputed to the Haplotype Reference Consortium reference panel. After increasing the number of analyzed variants in targeted protein-coding regions 10-fold via imputation, we confirmed associations for two previously known loci with corneal curvature (chr4q12, rs2114039; g.55092626T > C, β = -0.03 (95% confidence interval [CI]): -0.06, -0.01, P value = 0.01) and spherical equivalent (chr15q14, rs634990; g.35006073T > C, β = -0.27, 95% CI: -0.45, -0.09, P value = 3.79 × 10-3 ). Despite increased single nucleotide polymorphism (SNP) density, we did not detect any novel significant variants after correction for multiple comparisons. In summary, we confirmed two previous loci associated with corneal curvature and spherical equivalent in a European-American population highlighting the potential biological role of those regions in these traits.

Genomic Locus Modulating IOP in the BXD RI Mouse Strains

Intraocular pressure (IOP) is the primary risk factor for developing glaucoma, yet little is known about the contribution of genomic background to IOP regulation. The present study leverages an array of systems genetics tools to study genomic factors modulating normal IOP in the mouse. The BXD recombinant inbred (RI) strain set was used to identify genomic loci modulating IOP. We measured the IOP in a total of 506 eyes from 38 different strains. Strain averages were subjected to conventional quantitative trait analysis by means of composite interval mapping. Candidate genes were defined, and immunohistochemistry and quantitative PCR (qPCR) were used for validation. Of the 38 BXD strains examined the mean IOP ranged from a low of 13.2mmHg to a high of 17.1mmHg. The means for each strain were used to calculate a genome wide interval map. One significant quantitative trait locus (QTL) was found on Chr.8 (96 to 103 Mb). Within this 7 Mb region only 4 annotated genes were found: Gm15679, Cdh8, Cdh11 and Gm8730. Only two genes (Cdh8 and Cdh11) were candidates for modulating IOP based on the presence of non-synonymous SNPs. Further examination using SIFT (Sorting Intolerant From Tolerant) analysis revealed that the SNPs in Cdh8 (Cadherin 8) were predicted to not change protein function; while the SNPs in Cdh11 (Cadherin 11) would not be tolerated, affecting protein function. Furthermore, immunohistochemistry demonstrated that CDH11 is expressed in the trabecular meshwork of the mouse. We have examined the genomic regulation of IOP in the BXD RI strain set and found one significant QTL on Chr. 8. Within this QTL, there is one good candidate gene, Cdh11.

Caveolins as Regulators of Stress Adaptation

Caveolins have been recognized over the past few decades as key regulators of cell physiology. They are ubiquitously expressed and regulate a number of processes that ultimately impact efficiency of cellular processes. Though not critical to life, they are central to stress adaptation in a number of organs. The following review will focus specifically on the role of caveolin in stress adaptation in the heart, brain, and eye, three organs that are susceptible to acute and chronic stress and that show as well declining function with age. In addition, we consider some novel molecular mechanisms that may account for this stress adaptation and also offer potential to drive the future of caveolin research.

Genetic Risk Factors for Glaucoma and Exfoliation Syndrome Identified by Genome-wide Association Studies

BACKGROUND

Glaucoma is a neurodegenerative disease characterized by the progressive loss of retinal ganglion cells and optic nerve axons. According to its anatomical features, glaucoma is mainly subdivided into primary open-angle glaucoma (POAG) and primary angle-closure glaucoma (PACG). Exfoliation syndrome (XFS) and glaucoma (XFG) are characterized by the accumulation of extracellular materials in ocular tissues, particularly the lens surface and pupillary border. In addition to the two major forms of glaucoma, XFG is the most common cause of secondary open-angle glaucoma. Recent genome-wide association studies(GWASs) revealed genetic loci associated with each glaucoma subtype.

METHODS

Review of literatures regarding GWASs for POAG, PACG and XFS.

RESULTS

Several genetic loci were found to be independently associated with POAG, PACG, and XFS by large-scale GWASs.

CONCLUSIONS

Genetic studies may not only provide a better understanding of the pathophysiological mechanisms underlying the diseases, but also facilitate the development of new drugs or treatments.

A large multi-ethnic genome-wide association study identifies novel genetic loci for intraocular pressure

Elevated intraocular pressure (IOP) is a major risk factor for glaucoma, a leading cause of blindness. IOP heritability has been estimated to up to 67%, and to date only 11 IOP loci have been reported, accounting for 1.5% of IOP variability. Here, we conduct a genome-wide association study of IOP in 69,756 untreated individuals of European, Latino, Asian, and African ancestry. Multiple longitudinal IOP measurements were collected through electronic health records and, in total, 356,987 measurements were included. We identify 47 genome-wide significant IOP-associated loci (P < 5 × 10⁻⁸); of the 40 novel loci, 14 replicate at Bonferroni significance in an external genome-wide association study analysis of 37,930 individuals of European and Asian descent. We further examine their effect on the risk of glaucoma within our discovery sample. Using longitudinal IOP measurements from electronic health records improves our power to identify new variants, which together explain 3.7% of IOP variation.

Intraocular pressure (IOP) is a major risk factor for glaucoma. Here, Choquet and co-authors perform a multi-ethnic genome-wide association study of repeat IOP measurements in 69,756 individuals and identify 40 novel loci, 36 of which show directionally consistent effects in glaucoma.

Investigation of CAV1/CAV2 rs4236601 and CDKN2B-AS1 rs2157719 in primary open-angle glaucoma patients from Brazil

Large-scale genome-wide association studies have identified several susceptibility variants associated with the risk of primary open-angle glaucoma (POAG), among which rs4236601 (CAV1/CAV2) at chromosome 7q31 and rs2157719 at chromosome 9p21 (CDKN2B-AS1). The purpose of this study was to investigate whether these variants contribute to the incidence of POAG in a sample of the Brazilian Southeastern population and to determine the best-fitted genetic model for these single nucleotide polymorphisms (SNPs). A case-control study with 557 individuals, 310 with POAG, and 247 controls was conducted through PCR and direct sequencing. We observed a significant effect of the heterozygous genotype (G/A) of rs2157719 that occurred more frequently in the control group (p = 0.0004; OR: 0.517, CI 95%: 0.357-0.745). Allele frequencies also differed between cases and controls (p = 0.006; OR: 0.694, CI 95%: 0.522-0.922) with the best-fitted genetic model for rs2157719 being the codominant model. No differences were observed for genotype and allele distributions in relation to rs4236601 in the CAV1/CAV2 region. The association of rs2157719 (CDKN2B-AS1) with the POAG phenotype corroborates previously published results, reinforcing the importance of this variant in POAG etiology.

Polymorphism rs11656696 in GAS7 Is Not Associated with Primary Open Angle Glaucoma in a Saudi Cohort

AIMS

To conduct a case-control study to investigate the association between the polymorphism rs11656696 located in the growth arrest-specific 7 gene (GAS7)on human chromosome 17p13.1 and primary open angle glaucoma (POAG).

METHODS

The polymorphism rs11656696 was genotyped using the TaqMan^® assay in 187 subjects comprising 92 unrelated POAG cases and 95 controls of Saudi Arabian origin.

RESULTS

Association analysis between cases and controls revealed no significant genotype distribution under additive (p = 0.225), dominant (p = 0.635), or recessive (p = 0.085) models. Moreover, the allele frequency distribution was also nonsignificant (p = 0.70). The minor "A" allele frequency was 0.35 and 0.41 among POAG cases and controls, respectively. In addition, specific clinical indices used to assess severity of glaucoma such as intraocular pressure (IOP), cup/disk ratio, and number of antiglaucoma medications also did not show any significant genotype distribution in POAG cases. Moreover, a binary logistic regression analysis did not show any significant effect of age, sex, or genotype on disease outcome.

CONCLUSION

Polymorphism rs11656696 is not associated with POAG nor any of its endophenotypic traits such as IOP and cup/disk ratio and is thus not a risk factor for POAG in this Saudi cohort.

Major review: Molecular genetics of primary open-angle glaucoma

Glaucoma is a leading cause of irreversible blindness worldwide. Primary open-angle glaucoma (POAG), the most common type, is a complex inherited disorder that is characterized by progressive retinal ganglion cell death, optic nerve head excavation, and visual field loss. The discovery of a large, and growing, number of genetic and chromosomal loci has been shown to contribute to POAG risk, which carry implications for disease pathogenesis. Differential gene expression analyses in glaucoma-affected tissues as well as animal models of POAG are enhancing our mechanistic understanding in this common, blinding disorder. In this review we summarize recent developments in POAG genetics and molecular genetics research.

Caveolin-1 modulates intraocular pressure: implications for caveolae mechanoprotection in glaucoma

Polymorphisms in the CAV1/2 genes that encode signature proteins of caveolae are associated with glaucoma, the second leading cause of blindness worldwide, and with its major risk factor, intraocular pressure (IOP). We hypothesized that caveolin-1 (Cav-1) participates in IOP maintenance via modulation of aqueous humor drainage from the eye. We localize caveolae proteins to human and murine conventional drainage tissues and show that caveolae respond to mechanical stimulation. We show that Cav-1-deficient (Cav-1^−/−) mice display ocular hypertension explained by reduced pressure-dependent drainage of aqueous humor. Cav-1 deficiency results in loss of caveolae in the Schlemm’s canal (SC) and trabecular meshwork. However, their absence did not appear to impact development nor adult form of the conventional outflow tissues according to rigorous quantitative ultrastructural analyses, but did affect cell and tissue behavior. Thus, when IOP is experimentally elevated, cells of the Cav-1^−/− outflow tissues are more susceptible to plasma membrane rupture indicating that caveolae play a role in mechanoprotection. Additionally, aqueous drainage from Cav-1^−/− eyes was more sensitive to nitric oxide (NO) synthase inhibition than controls, suggesting that excess NO partially compensates for outflow pathway dysfunction. These results provide a functional link between a glaucoma risk gene and glaucoma-relevant pathophysiology.

Primary open-angle glaucoma

Glaucoma is an optic neuropathy that is characterized by the progressive degeneration of the optic nerve, leading to visual impairment. Glaucoma is the main cause of irreversible blindness worldwide, but typically remains asymptomatic until very severe. Open-angle glaucoma comprises the majority of cases in the United States and western Europe, of which, primary open-angle glaucoma (POAG) is the most common type. By contrast, in China and other Asian countries, angle-closure glaucoma is highly prevalent. These two types of glaucoma are characterized based on the anatomic configuration of the aqueous humour outflow pathway. The pathophysiology of POAG is not well understood, but it is an optic neuropathy that is thought to be associated with intraocular pressure (IOP)-related damage to the optic nerve head and resultant loss of retinal ganglion cells (RGCs). POAG is generally diagnosed during routine eye examination, which includes fundoscopic evaluation and visual field assessment (using perimetry). An increase in IOP, measured by tonometry, is not essential for diagnosis. Management of POAG includes topical drug therapies and surgery to reduce IOP, although new therapies targeting neuroprotection of RGCs and axonal regeneration are under development.

Caveolins and caveolae in ocular physiology and pathophysiology

Caveolae are specialized, invaginated plasma membrane domains that are defined morphologically and by the expression of signature proteins called, caveolins. Caveolae and caveolins are abundant in a variety of cell types including vascular endothelium, glia, and fibroblasts where they play critical roles in transcellular transport, endocytosis, mechanotransduction, cell proliferation, membrane lipid homeostasis, and signal transduction. Given these critical cellular functions, it is surprising that ablation of the caveolae organelle does not result in lethality suggesting instead that caveolae and caveolins play modulatory roles in cellular homeostasis. Caveolar components are also expressed in ocular cell types including retinal vascular cells, Müller glia, retinal pigment epithelium (RPE), conventional aqueous humor outflow cells, the corneal epithelium and endothelium, and the lens epithelium. In the eye, studies of caveolae and other membrane microdomains (i.e., "lipid rafts") have lagged behind what is a substantial body of literature outside vision science. However, interest in caveolae and their molecular components has increased with accumulating evidence of important roles in vision-related functions such as blood-retinal barrier homeostasis, ocular inflammatory signaling, pathogen entry at the ocular surface, and aqueous humor drainage. The recent association of CAV1/2 gene loci with primary open angle glaucoma and intraocular pressure has further enhanced the need to better understand caveolar functions in the context of ocular physiology and disease. Herein, we provide the first comprehensive review of literature on caveolae, caveolins, and other membrane domains in the context of visual system function. This review highlights the importance of caveolae domains and their components in ocular physiology and pathophysiology and emphasizes the need to better understand these important modulators of cellular function.

A Germline Variant in the PANX1 Gene Has Reduced Channel Function and Is Associated with Multisystem Dysfunction

Pannexin1 (PANX1) is probably best understood as an ATP release channel involved in paracrine signaling. Given its ubiquitous expression, PANX1 pathogenic variants would be expected to lead to disorders involving multiple organ systems. Using whole exome sequencing, we discovered the first patient with a homozygous PANX1 variant (c.650G→A) resulting in an arginine to histidine substitution at position 217 (p.Arg217His). The 17-year-old female has intellectual disability, sensorineural hearing loss requiring bilateral cochlear implants, skeletal defects, including kyphoscoliosis, and primary ovarian failure. Her consanguineous parents are each heterozygous for this variant but are not affected by the multiorgan syndromes noted in the proband. Expression of the p.Arg217His mutant in HeLa, N2A, HEK293T, and Ad293 cells revealed normal PANX1 glycosylation and cell surface trafficking. Dye uptake, ATP release, and electrophysiological measurements revealed p.Arg217His to be a loss-of-function variant. Co-expression of the mutant with wild-type PANX1 suggested the mutant was not dominant-negative to PANX1 channel function. Collectively, we demonstrate a PANX1 missense change associated with human disease in the first report of a "PANX1-related disorder."

Nailfold Capillary Abnormalities in Primary Open-Angle Glaucoma: A Multisite Study

PURPOSE

There is considerable evidence for systemic vascular dysfunction in primary open-angle glaucoma (POAG). We performed nailfold capillary video microscopy to observe directly the nature of nonocular microvasculature abnormalities in POAG.

METHODS

We enrolled 199 POAG patients and 124 control subjects from four sites. We used JH-1004 capillaroscopes to perform nailfold capillary video microscopy on the fourth and fifth digits of each subject's nondominant hand. Videos were evaluated for hemorrhages, dilated capillary loops > 50 μm, and avascular zones > 100 μm by graders masked to case status. Multivariable odds ratios (ORs) and 95% confidence intervals (CIs) for POAG were obtained by means of logistic regression analyses that were applied to data from all cases and controls. Corresponding estimates of moderate or severe POAG versus mild POAG (based on the Hodapp-Anderson-Parrish scale) were obtained among cases only.

RESULTS

After controlling for demographic factors, family history of glaucoma, systemic diseases, and use of anticoagulation and antiplatelet therapy, for each 100 nailfold capillaries assessed, all types of microvascular abnormalities were significantly associated with POAG. Specifically, the presence of any dilated capillaries (OR = 2.9; 95% CI, 1.6-5.6), avascular zones (OR = 4.4; 95% CI, 1.7-11.3) and hemorrhages (OR = 12.2; 95% CI, 5.9-25.1) were associated with POAG. Among cases, the frequency of microvascular abnormalities was not associated with glaucoma severity (P ≥ 0.43).

CONCLUSIONS

These data provided support for nonocular capillary bed abnormalities in POAG. Comparable vascular abnormalities in the optic nerve may render it susceptible to glaucomatous damage.

An Updated Review on the Genetics of Primary Open Angle Glaucoma

Epidemiological studies suggest that by 2020 the prevalence of primary open angle glaucoma (POAG) is estimated to increase to 76.0 million, and to 111.8 million by 2040 globally due to the population aging. The prevalence of POAG is the highest among those of African descent, followed by Asians, and the lowest in Europeans. POAG is a genetically complex trait with a substantial fraction exhibiting a significant heritability. Less than 10% of POAG cases in the general population are caused by specific gene mutations and the remaining cases are polygenic. Quantitative traits related to POAG pathogenesis such as intra-ocular pressure (IOP), vertical cup/disc ratio (VCDR), optic disc area, and central corneal thickness (CCT) are highly heritable, and likely to be influenced at least in part by genes and show substantial variation in human populations. Recent genome-wide association studies (GWAS) have identified several single nucleotide polymorphisms (SNPs) at different loci including CAV1/CAV2, TMCO1, CDKN2B-AS1, CDC7-TGFBR3, SIX1/SIX6, GAS7 and ATOH7 to be associated with POAG and its related quantitative traits (endophenotypes). The chapter provides a brief overview on the different GWAS and SNP association studies and their correlation with various clinical parameters important for POAG in the population worldwide, including the Middle East.

Expression Profiling of Human Schlemm's Canal Endothelial Cells From Eyes With and Without Glaucoma

PURPOSE

Ocular hypertension is a major risk factor for glaucoma and the inner wall of Schlemm's canal (SC) endothelia participates in the regulation of aqueous humor outflow resistance. This study aimed to identify differentially expressed genes in primary cultures of SC cells from glaucoma patients.

METHODS

This study examined SC samples from three glaucoma cases and four controls. Schlemm's canal cells were isolated from eight different postmortem human eyes. Total RNA was extracted, labeled, and hybridized to Illumina HumanWG-6 BeadChips containing probes for approximately 47,000 human transcripts. After extracting the data using Illumina GenomeStudio software, the data were normalized and analyzed using the R package limma in Bioconductor. Using Protein ANalysis THrough Evolutionary Relationships (PANTHER) software, gene ontology analysis of highly expressed genes was executed in controls and glaucoma groups separately. Pathway analysis was performed with differentially expressed genes using WebGestalt (WEB-based GEne SeT AnaLysis Toolkit). Selected genes were validated using droplet digital PCR (ddPCR).

RESULTS

Gene ontology analysis indicated similar functional categories in cases and controls. Differential analysis identified a total of 113 genes with at least 2-fold expression changes in cases. Pathway analysis indicated significant enrichment of genes in cell adhesion, heparin binding, glycosaminoglycan binding, filopodium, and extracellular matrix remodeling. Eighteen selected genes with differential expression were successfully validated using ddPCR.

CONCLUSIONS

This study represents the first genome-wide expression study of human primary SC cells from glaucoma patients and provides a potential list of targets regulating SC cell stiffness and pore formation, eventually the outflow resistance in glaucoma individuals.

Glaucoma Genes and Mechanisms

Genetic studies have yielded important genes contributing to both early-onset and adult-onset forms of glaucoma. The proteins encoded by the current collection of glaucoma genes participate in a broad range of cellular processes and biological systems. Approximately half the glaucoma-related genes function in the extracellular matrix, however proteins involved in cytokine signaling, lipid metabolism, membrane biology, regulation of cell division, autophagy, and ocular development also contribute to the disease pathogenesis. While the function of these proteins in health and disease are not completely understood, recent studies are providing insight into underlying disease mechanisms, a critical step toward the development of gene-based therapies. In this review, genes known to cause early-onset glaucoma or contribute to adult-onset glaucoma are organized according to the cell processes or biological systems that are impacted by the function of the disease-related protein product.

Discovery and Fine-Mapping of Glycaemic and Obesity-Related Trait Loci Using High-Density Imputation

Reference panels from the 1000 Genomes (1000G) Project Consortium provide near complete coverage of common and low-frequency genetic variation with minor allele frequency ≥0.5% across European ancestry populations. Within the European Network for Genetic and Genomic Epidemiology (ENGAGE) Consortium, we have undertaken the first large-scale meta-analysis of genome-wide association studies (GWAS), supplemented by 1000G imputation, for four quantitative glycaemic and obesity-related traits, in up to 87,048 individuals of European ancestry. We identified two loci for body mass index (BMI) at genome-wide significance, and two for fasting glucose (FG), none of which has been previously reported in larger meta-analysis efforts to combine GWAS of European ancestry. Through conditional analysis, we also detected multiple distinct signals of association mapping to established loci for waist-hip ratio adjusted for BMI (RSPO3) and FG (GCK and G6PC2). The index variant for one association signal at the G6PC2 locus is a low-frequency coding allele, H177Y, which has recently been demonstrated to have a functional role in glucose regulation. Fine-mapping analyses revealed that the non-coding variants most likely to drive association signals at established and novel loci were enriched for overlap with enhancer elements, which for FG mapped to promoter and transcription factor binding sites in pancreatic islets, in particular. Our study demonstrates that 1000G imputation and genetic fine-mapping of common and low-frequency variant association signals at GWAS loci, integrated with genomic annotation in relevant tissues, can provide insight into the functional and regulatory mechanisms through which their effects on glycaemic and obesity-related traits are mediated.

Human genetic studies have demonstrated that quantitative human anthropometric and metabolic traits, including body mass index, waist-hip ratio, and plasma concentrations of glucose and insulin, are highly heritable, and are established risk factors for type 2 diabetes and cardiovascular diseases. Although many regions of the genome have been associated with these traits, the specific genes responsible have not yet been identified. By making use of advanced statistical “imputation” techniques applied to more than 87,000 individuals of European ancestry, and publicly available “reference panels” of more than 37 million genetic variants, we have been able to identify novel regions of the genome associated with these glycaemic and obesity-related traits and localise genes within these regions that are most likely to be causal. This improved understanding of the biological mechanisms underlying glycaemic and obesity-related traits is extremely important because it may advance drug development for downstream disease endpoints, ultimately leading to public health benefits.

Expression-associated polymorphisms of CAV1-CAV2 affect intraocular pressure and high-tension glaucoma risk

PURPOSE

The human CAV1-CAV2 locus has been associated with susceptibility to primary open-angle glaucoma in four studies of Caucasian, Chinese, and Pakistani populations, although not in several other studies of non-Korean populations. In this study with Korean participants, the CAV1-CAV2 locus was investigated for associations with susceptibility to primary open-angle glaucoma accompanied by elevated intraocular pressure (IOP), namely, high-tension glaucoma (HTG), as well as with IOP elevation, which is a strong risk factor for glaucoma.

METHODS

Two single nucleotide polymorphisms (SNPs) were genotyped in 1,161 Korean participants including 229 patients with HTG and 932 healthy controls and statistically examined for association with HTG susceptibility and IOP. One SNP was rs4236601 G>A, which had been reported in the original study, and the other SNP was rs17588172 T>G, which was perfectly correlated (r2=1) with another reported SNP rs1052990. Expression quantitative trait loci (eQTL) analysis was performed using GENe Expression VARiation (Genevar) data.

RESULTS

Both SNPs were associated with HTG susceptibility, but the rs4236601 association disappeared when adjusted for the rs17588172 genotype and not vice versa. The minor allele G of rs17588172 was associated significantly with 1.5-fold increased susceptibility to HTG (p=0.0069) and marginally with IOP elevation (p=0.043) versus the major allele T. This minor allele was also associated with decreased CAV1 and CAV2 mRNA in skin and adipose according to the Genevar eQTL analysis.

CONCLUSIONS

The minor allele G of rs17588172 in the CAV1-CAV2 locus is associated with decreased expression of CAV1 and CAV2 in some tissues, marginally with IOP elevation, and consequently with increased susceptibility to HTG.