A Genome-Wide Association Study of Vertical Cup-Disc Ratio in a Latino Population

An Overview of Strategies for Detecting Genotype-Phenotype Associations Across Ancestrally Diverse Populations

Genome-wide association studies (GWAS) have been very successful at identifying genetic variants influencing a large number of traits. Although the great majority of these studies have been performed in European-descent individuals, it has been recognised that including populations with differing ancestries enhances the potential for identifying causal SNPs due to their differing patterns of linkage disequilibrium. However, when individuals from distinct ethnicities are included in a GWAS, it is necessary to implement a number of control steps to ensure that the identified associations are real genotype-phenotype relationships. In this Review, we discuss the analyses that are required when performing multi-ethnic studies, including methods for determining ancestry at the global and local level for sample exclusion, controlling for ancestry in association testing, and post-GWAS interrogation methods such as genomic control and meta-analysis. We hope that this overview provides a primer for those researchers interested in including distinct populations in their studies.

Glaucoma Heritability: Molecular Mechanisms of Disease

Glaucoma is one of the world’s leading causes of irreversible blindness. A complex, multifactorial disease, the underlying pathogenesis and reasons for disease progression are not fully understood. The most common form of glaucoma, primary open-angle glaucoma (POAG), was traditionally understood to be the result of elevated intraocular pressure (IOP), leading to optic nerve damage and functional vision loss. Recently, researchers have suggested that POAG may have an underlying genetic component. In fact, studies of genetic association and heritability have yielded encouraging results showing that glaucoma may be influenced by genetic factors, and estimates for the heritability of POAG and disease-related endophenotypes show encouraging results. However, the vast majority of the underlying genetic variants and their molecular mechanisms have not been elucidated. Several genes have been suggested to have molecular mechanisms contributing to alterations in key endophenotypes such as IOP (LMX1B, MADD, NR1H3, and SEPT9), and VCDR (ABCA1, ELN, ASAP1, and ATOH7). Still, genetic studies about glaucoma and its molecular mechanisms are limited by the multifactorial nature of the disease and the large number of genes that have been identified to have an association with glaucoma. Therefore, further study into the molecular mechanisms of the disease itself are required for the future development of therapies targeted at genes leading to POAG endophenotypes and, therefore, increased risk of disease.

Genetic association between CDKN2B/CDKN2B-AS1 gene polymorphisms with primary glaucoma in a North Indian cohort: an original study and an updated meta-analysis

BACKGROUND

Variants in CDKN2B/CDKN2B-AS1 have been reported to modulate glaucoma risk in several GWAS across different populations. CDKN2B/CDKN2A encodes tumor suppressor proteins p16INK4A/p15INK4B which influences cell proliferation/senescence in RGCs, the degeneration of which is a risk factor for glaucoma. CDKN2B-AS1 codes a long non-coding RNA in antisense direction and is involved in influencing nearby CDKN2A/CDKN2B via regulatory mechanisms.

METHODS

Current study investigated four SNPs (rs2157719, rs3217992, rs4977756, rs1063192) of aforementioned genes in a case-control study in a North Indian cohort. Genotyping was done with Taqman chemistry. In addition, an updated meta-analysis was performed.

RESULTS

Two SNPs, rs3217992 and rs2157719 were found to be significantly associated with the disease. The frequency of 'T' allele of rs3217992 was significantly lower in cases (POAG/PACG) [p = 0.045; OR = 0.80(CI = 0.65-0.99) and p = 0.024; OR = 0.73(CI = 0.55-0.96)], respectively than in controls. Genetic model analysis revealed that TT + CT genotype confers 0.73-fold protection against POAG [p = 0.047; OR = 0.73(CI = 0.54-0.99)] and trend assumed additive model gives 0.53 times higher protection against PACG progression. However the association of rs3217992 with POAG and PACG did not remain significant after Bonferroni correction. For rs2157719, the 'C' allele was found to be less prevalent among cases (POAG/PACG) with respect to controls. Cochran Armitage trend test assuming additive model revealed 0.77 and 0.64-fold protection against POAG and PACG respectively. Bonferroni correction (pcorr = 0.003) was applied and the association of rs2157719 remained significant in PACG cases but not among POAG cases (p = 0.024). The 'CC' genotype also confers protection against primary glaucoma (POAG/PACG) among males and female subjects. The frequency rs1063192 and rs4977756 did not vary significantly among subjects, however the haplotype 'CATA' was found to be associated with increased glaucoma risk. An updated meta-analysis conducted on pooled studies on POAG cases and controls revealed significant association between rs1063192, rs2157719, rs4977756 and POAG except rs3217992.

CONCLUSION

The study concludes significant association between INK4 variants and primary glaucoma in the targeted North Indian Punjabi cohort. We believe that deep-sequencing of INK4 locus may help in identifying novel variants modifying susceptibility to glaucoma. Functional studies can further delineate the role of CDKN2B and CDKN2B-AS1 in primary glaucoma for therapeutic intervention.

Molecular Genetics of Glaucoma: Subtype and Ethnicity Considerations

