Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Pharmacogenomics of statin-related myopathy: Meta-analysis of rare variants from whole-exome sequencing

  • James S. Floyd ,

    Contributed equally to this work with: James S. Floyd, Katarzyna M. Bloch, Jennifer A. Brody

    Roles Data curation, Funding acquisition, Investigation, Project administration, Validation, Writing – review & editing

    Affiliation Department of Medicine, University of Washington, Seattle, Washington, United States of America

  • Katarzyna M. Bloch ,

    Contributed equally to this work with: James S. Floyd, Katarzyna M. Bloch, Jennifer A. Brody

    Roles Formal analysis, Investigation, Project administration, Validation, Writing – original draft

    Affiliation MRC Centre for Drug Safety Science, Department of Molecular and Clinical Pharmacology, University of Liverpool, Liverpool, United Kingdom

  • Jennifer A. Brody ,

    Contributed equally to this work with: James S. Floyd, Katarzyna M. Bloch, Jennifer A. Brody

    Roles Data curation, Formal analysis, Methodology, Project administration, Validation, Visualization, Writing – review & editing

    Affiliation Department of Medicine, University of Washington, Seattle, Washington, United States of America

  • Cyrielle Maroteau,

    Roles Data curation, Formal analysis, Investigation, Software, Validation, Writing – review & editing

    Affiliation Division of Molecular and Clinical Medicine, University of Dundee, Dundee, United Kingdom

  • Moneeza K. Siddiqui,

    Roles Data curation, Investigation, Project administration, Software, Validation

    Affiliation Division of Molecular and Clinical Medicine, University of Dundee, Dundee, United Kingdom

  • Richard Gregory,

    Roles Data curation, Formal analysis, Methodology, Resources, Software, Supervision

    Affiliation Functional and Comparative Genomics, Institute of Integrative Biology, University of Liverpool, Liverpool, United Kingdom

  • Daniel F. Carr,

    Roles Conceptualization, Investigation, Methodology, Project administration, Supervision, Writing – review & editing

    Affiliation Department of Medicine, University of Washington, Seattle, Washington, United States of America

  • Mariam Molokhia,

    Roles Data curation, Funding acquisition, Project administration, Supervision, Validation

    Affiliation School of Population Health and Environmental Sciences, London, United Kingdom

  • Xiaoming Liu,

    Roles Data curation, Formal analysis, Investigation, Methodology, Resources, Software

    Affiliation Human Genetics Center, University of Texas Health Science Center, Houston, United States of America

  • Joshua C. Bis,

    Roles Data curation, Investigation, Project administration, Resources, Supervision, Visualization

    Affiliation Department of Medicine, University of Washington, Seattle, Washington, United States of America

  • Ammar Ahmed,

    Roles Data curation, Project administration, Validation, Visualization, Writing – review & editing

    Affiliation Medical School, University of Liverpool, Liverpool, United Kingdom

  • Xuan Liu,

    Roles Data curation, Formal analysis, Investigation, Methodology, Software, Validation

    Affiliation Functional and Comparative Genomics, Institute of Integrative Biology, University of Liverpool, Liverpool, United Kingdom

  • Pär Hallberg,

    Roles Data curation, Formal analysis, Funding acquisition, Investigation, Project administration, Resources, Software, Writing – review & editing

    Affiliation Department of Medical Sciences, Clinical Pharmacology and Science for Life Laboratory, Uppsala University, Uppsala University Hospital, Uppsala, Sweden

  • Qun-Ying Yue,

    Roles Data curation, Funding acquisition, Investigation, Methodology, Resources, Supervision, Writing – review & editing

    Affiliation Medical Products Agency, Uppsala, Sweden

  • Patrik K. E. Magnusson,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Project administration, Supervision, Writing – review & editing

    Affiliation Swedish Twin Registry, Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden

  • Diane Brisson,

    Roles Data curation, Formal analysis, Investigation, Methodology, Software, Supervision, Writing – review & editing

    Affiliation Clinical Lipidology and Rare Lipid Disorders Unit, Department of Medicine, Université de Montréal Community Gene Medicine Center, Lipid Clinic Chicoutimi Hospital and ECOGENE-21 Clinical and Translational Research Center, Chicoutimi, Quebec, Canada

  • Kerri L. Wiggins,

    Roles Data curation, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing

    Affiliation Department of Medicine, University of Washington, Seattle, Washington, United States of America

  • Alanna C. Morrison,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Project administration, Writing – review & editing

    Affiliation Human Genetics Center, Department of Epidemiology, Human Genetics, and Environmental Sciences, School of Public Health, The University of Texas Health Science Center at Houston, Houston, Texas, United States of America

  • Etienne Khoury,

    Roles Data curation, Formal analysis, Funding acquisition, Investigation, Software, Supervision, Writing – review & editing

    Affiliation Clinical Lipidology and Rare Lipid Disorders Unit, Department of Medicine, Université de Montréal Community Gene Medicine Center, Lipid Clinic Chicoutimi Hospital and ECOGENE-21 Clinical and Translational Research Center, Chicoutimi, Quebec, Canada

  • Paul McKeigue,

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Resources, Supervision, Writing – review & editing

    Affiliation Usher Institute of Population Health Sciences and Informatics, University of Edinburgh Medical School, Edinburgh, Scotland, United Kingdom

  • Bruno H. Stricker,

    Roles Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Writing – review & editing

    Affiliation Department of Epidemiology, Erasmus Medical Centre, Rotterdam, the Netherlands

  • Maryse Lapeyre-Mestre,

    Roles Conceptualization, Funding acquisition, Methodology, Project administration, Software, Supervision, Validation, Writing – review & editing

    Affiliation Paul Sabatier University - Toulouse III, UPS Toulouse, Laboratoire de Pharmacologie Medicale et Clinique, Toulouse, France

  • Susan R. Heckbert,

    Roles Conceptualization, Formal analysis, Funding acquisition, Project administration, Resources, Software

    Affiliation Department of Epidemiology, University of Washington, Seattle, Washington, United States of America

  • Arlene M. Gallagher,

    Roles Data curation, Formal analysis, Investigation, Methodology, Software, Validation

    Affiliation Clinical Practice Research Datalink (CPRD) Medicines and Healthcare Products Regulatory Agency, London, United Kingdom

  • Hector Chinoy,

    Roles Funding acquisition, Investigation, Methodology, Project administration, Resources, Validation

    Affiliation Rheumatology Department, Salford Royal NHS Foundation Trust, Manchester Academic Health Science Centre, Salford, United Kingdom

  • Richard A. Gibbs,

    Roles Data curation, Formal analysis, Funding acquisition, Methodology, Project administration, Supervision, Writing – review & editing

    Affiliation Human Genome Sequencing Center, Baylor College of Medicine, Houston, United States of America

  • Emmanuelle Bondon-Guitton,

    Roles Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Writing – review & editing

    Affiliation Centre Hospitalier Universitaire de Toulouse, CHU Toulouse, Centre de Pharmacovigilance, Toulouse, France

  • Russell Tracy,

    Roles Data curation, Formal analysis, Funding acquisition, Investigation, Project administration, Resources, Writing – review & editing

    Affiliation Departments of Pathology & Laboratory Medicine and Biochemistry, Larner College of Medicine, University of Vermont, Burlington, Vermont, United States of America

  • Eric Boerwinkle,

    Roles Conceptualization, Data curation, Funding acquisition, Methodology, Supervision, Writing – review & editing

    Affiliation Human Genetics Center, Department of Epidemiology, Human Genetics, and Environmental Sciences, School of Public Health, The University of Texas Health Science Center at Houston, Houston, Texas, United States of America

  • Daniel Gaudet,

    Roles Conceptualization, Data curation, Investigation, Methodology, Software, Writing – review & editing

    Affiliation Clinical Lipidology and Rare Lipid Disorders Unit, Department of Medicine, Université de Montréal Community Gene Medicine Center, Lipid Clinic Chicoutimi Hospital and ECOGENE-21 Clinical and Translational Research Center, Chicoutimi, Quebec, Canada

  • Anita Conforti,

    Roles Data curation, Formal analysis, Investigation, Methodology, Supervision, Validation

    Affiliation U.O. Farmacologia, Policlinico "Gb Rossi", Verona, Italy

  • Tjeerd van Staa,

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Resources, Software, Supervision, Writing – review & editing

    Affiliation Division of Informatics, Imaging & Data Sciences, University of Manchester, Manchester, United Kingdom

  • Colleen M. Sitlani,

    Roles Data curation, Formal analysis, Methodology, Validation, Writing – review & editing

    Affiliation Department of Medicine, University of Washington, Seattle, Washington, United States of America

  • Kenneth M. Rice,

    Roles Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Software, Supervision, Writing – review & editing

    Affiliation Department of Biostatistics, University of Washington, Seattle, Washington, United States of America

  • Anke-Hilse Maitland-van der Zee,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Project administration, Supervision, Validation, Writing – review & editing

    Affiliation Respiratory Medicine/Pediatric Respiratory Medicine, AMC, Amsterdam, the Netherlands

  • Mia Wadelius,

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – review & editing

    Affiliation Medical School, University of Liverpool, Liverpool, United Kingdom

  • Andrew P. Morris,

    Roles Conceptualization, Formal analysis, Methodology, Supervision, Validation, Visualization, Writing – original draft

    Affiliation Department of Biostatistics, University of Liverpool, Liverpool, United Kingdom

  • Munir Pirmohamed,

    Roles Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing

    Affiliation MRC Centre for Drug Safety Science, Department of Molecular and Clinical Pharmacology, University of Liverpool, Liverpool, United Kingdom

  • Colin A. N. Palmer,

    Roles Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Writing – review & editing

    Affiliation Division of Molecular and Clinical Medicine, University of Dundee, Dundee, United Kingdom

  • Bruce M. Psaty ,

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Resources, Supervision, Writing – review & editing

    ‡ These authors also contributed equally to this work.

    Affiliation Department of Medicine, University of Washington, Seattle, Washington, United States of America

  • Ana Alfirevic ,

    Roles Conceptualization, Data curation, Funding acquisition, Investigation, Supervision, Validation, Writing – original draft, Writing – review & editing

    Ana.Alfirevic@liverpool.ac.uk

    ‡ These authors also contributed equally to this work.

    Affiliation MRC Centre for Drug Safety Science, Department of Molecular and Clinical Pharmacology, University of Liverpool, Liverpool, United Kingdom

  •  [ ... ],
  • on behalf of the PREDICTION-ADR Consortium and EUDRAGENE

    The complete membership of the author groups can be found in the Acknowledgments

  • [ view all ]
  • [ view less ]

