Donate Help Contact The AHA Sign In Home
American Heart Association
Arteriosclerosis, Thrombosis, and Vascular Biology
Search: search_blue_button Advanced Search
Arteriosclerosis, Thrombosis, and Vascular Biology. 2008;28:2078-2084
Published online before print September 18, 2008, doi: 10.1161/ATVBAHA.108.172288
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Data Supplement
Right arrow All Versions of this Article:
28/11/2078    most recent
ATVBAHA.108.172288v1
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowRequest Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Burkhardt, R.
Right arrow Articles by Breslow, J. L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Burkhardt, R.
Right arrow Articles by Breslow, J. L.
Related Collections
Right arrow Functional genomics
Right arrow Gene expression
Right arrow Gene regulation
Right arrow Genomics
Right arrow Genetics of cardiovascular disease
Right arrow Lipid and lipoprotein metabolism
(Arteriosclerosis, Thrombosis, and Vascular Biology. 2008;28:2078.)
© 2008 American Heart Association, Inc.


Clinical and Population Studies

Common SNPs in HMGCR in Micronesians and Whites Associated With LDL-Cholesterol Levels Affect Alternative Splicing of Exon13

Ralph Burkhardt; Eimear E. Kenny; Jennifer K. Lowe; Andrew Birkeland; Rebecca Josowitz; Martha Noel; Jacqueline Salit; Julian B. Maller; Itsik Pe'er; Mark J. Daly; David Altshuler; Markus Stoffel; Jeffrey M. Friedman; Jan L. Breslow

From the Laboratory of Biochemical Genetics and Metabolism (R.B., E.E.K., J.K.L., A.B., R.J., M.N., J.L.B.), the Laboratory of Molecular Genetics (J.S., J.M.F.), Howard Hughes Medical Institute (J.M.F.), The Rockefeller University, New York; the Broad Institute of Harvard and MIT (J.K.L., J.B.M., M.J.D., D.A.), Cambridge, Mass; the Department of Computer Science (I.P.), Columbia University, New York; the Center for Human Genetics Research (J.B.M., M.J.D.) and the Department of Molecular Biology (D.A., J.K.L.), Massachusetts General Hospital, Boston; the Department of Medicine (M.J.D., D.A.) and the Department of Genetics (D.A.), Harvard Medical School, Boston, Mass; and the Institute of Molecular Systems Biology (M.S.), ETH Zurich, Switzerland.

Correspondence to Jan L. Breslow, Laboratory of Biochemical Genetics, The Rockefeller University, 1230 York Avenue, New York, NY 10065. E-mail breslow{at}rockefeller.edu


*    Abstract
up arrowTop
*Abstract
down arrowIntroduction
down arrowMethods
down arrowResults
down arrowDiscussion
down arrowReferences
 
Background— Variation in LDL-cholesterol (LDL-C) among individuals is a complex genetic trait involving multiple genes and gene-environment interactions.

Methods and Results— In a genome-wide association study (GWAS) to identify genetic variants influencing LDL-C in an isolated population from Kosrae, we observed associations for SNPs in the gene encoding 3hydroxy-3-methylglutaryl (HMG)-coenzyme A (CoA) reductase (HMGCR). Three of these SNPs (rs7703051, rs12654264, and rs3846663) met the statistical threshold of genome-wide significance when combined with data from the Diabetes Genetics Initiative GWAS. We followed up the association results and identified a functional SNP in intron13 (rs3846662), which was in linkage disequilibrium with the SNPs of genome-wide significance and affected alternative splicing of HMGCR mRNA. In vitro studies in human lymphoblastoid cells demonstrated that homozygosity for the rs3846662 minor allele was associated with up to 2.2-fold lower expression of alternatively spliced HMGCR mRNA lacking exon13, and minigene transfection assays confirmed that allele status at rs3846662 directly modulated alternative splicing of HMGCR exon13 (42.9±3.9 versus 63.7±1.0%{Delta}exon13/total HMGCR mRNA, P=0.02). Further, the alternative splice variant could not restore HMGCR activity when expressed in HMGCR deficient UT-2 cells.

Conclusion— We identified variants in HMGCR that are associated with LDL-C across populations and affect alternative splicing of HMGCR exon13.

In a genome-wide association study we identified variants in the HMG-CoA reductase gene that were associated with LDL-C across populations. Further, functional studies revealed that a common intronic SNP (rs384662) that is in linkage disequilibrium with the variants typed in the genome scan modulates alternative splicing of HMGCR mRNA.


Key Words: HMG-CoA reductase • SNP • genome-wide association study • LDL-C • alternative splicing


*    Introduction
up arrowTop
up arrowAbstract
*Introduction
down arrowMethods
down arrowResults
down arrowDiscussion
down arrowReferences
 
Elevated levels of LDL-cholesterol (LDL-C) are a primary risk factor for atherosclerotic cardiovascular disease, the major cause of morbidity and mortality in industrialized countries today.1 Variation in LDL-C among individuals is a complex genetic trait, involving multiple genes and significant gene-environment interactions.2 Candidate gene and linkage studies have identified some of the genetic factors contributing to the population variance in plasma lipoprotein levels,3,4 but these factors only explain a small fraction of the heritability, suggesting that additional variants influencing lipid levels remain to be identified.

