Email updates

Keep up to date with the latest news and content from BMC Genetics and BioMed Central.

Open Access Research article

Investigation of gene-by-sex interactions for lipid traits in diverse populations from the population architecture using genomics and epidemiology study

Kira C Taylor12*, Cara L Carty3, Logan Dumitrescu4, Petra Bůžková5, Shelley A Cole6, Lucia Hindorff7, Fred R Schumacher8, Lynne R Wilkens9, Ralph V Shohet10, P Miguel Quibrera1, Karen C Johnson11, Brian E Henderson128, Jeff Haessler3, Nora Franceschini1, Charles B Eaton13, David J Duggan14, Barbara Cochran15, Iona Cheng16, Chris S Carlson3, Kristin Brown-Gentry4, Garnet Anderson3, Jose Luis Ambite17, Christopher Haiman128, Loïc Le Marchand8, Charles Kooperberg3, Dana C Crawford184, Steven Buyske19, Kari E North120, Myriam Fornage2122 and for the PAGE Study

Author Affiliations

1 Department of Epidemiology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA

2 Department of Epidemiology and Population Health, School of Public Health and Information Sciences, University of Louisville, Louisville, KY, USA

3 Public Health Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA, USA

4 Center for Human Genetics Research, Vanderbilt University, Nashville, TN, USA

5 Department of Biostatistics, University of Washington, Seattle, WA, USA

6 Department of Genetics, Texas Biomedical Research Institute, San Antonio, TX, USA

7 Office of Population Genomics, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA

8 Department of Preventive Medicine, Keck School of Medicine, University of Southern California, Los Angeles, CA, USA

9 Epidemiology Program, University of Hawaii Cancer Center, University of Hawaii, Honolulu, HI, USA

10 Center of Cardiovascular Research, Department of Medicine, John A. Burns School of Medicine, University of Hawaii, Honolulu, HI, USA

11 Department of Preventive Medicine, University of Tennessee Health Science Center, Memphis, TN, USA

12 Zilkha Neurogenetic Institute, University of Southern California, Los Angeles, CA, USA

13 Department of Family Medicine and Community Health, Alpert Medical School of Brown University, Providence, RI, USA

14 Integrated Cancer Genomics Division, The Translational Genomics Research Institute, Phoenix, AZ, USA

15 Sponsored Programs, Baylor College of Medicine, Houston, TX, USA

16 Cancer Research Center, University of Hawaii, Honolulu, HI, USA

17 Information Sciences Institute, University of Southern California, Los Angeles, CA, USA

18 Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN, USA

19 Department of Statistics and Biostatistics, Rutgers University, Piscataway, NJ, USA

20 Carolina Center for Genome Sciences, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA

21 Division of Epidemiology, Human Genetics, and Environmental Sciences, School of Public Health, University of Texas Health Sciences Center at Houston, Houston, TX, USA

22 Institute of Molecular Medicine, University of Texas Health Sciences Center at Houston, Houston, TX, USA

For all author emails, please log on.

BMC Genetics 2013, 14:33  doi:10.1186/1471-2156-14-33


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2156/14/33


Received:14 August 2012
Accepted:17 April 2013
Published:1 May 2013

© 2013 Taylor et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

High-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglyceride (TG) levels are influenced by both genes and the environment. Genome-wide association studies (GWAS) have identified ~100 common genetic variants associated with HDL-C, LDL-C, and/or TG levels, mostly in populations of European descent, but little is known about the modifiers of these associations. Here, we investigated whether GWAS-identified SNPs for lipid traits exhibited heterogeneity by sex in the Population Architecture using Genomics and Epidemiology (PAGE) study.

Results

A sex-stratified meta-analysis was performed for 49 GWAS-identified SNPs for fasting HDL-C, LDL-C, and ln(TG) levels among adults self-identified as European American (25,013). Heterogeneity by sex was established when phet < 0.001. There was evidence for heterogeneity by sex for two SNPs for ln(TG) in the APOA1/C3/A4/A5/BUD13 gene cluster: rs28927680 (phet = 7.4x10-7) and rs3135506 (phet = 4.3x10-4), one SNP in PLTP for HDL levels (rs7679; phet = 9.9x10-4), and one in HMGCR for LDL levels (rs12654264; phet = 3.1x10-5). We replicated heterogeneity by sex in five of seventeen loci previously reported by genome-wide studies (binomial p = 0.0009). We also present results for other racial/ethnic groups in the supplementary materials, to provide a resource for future meta-analyses.

Conclusions

We provide further evidence for sex-specific effects of SNPs in the APOA1/C3/A4/A5/BUD13 gene cluster, PLTP, and HMGCR on fasting triglyceride levels in European Americans from the PAGE study. Our findings emphasize the need for considering context-specific effects when interpreting genetic associations emerging from GWAS, and also highlight the difficulties in replicating interaction effects across studies and across racial/ethnic groups.

Keywords:
Lipids; Genetics; Cardiovascular disease; Heterogeneity; Sex-specific effect; Association study

Background

The successes of genome-wide association studies (GWAS) in mapping over 100 loci associated with high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglyceride (TG) levels have advanced our understanding of genomic influences on common diseases [1-11]. However, the translation of such knowledge into clinical and public health applications requires exploration of the epidemiological architecture of these variants. Epidemiologic architecture describes the features of genetic associations in a population-based context that could act as modifiers of these associations [12]. Typical features include demographics (sex, age, and genetic ancestry) and various environmental exposures. The epidemiologic architecture of GWAS-identified genetic variants has largely been unexplored.

