Statins are the first-line lipid modifier for reducing cardiovascular morbidity and mortality (Ray et al., 2019; Michos et al., 2019). Statins have revolutionized the prevention and treatment of cardiovascular disease, and inspired the development of a range of effective interventions targeting the reduction of low-density lipoprotein (LDL)-cholesterol. Statins have long been suspected of having additional beneficial effects beyond lipid modulation (Schonbeck and Libby, 2004), such as on inflammation (Schonbeck and Libby, 2004), another potential target for reducing cardiovascular disease (Aday and Ridker, 2018). Meta-analysis of randomized controlled trials (RCTs) suggests statins are more effective at reducing mortality than proprotein convertase subtilisin/kexin type 9 (PCSK9) inhibitors or ezetimibe (Schmidt et al., 2017; Khan et al., 2018). However, these findings may be more apparent than real, stemming from differences in trial design, such as shorter duration of the PCSK9 inhibitor trials (Khan et al., 2018), the predominance of industry funded statin trials (Hobbs et al., 2016) or the difficulty of interpreting trials of ‘soft’ events when the treatment affects diagnostically relevant criteria, that is lipid levels (Schooling and Zhao, 2019). To investigate this anomaly, a previous study conducted a systematic agnostic scan of metabolic profile in a trial of statins compared to a PCSK9 inhibitor, which found few differences (Sliz et al., 2018). While characterization of the metabolic effects of statins suggested extensive effects on lipids and fatty acids (Würtz et al., 2016); these investigations were not able to include a factor which has previously been proposed as contributing to statin’s effectiveness, that is effects on male hormones (Schooling et al., 2013), although questions have been raised as to whether statins are as effective in women as men (Plakogiannis and Arif, 2016). However, meta-analysis of the available trial evidence suggests similar relative benefits of LDL-cholesterol reduction by statins for men and women (Fulcher et al., 2015), although the trials mainly concern men (73.2%) which may preclude detection of important sex differences. Men are also at substantially higher risk than women (Ezzati et al., 2015) giving larger absolute benefits in men than women at the same reduction in relative risk.

RCTs are not usually designed or powered to test mediating mechanisms. In addition, trials of statins on cardiovascular disease outcomes designed to be sex-specific are lacking. To assess a potential pathway by which statins might additionally operate, we used Mendelian randomization (MR), an observational study design that avoids confounding by taking advantage of the random allocation of genetic material at conception (Smith and Ebrahim, 2003), here specifically genetic variants mimicking effects of lipid modifiers. This random allocation at conception also avoids selection bias as long as few deaths have occurred between randomization and recruitment due to exposure, outcome, or other causes, that is competing risk, of the outcome (Schooling et al., 2020). So, here we focused on ischemic heart disease (IHD) (Kesteloot and Decramer, 2008), using the UK Biobank (Collins, 2012) to investigate whether testosterone mediated any of the effects of statins, PCSK9 inhibitors or ezetimibe on IHD in men or women using univariable and multivariable MR. As a further test, given some anti-inflammatories have also been shown to have opposite effects on male hormones compared to statins, specifically the interleukin one receptor antagonist (IL-1Ra), anakinra (Ebrahimi et al., 2018), we assessed whether the genetic variants mimicking effects of anakinra or tocilizumab targeting the interleukin six receptor (IL-6r) had opposite patterns of effects on testosterone and IHD (Aday and Ridker, 2018; Swerdlow et al., 2012; Interleukin 1 Genetics Consortium, 2015) to statins. Figure 1 illustrates the possible additional effects of statins, anakinra or tocilizumab on IHD via male hormones in the context of the well-established benefits of statins, PCSK9 inhibitors and ezetimibe acting via LDL-cholesterol and of anti-inflammatories in IHD.

Figure 1

Download asset Open asset

The six single neucleotide polymorphisms (SNPs) mimicking effects of statin (rs12916, rs5909, rs10066707, rs17238484, rs2006760 and rs2303152 from HMGCR) (Ference et al., 2019) were all correlated (r² >0.13). In the main analysis we used only the lead SNP, rs12916. Of the 7 SNPs mimicking effects of PCSK9 inhibitors (rs11206510, rs2149041, rs7552841, rs10888897, rs2479394, rs2479409 and rs562556 from PCSK9) (Ference et al., 2019), the three independently (r² <0.05) and most strongly associated with LDL-cholesterol (rs11206510, rs2149041 and rs7552841) were used in the main analysis. Of the 5 SNPs mimicking effects of ezetimibe (rs10260606 (proxy of rs2073547, r² = 0.99), rs2300414, rs10234070, rs7791240, rs217386 from NCP1L1) (Ference et al., 2019), two SNPs (rs2300414 and rs10234070) were discarded because their F-statistic for LDL-cholesterol was <10. The remaining three SNPs were all correlated at r² >0.05. rs2073547 was used in the main analysis because it had the strongest association with LDL-cholesterol. Supplementary file 1a shows the associations with LDL-cholesterol by sex for the independent SNPs used to mimic effects of statins, PCSK9 inhibitors and ezetimibe. Supplementary file 1b shows the associations of the two SNPs mimicking effects of the anti-inflammatory, anakinra, (rs6743376 and rs1542176 (r² <0.001)) with IL-1Ra and the associations of the SNP (rs7529229), mimicking effects of tocilizumab use, with IL-6r.

There were 179,918 men with 25,410 cases of IHD and 212,080 women with 12,511 cases of IHD in the UK Biobank.

Consistent with meta-analysis of RCTs this study provides genetic evidence that statins reduce testosterone in men (Schooling et al., 2013), and adds by showing that statins could partially operate on IHD, in men only, by reducing testosterone, as previously hypothesized (Schooling et al., 2013; Schooling et al., 2014) while having similar protective effects in men and women independent of testosterone. Previous Mendelian randomization studies have shown lower testosterone associated with lower risk of IHD, particularly in men (Schooling et al., 2018a; Luo et al., 2019; Mohammadi-Shemirani et al., 2019). Conversely, consistent with an RCT (Ebrahimi et al., 2018), this study also provides genetic validation that the anti-inflammatory anakinra, targeting IL-1Ra, increases testosterone in men, and is consistent with a previous Mendelian randomization study showing anakinra increases IHD (Interleukin 1 Genetics Consortium, 2015), but adds by showing why these associations might occur and that they may be specific to men.

A more marked association of testosterone with IHD in men than women (Table 2) is consistent with sex differences in biology, where testosterone is the main sex hormone in men and is much higher in men than in women. The associations of genetically mimicked effects of statins, PCSK9 inhibitors or ezetimibe on testosterone is consistent with the evidence available (Schooling et al., 2013; Ooi et al., 2015; Krysiak et al., 2015) and their mechanisms of action. Specifically, statins inhibit cholesterol synthesis, while PCSK9 inhibitors enable greater clearance of cholesterol, through increasing LDL-receptors, while ezetimibe reduces uptake of dietary cholesterol (Ray et al., 2019; Michos et al., 2019). However, some cells, such as Leydig cells, use de novo cholesterol synthesis to generate steroids, which can be reduced by statins (Shimizu-Albergine et al., 2016). Concerns about statins compromising androgen production pre-date the marketing of statins (Farnsworth et al., 1987; MacDonald et al., 1988). Beneficial immunosuppressive effects of androgens in rheumatoid arthritis have long been known (Cutolo et al., 1991), making androgen reduction a plausible mode of action for therapies, such as anakinra, whose primary indication is rheumatoid arthritis.

These findings may seem counter-intuitive given the essential role of testosterone in masculinity and reproduction. However, in 2015 the Food and Drug Administration in the United States required labelling changes for all testosterone prescriptions to warn of the risk of heart attacks and stroke on testosterone, although no sufficiently large RCT of testosterone administration has been conducted to confirm these effects (Onasanya et al., 2016). The Endocrine Society has also recommended caution in the use of testosterone (Bhasin et al., 2018). Meta-analysis of RCTs suggests androgen deprivation therapy reduces all-cause mortality, but is too small to quantify effects on specific diseases beyond prostate cancer (Nguyen et al., 2011). As such, our Mendelian randomization findings of the effects of testosterone have some consistency with the limited experimental evidence. In addition, our findings are consistent with well-established evolutionary biology theory, that is reproductive success may be at the expense of longevity, possibly in a sex-specific manner, impling that central drivers of the reproductive axis, as well as androgen production and catabolism, and their environmental cues may be relevant to IHD (Schooling, 2016; Schooling and Ng, 2019; Figure 2) encompassing the relations tested here (Figure 1). Notably, upregulation of indicators of plentiful living conditions, such as insulin, appear to cause IHD, particularly in men (Zhao et al., 2019), likely via gonadotropin releasing hormone (GnRH) (Schooling and Ng, 2019). Similarly, fatty acids may affect GnRH (Tran et al., 2016; Matsuyama and Kimura, 2015). In contrast, indicators of adversity, such as endotoxins promote an inflammatory response, involving interleukins, which suppresses the reproductive axis (Kalra et al., 1998) and thereby testosterone (Tremellen et al., 2018), which may be reversed by anakinra possibly outweighing the benefits for IHD of suppressing inflammation (Tardif et al., 2019). Statins, in contrast reduce androgen production, while agents that suppress androgen catabolism, such as rofecoxib, have also had unexpectedly adverse effects on IHD (Schooling, 2016). Mechanisms by which androgens might cause IHD have not been extensively investigated, but likely involve coagulation and red blood cell attributes. Several haemostatic and thrombotic factors, such as thromboxane A₂ (Pignatelli et al., 2012; Ajayi and Halushka, 2005), endothelin-1 (Polderman et al., 1993; van Kesteren et al., 1998; Sahebkar et al., 2015), nitric oxide (Pignatelli et al., 2012; Rosselli et al., 1998) and possibly thrombin (Orsi et al., 2019; Ferenchick et al., 1995), may be driven by androgens and likely play a role in IHD (Schooling et al., 2018b; Zhao, 2018; Nikpay et al., 2015). Von Willebrand factor (Sahebkar et al., 2016) and asymmetric dimethylarginine (Serban et al., 2015) are also modulated by statins and may cause IHD (Aday and Ridker, 2018; Au Yeung et al., 2016), whether they are driven by testosterone is unknown. Several red blood cell attributes are affected by androgens, from reticulocytes to hematocrit (Fernández-Balsells et al., 2010; Kanias et al., 2016), but exactly which causes IHD is unclear, although reticulocytes are a possibility (Astle et al., 2016). Currently, comprehensive genetic validation of these pathways is hampered by the lack of availability of large sex-specific genome wide association studies (GWAS) of cytokines and coagulation factors.

Figure 2

Download asset Open asset

Despite providing information that may be relevant to the performance of statins, and to the development of other therapies to protect against cardiovascular disease (Schooling, 2017), some limitations of this study exist. First, valid instruments should fulfill three assumptions, that is relate strongly to the exposure, not be associated with potential confounders and satisfy the exclusion restriction assumption. The F-statistics were >10. Despite high conditional F-statistics for testosterone the conditional F-statistics for the genetic mimics of statins were quite low and the Q-statistics for instrument validity were high suggesting pleiotropy, which we addressed by using multivariable MR-Egger. The associations with testosterone in women were not adjusted for factors, such as menopausal status, hormone use and history of oophorectomy, which could result in imprecision and weaker instruments. The SNPs used to mimic effects of statins, PCSK9 inhibitors and ezetimibe are well established (Ference et al., 2019), and in genes that harbor the target of each lipid modifier (HMGCR, PCSK9 and NCP1L1 respectively). We did not include body mass index (BMI) as a risk factor explaining the effect of statins on IHD, because statins increase BMI (Swerdlow et al., 2015) and decrease the risk of IHD, so including them in the multivariable analysis may inflate the effect of mimicking statins on IHD, rather than explaining part of their effect on IHD. The SNPs mimicking effects of the anti-inflammatory anakinra have been validated as increasing IL-1Ra (Interleukin 1 Genetics Consortium, 2015), and the SNP used to mimic effects of tocilizumab is well-established as affecting IL-6r (Swerdlow et al., 2012). Testosterone’s effects on IHD in men could be via adiposity, insulin or LDL-cholesterol rather than via testosterone. However, consistent with a previous MR study, we found testosterone did not affect BMI in men (Eriksson et al., 2017), we also found little evidence that testosterone in men affected LDL-cholesterol (data not shown). We could not test whether testosterone in men affects insulin because of the lack of an insulin GWAS including the X chromosome. Sex-specific genetic associations were used throughout with exception of the genetic mimics of effects of anakinra and tocilizumab on IL1Ra and IL-6r respectively. However, inflammation operating on the reproductive axis would be expected to have sex-specific effects not sex-specific drivers. We selected between correlated SNPs based on p-values which is relatively arbitrary, and the estimates could be sensitive to the choice of SNPs. Repeating the analysis using a larger number of correlated SNPs, where possible, taking into account their correlation, gave a similar interpretation. MR studies can be confounded by population stratification. However, we used genetic associations from GWAS mainly comprising people of white British ancestry with genomic control. Functions of each SNP predicting the exposures are not all fully understood, so we cannot rule out the possibility that the SNPs are linked with IHD through other pathways although we used sensitivity analysis.

