An investigation of the effects of lipid-lowering medications: genome-wide linkage analysis of lipids in the HyperGEN study

Abstract

Background

Use of anti-hyperlipidemic medications compromises genetic analysis because of altered lipid profiles. We propose an empirical method to adjust lipid levels for medication effects so that the adjusted lipid values substitute the unmedicated lipid values in the genetic analysis.

Results

Published clinical trials were reviewed for HMG-CoA reductase inhibitors and fibric acid derivatives as mono-drug therapy. HMG-CoA reductase inhibitors showed similar effects in African Americans (AA) and non-African Americans (non-AA) for lowering total cholesterol (TC, -50.7 mg/dl), LDL cholesterol (LDL-C, -48.1 mg/dl), and triglycerides (TG, -19.7 mg/dl). Their effect on increasing HDL cholesterol (HDL-C) in AA (+0.4 mg/dl) was lower than in Non-AA (+2.3 mg/dl). The effects of fibric acid derivatives were estimated as -46.1 mg/dl for TC, -40.1 mg/dl for LDL-C, and +5.9 mg/dl for HDL-C in non-AA. The corresponding effects in AA were less extreme (-20.1 mg/dl, -11.4 mg/dl, and +3.1 mg/dl). Similar effect for TG (59.0 mg/dl) was shown in AA and non-AA. The above estimated effects were applied to a multipoint variance components linkage analysis on the lipid levels in 2,403 Whites and 2,214 AA in the HyperGEN study. The familial effects did vary depending on whether the lipids were adjusted for medication use. For example, the heritabilities increased after medication adjustment for TC and LDL-C, but did not change significantly for HDL-C and TG.

Conclusion

Ethnicity-specific medication adjustments using our empirical method can be employed in epidemiological and genetic analysis of lipids.

Background

There is extensive epidemiological evidence showing that a dyslipidemic profile, characterized by elevated plasma levels of total cholesterol (TC), triglycerides (TG) and LDL cholesterol (LDL-C), and reduced levels of HDL cholesterol (HDL-C), is strongly associated with an increased risk of atherosclerosis and coronary heart disease [1]. Genetic contributions to plasma lipid levels also have been documented, with the magnitude as assessed by heritability being quite variable (26% to 83%) [2–4]. Similarly, genome scans from different populations have reported a variety of quantitative trait loci (QTLs) related to lipid profiles [5–13].

However, it is known that use of anti-hyperlipidemic medications by participants in such family studies distorts their blood lipid levels. Consequently, investigators may exclude subjects on medications during recruitment, or remove them at the time of the genetic analyses. Not only do such exclusions results in a smaller sample size, they also may discard the participants who are most likely to carry causative variants of dyslipidemia. To increase the power of the analysis by including these potentially informative subjects, some studies were designed to interrupt lipid lowering medication for specific time intervals before the blood samples were taken. Other methods that have been implemented are to adjust the lipid levels for medication by using medications as covariates in a linear regression model. This method tends to minimize the differences in lipid levels between medicated and unmedicated groups, but not necessarily adjust lipid levels toward their original values. To our knowledge, only one study has performed a genome-wide scan for LDL-C after increasing the measured levels by 25% for medicated subjects [13]. The lod score for their highest peak increased from 3.27, observed when analyzing 1977 white subjects not on lipid lowering medications, to 3.72 after including 350 medicated subjects.

The aim of the present study is to recover or impute the original unmedicated lipid levels using an empirical approach. The measured (biased) lipid values of individuals on lipid-lowering drugs are "adjusted" by increasing (or decreasing) the measured value by the same amount that clinical trials studies report decreases (or increases) for a given class of drugs using ethnic-specific data. That is, the approach statistically infers or projects what the original lipid value would have been prior to medication treatment. The impact of the medication adjustments will be evaluated by computing heritabilities and by conducting genome-wide linkage analyses in three samples: by excluding medicated subjects; by including medicated subjects after adjusting lipid values for medication use; and by including medicated subjects without any medication adjustment. We hypothesize that imputing the original lipid measurement by adjusting for medication effects can provide more precise estimates of the medication-free values of the phenotypes, with an expected increase in power to detect genes involved in the lipid metabolism.

Methods

Subject description