Differences in lipid profiles and metabolism between men and women have been well documented [13,14]. Pre-menopausal women have more favorable plasma lipid profiles than men, with lower levels of TG, total cholesterol (TC), and LDL-C, and higher of HDL-C [15-18]. Men and women also differ in lipoprotein particle concentration, subclass distribution, and sizes. For example, women have a larger average size of LDL and HDL particles than men [19,20]. Sex differences in lipid metabolism and lipoprotein kinetics are also well described [21]. However, the molecular basis of sexual dimorphism in lipid metabolism is poorly understood. Differences in sex hormones, body size and composition, and underlying genetic factors may each partially be implicated [17,18,22].

Sex-specific heritabilities of lipid traits have been previously reported [23,24]. Despite considerable advances in the identification of genetic variants influencing plasma lipid levels, few studies have examined the role of sex as a potential modifier of the effects of genetic variation on lipids. Sex-specific genetic associations may provide valuable insight into the factors responsible for the recognized sexual dimorphism in the plasma lipid profile, a major risk factor for cardiovascular disease, and, thus, may have significant public health and clinical relevance [25]. Accounting for the modifying effects of sex in genetic associations of lipid traits may also help in replication of results across studies and in generalization of findings across different populations.

We, therefore, investigated the sex-specific associations of 49 GWAS-identified SNPs with three common lipid traits (LDL-C, HDL-C, and TG) in the diverse cohorts of the Population Architecture using Genomics and Epidemiology (PAGE) study, established in 2008 by the National Human Genome Research Institute to investigate and characterize the epidemiologic architecture of GWAS-identified variants in diverse racial/ethnic groups [12]. We find significant heterogeneity by sex in three previously reported loci, and replicated heterogeneity by sex at nominal levels in three additional loci.

Results

Study populations

Table  1 illustrates the diversity of the PAGE study by racial/ethnic group or population, sex, age range, and years of data collection. There were 25,013 European Americans, 10,643 African Americans, 6,134 American Indians, 3,422 Mexican Americans/Hispanics, 827 Japanese/East Asians, and 200 Pacific Islanders/Native Hawaiians. Overall, there were more women (29,330) than men (16,909). Results are presented for the four largest racial/ethnic groups to maintain adequate sample sizes when investigating interactions. HDL-C, LDL-C, and ln(TG), stratified by sex, and racial/ethnic group, are shown in Table  2. Lipid levels are further stratified by cohort in Additional file 1: Table S1. Females had higher mean HDL-C levels than males across all four racial/ethnic groups, though LDL-C and ln(TG) profiles were similar (Table  2).

Table 1. Characteristics of PAGE study populations

Table 2. Mean lipid levels by sex and racial/ethnic group in the PAGE study

Additional file 1. This file contains Supplementary Tables 1-4 and Supplementary Figure 1 as described in the text.

Format: PDF Size: 1.4MB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Heterogeneity by sex in PAGE European Americans

Four SNPs across three different loci showed significant evidence of heterogeneity by sex at our established level of significance (phet < 0.001) in EAs, in HDL-C, ln(TG), and LDL-C levels (Table  3). The complete sex-stratified results for all 49 SNPs in the four largest PAGE racial/ethnic groups are shown in Additional file 1: Tables S2–S4.

Table 3. Loci with evidence of heterogeneity by sex (Phet< 0.001) in PAGE European Americans

In the APOA1/C3/A4/A5/BUD13 gene cluster, the C allele of rs28927680 was associated with increased ln(TG), but more strongly in males than females (Phet = 7.4×10-7). This SNP is in linkage disequilibrium with the other APOA1/C3/A4/A5 SNP presented in Table  3, rs3135506 (r2 = 0.83). In PLTP, rs7679 was strongly associated with HDL-C levels in females (p = 6.3×10-6), but there was no effect in males (p = 0.73) (phet = 9.9×10-4). In HMGCR, rs12654264 had a stronger effect on LDL-C levels in males (p = 8.3×10-14) than in females (p = 6.0×10-3) (phet = 3.1×10-5).

The only group with adequate power to detect interactions was EAs, and then only with SNPs that have larger allele frequencies and interaction effect sizes. For SNPs with a minor allele frequency of 0.15, the minimum detectable interaction effects with 80% power for HDL-C were 1.1 mg/dL in EAs, 1.8 mg/dL in AAs, 2.2 mg/dL in AIs, and 2.9 mg/dL in the Mexican population. These numbers represent the minimum detectable difference, comparing the effect size of a SNP on HDL-C levels in males vs. the effect size in females. The mean interaction effect for EAs across all SNPs for HDL-C was 0.35 mg/dL; only two SNPs with MAF > 0.15 had an interaction effect greater than 1.1 mg/dL (Additional file 1: Table S4). Additional file 1: Figure S1 shows the sample sizes required for minimum detectable interactions observable (in standard deviation units) for SNPs of varying allele frequencies.

Although not sufficiently powered to detect interaction in the other ethnic groups we include these results in the supplementary materials as a resource to build sample sizes large enough across these minority groups through meta-analysis. LDLR SNP rs6511720 met the criterion for heterogeneity in the Mexican American/Hispanic group (phet = 3.5×10-4), with a significant positive effect on ln(TG) levels in males and a significant negative effect in females (Additional file 1: Table S2). The same trend was also observed for American Indians for this SNP. In HNF4A, rs1800961 had a significant negative effect on HDL-C levels in AA males but no effect in AA females (phet = 2.6×10-4), although the minor allele frequency was only 0.01 (Additional file 1: Table S4). No other SNPs met the pre-specified criterion for declaring significant heterogeneity by sex for any of the three phenotypes in these racial/ethnic groups.

