Skip to main content

A phenome-wide association study of genetically mimicked statins



Beyond their success in cardiovascular disease prevention, statins are increasingly recognized to have sex-specific pleiotropic effects. To gain additional insight, we characterized associations of genetically mimicked statins across the phenotype sex-specifically. We also assessed whether any apparently non-lipid effects identified extended to genetically mimicking other widely used lipid modifiers (proprotein convertase subtilisin/kexin type 9 (PCSK9) inhibitors and ezetimibe) or were a consequence of low-density lipoprotein cholesterol (LDL-c).


We performed a sex-specific phenome-wide association study assessing the association of genetic variants in HMGCR, mimicking statins, with 1701 phenotypes. We used Mendelian randomization (MR) to assess if any non-lipid effects found were evident for genetically mimicked PCSK9 inhibitors and ezetimibe or for LDL-c.


As expected, genetically mimicking statins was inversely associated with LDL-c, apolipoprotein B (ApoB), and total cholesterol (TC) and positively associated with glycated hemoglobin (HbA1c) and was related to body composition. Genetically mimicking statins was also inversely associated with serum calcium, sex hormone-binding globulin (SHBG), and platelet count and positively associated with basal metabolic rate (BMR) and mean platelet volume. Stronger associations with genetically mimicked statins were evident for women than men for lipid traits (LDL-c, ApoB, and TC), calcium, and SHBG, but not for platelet attributes, body composition, or BMR. Genetically mimicking PCSK9 inhibitors or ezetimibe was also associated with lower lipids, but was not related to calcium, SHBG, BMR, or body composition. Genetically higher LDL-c increased lipids and decreased BMR, but did not affect calcium, HbA1c, platelet attributes, or SHBG with minor effects on body composition.


Similar inverse associations were found for genetically mimicking statins on lipid traits in men and women as for other lipid modifiers. Besides the positive associations with HbA1c, BMI (which may explain the higher BMR), and aspects of body composition in men and women, genetically mimicking statins was additionally associated with platelet attributes in both sexes and was inversely associated with serum calcium and SHBG in women. This genetic evidence suggests potential pathways that contribute to the effects of statins particularly in women. Further investigation is needed to confirm these findings and their implications for clinical practice.

Peer Review reports


Statins, 3-hydroxy-3-methyl-glutaryl coenzyme A reductase (HMGCR) inhibitors, reduce low-density lipoprotein cholesterol (LDL-c) significantly, resulting in a commeasurable reduction in morbidity and mortality from cardiovascular disease (CVD) [1, 2]. Beyond their effectiveness as a cardiovascular intervention via lipid modification, pleiotropic effects of statins have long been suggested [3,4,5]. Cardiovascular-related pleiotropic effects of statins may include beneficial effects via a range of potentially inter-related factors, including inflammatory responses [6, 7], endothelial function [5], possibly prothrombin time, and sex hormones [8]. More recently, Mendelian randomization (MR) studies have indicated that statins may reduce the risk of cancer by a lipid-independent pathway [9] as well as specifically reducing epithelial ovarian cancer [10]. Taken together, these studies highlight potential sex differences in the mechanisms underlying statins’ protective effects on CVD and overall mortality. However, exactly what mechanisms underlie statins’ complex sex-specific effects and how they might affect the positioning of statins in disease prevention and treatment remains unclear. Randomized trials are rarely designed or powered to investigate mechanisms or pleiotropic effects, although effects of statins on body weight and diabetes have been identified from a meta-analysis of trials [11, 12].

To identify statins’ pleiotropic effects, previous studies have used genetically mimicked statins to assess statins’ metabolomic profile [13], compared it with that of proprotein convertase subtilisin/kexin type 9 (PCSK9) inhibitors [14], and compared statin’s lipoprotein signature with that of cholesteryl ester transfer protein (CETP) inhibitors [15]. However, none of these studies has been comprehensive across the phenotype or sex-specific, when differences by sex are evident for the incidence of CVD [16] and some cancers [17, 18], highlighting the possibility of sex-specific pathways and sex-specific impacts. To fill this gap, we conducted a phenome-wide association study (PheWAS), a summary statistics-based [19] genotype-to-phenotype approach [20, 21], to assess systematically the sex-specific associations of genetically mimicking statins with a wide range of conditions and related phenotypes, using the largest available genome-wide association studies (GWAS), with validation where possible. The use of genetic mimics largely avoids confounding owing to the random allocation of alleles at conception [22]. To assess if any pleiotropic effects identified from the PheWAS were unique to statins, we also assessed whether these pleiotropic effects were evident for (a) genetically mimicking the major lipid modifiers in common use, i.e., PCSK9 inhibitors and ezetimibe, and (b) LDL-c.


Genetic mimics of statins, and other lipid modifiers

A well-established genetic mimic of statins is the T allele of rs12916 in the HMGCR gene [13]. rs12916 downregulates hepatic HMGCR expression and hence lowers serum LDL-c [12]. We used the lead SNPs rs11591147 [23] and rs10260606 [8, 24] to genetically mimic PCSK9 inhibitors and ezetimibe respectively (Additional file 1: Table S1). In sensitivity analysis for significant phenotypes for rs12916, we used all genetic variants mimicking statins taking into account their correlations, and similarly for PCSK9 inhibitors and ezetimibe (correlation coefficient matrixes shown at Additional file 1: Table S2-4). Finally, as previously, we used 56 independent single nucleotide polymorphisms (SNPs) from different genomic regions (r2 < 0.01) to genetically predict LDL-c (Additional file 1: Table S5) [25].

