1. Please check and comment entries here.
Table of Contents

    Topic review

    Glaucoma Heritability

    View times: 14
    Submitted by: Thomas Ciulla


    Glaucoma, a leading cause of irreversible blindness globally, is a degenerative optic neuropathy characterized by progressive visual field defects corresponding to retinal ganglion cell and retinal nerve fiber layer (RNFL) degeneration.

    1. Introduction

    Glaucoma, a leading cause of irreversible blindness globally, is a degenerative optic neuropathy characterized by progressive visual field defects corresponding to retinal ganglion cell and retinal nerve fiber layer (RNFL) degeneration [1][2]. Primary open-angle glaucoma (POAG) represents the most common form of glaucoma worldwide, responsible for nearly three-quarters of all cases of glaucoma [1]. POAG is a form of disease characterized by an anatomically open angle; in contrast to primary angle-closure glaucoma (PACG) where anatomic changes obstruct aqueous humor outflow from the anterior chamber, in POAG, drainage of aqueous humor out of the anterior chamber is obstructed despite an anatomy that should allow for aqueous drainage [3]. An elevated intraocular pressure (IOP) represents the major and only currently approved modifiable risk factor for glaucomatous disease onset and disease progression [4]. However, not all patients with glaucoma have elevated IOP, and not all patients with elevated IOP develop glaucoma. Therefore, the pathogenesis of glaucoma is still not fully understood and likely differs among persons and populations.
    As researchers have tried to explore and better understand the multifactorial nature of glaucoma, many studies have sought a genetic explanation for the disease. Genetic and genomic studies of single-gene variants, whole-exome sequencing, genome-wide association studies (GWAS), and genetic/polygenic risk scoring have all been used to describe genetic contributions to glaucoma and to evaluate the role of genetics as a risk factor for disease [5]. In addition to these studies, scientists have tried to evaluate the heritability of glaucoma.
    Heritability, represented by h2, is a complex statistical measurement ranging from zero to one that describes what percent of variability in a trait is due to genetic factors versus environmental factors [6]. For example, characteristics such as language and political ideology have a heritability close to zero—indicating that these characteristics are almost entirely influenced by environmental factors with no genetic influence. Meanwhile, traits caused by single-gene mutations, such as phenylketonuria, have a heritability close to one, indicating that variability in this trait is almost entirely the result of genetic influence [7]. Glaucoma represents a complex, multifactorial trait with a quantifiable genetic influence, and therefore h2 estimates for glaucoma fall between zero and one [8].
    In 2019, Asefa et al. conducted a systematic review and meta-analysis of glaucoma heritability, and the heritability of glaucoma-related endophenotypes. For glaucoma specifically, their search yielded six studies—three specifically for POAG, two for unspecified glaucoma, and one for PACG—and their pooled h2 for POAG ranged from 0.17 to 0.81 [8]. Meanwhile, their pooled h2 estimates for POAG endophenotypes included h2 estimates for IOP (h2 = 0.43), anterior chamber size (h2 = 0.67), central corneal thickness (CCT) (h2 = 0.81), cup-to-disc ratio (CDR) (h2 = 0.56), disc size (h2 = 0.61), cup size (h2 = 0.58), corneal hysteresis (h2 = 0.40), RNFL thickness (h2 = 0.73), cup shape (h2 = 0.62), and peripapillary atrophy (h2 = 0.73) [8].
    Since glaucoma is a multifactorial disease, heritability estimates of specific glaucoma endophenotypes and risk factors may be more important clues for researchers as to which endophenotypes, and in turn molecular mechanisms, are important to the pathogenesis and progression of disease. While these estimates for glaucoma and endophenotype heritability are generally strong, the majority of genetic variants underlying these traits are still unknown. Additionally, the molecular etiology of these endophenotypes, and their role in the pathogenesis of glaucoma, is similarly unclear. For instance, it has been suggested that IOP and vertical CDR (VCDR) may represent the true POAG endophenotypes, while CCT may not be physiologically significant, but the molecular mechanisms underlying genetic contributions to these traits are not well understood [9][10].

    2. Intraocular Pressure

    Elevated intraocular pressure represents the only currently approved modifiable risk factor for POAG. While the exact molecular mechanisms underlying elevated IOP and its relationship to glaucoma are not fully understood, it is well known that elevated IOP is strongly associated with both structural and functional glaucomatous change [3][11]. Currently, this relationship is described by mechanical strain that elevated IOP places on the lamina cribrosa and adjacent optic nerve tissues, ultimately resulting in structural damage that disrupts axonal transport [3][12]. It has also been suggested that there may be a pressure-induced metabolic stress effect, and even circulatory, immune, and oxidative stress mechanisms associated with glaucoma and related to IOP [3][13][14].
    Interestingly, elevated IOP has shown a relatively strong heritability estimate. When Asefa et al. conducted a systematic review and meta-analysis for the heritability of glaucoma, they also evaluated pooled h2 estimates for endophenotypes, including IOP, which they estimated at h2 = 0.43 (95% confidence interval (CI) 0.38–0.48). Importantly, in their analysis, IOP was the most commonly reported endophenotype, reported in 41 different studies. This large sample also allowed for subgroup analysis, showing that there was heterogeneity between age and ethnicity. Findings showing different heritability estimates between White Europeans (h2 = 0.47, 95% CI 0.37–0.57, p = 0.006), East Asians (h2 = 0.49, 95% CI 0.40–0.59), and Other Populations (h2 = 0.30, 95% CI 0.22–0.39) are particularly interesting given suspected genetic contributions to disease and genetic differences among different racial groups [5][8]. Additionally, it should be noted that IOP has also shown a genetic correlation with POAG (odds ratio 1.18, 95% CI 1.14–1.21, p = 1.8 × 10−27) [15].
    To date, a wide variety of genes have been associated with elevated IOP, including, but not limited to: ABCA1, ABO, ADAMTS8, ADAMTS17, ADAMTS18-NUDT7, AFAP1, ANGPT1, ANTXR1, ARHGEF12, ARID5B, ATXN2, CAV1-CAV2, CDKN2B-AS1, CELF1, CYP26A1-MYOF, FAM125B, FNDC3B, FOXC1, FOXP1, GAS7, GLCCI1-ICA1, GLIS3, GMDS, HIVEP3, INCA1, LMX1B, LOC171391, MADD, MIR548F3, MYBPC3, NDUFS3, NR1H3, PDDC1, PKHD1, PTPRJ, RAPSN, RPLP2-PNPLA2, SIX1/SIX6, SEPT9, SP11, TFEC-TES, TMCO1, and TXNRD2 [5][16]. It should be noted that all of these genes have not also been independently associated with POAG. Similarly, molecular mechanisms of these genes are generally not well understood (Table 1).
    Table 1. Selection of genes associated with IOP and suggested molecular mechanisms.
    Gene GWAS Identified Suggested Mechanism Reference
    LMX1B Gao et al., 2018 [16] Altered Anterior Segment Development; Altered Aqueous Humor Dynamics Pressman et al., 2000 [17]
    MADD Gao et al., 2018 [16] TNF-a-Mediated Microglial Activation Tezel et al., 2000 [18]
    NR1H3 Gao et al., 2018 [16] Alteration of IOP via-ABCA1-Regulated Aqueous Humor Dynamic Alterations Hu et al., 2020 [19]
    SEPT9 Gao et al., 2018 [16] Cytoskeletal Alterations Ou et al., 1998 [20]
    In 2018, Gao et al. conducted the largest known GWAS for IOP using patients of European Descent (ED) from the UK Biobank (n = 115,486). They identified 671 significant genetic variants associated with IOP, mapping to 149 loci, significantly expanding the known literature and replicating findings from previous studies. Importantly, they specifically identified four significant novel genes: LMX1B, MADD, NR1H3, and SEPT9 [16].
    LMX1B is a gene that is directly implicated in the development of serotonergic (5-HT) neurons [21]. According to Ding et al., LMX1B is a LIM homeodomain-containing gene, which they suggest is a key intermediate factor necessary for the cascade of specification and differentiation of 5-HT neurons [21]. In mouse models, the LMX1B transcription factor has been associated with the regulation of anterior segment development, with LMX1B transcripts being identified in the ciliary body and trabecular meshwork (Table 2) [17]. In fact, Pressman et al. described that mice with targeted LMX1B mutations have structural deficiencies in the ciliary body, corneal stroma, and iris [17]. If LMX1B mutations are associated with alterations in the ciliary body and anterior chamber structures, it is logical to consider that aqueous humor dynamics may be altered, leading to elevated IOP and, potentially, POAG, though further evaluation is surely warranted (Figure 1).
    Figure 1. Schematic representation of the proposed location of the genes’ mechanism of action (adapted with permission from: Prada D, Harris A, Guidoboni G, Siesky B, Huang AM, Arciero J. Autoregulation and neurovascular coupling in the optic nerve head. Survey of Ophthalmology 1 March 2016; 61(2): 164–186).
    Table 2. List of IOP and POAG genes under investigation.
    Gene Study Environment
    ABCA1 Hu et al., 2020 [19] Ex vivo Human Model
    ABCA1 Li et al., 2018 [22] In vivo Mouse Model
    ABCA1 Yeghiazaryan et al., 2005 [23] In vivo Human Model
    ATOH7 Song et al., 2015 [24] In vivo Rat Model
    ATOH7 Miesfeld et al., 2020 [25] In vivo Mouse Model
    ELN Gelman et al., 2010 [26] Ex vivo Mouse Model
    LMX1B Cross et al., 2014 [27] In vivo Mouse Model
    LMX1B Pressman et al., 2000 [17] In vivo Mouse Model
    MADD Schievella et al., 1997 [28] In vitro Yeast Interaction Trap
    NR1H3 Wang et al., 2002 [29] In vivo Mouse Model
    NR1H3 Yang et al., 2014 [30] In vivo Mouse Model
    NR1H3 Zheng et al., 2015 [31] In vivo Mouse Model
    NR1H3 Song et al., 2019 [32] Ex vivo Mouse Model
    SEPT9 Ghossoub et al., 2013 [33] Ex vivo Human Model
    SEPT9 Ou et al., 1998 [20] In vivo Guinea Pig Model
    It should be noted that LMX1B mutations have been associated with families affected by nail–patella syndrome (NPS), a rare autosomal dominant developmental disorder with variable presentation, and POAG, suggesting an explanation for the association between NPS and glaucoma (approximately 1/3 of patients with NPS develop glaucoma), though not explaining a mechanism for causation [34][35][36]. In 2018, Choquet et al. suggested that the variability in glaucoma phenotypes associated with LMX1B may be explained by different mutations. By inducing different LMX1B mutations in a murine model, they found that different LMX1B mutations can result in elevated IOP and glaucomatous damage independent of developmental defects, suggesting that LMX1B is a significant glaucoma-susceptibility locus even in the absence of developmental abnormalities [37]. Still, further research is necessary to better understand the role of LMX1B in glaucoma pathogenesis.
    NR1H3 is a nuclear receptor also referred to as the liver X receptor α. Liver X receptors regulate lipid homeostasis, with expression in both the liver and the brain, and play a key role neurodegenerative disease [29]. According to Yang et al., liver X receptor α may also play an anti-inflammatory role, as demonstrated by a mouse model of autoimmune uveitis (Table 2) [30]. In fact, they suggested that liver X receptor α may inhibit the NF-kB signaling pathway, resulting in the reduction in ocular inflammation. These findings reflect those of Zheng at al., which suggest that the activation of liver X receptor β induces a protective effect on the retina perhaps by decreasing retinal amyloid β formation and by inhibiting the NF-kB signaling pathway [31]. While Yang et al. investigated this anti-inflammatory effect in the context of autoimmune uveitis therapeutics, recent research has demonstrated the utility of the anti-inflammatory effect of liver X receptors in relation to glaucoma.
    In 2019, Song et al. demonstrated that liver X receptors α and β were expressed in the retina and optic nerve (Figure 1) [32]. They also demonstrated that there was an accumulation of amyloid deposition in the retinal ganglion cells of liver X receptor β knockout mice, and that there was activation of microglia and loss of aquaporin 4 (a biomarker of optic nerve function) in the optic nerve. Ultimately, they reported that the loss of liver X receptor β results in retinal ganglion cell loss and optic nerve degeneration (Table 2) [32]. Interestingly, they noted that aquaporin 4 expression decreased without the presence of aquaporin 4 antibodies, suggesting that optic nerve degeneration in this model is caused by the loss of aquaporin 4 expression [32]. While these findings are encouraging, they do not directly explain the association between NR1H3 and IOP as suggested by the 2018 GWAS from Gao et al. [16].
    The association between NR1H3 and IOP may actually be explained by ABCA1, a POAG-susceptibility locus associated with ED populations as well as Asian populations [5][38][39]. ABCA1, the ATP-binding cassette transporter A1, is a transmembrane transporter involved in cholesterol homeostasis through the formation of high-density lipoprotein [40]. Importantly, ABCA1 is regulated by the liver X receptor [41]. While ABCA1 has been implicated in retinal inflammation and retinal ganglion cell apoptosis, recent research has also demonstrated that ABCA1 has a role in regulating IOP through aqueous humor dynamics [19][22]. In a 2020 study, Hu et al. proposed a pathophysiologic mechanism for the role of ABCA1 in regulating IOP, and an alternate function to its role in cholesterol homeostasis. They suggested that ABCA1 changes Cav1 expression, resulting in a decreased inhibitory influence of Cav1 on eNos, resulting in the increased nitric oxide production by eNOS and, ultimately, reduced IOP (Table 2) [19]. Additionally, they confirmed that ABCA1 had higher expression in the trabecular meshwork of POAG patients versus controls, confirming previous work suggesting the enhanced expression of ABCA1 as a diagnostic marker for glaucoma [19][23]. Therefore, it is possible that variants in ABCA1 and NR1H3 result in altered IOP and aqueous outflow regulation.
    The third most significant novel gene identified by Gao et al. in their 2018 GWAS was MADD, a gene from the MAP kinase activating death domain [16]. MADD plays a role in the tumor necrosis factor α (TNF-a) signaling pathway and, therefore, apoptotic signal transduction, by way of mitogen-activated protein kinase (MAPK) [28]. MADD has been found to be upregulated in the retinas of POAG patients (Figure 1) [28][42]. This finding reflects those of increased TNF-a and TNF-a receptor 1 in glaucomatous optic nerve heads as measured by immunohistochemistry, suggesting a role for this signaling pathway in the neurodegenerative pathophysiology for POAG [42][43][44].
    Currently, there are no studies specifically linking MADD expression and IOP. However, TNF-a and TNF-a receptor expression increased by nearly eightfold has been demonstrated in rats with elevated IOP compared to other animal models [45]. Additionally, it has been suggested that TNF-a may indirectly mediate the cytotoxicity of elevated IOP via-microglial activation, and that stimulated glial cells may secrete TNF-a and nitric oxide to initially attenuate the effects of elevated IOP—despite the detrimental downstream effects of microglial activation such as the loss of aquaporin 4 [18][46][47][48]. Therefore, there may be an association between MADD expression and IOP, as described by Gao et al. [16]. Still, a 2010 analysis of TNF-a concentrations in the aqueous humor of glaucoma patients found a significant association between TNF-a concentrations and glaucoma patients (mean concentration 15.9 ± 3.6 pg/mL, range 1.7–57.6 pg/mL) but was unable to find a relationship between IOP and TNF-a concentration (r = 0.163) [49]. Therefore, the relationship between MADD, and by extension TNF-a, and IOP remains unclear.
    The final significant novel gene identified by Gao et al. was SEPT9, a septin protein [16]. Septins are a group of GTP-binding proteins that function in cell division and cytoskeleton formation; they have also been associated with neurodegenerative diseases [50]. In 2013, Ghossoub et al. described the role of septins in retinal pigment epithelial cells, noting how SEPT2 forms a complex with SEPT7 and SEPT9, with the complex controlling ciliogenesis and, more specifically, the elongation of the primary cilium in cells of the retinal pigment epithelium (Table 2) [33].
    The molecular mechanism linking septin activity and intraocular pressure has not been elucidated. It should be noted, though, that microtubule-associated protein 4 appears to play an antagonistic role with the activity of the SEPT2/SEPT7/SEPT9 complex in retinal pigment epithelial cells (Figure 1) [33]. Although not the same protein, one previous immunohistochemical study have shown that acutely increased IOP leads to a decrease in microtubule-associated protein 1 expression and a resulting decrease in the number of microtubules in optic nerve tissues, suggesting a possible mechanism by which cytoskeletal structures may be linked to IOP changes [20].
    While these four genes were identified as the most significant by Gao et al. in their 2018 GWAS, it is important to note that these genes are not conclusively the most important genetic markers of IOP. In fact, when Gao et al. calculated a heritability estimate for the UK Biobank sample limited to only significant genetic variants (n = 671), the significant variants only explained 7.2% of variance [16]. Additionally, when evaluated for pleiotropy, the 671 significant genetic variants also directly matched SNPs in the GWAS Catalog associated with PACG; age-related macular degeneration; other ocular characteristics such as central corneal thickness, axial length, and iris characteristics; and other systemic traits and conditions such as blood pressure, body mass index, type 2 diabetes mellitus, cancers, and immune disorders [16].
    Additionally, it is important to note, that IOP is not a perfect indicator for glaucoma, and therefore it cannot fully explain the heritability, nor the underlying molecular mechanisms, of all cases of POAG. In fact, a large subset of POAG patients do not have elevated IOP. These patients are referred to as having normal-tension glaucoma (NTG)—a unique classification of glaucoma where glaucomatous damage occurs independent of elevated IOP [51]. In fact, Park et al. even suggested that the alteration of LMX1B transcription factor function may lead to glaucomatous damage independent of IOP alterations, as they noted that single-nucleotide polymorphisms (SNPs) of LMX1B were also associated with NTG—and highlighting the danger of defining and evaluating POAG singularly through a lens of elevated IOP [36].

    3. Vertical Cup-to-Disc Ratio

    In addition to IOP, CDR has been similarly suggested to be a “true” POAG endophenotype [9]. VCDR is a reflection of the loss of retinal ganglion cells and axons, which results in the thinning of the neuroretinal rim and is represented clinically by alterations to the VCDR [52]. Similar to IOP, VCDR is an imperfect indicator of disease. According to Tatham et al., VCDR has limited ability to detect glaucoma given the wide variability of VCDR in healthy eyes [53]. Additionally, even though VCDR reflects retinal ganglion cell loss, VCDR is not a strong method for estimating quantitative measurements of retinal ganglion cell loss, largely due to wide variability of optic nerve head parameters [53]. In fact, in a large population-based study of healthy eyes (n = 6616), Hopley et al. found that the wide variability in normal CDR is influenced at least in part by the size of the optic disc [54]. Still, according to a large systematic review from Hollands et al., increased CDR is a key risk factor for POAG, as the likelihood for POAG increases as CDR increases. [55].
    In their 2019 systematic review and meta-analysis, Asefa et al. calculated a heritability estimate for CDR at h2 = 0.56 (0.44–0.68) [8]. Importantly, they suggested that the wide variability in heritability estimates of CDR may be a reflection of the variability in cup size in the general population. Still, CDR shows relatively high heritability as a trait.
    Genes significantly associated with CDR include, but are not limited to: ABCA1, ABG, ADAMTS8, ASAP1, ASB7, ATOH7, ATOH7-PBLD, BMP2, CARD10, CDC7-TGFBR3, CDKN2B, CDKN2B-CDKN2BAS, CHEK2, COL8A1, CRISPLD1, DCLK1, DGKB, DUSP1, ELN, ENO4, EXOC2, F5, FAM101A, GAS7, HSF2, PDZD2, PLCE1, PSCA, RARB, RERE, RPAP3, RPE65, RREB1, SALL1, SCYL1, SIX1, SIX6, SSSCA1, TMTC2, and VCAN [56][57][58][59][60][61][62]. Several of these genes, including ABCA1, ABG, AFAP1, CAV1, GAS7, and LMX1B have been independently associated with both IOP and CDR [59][61]. Importantly, however, not all of these genes have also been independently associated with POAG.
    In their 2017 GWAS and meta-analysis, Springelkamp et al. identified a variety of novel genes influencing CDR and other POAG endophenotypes. In order to better understand the molecular mechanisms behind these genes, they conducted a pathway analysis using DEPICT, a genetic association tool intended to predict gene function [63]. Their findings were diverse; for example, they noted that genes associated with CDR, cup area and disc area were associated with biologic pathways associated with metabolic processes [59]. Given redundancy between significant pathways, they clustered pathways into a series of meta-pathways, noting that gene sets centered on pathways involving cell differentiation and other developmental pathways [59].
    Interestingly, the meta-pathways Springelkamp et al. developed singled out pathways associated with abnormal fat cell and liver morphologies, noting that these pathways centered on ABCA1 [59]. As described above, ABCA1 is a transmembrane transporter with a role in cholesterol processing that has recently been linked to IOP regulation and aqueous humor dynamics (Table 3) [19]. So, while this recent development explains the linkage between ABCA1, IOP, and the meta-pathways described by Springelkamp et al., the link to CDR is still unclear.
    Table 3. Selection of genes associated with CDR and suggested molecular mechanisms.
    Gene GWAS Identified Suggested Mechanism Reference
    ABCA1 Springelkamp et al., 2017 [59] Normal Function and Cell Death of Retinal Ganglion Cells Chen et al., 2014 [64]
    ELN Han et al., 2021 [61] Alteration to Normal Function of Elastin, Leading to Optic Nerve Head Degeneration Gelman et al., 2010 [26]
    ASAP1 Alipanahi et al., 2021 [65] Glial Cell-Mediated Retinal Ganglion Cell Loss García-Bermúdez et al., 2021 [66]
    ATOH7 Nannini et al., 2018 [60] Alteration to Müller Cel Differentiation and Retinal Ganglion Cell Genesis Miesfeld et al., 2020 [25]
    In 2014, Chen et al. found that ABCA1 was highly expressed in the retinal ganglion cell layer and is significantly higher in POAG patients. They suggested that this high expression indicates that ABCA1 plays a role in the normal functioning and cell death of retinal ganglion cells [64]. Given that CDR is a reflection of retinal ganglion cell death, it is likely that, in addition to its suggested role in regulating IOP and aqueous humor dynamics, ABCA1 is linked to CDR through a molecular mechanism involving retinal ganglion cell death (Figure 1). It is important to note that in 2018 Li et al. demonstrated that the liver X receptor/ABCA1 pathway may actually be protective against retinal ganglion cell apoptosis (Table 2) [22]. Their model, however, noted reduced ABCA1 retinal expression in response to acute IOP elevation, which may not be reflective of the role of ABCA1 protein in POAG and its link to CDR. Additionally, it is possible that associations between CDR and ABCA1 may be reflective of IOP-mediated retinal ganglion cell death. Therefore, further research is surely warranted to better understand the true mechanistic relationship between these closely intertwined variables.
    In 2021, Han et al. conducted the largest known GWAS for VCDR and vertical disc diameter (VDD) using by utilizing artificial intelligence (AI) in the form of a convolutional neural network model. Using this technology they were able to analyze VCDR and VDD from 282,100 images from the UK Biobank and another independent study and then were able to systematically use these studies to conduct a stronger GWAS of optic nerve head parameters [61]. Using this data, they compared heritability estimates for VCDR as calculated from clinician gradings of CDR versus calculated from the AI-based GWAS. They found that the SNP-based heritability estimates were nearly 50% higher (VCDR h2 = 0.22 vs. 0.35), indicating a future role for AI in heritability studies. In fact, they suggested that AI techniques may lead to higher heritability estimates in twin-studies [61].
    In addition to stronger heritability studies, the AI-assisted GWAS more than doubled the number of identified genetic loci for VCDR and VDD, similarly uncovering new mechanistic pathways. For example, Han et al. described new VCDR genes with possible links to retinal ganglion cell biology. Among these, a missense variant in the elastin gene, ELN, was particularly noteworthy [61].
    Elastin is an extracellular matrix protein that has a significant role in tissue strength and elasticity [67]. In the eye, elastin has been found extensively throughout the choroid, conjunctiva, meninges, muscle tendons and, most densely, in the sclera surrounding the optic nerve head (Figure 1) [26]. Given its function, elastin is a key component to normal ocular physiology as the sclera must be able to withstand the constant mechanical strain and stress induced by normal IOP and ocular blood flow pulsations (Table 2) [26]. Therefore, an alteration to the function of elastin and its ability to permit elastic structural movement and recovery, could have a dramatic degenerative effect. In the case of ELN and VCDR, it is logical that alterations to the connective tissue of the lamina cribrosa could predispose a patient to glaucomatous damage and, as a result, increased VCDR.
    Interestingly, in 2021, Alipanahi et al. ran a similar study using AI to enhance the power of a GWAS. Using a machine learning model based on 81,830 ophthalmologist-labeled fundus images to phenotype VCDR, they predicted VCDR on 65,680 participants from the UK Biobank. Then, using the machine learning-generated VCDR, they performed a GWAS to evaluate for genetic associations with VCDR. They uncovered 93 novel loci associated with VCDR [65].
    Among the novel genetic loci they uncovered was a SNP near ASAP1, a locus previously associated with glioma [65][68]. Recent research has demonstrated that glaucomatous retinal ganglion cell loss, and therefore clinical measurements of VCDR, may in some cases be mediated by glial cells [66]. For example, the interplay of microglial activation and nitric oxide, as described previously, can be a beneficial response to elevations in IOP but may also lead to the downstream loss of aquaporin channels, leading to optic nerve degeneration [46]. In the case of ASAP1, Alipanahi et al. suggested that this association may indicate a mechanism by which glial cells mediate VCDR (Table 2) [65].
    An important example of a glial cell-mediated mechanism that would lead to glaucomatous damage reflected clinically by VCDR is the gene ATOH7, which has been significantly associated with VCDR [56][57][69]. The ATHO7 protein, Atonal BHLH Transcription Factor 7, is a key transcription factor for the differentiation of Müller cells, the key retinal glial cell [24][70]. According to Song et al., ATOH7 may promote the differentiation of retinal Müller cell-derived stem cells into retinal ganglion cells (Table 2) [24]. In fact, recent studies have shown that ATOH7 is a critical regulator of retinal ganglion cell genesis (Figure 1). In 2020, Miesfeld et al. demonstrated that retinal ganglion cell genesis was strictly dependent on ATOH7 expression (Table 2) [25]. Using a mouse model, they demonstrated how the number of adult retinal ganglion cells was directly correlated with ATOH7 expression. Additionally, they noted how ATOH7 mutants with decreased retinal ganglion cells experienced a variety of downstream effects and secondary malformations, including alterations to the retinal vasculature and optic nerve head parameters [25]. These findings demonstrate that the association between ATOH7 and VCDR may be mediated by the effect of the ATOH7 protein on glial cell differentiation into retinal ganglion cells, indicating one mechanism by which retinal glial cell alterations secondary to genetic alterations may predispose patients to alterations in VCDR.

    The entry is from 10.3390/genes12081135


    1. Quigley, H.A.; Broman, A.T. The number of people with glaucoma worldwide in 2010 and 2020. Br. J. Ophthalmol. 2006, 90, 262–267.
    2. Tham, Y.C.; Li, X.; Wong, T.Y.; Quigley, H.A.; Aung, T.; Cheng, C.Y. Global prevalence of glaucoma and projections of glaucoma burden through 2040: A systematic review and meta-analysis. Ophthalmology 2014, 121, 2081–2090.
    3. Weinreb, R.N.; Aung, T.; Medeiros, F.A. The Pathophysiology and Treatment of Glaucoma: A Review. JAMA 2014, 311, 1901–1911.
    4. Boland, M.V.; Ervin, A.M.; Friedman, D.S.; Jampel, H.D.; Hawkins, B.S.; Vollenweider, D.; Chelladurai, Y.; Ward, D.; Suarez-Cuervo, C.; Robinson, K.A. Comparative effectiveness of treatments for open-angle glaucoma: A systematic review for the U.S. Preventive Services Task Force. Ann. Intern. Med. 2013, 158, 271–279.
    5. Zukerman, R.; Harris, A.; Verticchio Vercellin, A.; Siesky, B.; Pasquale, L.R.; Ciulla, T.A. Molecular Genetics of Glaucoma: Subtype and Ethnicity Considerations. Genes 2021, 12, 55.
    6. Tenesa, A.; Haley, C.S. The heritability of human disease: Estimation, uses and abuses. Nat. Rev. Genet. 2013, 14, 139–149.
    7. Haworth, C.M.; Davis, O.S. From observational to dynamic genetics. Front. Genet. 2014, 5, 6.
    8. Asefa, N.G.; Neustaeter, A.; Jansonius, N.M.; Snieder, H. Heritability of glaucoma and glaucoma-related endophenotypes: Systematic review and meta-analysis. Surv. Ophthalmol. 2019, 64, 835–851.
    9. Charlesworth, J.; Kramer, P.L.; Dyer, T.; Diego, V.; Samples, J.R.; Craig, J.E.; Mackey, D.A.; Hewitt, A.W.; Blangero, J.; Wirtz, M.K. The path to open-angle glaucoma gene discovery: Endophenotypic status of intraocular pressure, cup-to-disc ratio, and central corneal thickness. Invest. Ophthalmol. Vis. Sci. 2010, 51, 3509–3514.
    10. Bonnemaijer, P.W.M.; Leeuwen, E.M.V.; Iglesias, A.I.; Gharahkhani, P.; Vitart, V.; Khawaja, A.P.; Simcoe, M.; Höhn, R.; Cree, A.J.; Igo, R.P.; et al. Multi-trait genome-wide association study identifies new loci associated with optic disc parameters. Commun. Biol. 2019, 2, 435.
    11. Medeiros, F.A.; Alencar, L.M.; Zangwill, L.M.; Sample, P.A.; Weinreb, R.N. The Relationship between intraocular pressure and progressive retinal nerve fiber layer loss in glaucoma. Ophthalmology 2009, 116, 1125–1133.
    12. Quigley, H.A.; Addicks, E.M.; Green, W.R.; Maumenee, A.E. Optic nerve damage in human glaucoma. II. The site of injury and susceptibility to damage. Arch. Ophthalmol. 1981, 99, 635–649.
    13. Ju, W.K.; Kim, K.Y.; Lindsey, J.D.; Angert, M.; Duong-Polk, K.X.; Scott, R.T.; Kim, J.J.; Kukhmazov, I.; Ellisman, M.H.; Perkins, G.A.; et al. Intraocular pressure elevation induces mitochondrial fission and triggers OPA1 release in glaucomatous optic nerve. Invest. Ophthalmol. Vis. Sci. 2008, 49, 4903–4911.
    14. Moore, D.; Harris, A.; Wudunn, D.; Kheradiya, N.; Siesky, B. Dysfunctional regulation of ocular blood flow: A risk factor for glaucoma? Clin. Ophthalmol. 2008, 2, 849–861.
    15. Aschard, H.; Kang, J.H.; Iglesias, A.I.; Hysi, P.; Cooke Bailey, J.N.; Khawaja, A.P.; Allingham, R.R.; Ashley-Koch, A.; Lee, R.K.; Moroi, S.E.; et al. Genetic correlations between intraocular pressure, blood pressure and primary open-angle glaucoma: A multi-cohort analysis. Eur. J. Hum. Genet. 2017, 25, 1261–1267.
    16. Gao, X.R.; Huang, H.; Nannini, D.R.; Fan, F.; Kim, H. Genome-wide association analyses identify new loci influencing intraocular pressure. Hum. Mol. Genet. 2018, 27, 2205–2213.
    17. Pressman, C.L.; Chen, H.; Johnson, R.L. LMX1B, a LIM homeodomain class transcription factor, is necessary for normal development of multiple tissues in the anterior segment of the murine eye. Genesis 2000, 26, 15–25.
    18. Tezel, G.; Wax, M.B. Increased production of tumor necrosis factor-α by glial cells exposed to simulated ischemia or elevated hydrostatic pressure induces apoptosis in cocultured retinal ganglion cells. J. Neurosci. 2000, 20, 8693–8700.
    19. Hu, C.; Niu, L.; Li, L.; Song, M.; Zhang, Y.; Lei, Y.; Chen, Y.; Sun, X. ABCA1 Regulates IOP by Modulating Cav1/eNOS/NO Signaling Pathway. Invest. Ophthalmol. Vis. Sci. 2020, 61, 33.
    20. Ou, B.; Ohno, S.; Tsukahara, S. Ultrastructural changes and immunocytochemical localization of microtubule-associated protein 1 in guinea pig optic nerves after acute increase in intraocular pressure. Invest. Ophthalmol. Vis. Sci. 1998, 39, 963–971.
    21. Ding, Y.Q.; Marklund, U.; Yuan, W.; Yin, J.; Wegman, L.; Ericson, J.; Deneris, E.; Johnson, R.L.; Chen, Z.F. Lmx1b is essential for the development of serotonergic neurons. Nat. Neurosci. 2003, 6, 933–938.
    22. Li, L.; Xu, L.; Chen, W.; Li, X.; Xia, Q.; Zheng, L.; Duan, Q.; Zhang, H.; Zhao, Y. Reduced Annexin A1 Secretion by ABCA1 Causes Retinal Inflammation and Ganglion Cell Apoptosis in a Murine Glaucoma Model. Front. Cell Neurosci. 2018, 12, 347.
    23. Yeghiazaryan, K.; Flammer, J.; Wunderlich, K.; Schild, H.H.; Orgul, S.; Golubnitschaja, O. An enhanced expression of ABC 1 transporter in circulating leukocytes as a potential molecular marker for the diagnostics of glaucoma. Amino Acids 2005, 28, 207–211.
    24. Song, W.T.; Zhang, X.Y.; Xia, X.B. Atoh7 promotes the differentiation of Müller cells-derived retinal stem cells into retinal ganglion cells in a rat model of glaucoma. Exp. Biol. Med. 2015, 240, 682–690.
    25. Miesfeld, J.B.; Ghiasvand, N.M.; Marsh-Armstrong, B.; Marsh-Armstrong, N.; Miller, E.B.; Zhang, P.; Manna, S.K.; Zawadzki, R.J.; Brown, N.L.; Glaser, T. The Atoh7 remote enhancer provides transcriptional robustness during retinal ganglion cell development. Proc. Natl. Acad. Sci. USA 2020, 117, 21690–21700.
    26. Gelman, S.; Cone, F.E.; Pease, M.E.; Nguyen, T.D.; Myers, K.; Quigley, H.A. The presence and distribution of elastin in the posterior and retrobulbar regions of the mouse eye. Exp. Eye Res. 2010, 90, 210–215.
    27. Cross, S.H.; Macalinao, D.G.; McKie, L.; Rose, L.; Kearney, A.L.; Rainger, J.; Thaung, C.; Keighren, M.; Jadeja, S.; West, K.; et al. A Dominant-Negative Mutation of Mouse Lmx1b Causes Glaucoma and Is Semi-lethal via LBD1-Mediated Dimerisation. PLoS Genet 2014, 10, e1004359.
    28. Schievella, A.R.; Chen, J.H.; Graham, J.R.; Lin, L.L. MADD, a novel death domain protein that interacts with the type 1 tumor necrosis factor receptor and activates mitogen-activated protein kinase. J. Biol. Chem. 1997, 272, 12069–12075.
    29. Wang, L.; Schuster, G.U.; Hultenby, K.; Zhang, Q.; Andersson, S.; Gustafsson, J.A. Liver X receptors in the central nervous system: From lipid homeostasis to neuronal degeneration. Proc. Natl. Acad. Sci. USA 2002, 99, 13878–13883.
    30. Yang, H.; Zheng, S.; Qiu, Y.; Yang, Y.; Wang, C.; Yang, P.; Li, Q.; Lei, B. Activation of liver X receptor alleviates ocular inflammation in experimental autoimmune uveitis. Invest. Ophthalmol. Vis. Sci. 2014, 55, 2795–2804.
    31. Zheng, S.; Yang, H.; Chen, Z.; Zheng, C.; Lei, C.; Lei, B. Activation of Liver X Receptor Protects Inner Retinal Damage Induced by N-methyl-D-aspartate. Invest. Ophthalmol. Vis. Sci. 2015, 56, 1168–1180.
    32. Song, X.Y.; Wu, W.F.; Gabbi, C.; Dai, Y.B.; So, M.; Chaurasiya, S.P.; Wang, L.; Warner, M.; Gustafsson, J. Retinal and optic nerve degeneration in liver X receptor β knockout mice. Proc. Natl. Acad. Sci. USA 2019, 116, 16507–16512.
    33. Ghossoub, R.; Hu, Q.; Failler, M.; Rouyez, M.C.; Spitzbarth, B.; Mostowy, S.; Wolfrum, U.; Saunier, S.; Cossart, P.; Jamesnelson, W.; et al. Septins 2, 7 and 9 and MAP4 colocalize along the axoneme in the primary cilium and control ciliary length. J. Cell Sci. 2013, 126, 2583–2594.
    34. Mimiwati, Z.; Mackey, D.A.; Craig, J.E.; Mackinnon, J.R.; Rait, J.L.; Liebelt, J.E.; Ayala-Lugo, R.; Vollrath, D.; Richards, J.E. Nail-patella syndrome and its association with glaucoma: A review of eight families. Br. J. Ophthalmol. 2006, 90, 1505–1509.
    35. Vollrath, D.; Jaramillo-Babb, V.L.; Clough, M.V.; McIntosh, I.; Scott, K.M.; Lichter, P.R.; Richards, J.E. Loss-of-function mutations in the LIM-homeodomain gene, LMX1B, in nail-patella syndrome. Hum. Mol. Genet. 1998, 7, 1091–1098.
    36. Park, S.; Jamshidi, Y.; Vaideanu, D.; Bitner-Glindzicz, M.; Fraser, S.; Sowden, J.C. Genetic risk for primary open-angle glaucoma determined by LMX1B haplotypes. Invest. Ophthalmol. Vis. Sci. 2009, 50, 1522–1530.
    37. Choquet, H.; Paylakhi, S.; Kneeland, S.C.; Thai, K.K.; Hoffmann, T.J.; Yin, J.; Kvale, M.N.; Banda, Y.; Tolman, N.G.; Williams, P.A.; et al. A multiethnic genome-wide association study of primary open-angle glaucoma identifies novel risk loci. Nat. Commun. 2018, 9, 2278.
    38. Gharahkhani, P.; Burdon, K.P.; Fogarty, R.; Sharma, S.; Hewitt, A.W.; Martin, S.; Law, M.H.; Cremin, K.; Bailey, J.N.C.; Loomis, S.J.; et al. Common variants near ABCA1, AFAP1 and GMDS confer risk of primary open-angle glaucoma. Nat. Genet. 2014, 46, 1120–1125.
    39. Hysi, P.G.; Cheng, C.Y.; Springelkamp, H.; Macgregor, S.; Bailey, J.N.C.; Wojciechowski, R.; Vitart, V.; Nag, A.; Hewitt, A.W.; Höhn, R.; et al. Genome-wide analysis of multi-ancestry cohorts identifies new loci influencing intraocular pressure and susceptibility to glaucoma. Nat. Genet. 2014, 46, 1126–1130.
    40. Wang, S.; Smith, J.D. ABCA1 and nascent HDL biogenesis. Biofactors 2014, 40, 547–554.
    41. Hozoji, M.; Munehira, Y.; Ikeda, Y.; Makishima, M.; Matsuo, M.; Kioka, N.; Ueda, K. Direct interaction of nuclear liver X receptor-β with ABCA1 modulates cholesterol efflux. J. Biol. Chem. 2008, 283, 30057–30063.
    42. Yang, X.; Luo, C.; Cai, J.; Powell, D.W.; Yu, D.; Kuehn, M.H.; Tezel, G. Neurodegenerative and inflammatory pathway components linked to TNF-α/TNFR1 signaling in the glaucomatous human retina. Invest. Ophthalmol. Vis. Sci. 2011, 52, 8442–8454.
    43. Yan, X.; Tezel, G.; Wax, M.B.; Edward, D.P. Matrix metalloproteinases and tumor necrosis factor α in glaucomatous optic nerve head. Arch. Ophthalmol. 2000, 118, 666–673.
    44. Tezel, G.; Li, L.Y.; Patil, R.V.; Wax, M.B. TNF-α and TNF-α receptor-1 in the retina of normal and glaucomatous eyes. Invest. Ophthalmol. Vis. Sci. 2001, 42, 1787–1794.
    45. Yang, Z.; Quigley, H.A.; Pease, M.E.; Yang, Y.; Qian, J.; Valenta, D.; Zack, D.J. Changes in gene expression in experimental glaucoma and optic nerve transection: The equilibrium between protective and detrimental mechanisms. Investig. Ophthalmol. Vis. Sci. 2007, 48, 5539–5548.
    46. Yuan, L.; Neufeld, A.H. Activated microglia in the human glaucomatous optic nerve head. J. Neurosci. Res. 2001, 64, 523–532.
    47. Nakazawa, T.; Nakazawa, C.; Matsubara, A.; Noda, K.; Hisatomi, T.; She, H.; Michaud, N.; Hafezi-Moghadam, A.; Miller, J.W.; Benowitz, L.I. Tumor necrosis factor-α mediates oligodendrocyte death and delayed retinal ganglion cell loss in a mouse model of glaucoma. J. Neurosci. 2006, 26, 12633–12641.
    48. Agarwal, R.; Agarwal, P. Glaucomatous neurodegeneration: An eye on tumor necrosis factor-α. Indian J. Ophthalmol. 2012, 60, 255–261.
    49. Sawada, H.; Fukuchi, T.; Tanaka, T.; Abe, H. Tumor Necrosis Factor-α Concentrations in the Aqueous Humor of Patients with Glaucoma. Invest. Ophthalmol. Vis. Sci. 2010, 51, 903–906.
    50. Bridges, A.A.; Gladfelter, A.S. Septin Form and Function at the Cell Cortex. J. Biol. Chem. 2015, 290, 17173–17180.
    51. Anderson, D.R. Normal-tension glaucoma (Low-tension glaucoma). Indian J. Ophthalmol. 2011, 59 (Suppl. 1), S97–S101.
    52. Weisman, R.L.; Asseff, C.F.; Phelps, C.D.; Podos, S.M.; Becker, B. Vertical elongation of the optic cup in glaucoma. Trans. Am. Acad. Ophthalmol. Otolaryngol. 1973, 77, OP157–OP161.
    53. Tatham, A.J.; Weinreb, R.N.; Zangwill, L.M.; Liebmann, J.M.; Girkin, C.A.; Medeiros, F.A. The relationship between cup-to-disc ratio and estimated number of retinal ganglion cells. Invest. Ophthalmol. Vis. Sci. 2013, 54, 3205–3214.
    54. Hopley, C.R.; Crowston, J.; Healey, P.R.; Lee, A.; Rochtchina, E.; Mitchell, P.; Study, B.M.E. Vertical Cup: Disc Ratio in a Non-glaucomatous Population-based Sample: The Blue Mountains Eye Study. Invest. Ophthalmol. Vis. Sci. 2003, 44, 3408.
    55. Hollands, H.; Johnson, D.; Hollands, S.; Simel, D.L.; Jinapriya, D.; Sharma, S. Do findings on routine examination identify patients at risk for primary open-angle glaucoma? The rational clinical examination systematic review. JAMA 2013, 309, 2035–2042.
    56. Nannini, D.R.; Torres, M.; Chen, Y.I.; Taylor, K.D.; Rotter, J.I.; Varma, R.; Gao, X. A Genome-Wide Association Study of Vertical Cup-Disc Ratio in a Latino Population. Invest. Ophthalmol. Vis. Sci. 2017, 58, 87–95.
    57. Ramdas, W.D.; van Koolwijk, L.M.; Ikram, M.K.; Jansonius, N.M.; de Jong, P.T.; Bergen, A.A.; Isaacs, A.; Amin, N.; Aulchenko, Y.S.; Wolfs, R.C.; et al. A genome-wide association study of optic disc parameters. PLoS Genet. 2010, 6, e1000978.
    58. Springelkamp, H.; Höhn, R.; Mishra, A.; Hysi, P.G.; Khor, C.-C.; Loomis, S.J.; Bailey, J.N.C.; Gibson, J.; Thorleifsson, G.; Janssen, S.F.; et al. Meta-analysis of genome-wide association studies identifies novel loci that influence cupping and the glaucomatous process. Nat. Commun. 2014, 5, 4883.
    59. Springelkamp, H.; Iglesias, A.I.; Mishra, A.; Höhn, R.; Wojciechowski, R.; Khawaja, A.P.; Nag, A.; Wang, Y.X.; Wang, J.J.; Cuellar-Partida, G.; et al. New insights into the genetics of primary open-angle glaucoma based on meta-analyses of intraocular pressure and optic disc characteristics. Hum. Mol. Genet. 2017, 26, 438–453.
    60. Nannini, D.R.; Kim, H.; Fan, F.; Gao, X. Genetic Risk Score Is Associated with Vertical Cup-to-Disc Ratio and Improves Prediction of Primary Open-Angle Glaucoma in Latinos. Ophthalmology 2018, 125, 815–821.
    61. Han, X.; Steven, K.; Qassim, A.; Marshall, H.N.; Bean, C.; Tremeer, M.; An, J.; Siggs, O.M.; Gharahkhani, P.; Craig, J.E.; et al. Automated AI labeling of optic nerve head enables insights into cross-ancestry glaucoma risk and genetic discovery in >280,000 images from UKB and CLSA. Am. J. Hum. Genet. 2021, 108, 1204–1216.
    62. Han, X.; Qassim, A.; An, J.; Marshall, H.; Zhou, T.; Ong, J.S.; Hassall, M.M.; Hysi, P.G.; Foster, P.J.; Khaw, P.T.; et al. Genome-wide association analysis of 95 549 individuals identifies novel loci and genes influencing optic disc morphology. Hum. Mol. Genet. 2019, 28, 3680–3690.
    63. Pers, T.H.; Karjalainen, J.M.; Chan, Y.; Westra, H.J.; Wood, A.R.; Yang, J.; Lui, J.C.; Vedantam, S.; Gustafsson, S.; Esko, T.; et al. Biological interpretation of genome-wide association studies using predicted gene functions. Nat. Commun. 2015, 6, 5890.
    64. Chen, Y.; Lin, Y.; Vithana, E.N.; Jia, L.; Zuo, X.; Wong, T.Y.; Chen, L.J.; Zhu, X.; Tam, P.O.S.; Gong, B.; et al. Common variants near ABCA1 and in PMM2 are associated with primary open-angle glaucoma. Nat. Genet. 2014, 46, 1115–1119.
    65. Alipanahi, B.; Hormozdiari, F.; Behsaz, B.; Cosentino, J.; McCaw, Z.R.; Schorsch, E.; Sculley, D.; Dorfman, E.H.; Foster, P.J.; Peng, L.H.; et al. Large-scale machine learning-based phenotyping significantly improves genomic discovery for optic nerve head morphology. Am. J. Hum. Genet. 2021, 108, 1217–1230.
    66. García-Bermúdez, M.Y.; Freude, K.K.; Mouhammad, Z.A.; van Wijngaarden, P.; Martin, K.K.; Kolko, M. Glial Cells in Glaucoma: Friends, Foes, and Potential Therapeutic Targets. Front. Neurol. 2021, 12, 624983.
    67. Mithieux, S.M.; Weiss, A.S. Elastin. Adv. Protein Chem. 2005, 70, 437–461.
    68. Melin, B.S.; Barnholtz-Sloan, J.S.; Wrensch, M.R.; Johansen, C.; Il’yasova, D.; Kinnersley, B.; Ostrom, Q.T.; Labreche, K.; Chen, Y.; Armstrong, G.; et al. Genome-wide association study of glioma subtypes identifies specific differences in genetic susceptibility to glioblastoma and non-glioblastoma tumors. Nat. Genet. 2017, 49, 789–794.
    69. Fan, B.J.; Wang, D.Y.; Pasquale, L.R.; Haines, J.L.; Wiggs, J.L. Genetic variants associated with optic nerve vertical cup-to-disc ratio are risk factors for primary open angle glaucoma in a US Caucasian population. Invest. Ophthalmol. Vis. Sci. 2011, 52, 1788–1792.
    70. Newman, E.; Reichenbach, A. The Müller cell: A functional element of the retina. Trends Neurosci. 1996, 19, 307–312.