Replication of previously published reports of heterogeneity

Seventeen SNPs in our dataset have shown evidence of heterogeneity by sex for LDL-C, HDL-C, or TG in previously published genome-wide studies, or are in linkage disequilibrium (R2 > 0.2) with previously reported SNPs (Table  4). Heterogeneity by sex was replicated for PLTP and the BUD13 locus at our established level of significance (p < 0.001), as also shown in Table  3. (Note, heterogeneity by sex has been previously reported in HMGCR for total cholesterol levels, but we did not have results for total cholesterol and therefore this locus did not meet the criteria for inclusion in Table  4). Heterogeneity by sex was replicated for three additional loci (LPL, TRIB1, and GCKR) at a nominal p < 0.05 level. The binomial p-value for replicating 5 of 17 findings at p < 0.05 is 0.0009. For all five of these loci, the direction of the interaction was consistent with the previous report [26]. The effect was larger in males for four of the five loci, and stronger in females for one locus.

Table 4. Replication of previously reported heterogeneity by sex in PAGE

We also present the sample size that would have been required in order to have 80% power to replicate the previously reported interaction effect, given the allele frequency in PAGE and the R2 between the PAGE SNP and the prior SNP (Table  4). Notably, the required sample sizes to achieve 80% power were >100,000 for 14 of the 17 SNPs. For the one SNP where we were adequately powered to detect interaction (rs10401969 in CILP2), we observed a stronger effect in males than in females, consistent with the previous report, though the directions of effect for the T allele appear to be opposite (positive for both males and females in the prior study, but negative for both sexes in ours) [8].

Discussion

We have examined the sex-specific effects of 49 selected SNPs on circulating lipids in 46,349 PAGE study participants from four racial/ethnic groups. Four SNPs in three loci (HMGCR, PLTP, and APOA5/BUD13) displayed evidence of sex-SNP interactions for fasting lipid levels according to our pre-specified significance criterion. Heterogeneity by sex for lipid levels had been previously reported in other contexts for each of these loci. We were also able to replicate previously reported heterogeneity by sex at three additional loci (TRIB1, LPL, and GCKR) at nominal levels of significance, despite being underpowered to detect these interactions, with required sample sizes exceeding 100,000 in most cases.

Although we also analyzed these data in African Americans, American Indians, and Mexican American/Hispanics, we were underpowered to detect interactions in these groups, though two associations reached statistical significance. Sex-SNP interactions that are consistent across populations likely reflect the effects of the biological differences between men and women that are expected to be shared by all population groups. However, differences in power and LD patterns between population groups, as well as possible unrecognized racial/ethnic-specific effects on the sexual dimorphism of lipid levels may have obscured consistent sex-SNP interactions effects across population groups.

Two published GWAS and one gene-centric genome-wide study have tested for heterogeneity by sex for lipid levels among their associated loci. Three of 22 lipid-associated loci (HMGCR, LPL, and NCAN) exhibited evidence of heterogeneity by sex for either TC or HDL-C levels in a meta-analysis of >20,000 European individuals, though the criterion used for heterogeneity was not stated explicitly [1]. Four of 95 associated loci (LPL, CILP2, APOE, and ZNF664) exhibited evidence of heterogeneity by sex (Phet < 0.0005) for TC, TG, LDL-C, or HDL-C levels in a recent meta-analysis of >100,000 individuals of European ancestry, and seven additional loci were genome-wide significant in one sex but not the other [8]. In the most recent gene-centric genome-wide study, 44 SNPs of ~50,000 were defined as exhibiting heterogeneity by sex for TC, TG, HDL-C, or LDL-C [26].

Association data for seventeen of the previously reported loci from GWAS were available in PAGE. We attempted to replicate the prior findings for TG, HDL and LDL if we had a SNP in linkage disequilibrium (R2 > 0.20) with the reported SNP. We did not have association results for total cholesterol, and therefore did not attempt to replicate those findings from any of the prior studies. Five of the previously reported GWAS-identified sex-specific effects consistently replicated in the PAGE study for the same phenotype. There was some overlap in the sample between PAGE and the other meta-analyses; for example PAGE cohorts also analyzed in the gene-centric meta-analysis by Asselbergs et al. include ARIC, WHI, CARDIA, and CHS [26].

Power impacted our ability to detect significant sex effects, with sample sizes necessary for 80% power exceeding 100,000 in 14 of the 17 loci. Although we did not have large sample sizes, we did have the advantage of having previously unreported data for these loci in additional racial/ethnic groups. Of these five loci that were replicated in EAs, the four TG loci (APO gene cluster, GCKR, LPL, and TRIB1) had consistent directions of interaction in Mexican Americans/Hispanics, though effects were less consistent for the other racial/ethnic groups (Additional file 1: Table S2). It would be worthwhile to continue this area of investigation to determine whether interaction effects discovered in EAs tend to generalize to other racial/ethnic groups, as this evidence points to generalization in MA/H but not the other groups.

Several candidate gene studies have reported heterogeneity by sex for rs3135506 [27-30], and a nearby SNP was recently reported with heterogeneity by sex (rs11820589; R2 = 0.65 with rs28927680) [26]. In PAGE, this variant was associated with ln(TG) in all four major populations tested in the sex-combined meta-analysis [12]. Significant evidence of heterogeneity by sex in PAGE was observed for European Americans and Mexican Americans/Hispanics. In these two groups, significantly stronger effects in males compared with females were observed. A British study reported results consistent with the PAGE study [30]. A Turkish study reported significant associations between triglyceride levels and rs3135506, with a stronger association in women [29]. Sex-differentiated effects have also been reported for triglyceride levels and rs3135506 in a Brazilian population of European descent, again with the female effects stronger than male effects [27]. Klos et al. [28] had previously reported heterogeneity by sex for this SNP in the CARDIA study (data not included in the present study), where it was significantly associated with plasma TG levels in African-American females, but not males. There is also evidence that serum ApoA5 levels are correlated with triglyceride levels and HDL-C levels more strongly in females than in males [31]. Different results observed across different cohorts that represent various genetic ancestries highlight the complexities in replicating and ultimately interpreting sex differences in genetic association studies.

Aulchenko et al. [1] reported heterogeneity by sex for HMGCR rs3846662 for total cholesterol. We did not test for heterogeneity by sex for total cholesterol in the PAGE study. We did note, however, that HMGCR rs12654264, which is in strong LD with previously reported rs3846662 (r2 = 0.87), displayed significant heterogeneity by sex for LDL-C (Phet = 3.1×10-5), a trait highly correlated with total cholesterol. In the PAGE study, however, the genetic effect was greater in males (β = −4.03) compared with females (β = −1.11), which does not replicate the findings of Aulchenko et al. [1]. The product of the HMGCR gene (3-hydroxy-3-methylglutaryl-CoA reductase) is the rate-limiting enzyme of the cholesterol biosynthesis pathway and the target of statins, a class of drugs widely used for the treatment of high cholesterol. Sequence variants of this gene have been associated with variation in response to statin therapy [32]. Among patients with asymptomatic plaques in the carotid artery from the Malmö Diet and Cancer-Cardiovascular Cohort, rs12654264 was associated with reduction in LDL-C levels in response to fluvastatin treatment in men but not women [33].

Heterogeneity by sex for HDL-C levels for PLTP (phospholipid transfer protein) has been previously reported, with consistent findings [26]. Studies have provided evidence that PLTP activity may affect HDL particle size [34]. In PAGE, the major allele of rs7679 was associated with higher HDL-C levels in women only.

The locus with the most consistent evidence for heterogeneity by sex across the studies is LPL, or lipoprotein lipase. Different SNPs in this gene exhibited heterogeneity by sex for HDL levels in two prior studies, with a larger effect in males [1,8]. In PAGE and in Asselbergs et al., LPL exhibited heterogeneity by sex for TG levels, also with a stronger effect in males [26]. LPL (lipoprotein lipase) is the rate-limiting enzyme for hydrolysis of triglycerides in lipoproteins. Polymorphisms and mutations in LPL have been associated with lipid metabolism disorders. Hormone levels have been shown to affect regulation of LPL, including thyroid hormone, estrogen, and testosterone [35].

Strengths and limitations

The diversity of the PAGE study potentially enabled us to examine potential sex-effects across populations, though we were underpowered to detect interactions in three of the four racial/ethnic groups. Physiological, anatomical, or even behavioral differences between men and women that may modify the effects of SNPs on lipid metabolism are expected to be largely shared across racial/ethnic groups. Hence, the consistent effects across racial/ethnic groups described above provide added support to the findings reported here.

Some limitations must be acknowledged, including sample size. Power to detect interaction effects typically requires substantially larger samples than those for main effects [36,37]. Sample size was greatest for European-Americans and, not unexpectedly, most evidence of sex-specific effects in this study was observed in this population. The required sample sizes to detect interaction for these loci (generally exceeding 100,000, and sometimes exceeding 1,000,000) should alert other investigators to the difficulty of replicating interaction effects.

The PAGE study cohorts differ in many aspects, including study design, period of collection, demographics and cardiovascular and metabolic risk factors of the participants. These differences may have further reduced our power to detect significant modifying effects of sex on SNPs-lipids associations. Indeed, sex differences not only represent biological differences between men and women but also encompass or are confounded by social and behavioral differences between the two sexes. Difficulties in assessing and accounting for such factors in a consistent and accurate manner across the multiple cohorts likely further reduced our power to detect interaction with genetic factors.

This difficulty of replicating heterogeneity by sex is compounded by the nature of the tested SNPs, which are likely to be in LD with the causal variant(s) but not themselves causal. Differences in LD patterns across studies and racial/ethnic groups may have hampered our ability to detect consistent sex modification effects in some population groups.

Claims of modification by sex have been difficult to replicate for most complex diseases and many studies lack the proper documentation for the claim of significant sex-effects [38]. In this study, we note that the three loci which met our a priori criterion for significant interaction had been previously reported to display heterogeneity by sex, suggesting that the initial reports for these loci were not attributable to type I error.

Differences in sex hormone levels have been hypothesized to play a role in the sexual dimorphism of circulating lipids. In this study, we did not examine whether menopausal status modified the association of these SNPs with lipid levels in women. Future such investigations may help shed light on the biological basis of the sex-specific associations reported here.

Conclusions

Using a rigorous methodology and the diverse populations of the PAGE study, we have confirmed previously reported heterogeneity by sex for lipid levels for six loci. Genotype-sex interactions may represent an important source of genetic variation that may contribute to the “missing heritability” of complex traits. Although challenging, assessment of sex-specific associations should be more widely considered in order to characterize the genetic architecture of complex, sexually dimorphic traits, such as lipids.

Methods

Study populations and phenotypes