We used SNPs predicting testosterone, but not other exposures, obtained from the same study as the genetic effects on IHD. However given the estimates for testosterone were largely obtained from non-cases, the overlap unlikely introduced substantial bias (Burgess et al., 2016). Canalization, that is buffering of genetic factors during development, may occur however; whether it does so is unknown. Our findings, largely in Europeans, may not be applicable to other populations. However, causes are unlikely to act differently in different populations, although the causal mechanisms may not be as relevant in all settings (Lopez et al., 2019). The SNPs mimicking effects of statins, PCSK9 inhibitors and ezetimibe were previously selected for their relations with LDL-cholesterol and on functional grounds (Ference et al., 2019), assuming the lipid modifiers act by action on lipids (Ference et al., 2019), so it is possible that relevant SNPs might have been discarded if they work through other mechanisms independent of lipids. It is also possible that the SNPs mimicking lipid modifiers might act via a different lipid trait, such as apoB (Richardson et al., 2020). The SNPs mimicking effects of anakinra and tocilizumab were similarly selected. Replication based on genetic instruments functionally relevant to all the exposures would be ideal. However, we used the most recent, published genetic instruments for testosterone (Ruth et al., 2020). Replication based on another large sex-specific IHD GWAS where the IHD cases are not from the same study as the testosterone instruments, would be ideal. However, sex-specific summary statistics are not available for large existing IHD GWAS, such as CARDIoGRAM (Nikpay et al., 2015). Moreover, the UK Biobank has the advantage of being very intensively genotyped and including the X chromosome, which is important for testosterone, but is not usually included in publicly available summary statistics. Lack of replication is a limitation of this study. Lastly, Mendelian randomization assesses the lifelong effects of an endogenous exposure rather than short-term effects of an interventions assessed in an RCT. Our estimates give an indication of the role of the exposures rather than the exact effects of the corresponding interventions. Nevertheless our estimates for statins on IHD are comparable with meta-analyses of statin trials considering similar outcomes (Fulcher et al., 2015).

Here, we present a hypothesis driven study examining the role of testosterone in mediating the effect of specifically statins in IHD, particularly in men. Future work could encompass a comprehensive sex-specific multivariable MR to confirm the role of sex hormones and sex hormone binding globulin in IHD as well as any mediation of their effects by key lipids, such as LDL-cholesterol or apoB. This work would be facilitated by the development of published genetic instruments for estrogen in women. Future work could also encompass assessing whether any other drugs that reduce cardiovascular disease, such as canakinumab (Ridker et al., 2017), also impact testosterone.

Taken together these complimentary findings for statins and anakinra raise the possibility that modulating testosterone, by whatever means, is a relevant feature for modulating IHD in men, with potential relevance to the development of new interventions, side-effects of existing interventions, re-purposing and appropriate use. Statins lowering testosterone could also be relevant to the muscle weakness or pain experienced on statins (Collins et al., 2016). Recognition that statins lower testosterone might also provide greater impetus for investigation of their role in other relevant conditions, such as prostate cancer (Alfaqih et al., 2017). Conversely, statins and anakinra did not clearly affect testosterone in women (Table 1) nor did testosterone mediate the effect of statins on IHD in women (Table 3). These differences by sex highlight the need for sex-specific approaches to IHD prevention and management, specifically in terms of the use of statins and investigation more broadly of causes of IHD.

All data generated or analysed during this study are included in the manuscript and supporting files.

1. 1. Aday AW
   2. Ridker PM

   (2018) Antiinflammatory therapy in clinical care: the CANTOS trial and beyond

   Frontiers in Cardiovascular Medicine 5:62.

   https://doi.org/10.3389/fcvm.2018.00062

   • PubMed
   • Google Scholar
2. 1. Ajayi AAL
   2. Halushka PV

   (2005) Castration reduces platelet thromboxane A2 receptor density and aggregability

   QJM: An International Journal of Medicine 98:349–356.

   https://doi.org/10.1093/qjmed/hci054

   • Google Scholar
3. 1. Alfaqih MA
   2. Allott EH
   3. Hamilton RJ
   4. Freeman MR
   5. Freedland SJ

   (2017) The current evidence on statin use and prostate Cancer prevention: are we there yet?

   Nature Reviews Urology 14:107–119.

   https://doi.org/10.1038/nrurol.2016.199

   • PubMed
   • Google Scholar
4. 1. Astle WJ
   2. Elding H
   3. Jiang T
   4. Allen D
   5. Ruklisa D
   6. Mann AL
   7. Mead D
   8. Bouman H
   9. Riveros-Mckay F
   10. Kostadima MA
   11. Lambourne JJ
   12. Sivapalaratnam S
   13. Downes K
   14. Kundu K
   15. Bomba L
   16. Berentsen K
   17. Bradley JR
   18. Daugherty LC
   19. Delaneau O
   20. Freson K
   21. Garner SF
   22. Grassi L
   23. Guerrero J
   24. Haimel M
   25. Janssen-Megens EM
   26. Kaan A
   27. Kamat M
   28. Kim B
   29. Mandoli A
   30. Marchini J
   31. Martens JHA
   32. Meacham S
   33. Megy K
   34. O’Connell J
   35. Petersen R
   36. Sharifi N
   37. Sheard SM
   38. Staley JR
   39. Tuna S
   40. van der Ent M
   41. Walter K
   42. Wang S-Y
   43. Wheeler E
   44. Wilder SP
   45. Iotchkova V
   46. Moore C
   47. Sambrook J
   48. Stunnenberg HG
   49. Di Angelantonio E
   50. Kaptoge S
   51. Kuijpers TW
   52. Carrillo-de-Santa-Pau E
   53. Juan D
   54. Rico D
   55. Valencia A
   56. Chen L
   57. Ge B
   58. Vasquez L
   59. Kwan T
   60. Garrido-Martín D
   61. Watt S
   62. Yang Y
   63. Guigo R
   64. Beck S
   65. Paul DS
   66. Pastinen T
   67. Bujold D
   68. Bourque G
   69. Frontini M
   70. Danesh J
   71. Roberts DJ
   72. Ouwehand WH
   73. Butterworth AS
   74. Soranzo N

   (2016) The allelic landscape of human blood cell trait variation and links to common complex disease

   Cell 167:1415–1429.

   https://doi.org/10.1016/j.cell.2016.10.042

   • Google Scholar
5. 1. Au Yeung SL
   2. Lin SL
   3. Lam HS
   4. Schooling CM

   (2016) Effect of l-arginine, asymmetric dimethylarginine, and symmetric dimethylarginine on ischemic heart disease risk: a mendelian randomization study

   American Heart Journal 182:54–61.

   https://doi.org/10.1016/j.ahj.2016.07.021

   • PubMed
   • Google Scholar
6. 1. Bhasin S
   2. Brito JP
   3. Cunningham GR
   4. Hayes FJ
   5. Hodis HN
   6. Matsumoto AM
   7. Snyder PJ
   8. Swerdloff RS
   9. Wu FC
   10. Yialamas MA

   (2018) Testosterone therapy in men with hypogonadism: an endocrine society clinical practice guideline

   The Journal of Clinical Endocrinology & Metabolism 103:1715–1744.

   https://doi.org/10.1210/jc.2018-00229

   • PubMed
   • Google Scholar
7. 1. Bowden J
   2. Davey Smith G
   3. Burgess S

   (2015) Mendelian randomization with invalid instruments: effect estimation and Bias detection through egger regression

   International Journal of Epidemiology 44:512–525.

   https://doi.org/10.1093/ije/dyv080

   • PubMed
   • Google Scholar
8. 1. Bowden J
   2. Del Greco M F
   3. Minelli C
   4. Davey Smith G
   5. Sheehan NA
   6. Thompson JR

   (2016a) Assessing the suitability of summary data for two-sample mendelian randomization analyses using MR-Egger regression: the role of the I2 statistic

   International Journal of Epidemiology 45:1961–1974.

   https://doi.org/10.1093/ije/dyw220

   • PubMed
   • Google Scholar
9. 1. Bowden J
   2. Davey Smith G
   3. Haycock PC
   4. Burgess S

   (2016b) Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator

   Genetic Epidemiology 40:304–314.

   https://doi.org/10.1002/gepi.21965

   • PubMed
   • Google Scholar
10. 1. Burgess S
    2. Butterworth A
    3. Thompson SG

    (2013) Mendelian randomization analysis with multiple genetic variants using summarized data

    Genetic Epidemiology 37:658–665.

    https://doi.org/10.1002/gepi.21758

    • PubMed
    • Google Scholar
11. 1. Burgess S
    2. Davies NM
    3. Thompson SG

    (2016) Bias due to participant overlap in two-sample mendelian randomization

    Genetic Epidemiology 40:597–608.

    https://doi.org/10.1002/gepi.21998

    • PubMed
    • Google Scholar
12. 1. Burgess S
    2. Thompson SG

    (2017) Interpreting findings from mendelian randomization using the MR-Egger method

    European Journal of Epidemiology 32:377–389.

    https://doi.org/10.1007/s10654-017-0255-x

    • PubMed
    • Google Scholar
13. 1. Collins R

    (2012) What makes UK biobank special?

    The Lancet 379:1173–1174.

    https://doi.org/10.1016/S0140-6736(12)60404-8

    • Google Scholar
14. 1. Collins R
    2. Reith C
    3. Emberson J
    4. Armitage J
    5. Baigent C
    6. Blackwell L
    7. Blumenthal R
    8. Danesh J
    9. Smith GD
    10. DeMets D
    11. Evans S
    12. Law M
    13. MacMahon S
    14. Martin S
    15. Neal B
    16. Poulter N
    17. Preiss D
    18. Ridker P
    19. Roberts I
    20. Rodgers A
    21. Sandercock P
    22. Schulz K
    23. Sever P
    24. Simes J
    25. Smeeth L
    26. Wald N
    27. Yusuf S
    28. Peto R

    (2016) Interpretation of the evidence for the efficacy and safety of statin therapy

    The Lancet 388:2532–2561.

    https://doi.org/10.1016/S0140-6736(16)31357-5

    • Google Scholar
15. 1. Cutolo M
    2. Balleari E
    3. Giusti M
    4. Intra E
    5. Accardo S

    (1991) Androgen replacement therapy in male patients with rheumatoid arthritis

    Arthritis & Rheumatism 34:1–5.

    https://doi.org/10.1002/art.1780340102

    • Google Scholar
16. 1. Ebrahimi F
    2. Urwyler SA
    3. Straumann S
    4. Doerpfeld S
    5. Bernasconi L
    6. Neyer P
    7. Schuetz P
    8. Mueller B
    9. Donath MY
    10. Christ-Crain M

    (2018) IL-1 antagonism in men with metabolic syndrome and low testosterone: a randomized clinical trial

    The Journal of Clinical Endocrinology & Metabolism 103:3466–3476.

    https://doi.org/10.1210/jc.2018-00739

    • PubMed
    • Google Scholar
17. 1. Eriksson J
    2. Haring R
    3. Grarup N
    4. Vandenput L
    5. Wallaschofski H
    6. Lorentzen E
    7. Hansen T
    8. Mellström D
    9. Pedersen O
    10. Nauck M
    11. Lorentzon M
    12. Nystrup Husemoen LL
    13. Völzke H
    14. Karlsson M
    15. Baumeister SE
    16. Linneberg A
    17. Ohlsson C

    (2017) Causal relationship between obesity and serum testosterone status in men: a bi-directional mendelian randomization analysis

    PLOS ONE 12:e0176277.

    https://doi.org/10.1371/journal.pone.0176277

    • PubMed
    • Google Scholar
18. 1. Ezzati M
    2. Obermeyer Z
    3. Tzoulaki I
    4. Mayosi BM
    5. Elliott P
    6. Leon DA

    (2015) Contributions of risk factors and medical care to cardiovascular mortality trends

    Nature Reviews Cardiology 12:508–530.

    https://doi.org/10.1038/nrcardio.2015.82

    • PubMed
    • Google Scholar
