Integrative genomic analyses identify candidate causal genes for calcific aortic valve stenosis involving tissue-specific regulation.

There is currently no medical therapy to prevent calcific aortic valve stenosis (CAVS). Multi-omics approaches could lead to the identification of novel molecular targets. Here, we perform a genome-wide association study (GWAS) meta-analysis including 14,819 cases among 941,863 participants of European ancestry. We report 32 genomic loci, among which 20 are novel. RNA sequencing of 500 human aortic valves highlights an enrichment in expression regulation at these loci and prioritizes candidate causal genes. Homozygous genotype for a risk variant near TWIST1, a gene involved in endothelial-mesenchymal transition, has a profound impact on aortic valve transcriptomics. We identify five genes outside of GWAS loci by combining a transcriptome-wide association study, colocalization, and Mendelian randomization analyses. Using cross-phenotype and phenome-wide approaches, we highlight the role of circulating lipoproteins, blood pressure and inflammation in the disease process. Our findings pave the way for the development of novel therapies for CAVS.

A Large-Scale Genome-Wide Study of Gene-Sleep Duration Interactions for Blood Pressure in 811,405 Individuals from Diverse Populations.

Although both short and long sleep duration are associated with elevated hypertension risk, our understanding of their interplay with biological pathways governing blood pressure remains limited. To address this, we carried out genome-wide cross-population gene-by-short-sleep and long-sleep duration interaction analyses for three blood pressure traits (systolic, diastolic, and pulse pressure) in 811,405 individuals from diverse population groups. We discover 22 novel gene-sleep duration interaction loci for blood pressure, mapped to genes involved in neurological, thyroidal, bone metabolism, and hematopoietic pathways. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausibility of distinct influences of both sleep duration extremes in cardiovascular health. With several of our loci reflecting specificity towards population background or sex, our discovery sheds light on the importance of embracing granularity when addressing heterogeneity entangled in gene-environment interactions, and in therapeutic design approaches for blood pressure management.

Genetic Associations and Architecture of Asthma-COPD Overlap.

BACKGROUND: Some people have characteristics of both asthma and COPD (asthma-COPD overlap), and evidence suggests they experience worse outcomes than those with either condition alone. RESEARCH QUESTION: What is the genetic architecture of asthma-COPD overlap, and do the determinants of risk for asthma-COPD overlap differ from those for COPD or asthma? STUDY DESIGN AND METHODS: We conducted a genome-wide association study in 8,068 asthma-COPD overlap case subjects and 40,360 control subjects without asthma or COPD of European ancestry in UK Biobank (stage 1). We followed up promising signals (P < 5 × 10-6) that remained associated in analyses comparing (1) asthma-COPD overlap vs asthma-only control subjects, and (2) asthma-COPD overlap vs COPD-only control subjects. These variants were analyzed in 12 independent cohorts (stage 2). RESULTS: We selected 31 independent variants for further investigation in stage 2, and discovered eight novel signals (P < 5 × 10-8) for asthma-COPD overlap (meta-analysis of stage 1 and 2 studies). These signals suggest a spectrum of shared genetic influences, some predominantly influencing asthma (FAM105A, GLB1, PHB, TSLP), others predominantly influencing fixed airflow obstruction (IL17RD, C5orf56, HLA-DQB1). One intergenic signal on chromosome 5 had not been previously associated with asthma, COPD, or lung function. Subgroup analyses suggested that associations at these eight signals were not driven by smoking or age at asthma diagnosis, and in phenome-wide scans, eosinophil counts, atopy, and asthma traits were prominent. INTERPRETATION: We identified eight signals for asthma-COPD overlap, which may represent loci that predispose to type 2 inflammation, and serious long-term consequences of asthma.

Utility of Genetically Predicted Lp(a) (Lipoprotein [a]) and ApoB Levels for Cardiovascular Risk Assessment.

