Whole-exome sequencing identifies rare variants in STAB2 associated with venous thromboembolic disease.

Deep vein thrombosis and pulmonary embolism, collectively defined as venous thromboembolism (VTE), are the third leading cause of cardiovascular death in the United States. Common genetic variants conferring increased varying degrees of VTE risk have been identified by genome-wide association studies (GWAS). Rare mutations in the anticoagulant genes PROC, PROS1 and SERPINC1 result in perinatal lethal thrombosis in homozygotes and markedly increased VTE risk in heterozygotes. However, currently described VTE variants account for an insufficient portion of risk to be routinely used for clinical decision making. To identify new rare VTE risk variants, we performed a whole-exome study of 393 individuals with unprovoked VTE and 6114 controls. This study identified 4 genes harboring an excess number of rare damaging variants in patients with VTE: PROS1, STAB2, PROC, and SERPINC1. At STAB2, 7.8% of VTE cases and 2.4% of controls had a qualifying rare variant. In cell culture, VTE-associated variants of STAB2 had a reduced surface expression compared with reference STAB2. Common variants in STAB2 have been previously associated with plasma von Willebrand factor and coagulation factor VIII levels in GWAS, suggesting that haploinsufficiency of stabilin-2 may increase VTE risk through elevated levels of these procoagulants. In an independent cohort, we found higher von Willebrand factor levels and equivalent propeptide levels in individuals with rare STAB2 variants compared with controls. Taken together, this study demonstrates the utility of gene-based collapsing analyses to identify loci harboring an excess of rare variants with functional connections to a complex thrombotic disease.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture.

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF ≤ 1%) variants contribute to complex traits and disease in the general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 × 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 × 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 × 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Polygenic risk for coronary heart disease acts through atherosclerosis in type 2 diabetes.

BACKGROUND: Type 2 diabetes increases the risk of coronary heart disease (CHD), yet the mechanisms involved remain poorly described. Polygenic risk scores (PRS) provide an opportunity to understand risk factors since they reflect etiologic pathways from the entire genome. We therefore tested whether a PRS for CHD influenced risk of CHD in individuals with type 2 diabetes and which risk factors were associated with this PRS. METHODS: We tested the association of a CHD PRS with CHD and its traditional clinical risk factors amongst individuals with type 2 diabetes in UK Biobank (N = 21,102). We next tested the association of the CHD PRS with atherosclerotic burden in a cohort of 352 genome-wide genotyped participants with type 2 diabetes who had undergone coronary angiograms. RESULTS: In the UK Biobank we found that the CHD PRS was strongly associated with CHD amongst individuals with type 2 diabetes (OR per standard deviation increase = 1.50; p = 1.5 × 10- 59). But this CHD PRS was, at best, only weakly associated with traditional clinical risk factors, such as hypertension, hyperlipidemia, glycemic control, obesity and smoking. Conversely, in the angiographic cohort, the CHD PRS was strongly associated with multivessel stenosis (OR = 1.65; p = 4.9 × 10- 4) and increased number of major stenotic lesions (OR = 1.35; p = 9.4 × 10- 3). CONCLUSIONS: Polygenic predisposition to CHD is strongly associated with atherosclerotic burden in individuals with type 2 diabetes and this effect is largely independent of traditional clinical risk factors. This suggests that genetic risk for CHD acts through atherosclerosis with little effect on most traditional risk factors, providing the opportunity to explore new biological pathways.

Large differences in adiponectin levels have no clear effect on multiple sclerosis risk: A Mendelian randomization study.