Glaucoma, the world’s leading cause of irreversible blindness, is a complex disease, with differential presentation as well as ethnic and geographic disparities. The multifactorial nature of glaucoma complicates the study of genetics and genetic involvement in the disease process. This review synthesizes the current literature on glaucoma and genetics, as stratified by glaucoma subtype and ethnicity. Primary open-angle glaucoma (POAG) is the most common cause of glaucoma worldwide, with the only treatable risk factor (RF) being the reduction of intraocular pressure (IOP). Genes associated with elevated IOP or POAG risk include: ABCA1, AFAP1, ARHGEF12, ATXN2, CAV1, CDKN2B-AS1, FOXC1, GAS7, GMDS, SIX1/SIX6, TMCO1, and TXNRD2. However, there are variations in RF and genetic factors based on ethnic and geographic differences; it is clear that unified molecular pathways accounting for POAG pathogenesis remain uncertain, although inflammation and senescence likely play an important role. There are similar ethnic and geographic complexities in primary angle closure glaucoma (PACG), but several genes have been associated with this disorder, including MMP9, HGF, HSP70, MFRP, and eNOS. In exfoliation glaucoma (XFG), genes implicated include LOXL1, CACNA1A, POMP, TMEM136, AGPAT1, RBMS3, and SEMA6A. Despite tremendous progress, major gaps remain in resolving the genetic architecture for the various glaucoma subtypes across ancestries. Large scale carefully designed studies are required to advance understanding of genetic loci as RF in glaucoma pathophysiology and to improve diagnosis and treatment options.

Insights into the regulatory molecules involved in glaucoma pathogenesis

Glaucoma is an important cause of irreversible blindness, characterized by optic nerve anomalies. Increased intraocular pressure (IOP) and aging are major risk factors. Retinal ganglion cells and trabecular meshwork cells are certainly involved in the etiology of glaucoma. Glaucoma is usually a complex disease, and various genes and functions may contribute to its etiology. Among these may be genes that encode regulatory molecules. In this review, regulatory molecules including 18 transcription factors (TFs), 195 microRNAs (miRNAs), 106 long noncoding RNAs (lncRNAs), and two circular RNAs (circRNAs) that are reasonable candidates for having roles in glaucoma pathogenesis are described. The targets of the regulators are reported. Glaucoma-related features including apoptosis, stress responses, immune functions, ECM properties, IOP, and eye development are affected by the targeted genes. The targeted genes that are frequently targeted by multiple regulators most often affect apoptosis and the related features of cell death and cell survival. BCL2, CDKN1A, and TP53 are among the frequent targets of three types of glaucoma-relevant regulators, TFs, miRNAs, and lncRNAs. TP53 was itself identified as a glaucoma-relevant TF. Several of the glaucoma-relevant TFs are themselves among frequent targets of regulatory molecules, which is consistent with existence of a complex network involved in glaucoma pathogenesis.

The Atoh7 remote enhancer provides transcriptional robustness during retinal ganglion cell development

The retinal ganglion cell (RGC) competence factor ATOH7 is dynamically expressed during retinal histogenesis. ATOH7 transcription is controlled by a promoter-adjacent primary enhancer and a remote shadow enhancer (SE). Deletion of the ATOH7 human SE causes nonsyndromic congenital retinal nonattachment (NCRNA) disease, characterized by optic nerve aplasia and total blindness. We used genome editing to model NCRNA in mice. Deletion of the murine SE reduces Atoh7 messenger RNA (mRNA) fivefold but does not recapitulate optic nerve loss; however, SEdel/knockout (KO) trans heterozygotes have thin optic nerves. By analyzing Atoh7 mRNA and protein levels, RGC development and survival, and chromatin landscape effects, we show that the SE ensures robust Atoh7 transcriptional output. Combining SE deletion and KO and wild-type alleles in a genotypic series, we determined the amount of Atoh7 needed to produce a normal complement of adult RGCs, and the secondary consequences of graded reductions in Atoh7 dosage. Together, these data reveal the workings of an evolutionary fail-safe, a duplicate enhancer mechanism that is hard-wired in the machinery of vertebrate retinal ganglion cell genesis.

The Genetic and Endoplasmic Reticulum-Mediated Molecular Mechanisms of Primary Open-Angle Glaucoma

Glaucoma is a heterogenous, chronic, progressive group of eye diseases, which results in irreversible loss of vision. There are several types of glaucoma, whereas the primary open-angle glaucoma (POAG) constitutes the most common type of glaucoma, accounting for three-quarters of all glaucoma cases. The pathological mechanisms leading to POAG pathogenesis are multifactorial and still poorly understood, but it is commonly known that significantly elevated intraocular pressure (IOP) plays a crucial role in POAG pathogenesis. Besides, genetic predisposition and aggregation of abrogated proteins within the endoplasmic reticulum (ER) lumen and subsequent activation of the protein kinase RNA-like endoplasmic reticulum kinase (PERK)-dependent unfolded protein response (UPR) signaling pathway may also constitute important factors for POAG pathogenesis at the molecular level. Glaucoma is commonly known as a ‘silent thief of sight’, as it remains asymptomatic until later stages, and thus its diagnosis is frequently delayed. Thereby, detailed knowledge about the glaucoma pathophysiology is necessary to develop both biochemical and genetic tests to improve its early diagnosis as well as develop a novel, ground-breaking treatment strategy, as currently used medical therapies against glaucoma are limited and may evoke numerous adverse side-effects in patients.