The Hypertension Genetic Epidemiology Network (HyperGEN) study is part of the National Heart, Lung and Blood Institute (NHLBI) Family Blood Pressure Program. The design and methods have been described in detail elsewhere [14]. In the present report 2,403 White and 2,214 African American participants represent 1,924 pairs of White hypertensive siblings in 1,242 nuclear families and 1,889 pairs of African American hypertensive siblings in 1,801 nuclear families. Among all participants with measured serum lipid levels, 402 subjects reported that they were taking medications to lower blood lipid levels. This included 385 subjects taking one medication, 16 subjects on combination drug therapy of two medications, and 1 subject reporting the use of three drugs. IRB approvals were obtained for all participating centers, and all participants gave their written informed consent.

Measurement of Lipids and Other Phenotypes

All blood assays on samples were performed at the Central Biochemistry Laboratory at the University of Minnesota. TG levels were colorimetrically determined by using a peroxidase-coupled method on a centrifugal analyzer [15]. TC was measured by using a commercial cholesterol oxidase method [16]. LDL-C was estimated by the Friedewald equation [17], or by ultracentrifugation for subjects with TG above 400 mg/dl [18]. HDL-C was quantified by COBAS FARA (Roche Diagnostic Systems) after precipitation of the other lipoprotein fractions by dextran sulfate [19].

Anthropometric measurements were collected with subjects wearing lightweight suits. Body mass index (BMI) was subsequently calculated as weight in kilograms divided by height in meters squared. Waist circumference was measured to the nearest centimeter at the level of the umbilicus, and hip circumference was measured at the level of the maximal circumference of the gluteus. All other measurements were collected using a common, standardized protocol administered by trained personnel.

Genotyping

The NHLBI Mammalian Genotyping Service (Marshfield, Wisconsin) genotyped 391 markers at approximate 10-cM intervals throughout the genome, with an average marker heterozygosity of 76%. Detailed information on the genotyping methods, instruments, and software used can be found at the Marshfield Laboratory World Wide Web site [20]. Quality control of the markers and the relationships in the pedigrees were carried out at the HyperGEN Data Coordinating Center at Washington University in Saint Louis, by using ASPEX [21], GRR [22], MAPMAKER/SIBS [23], and PEDCHECK [24].

Adjustment for the effects of lipid-lowering medications

The majority of the participants in the HyperGEN study, that were using anti-hyperlipidemic medications, were prescribed two groups of medications: HMG-CoA reductase inhibitors (statins) and/or fibric acid derivatives (Table 1). Therefore, we only focused on lipid-lowering effects of these two groups of medications. Clinical trials on the long-term effects of mono-drug therapy using these anti-hyperlipidemic medications published in English during 1989–2003 were searched online through PubMed. The search terms used included lipids, efficacy, mono-drug therapy, combination drug therapy, as well as the names of medication classes. We excluded clinical trials with duration of less than 8 weeks and those without placebo control. Clinical trials for medication combinations were also excluded unless they reported the data for mono-drug therapy as well. Tables A1 and A2 in Appendix A [see Additional File 1] summarize 32 clinical trials that have investigated the effects of HMG-CoA inhibitors and fibric acid derivatives as mono-drug therapy with commonly used dosages. A total of 19,384 participants, aged 18 to 80 years, with high blood lipid levels were pooled. The clinical trials were double-blinded, randomized with a course which varied from 8 to 280 weeks. Gemfibrozil was the only fibric acid derivative used by the subjects in the HyperGEN study; this medication was well represented in non-African Americans of the HyperGEN, but only two African Americans used the drug (Table 1). Therefore we reviewed only the clinical trials for the effects of this medication in non-African Americans (see Additional File 1 Table A2).

Table 1 Mechanisms of action for the major groups of antihyperlipidemic medications and their frequency of use by the subjects in HyperGEN

The weighted (by sample size) average absolute change for each lipid variables, which took into account placebo effect, was calculated by medication groups based on ethnicity (African Americans vs. non-African Americans, and vs. all ethnicities) (Table 2). Each phenotype of lipids responded to different classes of drugs differently, and the response also varied considerably by ethnicity. Non-African Americans included mainly Whites and Hispanics. The clinical trials for all ethnicities include more trials, because some studies had samples that could not be categorized into either African Americans or non-African Americans. HMG-CoA inhibitors showed similar reductions of TC, LDL-C and TG in African Americans and non-African Americans. Thus, the following effect estimates, -50.7 mg/dl, -48.1 mg/dl, and -19.7 mg/dl, were used for all ethnicities in performing the adjustments, respectively. However, non-African Americans appeared to be more responsive than African Americans to HMG-CoA inhibitors effecting HDL-C. The ethnicity specific average effects (+2.3 mg/dl for non-African Americans and +0.4 mg/dl for African Americans) were used to obtain the adjusted HDL-C. Gemfibrozil appeared to be more effective as an anti-hyperlipidemic agent for TC, LDL-C, and HDL-C in non-African Americans than in populations with all pooled ethnicities. For non-African Americans subjects receiving Gemfibrozil, the measured levels of TC, LDL-C, and HDL-C were adjusted by the following respective estimates: -46.1 mg/dl, -40.1 mg/dl, and +5.9 mg/dl. Instead, the measured lipid levels of two African Americans who took Gemfibrozil were adjusted by the effect for all ethnicities: -20.1 mg/dl for TC, -11.4 mg/dl for LDL-C, and + 3.1 mg/dl for HDL-C. Both African Americans and non-African Americans were adjusted by -59.0 mg/dl for TG.

Table 2 Average effects estimated from published clinical trails by ethnicity for HMG-CoA reductase inhibitors and fibric acid derivatives

For the participants receiving mono-drug therapy of anti-hyperlipidemics, the levels of TC, HDL-C, LDL-C, and TG were adjusted with the effect estimates, where the sign (+/-) indicates the increasing or lowering effect of a medication to a specific lipid component. When an individual took a combination therapy of two drugs, the total effect was estimated by selecting the medication with higher effect on the lipids and conservatively adding up half effect of the other medication. Unadjusted measurements were used for 12 African Americans and 2 Whites who took Bile Sequestrants or Nicotinic Acid Derivatives. The third medication effect, which was presented for only one participant, was not considered in our study.

Statistical Methods and Linkage Analysis

TG with and without medication adjustments were transformed using natural logarithm to an approximate normal distribution. Each phenotype was then adjusted for its mean and variance to remove the effects of age as well as the effects of BMI, current smoking and drinking status, the use of estrogen replacement therapy, diabetes status, and the field center. Stepwise multiple regressions were performed for the adjustment within each ethnicity and sex group, retaining only the terms significant at the 5% level. The adjusted variables were finally standardized to a mean of 0 and a standard deviation (SD) of 1. Four to ten outliers greater than 4 SD-s from the mean and at least 1 SD from the nearest data point were excluded within each ethnic group for the lipid phenotypes.

Heritability (h²), defined as the percent of variance due to additive familial effects, was estimated through the maximum-likelihood methods available in the software SEGPATH under the most parsimonious model using the following equation [25]:

$$h^2=\frac{(r_{parent-offspring}+r_{sibling})(1+r_{spouse})}{1+r_{spouse}+2r_{spouse}{r}_{parent-offspring}}$$

where r_{parent-offspring}, r_sibling, and r_spouse denote the average parent-offspring, the average sibling and the spouse correlations, respectively. Two heritability estimates were regarded as significantly different from each other if their difference was greater than the sum of their standard errors.

Genome-wide linkage analyses, both with and without adjustments for medication effects, were performed for each lipid phenotype by using variance component linkage analysis as implemented in SEGPATH [26]. Linkage analysis was carried out separately within ethnicity by using ethnicity-specific marker allele frequencies derived from the corresponding random sub-samples. Multipoint IBD estimates were used in the linkage analysis.

Results

Table 3 presents the selected characteristics of participants in the HyperGEN study by ethnic group and gender, and the effects of adjustment for medication for each lipid parameter. The average age was approximately 53 years in Whites and 46 years in African Americans. African American females had higher levels of BMI than the other groups. Alcohol consumption, smoking status, diabetes, estrogen use, and field center also showed significant ethnic differences. When lipid levels in subjects on lipid-lowering medications were adjusted for medication effects, there were observable elevations in the mean of TC and LDL-C, but very slight changes in the mean of HDL-C and in the log transformed TG.

Table 3 Characteristics of subjects by ethnicity and sex in the HyperGEN study

