Genome-wide association study identifies four pan-ancestry loci for suicidal ideation in the Million Veteran Program.

Suicidal ideation (SI) often precedes and predicts suicide attempt and death, is the most common suicidal phenotype and is over-represented in veterans. The genetic architecture of SI in the absence of suicide attempt (SA) is unknown, yet believed to have distinct and overlapping risk with other suicidal behaviors. We performed the first GWAS of SI without SA in the Million Veteran Program (MVP), identifying 99,814 SI cases from electronic health records without a history of SA or suicide death (SD) and 512,567 controls without SI, SA or SD. GWAS was performed separately in the four largest ancestry groups, controlling for sex, age and genetic substructure. Ancestry-specific results were combined via meta-analysis to identify pan-ancestry loci. Four genome-wide significant (GWS) loci were identified in the pan-ancestry meta-analysis with loci on chromosomes 6 and 9 associated with suicide attempt in an independent sample. Pan-ancestry gene-based analysis identified GWS associations with DRD2, DCC, FBXL19, BCL7C, CTF1, ANNK1, and EXD3. Gene-set analysis implicated synaptic and startle response pathways (q's<0.05). European ancestry (EA) analysis identified GWS loci on chromosomes 6 and 9, as well as GWS gene associations in EXD3, DRD2, and DCC. No other ancestry-specific GWS results were identified, underscoring the need to increase representation of diverse individuals. The genetic correlation of SI and SA within MVP was high (rG = 0.87; p = 1.09e-50), as well as with post-traumatic stress disorder (PTSD; rG = 0.78; p = 1.98e-95) and major depressive disorder (MDD; rG = 0.78; p = 8.33e-83). Conditional analysis on PTSD and MDD attenuated most pan-ancestry and EA GWS signals for SI without SA to nominal significance, with the exception of EXD3 which remained GWS. Our novel findings support a polygenic and complex architecture for SI without SA which is largely shared with SA and overlaps with psychiatric conditions frequently comorbid with suicidal behaviors.

Rare protective variants and glaucoma-relevant cell stressors modulate Angiopoietin-like 7 expression.

Rare missense and nonsense variants in the Angiopoietin-like 7 (ANGPTL7) gene confer protection from primary open-angle glaucoma (POAG), though the functional mechanism remains uncharacterized. Interestingly, a larger variant effect size strongly correlates with in silico predictions of increased protein instability (r = -0.98), suggesting that protective variants lower ANGPTL7 protein levels. Here, we show that missense and nonsense variants cause aggregation of mutant ANGPTL7 protein in the endoplasmic reticulum (ER) and decreased levels of secreted protein in human trabecular meshwork (TM) cells; a lower secreted:intracellular protein ratio strongly correlates with variant effects on intraocular pressure (r = 0.81). Importantly, accumulation of mutant protein in the ER does not increase expression of ER stress proteins in TM cells (P > 0.05 for all variants tested). Cyclic mechanical stress, a glaucoma-relevant physiologic stressor, also significantly lowers ANGPTL7 expression in primary cultures of human Schlemm's canal (SC) cells (-2.4-fold-change, P = 0.01). Collectively, these data suggest that the protective effects of ANGPTL7 variants in POAG stem from lower levels of secreted protein, which may modulate responses to physiologic and pathologic ocular cell stressors. Downregulation of ANGPTL7 expression may therefore serve as a viable preventative and therapeutic strategy for this common, blinding disease.

Lysyl oxidase-like 1-antisense 1 (LOXL1-AS1) lncRNA differentially regulates gene and protein expression, signaling and morphology of human ocular cells.

Pseudoexfoliation glaucoma (PEXG) is characterized by dysregulated extracellular matrix (ECM) homeostasis that disrupts conventional outflow function and increases intraocular pressure (IOP). Prolonged IOP elevation results in optic nerve head damage and vision loss. Uniquely, PEXG is a form of open angle glaucoma that has variable penetrance, is difficult to treat and does not respond well to common IOP-lowering pharmaceuticals. Therefore, understanding modulators of disease severity will aid in targeted therapies for PEXG. Genome-wide association studies have identified polymorphisms in the long non-coding RNA lysyl oxidase-like 1-antisense 1 (LOXL1-AS1) as a risk factor for PEXG. Risk alleles, oxidative stress and mechanical stretch all alter LOXL1-AS1 expression. As a long non-coding RNA, LOXL1-AS1 binds hnRNPL and regulates global gene expression. In this study, we focus on the role of LOXL1-AS1 in the ocular cells (trabecular meshwork and Schlemm's canal) that regulate IOP. We show that selective knockdown of LOXL1-AS1 leads to cell-type-specific changes in gene expression, ECM homeostasis, signaling and morphology. These results implicate LOXL1-AS1 as a modulator of cellular homeostasis, altering cell contractility and ECM turnover, both of which are well-known contributors to PEXG. These findings support LOXL1-AS1 as a key target for modifying the disease.

A genome-wide association study of suicide attempts in the million veterans program identifies evidence of pan-ancestry and ancestry-specific risk loci.

To identify pan-ancestry and ancestry-specific loci associated with attempting suicide among veterans, we conducted a genome-wide association study (GWAS) of suicide attempts within a large, multi-ancestry cohort of U.S. veterans enrolled in the Million Veterans Program (MVP). Cases were defined as veterans with a documented history of suicide attempts in the electronic health record (EHR; N = 14,089) and controls were defined as veterans with no documented history of suicidal thoughts or behaviors in the EHR (N = 395,064). GWAS was performed separately in each ancestry group, controlling for sex, age and genetic substructure. Pan-ancestry risk loci were identified through meta-analysis and included two genome-wide significant loci on chromosomes 20 (p = 3.64 × 10-9) and 1 (p = 3.69 × 10-8). A strong pan-ancestry signal at the Dopamine Receptor D2 locus (p = 1.77 × 10-7) was also identified and subsequently replicated in a large, independent international civilian cohort (p = 7.97 × 10-4). Additionally, ancestry-specific genome-wide significant loci were also detected in African-Americans, European-Americans, Asian-Americans, and Hispanic-Americans. Pathway analyses suggested over-representation of many biological pathways with high clinical significance, including oxytocin signaling, glutamatergic synapse, cortisol synthesis and secretion, dopaminergic synapse, and circadian rhythm. These findings confirm that the genetic architecture underlying suicide attempt risk is complex and includes both pan-ancestry and ancestry-specific risk loci. Moreover, pathway analyses suggested many commonly impacted biological pathways that could inform development of improved therapeutics for suicide prevention.

Identification and activity of the functional complex between hnRNPL and the pseudoexfoliation syndrome-associated lncRNA, LOXL1-AS1.

Individuals with pseudoexfoliation (PEX) syndrome exhibit various connective tissue pathologies associated with dysregulated extracellular matrix homeostasis. PEX glaucoma is a common, aggressive form of open-angle glaucoma resulting from the deposition of fibrillary material in the conventional outflow pathway. However, the molecular mechanisms that drive pathogenesis and genetic risk remain poorly understood. PEX glaucoma-associated single-nucleotide polymorphisms are located in and affect activity of the promoter of LOXL1-AS1, a long non-coding RNA (lncRNA). Nuclear and non-nuclear lncRNAs regulate a host of biological processes, and when dysregulated, contribute to disease. Here we report that LOXL1-AS1 localizes to the nucleus where it selectively binds to the mRNA processing protein, heterogeneous nuclear ribonucleoprotein-L (hnRNPL). Both components of this complex are critical for the regulation of global gene expression in ocular cells, making LOXL1-AS1 a prime target for investigation in PEX syndrome and glaucoma.

Integral role for lysyl oxidase-like-1 in conventional outflow tissue function and behavior.

Lysyl oxidase-like-1 (LOXL1), a vital crosslinking enzyme in elastin fiber maintenance, is essential for the stability and strength of elastic vessels and tissues. Variants in the LOXL1 locus associate with a dramatic increase in risk of exfoliation syndrome (XFS), a systemic fibrillopathy, which often presents with ocular hypertension and exfoliation glaucoma (XFG). We examined the role of LOXL1 in conventional outflow function, the prime regulator of intraocular pressure (IOP). Using Loxl1-/- , Loxl1+/- , and Loxl1+/+ mice, we observed an inverse relationship between LOXL1 expression and IOP, which worsened with age. Elevated IOP in Loxl1-/- mice was associated with a larger globe, decreased ocular compliance, increased outflow facility, extracellular matrix (ECM) abnormalities, and dilated intrascleral veins, yet, no dilation of arteries or capillaries. Interestingly, in living Loxl1-/- mouse eyes, Schlemm's canal (SC) was less susceptible to collapse when challenged with acute elevations in IOP, suggesting elevated episcleral venous pressure (EVP). Thus, LOXL1 expression is required for normal IOP control, while ablation results in altered ECM repair/homeostasis and conventional outflow physiology. Dilation of SC and distal veins, but not arteries, is consistent with key structural and functional roles for elastin in low-pressure vessels subjected to cyclical mechanical stress.

Genomic locus modulating corneal thickness in the mouse identifies POU6F2 as a potential risk of developing glaucoma.

Central corneal thickness (CCT) is one of the most heritable ocular traits and it is also a phenotypic risk factor for primary open angle glaucoma (POAG). The present study uses the BXD Recombinant Inbred (RI) strains to identify novel quantitative trait loci (QTLs) modulating CCT in the mouse with the potential of identifying a molecular link between CCT and risk of developing POAG. The BXD RI strain set was used to define mammalian genomic loci modulating CCT, with a total of 818 corneas measured from 61 BXD RI strains (between 60-100 days of age). The mice were anesthetized and the eyes were positioned in front of the lens of the Phoenix Micron IV Image-Guided OCT system or the Bioptigen OCT system. CCT data for each strain was averaged and used to QTLs modulating this phenotype using the bioinformatics tools on GeneNetwork (www.genenetwork.org). The candidate genes and genomic loci identified in the mouse were then directly compared with the summary data from a human POAG genome wide association study (NEIGHBORHOOD) to determine if any genomic elements modulating mouse CCT are also risk factors for POAG.This analysis revealed one significant QTL on Chr 13 and a suggestive QTL on Chr 7. The significant locus on Chr 13 (13 to 19 Mb) was examined further to define candidate genes modulating this eye phenotype. For the Chr 13 QTL in the mouse, only one gene in the region (Pou6f2) contained nonsynonymous SNPs. Of these five nonsynonymous SNPs in Pou6f2, two resulted in changes in the amino acid proline which could result in altered secondary structure affecting protein function. The 7 Mb region under the mouse Chr 13 peak distributes over 2 chromosomes in the human: Chr 1 and Chr 7. These genomic loci were examined in the NEIGHBORHOOD database to determine if they are potential risk factors for human glaucoma identified using meta-data from human GWAS. The top 50 hits all resided within one gene (POU6F2), with the highest significance level of p = 10-6 for SNP rs76319873. POU6F2 is found in retinal ganglion cells and in corneal limbal stem cells. To test the effect of POU6F2 on CCT we examined the corneas of a Pou6f2-null mice and the corneas were thinner than those of wild-type littermates. In addition, these POU6F2 RGCs die early in the DBA/2J model of glaucoma than most RGCs. Using a mouse genetic reference panel, we identified a transcription factor, Pou6f2, that modulates CCT in the mouse. POU6F2 is also found in a subset of retinal ganglion cells and these RGCs are sensitive to injury.

Identification of Estrogen Signaling in a Prioritization Study of Intraocular Pressure-Associated Genes.

Elevated intraocular pressure (IOP) is the only modifiable risk factor for primary open-angle glaucoma (POAG). Herein we sought to prioritize a set of previously identified IOP-associated genes using novel and previously published datasets. We identified several genes for future study, including several involved in cytoskeletal/extracellular matrix reorganization, cell adhesion, angiogenesis, and TGF-ß signaling. Our differential correlation analysis of IOP-associated genes identified 295 pairs of 201 genes with differential correlation. Pathway analysis identified ß-estradiol as the top upstream regulator of these genes with ESR1 mediating 25 interactions. Several genes (i.e., EFEMP1, FOXC1, and SPTBN1) regulated by ß-estradiol/ESR1 were highly expressed in non-glaucomatous human trabecular meshwork (TM) or Schlemm's canal (SC) cells and specifically expressed in TM/SC cell clusters defined by single-cell RNA-sequencing. We confirmed ESR1 gene and protein expression in human TM cells and TM/SC tissue with quantitative real-time PCR and immunofluorescence, respectively. 17ß-estradiol was identified in bovine, porcine, and human aqueous humor (AH) using ELISA. In conclusion, we have identified estrogen receptor signaling as a key modulator of several IOP-associated genes. The expression of ESR1 and these IOP-associated genes in TM/SC tissue and the presence of 17ß-estradiol in AH supports a role for estrogen signaling in IOP regulation.

Differential DNA methylation patterns in human Schlemm's canal endothelial cells with glaucoma.

Purpose: Schlemm's canal (SC) endothelial cells derived from donors with or without glaucoma showed different mechanical properties and gene expression. As an important contributor to the regulation of intraocular pressure (IOP) and pathogenesis of primary open-angle glaucoma (POAG), the heritable key epigenetic changes, methylation may play an important role in the physiologic function of SC cells. This study aims to identify differentially methylated CpG sites (DMSs) in primary cultures of human SC cells with or without glaucoma. Methods: We examined the methylation pattern of seven strains of primary human cells (two glaucoma and five normal SC cell samples), which were isolated and characterized using established protocols. DNA methylation was profiled using Illumina Human Methylation 450 BeadChip. Raw data were extracted and exported using Illumina GenomeStudio software. After quantile normalization, DNA methylation data were analyzed using R package RnBeads in Bioconductor. DMSs were filtered with p ≤ 1E-5, methylation change ≥ 0.1, and false discovery rate ≤ 0.05. The closest genes and the location of each CpG site were annotated using R package FDb.InfiniumMethylation.hg19. Gene Ontology and pathway analysis was performed using WebGestalt. Selected DMSs were validated using the Zymo qMethyl kit. Results: We used five non-glaucoma and two glaucomatous SC cell samples to profile genome-wide DNA methylation using Illumina Infinium Methylation BeadChips. Principle component analysis showed the separation between the glaucoma and control samples. After quality control and differential analysis, we identified 298 highly significant DMSs (p ≤ 1E-5). Among them, 221 DMSs were within 1 kb of a nearby gene. Gene Ontology analysis demonstrated significant enrichment in positive regulation of cell migration, negative regulation of endothelial cell proliferation, and stress fiber and actin filament bundles. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed enrichment in cell adhesion and gap junctions. Several glaucoma-related genes were identified, including TGFBR3, THBS1, PITX2, DAXX, TBX3, TNXB, ANGPT1, and PLEKHA7. We also examined differentially methylated regions (DMRs) near these CpG sites and identified significant DMRs in TBX3, TNXB1, DAXX, and PITX2. Conclusions: This study represents the first genome-wide DNA methylation profiling in cultured human primary SC cells. The DMSs were enriched in the pathways related to outflow resistance. Several DMRs were validated in glaucoma-associated genes, further suggesting the role of DNA methylation in glaucoma development. This study could provide comprehensive understanding of DNA methylation in glaucoma and its effect on aqueous humor outflow.

Update on the genetics of primary open-angle glaucoma.

Affecting nearly 80 million individuals, glaucoma is the number one cause of irreversible blindness in the world. This ocular disease describes a set of optic neuropathies of which primary open angle glaucoma (POAG) is the most common. POAG is associated with progressive visual field deterioration resulting from damage to the optic nerve and loss of retinal ganglion cells. Risk factors for POAG include elevated intraocular pressure, aging, African and Hispanic ancestry, and a positive family history of POAG. Multiple genes have been found to contribute to POAG. Much of POAG genetics and pathology has yet to be explained. Recent genome-wide association studies have identified a large number of novel loci associated with POAG and its endophenotypes. Genomic and proteomic profiling of biofluids has contributed to our knowledge of differential gene expression in POAG. Functional studies both in cell culture and animal models have confirmed the effects of variants and differential gene expression on ocular physiology while in silico analyses have increased our understanding of disease risk and progression so that we might better diagnose and treat this complex genetic illness.

Epistatic Gene-Based Interaction Analyses for Glaucoma in eMERGE and NEIGHBOR Consortium.

Primary open angle glaucoma (POAG) is a complex disease and is one of the major leading causes of blindness worldwide. Genome-wide association studies have successfully identified several common variants associated with glaucoma; however, most of these variants only explain a small proportion of the genetic risk. Apart from the standard approach to identify main effects of variants across the genome, it is believed that gene-gene interactions can help elucidate part of the missing heritability by allowing for the test of interactions between genetic variants to mimic the complex nature of biology. To explain the etiology of glaucoma, we first performed a genome-wide association study (GWAS) on glaucoma case-control samples obtained from electronic medical records (EMR) to establish the utility of EMR data in detecting non-spurious and relevant associations; this analysis was aimed at confirming already known associations with glaucoma and validating the EMR derived glaucoma phenotype. Our findings from GWAS suggest consistent evidence of several known associations in POAG. We then performed an interaction analysis for variants found to be marginally associated with glaucoma (SNPs with main effect p-value <0.01) and observed interesting findings in the electronic MEdical Records and GEnomics Network (eMERGE) network dataset. Genes from the top epistatic interactions from eMERGE data (Likelihood Ratio Test i.e. LRT p-value <1e-05) were then tested for replication in the NEIGHBOR consortium dataset. To replicate our findings, we performed a gene-based SNP-SNP interaction analysis in NEIGHBOR and observed significant gene-gene interactions (p-value <0.001) among the top 17 gene-gene models identified in the discovery phase. Variants from gene-gene interaction analysis that we found to be associated with POAG explain 3.5% of additional genetic variance in eMERGE dataset above what is explained by the SNPs in genes that are replicated from previous GWAS studies (which was only 2.1% variance explained in eMERGE dataset); in the NEIGHBOR dataset, adding replicated SNPs from gene-gene interaction analysis explain 3.4% of total variance whereas GWAS SNPs alone explain only 2.8% of variance. Exploring gene-gene interactions may provide additional insights into many complex traits when explored in properly designed and powered association studies.

Association of Genetic Variants With Primary Open-Angle Glaucoma Among Individuals With African Ancestry.

Importance: Primary open-angle glaucoma presents with increased prevalence and a higher degree of clinical severity in populations of African ancestry compared with European or Asian ancestry. Despite this, individuals of African ancestry remain understudied in genomic research for blinding disorders. Objectives: To perform a genome-wide association study (GWAS) of African ancestry populations and evaluate potential mechanisms of pathogenesis for loci associated with primary open-angle glaucoma. Design, Settings, and Participants: A 2-stage GWAS with a discovery data set of 2320 individuals with primary open-angle glaucoma and 2121 control individuals without primary open-angle glaucoma. The validation stage included an additional 6937 affected individuals and 14 917 unaffected individuals using multicenter clinic- and population-based participant recruitment approaches. Study participants were recruited from Ghana, Nigeria, South Africa, the United States, Tanzania, Britain, Cameroon, Saudi Arabia, Brazil, the Democratic Republic of the Congo, Morocco, Peru, and Mali from 2003 to 2018. Individuals with primary open-angle glaucoma had open iridocorneal angles and displayed glaucomatous optic neuropathy with visual field defects. Elevated intraocular pressure was not included in the case definition. Control individuals had no elevated intraocular pressure and no signs of glaucoma. Exposures: Genetic variants associated with primary open-angle glaucoma. Main Outcomes and Measures: Presence of primary open-angle glaucoma. Genome-wide significance was defined as P < 5 × 10-8 in the discovery stage and in the meta-analysis of combined discovery and validation data. Results: A total of 2320 individuals with primary open-angle glaucoma (mean [interquartile range] age, 64.6 [56-74] years; 1055 [45.5%] women) and 2121 individuals without primary open-angle glaucoma (mean [interquartile range] age, 63.4 [55-71] years; 1025 [48.3%] women) were included in the discovery GWAS. The GWAS discovery meta-analysis demonstrated association of variants at amyloid-ß A4 precursor protein-binding family B member 2 (APBB2; chromosome 4, rs59892895T>C) with primary open-angle glaucoma (odds ratio [OR], 1.32 [95% CI, 1.20-1.46]; P = 2 × 10-8). The association was validated in an analysis of an additional 6937 affected individuals and 14 917 unaffected individuals (OR, 1.15 [95% CI, 1.09-1.21]; P < .001). Each copy of the rs59892895*C risk allele was associated with increased risk of primary open-angle glaucoma when all data were included in a meta-analysis (OR, 1.19 [95% CI, 1.14-1.25]; P = 4 × 10-13). The rs59892895*C risk allele was present at appreciable frequency only in African ancestry populations. In contrast, the rs59892895*C risk allele had a frequency of less than 0.1% in individuals of European or Asian ancestry. Conclusions and Relevance: In this genome-wide association study, variants at the APBB2 locus demonstrated differential association with primary open-angle glaucoma by ancestry. If validated in additional populations this finding may have implications for risk assessment and therapeutic strategies.

Genome-wide association study of primary open-angle glaucoma in continental and admixed African populations.

Primary open angle glaucoma (POAG) is a complex disease with a major genetic contribution. Its prevalence varies greatly among ethnic groups, and is up to five times more frequent in black African populations compared to Europeans. So far, worldwide efforts to elucidate the genetic complexity of POAG in African populations has been limited. We conducted a genome-wide association study in 1113 POAG cases and 1826 controls from Tanzanian, South African and African American study samples. Apart from confirming evidence of association at TXNRD2 (rs16984299; OR[T] 1.20; P = 0.003), we found that a genetic risk score combining the effects of the 15 previously reported POAG loci was significantly associated with POAG in our samples (OR 1.56; 95% CI 1.26-1.93; P = 4.79 × 10-5). By genome-wide association testing we identified a novel candidate locus, rs141186647, harboring EXOC4 (OR[A] 0.48; P = 3.75 × 10-8), a gene transcribing a component of the exocyst complex involved in vesicle transport. The low frequency and high degree of genetic heterogeneity at this region hampered validation of this finding in predominantly West-African replication sets. Our results suggest that established genetic risk factors play a role in African POAG, however, they do not explain the higher disease load. The high heterogeneity within Africans remains a challenge to identify the genetic commonalities for POAG in this ethnicity, and demands studies of extremely large size.

Modeling Glaucoma: Retinal Ganglion Cells Generated from Induced Pluripotent Stem Cells of Patients with SIX6 Risk Allele Show Developmental Abnormalities.

Glaucoma represents a group of multifactorial diseases with a unifying pathology of progressive retinal ganglion cell (RGC) degeneration, causing irreversible vision loss. To test the hypothesis that RGCs are intrinsically vulnerable in glaucoma, we have developed an in vitro model using the SIX6 risk allele carrying glaucoma patient-specific induced pluripotent stem cells (iPSCs) for generating functional RGCs. Here, we demonstrate that the efficiency of RGC generation by SIX6 risk allele iPSCs is significantly lower than iPSCs-derived from healthy, age- and sex-matched controls. The decrease in the number of RGC generation is accompanied by repressed developmental expression of RGC regulatory genes. The SIX6 risk allele RGCs display short and simple neurites, reduced expression of guidance molecules, and immature electrophysiological signature. In addition, these cells have higher expression of glaucoma-associated genes, CDKN2A and CDKN2B, suggesting an early onset of the disease phenotype. Consistent with the developmental abnormalities, the SIX6 risk allele RGCs display global dysregulation of genes which map on developmentally relevant biological processes for RGC differentiation and signaling pathways such as mammalian target of rapamycin that integrate diverse functions for differentiation, metabolism, and survival. The results suggest that SIX6 influences different stages of RGC differentiation and their survival; therefore, alteration in SIX6 function due to the risk allele may lead to cellular and molecular abnormalities. These abnormalities, if carried into adulthood, may make RGCs vulnerable in glaucoma. Stem Cells 2017;35:2239-2252.

Transcriptome analysis of adult and fetal trabecular meshwork, cornea, and ciliary body tissues by RNA sequencing.

PURPOSE: To characterize the transcriptional landscape of human adult and fetal trabecular meshwork (TM), cornea, and ciliary body (CB) tissues, and to evaluate the expression level of candidate genes selected from genetic association studies of primary-open angle glaucoma, central corneal thickness, intraocular pressure, vertical cup to disc ratio, and optic nerve parameters. METHODS: Deep RNA sequencing was performed on the selected human tissues. Transcriptome analyses were performed to 1) characterize the total number of expressed genes, 2) identify the most highly expressed genes, 3) estimate the number of novel transcripts, and 4) evaluate the expression of candidate genes in each tissue. Finally, a differential gene expression analysis was conducted to compare the adult and fetal ocular tissues. RESULTS: There was an average of 12,362 protein coding genes and 3725 novel transcripts expressed in each tissue. The top most expressed genes in each tissue included SPARC (fetal cornea and TM), APOD (adult TM), CLU (adult cornea), and PTGDS (adult and fetal CB). Twenty-nine candidate genes selected from genetic association studies primarily showed high expression levels in the trabecular meshwork and cornea. Comparison of adult and fetal samples identified 2012 and 1261 differentially expressed protein-coding genes within the cornea and trabecular meshwork, respectively. CONCLUSIONS: This study has provided an unbiased glimpse into the transcriptome of three essential anterior ocular tissues, resulting in the development of several novel hypotheses. These data can be used in the future to better guide ocular research questions.

Genetic variants and cellular stressors associated with exfoliation syndrome modulate promoter activity of a lncRNA within the LOXL1 locus.

Exfoliation syndrome (XFS) is a common, age-related, systemic fibrillinopathy. It greatly increases risk of exfoliation glaucoma (XFG), a major worldwide cause of irreversible blindness. Coding variants in the lysyl oxidase-like 1 (LOXL1) gene are strongly associated with XFS in all studied populations, but a functional role for these variants has not been established. To identify additional candidate functional variants, we sequenced the entire LOXL1 genomic locus (∼40 kb) in 50 indigenous, black South African XFS cases and 50 matched controls. The variants with the strongest evidence of association were located in a well-defined 7-kb region bounded by the 3'-end of exon 1 and the adjacent region of intron 1 of LOXL1. We replicated this finding in US Caucasian (91 cases/1031 controls), German (771 cases/1365 controls) and Japanese (1484 cases/1188 controls) populations. The region of peak association lies upstream of LOXL1-AS1, a long non-coding RNA (lncRNA) encoded on the opposite strand of LOXL1. We show that this region contains a promoter and, importantly, that the strongly associated XFS risk alleles in the South African population are functional variants that significantly modulate the activity of this promoter. LOXL1-AS1 expression is also significantly altered in response to oxidative stress in human lens epithelial cells and in response to cyclic mechanical stress in human Schlemm's canal endothelial cells. Taken together, these findings support a functional role for the LOXL1-AS1 lncRNA in cellular stress response and suggest that dysregulation of its expression by genetic risk variants plays a key role in XFS pathogenesis.

Major review: Exfoliation syndrome; advances in disease genetics, molecular biology, and epidemiology.

Exfoliation syndrome (XFS) is a common age-related disorder that leads to deposition of extracellular fibrillar material throughout the body. The most recognized disease manifestation is exfoliation glaucoma (XFG), which is a common cause of blindness worldwide. Recent developments in XFS genetics, cell biology and epidemiology have greatly improved our understanding of the etiology of this complex inherited disease. This review summarizes current knowledge of XFS pathogenesis, identifies gaps in knowledge, and discusses areas for future research.

Discovery and functional annotation of SIX6 variants in primary open-angle glaucoma.

Glaucoma is a leading cause of blindness worldwide. Primary open-angle glaucoma (POAG) is the most common subtype and is a complex trait with multigenic inheritance. Genome-wide association studies have previously identified a significant association between POAG and the SIX6 locus (rs10483727, odds ratio (OR) = 1.32, p = 3.87×10(-11)). SIX6 plays a role in ocular development and has been associated with the morphology of the optic nerve. We sequenced the SIX6 coding and regulatory regions in 262 POAG cases and 256 controls and identified six nonsynonymous coding variants, including five rare and one common variant, Asn141His (rs33912345), which was associated significantly with POAG (OR = 1.27, p = 4.2×10(-10)) in the NEIGHBOR/GLAUGEN datasets. These variants were tested in an in vivo Danio rerio (zebrafish) complementation assay to evaluate ocular metrics such as eye size and optic nerve structure. Five variants, found primarily in POAG cases, were hypomorphic or null, while the sixth variant, found only in controls, was benign. One variant in the SIX6 enhancer increased expression of SIX6 and disrupted its regulation. Finally, to our knowledge for the first time, we have identified a clinical feature in POAG patients that appears to be dependent upon SIX6 genotype: patients who are homozygous for the SIX6 risk allele (His141) have a statistically thinner retinal nerve fiber layer than patients homozygous for the SIX6 non-risk allele (Asn141). Our results, in combination with previous SIX6 work, lead us to hypothesize that SIX6 risk variants disrupt the development of the neural retina, leading to a reduced number of retinal ganglion cells, thereby increasing the risk of glaucoma-associated vision loss.

ABCC5, a gene that influences the anterior chamber depth, is associated with primary angle closure glaucoma.

Anterior chamber depth (ACD) is a key anatomical risk factor for primary angle closure glaucoma (PACG). We conducted a genome-wide association study (GWAS) on ACD to discover novel genes for PACG on a total of 5,308 population-based individuals of Asian descent. Genome-wide significant association was observed at a sequence variant within ABCC5 (rs1401999; per-allele effect size =  -0.045 mm, P = 8.17 × 10(-9)). This locus was associated with an increase in risk of PACG in a separate case-control study of 4,276 PACG cases and 18,801 controls (per-allele OR = 1.13 [95% CI: 1.06-1.22], P = 0.00046). The association was strengthened when a sub-group of controls with open angles were included in the analysis (per-allele OR = 1.30, P = 7.45 × 10(-9); 3,458 cases vs. 3,831 controls). Our findings suggest that the increase in PACG risk could in part be mediated by genetic sequence variants influencing anterior chamber dimensions.

Epigenome-wide association of PTSD from heterogeneous cohorts with a common multi-site analysis pipeline.

Compelling evidence suggests that epigenetic mechanisms such as DNA methylation play a role in stress regulation and in the etiologic basis of stress related disorders such as Post traumatic Stress Disorder (PTSD). Here we describe the purpose and methods of an international consortium that was developed to study the role of epigenetics in PTSD. Inspired by the approach used in the Psychiatric Genomics Consortium, we brought together investigators representing seven cohorts with a collective sample size of N = 1147 that included detailed information on trauma exposure, PTSD symptoms, and genome-wide DNA methylation data. The objective of this consortium is to increase the analytical sample size by pooling data and combining expertise so that DNA methylation patterns associated with PTSD can be identified. Several quality control and analytical pipelines were evaluated for their control of genomic inflation and technical artifacts with a joint analysis procedure established to derive comparable data over the cohorts for meta-analysis. We propose methods to deal with ancestry population stratification and type I error inflation and discuss the advantages and disadvantages of applying robust error estimates. To evaluate our pipeline, we report results from an epigenome-wide association study (EWAS) of age, which is a well-characterized phenotype with known epigenetic associations. Overall, while EWAS are highly complex and subject to similar challenges as genome-wide association studies (GWAS), we demonstrate that an epigenetic meta-analysis with a relatively modest sample size can be well-powered to identify epigenetic associations. Our pipeline can be used as a framework for consortium efforts for EWAS.