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.

Hepatocyte mARC1 promotes fatty liver disease.

Background & Aims: Non-alcoholic fatty liver disease (NAFLD) has a prevalence of ∼25% worldwide, with significant public health consequences yet few effective treatments. Human genetics can help elucidate novel biology and identify targets for new therapeutics. Genetic variants in mitochondrial amidoxime-reducing component 1 (MTARC1) have been associated with NAFLD and liver-related mortality; however, its pathophysiological role and the cell type(s) mediating these effects remain unclear. We aimed to investigate how MTARC1 exerts its effects on NAFLD by integrating human genetics with in vitro and in vivo studies of mARC1 knockdown. Methods: Analyses including multi-trait colocalisation and Mendelian randomisation were used to assess the genetic associations of MTARC1. In addition, we established an in vitro long-term primary human hepatocyte model with metabolic readouts and used the Gubra Amylin NASH (GAN)-diet non-alcoholic steatohepatitis mouse model treated with hepatocyte-specific N-acetylgalactosamine (GalNAc)-siRNA to understand the in vivo impacts of MTARC1. Results: We showed that genetic variants within the MTARC1 locus are associated with liver enzymes, liver fat, plasma lipids, and body composition, and these associations are attributable to the same causal variant (p.A165T, rs2642438 G>A), suggesting a shared mechanism. We demonstrated that increased MTARC1 mRNA had an adverse effect on these traits using Mendelian randomisation, implying therapeutic inhibition of mARC1 could be beneficial. In vitro mARC1 knockdown decreased lipid accumulation and increased triglyceride secretion, and in vivo GalNAc-siRNA-mediated knockdown of mARC1 lowered hepatic but increased plasma triglycerides. We found alterations in pathways regulating lipid metabolism and decreased secretion of 3-hydroxybutyrate upon mARC1 knockdown in vitro and in vivo. Conclusions: Collectively, our findings from human genetics, and in vitro and in vivo hepatocyte-specific mARC1 knockdown support the potential efficacy of hepatocyte-specific targeting of mARC1 for treatment of NAFLD. Impact and implications: We report that genetically predicted increases in MTARC1 mRNA associate with poor liver health. Furthermore, knockdown of mARC1 reduces hepatic steatosis in primary human hepatocytes and a murine NASH model. Together, these findings further underscore the therapeutic potential of targeting hepatocyte MTARC1 for NAFLD.

MLX plays a key role in lipid and glucose metabolism in humans: Evidence from in vitro and in vivo studies.

BACKGROUND AND AIM: Enhanced hepatic de novo lipogenesis (DNL) has been proposed as an underlying mechanism for the development of NAFLD and insulin resistance. Max-like protein factor X (MLX) acts as a heterodimer binding partner for glucose sensing transcription factors and inhibition of MLX or downstream targets has been shown to alleviate intrahepatic triglyceride (IHTG) accumulation in mice. However, its effect on insulin sensitivity remains unclear. As human data is lacking, the aim of the present work was to investigate the role of MLX in regulating lipid and glucose metabolism in primary human hepatocytes (PHH) and in healthy participants with and without MLX polymorphisms. METHODS: PHH were transfected with non-targeting or MLX siRNA to assess the effect of MLX knockdown on lipid and glucose metabolism, insulin signalling and the hepatocellular transcriptome. A targeted association analysis on imputed genotype data for MLX on healthy individuals was undertaken to assess associations between specific MLX SNPs (rs665268, rs632758 and rs1474040), plasma biochemistry, IHTG content, DNL and gluconeogenesis. RESULTS: MLX knockdown in PHH altered lipid metabolism (decreased DNL (p < 0.05), increased fatty acid oxidation and ketogenesis (p < 0.05), and reduced lipid accumulation (p < 0.001)). Additionally, MLX knockdown increased glycolysis, lactate secretion and glucose production (p < 0.001) and insulin-stimulated pAKT levels (p < 0.01) as assessed by transcriptomic, steady-state and dynamic measurements. Consistent with the in vitro data, individuals with the rs1474040-A and rs632758-C variants had lower fasting plasma insulin (p < 0.05 and p < 0.01, respectively) and TG (p < 0.05 and p < 0.01, respectively). Although there was no difference in IHTG or gluconeogenesis, individuals with rs632758 SNP had notably lower hepatic DNL (p < 0.01). CONCLUSION: We have demonstrated using human in vitro and in vivo models that MLX inhibition favored lipid catabolism over anabolism and increased glucose production, despite increased glycolysis and phosphorylation of Akt, suggesting a metabolic mechanism that involves futile cycling.