BACKGROUND: Current lipid guidelines suggest measurement of Lp(a) (lipoprotein[a]) and ApoB (apolipoprotein B) for atherosclerotic cardiovascular disease risk assessment. Polygenic risk scores (PRSs) for Lp(a) and ApoB may identify individuals unlikely to have elevated Lp(a) or ApoB and thus reduce such suggested testing. METHODS: PRSs were developed using least absolute shrinkage and selection operator regression among 273 222 and 356 958 UK Biobank participants of white British ancestry for Lp(a) and ApoB, respectively, and validated in separate sets of 60 771 UK Biobank and 15 050 European Prospective Investigation into Cancer and Nutrition-Norfolk participants. We then assessed the proportion of participants who, based on these PRSs, were unlikely to benefit from Lp(a) or ApoB measurements, according to current lipid guidelines. RESULTS: In the UK Biobank and European Prospective Investigation into Cancer and Nutrition-Norfolk cohorts, the area under the receiver operating curve for the PRS-predicted Lp(a) and ApoB to identify individuals with elevated Lp(a) and ApoB was at least 0.91 (95% CI, 0.90-0.92) and 0.74 (95% CI, 0.73-0.75), respectively. The Lp(a) PRS and measured Lp(a) showed comparable association with atherosclerotic cardiovascular disease incidence, whereas the ApoB PRS was in general less predictive of atherosclerotic cardiovascular disease risk than measured ApoB. In the context of the European Society of Cardiology/European Atherosclerosis Society lipid guidelines, at a 95% sensitivity to identify individuals with elevated Lp(a) and ApoB levels, at least 54% of Lp(a) and 24% of ApoB testing could be reduced by prescreening with a PRS while maintaining a low false-negative rate. CONCLUSIONS: A substantial proportion of suggested testing for elevated Lp(a) and a modest proportion of testing for elevated ApoB could potentially be reduced by prescreening individuals with PRSs.

Genome-wide meta-analysis of muscle weakness identifies 15 susceptibility loci in older men and women.

Low muscle strength is an important heritable indicator of poor health linked to morbidity and mortality in older people. In a genome-wide association study meta-analysis of 256,523 Europeans aged 60 years and over from 22 cohorts we identify 15 loci associated with muscle weakness (European Working Group on Sarcopenia in Older People definition: n = 48,596 cases, 18.9% of total), including 12 loci not implicated in previous analyses of continuous measures of grip strength. Loci include genes reportedly involved in autoimmune disease (HLA-DQA1 p = 4 × 10-17), arthritis (GDF5 p = 4 × 10-13), cell cycle control and cancer protection, regulation of transcription, and others involved in the development and maintenance of the musculoskeletal system. Using Mendelian randomization we report possible overlapping causal pathways, including diabetes susceptibility, haematological parameters, and the immune system. We conclude that muscle weakness in older adults has distinct mechanisms from continuous strength, including several pathways considered to be hallmarks of ageing.

Genome-wide association analysis of type 2 diabetes in the EPIC-InterAct study.

Type 2 diabetes (T2D) is a global public health challenge. Whilst the advent of genome-wide association studies has identified >400 genetic variants associated with T2D, our understanding of its biological mechanisms and translational insights is still limited. The EPIC-InterAct project, centred in 8 countries in the European Prospective Investigations into Cancer and Nutrition study, is one of the largest prospective studies of T2D. Established as a nested case-cohort study to investigate the interplay between genetic and lifestyle behavioural factors on the risk of T2D, a total of 12,403 individuals were identified as incident T2D cases, and a representative sub-cohort of 16,154 individuals was selected from a larger cohort of 340,234 participants with a follow-up time of 3.99 million person-years. We describe the results from a genome-wide association analysis between more than 8.9 million SNPs and T2D risk among 22,326 individuals (9,978 cases and 12,348 non-cases) from the EPIC-InterAct study. The summary statistics to be shared provide a valuable resource to facilitate further investigations into the genetics of T2D.

Meta-analysis investigating the role of interleukin-6 mediated inflammation in type 2 diabetes.

