Mitochondrial TXNRD2 and ME3 Genetic Risk Scores Are Associated with Specific Primary Open-Angle Glaucoma Phenotypes.

PURPOSE: Genetic variants in regions that include the mitochondrial genes thioredoxin reductase 2 (TXNRD2) and malic enzyme 3 (ME3) are associated with primary open-angle glaucoma (POAG) in genome-wide association studies (GWASs). To assess their clinical impact, we investigated whether TXNRD2 and ME3 genetic risk scores (GRSs) are associated with specific glaucoma phenotypes. DESIGN: Cross-sectional study. PARTICIPANTS: A total of 2617 patients with POAG and 2634 control participants from the National Eye Institute Glaucoma Human Genetics Collaboration Hereditable Overall Operational Database (NEIGHBORHOOD) consortium. METHODS: All POAG-associated single nucleotide polymorphisms (SNPs) in the TXNRD2 and ME3 loci were identified using GWAS data (P < 0.05). Of these, 20 TXNRD2 and 24 ME3 SNPs were selected after adjusting for linkage disequilibrium. The correlation between SNP effect size and gene expression levels was investigated using the Gene-Tissue Expression database. Genetic risk scores were constructed for each individual using the unweighted sum of TXNRD2, ME3, and TXNRD2 + ME3 combined risk alleles. Age- and sex-adjusted odds ratios (ORs) for POAG diagnosis were calculated per decile for each GRS. Additionally, the clinical features of patients with POAG in the top 1%, 5%, and 10% of each GRS were compared with those in the bottom 1%, 5%, and 10%, respectively. MAIN OUTCOME MEASURES: Primary open-angle glaucoma OR per GRS decile, maximum treated intraocular pressure (IOP), and prevalence of paracentral visual field loss among patients with POAG with high versus low GRSs. RESULTS: A larger SNP effect size strongly correlated with higher TXNRD2 and lower ME3 expression levels (r = 0.95 and r = -0.97, respectively; P < 0.05 for both). Individuals in decile 10 of the TXNRD2 + ME3 GRS had the highest odds of POAG diagnosis (OR, 1.79 compared with decile 1; 95% confidence interval, 1.39-2.30; P < 0.001). Patients with POAG in the top 1% of the TXNRD2 GRS showed higher mean maximum treated IOP compared with the bottom 1% (19.9 mmHg vs. 15.6 mmHg; adjusted P = 0.03). Patients with POAG in the top 1% of the ME3 and TXNRD2 + ME3 GRS showed a higher prevalence of paracentral field loss than the bottom 1% (72.7% vs. 14.3% for ME3 GRS and 88.9% vs. 33.3% for TXNRD2+ME3 GRS; adjusted P = 0.03 for both). CONCLUSIONS: Patients with POAG with higher TXNRD2 and ME3 GRSs showed higher treated IOP and a greater prevalence of paracentral field loss. Functional studies exploring how these variants impact mitochondrial function in patients with glaucoma are warranted. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found after the references.

T-cell receptor diversity in minimal change disease in the NEPTUNE study.

BACKGROUND: Minimal change disease (MCD) is the major cause of childhood idiopathic nephrotic syndrome, which is characterized by massive proteinuria and debilitating edema. Proteinuria in MCD is typically rapidly reversible with corticosteroid therapy, but relapses are common, and children often have many adverse events from the repeated courses of immunosuppressive therapy. The pathobiology of MCD remains poorly understood. Prior clinical observations suggest that abnormal T-cell function may play a central role in MCD pathogenesis. Based on these observations, we hypothesized that T-cell responses to specific exposures or antigens lead to a clonal expansion of T-cell subsets, a restriction in the T-cell repertoire, and an elaboration of specific circulating factors that trigger disease onset and relapses. METHODS: To test these hypotheses, we sequenced T-cell receptors in fourteen MCD, four focal segmental glomerulosclerosis (FSGS), and four membranous nephropathy (MN) patients with clinical data and blood samples drawn during active disease and during remission collected by the Nephrotic Syndrome Study Network (NEPTUNE). We calculated several T-cell receptor diversity metrics to assess possible differences between active disease and remission states in paired samples. RESULTS: Median productive clonality did not differ between MCD active disease (0.0083; range: 0.0042, 0.0397) and remission (0.0088; range: 0.0038, 0.0369). We did not identify dominant clonotypes in MCD active disease, and few clonotypes were shared with FSGS and MN patients. CONCLUSIONS: While these data do not support an obvious role of the adaptive immune system T-cells in MCD pathogenesis, further study is warranted given the limited sample size. A higher resolution version of the Graphical abstract is available as Supplementary information.

