Marisa W. Medina

Variation in the 3-Hydroxyl-3-Methylglutaryl Coenzyme A Reductase Gene Is Associated With Racial Differences in Low-Density Lipoprotein Cholesterol Response to Simvastatin Treatment

Background— Use of 3-hydroxyl-3-methylglutaryl-3 coenzyme A reductase (HMGCR) inhibitors, or statins, reduces cardiovascular disease risk by lowering plasma low-density lipoprotein cholesterol (LDL-C) concentrations. However, LDL-C response is variable and influenced by many factors, including racial ancestry, with attenuated response in blacks compared with whites. We hypothesized that single nucleotide polymorphisms in the gene encoding HMGCR, a rate-limiting enzyme in cholesterol synthesis and the direct enzymatic target of statins, contribute to variation in statin response. Methods and Results— Genomic resequencing of HMGCR in 24 blacks and 23 whites identified 79 single nucleotide polymorphisms. Eleven single nucleotide polymorphisms were selected to tag common linkage disequilibrium clusters. These single nucleotide polymorphisms and the common haplotypes inferred from them were tested for association with plasma LDL-C and LDL-C response to simvastatin treatment (40 mg/d for 6 weeks) in 326 blacks and 596 whites. Black carriers of H7 and/or H2 had significantly lower baseline LDL-C (P=0.0006) and significantly attenuated LDL-C response compared with black participants who did not carry either haplotype as measured by absolute response (−1.23±0.04 mmol/L, n=209, versus −1.45±0.06 mmol/L, n=117; P=0.0008) and percent response (−36.9±1.0% versus −40.6±1.3%; P=0.02), but no haplotype effect was observed in whites. Percent LDL-C response was lowest in carriers of both H2 and H7, all but one of whom were black (−28.2±4.9%, n=12 H2+H7 carriers, versus −41.5±0.5%, n=650 H2/H7 noncarriers; P=0.001). LDL-C responses in H7 and/or H2 noncarriers were indistinguishable between blacks and whites. Conclusions— HMGCR gene polymorphisms are associated with reduced plasma LDL-C and LDL-C response to simvastatin, and these effects are most evident in blacks.

A statin-dependent QTL for GATM expression is associated with statin-induced myopathy.

Statins are prescribed widely to lower plasma low-density lipoprotein (LDL) concentrations and cardiovascular disease risk and have been shown to have beneficial effects in a broad range of patients. However, statins are associated with an increased risk, albeit small, of clinical myopathy and type 2 diabetes. Despite evidence for substantial genetic influence on LDL concentrations, pharmacogenomic trials have failed to identify genetic variations with large effects on either statin efficacy or toxicity, and have produced little information regarding mechanisms that modulate statin response. Here we identify a downstream target of statin treatment by screening for the effects of in vitro statin exposure on genetic associations with gene expression levels in lymphoblastoid cell lines derived from 480 participants of a clinical trial of simvastatin treatment. This analysis identified six expression quantitative trait loci (eQTLs) that interacted with simvastatin exposure, including rs9806699, a cis-eQTL for the gene glycine amidinotransferase (GATM) that encodes the rate-limiting enzyme in creatine synthesis. We found this locus to be associated with incidence of statin-induced myotoxicity in two separate populations (meta-analysis odds ratio = 0.60). Furthermore, we found that GATM knockdown in hepatocyte-derived cell lines attenuated transcriptional response to sterol depletion, demonstrating that GATM may act as a functional link between statin-mediated lowering of cholesterol and susceptibility to statin-induced myopathy.

Alternative splicing of 3-hydroxy-3-methylglutaryl coenzyme A reductase is associated with plasma low-density lipoprotein cholesterol response to simvastatin.

Background— HMGCR(3-Hydroxy-3-methylglutaryl coenzyme A reductase), the direct target of statin inhibition, undergoes alternative splicing of exon 13, which encodes part of the statin-binding domain of the enzyme. We hypothesized that HMGCR alternative splicing might be related to the interindividual variation in plasma low-density lipoprotein cholesterol response to statin treatment. Methods and Results— We measured mRNA expression of both the full-length and the alternatively spliced HMGCR transcript lacking exon 13 (HMGCRv_1) in 170 simvastatin-incubated immortalized lymphocyte cell lines derived from participants in the Cholesterol and Pharmacogenetics (CAP) study who were treated with simvastatin 40 mg/d for 6 weeks. Greater upregulation of HMGCRv_1 in vitro was significantly correlated (P≤0.0001) with smaller in vivo reductions of plasma total cholesterol, low-density lipoprotein cholesterol, apoprotein B, and triglycerides and explained 6% to 15% of the variation in their response to treatment. In contrast, no significant relationship was found between expression of the full-length HMGCR transcript and in vivo response. By siRNA knockdown of the full-length transcript, we found that HMGCR enzyme activity measured in cells enriched in HMGCRv_1 was relatively resistant to statin inhibition, consistent with the association of increased alternative splicing with reduced statin response in the CAP study. In addition, we found that a common HMGCR single-nucleotide polymorphism (rs3846662) located within intron 13 was associated with variation in the proportion of HMGCR mRNA that is alternatively spliced. Conclusions— Variation in the production of an HMGCR isoform with reduced statin sensitivity is a determinant of interindividual differences in low-density lipoprotein cholesterol, apolipoprotein B, and triglyceride response to statin treatment.

Pharmacogenomic discovery using cell-based models.

Quantitative variation in response to drugs in human populations is multifactorial; genetic factors probably contribute to a significant extent. Identification of the genetic contribution to drug response typically comes from clinical observations and use of classic genetic tools. These clinical studies are limited by our inability to control environmental factors in vivo and the difficulty of manipulating the in vivo system to evaluate biological changes. Recent progress in dissecting genetic contribution to natural variation in drug response through the use of cell lines has been made and is the focus of this review. A general overview of current cell-based models used in pharmacogenomic discovery and validation is included. Discussion includes the current approach to translate findings generated from these cell-based models into the clinical arena and the use of cell lines for functional studies. Specific emphasis is given to recent advances emerging from cell line panels, including the International HapMap Project and the NCI60 cell panel. These panels provide a key resource of publicly available genotypic, expression, and phenotypic data while allowing researchers to generate their own data related to drug treatment to identify genetic variation of interest. Interindividual and interpopulation differences can be evaluated because human lymphoblastoid cell lines are available from major world populations of European, African, Chinese, and Japanese ancestry. The primary focus is recent progress in the pharmacogenomic discovery area through ex vivo models.

Coordinately Regulated Alternative Splicing of Genes Involved in Cholesterol Biosynthesis and Uptake

Genes involved in cholesterol biosynthesis and uptake are transcriptionally regulated in response to cellular sterol content in a coordinated manner. A number of these genes, including 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGCR) and LDL receptor (LDLR), undergo alternative splicing, resulting in reductions of enzyme or protein activity. Here we demonstrate that cellular sterol depletion suppresses, and sterol loading induces, alternative splicing of multiple genes involved in the maintenance of cholesterol homeostasis including HMGCR and LDLR, the key regulators of cellular cholesterol biosynthesis and uptake, respectively. These changes were observed in both in vitro studies of the HepG2 human hepatoma derived cell line, as well as in vivo studies of St. Kitts vervets, also known as African green monkeys, a commonly used primate model for investigating cholesterol metabolism. These effects are mediated in part by sterol regulation of polypyrimidine tract binding protein 1 (PTBP1), since knock-down of PTBP1 eliminates sterol induced changes in alternative splicing of several of these genes. Single nucleotide polymorphisms (SNPs) that influence HMGCR and LDLR alternative splicing (rs3846662 and rs688, respectively), have been associated with variation in plasma LDL-cholesterol levels. Sterol-induced changes in alternative splicing are blunted in carriers of the minor alleles for each of these SNPs, indicating an interaction between genetic and non-genetic regulation of this process. Our results implicate alternative splicing as a novel mechanism of enhancing the robust transcriptional response to conditions of cellular cholesterol depletion or accumulation. Thus coordinated regulation of alternative splicing may contribute to cellular cholesterol homeostasis as well as plasma LDL levels.