The cholesterol-lowering effect of statins is modified by LILRB5 intolerance genotype: Results from a recruit-by-genotype clinical trial.

Background/Aims: Statin intolerance leads to poor adherence to statin therapy, resulting in a failure to achieve desired cholesterol reduction and adverse outcomes. The LILRB5 Asp247Gly genotype has been identified as being associated with statin intolerance and statin-induced myalgia. We conducted a randomized clinical trial to examine its role in immune response through T regulatory cell aggregation and in achieving cholesterol reduction targets. Methods: A double-blind, cross-over, recruit-by-genotype trial was undertaken. A total of 18 participants who had either the Asp247Asp (T/T) genotype or the Gly247Gly (C/C) genotype were recruited to the study. Participants were randomised to receive placebo or atorvastatin 80 mg daily for 28 days. Following a washout period of 3 weeks, they were then switched to the opposite treatment. Biochemical and immunological measurements as well as interviews were performed prior to and after both treatment periods. Within genotype group comparisons were performed using repeated measures Wilcoxon tests. Two-way repeated measures ANOVA with genotype and treatment as factors were used to compare changes in biochemical parameters between groups during placebo and atorvastatin periods. Results: Individuals with the Asp247Asp genotype had a greater increase in creatine kinase (CK) compared to those with Gly247Gly genotype in response to atorvastatin (p = 0.03). Those with Gly247Gly genotype had a mean non-HDL cholesterol reduction of 2.44 (95% CI:1.59 - 3.29) mmol/L while in Asp247Asp genotype group the mean reduction was 1.28 (95%CI: 0.48 - 2.07) mmol/L. The interaction between the genotype and atorvastatin treatment for total cholesterol (p = 0.007) and non-HDL cholesterol response was significant (p = 0.025). Immunological assessment showed no significant changes in aggregation of T regulatory cells by genotype. Conclusion: The Asp247Gly variant in LILRB5, previously associated with statin intolerance, was associated with differential increases in creatine kinase and total cholesterol and non-HDL cholesterol-lowering response to atorvastatin. Taken together, these results suggest that this variant could have utility in precision cardiovascular therapy.

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.

Common Statin Intolerance Variants in ABCB1 and LILRB5 Show Synergistic Effects on Statin Response: An Observational Study Using Electronic Health Records.

Background: Statin intolerance impacts approximately 10% of statin users, with side effects ranging from mild myalgia to extreme intolerance resulting in myopathy and rhabdomyolysis. Statin intolerance results in poor adherence to therapy and can impact statin efficacy. Many genetic variants are associated with statin intolerance. The effect of these variants on statin efficacy has not been systematically explored. Methods: Using longitudinal electronic health records and genetic biobank data from Tayside, Scotland, we examined the effect of seven genetic variants with previously reported associations with simvastatin or atorvastatin intolerance on the outcome of statin response. Statin response was measured by the reduction achieved when comparing pre- and post-statin non-high-density lipoprotein-cholesterol (non-HDL-C). Post-treatment statin response was limited to non-HDL-C measured within 6months of therapy initiation. Univariate and multivariable linear regression models were used to assess the main and adjusted effect of the variants on statin efficacy. Results: Around 9,401 statin users met study inclusion criteria, of whom 8,843 were first prescribed simvastatin or atorvastatin. The average difference in post-treatment compared to pre-treatment non-HDL-cholesterol was 1.45 (±1.04) mmol/L. In adjusted analyses, only two variants, one in the gene ATP-binding cassette transporter B1 (ABCB1; rs1045642), and one in leukocyte immunoglobulin like receptor B5 (LILRB5; rs12975366), were associated with statin efficacy. In ABCB1, homozygous carriers of the C allele at rs1045642 had 0.06mmol/L better absolute reduction in non-HDL-cholesterol than carriers of the T allele (95% CI: 0.01, 0.1). In LILRB5 (rs12975366), carriers of the C allele had 0.04mmol/L better absolute reduction compared to those homozygous for the T allele (95% CI: 0.004, 0.08). When combined into a two-variant risk score, individuals with both the rs1045642-CC genotype and the rs12975366-TC or CC genotype had a 0.11mmol/L greater absolute reduction in non-HDL-cholesterol compared to those with rs1045642-TC or TT genotype and the rs12975366-TT genotype (95% CI: 0.05, 0.16; p<0.001). Conclusion: We report two genetic variants for statin adverse drug reactions (ADRs) that are associated with statin efficacy. While the ABCB1 variant has been shown to have an association with statin pharmacokinetics, no similar evidence for LILRB5 has been reported. These findings highlight the value of genetic testing to deliver precision therapeutics to statin users.