19. 1. Farnsworth WH
    2. Hoeg JM
    3. Maher M
    4. Brittain EH
    5. Sherins RJ
    6. Brewer HB

    (1987) Testicular function in type II hyperlipoproteinemic patients treated with lovastatin (mevinolin) or neomycin

    The Journal of Clinical Endocrinology & Metabolism 65:546–550.

    https://doi.org/10.1210/jcem-65-3-546

    • PubMed
    • Google Scholar
20. 1. Ference BA
    2. Ray KK
    3. Catapano AL
    4. Ference TB
    5. Burgess S
    6. Neff DR
    7. Oliver-Williams C
    8. Wood AM
    9. Butterworth AS
    10. Di Angelantonio E
    11. Danesh J
    12. Kastelein JJP
    13. Nicholls SJ

    (2019) Mendelian randomization study of ACLY and Cardiovascular Disease

    New England Journal of Medicine 380:1033–1042.

    https://doi.org/10.1056/NEJMoa1806747

    • PubMed
    • Google Scholar
21. 1. Ferenchick GS
    2. Hirokawa S
    3. Mammen EF
    4. Schwartz KA

    (1995) Anabolic-androgenic steroid abuse in weight lifters: evidence for activation of the hemostatic system

    American Journal of Hematology 49:282–288.

    https://doi.org/10.1002/ajh.2830490405

    • PubMed
    • Google Scholar
22. 1. Fernández-Balsells MM
    2. Murad MH
    3. Lane M
    4. Lampropulos JF
    5. Albuquerque F
    6. Mullan RJ
    7. Agrwal N
    8. Elamin MB
    9. Gallegos-Orozco JF
    10. Wang AT
    11. Erwin PJ
    12. Bhasin S
    13. Montori VM

    (2010) Clinical review 1: adverse effects of testosterone therapy in adult men: a systematic review and meta-analysis

    The Journal of Clinical Endocrinology & Metabolism 95:2560–2575.

    https://doi.org/10.1210/jc.2009-2575

    • PubMed
    • Google Scholar
23. 1. Fulcher J
    2. O'Connell R
    3. Voysey M
    4. Emberson J
    5. Blackwell L
    6. Mihaylova B
    7. Simes J
    8. Collins R
    9. Kirby A
    10. Colhoun H
    11. Braunwald E
    12. La Rosa J
    13. Pedersen TR
    14. Tonkin A
    15. Davis B
    16. Sleight P
    17. Franzosi MG
    18. Baigent C
    19. Keech A
    20. Cholesterol Treatment Trialists' (CTT) Collaboration

    (2015) 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 385:1397–1405.

    https://doi.org/10.1016/S0140-6736(14)61368-4

    • PubMed
    • Google Scholar
24. 1. Haring R
    2. Hannemann A
    3. John U
    4. Radke D
    5. Nauck M
    6. Wallaschofski H
    7. Owen L
    8. Adaway J
    9. Keevil BG
    10. Brabant G

    (2012) Age-specific reference ranges for serum testosterone and androstenedione concentrations in women measured by liquid chromatography-tandem mass spectrometry

    The Journal of Clinical Endocrinology & Metabolism 97:408–415.

    https://doi.org/10.1210/jc.2011-2134

    • PubMed
    • Google Scholar
25. 1. Hemani G
    2. Tilling K
    3. Davey Smith G

    (2017) Orienting the causal relationship between imprecisely measured traits using GWAS summary data

    PLOS Genetics 13:e1007081.

    https://doi.org/10.1371/journal.pgen.1007081

    • PubMed
    • Google Scholar
26. 1. Hobbs FD
    2. Banach M
    3. Mikhailidis DP
    4. Malhotra A
    5. Capewell S

    (2016) Is statin-modified reduction in lipids the most important preventive therapy for cardiovascular disease? A pro/con debate

    BMC Medicine 14:4.

    https://doi.org/10.1186/s12916-016-0550-5

    • PubMed
    • Google Scholar
27. 1. Interleukin 1 Genetics Consortium

    (2015) Cardiometabolic effects of genetic upregulation of the interleukin 1 receptor antagonist: a mendelian randomisation analysis

    The Lancet Diabetes & Endocrinology 3:243–253.

    https://doi.org/10.1016/S2213-8587(15)00034-0

    • PubMed
    • Google Scholar
28. 1. Kalra PS
    2. Edwards TG
    3. Xu B
    4. Jain M
    5. Kalra SP

    (1998) The anti-gonadotropic effects of cytokines: the role of neuropeptides

    Domestic Animal Endocrinology 15:321–332.

    https://doi.org/10.1016/S0739-7240(98)00030-7

    • PubMed
    • Google Scholar
29. 1. Kanias T
    2. Sinchar D
    3. Osei-Hwedieh D
    4. Baust JJ
    5. Jordan A
    6. Zimring JC
    7. Waterman HR
    8. de Wolski KS
    9. Acker JP
    10. Gladwin MT

    (2016) Testosterone-dependent sex differences in red blood cell hemolysis in storage, stress, and disease

    Transfusion 56:2571–2583.

    https://doi.org/10.1111/trf.13745

    • Google Scholar
30. 1. Kesteloot H
    2. Decramer M

    (2008) Age at death from different diseases: the flemish experience during the period 2000-2004

    Acta Clinica Belgica 63:256–261.

    https://doi.org/10.1179/acb.2008.047

    • PubMed
    • Google Scholar
31. 1. Khan SU
    2. Talluri S
    3. Riaz H
    4. Rahman H
    5. Nasir F
    6. Bin Riaz I
    7. Sattur S
    8. Ahmed H
    9. Kaluski E
    10. Krasuski R

    (2018) A bayesian network meta-analysis of PCSK9 inhibitors, statins and ezetimibe with or without statins for cardiovascular outcomes

    European Journal of Preventive Cardiology 25:844–853.

    https://doi.org/10.1177/2047487318766612

    • PubMed
    • Google Scholar
32. 1. Krysiak R
    2. Kowalska B
    3. Żmuda W
    4. Okopień B

    (2015) The effect of ezetimibe-statin combination on steroid hormone production in men with coronary artery disease and low cholesterol levels

    Pharmacological Reports 67:305–309.

    https://doi.org/10.1016/j.pharep.2014.10.004

    • PubMed
    • Google Scholar
33. 1. Lopez PM
    2. Subramanian SV
    3. Schooling CM

    (2019) Effect measure modification conceptualized using selection diagrams as mediation by mechanisms of varying population-level relevance

    Journal of Clinical Epidemiology 113:123–128.

    https://doi.org/10.1016/j.jclinepi.2019.05.005

    • PubMed
    • Google Scholar
34. 1. Luo S
    2. Au Yeung SL
    3. Zhao JV
    4. Burgess S
    5. Schooling CM

    (2019) Association of genetically predicted testosterone with thromboembolism, heart failure, and myocardial infarction: mendelian randomisation study in UK biobank

    BMJ 364:l476.

    https://doi.org/10.1136/bmj.l476

    • PubMed
    • Google Scholar
35. 1. MacDonald JS
    2. Gerson RJ
    3. Kornbrust DJ
    4. Kloss MW
    5. Prahalada S
    6. Berry PH
    7. Alberts AW
    8. Bokelman DL

    (1988) Preclinical evaluation of lovastatin

    The American Journal of Cardiology 62:J16–J27.

    https://doi.org/10.1016/0002-9149(88)90003-3

    • PubMed
    • Google Scholar
36. 1. Matsuyama S
    2. Kimura K

    (2015) Regulation of gonadotropin secretion by monitoring energy availability

    Reproductive Medicine and Biology 14:39–47.

    https://doi.org/10.1007/s12522-014-0194-0

    • PubMed
    • Google Scholar
37. 1. Michos ED
    2. McEvoy JW
    3. Blumenthal RS

    (2019) Lipid management for the prevention of atherosclerotic cardiovascular disease

    New England Journal of Medicine 381:1557–1567.

    https://doi.org/10.1056/NEJMra1806939

    • PubMed
    • Google Scholar
38. Preprint

    1. Mohammadi-Shemirani PC
    2. Pigeyre M
    3. Morton M
    4. Gerstein RW

    (2019) Clinical benefits and adverse effects of genetically-elevated free testosterone levels: a mendelian randomization analysis

    medRxiv.

    https://doi.org/10.1101/19005132

    • Google Scholar
39. 1. Mohr BA
    2. Guay AT
    3. O'Donnell AB
    4. McKinlay JB

    (2005) Normal, bound and nonbound testosterone levels in normally ageing men: results from the massachusetts male ageing study

    Clinical Endocrinology 62:64–73.

    https://doi.org/10.1111/j.1365-2265.2004.02174.x

    • PubMed
    • Google Scholar
40. 1. Nguyen PL
    2. Je Y
    3. Schutz FA
    4. Hoffman KE
    5. Hu JC
    6. Parekh A
    7. Beckman JA
    8. Choueiri TK

    (2011) Association of androgen deprivation therapy with cardiovascular death in patients with prostate Cancer: a meta-analysis of randomized trials

    Jama 306:2359–2366.

    https://doi.org/10.1001/jama.2011.1745

    • PubMed
    • Google Scholar
41. 1. Nikpay M
    2. Goel A
    3. Won HH
    4. Hall LM
    5. Willenborg C
    6. Kanoni S
    7. Saleheen D
    8. Kyriakou T
    9. Nelson CP
    10. Hopewell JC
    11. Webb TR
    12. Zeng L
    13. Dehghan A
    14. Alver M
    15. Armasu SM
    16. Auro K
    17. Bjonnes A
    18. Chasman DI
    19. Chen S
    20. Ford I
    21. Franceschini N
    22. Gieger C
    23. Grace C
    24. Gustafsson S
    25. Huang J
    26. Hwang SJ
    27. Kim YK
    28. Kleber ME
    29. Lau KW
    30. Lu X
    31. Lu Y
    32. Lyytikäinen LP
    33. Mihailov E
    34. Morrison AC
    35. Pervjakova N
    36. Qu L
    37. Rose LM
    38. Salfati E
    39. Saxena R
    40. Scholz M
    41. Smith AV
    42. Tikkanen E
    43. Uitterlinden A
    44. Yang X
    45. Zhang W
    46. Zhao W
    47. de Andrade M
    48. de Vries PS
    49. van Zuydam NR
    50. Anand SS
    51. Bertram L
    52. Beutner F
    53. Dedoussis G
    54. Frossard P
    55. Gauguier D
    56. Goodall AH
    57. Gottesman O
    58. Haber M
    59. Han BG
    60. Huang J
    61. Jalilzadeh S
    62. Kessler T
    63. König IR
    64. Lannfelt L
    65. Lieb W
    66. Lind L
    67. Lindgren CM
    68. Lokki ML
    69. Magnusson PK
    70. Mallick NH
    71. Mehra N
    72. Meitinger T
    73. Memon FU
    74. Morris AP
    75. Nieminen MS
    76. Pedersen NL
    77. Peters A
    78. Rallidis LS
    79. Rasheed A
    80. Samuel M
    81. Shah SH
    82. Sinisalo J
    83. Stirrups KE
    84. Trompet S
    85. Wang L
    86. Zaman KS
    87. Ardissino D
    88. Boerwinkle E
    89. Borecki IB
    90. Bottinger EP
    91. Buring JE
    92. Chambers JC
    93. Collins R
    94. Cupples LA
    95. Danesh J
    96. Demuth I
    97. Elosua R
    98. Epstein SE
    99. Esko T
    100. Feitosa MF
    101. Franco OH
    102. Franzosi MG
    103. Granger CB
    104. Gu D
    105. Gudnason V
    106. Hall AS
    107. Hamsten A
    108. Harris TB
    109. Hazen SL
    110. Hengstenberg C
    111. Hofman A
    112. Ingelsson E
    113. Iribarren C
    114. Jukema JW
    115. Karhunen PJ
    116. Kim BJ
    117. Kooner JS
    118. Kullo IJ
    119. Lehtimäki T
    120. Loos RJF
    121. Melander O
    122. Metspalu A
    123. März W
    124. Palmer CN
    125. Perola M
    126. Quertermous T
    127. Rader DJ
    128. Ridker PM
    129. Ripatti S
    130. Roberts R
    131. Salomaa V
    132. Sanghera DK
    133. Schwartz SM
    134. Seedorf U
    135. Stewart AF
    136. Stott DJ
    137. Thiery J
    138. Zalloua PA
    139. O'Donnell CJ
    140. Reilly MP
    141. Assimes TL
    142. Thompson JR
    143. Erdmann J
    144. Clarke R
    145. Watkins H
    146. Kathiresan S
    147. McPherson R
    148. Deloukas P
    149. Schunkert H
    150. Samani NJ
    151. Farrall M

    (2015) A comprehensive 1,000 Genomes-based genome-wide association meta-analysis of coronary artery disease

    Nature Genetics 47:1121–1130.

    https://doi.org/10.1038/ng.3396

    • PubMed
    • Google Scholar