The effects of medication adjustments on the heritabilities of lipids are summarized in Table 4. Three groups of data were analyzed: 1) Excluding medicated participants; 2) All participants without any adjustment for medication effects; 3) All participants with adjustment for medication effects. Compared to the estimates without medication adjustment, the heritabilities in Whites improved significantly by 20% and 25% for TC and LDL-C, respectively, when medication adjustments were applied. In African Americans this significant improvement was 36% for TC and 27% for LDL-C. By contrast, the heritability estimates after medication adjustments for HDL-C and TG stayed similar in both races. The heritability estimates for the adjusted lipid levels were similar to those when excluding medicated participants, suggesting that the medication adjustment may recover the phenotypic correlations within pedigrees previously obscured by medication.

Table 4 The effect of medication correction on heritability of lipids: the HyperGEN study

A lod score of 1.75 (P < 0.0023) was set a priori as the threshold for promising linkages, which previously was shown to reduce the risk of type II errors to one false positive per genome scan (based on the assumption of discrete marker density with approximately 400 markers) [27]. Accordingly, all promising lod scores from the multipoint variance components linkage analysis for three samples, are presented in Table 5. Genome scan results of the four lipid traits for three groups of data are depicted in Figures A1-4 [see Additional File 1 Appendix A]. The strongest linkage signal (lod score = 3.08) was obtained for TG in African Americans excluding medicated subjects. This was localized at 131 cM from p-telomere (pter) on chromosome 11 (near marker AFM157XH6). The lod score dropped to 2.73 when medicated subjects were included without adjustment, and dropped to 2.71 when data were adjusted for the use of lipid-lowering medications. In contrast, the peak for TC on chromosome 6 at the region 53.81 cM p-ter (near marker GGAA15B08) increased in African Americans from 2.18 to 2.66 after the medication adjustment.

Table 5 Genome scan results (lod score ≥ 1.75) for lipids with and without medication correction by ethnicity: the HyperGEN study

Discussion

The purpose of this study was to propose and explore a method to adjust blood lipid measurements for the effects of anti-hyperlipidemic medications. To our knowledge, our empirical approach to adjust medication effects on lipids has not been applied before in the genetic epidemiology studies, although examples of other complex traits exist. Several approaches have been attempted to adjust blood pressure for the use of antihypertensive medications. We reported a similar empirical method to adjust blood pressure based on a summary of over 160 clinical trials on effects of anti-hypertensive medication [28]. A nonparametric algorithm by ranking (the adjusted) blood pressure was implemented by Levy et al. in the Framingham Heart Study [29]. Hunt et al. compared three methods, including the use of a threshold method assigning a fixed blood pressure (140/90 mm Hg) for all medicated participants, and a random assignment of systolic blood pressure between 140–160, and a random assignment of diastolic blood pressure between 90–100 mm Hg [30]. The investigators concluded that excluding medicated measurements or the use of unadjusted blood pressure both lowered lod score results as compared to adjusting blood pressure measurements. Cui et al. examined the adjustment effects by adding an arbitrary 10 mm Hg and 5 mm Hg to systolic and diastolic blood pressures respectively, according to a publication on clinical trials [31, 32]. The adjusted blood pressure resulted in an increase of both genetic and shared environmental variance components, while the unadjusted blood pressure for antihypertensive treatments obscured the relationship among phenotypes and genotypes.

The heritabilities of lipids, estimated in our study by excluding all the participants who took lipid-lowering medications, should be less biased heritabilities. However, by excluding subjects there is a smaller sample size and some valuable genetic information may be lost, especially alleles of the subjects affected more by dyslipidemia. Including medicated participants would introduce bias, because of a confounding of the medication effects with the lipid levels. That is, the measured lipid levels are artificial. This could cause a reduction in the genetic variance component as was seen for TC and LDL-C when we used the whole sample without performing medication adjustments. Heritabilities for TC and LDL-C were significantly improved after medication adjustments to the similar level as those obtained when excluding the medicated subjects. This indicates that the adjustment with our empirical method recovered relationships among phenotypes and genotypes obscured by the use of medications. Therefore, our method can provide the benefit of a greater power in genetic analysis by increasing the sample size as well as the information content.

The lod score for the TC peak in African Americans, centered at 53.81 cM on chromosome 6, increased from 2.03 to 2.66 when the medication adjustments were applied. This linkage signal is about 25 cM from a suggestive signal which was reported in a bivariate analysis of TG and LDL-C (LOD = 2.0) from 1,747 African American siblings in the Hypertension Genetic Epidemiology Network study [33]. Another study provided evidence of linkage (LOD = 2.14) at 102 cM on chromosome 6 for combined hyperlipidemia of TC and TG in 113 Caucasian families [5]. In contrast, the peak with a lod score of 3.08 dropped slightly to 2.71 when including those measures with medication adjustments for TG in African Americans on chromosome 11 at 131 cM from p-ter. This peak replicates a suggestive signal (LOD = 1.93) detected in French Canadians [6], and 30 cM from the region (LOD = 3.34) that have previously been identified in Turkish families [34]. We performed meta-analysis across African American and Caucasian families but there was no improvement in the LOD scores (results not shown).

Blood lipid levels are the result of a complex interaction among genetic and various environmental factors. The power of genome scans for detecting true signals is generally low with small-to-modest effect for complex traits like blood lipids. In addition, lipid-lowering treatment represents a pharmacologically induced alteration – typically of large magnitude – of the plasma lipid levels which further reduces the power to detect the true signals. One of the methods to overcome the confounding of medication effects with the measured lipid levels is to adjust the measured blood lipids by trying to recover the untreated levels as closely as possible. Although we noticed that the lod scores of blood lipids went both ways before and after medication adjustment in the present study, we expect lod scores for the medication adjusted analysis to be closest to their optimal (true) estimates. To evaluate our empirical approach of adjusting the treated lipid phenotypes, the heritability is a more appropriate parameter, since it reflects the effect size in linkage analysis. The higher heritabilities for TC and LDL-C when medication adjustments were made in the present study indicate an improvement in effect size of the linkage analysis with this empirical method. This is consistent with what is known about the HMG-CoA reductase inhibitors. First, as shown in Table 1 most of the medicated participants in the HyperGEN study took HMG-CoA reductase inhibitors and the effects is strongest on LDL-C and TC as shown in Table 3. Second, these inhibitors primarily target LDL-C (and by extension TC). For example, Table 2 shows larger mean and variance effects due to medication use for LDL-C and TC in the clinical trials literature.

Although our study has accounted for all participants and adjusted lipid measurements for aggregate medication effects, nevertheless there are limitations in the present study. The effects of anti-hyperlipidemics were estimated as group means rather than as individual responses to the medications, which are notoriously dose-dependent. Instead of producing precise untreated phenotype measurements, an adjustment with the average effects of medications may further serve as a shift by a constant and obscure the association between the phenotype and genotype. Also, we focused only on two groups of lipid-lowering medications, i.e., the most commonly used by subjects in the HyperGEN study. Under different circumstances an overall weighted average effect of other groups of anti-hyperlipidemic medications could be used. Another issue is that the present study did not ascertain the effect of the different patient groups which may influence antihyperlipidemic responses. Although dose-specific effects should be estimated when the information is available, the doses used in each of the source clinical trials were pooled to estimate the general efficacy for each group of medications because we do not have all the desired data in the HyperGEN study, such as dosage of anti-hyperlipidemics used by subjects. Similarly, we had to assume that the medicated subjects in the HyperGEN study had been prescribed, and were adhering, to the most commonly recommended doses. However, various strategies have been implemented to improve the compliance significantly during the clinical trials. The efficacy estimated from clinical trials is likely to overestimate the results observed in a general population of hyperlipidemic individuals under the treatment, such as in the HyperGEN study. Therefore, applying those estimates of drug effects would tend to over adjust the lipid values for studies using general populations.

Conclusion

In conclusion, this study provides a summary of a large number of clinical trials on the average effects of anti-hyperlipidemic medications on the actual lipid phenotypes, and an empirical method to adjust measured lipids for such medication effects. An illustrative example of how to apply this methodology is presented for the HyperGEN data. The medication adjustments using this method resulted in an improvement in the heritability estimates for TC, LDL-C, and TG. Genome scans with the imputed lipid measurements revealed some promising regions, such as a QTL with a peak lod score 2.66 at 53.81 cM from p-ter on chromosome 6 for TC and a region with lod score of 2.71 centered at 131 cM on chromosome 11 for TG in African Americans. The method described and the medication-use adjustments derived from a large body of clinical trials can be employed in other epidemiologic and genetic studies to obtain more accurate estimates of the lipid measurements as well to preserve the power of the study.