Abstract

Aims

Statin-related myopathy (SRM), which includes rhabdomyolysis, is an uncommon but important adverse drug reaction because the number of people prescribed statins world-wide is large. Previous association studies of common genetic variants have had limited success in identifying a genetic basis for this adverse drug reaction. We conducted a multi-site whole-exome sequencing study to investigate whether rare coding variants confer an increased risk of SRM.

Methods and results

SRM 3–5 cases (N = 505) and statin treatment-tolerant controls (N = 2047) were recruited from multiple sites in North America and Europe. SRM 3–5 was defined as symptoms consistent with muscle injury and an elevated creatine phosphokinase level >4 times upper limit of normal without another likely cause of muscle injury. Whole-exome sequencing and variant calling was coordinated from two analysis centres, and results of single-variant and gene-based burden tests were meta-analysed. No genome-wide significant associations were identified. Given the large number of cases, we had 80% power to identify a variant with minor allele frequency of 0.01 that increases the risk of SRM 6-fold at genome-wide significance.

Conclusions

In this large whole-exome sequencing study of severe statin-related muscle injury conducted to date, we did not find evidence that rare coding variants are responsible for this adverse drug reaction. Larger sample sizes would be required to identify rare variants with small effects, but it is unclear whether such findings would be clinically actionable.

