SLCO1B1 rs4149056 polymorphism associated with statin-induced myopathy is differently distributed according to ethnicity in the Brazilian general population: Amerindians as a high risk ethnic group
© Santos et al; licensee BioMed Central Ltd. 2011
Received: 1 March 2011
Accepted: 12 October 2011
Published: 12 October 2011
Recent studies reported the association between SLCO1B1 polymorphisms and the development of statin-induced myopathy. In the scenario of the Brazilian population, being one of the most heterogeneous in the world, the main aim here was to evaluate SLCO1B1 polymorphisms according to ethnic groups as an initial step for future pharmacogenetic studies.
One hundred and eighty-two Amerindians plus 1,032 subjects from the general urban population were included. Genotypes for the SLCO1B1 rs4149056 (c.T521C, p.V174A, exon 5) and SLCO1B1 rs4363657 (g.T89595C, intron 11) polymorphisms were detected by polymerase chain reaction followed by high resolution melting analysis with the Rotor Gene 6000® instrument.
The frequencies of the SLCO1B1 rs4149056 and rs4363657 C variant allele were higher in Amerindians (28.3% and 26.1%) and were lower in African descent subjects (5.7% and 10.8%) compared with Mulatto (14.9% and 18.2%) and Caucasian descent (14.8% and 15.4%) ethnic groups (p < 0.001 and p < 0.001, respectively). Linkage disequilibrium analysis show that these variant alleles are in different linkage disequilibrium patterns depending on the ethnic origin.
Our findings indicate interethnic differences for the SLCO1B1 rs4149056 C risk allele frequency among Brazilians. These data will be useful in the development of effective programs for stratifying individuals regarding adherence, efficacy and choice of statin-type.
KeywordsSLCO1B1 rs4149056 statins myopathy Amerindian pharmacogenetic
Statins, inhibitors of 3-hydroxy-3-methylglutaryl-coenzyme A (HMG-CoA), are widely prescribed drugs to reduce cardiovascular morbidity and mortality and their use is well tolerated by the majority of patients [1, 2]. Nonetheless, statins can lead to myopathy (i.e., muscle pain) with symptoms ranging from mild myalgia without elevated serum creatine kinase to fatal rhabdomyolysis (i.e., muscle breakdown and myoglobin release) [2, 3]. In clinical practice, the incidence of mild statin-induced side effects (about 5-10%) appears to be greater than in controlled trials and the mechanisms by which statins lead to myopathy remain unclear. However, studies have suggested that blood statin concentrations, daily dose, female gender, hepatic or renal dysfunction, and concomitant drugs are all predictors of statin-induced myopathy [4–7].
In addition, some recent studies reported the association between SLCO1B1 polymorphisms and the development of musculoskeletal side effects [8–10]. The SLCO1B1 gene (Gene ID: 10599) encodes the organic anion-transporting polypeptide 1B1 (OATP1B1) that plays a crucial role in the hepatic uptake and clearance of albumin-bound organic compounds . In particular, the SLCO1B1 rs4149056 functional polymorphism (c.T521C, p.V174A) in exon 5, leads to higher statin circulating concentration [12, 13]. SLCO1B1 rs4149056 was associated with an odds ratio for myopathy of 4.5 (95% CI, 2.6 to 7.7) per copy of the C allele and 16.9 (95% CI, 4.7 to 61.1) among CC homozygotes as compared with the TT genotype, in a genome-wide association study plus replication in a trial of simvastatin 40 mg daily involving 20,000 participants . The frequencies of SLCO1B1 polymorphisms vary significantly among different world-wide populations, leaving the impact of pharmacogenetic testing for these variants extremely dependent on a particular population genetic architecture.
The Brazilian population is one of the most heterogeneous in the world, being a mixture of different ethnic groups, composed mainly of European (Caucasian descent), African descent and Amerindians [14, 15]. Thus, the main aim of the present study was to evaluate the SLCO1B1 polymorphisms according to ethnic groups as an initial step for future pharmacogenetic studies and screening programmes.
This study included 1,032 subjects of the general urban population selected from the Hearts of Brazil Project (HBP) . The universe of the HBP consisted in the set of inhabitants of Brazilian urban centers with more than 100,000 inhabitants. The HBP sample plan was calculated as 2,500 interviews, distributed in 72 cities from the 5 regions of the country proportionally to the number of inhabitants, per sex and age range, based on data from IBGE (Brazilian Institute of Geography and Statistic). In the selected cities, the "households" constituted the second-stage units, with one interview per household. Subjects were separated in self-identified sub-groups according to ethnicity, as Caucasian descent, African descent, or Mulattos (considered mixed ethnic subjects) . In addition, one hundred and eighty-two Amerindians derived from two different groups (Guarani and Tupinikin; Aracruz Indian Reserve, Espirito Santo State in the Southeast Brazilian coast) were also analyzed. The study protocol was approved by the involved Institutional Ethics Committees and National Ethic Committee for Human Research (CONEP Register Number 4599), and written informed consent was obtained from all participants prior to entering the study.
Genomic DNA from subjects was extracted from peripheral blood following standard salting-out procedure . Genotypes for the SLCO1B1 rs4149056 (c.T521C, p.V174A, exon 5) and SLCO1B1 rs4363657 (g.T89595C, intron 11) polymorphisms were detected by polymerase chain reaction (PCR) followed by high resolution melting (HRM) analysis with the Rotor Gene 6000® instrument (Qiagen, Courtaboeuf, France) [18–20]. The QIAgility® (Qiagen, Courtaboeuf, France), an automated instrument, was used according to instructions to optimize the sample preparation step (approximately 30 minutes). One specific disc is able to genotype 96 samples for these polymorphisms.
Amplification of the fragments was performed using the primer sense 5'-TTGTTTAAAGGAATCTGGGTCA-3' and antisense 5'-GAGTCTCCCCTATTCCACGA-3' (77 pairs base) for the rs4149056; and sense 5'-GAGTCCTTCTTTCTCAATTTTTCA-3' and antisense 5'-AAAAGCCCTAGACCAAATGC-3' (107 pairs base) for the rs4363657. A 35-cycle PCR was carried out with the following conditions: denaturation of the template DNA for first cycle of 94°C for 120 s, denaturation of 94°C for 20 s, annealing of 53.4°C (rs4149056) or of 50.5°C (rs4363657) for 20 s, and extension of 72°C for 22 s. PCR was performed using a 10 μL reactive solution (10 mM Tris-HCl, 50 mM KCl, pH 9.0; 2.0 mM MgCl2; 200 μM of each dNTP; 0.5 U Taq DNA Polymerase; 200 nM of each primer; 10 ng of genomic DNA template) with addition of fluorescent DNA-intercalating SYTO9® (1.5 μM; Invitrogen, Carlsbad, USA).
Chi-square test was performed for comparative analysis of the SLCO1B1 polymorphism frequencies according to ethnicity (African descent, Mulatto, Caucasian descent, and Amerindian). Chi-square test was also performed for the comparative analysis between our data and previous reports. Hardy-Weinberg equilibrium and linkage disequilibrium analyses were conducted with Haploview 4.0. All statistical analyses were carried out using the SPSS software (v. 16.0), with the level of significance set at p < 0.05.
General data of the studied population sample
Of the 1,214 subjects (mean age 43.1 ± 14.2), 647 (53.3%) were female and 567 (46.7%) male. Ethnic group distribution was: African descent 8.0% (n = 97), Mulatto 27.3% (n = 332), Caucasian descent 49.7% (n = 603), and Amerindians 15.0% (n = 182).
Frequencies of the SLCO1B1polymorphisms
Allelic and genotypic frequencies for the SLCO1B1 polymorphisms according to ethnic groups
Nucleotide Change and Genotype
The frequencies of the SLCO1B1 rs4363657 C variant allele and of the homozygous genotype (CC) were also higher in Amerindians (26.1% and 6.6%) and were lower in African descent subjects (10.8% and 0%) compared with Mulatto (18.2% and 1.5%) and Caucasian descent (15.4% and 2.7%) ethnic groups (p < 0.001 and p < 0.001, respectively) (Table 1).
The genotypic distributions for the SLCO1B1 rs4149056 and rs4363657 polymorphisms were in Hardy-Weinberg equilibrium according to ethnic groups.
Linkage disequilibrium analysis between SLCO1B1rs4149056 and rs4363657 variant alleles according to ethnic groups
The main finding of the present study was to observe significant differences between the allelic and genotypic frequencies of SLCO1B1 rs4149056 polymorphism associated with statin-induced myopathy, according to ethnicity in the Brazilian general population. To our knowledge this is the first study indicating that Brazilian Amerindians present higher frequencies of rs4149056 C variant allele (28.3%) and probably higher risk for statin-induced myopathy. These findings could be useful in the strategic planning of the implementation of pharmacogenetic testing according to ethnicity.
The SLCO1B1 rs4149056 C variant allele frequency found in our study was higher (p < 0.001) in Amerindian subjects (28.3%) compared with other ethnic groups and with Turkish  (12.2%); English  (13.0%); European American  (14.3%); and German  (15.0%) subjects. However, the Caucasian descent (14.8%) and Mulatto (14.9%) ethnic groups of our study presented similar (p > 0.050) rs4149056 C variant allele frequency compared with Caucasian descent subjects from Europe and USA. In addition, our African descent (5.7%) subject group had similar (p > 0.050) frequency compared with African American  (2.3%), Uganda African  (3.9%), and Tanzania  (6.0%) subjects. Interestingly, the frequency observed for the Mulatto subjects indicates that the Amerindian ancestry may be present in this multi-ethnic group, besides the presence of known African ancestry.
For the SLCO1B1 rs4363657 polymorphism, C allele frequency found in Caucasian descent in our study (15.4%) was similar (p > 0.050) to that reported in English subjects (13.0%) who participated as controls in a genome-wide association study (GWAS) . The non-coding rs4363657 polymorphism was in nearly complete linkage disequilibrium (LD = 97) with the non-synonymous rs4149056 polymorphism, in the population sample in which the cited GWAS was performed . In our population sample plus Amerindians, it is possible to observe that the SLCO1B1 rs4149056 and rs4363657 are in different linkage disequilibrium patterns depending on the ethnic origin. The African descent (LD = 89) and Mulatto (LD = 81) ethnic groups present higher degree of linkage disequilibrium. This highlights two main points: first, even with the different linkage disequilibrium among ethnic groups, the haplotype analysis are not need because only the rs4149056 polymorphism have been described as functional [11, 24]. Second, the different frequencies among ethnic groups for the polymorphisms used as markers in the GWAS platforms are limitations/peculiarities of the GWAS methodology , which may be acting as interferents in replication studies, mainly in countries with admixed population.
The association between SLCO1B1 rs4149056 with myopathy (odds ratio: 4.5 - 95% CI 2.6 to 7.7, per copy of the C allele; and odds ratio: 16.9 - 95% CI 4.7 to 61.1, among CC homozygotes as compared with TT genotypes) was observed in 85 subjects with myopathy taking 80 mg of simvastatin daily as part of a trial involving 12,000 participants; plus replication in a trial of 40 mg of simvastatin daily involving 20,000 participants . The authors of this GWAS plus replication study concluded that common variants in the SLCO1B1 gene are strongly associated with an increased risk of statin-induced myopathy and genotyping may help to obtain benefits of statin therapy more safely and effectively .
Other studies were performed to replicate this association. Brunham et al observed that rs4149056 genotype was significantly associated (odds ratio 2.3 per C allele) with myopathy in patients receiving simvastatin (mean ~30 mg daily), but not in patients who received atorvastatin . In this way, genotype CC individuals compared with wild-type TT individuals, the area under the plasma concentration time curve for simvastatin acid was increased 221% and for atorvastatin was increased 145% [12, 13]. Voora et al confirmed, in 99 subjects with statin-related adverse events, the association previously identified, but also expanded this finding to the most commom statin-induced side effects (e.g., myalgia or muscle ache without significant creatine kinase elevations) . Puccetti et al reported the same association, but in a group of 46 subjects affected by familial hypercholesterolemia treated with atorvastatin and developing myopathy .
There are potential limitations in our study. First, the self-referred ethnicity could be creating wrong separation. However, this form of classification is the one most probable to be encountered in real-life situations in the clinical practice. In addition, even self-referred ethnicity was able to clearly differentiate groups of individuals with rather different allele and genotype frequencies. Second, in order to implement a genotyping program in this context, some points could still be evaluated, for example: other genetic markers, primary disease, use of concomitant medications, types of statins, age and gender .
Our findings indicate interethnic differences for the SLCO1B1 rs4149056 C risk allele frequency among Brazilian general population plus Amerindians. Recent studies indicated the genotyping for this functional polymorphism as an important clinical tool in future patient-tailored programs. Data of the present study will be useful in the development of effective programs for stratifying individuals regarding adherence, efficacy and choice of statin-type for patients who need statin, specially in high-dose (e.g., familial hypercholesterolemia).
PCJLS is recipient from fellowship from FAPESP, Proc. 2010-17465-8, Brazil. The technical assistance of the Laboratory of Genetics and Molecular Cardiology group, Heart Institute group is gratefully acknowledged.
- Baigent C, Keech A, Kearney PM, Blackwell L, Buck G, Pollicino C, Kirby A, Sourjina T, Peto R, Collins R, et al: Efficacy and safety of cholesterol-lowering treatment: prospective meta-analysis of data from 90,056 participants in 14 randomised trials of statins. Lancet. 2005, 366 (9493): 1267-1278.View ArticlePubMed
- Ghatak A, Faheem O, Thompson PD: The genetics of statin-induced myopathy. Atherosclerosis. 210 (2): 337-343.
- Thompson PD, Clarkson PM, Rosenson RS: An assessment of statin safety by muscle experts. Am J Cardiol. 2006, 97 (8A): 69C-76C.View ArticlePubMed
- Armitage J: The safety of statins in clinical practice. Lancet. 2007, 370 (9601): 1781-1790. 10.1016/S0140-6736(07)60716-8.View ArticlePubMed
- Gotto AM: Statins, cardiovascular disease, and drug safety. Am J Cardiol. 2006, 97 (8A): 3C-5C.View ArticlePubMed
- Law M, Rudnicka AR: Statin safety: a systematic review. Am J Cardiol. 2006, 97 (8A): 52C-60C.View ArticlePubMed
- Nichols GA, Koro CE: Does statin therapy initiation increase the risk for myopathy? An observational study of 32,225 diabetic and nondiabetic patients. Clin Ther. 2007, 29 (8): 1761-1770. 10.1016/j.clinthera.2007.08.022.View ArticlePubMed
- Brunham LR, Lansberg PJ, Zhang L, Miao F, Carter C, Hovingh GK, Visscher H, Jukema JW, Stalenhoef AF, Ross CJ, et al: Differential effect of the rs4149056 variant in SLCO1B1 on myopathy associated with simvastatin and atorvastatin. Pharmacogenomics J. 2011.
- Link E, Parish S, Armitage J, Bowman L, Heath S, Matsuda F, Gut I, Lathrop M, Collins R: SLCO1B1 variants and statin-induced myopathy--a genomewide study. N Engl J Med. 2008, 359 (8): 789-799.View ArticlePubMed
- Voora D, Shah SH, Spasojevic I, Ali S, Reed CR, Salisbury BA, Ginsburg GS: The SLCO1B1*5 genetic variant is associated with statin-induced side effects. J Am Coll Cardiol. 2009, 54 (17): 1609-1616. 10.1016/j.jacc.2009.04.053.PubMed CentralView ArticlePubMed
- Niemi M, Pasanen MK, Neuvonen PJ: Organic Anion Transporting Polypeptide 1B1: a Genetically Polymorphic Transporter of Major Importance for Hepatic Drug Uptake. Pharmacol Rev. 63 (1): 157-181.
- Pasanen MK, Fredrikson H, Neuvonen PJ, Niemi M: Different effects of SLCO1B1 polymorphism on the pharmacokinetics of atorvastatin and rosuvastatin. Clin Pharmacol Ther. 2007, 82 (6): 726-733. 10.1038/sj.clpt.6100220.View ArticlePubMed
- Pasanen MK, Neuvonen M, Neuvonen PJ, Niemi M: SLCO1B1 polymorphism markedly affects the pharmacokinetics of simvastatin acid. Pharmacogenet Genomics. 2006, 16 (12): 873-879. 10.1097/01.fpc.0000230416.82349.90.View ArticlePubMed
- de Lima Santos PC, de Oliveira Alvim R, Ferreira NE, de Sa Cunha R, Krieger JE, Mill JG, Pereira AC: Ethnicity and arterial stiffness in Brazil. Am J Hypertens. 24 (3): 278-284.
- Santos PC, Pereira AC, Cancado RD, Schettert IT, Sobreira TJ, Oliveira PS, Hirata RD, Hirata MH, Figueiredo MS, Chiattone CS, et al: HFE gene mutations in patients with primary iron overload: is there a significant improvement in molecular diagnosis yield with HFE sequencing?. Blood Cells Mol Dis. 45 (4): 302-307.
- Makdisse M, Pereira Ada C, Brasil Dde P, Borges JL, Machado-Coelho GL, Krieger JE, Nascimento Neto RM, Chagas AC: Prevalence and risk factors associated with peripheral arterial disease in the Hearts of Brazil Project. Arq Bras Cardiol. 2008, 91 (6): 370-382.View ArticlePubMed
- Miller SA, Dykes DD, Polesky HF: A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988, 16 (3): 1215-10.1093/nar/16.3.1215.PubMed CentralView ArticlePubMed
- Alvim RO, Freitas SR, Ferreira NE, Santos PC, Cunha RS, Mill JG, Krieger JE, Pereira AC: APOE polymorphism is associated with lipid profile, but not with arterial stiffness in the general population. Lipids Health Dis. 9: 128-
- Santos PC, Soares RA, Krieger JE, Guerra-Shinohara EM, Pereira AC: Genotyping of the hemochromatosis HFE p.H63D and p.C282Y mutations by high-resolution melting with the Rotor-Gene 6000((R)) instrument. Clin Chem Lab Med. 2011.
- Santos PC, Soares RA, Santos DB, Nascimento RM, Coelho GL, Nicolau JC, Mill JG, Krieger JE, Pereira AC: CYP2C19 and ABCB1 gene polymorphisms are differently distributed according to ethnicity in the Brazilian general population. BMC Med Genet. 12: 13-
- Mwinyi J, Kopke K, Schaefer M, Roots I, Gerloff T: Comparison of SLCO1B1 sequence variability among German, Turkish, and African populations. Eur J Clin Pharmacol. 2008, 64 (3): 257-266. 10.1007/s00228-007-0409-y.View ArticlePubMed
- Tirona RG, Leake BF, Merino G, Kim RB: Polymorphisms in OATP-C: identification of multiple allelic variants associated with altered transport activity among European- and African-Americans. J Biol Chem. 2001, 276 (38): 35669-35675. 10.1074/jbc.M103792200.View ArticlePubMed
- Aklillu E, Mugusi S, Ngaimisi E, Hoffmann MM, Konig S, Ziesenitz V, Mikus G, Haefeli WE, Weiss J: Frequency of the SLCO1B1 388A > G and the 521T > C polymorphism in Tanzania genotyped by a new LightCycler(R)-based method. Eur J Clin Pharmacol. 2011.
- Ho RH, Choi L, Lee W, Mayo G, Schwarz UI, Tirona RG, Bailey DG, Michael Stein C, Kim RB: Effect of drug transporter genotypes on pravastatin disposition in European- and African-American participants. Pharmacogenet Genomics. 2007, 17 (8): 647-656. 10.1097/FPC.0b013e3280ef698f.PubMed CentralView ArticlePubMed
- Haiman CA, Stram DO: Exploring genetic susceptibility to cancer in diverse populations. Curr Opin Genet Dev. 20 (3): 330-335.
- Puccetti L, Ciani F, Auteri A: Genetic involvement in statins induced myopathy. Preliminary data from an observational case-control study. Atherosclerosis. 211 (1): 28-29.
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/12/136/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.