LTA4H rs2660845 association with montelukast response in early and late-onset asthma.

Leukotrienes play a central pathophysiological role in both paediatric and adult asthma. However, 35% to 78% of asthmatics do not respond to leukotriene inhibitors. In this study we tested the role of the LTA4H regulatory variant rs2660845 and age of asthma onset in response to montelukast in ethnically diverse populations. We identified and genotyped 3,594 asthma patients treated with montelukast (2,514 late-onset and 1,080 early-onset) from seven cohorts (UKBiobank, GoSHARE, BREATHE, Tayside RCT, PAGES, GALA II and SAGE). Individuals under montelukast treatment experiencing at least one exacerbation in a 12-month period were compared against individuals with no exacerbation, using logistic regression for each cohort and meta-analysis. While no significant association was found with European late-onset subjects, a meta-analysis of 523 early-onset individuals from European ancestry demonstrated the odds of experiencing asthma exacerbations by carriers of at least one G allele, despite montelukast treatment, were increased (odds-ratio = 2.92, 95%confidence interval (CI): 1.04-8.18, I2 = 62%, p = 0.0412) compared to those in the AA group. When meta-analysing with other ethnic groups, no significant increased risk of asthma exacerbations was found (OR = 1.60, 95% CI: 0.61-4.19, I2 = 85%, p = 0.342). Our study demonstrates that genetic variation in LTA4H, together with timing of asthma onset, may contribute to variability in montelukast response. European individuals with early-onset (≤18y) carrying at least one copy of rs2660845 have increased odd of exacerbation under montelukast treatment, presumably due to the up-regulation of LTA4H activity. These findings support a precision medicine approach for the treatment of asthma with montelukast.

Polymorphism in INSR Locus Modifies Risk of Atrial Fibrillation in Patients on Thyroid Hormone Replacement Therapy.

AIMS: Atrial fibrillation (AF) is a risk for patients receiving thyroid hormone replacement therapy. No published work has focused on pharmacogenetics relevant to thyroid dysfunction and AF risk. We aimed to assess the effect of L-thyroxine on AF risk stratified by a variation in a candidate gene. METHODS AND RESULTS: A retrospective follow-up study was done among European Caucasian patients from the Genetics of Diabetes Audit and Research in Tayside Scotland cohort (Scotland, United Kingdom). Linked data on biochemistry, prescribing, hospital admissions, demographics, and genetic biobank were used to ascertain patients on L-thyroxine and diagnosis of AF. A GWAS-identified insulin receptor-INSR locus (rs4804416) was the candidate gene. Cox survival models and sensitivity analyses by taking competing risk of death into account were used. Replication was performed in additional sample (The Genetics of Scottish Health Research register, GoSHARE), and meta-analyses across the results of the study and replication cohorts were done. We analyzed 962 exposed to L-thyroxine and 5,840 unexposed patients who were rs4804416 genotyped. The rarer G/G genotype was present in 18% of the study population. The total follow-up was up to 20 years, and there was a significant increased AF risk for patients homozygous carriers of the G allele exposed to L-thyroxine (RHR = 2.35, P = 1.6e-02). The adjusted increased risk was highest within the first 3 years of exposure (RHR = 9.10, P = 8.5e-04). Sensitivity analysis yielded similar results. Effects were replicated in GoSHARE (n = 3,190). CONCLUSION: Homozygous G/G genotype at the INSR locus (rs4804416) is associated with an increased risk of AF in patients on L-thyroxine, independent of serum of free thyroxine and thyroid-stimulating hormone serum concentrations.

Genome-wide association study of asthma exacerbations despite inhaled corticosteroid use.

RATIONALE: Substantial variability in response to asthma treatment with inhaled corticosteroids (ICS) has been described among individuals and populations, suggesting the contribution of genetic factors. Nonetheless, only a few genes have been identified to date. We aimed to identify genetic variants associated with asthma exacerbations despite ICS use in European children and young adults and to validate the findings in non-Europeans. Moreover, we explored whether a gene-set enrichment analysis could suggest potential novel asthma therapies. METHODS: A genome-wide association study (GWAS) of asthma exacerbations was tested in 2681 children of European descent treated with ICS from eight studies. Suggestive association signals were followed up for replication in 538 European asthma patients. Further evaluation was performed in 1773 non-Europeans. Variants revealed by published GWAS were assessed for replication. Additionally, gene-set enrichment analysis focused on drugs was performed. RESULTS: 10 independent variants were associated with asthma exacerbations despite ICS treatment in the discovery phase (p≤5×10-6). Of those, one variant at the CACNA2D3-WNT5A locus was nominally replicated in Europeans (rs67026078; p=0.010), but this was not validated in non-European populations. Five other genes associated with ICS response in previous studies were replicated. Additionally, an enrichment of associations in genes regulated by trichostatin A treatment was found. CONCLUSIONS: The intergenic region of CACNA2D3 and WNT5A was revealed as a novel locus for asthma exacerbations despite ICS treatment in European populations. Genes associated were related to trichostatin A, suggesting that this drug could regulate the molecular mechanisms involved in treatment response.

