Modulation of Genetic Associations with Serum Urate Levels by Body-Mass-Index in Humans

Full lists of genome-wide association results for urate level as analysed in the PLOS-ONE manuscript "Modulation of Genetic Associations with Serum Urate Levels by Body-Mass-Index in Humans”. Please read the READ_DATASHARE.txt for full description

Genome-wide association study identifies _FUT8_ and _ESR2_ as co-regulators of a bi-antennary N-linked glycan A2 (GlcNAc~2~Man~3~GlcNAc~2~) in human plasma proteins

HPLC analysis of N-glycans quantified levels of the biantennary glycan (A2) in plasma proteins of 924 individuals. Subsequent genome-wide association study (GWAS) using 317,503 single nucleotide polymorphysms (SNP) identified two genetic loci influencing variation in A2: FUT 8 and ESR2. We demonstrate that human glycans are amenable to GWAS and their genetic regulation shows sex-specific effects with _FUT 8_ variants explaining 17.3% of the variance in pre-menopausal women, while _ESR2_ variants explained 6.0% of the variance in post-menopausal women

Colocalization of corneal resistance factor GWAS loci with GTEx e/sQTLs highlights plausible candidate causal genes for keratoconus postnatal corneal stroma weakening

Background: Genome-wide association studies (GWAS) for corneal resistance factor (CRF) have identified 100s of loci and proved useful to uncover genetic determinants for keratoconus, a corneal ectasia of early-adulthood onset and common indication of corneal transplantation. In the current absence of studies to probe the impact of candidate causal variants in the cornea, we aimed to fill some of this knowledge gap by leveraging tissue-shared genetic effects.Methods: 181 CRF signals were examined for evidence of colocalization with genetic signals affecting steady-state gene transcription and splicing in adult, non-eye, tissues of the Genotype-Tissue Expression (GTEx) project. Expression of candidate causal genes thus nominated was evaluated in single cell transcriptomes from adult cornea, limbus and conjunctiva. Fine-mapping and colocalization of CRF and keratoconus GWAS signals was also deployed to support their sharing causal variants.Results and discussion: 26.5% of CRF causal signals colocalized with GTEx v8 signals and nominated genes enriched in genes with high and specific expression in corneal stromal cells amongst tissues examined. Enrichment analyses carried out with nearest genes to all 181 CRF GWAS signals indicated that stromal cells of the limbus could be susceptible to signals that did not colocalize with GTEx’s. These cells might not be well represented in GTEx and/or the genetic associations might have context specific effects. The causal signals shared with GTEx provide new insights into mediation of CRF genetic effects, including modulation of splicing events. Functionally relevant roles for several implicated genes’ products in providing tensile strength, mechano-sensing and signaling make the corresponding genes and regulatory variants prime candidates to be validated and their roles and effects across tissues elucidated. Colocalization of CRF and keratoconus GWAS signals strengthened support for shared causal variants but also highlighted many ways into which likely true shared signals could be missed when using readily available GWAS summary statistics

Uncovering Networks from Genome-Wide Association Studies via Circular Genomic Permutation

Genome-wide association studies (GWAS) aim to detect single nucleotide polymorphisms (SNP) associated with trait variation. However, due to the large number of tests, standard analysis techniques impose highly stringent significance thresholds, leaving potentially associated SNPs undetected, and much of the trait genetic variation unexplained. Pathway- and network-based methodologies applied to GWAS aim to detect associations missed by standard single-marker approaches. The complex and non-random architecture of the genome makes it a challenge to derive an appropriate testing framework for such methodologies. We developed a rapid and simple permutation approach that uses GWAS SNP association results to establish the significance of pathway associations while accounting for the linkage disequilibrium structure of SNPs and the clustering of functionally related elements in the genome. All SNPs used in the GWAS are placed in a “circular genome” according to their location. Then the complete set of SNP association P values are permuted by rotation with respect to the genomic locations of the SNPs. Once these “simulated” P values are assigned, the joint gene P values are calculated using Fisher’s combination test, and the association of pathways is tested using the hypergeometric test. The circular genomic permutation approach was applied to a human genome-wide association dataset. The data consists of 719 individuals from the ORCADES study genotyped for ∼300,000 SNPs and measured for 51 traits ranging from physical to biochemical measurements. KEGG pathways (n = 225) were used as the sets of pathways to be tested. Our results demonstrate that the circular genomic permutations provide robust association P values. The non-permuted hypergeometric analysis generates ∼1400 pathway-trait combination results with an association P value more significant than P ≤ 0.05, whereas applying circular genomic permutation reduces the number of significant results to a more credible 40% of that value. The circular permutation software (“genomicper”) is available as an R package at http://cran.r-project.org/