42. 1. Onasanya O
    2. Iyer G
    3. Lucas E
    4. Lin D
    5. Singh S
    6. Alexander GC

    (2016) Association between exogenous testosterone and cardiovascular events: an overview of systematic reviews

    The Lancet Diabetes & Endocrinology 4:943–956.

    https://doi.org/10.1016/S2213-8587(16)30215-7

    • PubMed
    • Google Scholar
43. 1. Ooi TC
    2. Raymond A
    3. Cousins M
    4. Favreau C
    5. Taljaard M
    6. Gavin C
    7. Jolly EE
    8. Malone S
    9. Eapen L
    10. Chretien M
    11. Mbikay M
    12. Mayne J

    (2015) Relationship between testosterone, estradiol and circulating PCSK9: cross-sectional and interventional studies in humans

    Clinica Chimica Acta 446:97–104.

    https://doi.org/10.1016/j.cca.2015.03.036

    • Google Scholar
44. 1. Orsi FA
    2. Biedermann JS
    3. Kruip M
    4. van der Meer FJ
    5. Rosendaal FR
    6. van Hylckama Vlieg A
    7. Bos MHA
    8. Leebeek FWG
    9. Cannegieter SC
    10. Lijfering WM

    (2019) Rosuvastatin use reduces thrombin generation potential in patients with venous thromboembolism: a randomized controlled trial

    Journal of Thrombosis and Haemostasis 17:319–328.

    https://doi.org/10.1111/jth.14364

    • PubMed
    • Google Scholar
45. 1. Pignatelli P
    2. Carnevale R
    3. Pastori D
    4. Cangemi R
    5. Napoleone L
    6. Bartimoccia S
    7. Nocella C
    8. Basili S
    9. Violi F

    (2012) Immediate antioxidant and antiplatelet effect of atorvastatin via inhibition of Nox2

    Circulation 126:92–103.

    https://doi.org/10.1161/CIRCULATIONAHA.112.095554

    • PubMed
    • Google Scholar
46. 1. Plakogiannis R
    2. Arif SA

    (2016) Women versus men: is there equal benefit and safety from statins?

    Current Atherosclerosis Reports 18:6.

    https://doi.org/10.1007/s11883-016-0562-9

    • PubMed
    • Google Scholar
47. 1. Polderman KH
    2. Stehouwer CD
    3. van Kamp GJ
    4. Dekker GA
    5. Verheugt FW
    6. Gooren LJ

    (1993) Influence of sex hormones on plasma endothelin levels

    Annals of Internal Medicine 118:429–432.

    https://doi.org/10.7326/0003-4819-118-6-199303150-00006

    • PubMed
    • Google Scholar
48. 1. Rafiq S
    2. Frayling TM
    3. Murray A
    4. Hurst A
    5. Stevens K
    6. Weedon MN
    7. Henley W
    8. Ferrucci L
    9. Bandinelli S
    10. Corsi AM
    11. Guralnik JM
    12. Melzer D

    (2007) A common variant of the interleukin 6 receptor (IL-6r) gene increases IL-6r and IL-6 levels, without other inflammatory effects

    Genes & Immunity 8:552–559.

    https://doi.org/10.1038/sj.gene.6364414

    • PubMed
    • Google Scholar
49. 1. Ray KK
    2. Corral P
    3. Morales E
    4. Nicholls SJ

    (2019) Pharmacological lipid-modification therapies for prevention of ischaemic heart disease: current and future options

    The Lancet 394:697–708.

    https://doi.org/10.1016/S0140-6736(19)31950-6

    • Google Scholar
50. 1. Rees JMB
    2. Wood AM
    3. Burgess S

    (2017) Extending the MR-Egger method for multivariable mendelian randomization to correct for both measured and unmeasured pleiotropy

    Statistics in Medicine 36:4705–4718.

    https://doi.org/10.1002/sim.7492

    • PubMed
    • Google Scholar
51. 1. Richardson TG
    2. Sanderson E
    3. Palmer TM
    4. Ala-Korpela M
    5. Ference BA
    6. Davey Smith G
    7. Holmes MV

    (2020) Evaluating the relationship between circulating lipoprotein lipids and apolipoproteins with risk of coronary heart disease: a multivariable mendelian randomisation analysis

    PLOS Medicine 17:e1003062.

    https://doi.org/10.1371/journal.pmed.1003062

    • PubMed
    • Google Scholar
52. 1. Ridker PM
    2. Everett BM
    3. Thuren T
    4. MacFadyen JG
    5. Chang WH
    6. Ballantyne C
    7. Fonseca F
    8. Nicolau J
    9. Koenig W
    10. Anker SD
    11. Kastelein JJP
    12. Cornel JH
    13. Pais P
    14. Pella D
    15. Genest J
    16. Cifkova R
    17. Lorenzatti A
    18. Forster T
    19. Kobalava Z
    20. Vida-Simiti L
    21. Flather M
    22. Shimokawa H
    23. Ogawa H
    24. Dellborg M
    25. Rossi PRF
    26. Troquay RPT
    27. Libby P
    28. Glynn RJ
    29. CANTOS Trial Group

    (2017) Antiinflammatory therapy with canakinumab for atherosclerotic disease

    New England Journal of Medicine 377:1119–1131.

    https://doi.org/10.1056/NEJMoa1707914

    • PubMed
    • Google Scholar
53. 1. Rosselli M
    2. Keller PJ
    3. Dubey RK

    (1998) Role of nitric oxide in the biology, physiology and pathophysiology of reproduction

    Human Reproduction Update 4:3–24.

    https://doi.org/10.1093/humupd/4.1.3

    • PubMed
    • Google Scholar
54. 1. Ruth KS
    2. Day FR
    3. Tyrrell J
    4. Thompson DJ
    5. Wood AR
    6. Mahajan A
    7. Beaumont RN
    8. Wittemans L
    9. Martin S
    10. Busch AS
    11. Erzurumluoglu AM
    12. Hollis B
    13. O'Mara TA
    14. McCarthy MI
    15. Langenberg C
    16. Easton DF
    17. Wareham NJ
    18. Burgess S
    19. Murray A
    20. Ong KK
    21. Frayling TM
    22. Perry JRB
    23. Endometrial Cancer Association Consortium

    (2020) Using human genetics to understand the disease impacts of testosterone in men and women

    Nature Medicine 26:252–258.

    https://doi.org/10.1038/s41591-020-0751-5

    • PubMed
    • Google Scholar
55. 1. Sahebkar A
    2. Kotani K
    3. Serban C
    4. Ursoniu S
    5. Mikhailidis DP
    6. Jones SR
    7. Ray KK
    8. Blaha MJ
    9. Rysz J
    10. Toth PP
    11. Muntner P
    12. Lip GY
    13. Banach M
    14. Lipid and Blood Pressure Meta-analysis Collaboration (LBPMC) Group

    (2015) Statin therapy reduces plasma endothelin-1 concentrations: a meta-analysis of 15 randomized controlled trials

    Atherosclerosis 241:433–442.

    https://doi.org/10.1016/j.atherosclerosis.2015.05.022

    • PubMed
    • Google Scholar
56. 1. Sahebkar A
    2. Serban C
    3. Ursoniu S
    4. Mikhailidis DP
    5. Undas A
    6. Lip GY
    7. Bittner V
    8. Ray K
    9. Watts GF
    10. Hovingh GK
    11. Rysz J
    12. Kastelein JJ
    13. Banach M
    14. Lipid and Blood Pressure Meta-analysis Collaboration (LBPMC) Group

    (2016) The impact of statin therapy on plasma levels of von willebrand factor antigen. systematic review and meta-analysis of randomised placebo-controlled trials

    Thrombosis and Haemostasis 115:520–532.

    https://doi.org/10.1160/TH15-08-0620

    • PubMed
    • Google Scholar
57. 1. Sanderson E
    2. Davey Smith G
    3. Windmeijer F
    4. Bowden J

    (2019) An examination of multivariable mendelian randomization in the single-sample and two-sample summary data settings

    International Journal of Epidemiology 48:713–727.

    https://doi.org/10.1093/ije/dyy262

    • PubMed
    • Google Scholar
58. 1. Sanderson E
    2. Windmeijer F

    (2016) A weak instrument [Formula: see text]-test in linear IV models with multiple endogenous variables

    Journal of Econometrics 190:212–221.

    https://doi.org/10.1016/j.jeconom.2015.06.004

    • PubMed
    • Google Scholar
59. 1. Schmidt AF
    2. Pearce LS
    3. Wilkins JT
    4. Overington JP
    5. Hingorani AD
    6. Casas JP

    (2017) PCSK9 monoclonal antibodies for the primary and secondary prevention of cardiovascular disease

    Cochrane Database of Systematic Reviews 4:Cd011748.

    https://doi.org/10.1002/14651858.CD011748.pub2

    • PubMed
    • Google Scholar
60. 1. Schonbeck U
    2. Libby P

    (2004) Inflammation, immunity, and HMG-CoA reductase inhibitors: statins as antiinflammatory agents?

    Circulation 109:23.

    https://doi.org/10.1161/01.CIR.0000129505.34151.23

    • PubMed
    • Google Scholar
61. 1. Schooling CM
    2. Au Yeung SL
    3. Freeman G
    4. Cowling BJ

    (2013) The effect of statins on testosterone in men and women, a systematic review and meta-analysis of randomized controlled trials

    BMC Medicine 11:57.

    https://doi.org/10.1186/1741-7015-11-57

    • Google Scholar
62. 1. Schooling CM
    2. Au Yeung SL
    3. Leung GM

    (2014) Why do statins reduce cardiovascular disease more than other lipid modulating therapies?

    European Journal of Clinical Investigation 44:1135–1140.

    https://doi.org/10.1111/eci.12342

    • PubMed
    • Google Scholar
63. 1. Schooling CM

    (2016) Could androgens be relevant to partly explain why men have lower life expectancy than women?

    Journal of Epidemiology and Community Health 70:324–328.

    https://doi.org/10.1136/jech-2015-206336

    • PubMed
    • Google Scholar
64. 1. Schooling CM

    (2017) Tachykinin neurokinin 3 receptor antagonists: a new treatment for cardiovascular disease?

    The Lancet 390:709–711.

    https://doi.org/10.1016/S0140-6736(16)31648-8

    • PubMed
    • Google Scholar
65. 1. Schooling CM
    2. Luo S
    3. Au Yeung SL
    4. Thompson DJ
    5. Karthikeyan S
    6. Bolton TR
    7. Mason AM
    8. Ingelsson E
    9. Burgess S

    (2018a) Genetic predictors of testosterone and their associations with cardiovascular disease and risk factors: a mendelian randomization investigation

    International Journal of Cardiology 267:171–176.

    https://doi.org/10.1016/j.ijcard.2018.05.051

    • PubMed
    • Google Scholar
66. 1. Schooling CM
    2. Kodali H
    3. Li S
    4. Borrell LN

    (2018b) ET (Endothelin)-1 and ischemic heart disease: a mendelian randomization study

    Circulation. Genomic and Precision Medicine 11:e002026.

    https://doi.org/10.1161/CIRCGENETICS.117.002026

    • PubMed
    • Google Scholar
67. Preprint

    1. Schooling CM
    2. Lopez PM
    3. Yang Z
    4. Zhao JV
    5. Au Yeung SL
    6. Huang JV

    (2020) Use of multivariable mendelian randomization to address biases due to competing risk before recruitment

    bioRxiv.

    https://doi.org/10.1101/716621

    • Google Scholar
68. 1. Schooling CM
    2. Ng JCM

    (2019) Reproduction and longevity: a mendelian randomization study of gonadotropin-releasing hormone and ischemic heart disease

    SSM - Population Health 8:100411.

    https://doi.org/10.1016/j.ssmph.2019.100411

    • PubMed
    • Google Scholar
69. 1. Schooling CM
    2. Zhao JV

    (2019) How might bromodomain and Extra-Terminal (BET) Inhibitors operate in cardiovascular disease?

    American Journal of Cardiovascular Drugs 19:107–111.

    https://doi.org/10.1007/s40256-018-00315-3

    • PubMed
    • Google Scholar