Combined Influence of LDLR and HMGCR Sequence Variation on Lipid-Lowering Response to Simvastatin

Objective—Although statins are efficacious for lowering low-density lipoprotein cholesterol, there is wide interindividual variation in response. We tested the extent to which combined effects of common alleles of LDLR and HMGCR can contribute to this variability. Methods and Results—Haplotypes in the LDLR 3′-untranslated region (3-UTR) were tested for association with lipid-lowering response to simvastatin treatment in the Cholesterol and Pharmacogenetics trial (335 blacks and 609 whites). LDLR haplotype 5 (LDLR L5) was associated with smaller simvastatin-induced reductions in low-density lipoprotein cholesterol, total cholesterol, non-high-density lipoprotein cholesterol, and apolipoprotein B (P=0.0002 to 0.03) in blacks but not whites. The combined presence of LDLR L5 and previously described HMGCR haplotypes in blacks was associated with significantly attenuated apolipoprotein B reduction (−22.4±1.5%, N=89) compared with both noncarriers (−30.6±1.5%, N=78, P=0.0001) and carriers of either individual haplotype (−28.2±1.1%, N=158, P=0.001). We observed similar differences when measuring simvastatin-mediated induction of low-density lipoprotein receptor surface expression using lymphoblast cell lines (P=0.03). Conclusion—We have identified a common LDLR 3-UTR haplotype that is associated with attenuated lipid-lowering response to simvastatin treatment. Response was further reduced in individuals with both LDLR and previously described HMGCR haplotypes. Previously identified racial differences in statin efficacy were partially explained by the greater prevalence of these combined haplotypes in blacks.

The role of HMGCR alternative splicing in statin efficacy

Statins, or 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGCR) inhibitors, are widely prescribed to lower plasma cholesterol levels and reduce cardiovascular disease risk. Despite the well-documented efficacy of statins, there is large interindividual variation in response. Using a panel of immortalized lymphocyte cell lines incubated with simvastatin, we recently found that the magnitude of expression of an alternatively spliced HMGCR transcript lacking exon 13 was inversely correlated with in vivo reductions of total cholesterol, low-density lipoprotein cholesterol, apoB, and triglycerides after statin treatment of the individuals from whom the cells were derived. This review will discuss the potential significance of alternative splicing as a mechanism contributing to variation in statin efficacy as well as the use of immortalized lymphocyte cell lines for identifying pharmacogenetically relevant polymorphisms and molecular mechanisms.

HNRNPA1 regulates HMGCR alternative splicing and modulates cellular cholesterol metabolism

3-hydroxy-3-methylglutaryl-Coenzyme A reductase (HMGCR) encodes the rate-limiting enzyme in the cholesterol biosynthesis pathway and is inhibited by statins, a class of cholesterol-lowering drugs. Expression of an alternatively spliced HMGCR transcript lacking exon 13, HMGCR13(-), has been implicated in the variation of plasma LDL-cholesterol (LDL-C) and is the single most informative molecular marker of LDL-C response to statins. Given the physiological importance of this transcript, our goal was to identify molecules that regulate HMGCR alternative splicing. We recently reported gene expression changes in 480 lymphoblastoid cell lines (LCLs) after in vitro simvastatin treatment, and identified a number of statin-responsive genes involved in mRNA splicing. Heterogeneous nuclear ribonucleoprotein A1 (HNRNPA1) was chosen for follow-up since rs3846662, an HMGCR SNP that regulates exon 13 skipping, was predicted to alter an HNRNPA1 binding motif. Here, we not only demonstrate that rs3846662 modulates HNRNPA1 binding, but also that sterol depletion of human hepatoma cell lines reduced HNRNPA1 mRNA levels, an effect that was reversed with sterol add-back. Overexpression of HNRNPA1 increased the ratio of HMGCR13(-) to total HMGCR transcripts by both directly increasing exon 13 skipping in an allele-related manner and specifically stabilizing the HMGCR13(-) transcript. Importantly, HNRNPA1 overexpression also diminished HMGCR enzyme activity, enhanced LDL-C uptake and increased cellular apolipoprotein B (APOB). rs1920045, an SNP associated with HNRNPA1 exon 8 alternative splicing, was also associated with smaller statin-induced reduction in total cholesterol from two independent clinical trials. These results suggest that HNRNPA1 plays a role in the variation of cardiovascular disease risk and statin response.