Phenotype data sources

To date, few GWAS have focused on sex-specific analysis. Sex-specific summary genetic associations for many diseases and phenotypes are available from the UK Biobank study of 502,665 people recruited in 2006 to 2010 from Great Britain intended to be aged 40 to 69 years [26]. The UK Biobank study involved a comprehensive baseline assessment including questionnaire, interview assessment, physical measures, and sample collection. Follow-up of the participants for health-related outcomes by record linkage to hospital and death records is ongoing [26]. Here, we used the UK Biobank to perform the primary analysis and replicated the results using other publicly available consortia data, where possible.

Currently, UK Biobank summary statistics curated in the MRC Integrative Epidemiology Unit (IEU) OpenGWAS database [27] include 3948 phenotypes including diagnoses, current health status, treatment records, biochemical assays, anthropometrics, physical measurements, family history, lifestyle, and psychological health, with detailed information on their derivation (source, original questionnaire, or measurement) available on the UK Biobank website ( keyed on the phenotype or ID (Additional file 1: Table S6). Sex-specific summary statistics for the UK Biobank were provided by Neale Lab ( based on up to 361194 white-British people.

Phenotype categorization

Binary outcomes were considered in groups, corresponding to selected International Classification of Diseases (ICD)-9/10 chapters, i.e., circulatory, endocrine, respiratory, neoplasms, digestive, neurological, musculoskeletal, gynecologic and obstetric, hematopoietic, dermatologic, genitourinary, mental health, infectious diseases, sense organs, injuries and poisonings, symptoms, and others. Categorical and ordered phenotypes were classified in groups using the original UK Biobank categories (, i.e., blood count, blood biochemistry, and physical measures.

Inclusion and exclusion criteria

Duplicate phenotypes were excluded. To ensure power, binary phenotypes with less than 100 cases and continuous or categorical ordered phenotypes with a sample size less than 10,000 were excluded [28]. ICD-coded binary phenotypes without main ICD codes or with external causes (codes as Z00–Z99) were also excluded. The original categories provided by the UK Biobank study ( were used. Subcategories unlikely to reflect effects of statins, such as the subcategory covering family history, hospital administration information, household attributes, lifestyle (smoking, diet, physical activity), environmental exposures (pollutant and sun exposure), employment, socio-demographic factors (education, ethnicity), early life factors, maternity experiences, or other factors unrelated to health status, were excluded. Figure 1 shows the selection of phenotypes with detailed lists of the subcategories excluded.

Fig. 1
figure 1

Flow chart of phenotype inclusion for the PheWAS of genetically mimicking statins

Statistical analysis

MR was used to obtain sex-specific effects of genetically mimicking lipid modifiers and of genetically predicted LDL-c scaled to a standard deviation (SD) unit difference in LDL-c [13]. MR estimates were based on Wald estimates (SNP on outcome divided by SNP on exposure) which were meta-analyzed using inverse-variance weighting (IVW) with multiplicative random effects when 3 or more SNPs were available, taking any correlations between SNPs into account obtained from the correlation coefficient matrix of corresponding SNPs. Estimates for binary traits extracted from GWAS using linear regression were converted to odds ratios (OR), for presentation, as necessary, using an established approximation [29]. MR-Egger was used as a sensitivity analysis for LDL-c because it can detect violations of the exclusion restriction assumption and provide a robust estimate when the genetic variants have pleiotropic effects on outcomes [30]. Differences by sex were assessed using a two-sided z-test [31]. The level of statistical significance was adjusted for multiple comparisons using a Bonferroni correction: P = 0.05/Np, where Np is the number of phenotypes tested [32], giving a corrected P value of 2.9 × 10−5 (0.05/1701).


Where possible, associations found in the UK Biobank were validated using appropriate consortia. A two-sided z-test was used to compare the primary and replication results.

Statistical analyses were performed using R (version 4.0.3; R Foundation for Statistical Computing, Vienna, Austria). Associations of selected SNPs with available outcomes were extracted using the TwoSampleMR R package (version 0.5.5). MR analysis was performed using the MendelianRandomization R package (version 0.4.2). Ethical approval is not necessary because this study was conducted using publicly available summary statistics.


In total, 1384 binary, 70 categorical ordered, and 247 continuous phenotypes were included (Additional file 1: Table S6). Of the 1384 binary phenotypes considered, 135 were circulatory, 55 endocrine, 80 respiratory, 69 neoplasms, 143 digestive, and the remaining 902 were neurological, musculoskeletal, gynecologic and obstetric, hematopoietic, dermatologic, genitourinary, mental health, infectious diseases, sense organs, injuries and poisonings, description of symptoms, and others.

Sex-specific PheWAS of genetically mimicked statins

A Manhattan plot of rs12916 (Additional file 2: Figure I) shows −log10 transformed P-values of all phenotypes included by category. Rs12916 had most associations in the endocrine category for binary phenotypes, and in the blood count, blood biochemistry, and physical measures for continuous phenotypes. Correspondingly, the 49 associations found for rs12916 were mainly with endocrine factors, blood counts, blood biochemistry, and anthropometrics (Additional file 1: Table S7).

For the binary phenotypes considered, the rs12916 T allele was mainly associated with a lower risk of high cholesterol and less use of cholesterol-lowering drugs, including atorvastatin, ezetimibe, and simvastatin. The rs12916 T allele was also negatively associated with snoring. These associations were largely similar by sex.

For the continuous phenotypes considered, the rs12916 T allele was associated with lower levels of several lipid traits including LDL-c, apolipoprotein B (ApoB), and total cholesterol, with stronger associations in women than men. The rs12916 T allele was associated with higher glycated hemoglobin (HbA1c), mean platelet volume (MPV), and BMR, and lower values of some other platelet attributes including platelet count and platelet crit, with similar associations by sex. In addition, the rs12916 T allele was also associated with lower serum calcium and lower sex hormone-binding globulin (SHBG) particularly in women (Table 1). The rs12916 T allele was also associated with a wide range of related anthropometric measures, shown in Table 1 for anthropometric measures (higher weight, body mass index (BMI), hip circumference, and waist circumference) and Additional file 1: Table S8 for body composition measures (higher fat percentage, fat mass, and fat-free mass of the arms, legs, trunk, and whole body; higher levels of impedance of the arms, legs, and whole body; and higher basal metabolic rate). These associations were largely similar by sex. Sensitivity analysis also showed consistent directions and similar effect sizes for the lead SNP or using all 6 SNPs with a correlation matrix for statins (Additional file 1: Table S11 and Additional file 2: Figure II).

Table 1 Sex-specific associations of genetically mimicked statins, PCSK9 inhibitors, and ezetimibe (based on rs12916, rs11591147, and rs10260606, respectively) with continuous phenotypes strongly associated with rs12916 (excluding some body composition traits)

Associations for PCSK9 and ezetimibe

Of the 49 associations found for rs12916, only the associations with lipids were evident for PCSK9 inhibitors (based on the rs11591147 T allele), while associations with lower lipids, higher HbA1c, and lower MPV were evident for ezetimibe (based on the rs10260606 Gallele), as shown in Table 1. However, other associations were not evident for either PCSK9 inhibitors or ezetimibe, as shown in Table 1 and Additional file 1: Table S8. Sensitivity analysis also showed consistent directions and similar effect sizes for the lead SNPs or using alternative SNPs with correlation matrixes for PCSK9 inhibitors and ezetimibe (Additional file 1: Table S12-13 and Additional file 2: Figure III and IV).

Associations for LDL-c

Of the 6 associations with binary phenotypes found for rs12916 (Additional file 1: Table S9), genetically predicted LDL-c was associated with higher self-reported high cholesterol and four cholesterol-lowering treatment phenotypes, but not with snoring, where rs12916 was more strongly associated with ezetimibe use in men than women for (P value = 0.003). Of the 43 associations with continuous phenotypes, lipid and body composition traits were a consequence of genetically predicted LDL-c. Specifically, LDL-c was positively associated with ApoB and total cholesterol. LDL-c was also and negatively associated with BMR and with most anthropometric and body composition traits with much smaller magnitude of effects than for statins. LDL-c was not associated with HbA1c, calcium, SHBG, platelet count, platelet crit, or MPV (P values > 0.05). Here, the univariable MR-Egger intercepts for LDL-c on ApoB and waist circumference were significant in both sexes (P = 0.022, 0.038, respectively), suggesting that the IVW estimate was invalid due to potential pleiotropic effects. The MR-Egger estimates for the association of LDL-c with ApoB and waist circumference were similar to the IVW estimates. These associations with continuous phenotypes were generally similar by sex, as shown in Table 2 and Additional file 1: Table S9.

Table 2 Sex-specific MR estimates of genetically predicted LDL-c on continuous phenotypes strongly associated with rs12916 (excluding some body composition traits)


Most of the associations for rs12916 replicated using different GWAS studies (Additional file 1: Table S10), where available, with the same direction of associations as the primary analysis. Most replication results had similar effect sizes to the primary ones, particularly for rs12916 T reducing ApoB, LDL-c, total cholesterol, calcium, platelet count, and platelet crit, and for the positive associations of rs12916 with HbA1c, MPV, BMI, hip circumference, and waist circumference. However, sex-specific GWAS for BMR and most body composition traits is currently not available in other publicly available consortia.


Consistent with previous findings about the effects of statins [13], the rs12916 T allele was associated with lower LDL-c, total cholesterol, and ApoB; higher HbA1c [33]; and greater adiposity [12, 34, 35], particularly higher BMI, with little difference by sex. Our study adds by assessing a wider range of phenotypes sex-specifically. We found the rs12916 T allele was associated with higher BMR, and with platelet attributes in both sexes and with lower SHGB and calcium particularly in women, which did not extend to genetically mimicking the major lipid modifiers in current use and largely were not a consequence of LDL-c.

Few RCTs have assessed the effects of statins on BMR, SHBG, or serum calcium, but statins reducing SHBG in women have been previously reported [36, 37]. Previous MR studies have suggested that calcium increases the risk of ischemic heart disease (IHD) [38,39,40] and SHBG reduces it [40], so these mechanisms together might have a relatively neutral effect on IHD in women. Sex-specific effects of calcium and SHBG on IHD have not been fully assessed, so how these effects would affect specifically women is unknown, although broadly statins appear to have the same effects on IHD in men and women after accounting for testosterone [8]. SHBG inactivates sex hormones, so lowering SHBG might increase the availability of sex hormones and increase the risk of any related conditions. A recent MR study showed that lower SHBG increases the risk of estrogen-positive breast cancer [41]. Effects of SHBG on other cancers in women have not been systematically examined. Concerns were raised about the possibility of statins increasing the risk of breast cancer in women more than a quarter of century ago [42], but are not feasible to investigate in trials. However, a recent MR study did not suggest that genetically mimicking statins increases breast cancer risk [9]. An MR study has suggested that higher BMR increases the risk of colorectal cancer [43]. Recent MR studies have shown that genetically mimicking statins reduces cancer overall, but did not provide sex-specific estimates [9]. Correspondingly, genetically mimicking effects of statins reduced epithelial ovarian cancer [10]. The role of lowering calcium in cancer is unclear and has not been extensively examined.

Few RCTs have examined the effects of statins on platelet attributes. An RCT indicated that statins might reduce platelet count [44], consistent with our findings. We also found that rs12916 T reduced platelet crit, and increased MPV. A cohort study found platelet count was positively associated with CVD risk and mortality [45], as did an MR study [46]. Few RCTs have examined these questions, but MPV is increasingly realized to be important to CVD [47, 48], and platelet crit may be associated with stroke [49], so these may be additional mechanisms by which protective effects of statins are actuated.

Statins increasing BMI [12, 34] and HbA1c [33, 50] have been reported previously, consistent with our findings. We also found that the rs12916 T allele affected BMR, fat-mass, and fat-free mass. Genetically mimicked PCSK9 inhibitors and ezetimibe were not associated with BMI or BMR, and genetically mimicked PCSK9 inhibitors were not associated with HbA1c, consistent with previous findings [51]. LDL-c had minor effects on BMR, anthropometrics, and most body composition traits. Given statins adversely impact body composition and glycemic traits more strongly than other lipid modifiers but have similar effects on IHD per unit change in LDL-c, it suggests that statins have greater effects via lipid and/or non-lipid mechanisms than other lipid modifiers.

Differences in pleiotropic effects of genetically mimicked statins compared to PCSK9 inhibitors and ezetimibe may be related to the differences of their mechanisms of action. PCSK9 inhibitors reduce the degradation of LDL receptors by blocking PCSK9 [52], resulting in lower levels of circulating LDL-c [53]. Ezetimibe only inhibits the absorption of cholesterol [54]. Statins target cholesterol synthesis [55], including de novo synthesis in Leydig cells, so statins specifically affect hormones, although the mechanism by which statins might affect SHBG in women is less clear.

This study aimed to identify pleiotropic effects of statins, and whether they might be mediated by the target of statins, LDL-c, by testing the associations of LDL-c with any pleiotropic effects. We found the pleiotropic effects of statins did not appear to be driven by LDL-c, and so are specific effects of statins. Exactly, what drives these pleiotropic effects of statins has not been definitively established, nor has their interrelationships, which may well be complex. A possibility is that the pleiotropic effects of statins are driven by effects of statins on BMI. However, statins raise BMI in both sexes (Table 1), while genetically mimicked effects on calcium and SHBG were specific to women, suggesting a more complex explanation. A trial of gastric bypass suggested that BMI increases platelet count [56] but genetically predicted BMI does not appear to affect platelet count in women (Additional file 1: Table S14). However, BMI is well-known to play a crucial role in BMR in both sexes [57]. Using multivariable MR, the pleiotropic effects of genetically mimicking statins on BMR were not independent of BMI (Additional file 1: Table S15), suggesting that any effects of statins on BMR are due to statins raising BMI.

Despite a comprehensive sex-specific scan in the largest available studies, this study has some limitations. First, not all phenotypes of interest were available for the main analysis or for replication, such as very-low-density lipoprotein and some other lipid sub-fractions, and some body composition traits; however, effects of statins on lipid fractions have been examined before [13,14,15]. Second, this study is systematic and comprehensive but is also agnostic. As such, it uses a stringent test for significance to avoid chance findings, so may not replicate all known effects of statins. Instead, this study may provide information about unknown or overlooked potential effects of statins, which is important because of the very widespread use of statins globally. The multiple comparison cut-off using a Bonferroni correction is suitable for agnostic studies, as here. An agnostic study design is most appropriate for identifying pleiotropic effects which have not been considered before rather than for replicating findings based on known physiological pathways, which can be evaluated on different criteria in the context of all the other evidence. The largest available GWAS of IHD does not show rs12916 associated with IHD at genome-wide significance [58], which does not invalidate the GWAS or the role of statins in preventing IHD; instead, the GWAS provides information about other, possibly overlooked, factors that could be relevant to IHD. Similarly, here, this study does not invalidate known relations of statins with IHD or testosterone, which have been demonstrated in a meta-analysis of RCTs [59, 60], but provides additional insight about potentially relevant pleiotropic effects of statins. Overall, this study design is most appropriate for identifying pleiotropic effects of statins rather than replicating known effects. We cannot exclude the possibility that some novel effects of statins have been missed, which could be addressed by repeating this study when larger sex-specific genetic studies are available. Third, we used rs12916 and associated genetic variants as a surrogate for the pharmacological effects of statins [12], which mimics a life-long small dose of endogenous statins [8, 10, 13], so the MR estimates represent life-long inhibition of HMGCR [8, 10] and do not necessarily reflect the effects of statin treatment which generally starts in middle age [10, 61]. These estimates are usually different in magnitude from the short-term effects of pharmacologic interventions in an RCT [62] although similar effect sizes have been seen for genetically mimicked statins and use of statins [13]. Fourth, the UK Biobank is not representative of the UK population. However, no confounding and no selection bias are the criteria for an internally valid study of associations, not population representativeness [63]. The UK Biobank has shown similar associations to an equivalent population-representative study [64]. However, as with any study recruited in middle to older age, the UK Biobank is missing people who died before recruitment from their genetic make-up, from a condition of interest, or from a competing risk of such a condition, which may generate selection bias particularly for conditions that share etiology with diseases that cause death before recruitment [65], so effects of statins on late-onset diseases may have been missed. The UK Biobank study undoubtedly comprises healthy volunteers, less vulnerable to disease, which may bias associations with disease towards null. Fifth, the underlying studies mostly concern populations of European descent, due to data availability, which may limit the generalizability of these findings to other populations. It would be extremely beneficial to validate these findings in consortia more representative of the global population.


After systematic examination, we found that genetically mimicking statins was similarly associated with lower lipids as other major lipid modifiers and was positively associated with adiposity and HbA1c as expected. We also found genetically mimicking statins was associated with higher BMR, consistent with effects on adiposity. In addition, genetically mimicking statins was also associated with lower serum calcium and SHBG and was related to platelet attributes. Associations with SHBG and calcium were specific to women. Besides lipid modulation, inferred from the genetic evidence obtained, lower platelet count may contribute to the benefits of statins, whereas serum calcium could be in a potential pathway that contributes to the benefits of statins on CVD in specifically women. The consequences of statins possibly lowering SHBG on hormone-related conditions in women require further investigation, initially using study designs, such as Mendelian randomization which can provide high-quality evidence expeditiously without any risk. Such information could identify whether statins’ effects on SHBG in women are likely to make any material difference to the net benefits of statins, so as to further inform the level of future risk at which statin use should be initiated in women for primary prevention.

Availability of data and materials

The datasets generated and analyzed during the current study are available in the UK Biobank repository with sex-specific datasets available, [], and the TwoSampleMR R package, [].



Apolipoprotein B


Body mass index


Basal metabolic rate


Coronary Artery Disease Genome wide Replication and Meta-analysis plus The Coronary Artery Disease Genetics


Cholesteryl ester transfer protein


Cardiovascular disease


The Global Lipids Genetics Consortium


Genome-wide association studies


Glycated hemoglobin


International Classification of Diseases


Integrative Epidemiology Unit


Ischemic heart disease


Inverse-variance weighting


Low-density lipoprotein cholesterol


Mean platelet volume


Mendelian randomization


Odds ratio


Proprotein convertase subtilisin/kexin type 9


Phenome-wide association study


Standard deviation


Sex hormone-binding globulin


Single nucleotide polymorphism


Total cholesterol


  1. Michos ED, McEvoy JW, Blumenthal RS. Lipid management for the prevention of atherosclerotic cardiovascular disease. N Engl J Med. 2019;381(16):1557–67.

    Article  CAS  PubMed  Google Scholar 

  2. Ray KK, Corral P, Morales E, Nicholls SJ. Pharmacological lipid-modification therapies for prevention of ischaemic heart disease: current and future options. Lancet. 2019;394(10199):697–708.

    Article  CAS  Google Scholar 

  3. Mitropoulos KA, Armitage JM, Collins R, Meade TW, Reeves BE, Wallendszus KR, et al. Randomized placebo-controlled study of the effects of simvastatin on haemostatic variables, lipoproteins and free fatty acids. The Oxford Cholesterol Study Group. Eur Heart J. 1997;18(2):235–41.

    Article  CAS  PubMed  Google Scholar 

  4. Rosenson RS, Tangney CC. Antiatherothrombotic properties of statins: implications for cardiovascular event reduction. JAMA. 1998;279(20):1643–50.

    Article  CAS  PubMed  Google Scholar 

  5. O'Driscoll G, Green D, Taylor RR. Simvastatin, an HMG-coenzyme A reductase inhibitor, improves endothelial function within 1 month. Circulation. 1997;95(5):1126–31.

    Article  CAS  PubMed  Google Scholar 

  6. Ridker PM, Rifai N, Pfeffer MA, Sacks FM, Moye LA, Goldman S, et al. Inflammation, pravastatin, and the risk of coronary events after myocardial infarction in patients with average cholesterol levels. Cholesterol and Recurrent Events (CARE) Investigators. Circulation. 1998;98(9):839–44.

    Article  CAS  PubMed  Google Scholar 

  7. Albert MA, Danielson E, Rifai N, Ridker PM. Effect of statin therapy on C-reactive protein levels: the pravastatin inflammation/CRP evaluation (PRINCE): a randomized trial and cohort study. JAMA. 2001;286(1):64–70.

    Article  CAS  PubMed  Google Scholar 

  8. Schooling CM, Zhao JV, Au Yeung SL, Leung GM. Investigating pleiotropic effects of statins on ischemic heart disease in the UK Biobank using Mendelian randomisation. Elife. 2020;9:e58567.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Carter P, Vithayathil M, Kar S, Potluri R, Mason AM, Larsson SC, et al. Predicting the effect of statins on cancer risk using genetic variants from a Mendelian randomization study in the UK Biobank. Elife. 2020;9:e57191.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Yarmolinsky J, Bull CJ, Vincent EE, Robinson J, Walther A, Smith GD, et al. Association between genetically proxied inhibition of HMG-CoA reductase and epithelial ovarian cancer. JAMA. 2020;323(7):646–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Liao JK, Laufs U. Pleiotropic effects of statins. Annu Rev Pharmacol Toxicol. 2005;45(1):89–118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Swerdlow DI, Preiss D, Kuchenbaecker KB, Holmes MV, Engmann JE, Shah T, et al. HMG-coenzyme A reductase inhibition, type 2 diabetes, and bodyweight: evidence from genetic analysis and randomised trials. Lancet. 2015;385(9965):351–61.

    Article  CAS  Google Scholar 

  13. Würtz P, Wang Q, Soininen P, Kangas AJ, Fatemifar G, Tynkkynen T, et al. Metabolomic profiling of statin use and genetic inhibition of HMG-CoA reductase. J Am Coll Cardiol. 2016;67(10):1200–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Sliz E, Kettunen J, Holmes MV, Williams CO, Boachie C, Wang Q, et al. Metabolomic consequences of genetic inhibition of PCSK9 compared with statin treatment. Circulation. 2018;138(22):2499–512.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Kettunen J, Holmes MV, Allara E, Anufrieva O, Ohukainen P, Oliver-Williams C, et al. Lipoprotein signatures of cholesteryl ester transfer protein and HMG-CoA reductase inhibition. Plos Biol. 2019;17(12):e3000572.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Sytkowski PA, D'Agostino RB, Belanger A, Kannel WB. Sex and time trends in cardiovascular disease incidence and mortality: the Framingham Heart Study, 1950-1989. Am J Epidemiol. 1996;143(4):338–50.

    Article  CAS  PubMed  Google Scholar 

  17. Lou L, Wang L, Zhang Y, Chen G, Lin L, Jin X, et al. Sex difference in incidence of gastric cancer: an international comparative study based on the Global Burden of Disease Study 2017. BMJ Open. 2020;10(1):e033323.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Koti M, Ingersoll MA, Gupta S, Lam CM, Li X, Kamat AM, et al. Sex differences in bladder cancer immunobiology and outcomes: a collaborative review with implications for treatment. Eur Urol Oncol. 2020;3(5):622–30.

    Article  PubMed  Google Scholar 

  19. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7:e34408.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Ye Z, Mayer J, Ivacic L, Zhou Z, He M, Schrodi SJ, et al. Phenome-wide association studies (PheWASs) for functional variants. Eur J Hum Genet. 2015;23(4):523–9.

    Article  CAS  PubMed  Google Scholar 

  21. Arguinano AA, Ndiaye NC, Masson C, Visvikis-Siest S. Pleiotropy of ABO gene: correlation of rs644234 with E-selectin and lipid levels. Clin Chem Lab Med. 2018;56(5):748–54.

    Article  CAS  PubMed  Google Scholar 

  22. Smith GD, Ebrahim S. ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32(1):1–22.

    Article  PubMed  Google Scholar 

  23. Ference BA, Robinson JG, Brook RD, Catapano AL, Chapman MJ, Neff DR, et al. Variation in PCSK9 and HMGCR and risk of cardiovascular disease and diabetes. N Engl J Med. 2016;375(22):2144–53.

    Article  CAS  PubMed  Google Scholar 

  24. Klarin D, O'Donnell CJ, Kathiresan S. Mendelian randomization study of ACLY and cardiovascular disease. N Engl J Med. 2020;383(7):e50.

    Article  PubMed  Google Scholar 

  25. Lotta LA, Stewart ID, Sharp SJ, Day FR, Burgess S, Luan J, et al. Association of genetically enhanced lipoprotein lipase-mediated lipolysis and low-density lipoprotein cholesterol-lowering alleles with risk of coronary disease and type 2 diabetes. JAMA Cardiol. 2018;3(10):957–66.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. Plos Med. 2015;12(3):e1001779.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Elsworth B, Lyon M, Alexander T, Liu Y, Matthews P, Hallett J, Bates P, Palmer T, Haberland V, Smith GD et al. The MRC IEU OpenGWAS data infrastructure. bioRxiv. 2020.08.10.244293v1.

  28. Delongchamp R, Faramawi MF, Feingold E, Chung D, Abouelenein S. The association between SNPs and a quantitative trait: power calculation. Eur J Environ Public Health. 2018;2(2):10.

    PubMed  PubMed Central  Google Scholar 

  29. Lloyd-Jones LR, Robinson MR, Yang J, Visscher PM. Transformation of summary statistics from linear mixed model association on all-or-none traits to odds ratio. Genetics. 2018;208(4):1397–408.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Altman DG, Bland JM. Interaction revisited: the difference between two estimates. BMJ. 2003;326(7382):219.

    Article  Google Scholar 

  32. Read RW, Schlauch KA, Elhanan G, Metcalf WJ, Slonim AD, Aweti R, et al. GWAS and PheWAS of red blood cell components in a Northern Nevadan cohort. Plos One. 2019;14(6):e0218078.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Cui JY, Zhou RR, Han S, Wang TS, Wang LQ, Xie XH. Statin therapy on glycemic control in type 2 diabetic patients: a network meta-analysis. J Clin Pharm Ther. 2018;43(4):556–70.

    Article  CAS  PubMed  Google Scholar 

  34. Ferrières J, Lautsch D, Gitt AK, De Ferrari G, Toplak H, Elisaf M, et al. Body mass index impacts the choice of lipid-lowering treatment with no correlation to blood cholesterol - findings from 52 916 patients in the Dyslipidemia International Study (DYSIS). Diab Obes Metab. 2018;20(11):2670–4.

    Article  CAS  Google Scholar 

  35. Mogadam M. Do statins increase the risk of diabetes or is it guilt by association? Curr Opin Endocrinol Diab Obes. 2014;21(2):140–5.

    Article  CAS  Google Scholar 

  36. Banaszewska B, Pawelczyk L, Spaczynski RZ, Duleba AJ. Comparison of simvastatin and metformin in treatment of polycystic ovary syndrome: prospective randomized trial. J Clin Endocrinol Metab. 2009;94(12):4938–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Dobs AS, Schrott H, Davidson MH, Bays H, Stein EA, Kush D, et al. Effects of high-dose simvastatin on adrenal and gonadal steroidogenesis in men with hypercholesterolemia. Metabolism. 2000;49(9):1234–8.

    Article  CAS  PubMed  Google Scholar 

  38. Larsson SC, Burgess S, Michaëlsson K. Association of genetic variants related to serum calcium levels with coronary artery disease and myocardial infarction. JAMA. 2017;318(4):371–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Xu L, Lin SL, Schooling CM. A Mendelian randomization study of the effect of calcium on coronary artery disease, myocardial infarction and their risk factors. Sci Rep. 2017;7(1):42691.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Schooling CM, Luo S, Au Yeung SL, Thompson DJ, Karthikeyan S, Bolton TR, et al. Genetic predictors of testosterone and their associations with cardiovascular disease and risk factors: a Mendelian randomization investigation. Int J Cardiol. 2018;267:171–6.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Dimou NL, Papadimitriou N, Gill D, Christakoudi S, Murphy N, Gunter MJ, et al. Sex hormone binding globulin and risk of breast cancer: a Mendelian randomization study. Int J Epidemiol. 2019;48(3):807–16.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Sacks FM, Pfeffer MA, Moye LA, Rouleau JL, Rutherford JD, Cole TG, et al. The effect of pravastatin on coronary events after myocardial infarction in patients with average cholesterol levels. Cholesterol and Recurrent Events Trial investigators. N Engl J Med. 1996;335(14):1001–9.

    Article  CAS  PubMed  Google Scholar 

  43. Cornish AJ, Law PJ, Timofeeva M, Palin K, Farrington SM, Palles C, et al. Modifiable pathways for colorectal cancer: a mendelian randomisation analysis. Lancet Gastroenterol Hepatol. 2020;5(1):55–62.

    Article  PubMed  Google Scholar 

  44. Keech A, Collins R, MacMahon S, Armitage J, Lawson A, Wallendszus K, et al. Three-year follow-up of the Oxford Cholesterol Study: assessment of the efficacy and safety of simvastin in preparation for a large mortality study. Eur Heart J. 1994;15:255–69.

  45. Vinholt PJ, Hvas AM, Frederiksen H, Bathum L, Jørgensen MK, Nybo M. Platelet count is associated with cardiovascular disease, cancer and mortality: a population-based cohort study. Thromb Res. 2016;148:136–42.

    Article  CAS  PubMed  Google Scholar 

  46. Gill D, Monori G, Georgakis MK, Tzoulaki I, Laffan M. Genetically determined platelet count and risk of cardiovascular disease. Arterioscler Thromb Vasc Biol. 2018;38(12):2862–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. He S, Lei W, Li J, Yu K, Yu Y, Zhou L, et al. Relation of platelet parameters with incident cardiovascular disease (the Dongfeng-Tongji Cohort Study). Am J Cardiol. 2019;123(2):239–48.

    Article  PubMed  Google Scholar 

  48. Pafili K, Penlioglou T, Mikhailidis DP, Papanas N. Mean platelet volume and coronary artery disease. Curr Opin Cardiol. 2019;34(4):390–8.

    Article  PubMed  Google Scholar 

  49. Harshfield EL, Sims MC, Traylor M, Ouwehand WH, Markus HS. The role of haematological traits in risk of ischaemic stroke and its subtypes. Brain. 2020;143(1):210–21.

    Article  PubMed  Google Scholar 

  50. Mach F, Ray KK, Wiklund O, Corsini A, Catapano AL, Bruckert E, et al. Adverse effects of statin therapy: perception vs. the evidence - focus on glucose homeostasis, cognitive, renal and hepatic function, haemorrhagic stroke and cataract. Eur Heart J. 2018;39(27):2526–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Schmidt AF, Swerdlow DI, Holmes MV, Patel RS, Fairhurst-Hunter Z, Lyall DM, et al. PCSK9 genetic variants and risk of type 2 diabetes: a Mendelian randomisation study. Lancet Diab Endocrinol. 2017;5(2):97–105.

    Article  CAS  Google Scholar 

  52. Page MM, Watts GF. PCSK9 inhibitors - mechanisms of action. Aust Prescr. 2016;39(5):164–7.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Handelsman Y, Lepor NE. PCSK9 inhibitors in lipid management of patients with diabetes mellitus and high cardiovascular risk: a review. J Am Heart Assoc. 2018;7(13):e008953.

    Article  CAS  Google Scholar 

  54. Phan BAP, Dayspring TD, Toth PP. Ezetimibe therapy: mechanism of action and clinical update. Vasc Health Risk Manag. 2012;8:415–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Stancu C, Sima A. Statins: mechanism of action and effects. J Cell Mol Med. 2001;5(4):378–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Johansson HE, Wåhlén A, Aldenbäck E, Haenni A. Platelet counts and liver enzymes after gastric bypass surgery. Obes Surg. 2018;28(6):1526–31.

    Article  PubMed  Google Scholar 

  57. Cole TJ, Henry CJ. The Oxford Brookes basal metabolic rate database--a reanalysis. Public Health Nutr. 2005;8(7a):1202–12.

    Article  CAS  PubMed  Google Scholar 

  58. Nikpay M, Goel A, Won H-H, Hall LM, Willenborg C, Kanoni S, et al. A comprehensive 1,000 Genomes-based genome-wide association meta-analysis of coronary artery disease. Nat Genet. 2015;47(10):1121–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Schooling CM, Au Yeung SL, Freeman G, Cowling BJ. The effect of statins on testosterone in men and women, a systematic review and meta-analysis of randomized controlled trials. BMC Med. 2013;11(1):57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Fulcher J, O'Connell R, Voysey M, Emberson J, Blackwell L, Mihaylova B, et al. Efficacy and safety of LDL-lowering therapy among men and women: meta-analysis of individual data from 174,000 participants in 27 randomised trials. Lancet. 2015;385(9976):1397–405.

    Article  Google Scholar 

  61. Hope HF, Binkley GM, Fenton S, Kitas GD, Verstappen SMM, Symmons DPM. Systematic review of the predictors of statin adherence for the primary prevention of cardiovascular disease. Plos One. 2019;14(1):e0201196.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Ference BA, Yoo W, Alesh I, Mahajan N, Mirowska KK, Mewada A, et al. Effect of long-term exposure to lower low-density lipoprotein cholesterol beginning early in life on the risk of coronary heart disease: a Mendelian randomization analysis. J Am Coll Cardiol. 2012;60(25):2631–9.

    Article  CAS  PubMed  Google Scholar 

  63. Bareinboim E, Pearl J. Causal inference and the data-fusion problem. Proc Natl Acad Sci USA. 2016;113(27):7345–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Batty GD, Gale CR, Kivimäki M, Deary IJ, Bell S. Comparison of risk factor associations in UK Biobank against representative, general population based studies with conventional response rates: prospective cohort study and individual participant meta-analysis. BMJ. 2020;368:m131.

    Article  Google Scholar 

  65. Schooling CM, Lopez PM, Yang Z, Zhao JV, Au Yeung SL, Huang JV. Use of multivariable Mendelian randomization to address biases due to competing risk before recruitment. Front Genet. 2021;11(1683):610852.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We thank Neale lab ( for providing publicly available sex-specific summary data from the UK Biobank and the University of Bristol for providing the MRC IEU OpenGWAS database for performing a PheWAS [19, 27]. We also thank Guoyi Yang for checking the study results.


Not applicable.

Author information

Authors and Affiliations



SL analyzed and interpreted the data and was a major contributor in writing the manuscript. CMS designed and supervised the study and revised the manuscript. The authors read and approved the final manuscript.

Corresponding author

Correspondence to C. M. Schooling.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Established genetic variants in the HMGCR, PCSK9 and NCP1L1 gene that mimic statins, PCSK9 inhibitors and ezetimibe respectively]. Table S2. Correlation matrix for genetic variants mimicking statins. Table S3. Correlation matrix for genetic variants mimicking PCSK9 inhibitors. Table S4. Correlation matrix for genetic variants mimicking ezetimibe. Table S5. 56 established independent SNPs genetically predicting LDL-c. Table S6. PheWAS table of rs12916. Table S7. Number of significant associations of rs12916 with phenotypes by category. Table S8. Sex-specific associations of genetically mimicked statins, PCSK9 inhibitors and ezetimibe (based on rs12916, rs11591147 and rs10260606, respectively) with body composition traits. Table S9. Sex-specific univariable MR of genetically predicted LDL-c on binary traits and body composition traits. Table S10. Replication table. Table S11. Sensitivity analysis for genetically mimicking statins. Table S12. Sensitivity analysis for genetically mimicking PCSK9 inhibitors. Table S13. Sensitivity analysis for genetically mimicking ezetimibe. Table S14. Univariable MR results of BMI on platelet count and BMR. Table S15. Multivariable MR results of genetically mimicked statins on BMR adjusting for BMI.

Additional file 2: Figure I.

Manhattan plots for rs12916 with all phenotypes included by category. The x-axis shows the phenotypes by category, and the y-axis shows the -log10 transformed P values. The blue line indicates the corrected statistical significance level, P = 2.9 × 10-5. The categories of phenotypes (from left to right on the x-axis) are blood biochemistry, circulatory, blood count, digestive, endocrine, neoplasms, respiratory, symptoms, physical measures, other categories (continuous) and other categories (binary). Other categories (continuous) include gynaecologic and obstetric, mental health, sense organs and others. Other categories (binary) include neurological, musculoskeletal, gynaecologic and obstetric, hematopoietic, dermatologic, genitourinary, mental health, infectious diseases, sense organs, injuries and poisonings and others. Figure II. Comparison between primary and sensitivity analysis for genetically mimicking statins. Figure III. Comparison between primary and sensitivity analysis for genetically mimicking PCSK9 inhibitors. Figure IV. Comparison between primary and sensitivity analysis for genetically mimicking ezetimibe.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, S., Schooling, C.M. A phenome-wide association study of genetically mimicked statins. BMC Med 19, 151 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: