|Home | About | Journals | Submit | Contact Us | Français|
Skeletal muscle toxicity is the primary adverse effect for statins. In this review, we summarize current knowledge regarding the genetic and non-genetic determinants of risk for statin induced myopathy. Many genetic factors were initially identified through candidate-gene association studies (CGAS) limited to pharmacokinetic (PK) targets. Through genome-wide association studies (GWAS), it has become clear that SLCO1B1 is among the strongest PK predictors of myopathy risk. GWAS have also expanded our understanding of pharmacodynamic (PD) candidate genes, including RYR2. It is anticipated that deep re-sequencing efforts will define new loci with rare variants that contribute as well, and sophisticated computational approaches will be needed to characterize gene-gene (GxG) and gene-environment (GxE) interactions. Beyond environment, race is a critical covariate, and its influence is only partly explained by geographic differences in the frequency of known PD and PK variants. As such, admixture analyses will be essential to a full understanding of statin-induced myopathy.
Statins are among the most commonly prescribed drugs in the industrialized world . They inhibit HMG Coenzyme A reductase (HMGCR), the rate-limiting enzyme in cholesterol biosynthesis (Figure 1) [2, 3], and they reduce the frequency of coronary artery disease in patients at risk . There are seven statins currently available within the class. The first statin to be approved by the U.S. Food and Drug Administration was lovastatin in 1987, followed by simvastatin, 1988; pravastatin, 1991; fluvastatin, 1994; atorvastatin, 1997; rosuvastatin, 2003; and pitavastatin, 2009 [5, 6]. One additional agent, cerivastatin, was released (in 1998) and subsequently withdrawn from the market due to a markedly increased frequency of muscle toxicity. Between 1998 and 2001, forty reported cases of muscle toxicity due to cerivastatin use were reported to be fatal .
In general, statins are considered safe . Statin-related adverse drug reactions (ADRs) include neurocognitive complaints, musculoskeletal complaints, and very rarely nephrotoxicity or hepatotoxicity . Musculoskeletal problems represent the most common form of statin intolerance. Over the past two decades, statin use has expanded dramatically, and statin-induced muscle toxicity is becoming more fully characterized [6, 10–12]. The clinical presentation of this ADR varies widely, from mild myalgias (focal or diffuse) to rhabdomyolysis (severe skeletal muscle damage accompanied by acute kidney injury) . In the context of clinical practice, statin-induced muscle toxicity is diagnosed based upon two variables: muscle pain and circulating levels of a relatively non-specific muscle enzyme, creatine kinase (CK). CK level is often used as a marker for severity for statin-related muscle damage [9, 14]. Because CK levels can occasionally be elevated in the absence of pain, accurate risk prediction models are needed for preventing rhabdomyolysis .
The true frequency of statin-related muscle damage is difficult to quantify. Randomized clinical trials often underestimate the frequency of this ADR because patients with symptoms of intolerance are typically excluded during the run-in period [16–18]. Frequency estimates derived from databases maintained by regulatory agencies also tend to underestimate the problem because such event-reporting is voluntary . Thus, large observational databases linked to electronic medical records may represent the most accurate way of quantifying statin-induced muscle toxicity within the community . Mild myalgias related to statin use are quite common, reported to occur at frequencies ranging from 1% to 10% . Myalgias accompanied by elevation in serum CK level occur at a much lower frequency [21, 22]. The Health Improvement Network (THIN) and MediPlus databases report an annual incidence approaching 700 per million exposed per year, for intermediate myotoxicity . Rhabdomyolysis represents the most severe and potentially lethal form of this ADR [9, 14], and many clinicians consider an elevation in serum CK level >50-fold upper limit of normal (~10,000 Units per liter) both necessary and sufficient to make the diagnosis. Graham and colleagues surveyed more than 250,000 statin-exposed patients, and reported rhabdomyolysis rates of 0.000044 events per person-year . Similar rates have been observed for more than 100,000 first-time statin users followed in the United Kingdom over a course of 20 months . In general, the frequency of severe muscle damage is low in the context of monotherapy (i.e., when statins are used in the absence of other medications).
Risk and severity increase, however, in the presence of co-medications known to influence statin disposition [9, 22, 26]. Examples include medications altering the expression or activity of CYP3A4/5, UGT1A1, ABCB1, and OATP1B1 (gene name SLCO1B1). Each of these determinants of statin myopathy risk is discussed further below, within the context of pharmacokinetic candidate genes. Additional clinical factors influencing the severity of statin-induced myopathy include advanced age, small body mass index, female gender, metabolic co-morbidities, and vigorous physical exercise [20, 27–31] (Table 1). Age is a particularly strong contributor. By studying the medical records of more than 250,000 patients, Schech et al. reported that statin users ≥65 years of age have four times the risk of hospitalization due to rhabdomyolysis than younger statin users . Some of this difference can be attributed to age-related changes in body composition, and increased frequency of relevant co-morbidities (e.g., chronic liver disease or chronic kidney disease). Female gender also appears to increase risk, approximately 2-fold, perhaps due to smaller vascular volumes and reduced muscle mass resulting in greater tissue drug exposure per statin dose .
McClure and colleagues have shown that dose is also a very strong predictor of risk; the incidence rate for myotoxicity was roughly 10-fold greater in patients on high-dose statin therapy (i.e., defined as a dose equivalent to 40 mg of lovastatin daily or greater)  This observation has been replicated in an independent clinical practice-based cohort . From the records of nearly 2,000,000 unique individuals served by a single comprehensive system of care, 213 validated cases of statin-induced muscle toxicity were enrolled in a population-based study of genetic risk determinants (http://www.pharmgkb.org/contributors/pgrn/parc_profile.jsp). Within this observational cohort, the relationship between simvastatin dose and severity of myotoxicity was dose-dependent [20, 27, 34]. In 2011, the U.S. Food and Drug Administration (FDA) announced a new safety recommendation for high-dose simvastatin. FDA recommended that the 80 mg dose of simvastatin should be limited to patients “…who have been taking this dose for >12 months and have not experience any muscle toxicity” (www.fda.org).
The FDA has also recommended limiting the dose of some statins based upon major continental race (FDA Public Health Advisory on rosuvastatin; Media release March 2, 2005; http://www.accessdata.fda.gov/drugsatfda_docs/label/2005/21366slr005lbl.pdf). This recommendation was based upon (1) the observation that patients of Asian ancestry exhibit exhibit2-fold2-fold increase in AUC for rosuvastatin, compared to patients of European ancestry, following single dose exposure , and (2) the observation that patients of Asian ancestry have greater lipid lowering efficacy at lower doses of rosuvastatin, compared to patients of European ancestry . After reviewing the evidence from these and similar studies, the FDA concluded that Asian Americans were one of three important groups with an elevated risk/benefit ratio (the others were patients on immune suppression and patients with severe kidney failure). The FDA now recommends limiting patients of Asian ancestry to a 5 mg starting dose for rosuvastatin.
Differences in the overall genetic architecture between major continental races can cause variability within the clinical response to statins [29–32]. As noted below, some of the variability in myopathy risk is due to race-specific differences in the frequency of pharmacodynamic and pharmacokinetic candidate gene variants. However, a large fraction of the race-dependent variance in this trait remains unexplained. Failure to account for population confounding in the characterization of statin-induced myopathy may therefore lead to errors in inferring true association. Practice-based cohorts of diverse ancestry will be essential to our full understanding of myopathy risk in the general population [22, 37].
Within the general community, population confounding is introduced by groups of individuals that have some degree of reproductive isolation from the rest of the population, and for which allele frequencies are likely to be different from the population as a whole . Any given study may therefore be subjected to population confounding due to a) population stratification, or b) admixture. Stratification occurs when cases and controls are unintentionally drawn from two or more population subgroups. Polymorphisms that genetically mark the high-risk subgroup (i.e., found by chance at a higher frequency in that subgroup) can erroneously appear to be associated with the trait . Further, gene variants occurring at a very low frequency are typically not shared across divergent populations because they have either arisen relatively recently or because their frequencies have been influenced by population history (e.g., the out-of-Africa expansion or natural selection). Admixture is a common form of gene flow between populations. It refers to the process in which two or more genetically and phenotypically diverse populations with different allele frequencies mate and form a new, mixed or ‘hybrid’ population . A classic example of an admixed population in humans is the African-American population. As a result of the genetic admixture, the African-American population contains stretches of DNA as large as 20–30 centiMorgan (cM) that resemble mosaics of chromosomal segments, or ancestral blocks.
Markers with large allele frequency differences between ancestral populations, known as ancestry informative markers (AIMs), can be used to control for population confounding arising due to variations in background ancestry . It will therefore be critical to adequately control for population stratification and admixture when analyzing pharmacogenetic variants . Although three different methods (genomic control, structured association and principle component analysis) are used to correct for confounding, good epidemiology study design is the most efficient way to avoid confounding in population-based cohorts. No amount of adjustment can correct for poor study design where population subgroups for cases and controls do not overlap.
The degree to which genetic factors contribute to inter-individual variability in myopathy risk has been an active area of investigation for more than a decade. Much of our initial understanding came from candidate gene association studies (CGAS), particularly within the context of pharmacokinetic (PK) targets. In general, candidate gene studies are based on prior knowledge . Because statin disposition varies drug-by-drug, reflecting subtle differences in uptake, oxidation, conjugation, and efflux, models used to predict the impact of PK genes on myopathy risk require flexibility.
The clinical severity of statin-induced muscle toxicity is clearly influenced by variability in enzymes modulating statin disposition (absorption, distribution, metabolism and elimination, ADME) (Figure 2) (Table 2) . While many statins undergo phase I oxidation (atorvastatin, fluvastatin, lovastatin, simvastatin), the impact of phase I oxidation on others (pitavastatin, pravastatin, rosuvastatin) is very limited . Atorvastatin and lovastatin are oxidized primarily by the cytochrome P450 (CYP) 3A4 and 3A5 enzymes. Although the same enzymes are known to be responsible for the metabolism of simvastatin, fluvastatin, and cerivastatin, the oxidation of fluvastatin (and possibly pitavastatin) is influenced by CYP2C9 , whereas both simvastatin and cerivastatin metabolites are further oxidized by CYP2C8 [43–46]. Although controversial, the oxidation of simvastatin metabolites may also be influenced by the highly polymorphic CYP2D6 [45, 47–50]. Each of these genes - CY3A4/5, CYP2C8/9, and CYP2D6 – is polymorphic, and variability in phase I drug metabolizing enzyme genes might therefore account for patient-to-patient differences in muscle related ADRs.
The potential effect of CYP2D6 polymorphisms on statin intolerance was first explored in a cohort of 88 participants . Frudakis et al. demonstrated that CYP2D6*4 was associated with the frequency of statin induced muscle events (p = 0.001), independent of demographic variables . Kaspera et al. sequenced CYP2C8 in 126 rhabdomyolysis cases, and identified 12 novel single nucleotide polymorphisms (SNPs) with a potential to alter CYP2C8 enzyme function . A common splice variant in CYP3A5 has been associated with the magnitude of CK elevation by our group, specifically within the context of atorvastatin . The strength of this latter association was dependent upon the presence of concomitant medications known to interact with statins through processes other than phase I oxidation (e.g., phase II conjugation) .
Many statins and hydroxy-statin derivatives undergo further modification, through UDP-glucuronosyl transferase 1 (UGT1) - dependent processes (Figure 2) . It is therefore likely that genetic variability in the UGT1 enzyme family would contribute to myopathy risk as well. The entire family of UGT1 gene products (UGT1A1-12) is derived from the same locus. Because atorvastatin δ-lactone is associated with toxicity, Riedmaier et al. studied the role of UGTs in atorvastatin lactonization . After analyzing 150 human liver samples, they showed that atorvastatin lactonization is associated with both UGT1A3 immunoreactive protein levels and mRNA levels. Genetic analysis UGT1A3 mRNA and protein levels are altered by the UGT1A3*2 allele, a variant also shown to influence the rate of atorvastatin lactonization. Interestingly, expression level of UGT1A3 mRNA was also positively influenced by the well-defined UGT1A1 variant allele – UGT1A1*28 (p<0.001). This variant has previously been associated with clinical outcome within the context of a number of drug classes .
Beyond phase I (oxidative) and phase II (conjugative) statin metabolism, variability in membrane transport also contributes strongly to myopathy risk. The organic anion transporting polypeptide OATP-1B1 (gene name SLCO1B1) is expressed on the sinusoidal membrane of human hepatocytes and facilitates the hepatic uptake of most statins. Other relevant hepatic uptake transporters include OATP1B3, OATP2B1, OATP1A2 and the sodium-dependent taurocholate co-transporting polypeptide, NTCP [56, 57]. Genetic variability in membrane transport clearly influences statin-related clinical outcome (Table 2). Polymorphisms in candidate solute transporter genes are associated with the altered hepatic uptake of simvastatin  and pravastatin . Much of this variability can be attributed to two coding variants in SLCO1B1 (Asn130Asp and Val174Ala) . As outlined in a later section, the latter variant has since been shown to be highly informative in determining risk for the development of toxicity to simvastatin.
As shown in Figure 2, other transporters can influence the development and severity of statin-induced muscle toxicity as well. Many statins are substrates for efflux transporters such as multidrug resistance protein MDR1 (gene name ABCB1) or multidrug resistance-associated protein MRP2 (gene name ABCC2) . Located on the canalicular membrane of hepatocytes, these ATP-binding cassette proteins mediate the final step in the hepatobiliary clearance of statins. It therefore seems likely that variability in the activity of these transporters would alter the course of statin related clinical events. Genotype-phenotype association study performed in a cohort of 116 hypercholesterolemic patients has revealed that ABCB1 variants influence the efficacy of simvastatin. The same analysis also revealed that ABCB1 variants (1236T, 2677 non-G and 3435T) were less frequent in patients with adverse muscle effects .
Change in the activity of another efflux protein from the same family (gene name ABCG2) further alters the pharmacokinetics of most statins [63, 64]. This is particularly true for atorvastatin and rosuvastatin, two of the most potent drugs in the class [63, 64]. In a recent study of 305 Chinese patients treated with 10 mg rosuvastatin daily, one SNP in ABCG2 (rs2231142) was strongly associated with statin efficacy . Participants carrying a CC genotype at rs2231142 had a 6.9% greater reduction in LDL cholesterol levels compared to those with AA genotypes. Because the frequency of this genotype differs widely by race, it may explain a significant portion of the increased myopathy risk observed in Asians [36, 65]. Other ABCG2 variants may be involved as well [66, 67].
The rate-limiting enzyme in cholesterol biosynthesis is HMGCR. Statins inhibit the activity of this enzyme. It has been demonstrated in multiple studies that genetic variants in the HMGCR gene are important determinants for statin efficacy [68–74]. Although it seems reasonable to assume that those polymorphisms would also alter risk of statin toxicity, Frudakis et al. failed to observe any association between HMGCR variants and statin myopathy in a well-designed case-control study (263 samples) . Thus, alteration in cholesterol biosynthesis alone might not be sufficient to induce myopathy . Inhibition of HMGCR also attenuates the levels of many distal intermediates (Figure 1A) . After the generation of mevalonic acid, the pathway subsequently produces geranyl pyrophosphate (10 carbons), farnesyl pyrophosphate (15 carbons) and geranylgeranyl pyrophosphate (20 carbons). The isoprenoid side chains of these biosynthetic intermediates can transfer farnesyl or geranyl moieties to C-terminal cysteine(s) of target proteins, through a process call “protein prenylation.”
Because prenylation is necessary for synthesizing the side chain within uibiquinone (coenzyme Q10, CoQ10), statins may disturb the integrity of electron transport within the mitochondria (Figure 1A). As such, mitochondrial dysfunction due to altered levels of CoQ10 has been suggested as a potential mechanism for statin myopathy. The role of mitochondrial dysfunction in the pathogenesis of statin-induced myopathy is supported by extensive pathological evidence [76–78]. Vladutiu et al. have demonstrated that 52% of muscle biopsies from patients with statin related myalgias revealed mitochondrial abnormalities, and 31% of these biopsies revealed multiple defects . Further work by the same group identified variants in adenosine monophosphate deaminase (AMPD1), myophosphorylase (PYGM) and carnitine palmitoyltransferase II (CPT2) as contributors to risk . Additional pharmacodynamic variants contributing to myopathy risk have been reviewed by Peters et al. (including subclinical McArdle disease) (Table 2) . Oh et al. genotyped two SNPs in COQ2 (encoding an important enzyme in CoQ10 biosynthesis) in 133 statin-induced myopathy cases and 158 matched controls . Both SNPs were associated with increased risk of statin intolerance, and a haplotype based on these variants yielded an even stronger association (2.5-fold increase in risk) . These observations have led a number of investigators to explore the possibility that statin myopathy could be attenuated by coadministration of oral CoQ10 . The trials, however, have been small, and the results have been disappointing; e.g., Young et al. randomized 44 patients, who had previously failed statins due to muscle pain, to receive simvastatin with either placebo or CoQ10 supplementation (200 mg/day). No difference in myalgia score was observed between the treatment groups. Thus, oral CoQ10 did not improve statin tolerance .
Prenylation also influences the balance between myocyte viability and apoptosis Statin-induced apoptosis has been demonstrated in vitro, using myotubes , myoblasts , and differentiated primary human skeletal muscle cells . This effect can be reproduced by geranyl-geranyl-transferase inhibitors, and rescued by replacement of mevalonic acid . Compelling evidence suggests that statins cause apoptosis in skeletal muscle by disrupting the prenylation of small G proteins like Rho , Rab , and Rap . For example, statins induce apoptosis at concentrations that suppress the prenylation of Rap1a (a 21kD GTPase) , and Itagaki et al. have shown that this process is accompanied by the redistribution of small G proteins in myoblasts . It remains unclear, however, whether the altered prenylation of small G-proteins is necessary and sufficient to produce myopathy in vivo, or whether myocyte apoptosis is first activated by disrupted Ca2+ homeostasis following mitochondrial injury.
Although the candidate gene approach was widely applied in the identification of genes responsible for complex diseases, and evolutionarily important quantitative traits, the utility of this approach is largely limited by its reliance upon a priori knowledge . On the other hand, genome-wide approaches usually proceed without any presuppositions regarding the importance of specific functional features of the traits being investigated. Genome-wide approaches include linkage studies (in families) and genome wide association studies (in unrelated individuals). Both approaches represent unbiased hypothesis-free experiments that hold the potential to identify new biology .
Linkage studies represent the earliest type of whole-genome scanning. By constructing pedigrees, early linkage analyses tested for the joint transmission of chromosomal segments and complex phenotypic traits within families. Linkage is the method of choice to identify rare variants with a large impact on disease risk if the trait aggregates in families. The diseases caused by such variants show obvious inheritance patterns and are typically called Mendelian diseases . Although powerful, linkage analyses typically only localize ~10 to 100 cM (centiMorgans) intervals because of the limited number of recombination events within pedigrees [91, 92]. Furthermore, this approach has limited capacity for identifying genes with low penetrance and modest effect size. Thus, the main advantage of family-based studies is that they are not susceptible to false positives from racial admixture and population stratification. Linkage studies have typically not been applied routinely within a pharmacogenomic context, due to the difficulty in identifying families with multiple members exposed to the same drug at the same dose.
Conversely, genome-wide association studies (GWAS) are applied often within the context of pharmacogenomics, for large cohorts of unrelated individuals . GWAS conducted in randomized controlled trials (RCTs) can provide an unbiased survey of the genomic architecture underlying treatment outcome. It is now possible to examine large numbers of polymorphisms, on the order of 100000–1000000, across the entire genome using highly parallel genotyping arrays . In 2008, the SEARCH Collaborative Group applied a 317K SNP scan to 85 cases of incipient myopathy and 90 frequency-matched drug exposed controls, to identify markers of muscle toxicity specifically within the context of high dose simvastatin (80 mg daily) . This was the first published genome-wide association study of statin-induced muscle toxicity. A single variant survived statistical correction for multiple testing: a base substitution in the SLCO1B1 gene .
After genomic re-sequencing of SLCO1B1, the putative causative allele (Val174Ala) was re-tested for association in a subset of definite myopathy cases from the original SEARCH study cohort, revealing an odds ratio for myopathy of 4.5 per copy of the variant allele . This association has since been replicated in several independent study populations [95–98]. In the Heart Protection Study (HPS), 24 cases of myopathy were identified in 10,269 participants receiving primary prevention with a lower dose of simvastatin (40 mg daily); 21 were genotyped retrospectively for the variant identified in SEARCH , and the relative risk was 2.6 per copy of the variant allele. In a practice-based setting, where the definition of intolerance includes discontinuation of the drug for any reason, the relative risk appears to be closer to 1.5 [96–98]. Efforts are now being made to move this pharmacogenetic association into clinical practice through the application of novel decision-support mechanisms [99, 100]
GWAS using statin-induced myopathy cases may also provide deeper insight into the underlying mechanism of toxicity (i.e., leveraging the genetics to inform the biology) [101, 102]. In 2011, Marciante et al. published a combined CGAS - GWAS using a cohort of 185 confirmed cerivastatin-induced myopathy cases (CK >10xULN with pain) and 732 matched controls . In addition to replicating the well-established SLCO1B1 association for another statin (odds ratio 1.9, p = 0.002), Marciante et al. also leveraged GWAS to identify an association between cerivastatin-induced myopathy and an intronic SNP (rs2819742) in the ryanodine receptor 2 gene (RYR2) (odds ratio 0.48, p = 1.74−07) . Other GWAS cohorts are providing new candidates for statin-related myopathy Muscle specific genes (e.g., gene products modulating Ca2+ flux and excitation-contraction coupling) represent attractive targets for mechanism-based study in vitro.
GWAS have identified numerous common variants associated with complex diseases and provided valuable insights into their fundamental genetic architecture. The main strength of GWAS is an ability to discover truly novel candidate SNPs/genes. However, GWAS are vulnerable to selection bias for “top hits.” . Genes highly relevant in the context of interaction with other important variants will not be detectable at conventional levels of significance . Thus, variants contributing to myopathy risk must be considered within the context of gene-gene (GxG) and gene-environment (GxE) interactions as well. Grouping SNPs into functional units (either genes or groups of genes) might improve the signal for association with a trait. The gene-based approach estimates the combined effect of markers within a gene rather than each marker individually, and is potentially more powerful than individual SNP-based analysis .A pathway based approach (interactive sets based on gene ontology, and biological networks based on structural or functional similarity) may facilitate a deeper understanding of findings related to determinants of statin outcome, and account for genetic heterogeneity, as well as improve statistical power to detect relevant genes .
Accurate risk prediction models must be flexible enough to simultaneously integrate clinical and genetic factors influencing the absorption, distribution, metabolism, and elimination of statins (i.e., within processes impacting PK), and they must be robust enough to further interpret such determinants in the context of variables influencing mechanism (i.e., across PD and PK processes). One potential approach includes the application of the worlds rapidly growing annotated databases (e.g., GO, gene ontology and gene network analysis). The GO provides structured information on the properties of the products of genes, using three general domains: cellular components (the parts of a cell or its extracellular environment), molecular function (the activities of gene products at the molecular level), and biological process (sets of molecular events with a defined beginning and end) (http://www.geneontology.org/). Clinical responses to statins are complex traits that likely extend from many gene products highly interconnected within networks. For purposes of illustration, we uploaded PD and PK candidate genes into an IPA network analysis (http://www.ingenuity.com/), to uncover the biological and functional interactions among these candidates relevant to statin-related clinical outcome (Figures 1B, ,2B).2B). This software application resolves biological and chemical interactions between genes and gene products with functional commonalities. A similar approach was used by Baye et al (2011) to validate recursive partitioning (RP) based gene-gene interactions in childhood asthma .
The scientific community is also moving toward whole genome sequencing, in cohorts of increasing sample size, and it is therefore inevitable that risk prediction models for statin-induced myopathy will need to be sufficiently robust to accommodate huge numbers of functional variants that are relatively rare . Most SNP genotyping platforms used to date have only captured polymorphisms at a frequency of 5% or more. Thus, rare variants have been missed in most myopathy GWAS. Rare variants, however, are becoming increasingly catalogued within large populations through efforts like The 1000 Genomes Project (www.1000genomes.org) and next generation sequencing (NGS) [not only to analyze genomic DNA but also as a valuable tool to dissect whole transcriptomes (mRNA, ncRNA, microRNAs) and epigenomes (DNA methylation, protein-DNA interaction)]. While these populations could serve as controls (i.e., statin-exposed, non-myopathic individuals) for sequencing efforts directed at myopathic individuals, sophisticated analytical methodologies will be required [109, 110]. Collapsing methods based on functionality look promising . Using knowledge in annotated databases, functional elements can be weighted by a variety of strategies, such as the regulation of expression through prediction methods that assess transcription factor binding (TF search, ConSite and TRANSFAC), enhancer interactions (VISTA enhancer browser), and mRNA stability (miRBase, TargetScan).
Lastly, these genomic approaches will need to be supported by additional technologies related to rigorous phenotyping (e.g., transcriptomics, proteomics and lipidomics). Transcriptional profiling arrays now consider alternative splicing , tissue specific gene expression , and evolutionary aspects of gene expression . Laaksonen et al. analyzed the expression of over 46,000 genes in muscle biopsy samples obtained from 6 subjects receiving atorvastatin, 5 receiving simvastatin (1 of 6 cases yielded insufficient RNA), and 6 receiving placebo (limited to males, and frequency-matched according to age) . Simvastatin treatment resulted in expression change in 111 genes - 26 down-regulated and 85 up-regulated. More than twenty biological pathways were affected according to their bioinformatics analysis. The most significant up-regulated genes included ALOX5AP, CCL5, COL3A1 MYL5 and MYBPH. The same muscle biopsy specimens were then also characterized by lipidomics (LC tandem MS), quantifying 132 unique molecular lipid species [114, 115]. Regression of lipidomic data on gene expression data for pathway-based signaling networks confirmed the involvement of lipid-derived signaling pathways (e.g., prostanoid biosynthesis), and suggested a role for Ca2+-dependent pathways capable of modulating excitation contraction coupling and apoptosis (e.g., Phospholipase C, PLC) [114, 115]. Thus, comprehensive approaches linking genomics and transcriptional profiling, with novel phenotyping strategies in the context of larger populations hold the potential to define the genetic architecture of statin-induced myopathy with unprecedented power.
Like most complex traits, statin-induced muscle toxicity is highly variable in its clinical presentation, and strongly influenced by both GxG and GxE interactions. Although candidate genes studies have revealed a number of variants influencing statin PK, only one locus (SLCO1B1) has been validated using genome-wide technologies. Additional PD variants are emerging from GWAS (RYR2), reinforcing the claim that the field of genomics will inform the underlying biology. While these results move future mechanistic queries in the direction of Ca2+ homeostasis, additional population-based studies are needed to place the findings into context alongside huge numbers of rare variants in cohorts that are diverse enough to facilitate analyses adjusting for admixture.
As the translational genomics community continues moving toward prospective implementation of ancestry- and gene-based dosing models for a variety of commonly used medications, some groups  have argued that statins should be among the first drugs to benefit from automated decision support utilizing these data. The Clinical Pharmacogenomics Implementation Consortium (CPIC) routinely publishes guidelines for such an approach . To date, three such guidelines have been published (for warfarin, clopidogrel, and thiopurines), and it is anticipated that a guideline for gene-based statin dosing will be among the next publications to be posted: http://www.pharmgkb.org/contributors/consortia/cpic_gene-drug_pairs.jsp. As medical centers across diverse systems of care begin adopting this approach, the health care enterprise will soon be in a position to assess impact for pharmacogenomics on quality of care for the most commonly prescribed class of drugs in the industrialized world.