BACKGROUND: Evidence from animal models and observational epidemiology points to a role for chronic inflammation, in which interleukin 6 (IL-6) is a key player, in the pathophysiology of type 2 diabetes (T2D). However, it is unknown whether IL-6 mediated inflammation is implicated in the pathophysiology of T2D. METHODS: We performed a meta-analysis of 15 prospective studies to investigate associations between IL-6 levels and incident T2D including 5,421 cases and 31,562 non-cases. We also estimated the association of a loss-of-function missense variant (Asp358Ala) in the IL-6 receptor gene (IL6R), previously shown to mimic the effects of IL-6R inhibition, in a large trans-ethnic meta-analysis of six T2D case-control studies including 260,614 cases and 1,350,640 controls. FINDINGS: In a meta-analysis of 15 prospective studies, higher levels of IL-6 (per log pg/mL) were significantly associated with a higher risk of incident T2D (1·24 95% CI, 1·17, 1·32; P = 1 × 10-12). In a trans-ethnic meta-analysis of 260,614 cases and 1,350,640 controls, the IL6R Asp358Ala missense variant was associated with lower odds of T2D (OR, 0·98; 95% CI, 0·97, 0·99; P = 2 × 10-7). This association was not due to diagnostic misclassification and was consistent across ethnic groups. IL-6 levels mediated up to 5% of the association between higher body mass index and T2D. INTERPRETATION: Large-scale human prospective and genetic data provide evidence that IL-6 mediated inflammation is implicated in the etiology of T2D but suggest that the impact of this pathway on disease risk in the general population is likely to be small. FUNDING: The EPICNorfolk study has received funding from the Medical Research Council (MRC) (MR/N003284/1, MC-UU_12015/1 and MC_PC_13048) and Cancer Research UK (C864/A14136). The Fenland Study is funded by the MRC (MC_UU_12015/1 and MC_PC_13046).

The association between circulating 25-hydroxyvitamin D metabolites and type 2 diabetes in European populations: A meta-analysis and Mendelian randomisation analysis.

BACKGROUND: Prior research suggested a differential association of 25-hydroxyvitamin D (25(OH)D) metabolites with type 2 diabetes (T2D), with total 25(OH)D and 25(OH)D3 inversely associated with T2D, but the epimeric form (C3-epi-25(OH)D3) positively associated with T2D. Whether or not these observational associations are causal remains uncertain. We aimed to examine the potential causality of these associations using Mendelian randomisation (MR) analysis. METHODS AND FINDINGS: We performed a meta-analysis of genome-wide association studies for total 25(OH)D (N = 120,618), 25(OH)D3 (N = 40,562), and C3-epi-25(OH)D3 (N = 40,562) in participants of European descent (European Prospective Investigation into Cancer and Nutrition [EPIC]-InterAct study, EPIC-Norfolk study, EPIC-CVD study, Ely study, and the SUNLIGHT consortium). We identified genetic variants for MR analysis to investigate the causal association of the 25(OH)D metabolites with T2D (including 80,983 T2D cases and 842,909 non-cases). We also estimated the observational association of 25(OH)D metabolites with T2D by performing random effects meta-analysis of results from previous studies and results from the EPIC-InterAct study. We identified 10 genetic loci associated with total 25(OH)D, 7 loci associated with 25(OH)D3 and 3 loci associated with C3-epi-25(OH)D3. Based on the meta-analysis of observational studies, each 1-standard deviation (SD) higher level of 25(OH)D was associated with a 20% lower risk of T2D (relative risk [RR]: 0.80; 95% CI 0.77, 0.84; p < 0.001), but a genetically predicted 1-SD increase in 25(OH)D was not significantly associated with T2D (odds ratio [OR]: 0.96; 95% CI 0.89, 1.03; p = 0.23); this result was consistent across sensitivity analyses. In EPIC-InterAct, 25(OH)D3 (per 1-SD) was associated with a lower risk of T2D (RR: 0.81; 95% CI 0.77, 0.86; p < 0.001), while C3-epi-25(OH)D3 (above versus below lower limit of quantification) was positively associated with T2D (RR: 1.12; 95% CI 1.03, 1.22; p = 0.006), but neither 25(OH)D3 (OR: 0.97; 95% CI 0.93, 1.01; p = 0.14) nor C3-epi-25(OH)D3 (OR: 0.98; 95% CI 0.93, 1.04; p = 0.53) was causally associated with T2D risk in the MR analysis. Main limitations include the lack of a non-linear MR analysis and of the generalisability of the current findings from European populations to other populations of different ethnicities. CONCLUSIONS: Our study found discordant associations of biochemically measured and genetically predicted differences in blood 25(OH)D with T2D risk. The findings based on MR analysis in a large sample of European ancestry do not support a causal association of total 25(OH)D or 25(OH)D metabolites with T2D and argue against the use of vitamin D supplementation for the prevention of T2D.

