Gabriel Cuellar-Partida

Assessment of polygenic effects links primary open-angle glaucoma and age-related macular degeneration.

Primary open-angle glaucoma (POAG) and age-related macular degeneration (AMD) are leading causes of irreversible blindness. Several loci have been mapped using genome-wide association studies. Until very recently, there was no recognized overlap in the genetic contribution to AMD and POAG. At genome-wide significance level, only ABCA1 harbors associations to both diseases. Here, we investigated the genetic architecture of POAG and AMD using genome-wide array data. We estimated the heritability for POAG (h2g = 0.42 ± 0.09) and AMD (h2g = 0.71 ± 0.08). Removing known loci for POAG and AMD decreased the h2g estimates to 0.36 and 0.24, respectively. There was evidence for a positive genetic correlation between POAG and AMD (rg = 0.47 ± 0.25) which remained after removing known loci (rg = 0.64 ± 0.31). We also found that the genetic correlation between sexes for POAG was likely to be less than 1 (rg = 0.33 ± 0.24), suggesting that differences of prevalence among genders may be partly due to heritable factors.

Genome-wide analysis of multi-ancestry cohorts identifies new loci influencing intraocular pressure and susceptibility to glaucoma

Elevated intraocular pressure (IOP) is an important risk factor in developing glaucoma, and variability in IOP might herald glaucomatous development or progression. We report the results of a genome-wide association study meta-analysis of 18 population cohorts from the International Glaucoma Genetics Consortium (IGGC), comprising 35,296 multi-ancestry participants for IOP. We confirm genetic association of known loci for IOP and primary open-angle glaucoma (POAG) and identify four new IOP-associated loci located on chromosome 3q25.31 within the FNDC3B gene (P = 4.19 × 10−8 for rs6445055), two on chromosome 9 (P = 2.80 × 10−11 for rs2472493 near ABCA1 and P = 6.39 × 10−11 for rs8176693 within ABO) and one on chromosome 11p11.2 (best P = 1.04 × 10−11 for rs747782). Separate meta-analyses of 4 independent POAG cohorts, totaling 4,284 cases and 95,560 controls, showed that 3 of these loci for IOP were also associated with POAG.

Epigenetic priors for identifying active transcription factor binding sites

MOTIVATION Accurate knowledge of the genome-wide binding of transcription factors in a particular cell type or under a particular condition is necessary for understanding transcriptional regulation. Using epigenetic data such as histone modification and DNase I, accessibility data has been shown to improve motif-based in silico methods for predicting such binding, but this approach has not yet been fully explored. RESULTS We describe a probabilistic method for combining one or more tracks of epigenetic data with a standard DNA sequence motif model to improve our ability to identify active transcription factor binding sites (TFBSs). We convert each data type into a position-specific probabilistic prior and combine these priors with a traditional probabilistic motif model to compute a log-posterior odds score. Our experiments, using histone modifications H3K4me1, H3K4me3, H3K9ac and H3K27ac, as well as DNase I sensitivity, show conclusively that the log-posterior odds score consistently outperforms a simple binary filter based on the same data. We also show that our approach performs competitively with a more complex method, CENTIPEDE, and suggest that the relative simplicity of the log-posterior odds scoring method makes it an appealing and very general method for identifying functional TFBSs on the basis of DNA and epigenetic evidence. AVAILABILITY AND IMPLEMENTATION FIMO, part of the MEME Suite software toolkit, now supports log-posterior odds scoring using position-specific priors for motif search. A web server and source code are available at http://meme.nbcr.net. Utilities for creating priors are at http://research.imb.uq.edu.au/t.bailey/SD/Cuellar2011. CONTACT t.bailey@uq.edu.au SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

ARHGEF12 influences the risk of glaucoma by increasing intraocular pressure