Exome Sequencing Reveals Common and Rare Variants in F5 Associated With ACE Inhibitor and Angiotensin Receptor Blocker-Induced Angioedema.

Angioedema occurring in the head and neck region is a rare and sometimes life-threatening adverse reaction to angiotensin-converting enzyme inhibitors (ACEIs) and angiotensin receptor blockers (ARBs). Few studies have investigated the association of common variants with this extreme reaction, but none have explored the combined influence of rare variants yet. Adjudicated cases of ACEI-induced angioedema (ACEI-AE) or ARB-induced angioedema (ARB-AE) and controls were recruited at five different centers. Sequencing of 1,066 samples (408 ACEI-AE, ARB-AE, and 658 controls) was performed using exome-enriched sequence data. A common variant of the F5 gene that causes an increase in blood clotting (rs6025, p.Arg506Gln, also called factor V Leiden), was significantly associated with both ACEI-AE and ARB-AE (odds ratio: 2.85, 95% confidence interval (CI), 1.89-4.25). A burden test analysis of five rare missense variants in F5 was also found to be associated with ACEI-AE or ARB-AE, P = 2.09 × 10-3 . A combined gene risk score of these variants, and the common variants rs6025 and rs6020, showed that individuals carrying at least one variant had 2.21 (95% CI, 1.49-3.27, P = 6.30 × 10-9 ) times the odds of having ACEI-AE or ARB-AE. The increased risk due to the common Leiden allele was confirmed in a genome-wide association study from the United States. A high risk of angioedema was also observed for the rs6020 variant that is the main coagulation defect-causing variant in black African and Asian populations. We found that deleterious missense variants in F5 are associated with an increased risk of ACEI-AE or ARB-AE.

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.

Author Correction: A genome scan for milk production traits in dairy goats reveals two new mutations in Dgat1 reducing milk fat content.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Genome-wide association mapping for type and mammary health traits in French dairy goats identifies a pleiotropic region on chromosome 19 in the Saanen breed.

Type traits and mammary health traits are important to dairy ruminant breeding because they influence animal health, milking ability, and longevity, as well as the economic sustainability of farms. The availability of the genomic sequence and a single nucleotide polymorphism chip in goats has opened up new fields of investigation to better understand the genes and mechanisms that underlie such complex traits and to be able to select them. Our objective was to perform a genome-wide association study in dairy goats for 11 type traits and somatic cell count (SCC) as proxies for mastitis resistance. A genome-wide association study was implemented using a daughter design composed of 1,941 Alpine and Saanen goats sired by 20 artificial insemination bucks, genotyped with the Illumina GoatSNP50 BeadChip (Illumina Inc., San Diego, CA). This association study was based on both linkage analyses and linkage disequilibrium using QTLmap software (http://dga7.jouy.inra.fr/qtlmap/) interval mapping was performed with the likelihood ratio test using linear regressions. Breeds were analyzed together and separately. The study highlighted 37 chromosome-wide significant quantitative trait loci (QTL) with linkage analyses and 222 genome-wide significant QTL for linkage disequilibrium, for type and SCC traits in dairy goats. Genomic control of those traits was mostly polygenic and breed-specific, suggesting that within-breed selection would be favored for those traits. Of note, Capra hircus autosome (CHI) 19 appeared to be highly enriched in single nucleotide polymorphisms associated with type and SCC, with 2 highly significant regions in the Saanen breed. One region (33-42 Mb) was significantly associated with SCC and includes candidate genes associated with response to intramammary infections (RARA, STAT3, STAT5A, and STAT5B). Another region of the CHI 19 (24.5-27 Mb) exhibited an adverse pleiotropic effect on milk production (milk, fat yield, and protein yield) and udder traits (udder floor position and rear udder attachment) that agreed with the negative genetic correlations that exist between those 2 groups of traits. These QTL were not found in the Alpine breed. In Alpine, the 2 most significant regions were associated with chest depth on CHI 6 (45.8-46.0 Mb) and CHI 8 (80.7-81.1 Mb). These results will be helpful for goat selection in the future and could lead to identification of causal mutations.

A common missense variant of LILRB5 is associated with statin intolerance and myalgia.

Aims: A genetic variant in LILRB5 (leukocyte immunoglobulin-like receptor subfamily-B) (rs12975366: T > C: Asp247Gly) has been reported to be associated with lower creatine phosphokinase (CK) and lactate dehydrogenase (LDH) levels. Both biomarkers are released from injured muscle tissue, making this variant a potential candidate for susceptibility to muscle-related symptoms. We examined the association of this variant with statin intolerance ascertained from electronic medical records in the GoDARTS study. Methods and results: In the GoDARTS cohort, the LILRB5 Asp247 variant was associated with statin intolerance (SI) phenotypes; one defined as having raised CK and being non-adherent to therapy [odds ratio (OR) 1.81; 95% confidence interval (CI): 1.34-2.45] and the other as being intolerant to the lowest approved dose of a statin before being switched to two or more other statins (OR 1.36; 95% CI: 1.07-1.73). Those homozygous for Asp247 had increased odds of developing both definitions of intolerance. Importantly the second definition did not rely on CK elevations. These results were replicated in adjudicated cases of statin-induced myopathy in the PREDICTION-ADR consortium (OR1.48; 95% CI: 1.05-2.10) and for the development of myalgia in the JUPITER randomized clinical trial of rosuvastatin (OR1.35, 95% CI: 1.10-1.68). A meta-analysis across the studies showed a consistent association between Asp247Gly and outcomes associated with SI (OR1.34; 95% CI: 1.16-1.54). Conclusion: This study presents a novel immunogenetic factor associated with statin intolerance, an important risk factor for cardiovascular outcomes. The results suggest that true statin-induced myalgia and non-specific myalgia are distinct, with a potential role for the immune system in their development. We identify a genetic group that is more likely to be intolerant to their statins.

CKM Glu83Gly Is Associated With Blunted Creatine Kinase Variation, but Not With Myalgia.

BACKGROUND: To test the association of a recently reported variant in the creatine kinase (CK) muscle gene, CKM Glu83Gly (rs11559024) with constitutive creatine phosphokinase (CK) levels, CK variation, and inducibility. Given the diagnostic importance of CK in determining muscle damage, we tested the association of the variant with myalgia. METHODS AND RESULTS: Meta-analysis between longitudinal cohort GoDARTS (Genetics of Diabetes Audit and Research, Tayside Scotland), minor allele frequency (=0.02), and randomized clinical trial (JUPITER [Justification for the Use of Statins in Primary Prevention: An Intervention Trial Evaluating Rosuvastatin], minor allele frequency=0.018) was used to replicate the association with baseline CK measures. GoDARTS was used to study the relationship with CK variability. Myalgia was studied in JUPITER trial participants. Baseline and SDs of CK were on average 18% (P value=6×10-63) and 24% (P value=2×10-5) lower for carriers of the variant, respectively. The variant was not associated with myalgia (odds ratio, 0.84; 95% confidence interval, 0.52-1.38). CONCLUSIONS: This study highlights that a genetic factor known to be associated with constitutive CK levels is also associated with CK variability and inducibility. This is discussed in the context of evidence to suggest that the variant has an impact on inducibility of CK by trauma through a previously reported case of a homozygous carrier. However, the lack of association between the variant and myalgia suggests that it cannot reliably be used as a biomarker for muscle symptoms.

A genome scan for milk production traits in dairy goats reveals two new mutations in Dgat1 reducing milk fat content.

The quantity of milk and milk fat and proteins are particularly important traits in dairy livestock. However, little is known about the regions of the genome that influence these traits in goats. We conducted a genome wide association study in French goats and identified 109 regions associated with dairy traits. For a major region on chromosome 14 closely associated with fat content, the Diacylglycerol O-Acyltransferase 1 (DGAT1) gene turned out to be a functional and positional candidate gene. The caprine reference sequence of this gene was completed and 29 polymorphisms were found in the gene sequence, including two novel exonic mutations: R251L and R396W, leading to substitutions in the protein sequence. The R251L mutation was found in the Saanen breed at a frequency of 3.5% and the R396W mutation both in the Saanen and Alpine breeds at a frequencies of 13% and 7% respectively. The R396W mutation explained 46% of the genetic variance of the trait, and the R251L mutation 6%. Both mutations were associated with a notable decrease in milk fat content. Their causality was then demonstrated by a functional test. These results provide new knowledge on the genetic basis of milk synthesis and will help improve the management of the French dairy goat breeding program.