70. 1. Serban C
    2. Sahebkar A
    3. Ursoniu S
    4. Mikhailidis DP
    5. Rizzo M
    6. Lip GY
    7. Kees Hovingh G
    8. Kastelein JJ
    9. Kalinowski L
    10. Rysz J
    11. Banach M

    (2015) A systematic review and meta-analysis of the effect of statins on plasma asymmetric dimethylarginine concentrations

    Scientific Reports 5:9902.

    https://doi.org/10.1038/srep09902

    • PubMed
    • Google Scholar
71. 1. Shimizu-Albergine M
    2. Van Yserloo B
    3. Golkowski MG
    4. Ong SE
    5. Beavo JA
    6. Bornfeldt KE

    (2016) SCAP/SREBP pathway is required for the full steroidogenic response to cyclic AMP

    PNAS 113:E5685–E5693.

    https://doi.org/10.1073/pnas.1611424113

    • PubMed
    • Google Scholar
72. Preprint

    1. Sinnott-Armstrong N
    2. Tanigawa Y
    3. Amar D
    4. Mars NJ
    5. Aguirre M
    6. Venkataraman GR
    7. Wainberg M
    8. Ollila HN
    9. Pirruccello JP
    10. Qian J

    (2019) Genetics of 38 blood and urine biomarkers in the UK biobank

    bioRxiv.

    https://doi.org/10.1101/660506

    • Google Scholar
73. 1. Sliz E
    2. Kettunen J
    3. Holmes MV
    4. Williams CO
    5. Boachie C
    6. Wang Q
    7. Männikkö M
    8. Sebert S
    9. Walters R
    10. Lin K
    11. Millwood IY
    12. Clarke R
    13. Li L
    14. Rankin N
    15. Welsh P
    16. Delles C
    17. Jukema JW
    18. Trompet S
    19. Ford I
    20. Perola M
    21. Salomaa V
    22. Järvelin MR
    23. Chen Z
    24. Lawlor DA
    25. Ala-Korpela M
    26. Danesh J
    27. Davey Smith G
    28. Sattar N
    29. Butterworth A
    30. Würtz P

    (2018) Metabolomic consequences of genetic inhibition of PCSK9 compared with statin treatment

    Circulation 138:2499–2512.

    https://doi.org/10.1161/CIRCULATIONAHA.118.034942

    • PubMed
    • Google Scholar
74. 1. Smith GD
    2. Ebrahim S

    (2003) 'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease?

    International Journal of Epidemiology 32:1–22.

    https://doi.org/10.1093/ije/dyg070

    • PubMed
    • Google Scholar
75. 1. Swerdlow DI
    2. Holmes MV
    3. Kuchenbaecker KB
    4. Engmann JE
    5. Shah T
    6. Sofat R
    7. Guo Y
    8. Chung C
    9. Peasey A
    10. Pfister R
    11. Mooijaart SP
    12. Ireland HA
    13. Leusink M
    14. Langenberg C
    15. Li KW
    16. Palmen J
    17. Howard P
    18. Cooper JA
    19. Drenos F
    20. Hardy J
    21. Nalls MA
    22. Li YR
    23. Lowe G
    24. Stewart M
    25. Bielinski SJ
    26. Peto J
    27. Timpson NJ
    28. Gallacher J
    29. Dunlop M
    30. Houlston R
    31. Tomlinson I
    32. Tzoulaki I
    33. Luan J
    34. Boer JM
    35. Forouhi NG
    36. Onland-Moret NC
    37. van der Schouw YT
    38. Schnabel RB
    39. Hubacek JA
    40. Kubinova R
    41. Baceviciene M
    42. Tamosiunas A
    43. Pajak A
    44. Topor-Madry R
    45. Malyutina S
    46. Baldassarre D
    47. Sennblad B
    48. Tremoli E
    49. de Faire U
    50. Ferrucci L
    51. Bandenelli S
    52. Tanaka T
    53. Meschia JF
    54. Singleton A
    55. Navis G
    56. Mateo Leach I
    57. Bakker SJ
    58. Gansevoort RT
    59. Ford I
    60. Epstein SE
    61. Burnett MS
    62. Devaney JM
    63. Jukema JW
    64. Westendorp RG
    65. Jan de Borst G
    66. van der Graaf Y
    67. de Jong PA
    68. Mailand-van der Zee AH
    69. Klungel OH
    70. de Boer A
    71. Doevendans PA
    72. Stephens JW
    73. Eaton CB
    74. Robinson JG
    75. Manson JE
    76. Fowkes FG
    77. Frayling TM
    78. Price JF
    79. Whincup PH
    80. Morris RW
    81. Lawlor DA
    82. Smith GD
    83. Ben-Shlomo Y
    84. Redline S
    85. Lange LA
    86. Kumari M
    87. Wareham NJ
    88. Verschuren WM
    89. Benjamin EJ
    90. Whittaker JC
    91. Hamsten A
    92. Dudbridge F
    93. Delaney JA
    94. Wong A
    95. Kuh D
    96. Hardy R
    97. Castillo BA
    98. Connolly JJ
    99. van der Harst P
    100. Brunner EJ
    101. Marmot MG
    102. Wassel CL
    103. Humphries SE
    104. Talmud PJ
    105. Kivimaki M
    106. Asselbergs FW
    107. Voevoda M
    108. Bobak M
    109. Pikhart H
    110. Wilson JG
    111. Hakonarson H
    112. Reiner AP
    113. Keating BJ
    114. Sattar N
    115. Hingorani AD
    116. Casas JP
    117. Interleukin-6 Receptor Mendelian Randomisation Analysis (IL6R MR) Consortium

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

    Lancet 379:1214–1224.

    https://doi.org/10.1016/S0140-6736(12)60110-X

    • PubMed
    • Google Scholar
76. 1. Swerdlow DI
    2. Preiss D
    3. Kuchenbaecker KB
    4. Holmes MV
    5. Engmann JEL
    6. Shah T
    7. Sofat R
    8. Stender S
    9. Johnson PCD
    10. Scott RA
    11. Leusink M
    12. Verweij N
    13. Sharp SJ
    14. Guo Y
    15. Giambartolomei C
    16. Chung C
    17. Peasey A
    18. Amuzu A
    19. Li K
    20. Palmen J
    21. Howard P
    22. Cooper JA
    23. Drenos F
    24. Li YR
    25. Lowe G
    26. Gallacher J
    27. Stewart MCW
    28. Tzoulaki I
    29. Buxbaum SG
    30. van der A DL
    31. Forouhi NG
    32. Onland-Moret NC
    33. van der Schouw YT
    34. Schnabel RB
    35. Hubacek JA
    36. Kubinova R
    37. Baceviciene M
    38. Tamosiunas A
    39. Pajak A
    40. Topor-Madry R
    41. Stepaniak U
    42. Malyutina S
    43. Baldassarre D
    44. Sennblad B
    45. Tremoli E
    46. de Faire U
    47. Veglia F
    48. Ford I
    49. Jukema JW
    50. Westendorp RGJ
    51. de Borst GJ
    52. de Jong PA
    53. Algra A
    54. Spiering W
    55. der Zee AHM
    56. Klungel OH
    57. de Boer A
    58. Doevendans PA
    59. Eaton CB
    60. Robinson JG
    61. Duggan D
    62. Kjekshus J
    63. Downs JR
    64. Gotto AM
    65. Keech AC
    66. Marchioli R
    67. Tognoni G
    68. Sever PS
    69. Poulter NR
    70. Waters DD
    71. Pedersen TR
    72. Amarenco P
    73. Nakamura H
    74. McMurray JJV
    75. Lewsey JD
    76. Chasman DI
    77. Ridker PM
    78. Maggioni AP
    79. Tavazzi L
    80. Ray KK
    81. Seshasai SRK
    82. Manson JE
    83. Price JF
    84. Whincup PH
    85. Morris RW
    86. Lawlor DA
    87. Smith GD
    88. Ben-Shlomo Y
    89. Schreiner PJ
    90. Fornage M
    91. Siscovick DS
    92. Cushman M
    93. Kumari M
    94. Wareham NJ
    95. Verschuren WMM
    96. Redline S
    97. Patel SR
    98. Whittaker JC
    99. Hamsten A
    100. Delaney JA
    101. Dale C
    102. Gaunt TR
    103. Wong A
    104. Kuh D
    105. Hardy R
    106. Kathiresan S
    107. Castillo BA
    108. van der Harst P
    109. Brunner EJ
    110. Tybjaerg-Hansen A
    111. Marmot MG
    112. Krauss RM
    113. Tsai M
    114. Coresh J
    115. Hoogeveen RC
    116. Psaty BM
    117. Lange LA
    118. Hakonarson H
    119. Dudbridge F
    120. Humphries SE
    121. Talmud PJ
    122. Kivimäki M
    123. Timpson NJ
    124. Langenberg C
    125. Asselbergs FW
    126. Voevoda M
    127. Bobak M
    128. Pikhart H
    129. Wilson JG
    130. Reiner AP
    131. Keating BJ
    132. Hingorani AD
    133. Sattar N

    (2015) HMG-coenzyme A reductase inhibition, type 2 diabetes, and bodyweight: evidence from genetic analysis and randomised trials

    The Lancet 385:351–361.

    https://doi.org/10.1016/S0140-6736(14)61183-1

    • Google Scholar
77. 1. Tardif JC
    2. Kouz S
    3. Waters DD
    4. Bertrand OF
    5. Diaz R
    6. Maggioni AP
    7. Pinto FJ
    8. Ibrahim R
    9. Gamra H
    10. Kiwan GS
    11. Berry C
    12. López-Sendón J
    13. Ostadal P
    14. Koenig W
    15. Angoulvant D
    16. Grégoire JC
    17. Lavoie MA
    18. Dubé MP
    19. Rhainds D
    20. Provencher M
    21. Blondeau L
    22. Orfanos A
    23. L'Allier PL
    24. Guertin MC
    25. Roubille F

    (2019) Efficacy and safety of Low-Dose colchicine after myocardial infarction

    New England Journal of Medicine 381:2497–2505.

    https://doi.org/10.1056/NEJMoa1912388

    • PubMed
    • Google Scholar
78. 1. Tran DQ
    2. Ramos EH
    3. Belsham DD

    (2016) Induction of gnrh mRNA expression by the ω-3 polyunsaturated fatty acid docosahexaenoic acid and the saturated fatty acid palmitate in a GnRH-synthesizing neuronal cell model, mHypoA-GnRH/GFP

    Molecular and Cellular Endocrinology 426:125–135.

    https://doi.org/10.1016/j.mce.2016.02.019

    • PubMed
    • Google Scholar
79. 1. Tremellen K
    2. McPhee N
    3. Pearce K
    4. Benson S
    5. Schedlowski M
    6. Engler H

    (2018) Endotoxin-initiated inflammation reduces testosterone production in men of reproductive age

    American Journal of Physiology-Endocrinology and Metabolism 314:E206–E213.

    https://doi.org/10.1152/ajpendo.00279.2017

    • PubMed
    • Google Scholar
80. 1. van Kesteren PJ
    2. Kooistra T
    3. Lansink M
    4. van Kamp GJ
    5. Asscheman H
    6. Gooren LJ
    7. Emeis JJ
    8. Vischer UM
    9. Stehouwer CD

    (1998) The effects of sex steroids on plasma levels of marker proteins of endothelial cell functioning

    Thrombosis and Haemostasis 79:1029–1033.

    https://doi.org/10.1055/s-0037-1615115

    • PubMed
    • Google Scholar