Primary open-angle glaucoma (POAG) is a blinding disease. Two important risk factors for this disease are a positive family history and elevated intraocular pressure (IOP), which is also highly heritable. Genes found to date associated with IOP and POAG are ABCA1, CAV1/CAV2, GAS7 and TMCO1. However, these genes explain only a small part of the heritability of IOP and POAG. We performed a genome-wide association study of IOP in the population-based Rotterdam Study I and Rotterdam Study II using single nucleotide polymorphisms (SNPs) imputed to 1000 Genomes. In this discovery cohort (n = 8105), we identified a new locus associated with IOP. The most significantly associated SNP was rs58073046 (β = 0.44, P-value = 1.87 × 10(-8), minor allele frequency = 0.12), within the gene ARHGEF12. Independent replication in five population-based studies (n = 7471) resulted in an effect size in the same direction that was significantly associated (β = 0.16, P-value = 0.04). The SNP was also significantly associated with POAG in two independent case-control studies [n = 1225 cases and n = 4117 controls; odds ratio (OR) = 1.53, P-value = 1.99 × 10(-8)], especially with high-tension glaucoma (OR = 1.66, P-value = 2.81 × 10(-9); for normal-tension glaucoma OR = 1.29, P-value = 4.23 × 10(-2)). ARHGEF12 plays an important role in the RhoA/RhoA kinase pathway, which has been implicated in IOP regulation. Furthermore, it binds to ABCA1 and links the ABCA1, CAV1/CAV2 and GAS7 pathway to Mendelian POAG genes (MYOC, OPTN, WDR36). In conclusion, this study identified a novel association between IOP and ARHGEF12.

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.

LocusTrack: integrated visualization of GWAS results and genomic annotation

BackgroundGenome-wide association studies (GWAS) are an important tool for the mapping of complex traits and diseases. Visual inspection of genomic annotations may be used to generate insights into the biological mechanisms underlying GWAS-identified loci.ResultsWe developed LocusTrack, a web-based application that annotates and creates plots of regional GWAS results and incorporates user-specified tracks that display annotations such as linkage disequilibrium (LD), phylogenetic conservation, chromatin state, and other genomic and regulatory elements. Currently, LocusTrack can integrate annotation tracks from the UCSC genome-browser as well as from any tracks provided by the user.ConclusionLocusTrack is an easy-to-use application and can be accessed at the following URL: http://gump.qimr.edu.au/general/gabrieC/LocusTrack/. Users can upload and manage GWAS results and select from and/or provide annotation tracks using simple and intuitive menus. LocusTrack scripts and associated data can be downloaded from the website and run locally.

WNT10A exonic variant increases the risk of keratoconus by decreasing corneal thickness

Keratoconus is a degenerative eye condition which results from thinning of the cornea and causes vision distortion. Treatments such as ultraviolet (UV) cross-linking have proved effective for management of keratoconus when performed in early stages of the disease. The central corneal thickness (CCT) is a highly heritable endophenotype of keratoconus, and it is estimated that up to 95% of its phenotypic variance is due to genetics. Genome-wide association efforts of CCT have identified common variants (i.e. minor allele frequency (MAF) >5%). However, these studies typically ignore the large set of exonic variants whose MAF is usually low. In this study, we performed a CCT exome-wide association analysis in a sample of 1029 individuals from a population-based study in Western Australia. We identified a genome-wide significant exonic variant rs121908120 (P = 6.63 × 10(-10)) in WNT10A. This gene is 437 kb from a gene previously associated with CCT (USP37). We showed in a conditional analysis that the WNT10A variant completely accounts for the signal previously seen at USP37. We replicated our finding in independent samples from the Brisbane Adolescent Twin Study, Twin Eye Study in Tasmania and the Rotterdam Study. Further, we genotyped rs121908120 in 621 keratoconus cases and compared the frequency to a sample of 1680 unscreened controls from the Queensland Twin Registry. We found that rs121908120 increases the risk of keratoconus two times (odds ratio 2.03, P = 5.41 × 10(-5)).

Assessing the Genetic Predisposition of Education on Myopia: A Mendelian Randomization Study

Myopia is the largest cause of uncorrected visual impairments globally and its recent dramatic increase in the population has made it a major public health problem. In observational studies, educational attainment has been consistently reported to be correlated to myopia. Nonetheless, correlation does not imply causation. Observational studies do not tell us if education causes myopia or if instead there are confounding factors underlying the association. In this work, we use a two‐step least squares instrumental‐variable (IV) approach to estimate the causal effect of education on refractive error, specifically myopia. We used the results from the educational attainment GWAS from the Social Science Genetic Association Consortium to define a polygenic risk score (PGRS) in three cohorts of late middle age and elderly Caucasian individuals (N = 5,649). In a meta‐analysis of the three cohorts, using the PGRS as an IV, we estimated that each z‐score increase in education (approximately 2 years of education) results in a reduction of 0.92 ± 0.29 diopters (P = 1.04 × 10−3). Our estimate of the effect of education on myopia was higher (P = 0.01) than the observed estimate (0.25 ± 0.03 diopters reduction per education z‐score [∼2 years] increase). This suggests that observational studies may actually underestimate the true effect. Our Mendelian Randomization (MR) analysis provides new evidence for a causal role of educational attainment on refractive error.

