Genetic structure in the European American population reflects waves of migration and recent gene flow among different populations. This complex structure can introduce bias in genetic association studies. Using Principal Components Analysis (PCA), we analyze the structure of two independent European American datasets (1,521 individuals-307,315 autosomal SNPs). Individual variation lies across a continuum with some individuals showing high degrees of admixture with non-European populations, as demonstrated through joint analysis with HapMap data. The CEPH Europeans only represent a small fraction of the variation encountered in the larger European American datasets we studied. We interpret the first eigenvector of this data as correlated with ancestry, and we apply an algorithm that we have previously described to select PCA-informative markers (PCAIMs) that can reproduce this structure. Importantly, we develop a novel method that can remove redundancy from the selected SNP panels and show that we can effectively remove correlated markers, thus increasing genotyping savings. Only 150-200 PCAIMs suffice to accurately predict fine structure in European American datasets, as identified by PCA. Simulating association studies, we couple our method with a PCA-based stratification correction tool and demonstrate that a small number of PCAIMs can efficiently remove false correlations with almost no loss in power. The structure informative SNPs that we propose are an important resource for genetic association studies of European Americans. Furthermore, our redundancy removal algorithm can be applied on sets of ancestry informative markers selected with any method in order to select the most uncorrelated SNPs, and significantly decreases genotyping costs.

We conducted a genome-wide association analysis of 7 subfractions of low density lipoproteins (LDLs) and 3 subfractions of intermediate density lipoproteins (IDLs) measured by gradient gel electrophoresis, and their response to statin treatment, in 1868 individuals of European ancestry from the Pharmacogenomics and Risk of Cardiovascular Disease study. Our analyses identified four previously-implicated loci (SORT1, APOE, LPA, and CETP) as containing variants that are very strongly associated with lipoprotein subfractions (log(10)Bayes Factor > 15). Subsequent conditional analyses suggest that three of these (APOE, LPA and CETP) likely harbor multiple independently associated SNPs. Further, while different variants typically showed different characteristic patterns of association with combinations of subfractions, the two SNPs in CETP show strikingly similar patterns--both in our original data and in a replication cohort--consistent with a common underlying molecular mechanism. Notably, the CETP variants are very strongly associated with LDL subfractions, despite showing no association with total LDLs in our study, illustrating the potential value of the more detailed phenotypic measurements. In contrast with these strong subfraction associations, genetic association analysis of subfraction response to statins showed much weaker signals (none exceeding log(10)Bayes Factor of 6). However, two SNPs (in APOE and LPA) previously-reported to be associated with LDL statin response do show some modest evidence for association in our data, and the subfraction response proles at the LPA SNP are consistent with the LPA association, with response likely being due primarily to resistance of Lp(a) particles to statin therapy. An additional important feature of our analysis is that, unlike most previous analyses of multiple related phenotypes, we analyzed the subfractions jointly, rather than one at a time. Comparisons of our multivariate analyses with standard univariate analyses demonstrate that multivariate analyses can substantially increase power to detect associations. Software implementing our multivariate analysis methods is available at http://stephenslab.uchicago.edu/software.html.

Pleiotropic genetic variants have independent effects on different phenotypes. C-reactive protein (CRP) is associated with several cardiometabolic phenotypes. Shared genetic backgrounds may partially underlie these associations. We conducted a genome-wide analysis to identify the shared genetic background of inflammation and cardiometabolic phenotypes using published genome-wide association studies (GWAS). We also evaluated whether the pleiotropic effects of such loci were biological or mediated in nature. First, we examined whether 283 common variants identified for 10 cardiometabolic phenotypes in GWAS are associated with CRP level. Second, we tested whether 18 variants identified for serum CRP are associated with 10 cardiometabolic phenotypes. We used a Bonferroni corrected p-value of 1.1×10-04 (0.05/463) as a threshold of significance. We evaluated the independent pleiotropic effect on both phenotypes using individual level data from the Women Genome Health Study. Evaluating the genetic overlap between inflammation and cardiometabolic phenotypes, we found 13 pleiotropic regions. Additional analyses showed that 6 regions (APOC1, HNF1A, IL6R, PPP1R3B, HNF4A and IL1F10) appeared to have a pleiotropic effect on CRP independent of the effects on the cardiometabolic phenotypes. These included loci where individuals carrying the risk allele for CRP encounter higher lipid levels and risk of type 2 diabetes. In addition, 5 regions (GCKR, PABPC4, BCL7B, FTO and TMEM18) had an effect on CRP largely mediated through the cardiometabolic phenotypes. In conclusion, our results show genetic pleiotropy among inflammation and cardiometabolic phenotypes. In addition to reverse causation, our data suggests that pleiotropic genetic variants partially underlie the association between CRP and cardiometabolic phenotypes.

Aims

To utilize previously reported lead SNPs for low-density lipoprotein cholesterol (LDL-c) levels to find additional loci of importance to statin response, and examine whether genetic predisposition to LDL-c levels associates with differential statin response.

Methods

We investigated effects on statin response of 59 LDL-c SNPs, by combining summary level statistics from the Global Lipids Genetics and Genomic Investigation of Statin Therapy consortia.

Results

Lead SNPs for APOE, SORT1 and NPC1L1 were associated with a decreased LDL-c response to statin treatment, as was overall genetic predisposition for increased LDL-c levels as quantified with 59 SNPs, with a 5.4% smaller statin response per standard deviation increase in genetically raised LDL-c levels.

Conclusion

Genetic predisposition for increased LDL-c level may decrease efficacy of statin therapy.

Linear mixed models are a powerful statistical tool for identifying genetic associations and avoiding confounding. However, existing methods are computationally intractable in large cohorts and may not optimize power. All existing methods require time cost O(MN(2)) (where N is the number of samples and M is the number of SNPs) and implicitly assume an infinitesimal genetic architecture in which effect sizes are normally distributed, which can limit power. Here we present a far more efficient mixed-model association method, BOLT-LMM, which requires only a small number of O(MN) time iterations and increases power by modeling more realistic, non-infinitesimal genetic architectures via a Bayesian mixture prior on marker effect sizes. We applied BOLT-LMM to 9 quantitative traits in 23,294 samples from the Women's Genome Health Study (WGHS) and observed significant increases in power, consistent with simulations. Theory and simulations show that the boost in power increases with cohort size, making BOLT-LMM appealing for genome-wide association studies in large cohorts.

Statins are prescribed widely to lower plasma low-density lipoprotein (LDL) concentrations and cardiovascular disease risk and have been shown to have beneficial effects in a broad range of patients. However, statins are associated with an increased risk, albeit small, of clinical myopathy and type 2 diabetes. Despite evidence for substantial genetic influence on LDL concentrations, pharmacogenomic trials have failed to identify genetic variations with large effects on either statin efficacy or toxicity, and have produced little information regarding mechanisms that modulate statin response. Here we identify a downstream target of statin treatment by screening for the effects of in vitro statin exposure on genetic associations with gene expression levels in lymphoblastoid cell lines derived from 480 participants of a clinical trial of simvastatin treatment. This analysis identified six expression quantitative trait loci (eQTLs) that interacted with simvastatin exposure, including rs9806699, a cis-eQTL for the gene glycine amidinotransferase (GATM) that encodes the rate-limiting enzyme in creatine synthesis. We found this locus to be associated with incidence of statin-induced myotoxicity in two separate populations (meta-analysis odds ratio = 0.60). Furthermore, we found that GATM knockdown in hepatocyte-derived cell lines attenuated transcriptional response to sterol depletion, demonstrating that GATM may act as a functional link between statin-mediated lowering of cholesterol and susceptibility to statin-induced myopathy.

It remains unclear whether the increased risk of new-onset type 2 diabetes (T2D) seen in statin users is due to low LDL-C concentrations, or due to the statin-induced proportional change in LDL-C. In addition, genetic instruments have not been proposed before to examine whether liability to T2D might cause greater proportional statin-induced LDL-C lowering. Using summary-level statistics from the Genomic Investigation of Statin Therapy (GIST, n_max = 40,914) and DIAGRAM (n_max = 159,208) consortia, we found a positive genetic correlation between LDL-C statin response and T2D using LD score regression (r_genetic = 0.36, s.e. = 0.13). However, mendelian randomization analyses did not provide support for statin response having a causal effect on T2D risk (OR 1.00 (95% CI: 0.97, 1.03) per 10% increase in statin response), nor that liability to T2D has a causal effect on statin-induced LDL-C response (0.20% increase in response (95% CI: -0.40, 0.80) per doubling of odds of liability to T2D). Although we found no evidence to suggest that proportional statin response influences T2D risk, a definitive assessment should be made in populations comprised exclusively of statin users, as the presence of nonstatin users in the DIAGRAM dataset may have substantially diluted our effect estimate.

Purpose

Erythropoietic protoporphyria (EPP), characterized by painful cutaneous photosensitivity, results from pathogenic variants in ferrochelatase (FECH). For 96% of patients, EPP results from coinheriting a rare pathogenic variant in trans of a common hypomorphic variant c.315-48T>C (minor allele frequency 0.05). The estimated prevalence of EPP derived from the number of diagnosed individuals in Europe is 0.00092%, but this may be conservative due to underdiagnosis. No study has estimated EPP prevalence using large genetic data sets.