81. 1. Willer CJ
    2. Schmidt EM
    3. Sengupta S
    4. Peloso GM
    5. Gustafsson S
    6. Kanoni S
    7. Ganna A
    8. Chen J
    9. Buchkovich ML
    10. Mora S
    11. Beckmann JS
    12. Bragg-Gresham JL
    13. Chang HY
    14. Demirkan A
    15. Den Hertog HM
    16. Do R
    17. Donnelly LA
    18. Ehret GB
    19. Esko T
    20. Feitosa MF
    21. Ferreira T
    22. Fischer K
    23. Fontanillas P
    24. Fraser RM
    25. Freitag DF
    26. Gurdasani D
    27. Heikkilä K
    28. Hyppönen E
    29. Isaacs A
    30. Jackson AU
    31. Johansson Å
    32. Johnson T
    33. Kaakinen M
    34. Kettunen J
    35. Kleber ME
    36. Li X
    37. Luan J
    38. Lyytikäinen LP
    39. Magnusson PKE
    40. Mangino M
    41. Mihailov E
    42. Montasser ME
    43. Müller-Nurasyid M
    44. Nolte IM
    45. O'Connell JR
    46. Palmer CD
    47. Perola M
    48. Petersen AK
    49. Sanna S
    50. Saxena R
    51. Service SK
    52. Shah S
    53. Shungin D
    54. Sidore C
    55. Song C
    56. Strawbridge RJ
    57. Surakka I
    58. Tanaka T
    59. Teslovich TM
    60. Thorleifsson G
    61. Van den Herik EG
    62. Voight BF
    63. Volcik KA
    64. Waite LL
    65. Wong A
    66. Wu Y
    67. Zhang W
    68. Absher D
    69. Asiki G
    70. Barroso I
    71. Been LF
    72. Bolton JL
    73. Bonnycastle LL
    74. Brambilla P
    75. Burnett MS
    76. Cesana G
    77. Dimitriou M
    78. Doney ASF
    79. Döring A
    80. Elliott P
    81. Epstein SE
    82. Ingi Eyjolfsson G
    83. Gigante B
    84. Goodarzi MO
    85. Grallert H
    86. Gravito ML
    87. Groves CJ
    88. Hallmans G
    89. Hartikainen AL
    90. Hayward C
    91. Hernandez D
    92. Hicks AA
    93. Holm H
    94. Hung YJ
    95. Illig T
    96. Jones MR
    97. Kaleebu P
    98. Kastelein JJP
    99. Khaw KT
    100. Kim E
    101. Klopp N
    102. Komulainen P
    103. Kumari M
    104. Langenberg C
    105. Lehtimäki T
    106. Lin SY
    107. Lindström J
    108. Loos RJF
    109. Mach F
    110. McArdle WL
    111. Meisinger C
    112. Mitchell BD
    113. Müller G
    114. Nagaraja R
    115. Narisu N
    116. Nieminen TVM
    117. Nsubuga RN
    118. Olafsson I
    119. Ong KK
    120. Palotie A
    121. Papamarkou T
    122. Pomilla C
    123. Pouta A
    124. Rader DJ
    125. Reilly MP
    126. Ridker PM
    127. Rivadeneira F
    128. Rudan I
    129. Ruokonen A
    130. Samani N
    131. Scharnagl H
    132. Seeley J
    133. Silander K
    134. Stančáková A
    135. Stirrups K
    136. Swift AJ
    137. Tiret L
    138. Uitterlinden AG
    139. van Pelt LJ
    140. Vedantam S
    141. Wainwright N
    142. Wijmenga C
    143. Wild SH
    144. Willemsen G
    145. Wilsgaard T
    146. Wilson JF
    147. Young EH
    148. Zhao JH
    149. Adair LS
    150. Arveiler D
    151. Assimes TL
    152. Bandinelli S
    153. Bennett F
    154. Bochud M
    155. Boehm BO
    156. Boomsma DI
    157. Borecki IB
    158. Bornstein SR
    159. Bovet P
    160. Burnier M
    161. Campbell H
    162. Chakravarti A
    163. Chambers JC
    164. Chen YI
    165. Collins FS
    166. Cooper RS
    167. Danesh J
    168. Dedoussis G
    169. de Faire U
    170. Feranil AB
    171. Ferrières J
    172. Ferrucci L
    173. Freimer NB
    174. Gieger C
    175. Groop LC
    176. Gudnason V
    177. Gyllensten U
    178. Hamsten A
    179. Harris TB
    180. Hingorani A
    181. Hirschhorn JN
    182. Hofman A
    183. Hovingh GK
    184. Hsiung CA
    185. Humphries SE
    186. Hunt SC
    187. Hveem K
    188. Iribarren C
    189. Järvelin MR
    190. Jula A
    191. Kähönen M
    192. Kaprio J
    193. Kesäniemi A
    194. Kivimaki M
    195. Kooner JS
    196. Koudstaal PJ
    197. Krauss RM
    198. Kuh D
    199. Kuusisto J
    200. Kyvik KO
    201. Laakso M
    202. Lakka TA
    203. Lind L
    204. Lindgren CM
    205. Martin NG
    206. März W
    207. McCarthy MI
    208. McKenzie CA
    209. Meneton P
    210. Metspalu A
    211. Moilanen L
    212. Morris AD
    213. Munroe PB
    214. Njølstad I
    215. Pedersen NL
    216. Power C
    217. Pramstaller PP
    218. Price JF
    219. Psaty BM
    220. Quertermous T
    221. Rauramaa R
    222. Saleheen D
    223. Salomaa V
    224. Sanghera DK
    225. Saramies J
    226. Schwarz PEH
    227. Sheu WH
    228. Shuldiner AR
    229. Siegbahn A
    230. Spector TD
    231. Stefansson K
    232. Strachan DP
    233. Tayo BO
    234. Tremoli E
    235. Tuomilehto J
    236. Uusitupa M
    237. van Duijn CM
    238. Vollenweider P
    239. Wallentin L
    240. Wareham NJ
    241. Whitfield JB
    242. Wolffenbuttel BHR
    243. Ordovas JM
    244. Boerwinkle E
    245. Palmer CNA
    246. Thorsteinsdottir U
    247. Chasman DI
    248. Rotter JI
    249. Franks PW
    250. Ripatti S
    251. Cupples LA
    252. Sandhu MS
    253. Rich SS
    254. Boehnke M
    255. Deloukas P
    256. Kathiresan S
    257. Mohlke KL
    258. Ingelsson E
    259. Abecasis GR
    260. Global Lipids Genetics Consortium

    (2013) Discovery and refinement of loci associated with lipid levels

    Nature Genetics 45:1274–1283.

    https://doi.org/10.1038/ng.2797

    • PubMed
    • Google Scholar
82. 1. Würtz P
    2. Wang Q
    3. Soininen P
    4. Kangas AJ
    5. Fatemifar G
    6. Tynkkynen T
    7. Tiainen M
    8. Perola M
    9. Tillin T
    10. Hughes AD
    11. Mäntyselkä P
    12. Kähönen M
    13. Lehtimäki T
    14. Sattar N
    15. Hingorani AD
    16. Casas JP
    17. Salomaa V
    18. Kivimäki M
    19. Järvelin MR
    20. Davey Smith G
    21. Vanhala M
    22. Lawlor DA
    23. Raitakari OT
    24. Chaturvedi N
    25. Kettunen J
    26. Ala-Korpela M

    (2016) Metabolomic profiling of statin use and genetic inhibition of HMG-CoA reductase

    Journal of the American College of Cardiology 67:1200–1210.

    https://doi.org/10.1016/j.jacc.2015.12.060

    • PubMed
    • Google Scholar
83. 1. Zhao JV

    (2018) Coagulation factors and the risk of ischemic heart disease, a mendelian randomization study

    Circulation. Genomic and Precision Medicine 11:001956.

    https://doi.org/10.1161/CIRCGEN.117.001956

    • Google Scholar
84. 1. Zhao JV
    2. Luo S
    3. Schooling CM

    (2019) Sex-specific mendelian randomization study of genetically predicted insulin and cardiovascular events in the UK biobank

    Communications Biology 2:332.

    https://doi.org/10.1038/s42003-019-0579-z

    • PubMed
    • Google Scholar
85. 1. Zhao JV
    2. Schooling CM

    (2020) The role of testosterone in chronic kidney disease and kidney function in men and women: a bi-directional mendelian randomization study in the UK biobank

    BMC Medicine 18:122.

    https://doi.org/10.1186/s12916-020-01594-x

    • PubMed
    • Google Scholar

Acceptance summary:

The authors performed a Mendelian randomization analysis of statin, PCSK9 inhibitor and ezetimibe use on ischemic heart disease to address a clinically relevant question – could the apparent additional benefits of statins and other similar treatments, over and above lowering lipids, be mediated by favourable effects on testosterone levels? They found a potentially interesting link between statins and testosterone – genetic variants in the gene encoding the HMGcoR receptor that are known to alter LDL-cholesterol levels.

Decision letter after peer review:

Thank you for submitting your article "Pleiotropic effects of statins on ischemic heart disease: a Mendelian Randomization study in the UK Biobank" for consideration by eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Matthias Barton as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Timothy Frayling (Reviewer #3).

The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.

We would like to draw your attention to changes in our revision policy that we have made in response to COVID-19 (https://elifesciences.org/articles/57162). Specifically, when editors judge that a submitted work as a whole belongs in eLife but that some conclusions require a modest amount of additional new data, as they do with your paper, we are asking that the manuscript be revised to either limit claims to those supported by data in hand, or to explicitly state that the relevant conclusions require additional supporting data.

Our expectation is that the authors will eventually carry out the additional experiments and report on how they affect the relevant conclusions either in a preprint on bioRxiv or medRxiv, or if appropriate, as a Research Advance in eLife, either of which would be linked to the original paper.

Title: Please modify the title to comply with eLife requirements. The title should not exceed 120 characters. Two part titles containing a colon punctuation mark (":") are not allowed.

Summary:

This is a neat study using genetics to ask a clinically relevant question – could the apparent additional benefits of statins and other similar treatments, over and above lowering lipids, be mediated by favourable effects on testosterone levels? The question becomes a little more complicated in that higher levels of testosterone tend to be beneficial for men's metabolic health but adverse for women and so the authors have used UK Biobank data, where it is easy to stratify analyses by sex.

The authors performed a Mendelian randomization (MR) analysis of statin, PCSK9 inhibitor and ezetimibe use on ischemic heart disease (IHD), and if these effects were mediated by testosterone in men or women using univariable and multivariable MR.

They have found a potentially interesting link between statins and testosterone – genetic variants in the gene encoding the HMGcoR receptor that are known to alter LDL-cholesterol levels (And so we know they are very good mimics for the on target effects of statins) and testosterone levels in men. The direction is that the alleles associated with lower LDL-cholesterol and therefore mimic statins also lower testosterone. There were no links for the other genes encoding the other targets of lipid lowering therapy. In addition, they assessed whether the genetic variants corresponding to anti-inflammatory effects of anakinra or tocilizumab use had effects on testosterone and IHD. The topic is very interesting.

Revisions for this paper:

1) The main finding is that, when adjusting for testosterone levels, the protective benefit of statins via the HMGcoR on heart disease attenuates from 0.54 odds ratio to 0.73 in men (when using the MR egger approach. There are wide 95% CIs around these estimates but the results imply that the testosterone effect is similar to the lowering LDL-cholesterol effect, which is extremely counter intuitive given that we know about LDL-cholesterol and heart disease. The reliance on a sensitivity analysis – MR egger – when the main analysis showed evidence of pleiotropy means the results could be very sensitive to slight differences in the parameters used. (SNPs and models).

2) A second concern is that the genetic variants in HMGCoR are associated with BMI and adiposity (see Swerdlow et al. Lancet publication using the HMGcoR SNPs to show an association with diabetes) whereas the others are less so. The authors need to assess the potential for the effect being mediated via BMI/adiposity rather than testosterone. Such a mechanism would be consistent with sex differences because women are at lower risk of heart disease than men for a given BMI. Likewise, many of the testosterone SNPs are highly pleiotropic, with many primarily associated with adiposity and insulin resistance and lipid levels. My concern is that the pleiotropy will lead to false inferences. The authors have partially addressed this with some approaches such as MR-Egger, but these are not infallible. The multivariable MR needs to include the adiposity measures. multivariable MR with adiposity measures, steiger filtering of SNPs that have larger effects on metabolic potentially mediating traits than they do on lipids or testosterone could be alternative approaches to try.

3) It would be relevant to see the results replicated in a two sample MR setting where the IHD cases are not from the same dataset as the testosterone and LDL SNPs. The heart disease GWAS consortia have not analysed separately by sex, which makes this very difficult, but there may be other large studies where this is possible ? If this cannot be done readily it should be discussed and noted as a limitation.

4) All of this means that the conclusion that "statins partially operate on IHD by reducing testosterone in men" is likely overstated.

5) A further major concern is the selection of SNPs in vicinity of the HMGCR, NCP1L1 and PCSK9 genes as instruments for the statin, PCSK9 and ezetimibe use, respectively. For example, a lookup of the instrument included for ezetimibe use (rs10260606) in the publicly available sex-combined GWAS on ezetimibe treatment (http://www.nealelab.is/uk-biobank) resulted in an association p-value of 0.02 being a rather weak instrument. Also in the referenced publication of Ference et al., 2019, that was used as basis for the instrument selection, these SNPs were (to my understanding) selected to assess potential drug targets, and not to reflect medication use.

6) Which statistical model was used to calculate the sex-specific genetic associations with IHD (I assume logistic regression)? Please provide the number of cases and controls per sex-stratum that were finally included in this association analysis also in the main text (not only in the Abstract).