Localising Loci underlying Complex Trait Variation Using Regional Genomic Relationship Mapping

The limited proportion of complex trait variance identified in genome-wide association studies may reflect the limited power of single SNP analyses to detect either rare causative alleles or those of small effect. Motivated by studies that demonstrate that loci contributing to trait variation may contain a number of different alleles, we have developed an analytical approach termed Regional Genomic Relationship Mapping that, like linkage-based family methods, integrates variance contributed by founder gametes within a pedigree. This approach takes advantage of very distant (and unrecorded) relationships, and this greatly increases the power of the method, compared with traditional pedigree-based linkage analyses. By integrating variance contributed by founder gametes in the population, our approach provides an estimate of the Regional Heritability attributable to a small genomic region (e.g. 100 SNP window covering ca. 1 Mb of DNA in a 300000 SNP GWAS) and has the power to detect regions containing multiple alleles that individually contribute too little variance to be detectable by GWAS as well as regions with single common GWAS-detectable SNPs. We use genome-wide SNP array data to obtain both a genome-wide relationship matrix and regional relationship (“identity by state" or IBS) matrices for sequential regions across the genome. We then estimate a heritability for each region sequentially in our genome-wide scan. We demonstrate by simulation and with real data that, when compared to traditional (“individual SNP") GWAS, our method uncovers new loci that explain additional trait variation. We analysed data from three Southern European populations and from Orkney for exemplar traits – serum uric acid concentration and height. We show that regional heritability estimates are correlated with results from genome-wide association analysis but can capture more of the genetic variance segregating in the population and identify additional trait loci

Examining the link between 179 lipid species and 7 diseases using genetic predictors

Background Genome-wide association studies of lipid species have identified several loci shared with various diseases, however, the relationship between lipid species and disease risk remains poorly understood. Here we investigated whether the plasma levels of lipid species are causally linked to disease risk. Methods We built genetic predictors of 179 lipid species, measured in 7174 Finnish individuals, by utilising either 11 high-impact genomic loci or genome-wide polygenic scores (PGS). We assessed the impact of the lipid species on seven diseases by performing disease association across FinnGen (n = 500,348), UK Biobank (n = 420,531), and Generation Scotland (n = 20,032). We performed univariable Mendelian randomizsation (MR) and multivariable MR (MVMR) analyses to examine whether lipid species impact disease risk independently of standard lipids. Findings PGS explained > 4% of the variance for 34 lipid species but variants outside the high-impact loci had only a marginal contribution. Variants within the high-impact loci showed association with all seven diseases. MVMR supported a causal role of ApoB in ischemic heart disease after accounting for lipid species. Phosphatidylethanolamine-increasing LIPC variants seemed to lower age-related macular degeneration risk independently of HDL-cholesterol. MVMR suggested a protective effect of four lipid species containing arachidonic acid on cholelithiasis risk independently of Total Cholesterol. Interpretation Our study demonstrates how genetic predictors of lipid species can be utilised to gain insights into disease risk. We report potential links between lipid species and age-related macular degeneration and cholelithiasis risk, which can be explored for their utility in disease risk prediction and therapy

Local exome sequences facilitate imputation of less common variants and increase power of genome wide association studies