Genetic and environmental factors in conjunctival UV autofluorescence

IMPORTANCE Conjunctival UV autofluorescence (CUVAF) photography was developed to detect and characterize preclinical sunlight-induced ocular damage. Ocular sun exposure has been related to cases of pterygia and was recently negatively correlated with myopia. Hence, CUVAF has excellent potential as an objective biomarker of sun exposure. However, much variation in CUVAF has been observed, and the relative contributions of genes and environment to this variation have not yet been identified. OBJECTIVE To investigate sources of variation in CUVAF in relation to its potential clinical relevance. DESIGN, SETTING, AND PARTICIPANTS We performed a cross-sectional analysis of 3 population-based cohort studies in the general community, including the Twins Eye Study in Tasmania, the Brisbane Adolescent Twin Study, and the Western Australian Pregnancy Cohort (Raine) Study. The twin studies were conducted between 2001 and 2009, and the 20-year follow-up of the Raine Study was completed between March 2010 and February 2012. We included genotypic and phenotypic data from 295 Australian families in the Tasmanian and Brisbane twin studies and from 661 participants in the 20-year follow-up of the Raine Study. We compared CUVAF levels in the 3 cohorts and performed a classic twin study to partition variation in CUVAF. We also conducted a genome-wide association analysis to identify specific genetic variants associated with CUVAF. MAIN OUTCOMES AND MEASURES The total area of CUVAF, heritability of CUVAF, and single-nucleotide polymorphisms (SNPs) associated with CUVAF from the genome-wide association study. RESULTS Within twin cohorts, individuals living closer to the equator (latitude, 27.47° S) had higher levels of CUVAF compared with individuals from southern regions (latitude, 42.88° S) (median [interquartile range], 45.4 [26.8-68.5] vs 28.7 [15.0-42.3] mm2; P < .001). The variation in CUVAF explained by the additive genetic component was 0.37 (95% CI, 0.22-0.56), whereas the variation due to the common environment was 0.50 (95% CI; 0.29-0.71). The SNP rs1060043, located approximately 800 base pairs away from the SLC1A5 gene, a member of the solute carrier family 1, had a genome-wide significant association with a P value of 3.2 × 10-8. Gene-based analysis did not improve our power to detect association with other genes. CONCLUSIONS AND RELEVANCE Our findings confirm that, although a large environmental component to CUVAF (equivalent of sun exposure) exists, genes also play a significant role. We identified a SNP (rs1060043) as being significantly associated with CUVAF; replication of this finding in future studies is warranted.

Genetically low vitamin D concentrations and myopic refractive error: A Mendelian randomization study

Background Myopia prevalence has increased in the past 20 years, with many studies linking the increase to reduced time spent outdoors. A number of recent observational studies have shown an inverse association between vitamin D [25(OH)D] serum levels and myopia. However, in such studies it is difficult to separate the effects of time outdoors and vitamin D levels. In this work we use Mendelian randomization (MR) to assess if genetically determined 25(OH)D levels contribute to the degree of myopia. Methods We performed MR using results from a meta-analysis of refractive error (RE) genome-wide association study (GWAS) that included 37 382 and 8 376 adult participants of European and Asian ancestry, respectively, published by the Consortium for Refractive Error And Myopia (CREAM). We used single nucleotide polymorphisms (SNPs) in the DHCR7, CYP2R1, GC and CYP24A1 genes with known effects on 25(OH)D concentration as instrumental variables (IV). We estimated the effect of 25(OH)D on myopia level using a Wald-type ratio estimator based on the effect estimates from the CREAM GWAS. Results Using the combined effect attributed to the four SNPs, the estimate for the effect of 25(OH)D on refractive error was -0.02 [95% confidence interval (CI) -0.09, 0.04] dioptres (D) per 10 nmol/l increase in 25(OH)D concentration in Caucasians and 0.01 (95% CI -0.17, 0.19) D per 10 nmol/l increase in Asians. Conclusions The tight confidence intervals on our estimates suggest the true contribution of vitamin D levels to degree of myopia is very small and indistinguishable from zero. Previous findings from observational studies linking vitamin D levels to myopia were likely attributable to the effects of confounding by time spent outdoors.