Please provide the p-value level that was applied to declare significance of the results.

7) By conducting a sex-stratified analyses, have the authors considered the impact of collider bias. As if the genetic instruments are associated with sex and the outcome measures are also associated with sex, by stratifying based on a common cause of the exposure and outcome, a distorted/erroneous association may occur.

8) Please give more details about multivariate MR analysis. E.g. did you pool the genetic proxies for statin and testosterone all together and then extracted the associations of all these SNPs with LDL-cholesterol and testosterone and fitted it in one model?

Revisions expected in follow-up work:

1) As the authors aim to investigate sex-specific benefits of statin therapies, why only focus on testosterone. I think you should also look at 17β-oestradiol and SHBG (both available in UKBB), and conduct a multivariate MR to determine which would be the independent causal factor for IHD among the three (e.g. PMID: 32203549) and then forward the independent hormones/proteins in a multivariate model with LDL-cholesterol or apoB (drug targets of lipid lowering drugs) to determine whether these hormones would exert independent roles in contributing to IHD beyond the lipids. If this is the case, then it makes sense to further determine whether the genetic instruments of the drug targets show a sex-specific association with the sex hormones. This issue could be addressed in the Discussion.

Revisions for this paper:

1) The main finding is that, when adjusting for testosterone levels, the protective benefit of statins via the HMGcoR on heart disease attenuates from 0.54 odds ratio to 0.73 in men (when using the MR egger approach. There are wide 95% CIs around these estimates but the results imply that the testosterone effect is similar to the lowering LDL-cholesterol effect, which is extremely counter intuitive given that we know about LDL-cholesterol and heart disease. The reliance on a sensitivity analysis – MR egger – when the main analysis showed evidence of pleiotropy means the results could be very sensitive to slight differences in the parameters used. (SNPs and models).

Please accept our apologies for a counter-intuitive finding. We are in no way trying to question the role of LDL-cholesterol in heart disease. We are instead trying to address two important health issues. First, the higher rates of heart disease in men than women,¹ which are a factor contributing to shorter life expectancy in men than women. Second, the apparent pleiotropic effects of statins, because their explication might facilitate the discovery of much needed new treatments for cardiovascular disease.

As regards, the size of the estimates for statins via HMGcoR on ischemic heart disease (IHD). After allowing for testosterone, the estimates in men changed from being larger than the estimates in women to being very similar to the estimates in women, as expected. However, we agree that the magnitude of the estimates in Mendelian randomization (MR) is difficult to interpret because MR estimates capture lifetime effects rather than a possibly short-term exposure. We have clarified this important point to the Discussion, by making the following change:

From: “Lastly, Mendelian randomization assesses the lifelong effects of an endogenous exposure rather than short-term effects of an interventions assessed in an RCT, nevertheless our estimates for statins on IHD are comparable with meta-analyses of statin trials considering similar outcomes.²”

To: “Lastly, Mendelian randomization assesses the lifelong effects of an endogenous exposure rather than short-term effects of an interventions assessed in an RCT. Our estimates give an indication of the role of the exposures rather than the exact effects of the corresponding interventions. Nevertheless, our estimates for genetically mimicked statins on IHD are comparable with meta-analyses of statin trials considering similar outcomes.² ”

Please accept our apologies for being unclear about the use of the MR-Egger estimate. For univariable mendelian randomization (MR), MR-Egger is a sensitivity analysis. For the multivariable MR, when the analysis is multivariable by design with few genetic variants specifically predicting the exposures, as here for statins, then the MR-Egger estimate is most likely to be reliable.³ It is also recommended to use the MR-Egger multivariable estimate when the multivariable MR-Egger intercept is significant,³ as here. We have added this point to the Materials and methods:

“Given this analysis is multivariable by design with few genetic variants available to mimic the effects of statins, we used the multivariable MR-Egger estimates³.”

We also amended the Results to be clearer by making the following change:

From: “Correspondingly, the multivariable MR-Egger intercepts were significant in men and women, so we used the MR-Egger estimates ³.”

To: “Correspondingly, the multivariable MR-Egger intercepts were significant in men and women.”

2) A second concern is that the genetic variants in HMGCoR are associated with BMI and adiposity (see Swerdlow et al. Lancet publication using the HMGcoR SNPs to show an association with diabetes) whereas the others are less so. The authors need to assess the potential for the effect being mediated via BMI/adiposity rather than testosterone. Such a mechanism would be consistent with sex differences because women are at lower risk of heart disease than men for a given BMI. Likewise, many of the testosterone SNPs are highly pleiotropic, with many primarily associated with adiposity and insulin resistance and lipid levels. My concern is that the pleiotropy will lead to false inferences. The authors have partially addressed this with some approaches such as MR-Egger, but these are not infallible. The multivariable MR needs to include the adiposity measures. multivariable MR with adiposity measures, steiger filtering of SNPs that have larger effects on metabolic potentially mediating traits than they do on lipids or testosterone could be alternative approaches to try.

Thank you very much indeed for raising these concerns as to whether the pleiotropic effects of statins are via BMI rather than testosterone. We have considered this possibility in detail.

First, we checked, as expected from Swerdlow,⁴ that statins increase BMI in men and women (beta 0.33 effect size, 95% confidence interval (CI) 0.20 to 0.47, in men and 0.31, 95% CI 0.22 to 0.40, in women, Author response table 1 below), using published instruments for BMI⁵ and rs12916-T to mimic effects of statins, and the UK Biobank sex-specific genetic associations with BMI (http://www.nealelab.is/uk-biobank).

Second, we considered the effect of BMI⁵ and genetically mimicked statins on IHD in univariable and multivariable MR. In univariable MR, as expected, BMI increased IHD (odds ratio (OR) 1.33, 95% confidence interval (CI) 1.15 to 1.55, in men and OR 1.23, 95% CI 1.02 to 1.48, in women) and genetically mimicked statins decreased IHD (OR 0.54, 95% CI 0.38 to 0.79 in men and had an estimate in the same direction in women OR 0.87, 95% CI 0.59 to 1.27, Author response table 1). In multivariable MR, considering the effects of BMI and genetically mimicked statins (rs12916-T) together on IHD, BMI had very similar associations with IHD as in the univariable analysis (OR 1.35, 95% CI 1.15 to 1.57 in men and women OR 1.27, 95% CI 1.07 to 1.52). In contrast, in the multivariable analysis including BMI genetically mimicked statins had a larger decreasing effect on IHD (OR 0.46, 95% CI 0.35 to 0.65 in men and OR 0.68, 95% CI 0.53 to 0.86 in women). As such, the protective effect of statins on IHD was increased by additionally considering BMI, so statins protective effect on IHD does not appear to be explained away by statins increasing BMI, in fact including BMI makes the protective effects of statins on IHD more marked. The multivariable MR is probably giving a larger effect size for statins than the univariable MR because it is showing the effects of “statins” after removing the harmful aspect driven by raising BMI. As such, statins raising BMI is unlikely to be a pathway by which statins protect against IHD, and we did not consider it further in the multivariable MR to keep the paper as simple and focused as possible. We have added this point in the Discussion as follows:

“We did not include body mass index (BMI) as a risk factor explaining the effect of mimicking statins on IHD, because statins increase BMI⁴ and decrease the risk of IHD, so including them in the multivariable analysis may inflate the effect of mimicking statins on IHD, rather than explaining part of the effect”.

We agree that the testosterone SNPs may have pleiotropic effects possibly acting via adiposity, insulin resistance or lipids. Downstream effects of the exposures help explain how the exposure affects the outcome. Downstream effects of the exposures generate vertical pleiotropy which does not bias MR estimates, because it is a pathway through the exposure, rather than the genetic instrument operating by a pathway distinct from the exposure. Pleiotropy biases MR estimates, when it is horizontal pleiotropy, and the genetic instrument affects the outcome by a pathway distinct from the exposure. We also checked whether the testosterone SNPs operated via adiposity or lipids. The testosterone SNPs did not affect BMI in men, consistent with a previous Mendelian randomization study,⁶ as given below (Author response table 2). In contrast testosterone raised BMI in women and reduced LDL-cholesterol, however, testosterone in women does not affect IHD. We have added that these points in the Discussion.

“Testosterone’s effects on IHD in men could be via adiposity, insulin or LDL-cholesterol rather than via testosterone. […] We could not test whether testosterone in men affects insulin because of the lack of an insulin GWAS including the X chromosome.”

Author response table 1: Univariable Mendelian randomization estimates for genetically mimicked statins (rs21916-T) and genetically predicted BMI (95 variants*) on IHD in men and women in the UK Biobank.

*after removing one variant predicting BMI (rs2112347) that was correlated with rs12916

Author response table 2: Univariable sex-specific Mendelian randomization estimates for testosterone on body mass index and LDL-cholesterol in men and women in the UK Biobank.

As suggested, we used Steiger filtering to check the direction from testosterone to IHD,⁷ which confirmed that the direction of causality was from testosterone to IHD in men and women. We have added this point in the Materials and methods and the Results.

In the Materials and methods:

“Steiger filtering was used to check the directionality between testosterone and IHD.⁷”

In the Results;

“Steiger filtering indicated directionality from testosterone to IHD in men and women.”

3) It would be relevant to see the results replicated in a two sample MR setting where the IHD cases are not from the same dataset as the testosterone and LDL SNPs. The heart disease GWAS consortia have not analysed separately by sex, which makes this very difficult, but there may be other large studies where this is possible ? If this cannot be done readily it should be discussed and noted as a limitation.

Please accept our apologies for being unclear. The testosterone SNPs were derived from the UK Biobank, but were also replicated in three independent studies,⁸ so they do have some replication. The LDL SNPs for statins, PCSK9 inhibitors and ezetimibe do not come from the same dataset as the IHD cases. The LDL SNPs come from Ference et al., as used previously to mimic the effects of lipid modifiers⁹. Only the beta co-efficients for the LDL SNPs on LDL come from the same study as the IHD cases. The beta coefficients are simply scaling factors which are quite similar to those from the Global Lipids Genetic consortium (GLGC), as we show in Supplementary file 1. We have added a column to that table to make it clearer. However, the GLGC does not have sex-specific genetic associations with LDL-c, so we used UK Biobank which gives sex-specific effect of the SNPs mimicking statins on LDL-c. We have explained this point, as given below:

“Strong (p-value<5×10^-8), independent (r²<0.05), sex-specific genetic predictors of testosterone were extracted from a published genome wide association study (GWAS) based on the UK Biobank and replicated in three independent studies (CHARGE Consortium, Twins UK and EPIC-Norfolk) ^{10 11}.”

“Genetic effects mimicking statins, PCSK9 inhibitors and ezetimibe were expressed in sex-specific effect sizes of LDL-cholesterol reduction taken from the largest available sex-specific GWAS summary statistics, i.e., the UK Biobank (http://www.nealelab.is/uk-biobank).”

Thank you for this suggestion about replication, we would be very keen to do so. We cannot identify any other sex-specific GWAS of IHD. Most GWAS of IHD are composed largely of men, so one option would be to repeat the analysis using an IHD GWAS of mainly men. However, there are two issues with that approach. First a GWAS of mainly men does not provide replication by sex. Second, almost all previous GWAS are autosomal, so we are missing key items in the testosterone instrument from the X chromosome, which means this is not a true replication. We have clarified this point in the Discussion, by making the following change:

From: “Replication based on another large IHD GWAS, such as CARDIoGRAM¹², would be ideal, but it is less intensively genotyped than the UK Biobank and excludes the X chromosome, making it unsuitable.”

To: “Replication based on another large sex-specific IHD GWAS where the IHD cases are not from the same study as the testosterone instruments, would be ideal. However, sex-specific summary statistics are not available for large existing IHD GWAS, such as CARDIoGRAM.¹² Moreover, the UK Biobank has the advantage of being very intensively genotyped and including the X chromosome, which is important for testosterone, but is not usually included in publicly available summary statistics. Lack of replication is a limitation of this study.”

4) All of this means that the conclusion that "statins partially operate on IHD by reducing testosterone in men" is likely overstated.

We fully accept that no study is ever definitive, so we have amended this sentence to say

“statins could partially operate on IHD by reducing testosterone in men".