Variants associated with HHIP expression have sex-differential effects on lung function.

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P<5x10 -8) interactions with sex on lung function, and 21 showed suggestive interactions (P<1x10 -6). The strongest signal, from rs7697189 (chr4:145436894) on forced expiratory volume in 1 second (FEV 1) (P=3.15x10 -15), was replicated (P=0.016) in SpiroMeta. The C allele increased FEV 1 more in males (untransformed FEV 1 ß=0.028 [SE 0.0022] litres) than females (ß=0.009 [SE 0.0014] litres), and this effect was not accounted for by differential effects on height, smoking or pubertal age. rs7697189 resides upstream of the hedgehog-interacting protein ( HHIP) gene and was previously associated with lung function and HHIP lung expression. We found HHIP expression was significantly different between the sexes (P=6.90x10 -6), but we could not detect sex differential effects of rs7697189 on expression. Conclusions: We identified a novel genotype-by-sex interaction at a putative enhancer region upstream of the HHIP gene. Establishing the mechanism by which HHIP SNPs have different effects on lung function in males and females will be important for our understanding of lung health and diseases in both sexes.

Human Gain-of-Function MC4R Variants Show Signaling Bias and Protect against Obesity.

The melanocortin 4 receptor (MC4R) is a G protein-coupled receptor whose disruption causes obesity. We functionally characterized 61 MC4R variants identified in 0.5 million people from UK Biobank and examined their associations with body mass index (BMI) and obesity-related cardiometabolic diseases. We found that the maximal efficacy of ß-arrestin recruitment to MC4R, rather than canonical Gαs-mediated cyclic adenosine-monophosphate production, explained 88% of the variance in the association of MC4R variants with BMI. While most MC4R variants caused loss of function, a subset caused gain of function; these variants were associated with significantly lower BMI and lower odds of obesity, type 2 diabetes, and coronary artery disease. Protective associations were driven by MC4R variants exhibiting signaling bias toward ß-arrestin recruitment and increased mitogen-activated protein kinase pathway activation. Harnessing ß-arrestin-biased MC4R signaling may represent an effective strategy for weight loss and the treatment of obesity-related cardiometabolic diseases.

Ranking and characterization of established BMI and lipid associated loci as candidates for gene-environment interactions.

Phenotypic variance heterogeneity across genotypes at a single nucleotide polymorphism (SNP) may reflect underlying gene-environment (G×E) or gene-gene interactions. We modeled variance heterogeneity for blood lipids and BMI in up to 44,211 participants and investigated relationships between variance effects (Pv), G×E interaction effects (with smoking and physical activity), and marginal genetic effects (Pm). Correlations between Pv and Pm were stronger for SNPs with established marginal effects (Spearman's ρ = 0.401 for triglycerides, and ρ = 0.236 for BMI) compared to all SNPs. When Pv and Pm were compared for all pruned SNPs, only BMI was statistically significant (Spearman's ρ = 0.010). Overall, SNPs with established marginal effects were overrepresented in the nominally significant part of the Pv distribution (Pbinomial <0.05). SNPs from the top 1% of the Pm distribution for BMI had more significant Pv values (PMann-Whitney = 1.46×10-5), and the odds ratio of SNPs with nominally significant (<0.05) Pm and Pv was 1.33 (95% CI: 1.12, 1.57) for BMI. Moreover, BMI SNPs with nominally significant G×E interaction P-values (Pint<0.05) were enriched with nominally significant Pv values (Pbinomial = 8.63×10-9 and 8.52×10-7 for SNP × smoking and SNP × physical activity, respectively). We conclude that some loci with strong marginal effects may be good candidates for G×E, and variance-based prioritization can be used to identify them.

Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits.

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.

Genome-Wide Association Study of the Modified Stumvoll Insulin Sensitivity Index Identifies BCL2 and FAM19A2 as Novel Insulin Sensitivity Loci.

Genome-wide association studies (GWAS) have found few common variants that influence fasting measures of insulin sensitivity. We hypothesized that a GWAS of an integrated assessment of fasting and dynamic measures of insulin sensitivity would detect novel common variants. We performed a GWAS of the modified Stumvoll Insulin Sensitivity Index (ISI) within the Meta-Analyses of Glucose and Insulin-Related Traits Consortium. Discovery for genetic association was performed in 16,753 individuals, and replication was attempted for the 23 most significant novel loci in 13,354 independent individuals. Association with ISI was tested in models adjusted for age, sex, and BMI and in a model analyzing the combined influence of the genotype effect adjusted for BMI and the interaction effect between the genotype and BMI on ISI (model 3). In model 3, three variants reached genome-wide significance: rs13422522 (NYAP2; P = 8.87 × 10(-11)), rs12454712 (BCL2; P = 2.7 × 10(-8)), and rs10506418 (FAM19A2; P = 1.9 × 10(-8)). The association at NYAP2 was eliminated by conditioning on the known IRS1 insulin sensitivity locus; the BCL2 and FAM19A2 associations were independent of known cardiometabolic loci. In conclusion, we identified two novel loci and replicated known variants associated with insulin sensitivity. Further studies are needed to clarify the causal variant and function at the BCL2 and FAM19A2 loci.

A genomic approach to therapeutic target validation identifies a glucose-lowering GLP1R variant protective for coronary heart disease.

Regulatory authorities have indicated that new drugs to treat type 2 diabetes (T2D) should not be associated with an unacceptable increase in cardiovascular risk. Human genetics may be able to guide development of antidiabetic therapies by predicting cardiovascular and other health endpoints. We therefore investigated the association of variants in six genes that encode drug targets for obesity or T2D with a range of metabolic traits in up to 11,806 individuals by targeted exome sequencing and follow-up in 39,979 individuals by targeted genotyping, with additional in silico follow-up in consortia. We used these data to first compare associations of variants in genes encoding drug targets with the effects of pharmacological manipulation of those targets in clinical trials. We then tested the association of those variants with disease outcomes, including coronary heart disease, to predict cardiovascular safety of these agents. A low-frequency missense variant (Ala316Thr; rs10305492) in the gene encoding glucagon-like peptide-1 receptor (GLP1R), the target of GLP1R agonists, was associated with lower fasting glucose and T2D risk, consistent with GLP1R agonist therapies. The minor allele was also associated with protection against heart disease, thus providing evidence that GLP1R agonists are not likely to be associated with an unacceptable increase in cardiovascular risk. Our results provide an encouraging signal that these agents may be associated with benefit, a question currently being addressed in randomized controlled trials. Genetic variants associated with metabolic traits and multiple disease outcomes can be used to validate therapeutic targets at an early stage in the drug development process.

Physical activity and all-cause mortality across levels of overall and abdominal adiposity in European men and women: the European Prospective Investigation into Cancer and Nutrition Study (EPIC).

BACKGROUND: The higher risk of death resulting from excess adiposity may be attenuated by physical activity (PA). However, the theoretical number of deaths reduced by eliminating physical inactivity compared with overall and abdominal obesity remains unclear. OBJECTIVE: We examined whether overall and abdominal adiposity modified the association between PA and all-cause mortality and estimated the population attributable fraction (PAF) and the years of life gained for these exposures. DESIGN: This was a cohort study in 334,161 European men and women. The mean follow-up time was 12.4 y, corresponding to 4,154,915 person-years. Height, weight, and waist circumference (WC) were measured in the clinic. PA was assessed with a validated self-report instrument. The combined associations between PA, BMI, and WC with mortality were examined with Cox proportional hazards models, stratified by center and age group, and adjusted for sex, education, smoking, and alcohol intake. Center-specific PAF associated with inactivity, body mass index (BMI; in kg/m²) (>30), and WC (≥102 cm for men, ≥88 cm for women) were calculated and combined in random-effects meta-analysis. Life-tables analyses were used to estimate gains in life expectancy for the exposures. RESULTS: Significant interactions (PA × BMI and PA × WC) were observed, so HRs were estimated within BMI and WC strata. The hazards of all-cause mortality were reduced by 16-30% in moderately inactive individuals compared with those categorized as inactive in different strata of BMI and WC. Avoiding all inactivity would theoretically reduce all-cause mortality by 7.35% (95% CI: 5.88%, 8.83%). Corresponding estimates for avoiding obesity (BMI >30) were 3.66% (95% CI: 2.30%, 5.01%). The estimates for avoiding high WC were similar to those for physical inactivity. CONCLUSION: The greatest reductions in mortality risk were observed between the 2 lowest activity groups across levels of general and abdominal adiposity, which suggests that efforts to encourage even small increases in activity in inactive individuals may be beneficial to public health.

Genome-wide meta-analysis identifies six novel loci associated with habitual coffee consumption.

Coffee, a major dietary source of caffeine, is among the most widely consumed beverages in the world and has received considerable attention regarding health risks and benefits. We conducted a genome-wide (GW) meta-analysis of predominately regular-type coffee consumption (cups per day) among up to 91,462 coffee consumers of European ancestry with top single-nucleotide polymorphisms (SNPs) followed-up in ~30 062 and 7964 coffee consumers of European and African-American ancestry, respectively. Studies from both stages were combined in a trans-ethnic meta-analysis. Confirmed loci were examined for putative functional and biological relevance. Eight loci, including six novel loci, met GW significance (log10Bayes factor (BF)>5.64) with per-allele effect sizes of 0.03-0.14 cups per day. Six are located in or near genes potentially involved in pharmacokinetics (ABCG2, AHR, POR and CYP1A2) and pharmacodynamics (BDNF and SLC6A4) of caffeine. Two map to GCKR and MLXIPL genes related to metabolic traits but lacking known roles in coffee consumption. Enhancer and promoter histone marks populate the regions of many confirmed loci and several potential regulatory SNPs are highly correlated with the lead SNP of each. SNP alleles near GCKR, MLXIPL, BDNF and CYP1A2 that were associated with higher coffee consumption have previously been associated with smoking initiation, higher adiposity and fasting insulin and glucose but lower blood pressure and favorable lipid, inflammatory and liver enzyme profiles (P<5 × 10(-8)).Our genetic findings among European and African-American adults reinforce the role of caffeine in mediating habitual coffee consumption and may point to molecular mechanisms underlying inter-individual variability in pharmacological and health effects of coffee.

Novel loci affecting iron homeostasis and their effects in individuals at risk for hemochromatosis.

Variation in body iron is associated with or causes diseases, including anaemia and iron overload. Here, we analyse genetic association data on biochemical markers of iron status from 11 European-population studies, with replication in eight additional cohorts (total up to 48,972 subjects). We find 11 genome-wide-significant (P<5 × 10(-8)) loci, some including known iron-related genes (HFE, SLC40A1, TF, TFR2, TFRC, TMPRSS6) and others novel (ABO, ARNTL, FADS2, NAT2, TEX14). SNPs at ARNTL, TF, and TFR2 affect iron markers in HFE C282Y homozygotes at risk for hemochromatosis. There is substantial overlap between our iron loci and loci affecting erythrocyte and lipid phenotypes. These results will facilitate investigation of the roles of iron in disease.

The interleukin-6 receptor as a target for prevention of coronary heart disease: a mendelian randomisation analysis.