Our group has previously used candidate gene5 and linkage6 approaches to identify genetic loci affecting plasma lipid and lipoprotein levels in a cohort from the Island of Kosrae, Federated States of Micronesia, a genetic isolate with significant founder effects and a high prevalence of traits related to the metabolic syndrome.7

Recently, genome-wide association studies (GWAS) have been shown to be successful in gene discovery for complex traits and offer a new approach to identify common genetic variants with modest effects.8 Using Affymetrix gene chip 500k arrays, we have performed a GWAS in {approx}2400 Kosraens for LDL-C and other metabolic traits (Lowe et al, submitted). The strongest association for LDL-C was found for 2 SNPs on chromosome 19q13 near APOE, a candidate gene with known common coding polymorphisms.9 The second best hits for LDL-C were SNPs that mapped to the HMG-CoA reductase (HMGCR) gene, the rate-limiting enzyme in cholesterol biosynthesis, and the target of LDL-C lowering statin drugs. During the preparation of this manuscript, 4 GWAS for LDL-C in white cohorts were published,10–13 and 1 of them13 included an association between the same SNPs in HMGCR and plasma LDL-C in whites.

Identifying a strong and replicable association signal in a GWAS is just the first step in elucidating the specific genetic variants involved in predisposing individuals to complex traits. In most cases, the underlying causal variant is not directly captured on the SNP array. Rather, one or more SNPs on the array are acting as a proxy for a functional nongenotyped SNP with which it is in linkage disequilibrium. Fine-mapping, resequencing, and hypothesis-driven approaches have been proposed to unearth the actual causal variants.

In this article, we report associations of SNPs in HMGCR with LDL-C in Kosraens that, in combination with similar findings from studies in whites, indicate that the same genetic variants at HMGCR contribute to differences in LDL-C across populations. To follow-up the association signal, we have implemented a hypothesis driven strategy and performed in vitro studies to identify a functional variant in intron13 of HMGCR that affects alternative splicing of exon13.


*    Methods
up arrowTop
up arrowAbstract
up arrowIntroduction
*Methods
down arrowResults
down arrowDiscussion
down arrowReferences
 