5) A further major concern is the selection of SNPs in vicinity of the HMGCR, NCP1L1 and PCSK9 genes as instruments for the statin, PCSK9 and ezetimibe use, respectively. For example, a lookup of the instrument included for ezetimibe use (rs10260606) in the publicly available sex-combined GWAS on ezetimibe treatment (http://www.nealelab.is/uk-biobank) resulted in an association p-value of 0.02 being a rather weak instrument.

Thank you for these comments. Please accept our apologies for being unclear, we did not obtain the instruments for statins, PCSK9 or ezetimibe use from the UK Biobank GWAS of self-reports of drug use. There were three reasons for this choice. First, for comparability with other studies we used instruments from a previous study.⁹ Second, self-reports are not always very reliable, so a GWAS based on self-reports may generate instruments that represent noise rather than signal. This may be why MR studies exploring the effects of medications usually use genetic predictors related to the drug target, based on biological pathways so as to minimise potential bias. Third, the SNPs for statins, PCSK9 or ezetimibe are to some extent selected on functional rather than statistical grounds from relevant genes which encode proteins targeted by these drugs. This approach has the advantage of being functionally relevant, and so hopefully is relevant to the drugs biological effect, and shows the effect of the drug. The SNPs we used all had F-statistics greater than 10. On reflection, we realize that the language used to describe exposure was unclear as “effects of statin use” and we have amended it throughout, from “effects of statin use”, to “mimic the effect of statins” or similar. Thank you for pointing this out.

Also in the referenced publication of Ference et al., 2019, that was used as basis for the instrument selection, these SNPs were (to my understanding) selected to assess potential drug targets, and not to reflect medication use.

Thank you for these comments. As noted in the previous comment, we found it very tricky to find exactly the right words to describe the instruments for statins, PCSK9 inhibitors and ezetimibe clearly. We have carefully reviewed the underlying publications to find more appropriate wording to describe use of these variants in similar situations. As given below the original publications use the words “mimic the effect of [the drug] or “as a proxy for [drug] treatment”.

Ference et al., 2019, used the wording “Genetic variants that mimic the effect of …statins”⁹.

Ference et al., 2016, used the wording “We constructed genetic scores that mimic the effect of PCSK9 inhibitors and the effect of statins”¹³.

Ference et al., 2015, used the wording “To compare the biological effect of lower LDL-cholesterol mediated by inhibition of NCP1L1, HMGCR, or both on the risk of CHD, and to provide a context for interpreting the results of IMPROVE-IT, we sought to compare the effect of naturally random allocation to lower LDL-cholesterol on the risk of CHD mediated by genetic polymorphisms in the NPC1L1 gene (as a proxy for ezetimibe treatment), the HMGCR gene (as a proxy for statin treatment)”.

We have amended the paper throughout to be consistent with the previous usage and to use the expression “mimic the effect of” or similar throughout.

6) Which statistical model was used to calculate the sex-specific genetic associations with IHD (I assume logistic regression)? Please provide the number of cases and controls per sex-stratum that were finally included in this association analysis also in the main text (not only in the Abstract).

Thank you for pointing out this oversight, we have added that the associations with IHD were “obtained using logistic regression”. Information about the number of cases is given in the main text.

Please provide the p-value level that was applied to declare significance of the results.

We have added this point in the text:

“Given this is a hypothesis driven study, with a positive control, we used a statistical significance level of 0.05.”

7) By conducting a sex-stratified analyses, have the authors considered the impact of collider bias. As if the genetic instruments are associated with sex and the outcome measures are also associated with sex, by stratifying based on a common cause of the exposure and outcome, a distorted/erroneous association may occur.

Thank you very much indeed for raising the issue of collider bias. Stratifying on a common cause of exposure and outcome is typically used to address confounding, i.e., it removes bias rather than adds bias.¹⁴ Collider bias (or selection bias) typically occurs if there is selection (or stratification) on a common effect of the exposure and outcome,¹⁴ i.e., here the genetic variants and IHD. We do not think that sex is a result of the genetic variants or of IHD and hence collider bias is very unlikely upon stratification by sex.

8) Please give more details about multivariate MR analysis. E.g. did you pool the genetic proxies for statin and testosterone all together and then extracted the associations of all these SNPs with LDL-cholesterol and testosterone and fitted it in one model?

Please accept our apologies for not explaining these points clearly. This is exactly the process we carried out. We have added this point, as follows:

“In the multivariable MR, we pooled the genetic instruments mimicking statins and the genetic instruments for testosterone together, extracted their associations with LDL-cholesterol and testosterone and fitted one multivariable model.”

Revisions expected in follow-up work:

1) As the authors aim to investigate sex-specific benefits of statin therapies, why only focus on testosterone. I think you should also look at 17β-oestradiol and SHBG (both available in UKBB), and conduct a multivariate MR to determine which would be the independent causal factor for IHD among the three (e.g. PMID: 32203549) and then forward the independent hormones/proteins in a multivariate model with LDL-cholesterol or apoB (drug targets of lipid lowering drugs) to determine whether these hormones would exert independent roles in contributing to IHD beyond the lipids. If this is the case, then it makes sense to further determine whether the genetic instruments of the drug targets show a sex-specific association with the sex hormones. This issue could be addressed in the Discussion.

Thank you very much indeed for asking why we focused on testosterone. We have for many years been examining the theory that statins operate specifically by testosterone,^{15 16} as presented explicitly in Figure 1 of the paper. In contrast, we have not been investigating the role of estrogen in IHD. Large randomized controlled trials have shown that hormone replacement therapy in women does not protect against IHD,¹⁷ and estrogen does not protect against IHD in men.¹⁸ Nevertheless, the possibility of a role of estrogen in IHD remains, because hormone replacement therapy may differ from naturally occurring estrogen in women and in men the Coronary Drug project estrogen arms were stopped early. Sex hormone binding globulin (SHBG) could be an independent influence on IHD or an important moderator of the role of testosterone, which deserves investigation. However, as shown in Author response table 3 (below) estrogen and SHBG do not seem relevant to specifically the effects of statins on IHD in men.

Author response table 3: Mendelian randomization estimates for effects of genetically mimicked statins (per effect sizes of LDL-cholesterol) on estrogen and SHBG in men and women from the UK Biobank summary statistics (http://www.nealelab.is/uk-biobank).

The unit of LDL-cholesterol is approximately 0.83mm/L.

We agree completely that it is a very important question to assess the independent role of testosterone, 17β-oestradiol and SHBG in IHD by sex and the extent to which they operate via LDL-cholesterol or apoB. Currently there is a dearth of comprehensive genetic instruments for estrogen. Ruth et al., by far the largest GWAS of sex hormones, only provided genetic instruments for estrogen in men, but not for estrogen in women.¹⁹ We have explained this future work in the Discussion as follows.

“Here, we present a hypothesis driven study examining the role of testosterone in mediating the effect of specifically statins in IHD, particularly in men. Future work could encompass a comprehensive sex-specific multivariable MR to confirm the role of sex hormones and sex hormone binding globulin in IHD as well as any mediation of their effects by key lipids, such as LDL-cholesterol or apoB. This work would be facilitated by the development of published genetic instruments for estrogen in women.

References:

https://wwwbiorxivorg/content/101101/660506v1

1) Moran AE, Tzong KY, Forouzanfar MH, et al. Variations in ischemic heart disease burden by age, country, and income: the Global Burden of Diseases, Injuries, and Risk Factors 2010 study. Glob Heart 2014;9(1):91-9. doi: 10.1016/j.gheart.2013.12.007 [published Online First: 2014/07/01]2) Fulcher J, O'Connell R, Voysey M, 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 (London, England) 2015;385(9976):1397-405. doi: 10.1016/s0140-6736(14)61368-4 [published Online First: 2015/01/13]3) Rees JMB, Wood AM, Burgess S. Extending the MR-Egger method for multivariable Mendelian randomization to correct for both measured and unmeasured pleiotropy. Stat Med 2017;36(29):4705-18. doi: 10.1002/sim.7492 [published Online First: 2017/09/30]4) Swerdlow DI, Preiss D, Kuchenbaecker KB, et al. HMG-coenzyme A reductase inhibition, type 2 diabetes, and bodyweight: evidence from genetic analysis and randomised trials. Lancet (London, England) 2015;385(9965):351-61. doi: 10.1016/s0140-6736(14)61183-1 [published Online First: 2014/09/30]5) Larsson SC, Bäck M, Rees JMB, et al. Body mass index and body composition in relation to 14 cardiovascular conditions in UK Biobank: a Mendelian randomization study. European heart journal 2020;41(2):221-26. doi: 10.1093/eurheartj/ehz388 [published Online First: 2019/06/14]6) Eriksson J, Haring R, Grarup N, et al. Causal relationship between obesity and serum testosterone status in men: A bi-directional mendelian randomization analysis. PLoS One 2017;12(4):e0176277. doi: 10.1371/journal.pone.0176277 [published Online First: 2017/04/28]7) Hemani G, Tilling K, Davey Smith G. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS genetics 2017;13(11):e1007081. doi: 10.1371/journal.pgen.1007081 [published Online First: 2017/11/18]8) Rafiq S, Frayling TM, Murray A, et al. A common variant of the interleukin 6 receptor (IL-6r) gene increases IL-6r and IL-6 levels, without other inflammatory effects. Genes Immun 2007;8(7):552-9. doi: 10.1038/sj.gene.6364414 [published Online First: 2007/08/03]9) Ference BA, Ray KK, Catapano AL, et al. Mendelian Randomization Study of ACLY and Cardiovascular Disease. The New England journal of medicine 2019;380(11):1033-42. doi: 10.1056/NEJMoa1806747 [published Online First: 2019/03/14]10) Ruth KS, Day FR, Tyrrell J, et al. Using human genetics to understand the disease impacts of testosterone in men and women. Nat Med 2020;26(2):252-58. doi: 10.1038/s41591-020-0751-511) Sinnott-Armstrong N, Tanigawa Y, Amar D, et al. Genetics of 38 blood and urine biomarkers in the UK Biobank. 201912) Nikpay M, Goel A, Won HH, et al. A comprehensive 1,000 Genomes-based genome-wide association meta-analysis of coronary artery disease. Nature genetics 2015;47(10):1121-30. doi: 10.1038/ng.3396 [published Online First: 2015/09/08]13) Ference BA, Robinson JG, Brook RD, et al. Variation in PCSK9 and HMGCR and Risk of Cardiovascular Disease and Diabetes. The New England journal of medicine 2016;375(22):2144-53. doi: 10.1056/NEJMoa1604304 [published Online First: 2016/12/14]14) Hernán MA, Hernández-Díaz S, Robins JM. A structural approach to selection bias. Epidemiology 2004;15(5):615-25. doi: 10.1097/01.ede.0000135174.63482.43 [published Online First: 2004/08/17]15) Schooling CM, Au Yeung SL, Freeman G, et al. The effect of statins on testosterone in men and women, a systematic review and meta-analysis of randomized controlled trials. BMC Med 2013;11:57. doi: 10.1186/1741-7015-11-57 [published Online First: 2013/03/02]16) Schooling CM, Au Yeung SL, Leung GM. Why do statins reduce cardiovascular disease more than other lipid modulating therapies? Eur J Clin Invest 2014;44(11):1135-40. doi: 10.1111/eci.12342 [published Online First: 2014/09/25]17) Manson JE, Chlebowski RT, Stefanick ML, et al. Menopausal hormone therapy and health outcomes during the intervention and extended poststopping phases of the Women's Health Initiative randomized trials. Jama 2013;310(13):1353-68. doi: 10.1001/jama.2013.278040 [published Online First: 2013/10/03]18) The Coronary Drug Project. Findings leading to discontinuation of the 2.5-mg day estrogen group. The coronary Drug Project Research Group. Jama 1973;226(6):652-7. [published Online First: 1973/11/05]19) Willer CJ, Schmidt EM, Sengupta S, et al. Discovery and refinement of loci associated with lipid levels. Nature genetics 2013;45(11):1274-83. doi: 10.1038/ng.2797 [published Online First: 2013/10/08]et al.

Author details

1. CM Schooling

   1. School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong, China
   2. City University of New York, Graduate School of Public Health and Health Policy, New York, United States

   Contribution

   Conceptualization, Formal analysis, Supervision, Validation, Methodology, Writing - original draft, Project administration

   For correspondence

   cms1@hku.hk

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-9933-5887

2. JV Zhao

   School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong, China

   Contribution

   Data curation, Validation, Writing - review and editing

   Competing interests

   No competing interests declared

3. SL Au Yeung

   School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong, China

   Contribution

   Resources, Data curation, Supervision, Investigation, Project administration, Writing - review and editing

   Competing interests

   No competing interests declared

4. GM Leung

   School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong, China

   Contribution

   Conceptualization, Methodology, Writing - review and editing

   Competing interests

   No competing interests declared

• 1,226

  Page views

• 128

  Downloads

• 18

  Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.