Phenomewide Association Study of Health Outcomes Associated With the Genetic Correlates of 25 Hydroxyvitamin D Concentration and Vitamin D Binding Protein Concentration.

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.

Cross-ancestry analyses identify new genetic loci associated with 25-hydroxyvitamin D.

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.

Polygenic risk prediction: why and when out-of-sample prediction R2 can exceed SNP-based heritability.

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.

Genetic correlates of vitamin D-binding protein and 25-hydroxyvitamin D in neonatal dried blood spots.

The vitamin D binding protein (DBP), encoded by the group-specific component (GC) gene, is a component of the vitamin D system. In a genome-wide association study of DBP concentration in 65,589 neonates we identify 26 independent loci, 17 of which are in or close to the GC gene, with fine-mapping identifying 2 missense variants on chromosomes 12 and 17 (within SH2B3 and GSDMA, respectively). When adjusted for GC haplotypes, we find 15 independent loci distributed over 10 chromosomes. Mendelian randomization analyses identify a unidirectional effect of higher DBP concentration and (a) higher 25-hydroxyvitamin D concentration, and (b) a reduced risk of multiple sclerosis and rheumatoid arthritis. A phenome-wide association study confirms that higher DBP concentration is associated with a reduced risk of vitamin D deficiency. Our findings provide valuable insights into the influence of DBP on vitamin D status and a range of health outcomes.

Estimation and implications of the genetic architecture of fasting and non-fasting blood glucose.

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 (h2 = 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.

A Comparison of Ten Polygenic Score Methods for Psychiatric Disorders Applied Across Multiple Cohorts.

BACKGROUND: Polygenic scores (PGSs), which assess the genetic risk of individuals for a disease, are calculated as a weighted count of risk alleles identified in genome-wide association studies. PGS methods differ in which DNA variants are included and the weights assigned to them; some require an independent tuning sample to help inform these choices. PGSs are evaluated in independent target cohorts with known disease status. Variability between target cohorts is observed in applications to real data sets, which could reflect a number of factors, e.g., phenotype definition or technical factors. METHODS: The Psychiatric Genomics Consortium Working Groups for schizophrenia and major depressive disorder bring together many independently collected case-control cohorts. We used these resources (31,328 schizophrenia cases, 41,191 controls; 248,750 major depressive disorder cases, 563,184 controls) in repeated application of leave-one-cohort-out meta-analyses, each used to calculate and evaluate PGS in the left-out (target) cohort. Ten PGS methods (the baseline PC+T method and 9 methods that model genetic architecture more formally: SBLUP, LDpred2-Inf, LDpred-funct, LDpred2, Lassosum, PRS-CS, PRS-CS-auto, SBayesR, MegaPRS) were compared. RESULTS: Compared with PC+T, the other 9 methods gave higher prediction statistics, MegaPRS, LDPred2, and SBayesR significantly so, explaining up to 9.2% variance in liability for schizophrenia across 30 target cohorts, an increase of 44%. For major depressive disorder across 26 target cohorts, these statistics were 3.5% and 59%, respectively. CONCLUSIONS: Although the methods that more formally model genetic architecture have similar performance, MegaPRS, LDpred2, and SBayesR rank highest in most comparisons and are recommended in applications to psychiatric disorders.

Genome-wide association study identifies 143 loci associated with 25 hydroxyvitamin D concentration.

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.

Minimal phenotyping yields genome-wide association signals of low specificity for major depression.

Minimal phenotyping refers to the reliance on the use of a small number of self-reported items for disease case identification, increasingly used in genome-wide association studies (GWAS). Here we report differences in genetic architecture between depression defined by minimal phenotyping and strictly defined major depressive disorder (MDD): the former has a lower genotype-derived heritability that cannot be explained by inclusion of milder cases and a higher proportion of the genome contributing to this shared genetic liability with other conditions than for strictly defined MDD. GWAS based on minimal phenotyping definitions preferentially identifies loci that are not specific to MDD, and, although it generates highly predictive polygenic risk scores, the predictive power can be explained entirely by large sample sizes rather than by specificity for MDD. Our results show that reliance on results from minimal phenotyping may bias views of the genetic architecture of MDD and impede the ability to identify pathways specific to MDD.

Effects of interleukin-6 receptor blockade on allergen-induced airway responses in mild asthmatics.

BACKGROUND: Interleukin (IL)-6 signalling has been implicated in allergic asthma by animal, genetic association and clinical studies. In this study, we tested the hypothesis that tocilizumab (TCZ), a human monoclonal antibody that blocks IL-6 signalling, can prevent the development of allergen-induced bronchoconstriction in humans. METHODS: We performed a randomised, double-blind, placebo-controlled study, with eligible participants completing two allergen inhalation challenge tests, conducted before and after treatment with a single dose of TCZ or placebo. The primary efficacy endpoint was the magnitude of the late asthmatic response recorded between 3 and 7 after allergen challenge. The secondary efficacy endpoint was the early asthmatic response, measured 20 min to 2 h after allergen challenge. RESULTS: A total of 66 patients enrolled between September 2014 and August 2017, when the trial was stopped for futility based on results from an interim analysis. Eleven patients fulfilled all eligibility criteria assessed at baseline and were subsequently randomised to the TCZ (n = 6) or placebo (n = 5) groups. Both the primary and secondary efficacy endpoints were not significantly different between the two groups. Five patients reported adverse events (AEs), three in the TCZ group (11 AEs) and two in the placebo group (four AEs). Only one AE was TCZ-related (mild neutropenia), and there were no serious AEs. Significant treatment effects were observed for serum levels of C-reactive protein, IL-6 and soluble IL-6R levels. CONCLUSION: In a small proof-of-concept clinical trial, we found no evidence that a single dose of tocilizumab was able to prevent allergen-induced bronchoconstriction. (Trial registered in the Australian New Zealand Clinical Trials Registry, number ACTRN12614000123640).

Shared genetic origin of asthma, hay fever and eczema elucidates allergic disease biology.

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-8), 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.

