A gene risk score using missense variants in SLCO1B1 is associated with earlier onset statin intolerance

Background and aims The efficacy of statin therapy is hindered by intolerance to the therapy, leading to discontinuation. Variants in SLCO1B1, which encodes the hepatic transporter OATB1B1, influence statin pharmacokinetics, resulting in altered plasma concentrations of the drug and its metabolites. Current pharmacogenetic guidelines require sequencing of the SLCO1B1 gene, which is more expensive and less accessible than genotyping. In this study, we aimed to develop an easy, clinically implementable functional gene risk score (GRS) of common variants in SLCO1B1 to identify patients at risk of statin intolerance. Methods and results A GRS was developed from four common variants in SLCO1B1. In statin users from Tayside, Scotland, UK, those with a high-risk GRS had increased odds across three phenotypes of statin intolerance [general statin intolerance (GSI): ORGSI 2.42; 95% confidence interval (CI): 1.29–4.31, P = 0.003; statin-related myopathy: ORSRM 2.51; 95% CI: 1.28–4.53, P = 0.004; statin-related suspected rhabdomyolysis: ORSRSR 2.85; 95% CI: 1.03–6.65, P = 0.02]. In contrast, using the Val174Ala genotype alone or the recommended OATP1B1 functional phenotypes produced weaker and less reliable results. A meta-analysis with results from adjudicated cases of statin-induced myopathy in the PREDICTION-ADR Consortium confirmed these findings (ORVal174Ala 1.99; 95% CI: 1.01–3.95, P = 0.048; ORGRS 1.76; 95% CI: 1.16–2.69, P = 0.008). For those requiring high-dose statin therapy, the high-risk GRS was more consistently associated with the time to onset of statin intolerance amongst the three phenotypes compared with Val174Ala (GSI: HRVal174Ala 2.49; 95% CI: 1.09–5.68, P = 0.03; HRGRS 2.44; 95% CI: 1.46–4.08, P < 0.001). Finally, sequence kernel association testing confirmed that rare variants in SLCO1B1 are associated with the risk of intolerance (P = 0.02). Conclusion We provide evidence that a GRS based on four common SLCO1B1 variants provides an easily implemented genetic tool that is more reliable than the current recommended practice in estimating the risk and predicting early-onset statin intolerance

Evidence of a Causal Relationship between Serum Thyroid-Stimulating Hormone and Osteoporotic Bone Fractures

OBJECTIVE: We aimed to validate the association of genome-wide association study (GWAS)-identified loci and polygenic risk score with serum thyroid-stimulating hormone (TSH) concentrations and the diagnosis of hypothyroidism. Then, the causal relationship between serum TSH and osteoporotic bone fracture risk was tested. METHODS: A cross-sectional study was done among patients of European Caucasian ethnicity recruited in Tayside (Scotland, UK). Electronic medical records (EMRs) were used to identify patients and average serum TSH concentration and linked to genetic biobank data. Genetic associations were performed by linear and logistic regression models. One-sample Mendelian randomization (MR) was used to test causality of serum TSH on bone fracture risk. RESULTS: Replication in 9,452 euthyroid individuals confirmed known loci previously reported. The 58 polymorphisms accounted for 11.08% of the TSH variation (p < 1e−04). TSH-GRS was directly associated with the risk of hypothyroidism with an odds ratio (OR) of 1.98 for the highest quartile compared to the first quartile (p = 2.2e−12). MR analysis of 5,599 individuals showed that compared with those in the lowest tertile of the TSH-GRS, men in the highest tertile had a decreased risk of osteoporotic bone fracture (OR = 0.59, p = 2.4e−03), while no difference in a similar comparison was observed in women (OR = 0.93, p = 0.61). Sensitivity analysis yielded similar results. CONCLUSIONS: EMRs linked to genomic data in large populations allow replication of GWAS discoveries without additional genotyping costs. This study suggests that genetically raised serum TSH concentrations are causally associated with decreased bone fracture risk in men

Genome-wide association study of angioedema induced by angiotensin-converting enzyme inhibitor and angiotensin receptor blocker treatment

Angioedema in the mouth or upper airways is a feared adverse reaction to angiotensin-converting enzyme inhibitor (ACEi) and angiotensin receptor blocker (ARB) treatment, which is used for hypertension, heart failure and diabetes complications. This candidate gene and genome-wide association study aimed to identify genetic variants predisposing to angioedema induced by these drugs. The discovery cohort consisted of 173 cases and 4890 controls recruited in Sweden. In the candidate gene analysis, ETV6, BDKRB2, MME, and PRKCQ were nominally associated with angioedema (p < 0.05), but did not pass Bonferroni correction for multiple testing (p < 2.89 × 10−5). In the genome-wide analysis, intronic variants in the calcium-activated potassium channel subunit alpha-1 (KCNMA1) gene on chromosome 10 were significantly associated with angioedema (p < 5 × 10−8). Whilst the top KCNMA1 hit was not significant in the replication cohort (413 cases and 599 ACEi-exposed controls from the US and Northern Europe), a meta-analysis of the replication and discovery cohorts (in total 586 cases and 1944 ACEi-exposed controls) revealed that each variant allele increased the odds of experiencing angioedema 1.62 times (95% confidence interval 1.05–2.50, p = 0.030). Associated KCNMA1 variants are not known to be functional, but are in linkage disequilibrium with variants in transcription factor binding sites active in relevant tissues. In summary, our data suggest that common variation in KCNMA1 is associated with risk of angioedema induced by ACEi or ARB treatment. Future whole exome or genome sequencing studies will show whether rare variants in KCNMA1 or other genes contribute to the risk of ACEi- and ARB-induced angioedema

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

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