A genome-wide association analysis of temozolomide response using lymphoblastoid cell lines shows a clinically relevant association with MGMT.

Objective Recently, lymphoblastoid cell lines (LCLs) have emerged as an innovative model system for mapping gene variants that predict the dose response to chemotherapy drugs. Methods In the current study, this strategy was expanded to the in-vitro genome-wide association approach, using 516 LCLs derived from a White cohort to assess the cytotoxic response to temozolomide. Results Genome-wide association analysis using ∼2.1 million quality-controlled single-nucleotide polymorphisms (SNPs) identified a statistically significant association (P<10−8) with SNPs in the O6-methylguanine-DNA methyltransferase (MGMT) gene. We also show that the primary SNP in this region is significantly associated with the differential gene expression of MGMT (P<10–26) in LCLs and differential methylation in glioblastoma samples from The Cancer Genome Atlas. Conclusion The previously documented clinical and functional relationships between MGMT and temozolomide response highlight the potential of well-powered genome-wide association studies of the LCL model system to identify meaningful genetic associations.

RHOA Is a Modulator of the Cholesterol-Lowering Effects of Statin

Although statin drugs are generally efficacious for lowering plasma LDL-cholesterol levels, there is considerable variability in response. To identify candidate genes that may contribute to this variation, we used an unbiased genome-wide filter approach that was applied to 10,149 genes expressed in immortalized lymphoblastoid cell lines (LCLs) derived from 480 participants of the Cholesterol and Pharmacogenomics (CAP) clinical trial of simvastatin. The criteria for identification of candidates included genes whose statin-induced changes in expression were correlated with change in expression of HMGCR, a key regulator of cellular cholesterol metabolism and the target of statin inhibition. This analysis yielded 45 genes, from which RHOA was selected for follow-up because it has been found to participate in mediating the pleiotropic but not the lipid-lowering effects of statin treatment. RHOA knock-down in hepatoma cell lines reduced HMGCR, LDLR, and SREBF2 mRNA expression and increased intracellular cholesterol ester content as well as apolipoprotein B (APOB) concentrations in the conditioned media. Furthermore, inter-individual variation in statin-induced RHOA mRNA expression measured in vitro in CAP LCLs was correlated with the changes in plasma total cholesterol, LDL-cholesterol, and APOB induced by simvastatin treatment (40 mg/d for 6 wk) of the individuals from whom these cell lines were derived. Moreover, the minor allele of rs11716445, a SNP located in a novel cryptic RHOA exon, dramatically increased inclusion of the exon in RHOA transcripts during splicing and was associated with a smaller LDL-cholesterol reduction in response to statin treatment in 1,886 participants from the CAP and Pravastatin Inflamation and CRP Evaluation (PRINCE; pravastatin 40 mg/d) statin clinical trials. Thus, an unbiased filter approach based on transcriptome-wide profiling identified RHOA as a gene contributing to variation in LDL-cholesterol response to statin, illustrating the power of this approach for identifying candidate genes involved in drug response phenotypes.