Gene-based analysis of regulatory variants identifies 4 putative novel asthma risk genes related to nucleotide synthesis and signaling.

BACKGROUND: Hundreds of genetic variants are thought to contribute to variation in asthma risk by modulating gene expression. Methods that increase the power of genome-wide association studies (GWASs) to identify risk-associated variants are needed. OBJECTIVE: We sought to develop a method that aggregates the evidence for association with disease risk across expression quantitative trait loci (eQTLs) of a gene and use this approach to identify asthma risk genes. METHODS: We developed a gene-based test and software package called EUGENE that (1) is applicable to GWAS summary statistics; (2) considers both cis- and trans-eQTLs; (3) incorporates eQTLs identified in different tissues; and (4) uses simulations to account for multiple testing. We applied this approach to 2 published asthma GWASs (combined n = 46,044) and used mouse studies to provide initial functional insights into 2 genes with novel genetic associations. RESULTS: We tested the association between asthma and 17,190 genes that were found to have cis- and/or trans-eQTLs across 16 published eQTL studies. At an empirical FDR of 5%, 48 genes were associated with asthma risk. Of these, for 37, the association was driven by eQTLs located in established risk loci for allergic disease, including 6 genes not previously implicated in disease cause (eg, LIMS1, TINF2, and SAFB). The remaining 11 significant genes represent potential novel genetic associations with asthma. The association with 4 of these replicated in an independent GWAS: B4GALT3, USMG5, P2RY13, and P2RY14, which are genes involved in nucleotide synthesis or nucleotide-dependent cell activation. In mouse studies, P2ry13 and P2ry14-purinergic receptors activated by adenosine 5-diphosphate and UDP-sugars, respectively-were upregulated after allergen challenge, notably in airway epithelial cells, eosinophils, and neutrophils. Intranasal exposure with receptor agonists induced the release of IL-33 and subsequent eosinophil infiltration into the lungs. CONCLUSION: We identified novel associations between asthma and eQTLs for 4 genes related to nucleotide synthesis/signaling and demonstrated the power of gene-based analyses of GWASs.

Lessons from ten years of genome-wide association studies of asthma.

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.

Allergen-induced IL-6 trans-signaling activates γδ T cells to promote type 2 and type 17 airway inflammation.

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.