Vitamin D deficiency is a candidate risk factor for a range of adverse health outcomes. In a genome-wide association study of 25 hydroxyvitamin D (25OHD) concentration in 417,580 Europeans we ...identify 143 independent loci in 112 1-Mb regions, providing insights into the physiology of vitamin D and implicating genes involved in lipid and lipoprotein metabolism, dermal tissue properties, and the sulphonation and glucuronidation of 25OHD. Mendelian randomization models find no robust evidence that 25OHD concentration has causal effects on candidate phenotypes (e.g. BMI, psychiatric disorders), but many phenotypes have (direct or indirect) causal effects on 25OHD concentration, clarifying the epidemiological relationship between 25OHD status and the health outcomes examined in this study.
Asthma, hay fever (or allergic rhinitis) and eczema (or atopic dermatitis) often coexist in the same individuals, partly because of a shared genetic origin. To identify shared risk variants, we ...performed a genome-wide association study (GWAS; n = 360,838) of a broad allergic disease phenotype that considers the presence of any one of these three diseases. We identified 136 independent risk variants (P < 3 × 10
), including 73 not previously reported, which implicate 132 nearby genes in allergic disease pathophysiology. Disease-specific effects were detected for only six variants, confirming that most represent shared risk factors. Tissue-specific heritability and biological process enrichment analyses suggest that shared risk variants influence lymphocyte-mediated immunity. Six target genes provide an opportunity for drug repositioning, while for 36 genes CpG methylation was found to influence transcription independently of genetic effects. Asthma, hay fever and eczema partly coexist because they share many genetic risk variants that dysregulate the expression of immune-related genes.
Vitamin D status-a complex trait influenced by environmental and genetic factors-is tightly associated with skin colour and ancestry. Yet very few studies have investigated the genetic underpinnings ...of vitamin D levels across diverse ancestries, and the ones that have, relied on small sample sizes, resulting in inconclusive results. Here, we conduct genome-wide association studies (GWAS) of 25 hydroxyvitamin D (25OHD)-the main circulating form of vitamin D-in 442,435 individuals from four broad genetically-determined ancestry groups represented in the UK Biobank: European (N = 421,867), South Asian (N = 9,983), African (N = 8,306) and East Asian (N = 2,279). We identify a new genetic determinant of 25OHD (rs146759773) in individuals of African ancestry, which was not detected in previous analysis of much larger European cohorts due to low minor allele frequency. We show genome-wide significant evidence of dominance effects in 25OHD that protect against vitamin D deficiency. Given that key events in the synthesis of 25OHD occur in the skin and are affected by pigmentation levels, we conduct GWAS of 25OHD stratified by skin colour and identify new associations. Lastly, we test the interaction between skin colour and variants associated with variance in 25OHD levels and identify two loci (rs10832254 and rs1352846) whose association with 25OHD differs in individuals of distinct complexions. Collectively, our results provide new insights into the complex relationship between 25OHD and skin colour and highlight the importance of diversity in genomic studies. Despite the much larger rates of vitamin D deficiency that we and others report for ancestry groups with dark skin (e.g., South Asian), our study highlights the importance of considering ancestral background and/or skin colour when assessing the implications of low vitamin D.
The genetic regulation of post-prandial glucose levels is poorly understood. Here, we characterise the genetic architecture of blood glucose variably measured within 0 and 24 h of fasting in 368,000 ...European ancestry participants of the UK Biobank. We found a near-linear increase in the heritability of non-fasting glucose levels over time, which plateaus to its fasting state value after 5 h post meal (h
= 11%; standard error: 1%). The genetic correlation between different fasting times is > 0.77, suggesting that the genetic control of glucose is largely constant across fasting durations. Accounting for heritability differences between fasting times leads to a ~16% improvement in the discovery of genetic variants associated with glucose. Newly detected variants improve the prediction of fasting glucose and type 2 diabetes in independent samples. Finally, we meta-analysed summary statistics from genome-wide association studies of random and fasting glucose (N = 518,615) and identified 156 independent SNPs explaining 3% of fasting glucose variance. Altogether, our study demonstrates the utility of random glucose measures to improve the discovery of genetic variants associated with glucose homeostasis, even in fasting conditions.
Abstract
Interleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have ...been associated with an increased risk of complex disease pathogenesis. We conducted a two-staged, discovery and replication meta genome-wide association study (GWAS) of circulating serum IL-6 levels comprising up to 67 428 (ndiscovery = 52 654 and nreplication = 14 774) individuals of European ancestry. The inverse variance fixed effects based discovery meta-analysis, followed by replication led to the identification of two independent loci, IL1F10/IL1RN rs6734238 on chromosome (Chr) 2q14, (Pcombined = 1.8 × 10−11), HLA-DRB1/DRB5 rs660895 on Chr6p21 (Pcombined = 1.5 × 10−10) in the combined meta-analyses of all samples. We also replicated the IL6R rs4537545 locus on Chr1q21 (Pcombined = 1.2 × 10−122). Our study identifies novel loci for circulating IL-6 levels uncovering new immunological and inflammatory pathways that may influence IL-6 pathobiology.
Background A variant in the IL-6 receptor (IL-6R) gene increases asthma risk and is predicted to decrease IL-6 classic signaling and increase IL-6 trans-signaling. This suggests that inhibition of ...IL-6 trans-signaling, but not classic signaling, might suppress allergic airway inflammation. Objectives We sought to determine whether IL-6 signaling contributes to (1) acute experimental asthma induced by clinically relevant allergens and (2) variation in asthma clinical phenotypes in asthmatic patients. Methods Mice were sensitized to house dust mite (HDM) or cockroach at day 0, treated with IL-6R inhibitors at day 13, and challenged with the same allergen at days 14 to 17. End points were measured 3 hours after the final challenge. IL-6 and soluble IL-6 receptor (sIL-6R) expression in induced sputum of asthmatic patients was correlated with asthma clinical phenotypes. Results Both HDM and cockroach induced a type 2/type 17 cytokine profile and mixed granulocytic inflammation in the airways. Both allergens increased IL-6 expression in the airways, but only cockroach induced sIL-6R expression. Therefore HDM challenge promoted IL-6 classic signaling but not trans-signaling; in this model treatment with anti–IL-6R did not suppress airway inflammation. In contrast, cockroach-induced inflammation involved activation of IL-6 trans-signaling and production of IL-17A by γδ T cells. Anti–IL-6R, selective blockade of sIL-6R, or γδ T-cell deficiency significantly attenuated cockroach-induced inflammation. Asthmatic patients with high airway IL-6 and sIL-6R levels were enriched for the neutrophilic and mixed granulocytic subtypes. Conclusion Experimental asthma associated with both high IL-6 and high sIL-6R levels in the airways is attenuated by treatment with IL-6R inhibitors.
Twenty‐five genome‐wide association studies (GWAS) of asthma were published between 2007 and 2016, the largest with a sample size of 157242 individuals. Across these studies, 39 genetic variants in ...low linkage disequilibrium (LD) with each other were reported to associate with disease risk at a significance threshold of P<5 × 10−8, including 31 in populations of European ancestry. Results from analyses of the UK Biobank data (n=380 503) indicate that at least 28 of the 31 associations reported in Europeans represent true‐positive findings, collectively explaining 2.5% of the variation in disease liability (median of 0.06% per variant). We identified 49 transcripts as likely target genes of the published asthma risk variants, mostly based on LD with expression quantitative trait loci (eQTL). Of these genes, 16 were previously implicated in disease pathophysiology by functional studies, including TSLP, TNFSF4, ADORA1, CHIT1 and USF1. In contrast, at present, there is limited or no functional evidence directly implicating the remaining 33 likely target genes in asthma pathophysiology. Some of these genes have a known function that is relevant to allergic disease, including F11R, CD247, PGAP3, AAGAB, CAMK4 and PEX14, and so could be prioritized for functional follow‐up. We conclude by highlighting three areas of research that are essential to help translate GWAS findings into clinical research or practice, namely validation of target gene predictions, understanding target gene function and their role in disease pathophysiology and genomics‐guided prioritization of targets for drug development.
Background A variant in the IL-6 receptor (IL-6R) gene increases asthma risk and is predicted to decrease IL-6 classic signaling and increase IL-6 trans-signaling. This suggests that inhibition of ...IL-6 trans-signaling, but not classic signaling, might suppress allergic airway inflammation. Objectives We sought to determine whether IL-6 signaling contributes to (1) acute experimental asthma induced by clinically relevant allergens and (2) variation in asthma clinical phenotypes in asthmatic patients. Methods Mice were sensitized to house dust mite (HDM) or cockroach at day 0, treated with IL-6R inhibitors at day 13, and challenged with the same allergen at days 14 to 17. End points were measured 3 hours after the final challenge. IL-6 and soluble IL-6 receptor (sIL-6R) expression in induced sputum of asthmatic patients was correlated with asthma clinical phenotypes. Results Both HDM and cockroach induced a type 2/type 17 cytokine profile and mixed granulocytic inflammation in the airways. Both allergens increased IL-6 expression in the airways, but only cockroach induced sIL-6R expression. Therefore HDM challenge promoted IL-6 classic signaling but not trans-signaling; in this model treatment with anti-IL-6R did not suppress airway inflammation. In contrast, cockroach-induced inflammation involved activation of IL-6 trans-signaling and production of IL-17A by gamma delta T cells. Anti-IL-6R, selective blockade of sIL-6R, or gamma delta T-cell deficiency significantly attenuated cockroach-induced inflammation. Asthmatic patients with high airway IL-6 and sIL-6R levels were enriched for the neutrophilic and mixed granulocytic subtypes. Conclusion Experimental asthma associated with both high IL-6 and high sIL-6R levels in the airways is attenuated by treatment with IL-6R inhibitors.
While it is known that vitamin D deficiency is associated with adverse bone outcomes, it remains unclear whether low vitamin D status may increase the risk of a wider range of health outcomes. We had ...the opportunity to explore the association between common genetic variants associated with both 25 hydroxyvitamin D (25OHD) and the vitamin D binding protein (DBP, encoded by the GC gene) with a comprehensive range of health disorders and laboratory tests in a large academic medical center. We used summary statistics for 25OHD and DBP to generate polygenic scores (PGS) for 66,482 participants with primarily European ancestry and 13,285 participants with primarily African ancestry from the Vanderbilt University Medical Center Biobank (BioVU). We examined the predictive properties of PGS25OHD, and two scores related to DBP concentration with respect to 1322 health-related phenotypes and 315 laboratory-measured phenotypes from electronic health records. In those with European ancestry: (a) the PGS25OHD and PGSDBP scores, and individual SNPs rs4588 and rs7041 were associated with both 25OHD concentration and 1,25 dihydroxyvitamin D concentrations; (b) higher PGS25OHD was associated with decreased concentrations of triglycerides and cholesterol, and reduced risks of vitamin D deficiency, disorders of lipid metabolism, and diabetes. In general, the findings for the African ancestry group were consistent with findings from the European ancestry analyses. Our study confirms the utility of PGS and two key variants within the GC gene (rs4588 and rs7041) to predict the risk of vitamin D deficiency in clinical settings and highlights the shared biology between vitamin D-related genetic pathways a range of health outcomes.
In polygenic score (PGS) analysis, the coefficient of determination (R2) is a key statistic to evaluate efficacy. R2 is the proportion of phenotypic variance explained by the PGS, calculated in a ...cohort that is independent of the genome-wide association study (GWAS) that provided estimates of allelic effect sizes. The SNP-based heritability (hSNP2, the proportion of total phenotypic variances attributable to all common SNPs) is the theoretical upper limit of the out-of-sample prediction R2. However, in real data analyses R2 has been reported to exceed hSNP2, which occurs in parallel with the observation that hSNP2 estimates tend to decline as the number of cohorts being meta-analyzed increases. Here, we quantify why and when these observations are expected. Using theory and simulation, we show that if heterogeneities in cohort-specific hSNP2 exist, or if genetic correlations between cohorts are less than one, hSNP2 estimates can decrease as the number of cohorts being meta-analyzed increases. We derive conditions when the out-of-sample prediction R2 will be greater than hSNP2 and show the validity of our derivations with real data from a binary trait (major depression) and a continuous trait (educational attainment). Our research calls for a better approach to integrating information from multiple cohorts to address issues of between-cohort heterogeneity.
Display omitted
SNP-based heritability estimates tend to decline and then plateau as the number of cohorts being meta-analyzed in a GWAS increase, and the out-of-sample prediction R2 in "target" cohorts can sometimes exceed its theoretical upper limit. Here, we provide theory to explain these observations that reflect heterogeneity between cohorts in meta-analyses.