Introduction

Lipid-lowering drugs that inhibit 5-hydroxy-3-methylglutaryl-coenzyme A (HMG-CoA), known as statins, are widely used for the primary and secondary prevention of cardiovascular disease (CVD). All statins can cause muscle toxicity that ranges in severity from mild muscle pain to severe myopathy or rhabdomyolysis, which can lead to kidney failure and death [1]. The most severe forms of statin-related muscle injury are uncommon, suggesting a genetic predisposition. For instance, severe statin-related myopathy (SRM) restricted to cases with creatine phosphokinase (CK) levels >10x the upper limit of normal (ULN) (SRM 4) requiring hospitalisation, occurs in approximately 1–10 out of every 10,000 individuals taking standard statin doses [2].

Previous genome-wide association studies (GWAS) have identified a common variant in the drug transporter gene SLCO1B1 that increases the risk of SRM defined by varying CK elevation criteria, 2- to 4-fold through a pharmacokinetic mechanism, but the underlying causes of this adverse drug reaction remain largely unexplained [35]. To investigate whether rare variants in coding regions of the genome may cause SRM, we conducted whole-exome sequencing on cases of SRM identified from multiple sites throughout North America and Europe, harmonized variant calling and analysis at two centres, and meta-analysed results from single-variant and gene-based burden tests across two large case-control studies.

Materials and methods

Study population and overall study design