Loss of foxc1 in zebrafish reduces optic nerve size and cell number in the retinal ganglion cell layer

Mutation of FOXC1 causes Axenfeld-Rieger Syndrome (ARS) with early onset or congenital glaucoma. We assessed retinal ganglion cell (RGC) number in zebrafish due to CRISPR-mediated mutation and antisense inhibition of two-forkhead box transcription factors, foxc1a and foxc1b. These genes represent duplicated homologues of human FOXC1. Using a CRISPR induced null mutation in foxc1b, in combination with antisense inhibition of foxc1a, we demonstrate reduced cell number in the retinal ganglion cell layer of developing zebrafish eyes. As early as 5 days post fertilization (dpf), fewer RGCs are found in foxc1b homozygous mutants injected with foxc1a morpholinos, and a thinner optic nerve results. Our data illustrates that foxc1 is required for the expression of atonal homolog 7 (atoh7), a gene that is necessary for RGC differentiation. As markers of differentiated RGCs (pou4f2) are downregulated in foxc1b-/- mutants injected with foxc1a morpholinos and no cell death is observed, our results are consistent with defects in the differentiation of RGCs leading to reduced cell number, as opposed to increased cell death of RGCs or off targets effects of morpholino injection. Our zebrafish model demonstrates that aberrant regulation of RGC number could act in concert with other known glaucoma risk factors to influence the development of congenital and early onset glaucoma due to FOXC1 mutation.

Lack of Association between Variant rs7916697 in ATOH7 and Primary Open Angle Glaucoma in a Saudi Cohort

A case-control genetic association study was performed to investigate whether variant rs7916697 in atonal bHLH transcription factor 7 (ATOH7), which has been previously reported to be associated with optic disc parameters and primary open angle glaucoma (POAG) in different ethnic groups, is a risk factor for POAG or any of its clinical phenotypes in a Saudi cohort. Genotyping of rs7916697 (G>A) variant was performed in 186 unrelated POAG cases and 171 unrelated nonglaucomatous controls of Saudi origin using real-time Taq-Man® assay. Genotypic and allelic association with POAG and its related clinical indices were evaluated. Demographic and systemic disease status did not differ significantly between POAG cases and controls. Association analysis between POAG cases and controls showed no significant genotype effect under additive (p=0.707), dominant (p=0.458), and recessive (p=0.554) models. Besides, the minor ‘A' allele frequency was 0.39 in POAG cases and 0.36 in controls with no significant distribution (p=0.406). In addition, there was no significant difference between genotypes and clinical phenotypes such as intraocular pressure and cup/disc ratio within the POAG group, or any age and sex adjusted genotype effect on the disease outcome in regression analysis. Variant rs7916697 in ATOH7 is not associated with POAG or its clinical indices such as IOP and cup/disc ratio in a Saudi cohort.

Genetic Risk Score Is Associated with Vertical Cup-to-Disc Ratio and Improves Prediction of Primary Open-Angle Glaucoma in Latinos

PURPOSE

Genome-wide association studies have identified multiple genetic variants associated with vertical cup-to-disc ratio (VCDR). Genetic risk scores (GRS) examine the aggregate genetic effect of individual variants on a trait by combining these separate genetic variants into a single measure. The purpose of this study was to construct GRS for VCDR and to determine whether the GRS are associated with VCDR and whether the GRS increase the discriminatory ability for primary open-angle glaucoma (POAG) in a Latino population.

DESIGN

Population-based genetic association study.

PARTICIPANTS

A total of 4018 Latino participants recruited from Los Angeles.

METHODS

Weighted and unweighted GRS were constructed using 68 previously reported VCDR single nucleotide polymorphisms (SNPs), as well as SNPs from our own genome-wide association data. Linear and logistic regression analyses examined the associations of GRS with VCDR and POAG, respectively. To evaluate the discriminatory ability of the GRS for POAG, we conducted receiver operating characteristic (ROC) analyses.

MAIN OUTCOME MEASURES

The relationship between GRS and VCDR in Latinos.

RESULTS

The GRS were associated significantly with VCDR (P < 0.0001), after adjusting for age, gender, central corneal thickness, intraocular pressure, and education. The weighted GRS explained an additional 2.74% of the variation in VCDR. Adding the weighted GRS derived from previously reported SNPs resulted in a moderate improvement in the discriminatory ability for POAG during ROC analyses, yielding an area under the ROC curve (AUC) of 0.735 (95% CI, 0.701-0.768). When our own SNPs were used, the AUC increased significantly to 0.809 (95% CI, 0.781-0.837; P < 0.0001). We obtained similar results for the unweighted GRS.

CONCLUSIONS

To our knowledge, we identified a novel association between GRS and VCDR and its improvement in the discriminatory ability of POAG in a Latino population.

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.