An expanded methods section is provided in the supplemental Methods (please see http://atvb.ahajournals.org).

Kosrae Study Subjects and GWAS
A genome-wide association study for LDL-C was carried out in 2346 people of the Island of Kosrae, Federated States of Micronesia using SNPs from the Affymetrix 500k platform. Details of this study will be described elsewhere (Lowe et al, submitted). All participants in the study provided written informed consent, and IRB approval was obtained from all participating institutions.

We used publicly available data from Saxena et al14 to validate our findings. The probability values of the Kosrae and DGI study analyses were combined using Fisher’s method15 to quantify the overall evidence for association. We set 5x10–8 as threshold for genome-wide statistical significance of the combined probability values.

Cell Lines
Lymphoblastoid cell lines from 18 white CEU individuals of the HapMap collection were obtained from the Coriell Institute for Medical Research. These individuals were homozygous for the major (n=9) or minor allele (n=9) at rs3846662 and the three proxy SNPs from the Affymetrix 500k array (rs3846663, rs7703051 and rs12654264). For HMGCR mRNA expression studies cells were seeded at a density of 250 000 cells/mL medium.

UT-2 cells (a gift from Drs Russell DeBose-Boyd and Michael Brown, UT Southwestern Medical Center, Dallas, Tx) are mutant Chinese hamster ovary (CHO) cells that lack HMG-CoA reductase and require mevalonate for growth. Stock cultures of UT-2 cells were grown in F12:MEM, 10% FCS, 0.2 mmol/L mevalonate and growth experiments were carried out as described.16

Quantitative Real-Time Polymerase Chain Reaction
RNA was extracted using TRIzol reagent (Invitrogen) and cDNA was synthesized using SuperScript III (Invitrogen). Quantitative RT-PCR was performed in an ABI PRISM 7700 Sequence Detector (Applied Biosystems). Specific primers and probes for full-length HMGCR, {Delta}exon13 HMGCR, and β-actin mRNA were selected to span exon junctions to avoid coamplification of genomic DNA. mRNA expression levels were normalized to β-actin as a housekeeping gene.

Minigene Constructs
Minigenes were used to assess the influence of SNP rs3846662 on alternative splicing of HMGCR exon13. PCR amplified genomic DNA for both alleles (rs3846662/A and rs3846662/G), containing exons 12–14, internal introns, and parts of the flanking introns, was cloned into the exon-trapping vector pSPL3 (kindly provided by Drs Woohyun Yoon and David B. Goldstein, Duke University, Durham, NC). Identities of the minigenes were confirmed by DNA sequencing. To confirm the causal polymorphism of the observed effects in the minigene systems, we converted SNP rs3846662 in the major form construct (allele A) to the minor form (allele G) by site directed mutagenesis (QuikChange Lightning Kit, Strategene).

Cell Transfection
HEK293 cells were transfected with HMGCR pSPL3 minigenes (A, G and A->G) and empty pSPL3 vector (negative control) using FuGene6 reagent (Roche Applied Bioscience). After 24 hours, RNA was isolated and reverse transcribed using the SA2 Primer (5'-ATCTCAGTGGTATTTGTGAGC-3'), corresponding to a transcribed exonic sequence in the pSPL3 vector and thus allowing analysis of only vector-specific HMGCR transcripts. HMGCR splicing pattern was analyzed by real-time PCR as described above.

Stable Expressing Cell Lines
The open reading frames of HMGCR full-length and {Delta}exon13 mRNA were PCR amplified and cloned into the pcDNA3.1 expression vector (Invitrogen). UT-2 cells with stable expression of human full-length (UT-2+FL) or {Delta}exon13 HMGCR (UT-2+ex13) were generated by G418 selection of FuGene6 transfected UT-2 cells.

Other Assays
Cell protein was measured using the BioRad DC kit. HMG-CoA reductase activity was measured in detergent solubilized cell extracts as described,17 except that mevalonolactone was separated by ion exchange chromatography.18


*    Results
up arrowTop
up arrowAbstract
up arrowIntroduction
up arrowMethods
*Results
down arrowDiscussion
down arrowReferences
 
Common Variants in HMGCR Are Associated With Plasma LDL-C Across Populations
We performed genome-wide association analysis for LDL-C using 2346 individuals from the Micronesian Island of Kosrae. The strongest association was found for SNPs on chromosome 19q13 in the APOE/C1/C4/C2 gene cluster (rs4420638, P=1.89x10–7). The effect of the APOE polymorphisms on LDL-C is well established in many ethnicities, and our result thus implicates that the same genetic factors are important in the Kosraen population.

The second best locus for LDL-C mapped to a region on chromosome 5q13 containing the HMGCR gene (Table). This gene encodes HMG-CoA reductase, the rate-limiting enzyme in cholesterol biosynthesis and thus represents an interesting candidate gene with high plausibility. The difference in LDL-C between homozygotes at this locus was 0.30 mmol/L (11.6 mg/dL), and the fraction of the population variance for LDL-C explained by this locus was 2.1%.


View this table:
[in this window]
[in a new window]

 
Table. Primary Associations and Meta-Analysis Results for LDL-C

However, the results did not surpass genome-wide significance when the Bonferroni correction was applied for multiple testing (most associated SNP rs3846663: P=1.28x10–6). To validate our findings in an independent cohort, we combined data from our analysis with public data for LDL-C of the Diabetes Genetics Initiative (DGI) GWAS. Because Kosraens (Micronesian) and DGI participants (European whites) are of different ancestry, we first investigated linkage disequilibrium (LD) patterns in a 1-Mb region around the locus. As shown in supplemental Figure I, pairwise LD (r2) was symmetrical between HapMap samples of European ancestry and Micronesians from Kosrae.

Combining probability values across both studies, we validated multiple SNPs at the HMGCR locus, 3 of which surpassed a genome-wide significance of P<5x10–8 (rs7703051, rs12654264, and rs3846663; Table and supplemental Figure II). A regional association plot for the combined results showed a peak of association signal over a 47-kb region containing the HMGCR gene (supplemental Figure II). In both study populations, the minor allele frequencies were comparable and the minor alleles were associated with an increase in LDL-C.

We also combined the probability values of the association analysis results for plasma total cholesterol from the Kosrae and DGI studies and found genome-wide significance for the same three SNPs at the HMGCR locus (supplemental Table I).

Association Between SNP rs3846662 and Alternative Splicing of HMGCR Exon13 in Human Lymphoblastoid Cell Lines
To follow-up the association results, we next aimed to discover functional variants at the HMGCR locus and study their molecular mode of action. We used existing resequencing data of the region containing the entire HMGCR gene from 23 whites—estimated to have >99% power to detect variants with a minor allele frequency of >5%19—to identify candidate functional SNPs. Under our hypothesis-driven model we focused on SNPs that would have strong potential for changing HMGCR function or levels.

All associated SNPs of the Kosrae and DGI studies were noncoding variants and in LD (r2 0.81 to 1), suggesting they represent the same association. The lack of LD (r2≤0.02) between these SNPs and the only known nonsynonymous SNP in HMGCR, rs5908 (I638V) in exon15, suggested that the association was not attributable to this protein coding mutation (Figure 1). Because the existence of a second HMGCR mRNA transcript resulting from alternative splicing had been reported in humans,20 we looked for SNPs in the vicinity of exon-intron borders. We detected that SNP rs3846662 was located 47 bp downstream of exon13 and in LD with the genotyped variants of genome-wide significance (r2:0.82 to 0.93, Figure 1). We hypothesized that this intronic SNP may be functional and modulate the splicing efficiency of exon13.


Figure 1
View larger version (45K):
[in this window]
[in a new window]

 
Figure 1. LD plot (r2 values) in HapMap CEU individuals for SNPs at the HMGCR gene locus demonstrating LD between Affymetrix 500k SNPs with genome-wide significance (green) and the functional intronic SNP rs3846662 (red). Coding exons are shown as blue boxes. Exon13, which is subject to alternative splicing, is double framed.

To analyze whether rs3846662 was associated with HMGCR splicing efficiency we obtained lymphoblastoid cell lines (LCLs) from white CEU individuals of the HapMap collection, who were either homozygous for the major (rs3846662/AA) or minor allele (rs3846662/GG). Likewise, these individuals were also homozygous for the 3 proxy SNPs with genome-wide significance in the association meta-analysis (rs7703051, rs12654264, and rs3846663), where these SNPs denote a haplotype. LCLs (n=9 per genotype) were seeded in medium supplemented with 10% fetal calf serum, and mRNA was harvested at various time points (0 hour, 10 hours, 24 hours, and 48 hours) for expression analysis. We first analyzed whether total HMGCR mRNA expression differed between both groups. As shown in Figure 2A, we did not detect significant differences in total HMGCR mRNA levels between the 2 groups at any time point, indicating that total HMGCR mRNA expression is not influenced by allele status at these SNPs.


Figure 2
View larger version (17K):
[in this window]
[in a new window]

 
Figure 2. HMGCR mRNA expression in lymphoblastoid cells of major and minor allele homozygotes at rs3846662 (n=9 per genotype). A, Total HMGCR mRNA expression normalized to β-actin over 48 hours. B, Alternatively spliced HMGCR mRNA ({Delta}ex13 mRNA) as percentage of total HMGCR mRNA expression. *P<0.01 and {dagger}P<0.05 between major and minor allele samples.

We went on to determine the amounts of full-length and alternatively spliced HMGCR ({Delta}exon13) mRNA separately. {Delta}exon13 HMGCR mRNA was detectable in all samples and showed significant variation along the time course (Figure 2B). We observed a distinct decrease in the percentage of {Delta}exon13 HMGCR mRNA in both genotype groups over the first 10 hours. However, the decrease in percentage of {Delta}exon13 HMGCR mRNA was significantly less pronounced in LCLs from homozygotes for the rs3846662 major allele. Hence, the percentage of {Delta}exon13 HMGCR mRNA per total HMGCR mRNA was significantly higher in homozygotes for the major allele as compared to homozygotes for the minor allele at 10 hours, 24 hours, and 48 hours (10 hours: 23.0±7.6 versus 10.4±3.4, 24 hours: 23.1±7.9 versus 11.6±3.8, 48 hours: 27.1±9.4 versus 15.7±9.1, % {Delta}exon13/total HMGCR mRNA, major versus minor allele, Figure 2B). In contrast to the other time points, the difference in {Delta}exon13 mRNA expression at 0 hours did not reach significance, possibly attributable to regulatory mechanisms that first need to be initiated during this stage of acute adaption.

Similar results were obtained when we studied expression patterns in LCLs which were incubated in medium supplemented with 10% lipoprotein deficient serum (supplemental Figure III).

Expression Pattern of Full-Length and {Delta}Exon13 HMGCR mRNA in Human Tissues
We next quantified the expression levels of both HMGCR transcripts in vivo, using cDNA samples from various human tissues. Both variants were expressed in cDNA pools from human liver, brain, spleen, lung, placenta, kidney, heart, ovary, peripheral blood leukocytes, skeletal muscle, and small intestine, but their relative amounts differed significantly. The percentage of {Delta}exon13 HMGCR mRNA per total HMGCR mRNA varied between 7 and 18%, with the exception of peripheral blood leukocytes. In peripheral blood leukocytes, {Delta}exon13 HMGCR mRNA accounted for 79% of total HMGCR mRNA transcripts (Figure 3).


Figure 3
View larger version (20K):
[in this window]
[in a new window]

 
Figure 3. Absolute amounts of full-length and {Delta}exon13 HMGCR mRNA expression in cDNA pools of various human tissues (Clontech). HMGCR transcripts were normalized to β-actin as housekeeping gene.

SNP rs3846662 Modulates Alternative Splicing of Exon13 in HMGCR Minigene Constructs
Our studies in LCLs demonstrated that the amount of {Delta}exon13 HMGCR mRNA was associated with allele status at SNP rs3846662. Therefore, we aimed to specifically evaluate the functionality of SNP rs3846662 in splicing efficiency. We created exon-trapping vectors containing the genomic DNA sequence of HMGCR from intron12 to intron14 of rs3846662/AA (major allele) and rs3846662/GG (minor allele) individuals, respectively, and transfected them into HEK293 cells.

In accordance with our previous results in human LCLs, we found significantly lower levels for {Delta}exon13 HMGCR mRNA in cells transfected with the minor allele minigene (rs3846662/G) as compared to cells transfected with the major allele minigene (rs3846662/A) (Figure 4). The difference in exon13 splicing efficiency between the two minigenes was 20.8% (42.9±3.9 versus 63.7±1.0%{Delta}exon13 HMGCR mRNA/total HMGCR mRNA, P=0.02). This difference in splicing efficiency was abolished when we transfected a construct in which we had used site directed mutagenesis to convert rs3846662/A to the minor G allele (Figure 4), further corroborating that allelic variants at rs3846662 directly modulate the efficiency of HMGCR exon13 splicing.


Figure 4
View larger version (20K):
[in this window]
[in a new window]

 
Figure 4. HMGCR minigene splicing in transfected HEK cells. Alternatively spliced RNA ({Delta}ex13) from minigene constructs is expressed as percentage of total HMGCR minigene RNA. rs3846662/G: minor allele minigene, rs3846662/A: major allele minigene, mut A->G: major allele minigene after site directed mutagenesis, *P<0.05 rs3846662/G vs rs3846662/A; {dagger}P<0.05 rs3846662/A vs mut A->G.

Full-Length and {Delta}Exon13 HMGCR Expression in Stable Transfected UT-2 Cells
Alternative splicing of HMGCR mRNA leads to an in-frame deletion of 53 amino acids in the catalytic domain of the protein. To investigate the effect of this deletion on enzyme activity, we stably expressed human full-length (UT-2+FL) and {Delta}exon13 (UT-2+ex13) HMGCR variants at comparable levels in UT-2 cells (Figure 5A), a CHO cell-line that lacks HMGCR activity and requires exogenous mevalonate for growth.16 UT-2+FL cells displayed 51% HMGCR enzyme activity of wild-type CHO cells, whereas UT-2+ex13 cells lacked enzyme activity and were indistinguishable from control UT-2 cells (Figure 5B). Further, UT-2+FL cells grew in the absence of mevalonate, whereas UT-2+ex13 and parental UT-2 cells died without mevalonate supplementation (Figure 5C), suggesting that the {Delta}exon13 HMGCR variant is unable to restore enzyme activity in these cells.


Figure 5
View larger version (82K):
[in this window]
[in a new window]

 
Figure 5. Characterization of UT-2 cells with stable expression of human full-length (UT-2+FL) and {Delta}exon13 (UT-2+ex13) HMGCR: (A) HMGCR mRNA expression; (B) HMGCR activity normalized to wild-type CHO cells; (C) Growth characteristics of UT-2, UT-2+ex13, and UT-2+FL cells in the presence and absence of 0.2 mmol/L mevalonate.


*    Discussion
up arrowTop
up arrowAbstract
up arrowIntroduction
up arrowMethods
up arrowResults
*Discussion
down arrowReferences
 
We identified variants in the HMGCR gene that were among our top hits for LDL-C in a GWAS in a population from the Island of Kosrae. We then conducted in vitro studies to follow-up the association signals and to identify a functional variant at the HMGCR locus. We present evidence that a common intronic SNP (rs384662) that is in linkage disequilibrium with the variants typed in the genome scan modulates alternative splicing of HMGCR mRNA. The resulting splice variant could not restore enzyme activity when expressed in HMGCR deficient UT-2 cells.

HMG-CoA reductase is a key enzyme in cholesterol homeostasis and catalyzes the rate-limiting step in cholesterol biosynthesis.21 In contrast to other well known determinants of cholesterol homeostasis, eg, LDL-receptor or Apolipoprotein E, associations between variants in HMGCR and LDL-C have only recently emerged in the context of GWAS. As in the Kosrae study, the initial results of the DGI GWAS in 2758 whites supported associations between SNPs in HMGCR and LDL-C but did not meet the statistical threshold of genome-wide significance by themselves (best associated SNP rs12654264: P=4.09x10–4).13 In this study, genome-wide significance was clearly established for HMGCR SNP rs12654264 after validation in 3 additional white cohorts, resulting in a combined probability value of 1x10–20 in a total of {approx}18 000 subjects.13 However, in a separate study, the associations between SNPs in HMGCR and LDL-C that were observed in the DGI study were not strengthened by a meta analysis approach, consisting of the DGI and 2 other GWAS in whites (best associated SNP rs3846663: P=2.79x10–4).12 This discrepancy might be attributable to some source of heterogeneity, eg, differences in sample ascertainment or the impact of nonadditive interactions with other genetic variants or unaccounted environmental exposures.22,23 Combining the association results from the Kosrae and DGI studies revealed 3 variants in LD (r2>0.81) with genome-wide significance at the HMGCR locus, including the 2 SNPs mentioned above (rs12654264, rs3846663) and SNP rs7703051. Our data obtained in the Kosrae isolate thereby adds important evidence about the generalizibility of genetic associations at the HMGCR locus, demonstrating that these associations also extend to other ancestries. Interestingly, 2 pharmacogenetic studies investigating whether genetic variants in HMGCR influence response to statin therapy demonstrated that common SNP haplotypes in HMGCR contribute to variation in statin response.24,25 These haplotypes included the SNPs that were associated with plasma LDL-C in the Kosrae and DGI studies, and it is possible that the same underlying mechanisms contribute to variation in LDL-C levels and variation in statin response.

A major aspect of our study was to follow-up the findings from the GWAS and to identify the putative functional variant at the HMGCR locus. To address this question we used human lymphoblastoid cells from the HapMap CEU collection, which have previously been established as a suitable model to study the regulation of cholesterol biosynthesis in normal subjects and subjects with genetic abnormalities in lipid metabolism.26 Our efforts were facilitated by a near complete inventory (99%) of all common (>5% minor allele frequency) regional sequence variations, resulting from resequencing of the complete HMGCR locus in 23 whites.19 Because the only known common coding SNP in HMGCR (rs5908, I638V) is not in LD with any of the genotyped SNPs, we consider it to be unlikely that this variant is responsible for the association signal. Likewise, because we did not detect significant differences in total HMGCR mRNA expression, we consider it to be unlikely that the causal SNP is located in a regulatory element affecting HMGCR transcription. On the other hand, we provide mutually supportive evidence that a common intronic variant (rs3846662) in LD with the genotyped variants is functional and alters the efficiency of HMGCR exon13 alternative splicing: We could demonstrate that (1) expression levels of alternatively spliced {Delta}exon13 HMGCR mRNA were significantly lower in lymphoblastoid cells from homozygotes for the rs3846662 minor allele, and (2) allele status at rs3846662 directly modulated alternative splicing of HMGCR mRNA in minigene constructs. Further, alternative splicing of HMGCR appeared to be regulated and was present in vivo, as we could detect {Delta}exon13 HMGCR mRNA in all 11 human tissues that we studied.

HMGCR mRNA lacking exon13 was described in a survey of alternative pre-mRNA splicing by Johnson et al,20 however its function and the underlying mechanisms remain unknown. The regulation of gene splicing in mammalians involves both cis- and transfactors, which are composed of auxiliary element sequences in the pre mRNA, known as splicing enhancers and silencers27 and cellular splicing factors which include several protein families.28 The most likely explanation for the observed differences in HMGCR mRNA splicing between major and minor allele homozygotes at rs3846662 is that this SNP is located in a binding motif for a splice auxiliary protein and allele status changes the binding affinity of this protein.

Homozygosity for the major allele at rs3846662 increased the proportion of HMGCR mRNA lacking exon13. Skipping of exon 13 (159 bp) does not change the reading frame, and the resulting protein lacks 53 amino acids in the catalytic domain. When we stably expressed both HMGCR variants in CHO cells deficient of endogenous HMGCR activity, the {Delta}exon13 variant appeared to be nonfunctional and was not able to restore cell growth in the absence of mevalonate. At present, we can only speculate about the exact underlying mechanisms of this observation. Exon13 encodes parts of the catalytic domain, and it contains the highly conserved sequence element ENVIGX3I/LP which is thought to mediate dimerization of the enzyme’s monomers.29 Thus, deletion of exon13 could potentially impact the stability of the enzyme, because experiments in which monomeric soluble proteins were fused to the HMGCR membrane domains illustrated that the protein was degraded faster when it was smaller than tetrameric.30 In addition, exon13 contains the E559 residue which is located at the front of the active site and was proposed to directly participate in the reduction of HMG-CoA by serving as a proton donor to mevaldehyde.29 Therefore, alternative splicing of HMGCR appears to result in altered enzymatic activity and could also lead to more rapid degradation of the protein. A decrease in HMGCR activity would lead to lower cellular cholesterol synthesis and subsequently a counter-regulatory increase of cholesterol uptake from the plasma via the LDL-receptor pathway to maintain intracellular cholesterol homeostasis. In accordance with this hypothesis the allele at rs3846662 that was causing higher levels of {Delta}exon13 HMGCR mRNA in our in vitro studies was sharing a haplotype with the alleles that were associated with lower LDL-C in the genome-wide association studies. HMG-CoA reductase activity is subject to multivalent control on transcriptional and posttranscriptional levels, and alternative splicing may be an additional regulatory mechanism.

Modulation of alternatively spliced HMGCR mRNA levels could be of pharmacological interest with regard to response to statin therapy or as target for antisense-mediated exon skipping. Recently, an antisense oligonucleotide (AON)-mediated skipping approach related to lowering plasma cholesterol levels was applied by Khoo et al.31 In their study, AON-mediated exon27 skipping of the Apolipoprotein B transcript specifically lowered the amount of functional ApoB100 protein while maintaining ApoB48 levels.31

Therefore, identification of specific factors that regulate HMGCR alternative splicing and elucidating the underlying mechanism may lead to a better understanding of its impact on regulating cellular cholesterol homeostasis and plasma cholesterol levels.


*    Acknowledgments
 
Addendum

While this manuscript was under review another study reported that alternative splicing of HMG-CoA Reductase exon13 is associated with plasma LDL-C response to simvastatin (Medina MW, et al. Circulation. 2008;118:355–362), further supporting the functional significance of HMGCR alternative splicing.

Sources of Funding

Ralph Burkhardt is a fellow of the Deutsche Forschungsgemeinschaft (DFG Bu2263/1-1).

Disclosures

None.


*    Footnotes
 
Original received June 17, 2008; final version accepted August 1, 2008.


*    References
up arrowTop
up arrowAbstract
up arrowIntroduction
up arrowMethods
up arrowResults
up arrowDiscussion
*References
 
1. Kuulasmaa K, Tunstall-Pedoe H, Dobson A, Fortmann S, Sans S, Tolonen H, Evans A, Ferrario M, Tuomilehto J. Estimation of contribution of changes in classic risk factors to trends in coronary-event rates across the WHO MONICA Project populations. Lancet. 2000; 355: 675–687.[CrossRef][Medline] [Order article via Infotrieve]

2. Heller DA, de Faire U, Pedersen NL, Dahlen G, McClearn GE. Genetic and environmental influences on serum lipid levels in twins. N Engl J Med. 1993; 328: 1150–1156.[Abstract/Free Full Text]

3. Breslow JL. Genetics of lipoprotein abnormalities associated with coronary artery disease susceptibility. Annu Rev Genet. 2000; 34: 233–254.[CrossRef][Medline] [Order article via Infotrieve]

4. Rader DJ, Cohen J, Hobbs HH. Monogenic hypercholesterolemia: new insights in pathogenesis and treatment. J Clin Invest. 2003; 111: 1795–1803.[CrossRef][Medline] [Order article via Infotrieve]

5. Han Z, Heath SC, Shmulewitz D, Li W, Auerbach SB, Blundell ML, Lehner T, Ott J, Stoffel M, Friedman JM, Breslow JL. Candidate genes involved in cardiovascular risk factors by a family-based association study on the island of Kosrae, Federated States of Micronesia. Am J Med Genet. 2002; 110: 234–242.[CrossRef][Medline] [Order article via Infotrieve]

6. Shmulewitz D, Heath SC, Blundell ML, Han Z, Sharma R, Salit J, Auerbach SB, Signorini S, Breslow JL, Stoffel M, Friedman JM. Linkage analysis of quantitative traits for obesity, diabetes, hypertension, and dyslipidemia on the island of Kosrae, Federated States of Micronesia. Proc Natl Acad Sci U S A. 2006; 103: 3502–3509.[Abstract/Free Full Text]

7. Shmulewitz D, Auerbach SB, Lehner T, Blundell ML, Winick JD, Youngman LD, Skilling V, Heath SC, Ott J, Stoffel M, Breslow JL, Friedman JM. Epidemiology and factor analysis of obesity, type II diabetes, hypertension, and dyslipidemia (syndrome X) on the Island of Kosrae, Federated States of Micronesia. Hum Hered. 2001; 51: 8–19.[Medline] [Order article via Infotrieve]

8. Kruglyak L. The road to genome-wide association studies. Nat Rev Genet. 2008; 9: 314–318.[CrossRef][Medline] [Order article via Infotrieve]

9. Breslow JL, McPherson J, Nussbaum AL, Williams HW, Lofquist-Kahl F, Karathanasis SK, Zannis VI. Identification and DNA sequence of a human apolipoprotein E cDNA clone. J Biol Chem. 1982; 257: 14639–14641.[Abstract/Free Full Text]

10. Wallace C, Newhouse SJ, Braund P, Zhang F, Tobin M, Falchi M, Ahmadi K, Dobson RJ, Marcano AC, Hajat C, Burton P, Deloukas P, Brown M, Connell JM, Dominiczak A, Lathrop GM, Webster J, Farrall M, Spector T, Samani NJ, Caulfield MJ, Munroe PB. Genome-wide association study identifies genes for biomarkers of cardiovascular disease: serum urate and dyslipidemia. Am J Hum Genet. 2008; 82: 139–149.[CrossRef][Medline] [Order article via Infotrieve]

11. Sandhu MS, Waterworth DM, Debenham SL, Wheeler E, Papadakis K, Zhao JH, Song K, Yuan X, Johnson T, Ashford S, Inouye M, Luben R, Sims M, Hadley D, McArdle W, Barter P, Kesaniemi YA, Mahley RW, McPherson R, Grundy SM, Bingham SA, Khaw KT, Loos RJ, Waeber G, Barroso I, Strachan DP, Deloukas P, Vollenweider P, Wareham NJ, Mooser V. LDL-cholesterol concentrations: a genome-wide association study. Lancet. 2008; 371: 483–491.[CrossRef][Medline] [Order article via Infotrieve]

12. Willer CJ, Sanna S, Jackson AU, Scuteri A, Bonnycastle LL, Clarke R, Heath SC, Timpson NJ, Najjar SS, Stringham HM, Strait J, Duren WL, Maschio A, Busonero F, Mulas A, Albai G, Swift AJ, Morken MA, Narisu N, Bennett D, Parish S, Shen H, Galan P, Meneton P, Hercberg S, Zelenika D, Chen WM, Li Y, Scott LJ, Scheet PA, Sundvall J, Watanabe RM, Nagaraja R, Ebrahim S, Lawlor DA, Ben-Shlomo Y, Davey-Smith G, Shuldiner AR, Collins R, Bergman RN, Uda M, Tuomilehto J, Cao A, Collins FS, Lakatta E, Lathrop GM, Boehnke M, Schlessinger D, Mohlke KL, Abecasis GR. Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nat Genet. 2008; 40: 161–169.[CrossRef][Medline] [Order article via Infotrieve]

13. Kathiresan S, Melander O, Guiducci C, Surti A, Burtt NP, Rieder MJ, Cooper GM, Roos C, Voight BF, Havulinna AS, Wahlstrand B, Hedner T, Corella D, Tai ES, Ordovas JM, Berglund G, Vartiainen E, Jousilahti P, Hedblad B, Taskinen MR, Newton-Cheh C, Salomaa V, Peltonen L, Groop L, Altshuler DM, Orho-Melander M. Six new loci associated with blood low-density lipoprotein cholesterol, high-density lipoprotein cholesterol or triglycerides in humans. Nat Genet. 2008; 40: 189–197.[CrossRef][Medline] [Order article via Infotrieve]

14. Saxena R, Voight BF, Lyssenko V, Burtt NP, de Bakker PI, Chen H, Roix JJ, Kathiresan S, Hirschhorn JN, Daly MJ, Hughes TE, Groop L, Altshuler D, Almgren P, Florez JC, Meyer J, Ardlie K, Bengtsson Bostrom K, Isomaa B, Lettre G, Lindblad U, Lyon HN, Melander O, Newton-Cheh C, Nilsson P, Orho-Melander M, Rastam L, Speliotes EK, Taskinen MR, Tuomi T, Guiducci C, Berglund A, Carlson J, Gianniny L, Hackett R, Hall L, Holmkvist J, Laurila E, Sjogren M, Sterner M, Surti A, Svensson M, Tewhey R, Blumenstiel B, Parkin M, Defelice M, Barry R, Brodeur W, Camarata J, Chia N, Fava M, Gibbons J, Handsaker B, Healy C, Nguyen K, Gates C, Sougnez C, Gage D, Nizzari M, Gabriel SB, Chirn GW, Ma Q, Parikh H, Richardson D, Ricke D, Purcell S. Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science. 2007; 316: 1331–1336.[Abstract/Free Full Text]

15. Fisher RA. Statistical Methods for Research Workers. London: Oliver and Boyd; 1932.

16. Mosley ST, Brown MS, Anderson RG, Goldstein JL. Mutant clone of Chinese hamster ovary cells lacking 3-hydroxy-3 -methylglutaryl coenzyme A reductase. J Biol Chem. 1983; 258: 13875–13881.[Abstract/Free Full Text]

17. Walli AK, Seidel D. Role of lipoprotein-X in the pathogenesis of cholestatic hypercholesterolemia. Uptake of lipoprotein-X and its effect on 3-hydroxy-3-methylglutaryl coenzyme A reductase and chylomicron remnant removal in human fibroblasts, lymphocytes, and in the rat. J Clin Invest. 1984; 74: 867–879.[CrossRef][Medline] [Order article via Infotrieve]

18. Edwards PA, Lemongello D, Fogelman AM. Improved methods for the solubilization and assay of hepatic 3-hydroxy-3-methylglutaryl coenzyme A reductase. J Lipid Res. 1979; 20: 40–46.[Abstract]

19. Crawford DC, Carlson CS, Rieder MJ, Carrington DP, Yi Q, Smith JD, Eberle MA, Kruglyak L, Nickerson DA. Haplotype diversity across 100 candidate genes for inflammation, lipid metabolism, and blood pressure regulation in two populations. Am J Hum Genet. 2004; 74: 610–622.[CrossRef][Medline] [Order article via Infotrieve]

20. Johnson JM, Castle J, Garrett-Engele P, Kan Z, Loerch PM, Armour CD, Santos R, Schadt EE, Stoughton R, Shoemaker DD. Genome-wide survey of human alternative pre-mRNA splicing with exon junction microarrays. Science. 2003; 302: 2141–2144.[Abstract/Free Full Text]

21. Goldstein JL, Brown MS. Regulation of the mevalonate pathway. Nature. 1990; 343: 425–430.[CrossRef][Medline] [Order article via Infotrieve]

22. McCarthy MI, Abecasis GR, Cardon LR, Goldstein DB, Little J, Ioannidis JP, Hirschhorn JN. Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nat Rev Genet. 2008; 9: 356–369.[CrossRef][Medline] [Order article via Infotrieve]

23. Ioannidis JP, Patsopoulos NA, Evangelou E. Heterogeneity in meta-analyses of genome-wide association investigations. PLoS ONE. 2007; 2: e841.[CrossRef]

24. Chasman DI, Posada D, Subrahmanyan L, Cook NR, Stanton VP Jr, Ridker PM. Pharmacogenetic study of statin therapy and cholesterol reduction. JAMA. 2004; 291: 2821–2827.[Abstract/Free Full Text]

25. Krauss RM, Mangravite LM, Smith JD, Medina MW, Wang D, Guo X, Rieder MJ, Simon JA, Hulley SB, Waters D, Saad M, Williams PT, Taylor KD, Yang H, Nickerson DA, Rotter JI. 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. Circulation. 2008; 117: 1537–1544.[Abstract/Free Full Text]

26. Kayden HJ, Hatam L, Beratis NG. Regulation of 3-hydroxy-3-methylglutaryl coenzyme A reductase activity and the esterification of cholesterol in human long term lymphoid cell lines. Biochemistry. 1976; 15: 521–528.[CrossRef][Medline] [Order article via Infotrieve]

27. Matlin AJ, Clark F, Smith CW. Understanding alternative splicing: towards a cellular code. Nat Rev Mol Cell Biol. 2005; 6: 386–398.[CrossRef][Medline] [Order article via Infotrieve]

28. Caceres JF, Kornblihtt AR. Alternative splicing: multiple control mechanisms and involvement in human disease. Trends Genet. 2002; 18: 186–193.[CrossRef][Medline] [Order article via Infotrieve]

29. Istvan ES, Palnitkar M, Buchanan SK, Deisenhofer J. Crystal structure of the catalytic portion of human HMG-CoA reductase: insights into regulation of activity and catalysis. EMBO J. 2000; 19: 819–830.[CrossRef][Medline] [Order article via Infotrieve]

30. Cheng HH, Xu L, Kumagai H, Simoni RD. Oligomerization state influences the degradation rate of 3-hydroxy-3-methylglutaryl-CoA reductase. J Biol Chem. 1999; 274: 17171–17178.[Abstract/Free Full Text]

31. Khoo B, Roca X, Chew SL, Krainer AR. Antisense oligonucleotide-induced alternative splicing of the APOB mRNA generates a novel isoform of APOB. BMC Mol Biol. 2007; 8: 3.[CrossRef][Medline] [Order article via Infotrieve]




This article has been cited by other articles:


Home page
Proc. Natl. Acad. Sci. USAHome page
E. E. Kenny, A. Gusev, K. Riegel, D. Lutjohann, J. K. Lowe, J. Salit, J. B. Maller, M. Stoffel, M. J. Daly, D. M. Altshuler, et al.
Systematic haplotype analysis resolves a complex plasma plant sterol locus on the Micronesian Island of Kosrae
PNAS, August 18, 2009; 106(33): 13886 - 13891.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Data Supplement
Right arrow All Versions of this Article:
28/11/2078    most recent
ATVBAHA.108.172288v1
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowRequest Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Burkhardt, R.
Right arrow Articles by Breslow, J. L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Burkhardt, R.
Right arrow Articles by Breslow, J. L.
Related Collections
Right arrow Functional genomics
Right arrow Gene expression
Right arrow Gene regulation
Right arrow Genomics
Right arrow Genetics of cardiovascular disease
Right arrow Lipid and lipoprotein metabolism