BACKGROUND: Mendelian randomization (MR) studies have demonstrated strong support for an association between genetically increased body mass index and risk of multiple sclerosis (MS). The adipokine adiponectin may be a potential mechanism linking body mass to risk of MS. OBJECTIVE: To evaluate whether genetically increased adiponectin levels influence risk of MS. METHODS: Using genome-wide significant single nucleotide polymorphisms (SNPs) for adiponectin, we undertook an MR study to estimate the effect of adiponectin on MS. This method prevents bias due to reverse causation and minimizes bias due to confounding. Sensitivity analyses were performed to evaluate the assumptions of MR. RESULTS: MR analyses did not support a role for genetically elevated adiponectin in risk of MS (odds ratio (OR) = 0.93 per unit increase in natural-log-transformed adiponectin, equivalent to a two-standard deviation increase in adiponectin on the absolute scale; 95% confidence interval (CI) = 0.66-1.33; p = 0.61). Further MR analysis suggested that genetic variation at the adiponectin gene, which influences adiponectin level, does not impact MS risk. Sensitivity analyses, including MR-Egger regression, suggested no bias due to pleiotropy. CONCLUSION: Lifelong genetically increased adiponectin levels in humans have no clear effect on risk of MS. Other biological factors driving the association between body mass and MS should be investigated.

Obesity and Multiple Sclerosis: A Mendelian Randomization Study.

BACKGROUND: Observational studies have reported an association between obesity, as measured by elevated body mass index (BMI), in early adulthood and risk of multiple sclerosis (MS). However, bias potentially introduced by confounding and reverse causation may have influenced these findings. Therefore, we elected to perform Mendelian randomization (MR) analyses to evaluate whether genetically increased BMI is associated with an increased risk of MS. METHODS AND FINDINGS: Employing a two-sample MR approach, we used summary statistics from the Genetic Investigation of Anthropometric Traits (GIANT) consortium and the International MS Genetics Consortium (IMSGC), the largest genome-wide association studies for BMI and MS, respectively (GIANT: n = 322,105; IMSGC: n = 14,498 cases and 24,091 controls). Seventy single nucleotide polymorphisms (SNPs) were genome-wide significant (p < 5 x 10-8) for BMI in GIANT (n = 322,105) and were investigated for their association with MS risk in the IMSGC. The effect of each SNP on MS was weighted by its effect on BMI, and estimates were pooled to provide a summary measure for the effect of increased BMI upon risk of MS. Our results suggest that increased BMI influences MS susceptibility, where a 1 standard deviation increase in genetically determined BMI (kg/m2) increased odds of MS by 41% (odds ratio [OR]: 1.41, 95% CI 1.20-1.66, p = 2.7 x 10-5, I2 = 0%, 95% CI 0-29). Sensitivity analyses, including MR-Egger regression, and the weighted median approach provided no evidence of pleiotropic effects. The main study limitations are that, while these sensitivity analyses reduce the possibility that pleiotropy influenced our results, residual pleiotropy is difficult to exclude entirely. CONCLUSION: Genetically elevated BMI is associated with risk of MS, providing evidence for a causal role for obesity in MS etiology. While obesity has been associated with many late-life outcomes, these findings suggest an important consequence of childhood and/or early adulthood obesity.

Mendelian randomisation applied to drug development in cardiovascular disease: a review.

Despite increased expenditure, productivity of the pharmaceutical industry has decreased and currently 90% of developed molecules entering phase II and phase III clinical trials fail to gain regulatory approval. Most of these failures are due to lack of therapeutic efficacy rather than lack of safety, suggesting that drug development failures may often be due to poor drug target validation. Currently, drug targets are largely validated using in vitro assays and animal models which may not translate well to human disease. Emerging methods from human genetics, such as Mendelian randomisation (MR), can enable the validation of putative biomarker drug targets in humans prior to the initiation of clinical trials. MR studies can provide evidence as to whether genetically determined levels of a biomarker influence disease aetiology, enabling investigators to infer whether the biomarker is causal. We review the extent to which MR techniques may be helpful in biomarker validation by assessing the concordance between the results from MR studies and phase III clinical trials for lipid therapy in cardiovascular disease. Our findings show that concordance is highest when MR provides evidence suggesting that a biomarker is not causal. In contrast, there are many examples of clinical trials that still failed despite targeting confirmed causal biomarkers. We discuss why such trials may not succeed, despite evidence for causality in MR studies, and outline important limitations when using MR for biomarker validation in drug development. Nonetheless, given the current inefficiencies in drug development, MR methods have potential to improve the success rate of drug development and ultimately the delivery of new molecules to clinical care.

Vitamin D and Risk of Multiple Sclerosis: A Mendelian Randomization Study.

BACKGROUND: Observational studies have demonstrated an association between decreased vitamin D level and risk of multiple sclerosis (MS); however, it remains unclear whether this relationship is causal. We undertook a Mendelian randomization (MR) study to evaluate whether genetically lowered vitamin D level influences the risk of MS. METHODS AND FINDINGS: We identified single nucleotide polymorphisms (SNPs) associated with 25-hydroxyvitamin D (25OHD) level from SUNLIGHT, the largest (n = 33,996) genome-wide association study to date for vitamin D. Four SNPs were genome-wide significant for 25OHD level (p-values ranging from 6 × 10-10 to 2 × 10-109), and all four SNPs lay in, or near, genes strongly implicated in separate mechanisms influencing 25OHD. We then ascertained their effect on 25OHD level in 2,347 participants from a population-based cohort, the Canadian Multicentre Osteoporosis Study, and tested the extent to which the 25OHD-decreasing alleles explained variation in 25OHD level. We found that the count of 25OHD-decreasing alleles across these four SNPs was strongly associated with lower 25OHD level (n = 2,347, F-test statistic = 49.7, p = 2.4 × 10-12). Next, we conducted an MR study to describe the effect of genetically lowered 25OHD on the odds of MS in the International Multiple Sclerosis Genetics Consortium study, the largest genetic association study to date for MS (including up to 14,498 cases and 24,091 healthy controls). Alleles were weighted by their relative effect on 25OHD level, and sensitivity analyses were performed to test MR assumptions. MR analyses found that each genetically determined one-standard-deviation decrease in log-transformed 25OHD level conferred a 2.0-fold increase in the odds of MS (95% CI: 1.7-2.5; p = 7.7 × 10-12; I2 = 63%, 95% CI: 0%-88%). This result persisted in sensitivity analyses excluding SNPs possibly influenced by population stratification or pleiotropy (odds ratio [OR] = 1.7, 95% CI: 1.3-2.2; p = 2.3 × 10-5; I2 = 47%, 95% CI: 0%-85%) and including only SNPs involved in 25OHD synthesis or metabolism (ORsynthesis = 2.1, 95% CI: 1.6-2.6, p = 1 × 10-9; ORmetabolism = 1.9, 95% CI: 1.3-2.7, p = 0.002). While these sensitivity analyses decreased the possibility that pleiotropy may have biased the results, residual pleiotropy is difficult to exclude entirely. CONCLUSIONS: A genetically lowered 25OHD level is strongly associated with increased susceptibility to MS. Whether vitamin D sufficiency can delay, or prevent, MS onset merits further investigation in long-term randomized controlled trials.

Interleukin-18 as a drug repositioning opportunity for inflammatory bowel disease: A Mendelian randomization study.

Support from human genetics increases the probability of success in drug development. However, few examples exist of successful genomically-driven drug repositioning. Given that a Mendelian form of severe enterocolitis is due to up-regulation of the interleukin-18 (IL18) signaling pathway, and pharmacologic inhibition of IL18 has been shown to reverse this enterocolitis, we undertook a Mendelian randomization study to test the causal effect of elevated IL18 levels on inflammatory bowel disease susceptibility (IBD) in 12,882 cases and 21,770 controls. Mendelian randomization is an established method to assess the role of biomarkers in disease etiology in a manner that minimizes confounding and prevents reverse causation. Using three SNPs that explained almost 7% of the variance in IL18 level, we found that each genetically predicted standard deviation increase in IL18 was associated with an increase in IBD susceptibility (odds ratio = 1.22, 95% CI = 1.11-1.34, P-value = 6 × 10-5). This association was further validated in 25,042 IBD cases and 34,915 controls (odds ratio = 1.13, 95% CI = 1.05-1.20). Recently, an anti-IL18 monoclonal antibody, which decreased free IL18 levels, was found to be safe, yet ineffective in a phase II trial for type 2 diabetes. Taken together, these genomic findings implicated IBD as an alternative indication for anti-IL18 therapy, which should be tested in randomized controlled trials.