The study population included 46,239 individuals from cohorts which are part of the PAGE study, a collaborative program across four large population-based studies or consortia, including EAGLE (Epidemiologic Architecture for Genes Linked to Environment), based on three National Health and Nutrition Examination Surveys (NHANES) [39]; the Multiethnic Cohort (MEC) [40]; the Women’s Health Initiative (WHI) [41]; and the Causal variants Across the Life Course (CALiCo) consortium, which encompasses five studies: Atherosclerosis Risk in Communities (ARIC) [42], Coronary Artery Risk in Young Adults (CARDIA) [43], the Cardiovascular Health Study (CHS) [44], Strong Heart Family Study (SHFS) and Strong Heart Study (SHS) [45]. Details about the design of the PAGE study have been previously published [12]. All participants were consented and all studies were approved by Institutional Review Boards at their respective sites.

Serum HDL-C, triglycerides, and total cholesterol were measured using standard enzymatic methods. LDL-C was calculated using the Friedewald equation, with missing values assigned for samples with triglyceride levels greater than 400 mg/dl. For PAGE cohorts with longitudinal data, measurements from the baseline examination were used in the analyses.

The PAGE study participants who were less than 18 years of age and those fasting for less than 8 hours prior to the blood draw were excluded from analyses. Participants with triglycerides values >1,000 mg/dl were excluded from analyses of that trait. The distribution of triglyceride levels was skewed and thus values were natural log transformed prior to analysis. A further description of the study design, methods of data collection and participants’ characteristics for each of the cohorts have been presented elsewhere [46].

SNP selection and genotyping

Detailed methods of SNP selection and genotyping have been described previously [12]. Briefly, a total of 52 SNPs previously associated with HDL-C, LDL-C, and/or triglycerides in published candidate gene and genome-wide association studies (through 2008) were targeted for genotyping in two or more PAGE cohorts. Of these, three (CETP rs1800775, APOE rs429358, and APOE rs7412) failed at all PAGE study sites that attempted genotyping; therefore, a total of 49 SNPs were included in this analysis. Genotyping was performed by each of the four PAGE studies using commercially available genotyping arrays (Affymetrix 6.0, Illumina 370CNV BeadChip), custom mid- and low-throughput assays (TaqMan, Sequenom, Illumina GoldenGate on the BeadXpress), or a combination thereof. Quality control was implemented at each study site independently. Only SNPs with high call rates (>95%) were included in the analyses. In addition, all PAGE study sites genotyped 360 DNA samples from the International HapMap Project for concordance analysis.

Statistical methods

Cohort-specific analyses

Statistical analyses were performed separately by each cohort following the same analysis plan. Within each racial/ethnic group and sex stratum, linear regression was used to evaluate the association of each SNP with HDL-C, LDL-C, or natural log-transformed TG levels, assuming an additive genetic model. Models were adjusted for age and field center (for multi-center studies). Previous PAGE study results for lipids demonstrated that further adjustment for body mass index, current smoking, type 2 diabetes, post-menopausal status, current hormone use, myocardial infarction, and ancestry using principal components, did not meaningfully impact the association of these SNPs with lipid levels [46]. Analyses were performed without regard to lipid lowering medication status given that relatively few participants (<10%) reported such medication use, and their inclusion did not appreciably alter results in previous PAGE lipids work [46].

Meta-analyses

Sex-specific beta coefficients were combined within each racial/ethnic group using an inverse variance-weighted fixed-effects meta-analysis, using METAL software [47]. For each SNP, significance of association within each stratum was evaluated using a 1-df Chi-square test. Heterogeneity of effects between the sexes was then evaluated using the following 1-df Chi-square test: (βmf)/(SEm2 + SEf2) ~ χ21, where βm and βf represent the meta-analyzed effect estimates among males and females, respectively. In addition, for each SNP, sex-differentiated tests of association were performed using a Chi-square test (2 df) as described by Magi et al. [48] in which the 1-df Chi-square statistics for the SNP-phenotype associations for males and females are summed. This method yields a single P-value of association but permits the effect estimates for males and females to differ. It also permits inclusion of single-sex studies such as WHI that otherwise would not be able to contribute to a traditional interaction meta-analysis. Reported P-values were not adjusted for multiple testing.

Reported P-values include Pcomb, the genotype-phenotype association p-value not allowing for a sex-specific association; Pdiff, the genotype-phenotype association P-value allowing for different effects by sex; and Phet, the P-value for heterogeneity of sex effects. The a priori criterion for heterogeneity was set at Phet < 0.001 (a Bonferroni correction for the number of SNPs examined; 0.05/49).

Sample size and power calculations