Multimodal MRI examination of structural and functional brain changes in older women with breast cancer in the first year of antiestrogen hormonal therapy.

PURPOSE: Cancer patients are concerned about treatment-related cognitive problems. We examined effects of antiestrogen hormonal therapy on brain imaging metrics in older women with breast cancer. METHODS: Women aged 60 + treated with hormonal therapy only and matched non-cancer controls (n = 29/group) completed MRI and objective and self-reported cognitive assessment at pre-treatment/enrollment and 12 months later. Gray matter was examined using voxel-based morphometry (VBM), FreeSurfer, and brain age calculations. Functional MRI (fMRI) assessed working memory-related activation. Analyses examined cross-sectional and longitudinal differences and tested associations between brain metrics, cognition, and days on hormonal therapy. RESULTS: The cancer group showed regional reductions over 12 months in frontal, temporal, and parietal gray matter on VBM, reduced FreeSurfer cortical thickness in prefrontal, parietal, and insular regions, and increased working memory-related fMRI activation in frontal, cingulate, and visual association cortex. Controls showed only reductions in fusiform gyrus on VBM and FreeSurfer temporal and parietal cortex thickness. Women with breast cancer showed higher estimated brain age and lower regional gray matter volume than controls at both time points. The cancer group showed a trend toward lower performance in attention, processing speed, and executive function at follow-up. There were no significant associations between brain imaging metrics and cognition or days on hormonal therapy. CONCLUSION: Older women with breast cancer showed brain changes in the first year of hormonal therapy. Increased brain activation during working memory processing may be a sign of functional compensation for treatment-related structural changes. This hypothesis should be tested in larger samples over longer time periods. GOV IDENTIFIER: NCT03451383.

Glaucoma Genetic Risk Scores in the Million Veteran Program.

PURPOSE: Primary open-angle glaucoma (POAG) is a degenerative eye disease for which early treatment is critical to mitigate visual impairment and irreversible blindness. POAG-associated loci individually confer incremental risk. Genetic risk score(s) (GRS) could enable POAG risk stratification. Despite significantly higher POAG burden among individuals of African ancestry (AFR), GRS are limited in this population. A recent large-scale, multi-ancestry meta-analysis identified 127 POAG-associated loci and calculated cross-ancestry and ancestry-specific effect estimates, including in European ancestry (EUR) and AFR individuals. We assessed the utility of the 127-variant GRS for POAG risk stratification in EUR and AFR Veterans in the Million Veteran Program (MVP). We also explored the association between GRS and documented invasive glaucoma surgery (IGS). DESIGN: Cross-sectional study. PARTICIPANTS: MVP Veterans with imputed genetic data, including 5830 POAG cases (445 with IGS documented in the electronic health record) and 64 476 controls. METHODS: We tested unweighted and weighted GRS of 127 published risk variants in EUR (3382 cases and 58 811 controls) and AFR (2448 cases and 5665 controls) Veterans in the MVP. Weighted GRS were calculated using effect estimates from the most recently published report of cross-ancestry and ancestry-specific meta-analyses. We also evaluated GRS in POAG cases with documented IGS. MAIN OUTCOME MEASURES: Performance of 127-variant GRS in EUR and AFR Veterans for POAG risk stratification and association with documented IGS. RESULTS: GRS were significantly associated with POAG (P < 5 × 10-5) in both groups; a higher proportion of EUR compared with AFR were consistently categorized in the top GRS decile (21.9%-23.6% and 12.9%-14.5%, respectively). Only GRS weighted by ancestry-specific effect estimates were associated with IGS documentation in AFR cases; all GRS types were associated with IGS in EUR cases. CONCLUSIONS: Varied performance of the GRS for POAG risk stratification and documented IGS association in EUR and AFR Veterans highlights (1) the complex risk architecture of POAG, (2) the importance of diverse representation in genomics studies that inform GRS construction and evaluation, and (3) the necessity of expanding diverse POAG-related genomic data so that GRS can equitably aid in screening individuals at high risk of POAG and who may require more aggressive treatment.

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.

Genome-wide association study identifies seven novel susceptibility loci for primary open-angle glaucoma.

Primary open-angle glaucoma (POAG) is the leading cause of irreversible blindness worldwide for which 15 disease-associated loci had been discovered. Among them, only 5 loci have been associated with POAG in Asians. We carried out a genome-wide association study and a replication study that included a total of 7378 POAG cases and 36 385 controls from a Japanese population. After combining the genome-wide association study and the two replication sets, we identified 11 POAG-associated loci, including 4 known (CDKN2B-AS1, ABCA1, SIX6 and AFAP1) and 7 novel loci (FNDC3B, ANKRD55-MAP3K1, LMX1B, LHPP, HMGA2, MEIS2 and LOXL1) at a genome-wide significance level (P < 5.0×10-8), bringing the total number of POAG-susceptibility loci to 22. The 7 novel variants were subsequently evaluated in a multiethnic population comprising non-Japanese East Asians (1008 cases, 591 controls), Europeans (5008 cases, 35 472 controls) and Africans (2341 cases, 2037 controls). The candidate genes located within the new loci were related to ocular development (LMX1B, HMGA2 and MAP3K1) and glaucoma-related phenotypes (FNDC3B, LMX1B and LOXL1). Pathway analysis suggested epidermal growth factor receptor signaling might be involved in POAG pathogenesis. Genetic correlation analysis revealed the relationships between POAG and systemic diseases, including type 2 diabetes and cardiovascular diseases. These results improve our understanding of the genetic factors that affect the risk of developing POAG and provide new insight into the genetic architecture of POAG in Asians.

New insights into the genetics of primary open-angle glaucoma based on meta-analyses of intraocular pressure and optic disc characteristics.

Primary open-angle glaucoma (POAG), the most common optic neuropathy, is a heritable disease. Siblings of POAG cases have a ten-fold increased risk of developing the disease. Intraocular pressure (IOP) and optic nerve head characteristics are used clinically to predict POAG risk. We conducted a genome-wide association meta-analysis of IOP and optic disc parameters and validated our findings in multiple sets of POAG cases and controls. Using imputation to the 1000 genomes (1000G) reference set, we identified 9 new genomic regions associated with vertical cup-disc ratio (VCDR) and 1 new region associated with IOP. Additionally, we found 5 novel loci for optic nerve cup area and 6 for disc area. Previously it was assumed that genetic variation influenced POAG either through IOP or via changes to the optic nerve head; here we present evidence that some genomic regions affect both IOP and the disc parameters. We characterized the effect of the novel loci through pathway analysis and found that pathways involved are not entirely distinct as assumed so far. Further, we identified a novel association between CDKN1A and POAG. Using a zebrafish model we show that six6b (associated with POAG and optic nerve head variation) alters the expression of cdkn1a. In summary, we have identified several novel genes influencing the major clinical risk predictors of POAG and showed that genetic variation in CDKN1A is important in POAG risk.

Mind the gap: resources required to receive, process and interpret research-returned whole genome data.

Most genotype-phenotype studies have historically lacked population diversity, impacting the generalizability of findings and thereby limiting the ability to equitably implement precision medicine. This well-documented problem has generated much interest in the ascertainment of new cohorts with an emphasis on multiple dimensions of diversity, including race/ethnicity, gender, age, socioeconomic status, disability, and geography. The most well known of these new cohort efforts is arguably All of Us, formerly known as the Precision Medicine Cohort Initiative Program. All of Us intends to ascertain at least one million participants in the United States representative of the multiple dimensions of diversity. As an incentive to participate, All of Us is offering the return of research results, including whole genome sequencing data, as well as the opportunity to contribute to the scientific process as non-scientists. The scale and scope of the proposed return of research results are unprecedented. Here, we briefly review possible return of genetic data models, including the likely data file formats and modes of data transfer or access. We also review the resources required to access and interpret the genetic or genomic data once received by the average participant, highlighting the nuanced anticipated barriers that will challenge both the digitally, computationally literate and illiterate participant alike. This inventory of resources required to receive, process, and interpret return of research results exposes the potential for access disparities and warns the scientific community to mind the gap so that all participants have equal access and understanding of the benefits of human genetic research.

Rare variants and loci for age-related macular degeneration in the Ohio and Indiana Amish.

Age-related macular degeneration (AMD) is a leading cause of blindness in the world. While dozens of independent genomic variants are associated with AMD, about one-third of AMD heritability is still unexplained. To identify novel variants and loci for AMD, we analyzed Illumina HumanExome chip data from 87 Amish individuals with early or late AMD, 79 unaffected Amish individuals, and 15 related Amish individuals with unknown AMD affection status. We retained 37,428 polymorphic autosomal variants across 175 samples for association and linkage analyses. After correcting for multiple testing (n = 37,428), we identified four variants significantly associated with AMD: rs200437673 (LCN9, p = 1.50 × 10-11), rs151214675 (RTEL1, p = 3.18 × 10-8), rs140250387 (DLGAP1, p = 4.49 × 10-7), and rs115333865 (CGRRF1, p = 1.05 × 10-6). These variants have not been previously associated with AMD and are not in linkage disequilibrium with the 52 known AMD-associated variants reported by the International AMD Genomics Consortium based on physical distance. Genome-wide significant linkage peaks were observed on chromosomes 8q21.11-q21.13 (maximum recessive HLOD = 4.03) and 18q21.2-21.32 (maximum dominant HLOD = 3.87; maximum recessive HLOD = 4.27). These loci do not overlap with loci previously linked to AMD. Through gene ontology enrichment analysis with ClueGO in Cytoscape, we determined that several genes in the 1-HLOD support interval of the chromosome 8 locus are involved in fatty acid binding and triglyceride catabolic processes, and the 1-HLOD support interval of the linkage region on chromosome 18 is enriched in genes that participate in serine-type endopeptidase inhibitor activity and the positive regulation of epithelial to mesenchymal transition. These results nominate novel variants and loci for AMD that require further investigation.

AMISH EYE STUDY: Baseline Spectral Domain Optical Coherence Tomography Characteristics of Age-Related Macular Degeneration.

PURPOSE: To describe spectral domain optical coherence tomography (SD-OCT) findings in an Amish cohort to assess SD-OCT markers for early age-related macular degeneration (AMD). METHODS: The authors performed a family-based prospective cohort study of 1,146 elderly Amish subjects (age range 50-99 years) (2,292 eyes) who had a family history of at least 1 individual with AMD. All subjects underwent complete ophthalmic examinations, SD-OCT using both Cirrus and Spectralis (20 × 20° scan area) instruments, fundus autofluorescence, infrared imaging, and color fundus photography. Spectral domain optical coherence tomography characteristics were analyzed in subjects with AMD (with and without subretinal drusenoid deposits [SDDs]) and normal healthy cohorts. RESULTS: Participants' mean age was 65.2 years (SD ± 11). Color fundus photographic findings in 596 (53%) subjects (1,009 eyes) were consistent with AMD; the remaining 478 (43%) subjects showed no signs of AMD. The choroid was significantly thinner on OCT (242 ± 76 µm, P < 0.001) in those with AMD compared with those without (263 ± 63 µm). Subretinal drusenoid deposits were found in 143 eyes (7%); 11 of the 143 eyes (8%) had no other manifestations of AMD. Drusen volume (P < 0.001) and area of geographic atrophy (P < 0.001) were significantly greater, and choroid was significantly (P < 0.001) thinner in subjects with SDDs versus those without SDDs. CONCLUSION: The authors describe spectral domain optical coherence tomography characteristics in an elderly Amish population with and without AMD, including the frequency of SDD. Although relatively uncommon in this population, the authors confirmed that SDDs can be found in the absence of other features of AMD and that eyes with SDDs have thinner choroids.

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.

Genomic loci modulating retinal ganglion cell death following elevated IOP in the mouse.

The present study was designed to identify genomic loci modulating the susceptibility of retinal ganglion cells (RGC) to elevated intraocular pressure (IOP) in the BXD recombinant inbred mouse strain set. IOP was elevated by injecting magnetic microspheres into the anterior chamber and blocking the trabecular meshwork using a handheld magnet to impede drainage. The IOP was then measured over the next 21 days. Only animals with IOP greater than 25 mmHg for two consecutive days or an IOP above 30 mmHg on a single day after microsphere-injection were used in this study. On day 21, mice were sacrificed and the optic nerve was processed for histology. Axons were counted for both the injected and the control eye in 49 BXD strains, totaling 181 normal counts and 191 counts associated with elevated IOP. The axon loss for each strain was calculated and the data were entered into genenetwork.org. The average number of normal axons in the optic nerve across all strains was 54,788 ±â€¯16% (SD), which dropped to 49,545 ±â€¯20% in animals with artificially elevated IOP. Interval mapping demonstrated a relatively similar genome-wide map for both conditions with a suggestive Quantitative Trait Locus (QTL) on proximal Chromosome 3. When the relative axon loss was used to generate a genome-wide interval map, we identified one significant QTL (p < 0.05) on Chromosome 18 between 53.6 and 57 Mb. Within this region, the best candidate gene for modulating axon loss was Aldh7a1. Immunohistochemistry demonstrated ALDH7A1 expression in mouse RGCs. ALDH7A1 variants were not significantly associated with glaucoma in the NEIGHBORHOOD GWAS dataset, but this enzyme was identified as part of the butanoate pathway previously associated with glaucoma risk. Our results suggest that genomic background influences susceptibility to RGC degeneration and death in an inducible glaucoma model.

Pathway analysis by randomization incorporating structure-PARIS: an update.

MOTIVATION: We present an update to the pathway enrichment analysis tool 'Pathway Analysis by Randomization Incorporating Structure (PARIS)' that determines aggregated association signals generated from genome-wide association study results. Pathway-based analyses highlight biological pathways associated with phenotypes. PARIS uses a unique permutation strategy to evaluate the genomic structure of interrogated pathways, through permutation testing of genomic features, thus eliminating many of the over-testing concerns arising with other pathway analysis approaches. RESULTS: We have updated PARIS to incorporate expanded pathway definitions through the incorporation of new expert knowledge from multiple database sources, through customized user provided pathways, and other improvements in user flexibility and functionality. AVAILABILITY AND IMPLEMENTATION: PARIS is freely available to all users at https://ritchielab.psu.edu/software/paris-download CONTACT: jnc43@case.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Meta-analysis of Genome-Wide Association Studies Identifies Novel Loci Associated With Optic Disc Morphology.

Primary open-angle glaucoma is the most common optic neuropathy and an important cause of irreversible blindness worldwide. The optic nerve head or optic disc is divided in two parts: a central cup (without nerve fibers) surrounded by the neuroretinal rim (containing axons of the retinal ganglion cells). The International Glaucoma Genetics Consortium conducted a meta-analysis of genome-wide association studies consisting of 17,248 individuals of European ancestry and 6,841 individuals of Asian ancestry. The outcomes of the genome-wide association studies were disc area and cup area. These specific measurements describe optic nerve morphology in another way than the vertical cup-disc ratio, which is a clinically used measurement, and may shed light on new glaucoma mechanisms. We identified 10 new loci associated with disc area (CDC42BPA, F5, DIRC3, RARB, ABI3BP, DCAF4L2, ELP4, TMTC2, NR2F2, and HORMAD2) and another 10 new loci associated with cup area (DHRS3, TRIB2, EFEMP1, FLNB, FAM101, DDHD1, ASB7, KPNB1, BCAS3, and TRIOBP). The new genes participate in a number of pathways and future work is likely to identify more functions related to the pathogenesis of glaucoma.

Whole exome sequencing of extreme age-related macular degeneration phenotypes.

PURPOSE: Demographic, environmental, and genetic risk factors for age-related macular degeneration (AMD) have been identified; however, a substantial portion of the variance in AMD disease risk and heritability remains unexplained. To identify AMD risk variants and generate hypotheses for future studies, we performed whole exome sequencing for 75 individuals whose phenotype was not well predicted by their genotype at known risk loci. We hypothesized that these phenotypically extreme individuals were more likely to carry rare risk or protective variants with large effect sizes. METHODS: A genetic risk score was calculated in a case-control set of 864 individuals (467 AMD cases, 397 controls) based on 19 common (≥1% minor allele frequency, MAF) single nucleotide variants previously associated with the risk of advanced AMD in a large meta-analysis of advanced cases and controls. We then selected for sequencing 39 cases with bilateral choroidal neovascularization with the lowest genetic risk scores to detect risk variants and 36 unaffected controls with the highest genetic risk score to detect protective variants. After minimizing the influence of 19 common genetic risk loci on case-control status, we targeted single variants of large effect and the aggregate effect of weaker variants within genes and pathways. Single variant tests were conducted on all variants, while gene-based and pathway analyses were conducted on three subsets of data: 1) rare (≤1% MAF in the European population) stop, splice, or damaging missense variants, 2) all rare variants, and 3) all variants. All analyses controlled for the effects of age and sex. RESULTS: No variant, gene, or pathway outside regions known to be associated with risk for advanced AMD reached genome-wide significance. However, we identified several variants with substantial differences in allele frequency between cases and controls with strong additive effects on affection status after controlling for age and sex. Protective effects trending toward significance were detected at two loci identified in single-variant analyses: an intronic variant in FBLN7 (the gene encoding fibulin 7) and at three variants near pyridoxal (pyridoxine, vitamin B6) kinase (PDXK). Aggregate rare-variant analyses suggested evidence for association at ASRGL1, a gene previously linked to photoreceptor cell death, and at BSDC1. In known AMD loci we also identified 29 novel or rare damaging missense or stop/splice variants in our sample of cases and controls. CONCLUSIONS: Identified variants and genes may highlight regions important in the pathogenesis of AMD and are key targets for replication.

Heritability of Choroidal Thickness in the Amish.

PURPOSE: To evaluate the heritability of choroidal thickness and its relationship to age-related macular degeneration (AMD). DESIGN: Cohort study. PARTICIPANTS: Six hundred eighty-nine individuals from Amish families with early or intermediate AMD. METHODS: Ocular coherence tomography was used to quantify choroidal thickness, and fundus photography was used to classify eyes into categories using a modified Clinical Age-Related Maculopathy Staging (CARMS) system. Repeatability and heritability of choroidal thickness and its phenotypic and genetic correlations with the AMD phenotype (CARMS category) were estimated using a generalized linear mixed model (GLMM) approach that accounted for relatedness, repeated measures (left and right eyes), and the effects of age, gender, and refraction. MAIN OUTCOME MEASURES: Heritability of choroidal thickness and its phenotypic and genetic correlation with the AMD phenotype (CARMS category). RESULTS: Phenotypic correlation between choroidal thickness and CARMS category was moderate (Spearman's rank correlation, rs = -0.24; n = 1313 eyes) and significant (GLMM posterior mean, -4.27; 95% credible interval [CI], -7.88 to -0.79; P = 0.02) after controlling for relatedness, age, gender, and refraction. Eyes with advanced AMD had thinner choroids than eyes without AMD (posterior mean, -73.8; 95% CI, -94.7 to -54.6; P < 0.001; n = 1178 eyes). Choroidal thickness was highly repeatable within individuals (repeatability, 0.78; 95% CI, 0.68 to 0.89) and moderately heritable (heritability, 0.40; 95% CI, 0.14 to 0.51), but did not show significant genetic correlation with CARMS category, although the effect size was moderate (genetic correlation, -0.18; 95% CI, -0.49 to 0.16). Choroidal thickness also varied with age, gender, and refraction. The CARMS category showed moderate heritability (heritability, 0.49; 95% CI, 0.26 to 0.72). CONCLUSIONS: We quantify the heritability of choroidal thickness for the first time, highlighting a heritable, quantitative trait that is measurable in all individuals regardless of AMD affection status, and moderately phenotypically correlated with AMD severity. Choroidal thickness therefore may capture variation not captured by the CARMS system. However, because the genetic correlation between choroidal thickness and AMD severity was not significant in our data set, genes associated with the 2 traits may not overlap substantially. Future studies should therefore test for genetic variation associated with choroidal thickness to determine the overlap in genetic basis with AMD.

Global and local ancestry in African-Americans: Implications for Alzheimer's disease risk.

INTRODUCTION: African-American (AA) individuals have a higher risk for late-onset Alzheimer's disease (LOAD) than Americans of primarily European ancestry (EA). Recently, the largest genome-wide association study in AAs to date confirmed that six of the Alzheimer's disease (AD)-related genetic variants originally discovered in EA cohorts are also risk variants in AA; however, the risk attributable to many of the loci (e.g., APOE, ABCA7) differed substantially from previous studies in EA. There likely are risk variants of higher frequency in AAs that have not been discovered. METHODS: We performed a comprehensive analysis of genetically determined local and global ancestry in AAs with regard to LOAD status. RESULTS: Compared to controls, LOAD cases showed higher levels of African ancestry, both globally and at several LOAD relevant loci, which explained risk for AD beyond global differences. DISCUSSION: Exploratory post hoc analyses highlight regions with greatest differences in ancestry as potential candidate regions for future genetic analyses.

Estimating cumulative pathway effects on risk for age-related macular degeneration using mixed linear models.

BACKGROUND: Age-related macular degeneration (AMD) is the leading cause of irreversible visual loss in the elderly in developed countries and typically affects more than 10% of individuals over age 80. AMD has a large genetic component, with heritability estimated to be between 45% and 70%. Numerous variants have been identified and implicate various molecular mechanisms and pathways for AMD pathogenesis but those variants only explain a portion of AMD's heritability. The goal of our study was to estimate the cumulative genetic contribution of common variants on AMD risk for multiple pathways related to the etiology of AMD, including angiogenesis, antioxidant activity, apoptotic signaling, complement activation, inflammatory response, response to nicotine, oxidative phosphorylation, and the tricarboxylic acid cycle. While these mechanisms have been associated with AMD in literature, the overall extent of the contribution to AMD risk for each is unknown. METHODS: In a case-control dataset with 1,813 individuals genotyped for over 600,000 SNPs we used Genome-wide Complex Trait Analysis (GCTA) to estimate the proportion of AMD risk explained by SNPs in genes associated with each pathway. SNPs within a 50 kb region flanking each gene were also assessed, as well as more distant, putatively regulatory SNPs, based on DNaseI hypersensitivity data from ocular tissue in the ENCODE project. RESULTS: We found that 19 previously associated AMD risk SNPs contributed to 13.3% of the risk for AMD in our dataset, while the remaining genotyped SNPs contributed to 36.7% of AMD risk. Adjusting for the 19 risk SNPs, the complement activation and inflammatory response pathways still explained a statistically significant proportion of additional risk for AMD (9.8% and 17.9%, respectively), with other pathways showing no significant effects (0.3% - 4.4%). DISCUSSION: Our results show that SNPs associated with complement activation and inflammation significantly contribute to AMD risk, separately from the risk explained by the 19 known risk SNPs. We found that SNPs within 50 kb regions flanking genes explained additional risk beyond genic SNPs, suggesting a potential regulatory role, but that more distant SNPs explained less than 0.5% additional risk for each pathway. CONCLUSIONS: From these analyses we find that the impact of complement SNPs on risk for AMD extends beyond the established genome-wide significant SNPs.

Advances in the genomics of common eye diseases.

Genome-wide association studies (GWAS) and other genomic technologies have accelerated the discovery of genes and genomic regions contributing to common human ocular disorders with complex inheritance. Age-related macular degeneration (AMD), diabetic retinopathy (DR), glaucoma and myopia account for the majority of visual impairment worldwide. Over 19 genes and/or genomic regions have been associated with AMD. Current investigations are assessing the clinical utility of risk score panels and therapies targeting disease-specific pathways. DR is the leading cause of blindness in the United States and globally is a major cause of vision loss. Genomic investigations have identified molecular pathways associated with DR in animal models which could suggest novel therapeutic targets. Three types of glaucoma, primary-open-angle glaucoma (POAG), angle-closure glaucoma and exfoliation syndrome (XFS) glaucoma, are common age-related conditions. Five genomic regions have been associated with POAG, three with angle-closure glaucoma and one with XFS. Myopia causes substantial ocular morbidity throughout the world. Recent large GWAS have identified >20 associated loci for this condition. In this report, we present a comprehensive overview of the genes and genomic regions contributing to disease susceptibility for these common blinding ocular disorders and discuss the next steps toward translation to effective gene-based screening tests and novel therapies targeting the molecular events contributing to disease.