Methods

Disease-associated FECH variants were identified in the UK Biobank, a data set of 500,953 individuals including 49,960 exome sequences. EPP prevalence was then estimated. The association of FECH variants with EPP-related traits was assessed.

Results

Analysis of pathogenic FECH variants in the UK Biobank provides evidence that EPP prevalence is 0.0059% (95% confidence interval [CI]: 0.0042-0.0076%), 1.7-3.0 times more common than previously thought in the UK. In homozygotes for the common c.315-48T>C FECH variant, there was a novel decrement in both erythrocyte mean corpuscular volume (MCV) and hemoglobin.

Conclusion

The prevalence of EPP has been underestimated secondary to underdiagnosis. The common c.315-48T>C allele is associated with both MCV and hemoglobin, an association that could be important both for those with and without EPP.

Background

Statins effectively lower total and plasma LDL-cholesterol, but the magnitude of decrease varies among individuals. To identify single nucleotide polymorphisms (SNPs) contributing to this variation, we performed a combined analysis of genome-wide association (GWA) results from three trials of statin efficacy.

Methods and principal findings

Bayesian and standard frequentist association analyses were performed on untreated and statin-mediated changes in LDL-cholesterol, total cholesterol, HDL-cholesterol, and triglyceride on a total of 3932 subjects using data from three studies: Cholesterol and Pharmacogenetics (40 mg/day simvastatin, 6 weeks), Pravastatin/Inflammation CRP Evaluation (40 mg/day pravastatin, 24 weeks), and Treating to New Targets (10 mg/day atorvastatin, 8 weeks). Genotype imputation was used to maximize genomic coverage and to combine information across studies. Phenotypes were normalized within each study to account for systematic differences among studies, and fixed-effects combined analysis of the combined sample were performed to detect consistent effects across studies. Two SNP associations were assessed as having posterior probability greater than 50%, indicating that they were more likely than not to be genuinely associated with statin-mediated lipid response. SNP rs8014194, located within the CLMN gene on chromosome 14, was strongly associated with statin-mediated change in total cholesterol with an 84% probability by Bayesian analysis, and a p-value exceeding conventional levels of genome-wide significance by frequentist analysis (P = 1.8 x 10(-8)). This SNP was less significantly associated with change in LDL-cholesterol (posterior probability = 0.16, P = 4.0 x 10(-6)). Bayesian analysis also assigned a 51% probability that rs4420638, located in APOC1 and near APOE, was associated with change in LDL-cholesterol.

Conclusions and significance

Using combined GWA analysis from three clinical trials involving nearly 4,000 individuals treated with simvastatin, pravastatin, or atorvastatin, we have identified SNPs that may be associated with variation in the magnitude of statin-mediated reduction in total and LDL-cholesterol, including one in the CLMN gene for which statistical evidence for association exceeds conventional levels of genome-wide significance.

Trial registration

PRINCE and TNT are not registered. CAP is registered at Clinicaltrials.gov NCT00451828.

Background: The pharmacogenetic effect on cardiovascular disease reduction in response to statin treatment has only been assessed in small studies. In a pharmacogenetic genome wide association study (GWAS) analysis within the Genomic Investigation of Statin Therapy (GIST) consortium, we investigated whether genetic variation was associated with the response of statins on cardiovascular disease risk reduction. Methods: The investigated endpoint was incident myocardial infarction (MI) defined as coronary heart disease death and definite and suspect non-fatal MI. For imputed single nucleotide polymorphisms (SNPs), regression analysis was performed on expected allelic dosage and meta-analysed with a fixed-effects model, inverse variance weighted meta-analysis. All SNPs with p-values <5.0 × 10^-4 in stage 1 GWAS meta-analysis were selected for further investigation in stage-2. As a secondary analysis, we extracted SNPs from the Stage-1 GWAS meta-analysis results based on predefined hypotheses to possibly modifying the effect of statin therapy on MI. Results: In stage-1 meta-analysis (eight studies, n = 10,769, 4,212 cases), we observed no genome-wide significant results (p < 5.0 × 10^-8). A total of 144 genetic variants were followed-up in the second stage (three studies, n = 1,525, 180 cases). In the combined meta-analysis, no genome-wide significant hits were identified. Moreover, none of the look-ups of SNPs known to be associated with either CHD or with statin response to cholesterol levels reached Bonferroni level of significance within our stage-1 meta-analysis. Conclusion: This GWAS analysis did not provide evidence that genetic variation affects statin response on cardiovascular risk reduction. It does not appear likely that genetic testing for predicting effects of statins on clinical events will become a useful tool in clinical practice.