Identification and recruitment of case and control subjects, whole-exome sequencing, variant calling, and statistical analyses were coordinated by two analysis centers: one based at the University of Washington that included 291 cases and 1540 controls from North America and the United Kingdom (“US-UK”) and one based in Liverpool that included 214 cases and 507 controls from Europe (“PREDICTION-ADR”). Within each of these two large case-control studies, participants were recruited from multiple sites. Because there were few non-European ancestry participants in either case-control study, and many sites had none, we restricted analyses to participants of European ancestry. The details of recruitment at each site are described in the Supporting Information files. Whole-exome sequencing and variant calling were harmonized within each of the two case-control studies, and investigators from both studies agreed prospectively to implement a coordinated statistical analysis plan, which included protocols for meta-analysis of results from both single-variant and gene-based burden tests. The University of Washington Human Subjects Division, East of Scotland Research Ethics Service REC1 and REC2, North West–Liverpool Central Research Ethics Committee, Regional Ethical Review Boards in Uppsala (2008/213 and 2010/231) and Stockholm (2007/644-31 and 2011/463-32) and the Medical Ethics Committee from the University Medical Center in Utrecht reviewed and approved this study.

Case definition

Using a previously published case definition (S1 Table), SRM 3–5 cases received statin treatment at the time of onset of muscle injury symptoms, had a peak CK level > 4x ULN (SRM 3) or supporting clinical evidence of rhabdomyolysis (SRM 5), and had no other likely cause of muscle injury other than statin treatment.[1] In each study, non-statin causes of muscle injury were excluded based on a review of medical records and/or interview of patients by physicians or study staff. Severe SRM cases (SRM 4) sometimes called severe myopathy and SRM 5, sometimes called rhabdomyolysis in some previous studies [6], met these criteria and had a peak CK level > 10x ULN (S1 Table). SRM cases with myopathy and evidence of anti-3-hydroxy-3-methyl-glutaryl-coenzyme A reductase (HMGCR) autoantibodies (SRM 6), which are suggestive of autoimmune myopathy, were excluded at all sites [7]. In both case-control studies, statin treatment-tolerant controls subjects were users of statins who had no evidence of muscle injury during follow up.

Whole-exome sequencing, variant calling, quality control

Exome sequencing, variant calling, and quality control of sequence data were performed separately in each case-control study (S1 Appendix and S2 Table). Each study conducted sequencing and called variants for cases and controls together, to eliminate batch effects across case-control sampling. We created a combined variant annotation file including all quality-controlled variant sites observed in either study. Variants were annotated using ANNOVAR [8] and dbNSFP v3.0 [9] according to the reference genome GRCh37.

Statistical analysis

Within each case-control study, which included only European ancestry participants, population structure was adjusted for using principal components derived from a genetic relationship matrix. To estimate associations between single variants and the risk of SRM, the Firth bias-corrected score test was used [10]. The primary single-variant analysis included all SRM cases and all control subjects. Single variant analyses were restricted to variants with at least 10 copies of the minor allele, but not further restricted by annotated variant function. Results from the two case-control studies were combined using a Z-based meta-analysis [11]. The following secondary analyses were also conducted: restriction to users of simvastatin, atorvastatin, and cerivastatin (the three most commonly-used statins in these studies); stratification of SRM cases into moderate (SRM 3: CK ≥ 4x ULN and <10x ULN) and severe (SRM 4 and 5: CK ≥ 10x ULN) categories; and inclusion of all missense variants. In the US-UK case-control study, information on fibrate use was collected. Because fibrates can cause drug-drug interactions with statins through inhibition of the metabolism of statins, and because strong drug-drug interactions can mask the presence of a genuine drug-gene interaction that may be less potent, we conducted secondary analyses restricted to non-users of fibrates.

Rare variants (MAF < 1%) were collapsed into gene-burden scores and tested for association with the risk of SRM [12]. Gene-burden results from the two case-control studies were combined using a Z-based meta-analysis. The primary gene-based analysis included variants with MAF < 1% that are missense, bioinformatically predicted to be damaging (MetaLR > .5) [13], stop-gain/loss, coding indels, or splice site variants. Secondary gene-burden analyses restricted to users of simvastatin or atorvastatin, stratified SRM cases into moderate (SRM 3: CK > = 4x ULN and <10x ULN) and severe (SRM 4 and 5: CK > = 10x ULN) categories, and included rare missense variants regardless of their prediction score for function.

For each analysis, the Bonferroni-corrected threshold for statistical significance was P = 0.05 / number of single variants or genes tested. For the primary single-variant analysis, we conducted post-hoc power calculations using QUANTO 1.1 [14, 15] to determine study power in the primary analysis population across a range of effect sizes and MAFs, based on α = 0.05 / number of single variants included in the primary analysis and an estimated disease probability among control subjects of 0.0001 [16].

The full results from each meta-analysis will be available on dbGaP within the next release of Cohort for Heart and Aging Research in Genomic Epidemiology Results CHARGE Consortium summary results (accession phs000930) (http://www.chargeconsortium.com/main/results). Information on how to access individual-level study data is provided in the Supporting Information files.

Results

Across the two case-control studies, there were 505 SRM cases and 2,047 treatment-tolerant controls that passed sample-level quality control (Table 1). Most of the US-UK cases (mean CK 159x ULN) were severe (65%) and used cerivastatin or simvastatin, while most of the PREDICTION-ADR cases (mean CK 32x ULN) were moderate (61%) and used simvastatin or atorvastatin.

thumbnail
Table 1. Characteristics of statin-related myopathy cases and statin-tolerant controls from two international case-control studies.

https://doi.org/10.1371/journal.pone.0218115.t001

Single-variant results

The mean sequencing depth in US-UK case-control study was 78x, and in PREDICTION-ADR it was 56x. After restricting to rare (MAF < 0.01) likely damaging variants, 162,813 variants passed quality control in one or both studies and were included in meta-analyses (Fig 1). In primary analyses that included all SRM cases and all statin-tolerant controls, there were no variants that met the Bonferroni-corrected threshold for statistical significance (P < 3.07 x 10−7) (Table 2). The widely replicated locus for statin-related muscle injury in SLCO1B1 (rs4149056) was associated with a 1.59-fold increased risk of SRM (95% confidence interval, 1.41–1.77), but this finding was not significant after correction for the number of single variants tested.

thumbnail
Fig 1. Single variant analysis.

Single variants included in primarily analysis by frequency (panel A) and type (panel B).

https://doi.org/10.1371/journal.pone.0218115.g001

thumbnail
Table 2. Top associations from meta-analysis of primary single-variant results from US-UK and PREDICTION-ADR case-control studies of statin-related myopathy.

https://doi.org/10.1371/journal.pone.0218115.t002

Secondary analyses restricted by statin type and stratified by moderate or severe cases of SRM did not yield any associations that met criteria for statistical significance (S3 Table). Among non-users of fibrates in the US-UK case-control study, the SLCO1B1 variant rs4149056 was associated with a 4.01-fold increased risk of SRM (95% confidence interval, 2.61–6.17), and this association met criteria for statistical significance (P = 5.46 x 10−11) (S4 Table).

Gene-burden results

In total, 10,244 genes passed quality control thresholds in one or both studies and were included in the meta-analyses of burden test results. In primary analyses of all SRM cases and all controls, no genes met the Bonferroni-corrected threshold for statistical significance (P < 4.9 x 10−6) (Table 3). Similarly, none of the secondary burden tests analyses yielded genome-wide significant findings (S5 Table).

thumbnail
Table 3. Top associations from meta-analysis of primary gene burden results from US-UK and PREDICTION-ADR case-control studies of statin-related myopathy.

https://doi.org/10.1371/journal.pone.0218115.t003

Post-hoc power calculations

We conducted post-hoc power calculations to estimate study power in the primary analysis population to identify single-variants associated with SRM given a range of MAFs and effect sizes (Fig 2). For a variant with a MAF of 0.01, we had 98% power to detect an OR of 5, and for a variant with a MAF of 0.005, we had 84% power to detect an OR of 6.

thumbnail
Fig 2. Power calculation.

Post-hoc power calculations to identify a single variant associated with statin-related myopathy in primary analysis population given effect size (odds ratio) and minor allele frequency (MAF), α = 0.05 / 162,813.

https://doi.org/10.1371/journal.pone.0218115.g002

Discussion

We recruited more than 500 patients with statin-related muscle injury from multiple sites in North America and Europe to conduct the largest whole-exome sequencing study of this adverse drug reaction to date. Using high-throughput sequencing technology, joint variant calling at two experienced analysis centres, and rigorous statistical methods, we did not identify any novel coding variants associated with SRM 3–5. We also evaluated hypotheses about the burden of rare coding variants within genes, muscle injury due to specific statin drugs, and the severity of muscle injury; the findings from these analyses were null as well.

In 2008, Link et al. found that a common nonsynonymous variant in the drug transporter SLCO1B1 (rs4149056) was associated with a 4-fold increase in the risk of myopathy among users of high dose simvastatin at genome-wide levels of statistical significance [4]. This finding has been widely replicated in other settings for other statin drugs [3, 5, 1719], and in our study this common coding variant was associated with a 4-fold increased risk of SRM 5 among non-users of fibrates. Interestingly, we did not find an association between the SLCO1B1 polymorphisms and mild forms of SRM. Attempts to discover additional common variants related to statin-related muscle injury have failed to yield replicable findings [2022]. A recent study from the Czech Republic conducted whole-exome sequencing on 88 patients with mild statin-associated muscle toxicity, reporting an increased burden of rare variation in 24 genes [23]. We did not replicate these findings in our study, perhaps due to different phenotypes as most of our patients have more severe SRM. The validity of these findings may be limited because of the use of publicly available sequencing repositories that do not have information on medication use as the referent group rather than statin-tolerant control subjects.

Strengths of our study include the recruitment of patients from multiple international sites; the use of a standardized case-definition for different forms of myopathy in our study (SRM 3–5) with the mildest being CK>4x ULN (SRM 3) and rhabdomyolysis ranging from CK>10 to >40xULN with end-organ impairment and muscle symptoms[1], the depth of whole-exome sequencing in both case-control studies; the use of statin exposed controls to minimize genetic associations with indications for statin therapy, the coordinated statistical analysis plan with correction for multiple testing to reduce the risk of false positive findings; and the prospective plan for a meta-analysis of results across two large-case control studies.

There were several limitations. As with any adverse drug reaction, it is difficult to attribute causality to a specific drug, even with prospectively developed case definitions and expert adjudication of cases. In addition, the presence of unknown drug-drug interactions with statin that increase the risk of SRM may have masked our ability to identify genuine drug-gene interactions. We recognize that our inability to recruit a large number of cases from diverse ancestries limits the generalizability of our findings to variants observed in European ancestry participants. Half of the SRM cases in the US-UK case-control study were due to the drug cerivastatin, and we were unable to identify a large number of controls subjects with WES who used this drug, since it was removed from the market in 2001. If there is heterogeneity in the effects of genetic variants on the risk of SRM by the type of statin drug or severity of SRM, this might have reduced our overall power to identify genetic loci for this adverse drug reaction in our primary analysis, which included all samples. Finally, using our stringent correction for multiple comparisons, we may have missed an association with rare, loss of function variants in biologically plausible metabolising enzyme or transporter genes that could have led to large increase in drug exposure and therefore muscle damage in some individuals [24].

While our post-hoc calculations demonstrated that our study was well-powered to identify rare coding variants (MAF 1%) with large effect sizes (OR 5), we had limited power to detect rare variants with small- or moderate-sized effects, or to detect even rarer variants with large effects. For instance, to have 90% power to identify a coding variant with a frequency of 0.001 among controls that increases the risk of SRM 6-fold, we would have to recruit and conduct whole-exome sequencing on approximately 2,700 cases with this rare adverse drug reaction, matched with 10,800 treatment-tolerant controls.

It is possible that rare variants exist that increase the risk of SRM, perhaps in non-coding regions of the genome, but they are unlikely to be detected with the sample sizes available in ongoing studies. Based on the null findings from our study, genotyping rare coding alleles for the prediction of severe statin-related muscle injury seems unlikely to yield clinically actionable findings.

Supporting information

S1 Appendix. Recruitment and study populations; genotyping and variant calling.

https://doi.org/10.1371/journal.pone.0218115.s001

(DOCX)

S1 Table. Definition of statin-related myopathy.

https://doi.org/10.1371/journal.pone.0218115.s002

(DOCX)

S2 Table. Whole exome sequence metrics, PREDICTION-ADR.

https://doi.org/10.1371/journal.pone.0218115.s003

(DOCX)

S3 Table. Top hits from single variant meta-analysis.

Primary analysis including all samples and analyses stratified by statin type and severity of statin-related myopathy.

https://doi.org/10.1371/journal.pone.0218115.s004

(XLSX)

S4 Table. Top hits from single variant meta-analysis, non-users of fibrates in the US-UK case-control study.

SLCO1B1 variant rs4149056 (12:21331549:T:C) was associated with a 4.01-fold increased risk of SRM (95% CI, 2.61–6.17), and this association met criteria for statistical significance (P = 5.46 x 10−11).

https://doi.org/10.1371/journal.pone.0218115.s005

(XLSX)

S5 Table. Top hits from gene burden meta-analysis.

Primary analysis including all samples and secondary analyses stratified by statin type and severity of statin-related myopathy.

https://doi.org/10.1371/journal.pone.0218115.s006

(XLSX)

S1 Fig. QQ plots of meta-analysed single-variant results from US-UK and PREDICTION-ADR case-control studies for primary and secondary analyses.

https://doi.org/10.1371/journal.pone.0218115.s007

(DOCX)

S2 Fig. QQ plots of meta-analysed gene burden test results from US-UK and PREDICTION-ADR case-control studies for primary and secondary analyses.

Columns: (1) predicted damaging variants, (2) all loss of function, nonsynonymous and missense. Rows: (1) all SRM, (2) simvastatin only, (3) atorvastatin only, (4) mild SRM, (5) severe SRM.

https://doi.org/10.1371/journal.pone.0218115.s008

(PDF)

Acknowledgments

We thank research nurses Anita Hanson, Ruth Westhead and Diane Van Eker, and administrators Joe Fahy, Emma McDonough and Angela Foxcroft in Liverpool, for recruiting patients and administering the clinical database. We thank the Centre for Genomic Research lead by Prof Christiane Hertz-Fowler, which is the MRC sequencing hub at the University of Liverpool, UK.

We thank research nurses Ulrica Ramqvist and Elisabeth Stjernberg and research assistants Sofie Collin, Eva Prado Lopez, Agnes Kataja Knight, Agnes Wadelius, and Martha Wadelius, Department of Medical Sciences, Clinical Pharmacology, Uppsala University, Uppsala, Sweden, for recruiting and interviewing cases and for administering the phenotype database. We acknowledge Barbro Sandin at the Swedish Twin Registry for access to data for the Swedish controls.

The lead author from PREDICTION-ADR is Prof Colin AN Palmer (c.n.a.palmer@dundee.ac.uk),

PREDICTION ADR consortium members: Colin AN Palmer, Ana Alfirevic, Mia Wadelius, Anke-Hilse Maitland-van der Zee, Ekaterina V Baranova, Niclas Eriksson, Cyrielle Maroteau, Katarzyna M Bloch, Moneeza K Siddiqui, Daniel F Carr, Pär Hallberg, Munir Pirmohamed, Anu Aaspollu, Alun McCarthy, Olivier Delrieu.

The lead author from EUDRAGENE is Dr Molokhia (mariam.molokhia@kcl.ac.uk),

EUDRAGEN consortium members: Mariam Molokhia, Paul McKeigue, Bruno H Stricker, Maryse Lapeyre-Mestre, Anita Conforti, Emmanuelle Bondon-Guitton.

References

  1. 1. Alfirevic A, Neely D, Armitage J, Chinoy H, Cooper RG, Laaksonen R, et al. Phenotype standardization for statin-induced myotoxicity. Clin Pharmacol Ther. 2014;96(4):470–6. pmid:24897241.
  2. 2. Graham DJ, Staffa JA, Shatin D, Andrade SE, Schech SD, La Grenade L, et al. Incidence of hospitalized rhabdomyolysis in patients treated with lipid-lowering drugs. JAMA. 2004;292(21):2585–90. pmid:15572716.
  3. 3. Carr DF, O’Meara H, Jorgensen AL, Campbell J, Hobbs M, McCann G, et al. SLCO1B1 Genetic Variant Associated With Statin-Induced Myopathy: A Proof-of-Concept Study Using the Clinical Practice Research Datalink. Clin Pharmacol Ther. 2013;94(6):695–701. pmid:23942138.
  4. 4. Link E, Parish S, Armitage J, Bowman L, Heath S, Matsuda F, et al. SLCO1B1 variants and statin-induced myopathy—a genomewide study. The New England journal of medicine. 2008;359(8):789–99. Epub 2008/07/25. pmid:18650507.
  5. 5. Marciante KD, Durda JP, Heckbert SR, Lumley T, Rice K, McKnight B, et al. Cerivastatin, genetic variants, and the risk of rhabdomyolysis. Pharmacogenetics and genomics. 2011;21(5):280–8. pmid:21386754.
  6. 6. Floyd JS, Kaspera R, Marciante KD, Weiss NS, Heckbert SR, Lumley T, et al. A screening study of drug-drug interactions in cerivastatin users: an adverse effect of clopidogrel. Clin Pharmacol Ther. 2012;91(5):896–904. pmid:22419147.
  7. 7. Mammen AL, Chung T, Christopher-Stine L, Rosen P, Rosen A, Doering KR, et al. Autoantibodies against 3-hydroxy-3-methylglutaryl-coenzyme A reductase in patients with statin-associated autoimmune myopathy. Arthritis Rheum. 2011;63(3):713–21. pmid:21360500.
  8. 8. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16):e164. pmid:20601685.
  9. 9. Liu X, Wu C, Li C, Boerwinkle E. dbNSFP v3.0: A One-Stop Database of Functional Predictions and Annotations for Human Nonsynonymous and Splice-Site SNVs. Human mutation. 2016;37(3):235–41. pmid:26555599.
  10. 10. Ma C, Blackwell T, Boehnke M, Scott LJ, Go TDi. Recommended joint and meta-analysis strategies for case-control association testing of single low-count variants. Genetic epidemiology. 2013;37(6):539–50. pmid:23788246.
  11. 11. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26(17):2190–1. Epub 2010/07/10. pmid:20616382.
  12. 12. Li B, Leal SM. Methods for detecting associations with rare variants for common diseases: application to analysis of sequence data. Am J Hum Genet. 2008;83(3):311–21. pmid:18691683.
  13. 13. Dong C, Wei P, Jian X, Gibbs R, Boerwinkle E, Wang K, et al. Comparison and integration of deleteriousness prediction methods for nonsynonymous SNVs in whole exome sequencing studies. Hum Mol Genet. 2015;24(8):2125–37. pmid:25552646.
  14. 14. Gauderman WJ. Sample size requirements for matched case-control studies of gene-environment interaction. Stat Med. 2002;21(1):35–50. Epub 2002/01/10. pmid:11782049.
  15. 15. Gauderman WJ, Morrison JM. QUANTO 1.1: A computer program for power and sample size calculations for genetic-epidemiology studies, http://hydra.usc.edu/gxe, 2006.
  16. 16. Graham DJ, Staffa JA, Shatin D, Andrade SE, Schech SD, La Grenade L, et al. Incidence of hospitalized rhabdomyolysis in patients treated with lipid-lowering drugs. JAMA: the journal of the American Medical Association. 2004;292(21):2585–90. Epub 2004/12/02. pmid:15572716.
  17. 17. Brunham LR, Lansberg PJ, Zhang L, Miao F, Carter C, Hovingh GK, et al. Differential effect of the rs4149056 variant in SLCO1B1 on myopathy associated with simvastatin and atorvastatin. The pharmacogenomics journal. 2012;12(3):233–7. pmid:21243006.
  18. 18. Linde R, Peng L, Desai M, Feldman D. The role of vitamin D and SLCO1B1*5 gene polymorphism in statin-associated myalgias. Dermato-endocrinology. 2010;2(2):77–84. pmid:21547103.
  19. 19. Voora D, Shah SH, Spasojevic I, Ali S, Reed CR, Salisbury BA, et al. The SLCO1B1*5 genetic variant is associated with statin-induced side effects. Journal of the American College of Cardiology. 2009;54(17):1609–16. pmid:19833260.
  20. 20. Carr DF, Alfirevic A, Johnson R, Chinoy H, van Staa T, Pirmohamed M. GATM gene variants and statin myopathy risk. Nature. 2014;513(7518):E1. Epub 2014/09/19. pmid:25230669.
  21. 21. Floyd JS, Bis JC, Brody JA, Heckbert SR, Rice K, Psaty BM. GATM locus does not replicate in rhabdomyolysis study. Nature. 2014;513(7518):E1–3. Epub 2014/09/19. pmid:25230668.
  22. 22. Mangravite LM, Engelhardt BE, Medina MW, Smith JD, Brown CD, Chasman DI, et al. A statin-dependent QTL for GATM expression is associated with statin-induced myopathy. Nature. 2013;502(7471):377–80. Epub 2013/09/03. pmid:23995691.
  23. 23. Neroldova M, Stranecky V, Hodanova K, Hartmannova H, Piherova L, Pristoupilova A, et al. Rare variants in known and novel candidate genes predisposing to statin-associated myopathy. Pharmacogenomics. 2016;17(13):1405–14. pmid:27296017.
  24. 24. Teslovich TM, Musunuru K, Smith AV, Edmondson AC, Stylianou IM, Koseki M, et al. Biological, clinical and population relevance of 95 loci for blood lipids. Nature. 2010;466(7307):707–13. pmid:20686565.