BACKGROUND: A high circulating concentration of interleukin 6 is associated with increased risk of coronary heart disease. Blockade of the interleukin-6 receptor (IL6R) with a monoclonal antibody (tocilizumab) licensed for treatment of rheumatoid arthritis reduces systemic and articular inflammation. However, whether IL6R blockade also reduces risk of coronary heart disease is unknown. METHODS: Applying the mendelian randomisation principle, we used single nucleotide polymorphisms (SNPs) in the gene IL6R to evaluate the likely efficacy and safety of IL6R inhibition for primary prevention of coronary heart disease. We compared genetic findings with the effects of tocilizumab reported in randomised trials in patients with rheumatoid arthritis. FINDINGS: In 40 studies including up to 133,449 individuals, an IL6R SNP (rs7529229) marking a non-synonymous IL6R variant (rs8192284; p.Asp358Ala) was associated with increased circulating log interleukin-6 concentration (increase per allele 9·45%, 95% CI 8·34-10·57) as well as reduced C-reactive protein (decrease per allele 8·35%, 95% CI 7·31-9·38) and fibrinogen concentrations (decrease per allele 0·85%, 95% CI 0·60-1·10). This pattern of effects was consistent with IL6R blockade from infusions of tocilizumab (4-8 mg/kg every 4 weeks) in patients with rheumatoid arthritis studied in randomised trials. In 25,458 coronary heart disease cases and 100,740 controls, the IL6R rs7529229 SNP was associated with a decreased odds of coronary heart disease events (per allele odds ratio 0·95, 95% CI 0·93-0·97, p=1·53×10(-5)). INTERPRETATION: On the basis of genetic evidence in human beings, IL6R signalling seems to have a causal role in development of coronary heart disease. IL6R blockade could provide a novel therapeutic approach to prevention of coronary heart disease that warrants testing in suitably powered randomised trials. Genetic studies in populations could be used more widely to help to validate and prioritise novel drug targets or to repurpose existing agents and targets for new therapeutic uses. FUNDING: UK Medical Research Council; British Heart Foundation; Rosetrees Trust; US National Heart, Lung, and Blood Institute; Du Pont Pharma; Chest, Heart and Stroke Scotland; Wellcome Trust; Coronary Thrombosis Trust; Northwick Park Institute for Medical Research; UCLH/UCL Comprehensive Medical Research Centre; US National Institute on Aging; Academy of Finland; Netherlands Organisation for Health Research and Development; SANCO; Dutch Ministry of Public Health, Welfare and Sports; World Cancer Research Fund; Agentschap NL; European Commission; Swedish Heart-Lung Foundation; Swedish Research Council; Strategic Cardiovascular Programme of the Karolinska Institutet; Stockholm County Council; US National Institute of Neurological Disorders and Stroke; MedStar Health Research Institute; GlaxoSmithKline; Dutch Kidney Foundation; US National Institutes of Health; Netherlands Interuniversity Cardiology Institute of the Netherlands; Diabetes UK; European Union Seventh Framework Programme; National Institute for Healthy Ageing; Cancer Research UK; MacArthur Foundation.

A latent variable partial least squares path modeling approach to regional association and polygenic effect with applications to a human obesity study.

Genetic association studies are now routinely used to identify single nucleotide polymorphisms (SNPs) linked with human diseases or traits through single SNP-single trait tests. Here we introduced partial least squares path modeling (PLSPM) for association between single or multiple SNPs and a latent trait that can involve single or multiple correlated measurement(s). Furthermore, the framework naturally provides estimators of polygenic effect by appropriately weighting trait-attributing alleles. We conducted computer simulations to assess the performance via multiple SNPs and human obesity-related traits as measured by body mass index (BMI), waist and hip circumferences. Our results showed that the associate statistics had type I error rates close to nominal level and were powerful for a range of effect and sample sizes. When applied to 12 candidate regions in data (N = 2,417) from the European Prospective Investigation of Cancer (EPIC)-Norfolk study, a region in FTO was found to have stronger association (rs7204609∼rs9939881 at the first intron P = 4.29×10(-7)) than single SNP analysis (all with P>10(-4)) and a latent quantitative phenotype was obtained using a subset sample of EPIC-Norfolk (N = 12,559). We believe our method is appropriate for assessment of regional association and polygenic effect on a single or multiple traits.