The analysis of less common variants in genome-wide association studies promises to elucidate complex trait genetics but is hampered by low power to reliably detect association. We show that addition of population-specific exome sequence data to global reference data allows more accurate imputation, particularly of less common SNPs (minor allele frequency 1–10%) in two very different European populations. The imputation improvement corresponds to an increase in effective sample size of 28–38%, for SNPs with a minor allele frequency in the range 1–3%

Enrichment of pathogenic alleles in the brittle cornea gene, ZNF469, in keratoconus

Keratoconus, a common inherited ocular disorder resulting in progressive corneal thinning, is the leading indication for corneal transplantation in the developed world. Genome-wide association studies have identified common SNPs 100 kb upstream of ZNF469 strongly associated with corneal thickness. Homozygous mutations in ZNF469 and PR domain-containing protein 5 (PRDM5) genes result in brittle cornea syndrome (BCS) Types 1 and 2, respectively. BCS is an autosomal recessive generalized connective tissue disorder associated with extreme corneal thinning and a high risk of corneal rupture. Some individuals with heterozygous PRDM5 mutations demonstrate a carrier ocular phenotype, which includes a mildly reduced corneal thickness, keratoconus and blue sclera. We hypothesized that heterozygous variants in PRDM5 and ZNF469 predispose to the development of isolated keratoconus. We found a significant enrichment of potentially pathologic heterozygous alleles in ZNF469 associated with the development of keratoconus (P = 0.00102) resulting in a relative risk of 12.0. This enrichment of rare potentially pathogenic alleles in ZNF469 in 12.5% of keratoconus patients represents a significant mutational load and highlights ZNF469 as the most significant genetic factor responsible for keratoconus identified to dat

Multicohort analysis of the maternal age effect on recombination

Several studies have reported that the number of crossovers increases with maternal age in humans, but others have found the opposite. Resolving the true effect has implications for understanding the maternal age effect on aneuploidies. Here, we revisit this question in the largest sample to date using single nucleotide polymorphism (SNP)-chip data, comprising over 6,000 meioses from nine cohorts. We develop and fit a hierarchical model to allow for differences between cohorts and between mothers. We estimate that over 10 years, the expected number of maternal crossovers increases by 2.1% (95% credible interval (0.98%, 3.3%)). Our results are not consistent with the larger positive and negative effects previously reported in smaller cohorts. We see heterogeneity between cohorts that is likely due to chance effects in smaller samples, or possibly to confounders, emphasizing that care should be taken when interpreting results from any specific cohort about the effect of maternal age on recombination

Novel Genetic Locus Influencing Retinal Venular Tortuosity Is Also Associated With Risk of Coronary Artery Disease

Objective: The retina may provide readily accessible imaging biomarkers of global cardiovascular health. Increasing evidence suggests variation in retinal vascular traits is highly heritable. This study aimed to identify the genetic determinants of retinal vascular traits. Approach and Results: We conducted a meta-analysis of genome-wide association studies for quantitative retinal vascular traits derived using semi-automatic image analysis of digital retinal photographs from the GoDARTS (Genetics of Diabetes Audit and Research in Tayside; N=1736) and ORCADES (Orkney Complex Disease Study; N=1358) cohorts. We identified a novel genome-wide significant locus at 19q13 (ACTN4/CAPN12) for retinal venular tortuosity (TortV), and one at 13q34 (COL4A2) for retinal arteriolar tortuosity (TortA); these 2 loci were subsequently confirmed in 3 independent cohorts (Ntotal=1413). In the combined analysis of discovery and replication cohorts, the lead single-nucleotide polymorphism in ACTN4/CAPN12 was rs1808382 (βs.d.=-0.109; SE=0.015; P=2.39×10-13) and in COL4A2 was rs7991229 (βs.d.=0.103; SE=0.015; P=4.66×10-12). Notably, the ACTN4/CAPN12 locus associated with TortV is also associated with coronary artery disease, heart rate, and atrial fibrillation.Conclusions: Genetic determinants of retinal vascular tortuosity are also linked to cardiovascular health. These findings provide a molecular pathophysiological foundation for the use of retinal vascular traits as biomarkers for cardiovascular diseases.</p