Mendelian Randomization Studies Do Not Support a Role for Vitamin D in Coronary Artery Disease.

BACKGROUND: Observational studies support a possible association between decreased vitamin D levels and risk of coronary artery disease (CAD); however, it remains unclear whether this relationship is causal. We aimed to evaluate whether genetically lowered vitamin D levels influence the risk of CAD using a Mendelian randomization approach. METHODS AND RESULTS: In this 2-stage Mendelian randomization study, we first identified single-nucleotide polymorphisms associated with 25-hydroxyvitamin D (25OHD) levels in the SUNLIGHT consortium (n=33 996), then tested them for possible violation of Mendelian randomization assumptions. A count of risk alleles was tested for association with 25OHD levels in a separate cohort (n=2347). Alleles were weighted by their relative effect on 25OHD and tested for their combined effect on CAD in the Coronary Artery Disease Genome-Wide Replication and Meta-Analysis (CARDIoGRAM) study (22 233 cases/64 762 controls). Four single-nucleotide polymorphisms were identified to be associated with 25OHD levels, all in or near genes implicated in 25OHD synthesis, transport or metabolism. A count of these risk alleles was strongly associated with 25OHD (n=2347, F-test statistic=49.7, P=2×10(-12)). None of the single-nucleotide polymorphisms associated with 25OHD levels were associated with CAD (all P values >0.6). The Mendelian randomization odds ratio (OR) for CAD was 0.99 (95% confidence interval, 0.84-1.17; P=0.93; I(2)=0) per SD decrease in log-transformed 25OHD levels. These results persisted after sensitivity analyses for population stratification and pleiotropy. CONCLUSIONS: Genetically lowered 25OHD levels were not associated with increased risk of CAD in a large, well-powered study, suggesting that previous associations between circulating 25OHD levels and CAD are possibly confounded or due to reverse causation.

Genetically decreased vitamin D and risk of Alzheimer disease.

OBJECTIVE: To test whether genetically decreased vitamin D levels are associated with Alzheimer disease (AD) using mendelian randomization (MR), a method that minimizes bias due to confounding or reverse causation. METHODS: We selected single nucleotide polymorphisms (SNPs) that are strongly associated with 25-hydroxyvitamin D (25OHD) levels (p < 5 × 10-8) from the Study of Underlying Genetic Determinants of Vitamin D and Highly Related Traits (SUNLIGHT) Consortium (N = 33,996) to act as instrumental variables for the MR study. We measured the effect of each of these SNPs on 25OHD levels in the Canadian Multicentre Osteoporosis Study (CaMos; N = 2,347) and obtained the corresponding effect estimates for each SNP on AD risk from the International Genomics of Alzheimer's Project (N = 17,008 AD cases and 37,154 controls). To produce MR estimates, we weighted the effect of each SNP on AD by its effect on 25OHD and meta-analyzed these estimates using a fixed-effects model to provide a summary effect estimate. RESULTS: The SUNLIGHT Consortium identified 4 SNPs to be genome-wide significant for 25OHD, which described 2.44% of the variance in 25OHD in CaMos. All 4 SNPs map to genes within the vitamin D metabolic pathway. MR analyses demonstrated that a 1-SD decrease in natural log-transformed 25OHD increased AD risk by 25% (odds ratio 1.25, 95% confidence interval 1.03-1.51, p = 0.021). After sensitivity analysis in which we removed SNPs possibly influenced by pleiotropy and population stratification, the results were largely unchanged. CONCLUSIONS: Our results provide evidence supporting 25OHD as a causal risk factor for AD. These findings provide further rationale to understand the effect of vitamin D supplementation on cognition and AD risk in randomized controlled trials.