Sample size and power for gene-environment interactions were calculated using Quanto software (http://hydra.usc.edu/gxe/ webcite). To calculate sample sizes required to replicate previous findings, we used the following assumptions: power = 80%; alpha = 0.05 (two-sided); interaction effect size = (male effect – female effect, in standard deviation units); minor allele frequency = the minor allele frequency in the PAGE population; additive genetic model. To calculate minimum detectable effect sizes, we set sample size equal to the PAGE sample size but kept the other assumptions. Because none of the other racial/ethnic groups in PAGE had adequate power to detect gene-environment interactions, only the European-American results are presented in the main text. Results for the other groups are presented in the supplementary materials.

Abbreviations

AA: African-American; AI: American Indian; ARIC: Atherosclerosis Risk In Communities; CAF: Coded allele frequency; CALICO: Causal Variants Across the Life Course; CARDIA: Coronary Artery Risk Development in Young Adults; CHS: Cardiovascular Health Study; df: Degrees of freedom; EA: European-American; EAGLE: Epidemiologic Architecture for Genes Linked to Environment; GWAS: Genome-wide association study; HDL-C: High density lipoprotein-cholesterol; LDL-C: Low density lipoprotein-cholesterol; MA/H: Mexican American/Hispanic; MEC: Multi-Ethnic Cohort; PAGE: Population Architecture using Genomics and Epidemiology; Pcomb: P-value for association, not allowing for different effects by sex; Pdiff: P-value for association, allowing different effect estimates for males and females; Phet: P-value for heterogeneity; R/E: Race/Ethnicity; TC: Total cholesterol; TG: Triglycerides; SE: Standard error; SHS: Strong Heart Study; SNP: Single Nucleotide Polymorphism; WHI: Women’s Health Initiative.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

KCT and MF designed the meta-analysis, conducted the meta-analysis and power analysis, and drafted and edited the manuscript. CLC, LD, PB, SAC, LH, and FRS assisted with planning the study, and interpretation of results, and with drafting and editing the manuscript. LRW, RVS, PMQ, KCJ, BEH, JH, NF, CBE, DJD, BC, and IC conducted the single-study analyses in the original cohorts and participated in discussions regarding the manuscript. CSC, KBG, GA, JLA, CH, LLM, and CK oversaw the single-study analyses and participated in discussions regarding the manuscript. DCC, SB, and KEN assisted with the design of the meta-analysis, and assisted with the drafting and editing the manuscript. All authors have read and approved the final manuscript.

Acknowledgements

The contents of this paper are solely the responsibility of the authors and do not necessarily represent the official views of the NIH. The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the Centers for Disease Control and Prevention. The opinions expressed in this paper are those of the author(s) and do not necessarily reflect the views of the Indian Health Service.

The PAGE consortium thanks the staff and participants of all PAGE studies for their important contributions. The complete list of PAGE members can be found at http://www.pagestudy.org.

The authors thank the WHI investigators and staff for their dedication, and the study participants for making the program possible. A full listing of WHI investigators can be found at: http://www.whiscience.org/publications/WHI_investigators_shortlist.pdf.

EAGLE would like to thank Dr. Geraldine McQuillan and Jody McLean for their help in accessing the Genetic NHANES data. The Vanderbilt University Center for Human Genetics Research, Computational Genomics Core, provided computational and/or analytical support for this work. The EAGLE/NHANES DNA samples are stored and plated by the Vanderbilt DNA Resources Core. Genotyping was performed by Ping Mayo, Melissa Allen, and Dr. Nathalie Schnetz-Boutaud in the laboratory of Dr. Jonathan Haines and Hailing Jin and Nila Gillani under the direction of Dr. Holli Dilks in the Vanderbilt DNA Resources Core.

References

  1. Aulchenko YS, Ripatti S, Lindqvist I, Boomsma D, Heid IM, Pramstaller PP, Penninx BW, Janssens AC, Wilson JF, Spector T: Loci influencing lipid levels and coronary heart disease risk in 16 European population cohorts.

    Nat Genet 2009, 41(1):47-55. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Chasman DI, Pare G, Mora S, Hopewell JC, Peloso G, Clarke R, Cupples LA, Hamsten A, Kathiresan S, Malarstig A: Forty-three loci associated with plasma lipoprotein size, concentration, and cholesterol content in genome-wide analysis.

    PLoS Genet 2009, 5(11):e1000730. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Kathiresan S, Manning AK, Demissie S, D’Agostino RB, Surti A, Guiducci C, Gianniny L, Burtt NP, Melander O, Orho-Melander M: A genome-wide association study for blood lipid phenotypes in the Framingham Heart Study.

    BMC Med Genet 2007, 8(1):S17. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Kathiresan S, Melander O, Guiducci C, Surti A, Burtt NP, Rieder MJ, Cooper GM, Roos C, Voight BF, Havulinna AS: Six new loci associated with blood low-density lipoprotein cholesterol, high-density lipoprotein cholesterol or triglycerides in humans.

    Nat Genet 2008, 40(2):189-197. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Kathiresan S, Willer CJ, Peloso GM, Demissie S, Musunuru K, Schadt EE, Kaplan L, Bennett D, Li Y, Tanaka T: Common variants at 30 loci contribute to polygenic dyslipidemia.

    Nat Genet 2009, 41(1):56-65. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Kooner JS, Chambers JC, Aguilar-Salinas CA, Hinds DA, Hyde CL, Warnes GR, Gomez Perez FJ, Frazer KA, Elliott P, Scott J: Genome-wide scan identifies variation in MLXIPL associated with plasma triglycerides.

    Nat Genet 2008, 40(2):149-151. PubMed Abstract | Publisher Full Text OpenURL

  7. Sabatti C, Service SK, Hartikainen AL, Pouta A, Ripatti S, Brodsky J, Jones CG, Zaitlen NA, Varilo T, Kaakinen M: Genome-wide association analysis of metabolic traits in a birth cohort from a founder population.

    Nat Genet 2009, 41(1):35-46. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Teslovich TM, Musunuru K, Smith AV, Edmondson AC, Stylianou IM, Koseki M, Pirruccello JP, Ripatti S, Chasman DI, Willer CJ: Biological, clinical and population relevance of 95 loci for blood lipids.

    Nature 2010, 466(7307):707-713. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Wallace C, Newhouse SJ, Braund P, Zhang F, Tobin M, Falchi M, Ahmadi K, Dobson RJ, Marcano AC, Hajat C: Genome-wide association study identifies genes for biomarkers of cardiovascular disease: serum urate and dyslipidemia.

    Am J Hum Genet 2008, 82(1):139-149. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Waterworth DM, Ricketts SL, Song K, Chen L, Zhao JH, Ripatti S, Aulchenko YS, Zhang W, Yuan X, Lim N: Genetic variants influencing circulating lipid levels and risk of coronary artery disease.

    Arterioscler Thromb Vasc Biol 2010, 30(11):2264-2276. PubMed Abstract | Publisher Full Text OpenURL

  11. Willer CJ, Sanna S, Jackson AU, Scuteri A, Bonnycastle LL, Clarke R, Heath SC, Timpson NJ, Najjar SS, Stringham HM: Newly identified loci that influence lipid concentrations and risk of coronary artery disease.

    Nat Genet 2008, 40(2):161-169. PubMed Abstract | Publisher Full Text OpenURL

  12. Matise T, Ambite J, Buyske S, Cole S, Crawford DC, Haiman CA, Heiss G, Kooperberg C, Le Marchand L, Manolio TA: The next page in understanding complex traits: Study design and analysis of Population Architecture using Genomics and Epidemiology.

    Am J Epidemiology 2011, 174(7):849-859. Publisher Full Text OpenURL

  13. Carroll MD, Lacher DA, Sorlie PD, Cleeman JI, Gordon DJ, Wolz M, Grundy SM, Johnson CL: Trends in serum lipids and lipoproteins of adults, 1960–2002.

    JAMA 2005, 294(14):1773-1781. PubMed Abstract | Publisher Full Text OpenURL

  14. Mittendorfer B: Sexual dimorphism in human lipid metabolism.

    J Nutr 2005, 135(4):681-686. PubMed Abstract | Publisher Full Text OpenURL

  15. Gardner CD, Winkleby MA, Fortmann SP: Population frequency distribution of non-high-density lipoprotein cholesterol (Third National Health and Nutrition Examination Survey [NHANES III], 1988–1994).

    Am J Cardiol 2000, 86(3):299-304. PubMed Abstract | Publisher Full Text OpenURL

  16. Gostynski M, Gutzwiller F, Kuulasmaa K, Doring A, Ferrario M, Grafnetter D, Pajak A: Analysis of the relationship between total cholesterol, age, body mass index among males and females in the WHO MONICA Project.

    Int J Obes Relat Metab Disord 2004, 28(8):1082-1090. PubMed Abstract | Publisher Full Text OpenURL

  17. Regitz-Zagrosek V, Lehmkuhl E, Mahmoodzadeh S: Gender aspects of the role of the metabolic syndrome as a risk factor for cardiovascular disease.

    Gend Med 2007, (4):S162-177. OpenURL

  18. Williams CM: Lipid metabolism in women.

    Proc Nutr Soc 2004, 63(1):153-160. PubMed Abstract | Publisher Full Text OpenURL

  19. Freedman DS, Otvos JD, Jeyarajah EJ, Shalaurova I, Cupples LA, Parise H, D’Agostino RB, Wilson PW, Schaefer EJ: Sex and age differences in lipoprotein subclasses measured by nuclear magnetic resonance spectroscopy: the Framingham Study.

    Clin Chem 2004, 50(7):1189-1200. PubMed Abstract | Publisher Full Text OpenURL

  20. Johnson JL, Slentz CA, Duscha BD, Samsa GP, McCartney JS, Houmard JA, Kraus WE: Gender and racial differences in lipoprotein subclass distributions: the STRRIDE study.

    Atherosclerosis 2004, 176(2):371-377. PubMed Abstract | Publisher Full Text OpenURL

  21. Wang X, Magkos F, Mittendorfer B: Sex differences in lipid and lipoprotein metabolism: it’s not just about sex hormones.

    J Clin Endocrinol Metab 2011, 96(4):885-893. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Seidell JC, Cigolini M, Charzewska J, Ellsinger BM, Bjorntorp P, Hautvast JG, Szostak W: Fat distribution and gender differences in serum lipids in men and women from four European communities.

    Atherosclerosis 1991, 87(2–3):203-210. PubMed Abstract OpenURL

  23. Miljkovic I, Yerges-Armstrong LM, Kuller LH, Kuipers AL, Wang X, Kammerer CM, Nestlerode CS, Bunker CH, Patrick AL, Wheeler VW: Association analysis of 33 lipoprotein candidate genes in multi-generational families of African ancestry.

    J Lipid Res 2010, 51(7):1823-1831. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Weiss LA, Pan L, Abney M, Ober C: The sex-specific genetic architecture of quantitative traits in humans.

    Nat Genet 2006, 38(2):218-222. PubMed Abstract | Publisher Full Text OpenURL

  25. Arnold AP: Promoting the understanding of sex differences to enhance equity and excellence in biomedical science.

    Biol Sex Differ 2010, 1(1):1. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  26. Asselbergs FW, Guo Y, van Iperen EP, Sivapalaratnam S, Tragante V, Lanktree MB, Lange LA, Almoguera B, Appelman YE, Barnard J: Large-scale gene-centric meta-analysis across 32 studies identifies multiple lipid loci.

    Am J Hum Genet 2012, 91(5):823-838. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. De Andrade FM, Maluf SW, Schuch JB, Voigt F, Barros AC, Lucatelli JF, Hutz MH: The influence of the S19W SNP of the APOA5 gene on triglyceride levels in southern Brazil: interactions with the APOE gene, sex and menopause status.

    Nutr Metab Cardiovasc Dis 2011, 21(8):584-590. PubMed Abstract | Publisher Full Text OpenURL

  28. Klos KL, Hamon S, Clark AG, Boerwinkle E, Liu K, Sing CF: APOA5 polymorphisms influence plasma triglycerides in young, healthy African Americans and whites of the CARDIA Study.

    J Lipid Res 2005, 46(3):564-571. PubMed Abstract | Publisher Full Text OpenURL

  29. Komurcu-Bayrak E, Onat A, Poda M, Humphries SE, Palmen J, Guclu F, Can G, Erginel-Unaltuna N: Gender-modulated impact of apolipoprotein A5 gene (APOA5) -1131T > C and c.56C > G polymorphisms on lipids, dyslipidemia and metabolic syndrome in Turkish adults.

    Clin Chem Lab Med 2008, 46(6):778-784. PubMed Abstract | Publisher Full Text OpenURL

  30. Olano-Martin E, Abraham EC, Gill-Garrison R, Valdes AM, Grimaldi K, Tang F, Jackson KG, Williams CM, Minihane AM: Influence of apoA-V gene variants on postprandial triglyceride metabolism: impact of gender.

    J Lipid Res 2008, 49(5):945-953. PubMed Abstract | Publisher Full Text OpenURL

  31. Zhao SP, Hu S, Li J, Hu M, Liu Q, Wu LJ, Zhang T: Association of human serum apolipoprotein A5 with lipid profiles affected by gender.

    Clin Chim Acta 2007, 376(1–2):68-71. PubMed Abstract | Publisher Full Text OpenURL

  32. Medina MW: The relationship between HMGCR genetic variation, alternative splicing, and statin efficacy.

    Discov Med 2010, 9(49):495-499. PubMed Abstract OpenURL

  33. Hamrefors V, Orho-Melander M, Krauss RM, Hedblad B, Almgren P, Berglund G, Melander O: A gene score of nine LDL and HDL regulating genes is associated with fluvastatin-induced cholesterol changes in women.

    J Lipid Res 2010, 51(3):625-634. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Albers JJ, Vuletic S, Cheung MC: Role of plasma phospholipid transfer protein in lipid and lipoprotein metabolism.

    Biochim Biophys Acta 2012, 1821(3):345-357. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Wang H, Eckel RH: Lipoprotein lipase: from gene to obesity.

    Am J Physiol Endocrinol Metab 2009, 297(2):E271-288. PubMed Abstract | Publisher Full Text OpenURL

  36. Greenland S: Tests for interaction in epidemiologic studies: a review and a study of power.

    Stat Med 1983, 2(2):243-251. PubMed Abstract | Publisher Full Text OpenURL

  37. Hunter DJ: Gene-environment interactions in human diseases.

    Nat Rev Genet 2005, 6(4):287-298. PubMed Abstract | Publisher Full Text OpenURL

  38. Patsopoulos NA, Tatsioni A, Ioannidis JP: Claims of sex differences: an empirical assessment in genetic associations.

    JAMA 2007, 298(8):880-893. PubMed Abstract | Publisher Full Text OpenURL

  39. Centers for Disease Control and Prevention: National Health and Nutrition Examination Survey (NHANES) DNA Samples: Guidelines for Proposals to Use Samples and Cost Schedule.

    Fed Regist 2010, (75):32191-32195. OpenURL

  40. Kolonel LN, Altshuler D, Henderson BE: The multiethnic cohort study: exploring genes, lifestyle and cancer risk.

    Nat Rev Cancer 2004, 4(7):519-527. PubMed Abstract | Publisher Full Text OpenURL

  41. The Women’s Health Initiative Study Group: Design of the Women’s Health Initiative clinical trial and observational study.

    Control Clin Trials 1998, 19(1):61-109. PubMed Abstract | Publisher Full Text OpenURL

  42. The ARIC investigators: The Atherosclerosis Risk in Communities (ARIC) Study: design and objectives.

    Am J Epidemiol 1989, 129(4):687-702. PubMed Abstract | Publisher Full Text OpenURL

  43. Friedman GD, Cutter GR, Donahue RP, Hughes GH, Hulley SB, Jacobs DR Jr, Liu K, Savage PJ: CARDIA: study design, recruitment, and some characteristics of the examined subjects.

    J Clin Epidemiol 1988, 41(11):1105-1116. PubMed Abstract | Publisher Full Text OpenURL

  44. Fried LP, Borhani NO, Enright P, Furberg CD, Gardin JM, Kronmal RA, Kuller LH, Manolio TA, Mittelmark MB, Newman A: The Cardiovascular Health Study: design and rationale.

    Ann Epidemiol 1991, 1(3):263-276. PubMed Abstract | Publisher Full Text OpenURL

  45. Lee ET, Welty TK, Fabsitz R, Cowan LD, Le NA, Oopik AJ, Cucchiara AJ, Savage PJ, Howard BV: The Strong Heart Study. A study of cardiovascular disease in American Indians: design and methods.

    Am J Epidemiol 1990, 132(6):1141-1155. PubMed Abstract | Publisher Full Text OpenURL

  46. Dumitrescu L, Carty CL, Taylor KC, Schumacher FR, Hindorff LA, Ambite JL, Anderson G, Best LG, Brown-Gentry K, Buzkova P: Genetic determinants of lipid traits in diverse populations from the population architecture using genomics and epidemiology (PAGE) study.

    PLoS Genet 2011, 7(6):e1002138. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. Willer CJ, Li Y, Abecasis GR: METAL: fast and efficient meta-analysis of genomewide association scans.

    Bioinformatics 2010, 26(17):2190-2191. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Magi R, Lindgren CM, Morris AP: Meta-analysis of sex-specific genome-wide association studies.

    Genet Epidemiol 2010, 34(8):846-853. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL