Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Gene-Gene Combination Effect and Interactions among ABCA1, APOA1, SR-B1, and CETP Polymorphisms for Serum High-Density Lipoprotein-Cholesterol in the Japanese Population

  • Akihiko Nakamura,

    Affiliation Department of International Island and Community Medicine, Kagoshima University Graduate School of Medical and Dental Sciences, Kagoshima, Japan

  • Hideshi Niimura,

    Affiliation Department of International Island and Community Medicine, Kagoshima University Graduate School of Medical and Dental Sciences, Kagoshima, Japan

  • Kazuyo Kuwabara,

    Affiliation Education Center for Doctors in Remote Islands and Rural Areas, Kagoshima University Graduate School of Medical and Dental Sciences, Kagoshima, Japan

  • Toshiro Takezaki ,

    takezaki@m.kufm.kagoshima-u.ac.jp

    Affiliation Department of International Island and Community Medicine, Kagoshima University Graduate School of Medical and Dental Sciences, Kagoshima, Japan

  • Emi Morita,

    Affiliation Department of Preventive Medicine, Nagoya University Graduate School of Medicine, Nagoya, Japan

  • Kenji Wakai,

    Affiliation Department of Preventive Medicine, Nagoya University Graduate School of Medicine, Nagoya, Japan

  • Nobuyuki Hamajima,

    Affiliation Department of Healthcare Administration, Nagoya University Graduate School of Medicine, Nagoya, Japan

  • Yuichiro Nishida,

    Affiliation Department of Preventive Medicine, Faculty of Medicine, Saga University, Saga, Japan

  • Tanvir Chowdhury Turin,

    Affiliations Department of Health Science, Shiga University of Medical Science Otsu, Japan, Department of Medicine, University of Calgary, Calgary, Canada

  • Sadao Suzuki,

    Affiliation Department of Public Health, Nagoya City University Graduate School of Medical Sciences, Nagoya, Japan

  • Keizo Ohnaka,

    Affiliation Department of Geriatric Medicine, Graduate School of Medical Sciences, Kyushu University, Fukuoka, Japan

  • Hirokazu Uemura,

    Affiliation Department of Preventive Medicine, Institute of Health Biosciences, University of Tokushima Graduate School, Tokushima, Japan

  • Etsuko Ozaki,

    Affiliation Department of Epidemiology for Community Health and Medicine, Kyoto Prefectural University of Medicine, Kyoto, Japan

  • Satoyo Hosono,

    Affiliation Division of Epidemiology and Prevention, Aichi Cancer Center Research Institute, Nagoya, Japan

  • Haruo Mikami,

    Affiliation Division of Cancer Registry, Prevention and Epidemiology, Chiba Cancer Center, Chiba, Japan

  • Michiaki Kubo,

    Affiliation Laboratory for Genotyping Development, Center for Genomic Medicine, RIKEN, Yokohama, Japan

  •  [ ... ],
  • Hideo Tanaka

    Affiliation Division of Epidemiology and Prevention, Aichi Cancer Center Research Institute, Nagoya, Japan

  • [ view all ]
  • [ view less ]

Abstract

Background/Objective

Gene-gene interactions in the reverse cholesterol transport system for high-density lipoprotein-cholesterol (HDL-C) are poorly understood. The present study observed gene-gene combination effect and interactions between single nucleotide polymorphisms (SNPs) in ABCA1, APOA1, SR-B1, and CETP in serum HDL-C from a cross-sectional study in the Japanese population.

Methods

The study population comprised 1,535 men and 1,515 women aged 35–69 years who were enrolled in the Japan Multi-Institutional Collaborative Cohort (J-MICC) Study. We selected 13 SNPs in the ABCA1, APOA1, CETP, and SR-B1 genes in the reverse cholesterol transport system. The effects of genetic and environmental factors were assessed using general linear and logistic regression models after adjusting for age, sex, and region.

Principal Findings

Alcohol consumption and daily activity were positively associated with HDL-C levels, whereas smoking had a negative relationship. The T allele of CETP, rs3764261, was correlated with higher HDL-C levels and had the highest coefficient (2.93 mg/dL/allele) among the 13 SNPs, which was statistically significant after applying the Bonferroni correction (p<0.001). Gene-gene combination analysis revealed that CETP rs3764261 was associated with high HDL-C levels with any combination of SNPs from ABCA1, APOA1, and SR-B1, although no gene-gene interaction was apparent. An increasing trend for serum HDL-C was also observed with an increasing number of alleles (p<0.001).

Conclusions

The present study identified a multiplier effect from a polymorphism in CETP with ABCA1, APOA1, and SR-B1, as well as a dose-dependence according to the number of alleles present.

Introduction

Coronary heart disease and cerebrovascular disease are the leading causes of mortality in both high- and low-income countries [1], [2]. A low level of high-density lipoprotein-cholesterol (HDL-C) is an important risk factor for these diseases [3], [4]. Environmental and genetic factors influence HDL-C levels; however, no specific treatments are available. A recent Mendelian randomization study reported that a genetic mechanism for raising serum HDL-C did not seem to lower the risk of myocardial infarction [5]. Further study is still required.

Several environmental factors have been reported to affect HDL-C levels. Alcohol consumption, habitual exercise, and high consumption of eggs and fish are positively associated with HDL-C levels [6][9]. Smoking leads to decreased levels of HDL-C [10]. Genetic factors, with or without environmental influence, also affect HDL-C levels. Several gene polymorphisms, especially in the genes encoding enzymes involved in the reverse cholesterol transport (RCT) system, have been reported to be associated with HDL-C levels [11][17]. The interaction between genetic and environmental factors has been investigated, but the reports have shown inconsistent results [18][20]. Furthermore, gene-gene interactions in the RCT system may modulate HDL-C levels, but few studies have been able clarify this effect, with the exception of one report from among the USA population [21].

HDL-C exerts an anti-atherogenic effect through several mechanisms such as anti-inflammation, anti-oxidation of low-density lipoprotein-cholesterol (LDL-C), and inhibition of vascular endothelial cell apoptosis. The RCT system also plays an important role in these processes [22]. The RCT system is involved in the transportation of cholesterol from the peripheral tissues to the liver, where the cholesterol is secreted into bile. ATP-binding-cassette A1 (ABCA1), apolipoprotein A-1 (ApoA-1), lecithin cholesterol acyltransferase (LCAT), cholesteryl ester transfer protein (CETP), and scavenger receptor class B1 (SR-B1) play important roles in the RCT system.

We are currently conducting the Japan Multi-Institutional Collaborative Cohort (J-MICC) study by compiling data from a relatively large number of subjects recruited from the Japanese population together with environmental data and DNA samples. To clarify the effects of gene-gene combination and interactions on HDL-C levels, we conducted a cross-sectional study by using data from the J-MICC study and focusing on gene polymorphisms in ABCA1, APOA1, CETP, and SR-B1 in the RCT system.

Subjects and Methods

Ethics Statement

All participants provided written informed consent. The ethics committees of all participating institutes and universities approved the protocol. All data and samples were sent to the Nagoya University School of Medicine as linkable anonymizing data and de-identified samples. New identification numbers were then given for the combined dataset, which was linked to the analyzed data of gene polymorphisms. These data were sent to each collaborator without identifying lists.

The participating institutes and universities included: (1) Department of Preventive Medicine, Nagoya University Graduate School of Medicine, Nagoya, Japan; (2) Department of International Island and Community Medicine, Kagoshima University Graduate School of Medical and Dental Sciences, Kagoshima, Japan; (3) Department of Preventive Medicine, Faculty of Medicine, Saga University, Saga, Japan; (4) Department of Health Science, Shiga University of Medical Science, Otsu, Japan; (5) Department of Public Health, Nagoya City University Graduate School of Medical Sciences, Nagoya, Japan; (6) Department of Geriatric Medicine, Graduate School of Medical Sciences, Kyushu University, Fukuoka, Japan; (7) Department of Preventive Medicine, Institute of Health Biosciences, University of Tokushima Graduate School, Tokushima, Japan; (8) Department of Epidemiology for Community Health and Medicine, Kyoto Prefectural University of Medicine, Kyoto, Japan; (9) Division of Epidemiology and Prevention, Aichi Cancer Center Research Institute, Nagoya, Japan; (10) Division of Cancer Registry, Prevention and Epidemiology, Chiba Cancer Center, Chiba, Japan; and (11) Laboratory for Genotyping Development, Center for Genomic Medicine, RIKEN, Yokohama, Japan.

Study Population

The J-MICC study has been conducted in 10 regions of Japan by 10 research institutes and universities since 2005, as described previously [23], [24]. In brief, the cohort participants were enrolled from the community through invitations mailed or leaflets distributed (3 regions) to patients on their first visit to a cancer hospital (1 region) or at health checkups (6 regions). First, we recruited 5,108 participants at the baseline of the J-MICC study (Figure 1). From these, we excluded the participants from whom we did not receive appropriate informed consent (n = 8), sufficient DNA (n = 442), questionnaire data (n = 9), or local government registration of residence in the study region (n = 7); anyone who had declined follow-up visits (n = 2); anyone who had withdrawn from the study (n = 1); and those who were under 35 or over 69 years of age (n = 120). Then, we proceeded with the SNPs analysis for the residual 4,519 subjects (2,124 men and 2,395 women). Furthermore, as the present study examined the association between HDL-C and SNPs, we also excluded those who had no HDL-C data (HDL-C examination was not included at one of the Cancer Center study region and at one of the community study regions; n = 1,088); those who had a history of liver cirrhosis (n = 8); those with low albumin levels as per blood examination (<3.5 g/dL); those with a low A/G ratio (<1.0; n = 5); those with a history of dyslipidemia with medication (n = 306); and those who had stopped drinking alcoholic beverages (excluded because most of them had stopped drinking owing to diseases associated with liver dysfunction; n = 62). Eventually, 3,050 subjects (1,535 men and 1,515 women) in 8 regions were deemed eligible for the present study. These regions were located in the western part of Japan, including the Amami Islands (Figure 2).

thumbnail
Figure 1. Diagram to list individuals who were excluded from the study sample.

https://doi.org/10.1371/journal.pone.0082046.g001

Genotyping

First, we selected 5 genes encoding enzymes in the RCT system. Using these genes, we then selected 15 single nucleotide polymorphisms (SNPs) that have been reported to be associated with HDL-C levels in previous studies [11][17]: ABCA1-565C>T (rs2422493), R1587K (rs2230808), -273G>C (rs1800976), V771M (rs2066718), -17C>G (rs2740483), and V825I (rs2066715); APOA1 A61T (rs12718465); LCAT (rs4986970); CETP Taq1B (rs708272), G/T (rs3764261), I405V (rs5882), and -629C>A (rs1800775); SR-B1 A>G (rs3782287), A350A (rs5888), and V135I (rs5891). The selected SNPs were genotyped by the multiplex polymerase chain reaction-based Invader assay at the Laboratory for Genotyping Development, Center for Genomic Medicine, RIKEN, as described previously [23]. Because no minor allele was found in the present study population, 1 SNP in LCAT (rs4986970) and 1 SNP in SR-B1 (rs5891) were excluded from further analysis.

Samples and Determination of Serum Lipid Levels

Venous blood samples were drawn during fasting. The mean duration of fasting was 9.8 hours and blood samples were drawn with the subjects in sitting position. The samples were separated into serum, plasma, and buffy coat, and were stored at −80°C on the day of sampling. Serum lipid levels were examined as part of the health checkup or for research purposes at the institutions affiliated with this study [25]. LDL-C values were estimated using the Friedewald formula when triglyceride levels were less than 400 mg/dL.

Environmental Exposure

Information on environmental factors was obtained using a structural questionnaire [22]. The duration and intensity of daily activity (hard work and walking) and the frequency and intensity of habitual exercise were used to estimate the metabolic equivalents (METs). Alcohol consumption was estimated for 8 kinds of alcoholic beverages on the basis of intake frequency and ingested amount. Smoking habits were recorded as smoking status, number of cigarettes smoked, age at which he/she started smoking, and age at which he/she quit smoking. Dietary habits of the study participants were categorized into 8 groups on the basis of intake frequencies for each food item [26], [27].

Statistical Analysis

Geographical differences in mean age and HDL-C values were confirmed by using one-way ANOVA between geographical regions. The characteristics of the clinical lipid profiles and the lifestyles of the study subjects were compared using a trend test and Chi-square test in 3 groups according to HDL-C levels, with almost an equal number of subjects in men and women.

The score variable (1–3) was used for all SNP analyses in the general linear model, because each SNP seemed to show a different effect on HDL-C levels depending on whether the additive or genotypic model was used. The coefficients of the genetic and environmental factors for the HDL-C value (the dependent variable) were estimated in the 2 models after adjusting for age (continuous), sex, smoking, drinking, daily activity (continuous METs), habitual exercise (continuous METs), intake of egg, meat, and fish, body mass index (BMI) and regions (site 1 to site 6). The Kyoto and Tokushima regions were combined into a single site because there was no difference in the distribution of SNPs between them, and the number of subjects in both regions was relatively small. We included each of the 6 regions in the model for adjustment variables except for the Amami Islands region, which was used as a reference site. METs per day were estimated from the information given in the questionnaire by using time and an intensity of 3.0 for walking and 4.5 for heavy work. Food intake was also estimated from the questionnaire on the basis of intake frequency; zero intake was included as the lowest category (“almost never”). To avoid any collinearity effect caused by gene linkage, 2 SNPs in ABCA1 (rs2422493 and rs1800976), which were in close proximity, were not included as adjusted factors at the same time.

HDL-C levels were also compared according to each genotype by sex and SNP by using the trend test to assess whether the effect of each SNP was from the threshold or additive model.

Gene-gene combination effect and interaction was evaluated using a logistic regression model by using the dummy variable (0, 1) for SNPs to compare the magnitude of their odd ratios (ORs). From each gene, a representative SNP was selected that showed the lowest p-value and highest coefficient in the general linear model. We evaluated gene-gene combination effects and interactions in both recessive and dominant models. Genotypes homozygous and heterozygous for the major allele were used as the reference in the recessive model, and genotypes homozygous for the minor allele were used as the case. The dominant model comprised the homozygous genotype of the major allele as the reference and other two genotypes as the case. As we selected 4 SNPs, and each SNP was categorized into 2 groups (reference and case), the number of combination was 16 (2×2×2×2), including one combination with all references of 4 SNPs. To minimize multicollinearity, p-values for the interaction were calculated after each variable was centralized. A dose-dependent effect was assessed using the ORs according to the number of minor alleles as score variables (0–5). The number of alleles was counted as zero with a homozygous genotype for the major allele; one with a heterozygous genotype; and 2 with a homozygous genotype for the minor allele. The reference group was categorized as those with 0–3 minor alleles because of the small number of subjects. High HDL-C levels in the logistic regression model were defined as a level ≥62 mg/dL at the subject median in men and women. ORs and their 95% confidence intervals (CIs) were estimated after adjusting for age (continuous), sex, and regions (site 1 to site 6).

P-values less than 0.05 were considered statistically significant. We applied the Bonferroni correction in the analysis of SNPs to decrease the potential of an alpha error by multiple hypothesis testing. P-values less than 0.0038, 0.0033, or 0.010 were considered nominally significant in the regression analysis of the 13 SNPs; the 15 combinations for the selected 4 SNPs; and the reference and 5 categories for the number of alleles, respectively. Statistical calculations were performed using the Stata Version 10 software (Stata Corp, College Station, TX). Genotypes with distributions from the Hardy–Weinberg equilibrium (HWE) were assessed using the Chi-square test with the “genhwi” command in Stata.

Results

The mean age and HDL-C values varied by region in both men and women (Table 1); the mean age ranged between 48.6 and 61.1 years in men and 46.4 and 60.2 years in women. The range of HDL-C values was between 55.0 and 64.7 mg/dL in men and 63.7 and 75.0 mg/dL in women. P-values for both the mean age and HDL-C values, determined using one-way ANOVA between the regions, were statistically significant by sex. The same degree of geographical variation in HDL-C values was not always observed between men and women.

thumbnail
Table 1. Geographical distribution of mean age and HDL-cholesterol values in the study subjects according to sex.

https://doi.org/10.1371/journal.pone.0082046.t001

The mean age did not differ according to HDL-C levels in men and was slightly lower in women with high HDL-C (Table 2). The BMI averages and levels of LDL-C and triglycerides tended to be lower with high HDL-C in both men and women; however, the level of total cholesterol was positively correlated with the level of HDL-C. The level of METs in daily activity was positively correlated with the level of HDL-C in men. Proportionately more male smokers had low HDL-C, and more current drinkers, both men and women, had high HDL-C. Neither increasing nor decreasing trends were observed for the subjects who reported frequent intake of egg, meat, and fish.

thumbnail
Table 2. Characteristics of clinical lipid profiles and lifestyles in the study subjects according to sex.

https://doi.org/10.1371/journal.pone.0082046.t002

The distribution of the minor allele frequency (MAF) from the selected 13 SNPs in the 4 genes is listed in Table 3. The MAF was distributed from 0.068 of APOA1 rs12718465 to 0.495 of CETP rs5882. The genotype frequencies were all in HWE, with the exception of CETP (rs3764261; p<0.001). The geographical variation in the MAF of the CETP rs3764261 allele ranged from 0.173 to 0.337 depending on the region. Genotype frequencies in the Shizuoka and Okazaki regions were not in HWE, but the frequencies of the other 6 regions were in HWE. We also compared the MAF between the studied 3,050 subjects and 1,469 subjects who were excluded from analysis due to various reasons to evaluate the selection bias for excluded subjects. The genotype frequencies of all 13 SNPs among excluded 1,469 subjects were in HWE.

thumbnail
Table 3. Distribution of minor allele frequency in the presented 13 single nucleotide polymorphisms in 4 genes (n = 3050).

https://doi.org/10.1371/journal.pone.0082046.t003

We evaluated the effect of genetic and environmental factors on high HDL-C levels using the general linear models. The SNPs in ABCA1 rs2422493, rs1800976, and rs2740483, and CETP rs708272, rs3764261, and rs1800775 were related to high HDL-C levels in model 1 after adjusting for age and sex (Table 4). SNPs (rs708272 and rs3764261) in the CETP gene were still statistically significant in model 2 after the Bonferroni correction (p = 0.003 and p<0.001, respectively), and the coefficient value (2.93 mg/dL/allele) for the CETP rs3764261 polymorphism was the highest among all the SNPs. Current drinking and daily activity were also associated with high HDL-C levels in both model 1 and 2 after adjusting for age, sex, environmental factors, and SNPs. Current smoking was negatively correlated with HDL levels. Egg intake had a borderline association with HDL-C levels in model 2.

thumbnail
Table 4. Coefficients of genetic and environmental factors in the general linear model for HDL-cholesterol levels.

https://doi.org/10.1371/journal.pone.0082046.t004

To evaluate the gene-gene combination effect and interaction, we selected one SNP for each of the 4 genes (ABCA1 rs2740483, APOA1 rs12718465, CETP rs3764261, and SR-B1 rs5888) that showed relatively stronger associations, with the lowest p-value and the highest coefficient for each gene in Table 4. At first, HDL-C levels were compared according to each genotype by sex and SNP. An apparent additive trend was observed in the CETP rs3764261 polymorphism in both men and women, and those with the CETP TT genotype had the highest HDL-C levels among men and women, respectively (Figure 3). The other 3 SNPs seemed to have either threshold or additive effects, although the increasing trend was not statistically significant.

thumbnail
Figure 3. Comparison of HDL-C levels according to each genotype by sex and SNP.

P-value was tested by the trend test.

https://doi.org/10.1371/journal.pone.0082046.g003

The ORs for the high levels of HDL-C in combination with the selected 4 SNPs suggested a gene-gene combination effect in the recessive model considering the homozygote and heterozygote genotypes as the reference. Increased ORs were observed in every combination with CETP rs3764261, and the combination was statistically significant after the Bonferroni correction (Table 5). The combinations of 3 SNPs (APOA1 rs12718465, CETP rs3764261, and SR-B1 rs5888) and 4 SNPs (adding ABCA1 rs2740483) had higher ORs (5.06, 2.72–9.39; and 4.74, 2.01–11.18, respectively) than any other combinations. The OR values in these combinations seemed to be higher than those obtained after multiplying each OR value, but their 95% CIs had a wide range. The association in the combination of the SNPs other than with CETP was not apparent. Gene-gene interactions between the 4 SNPs were not observed. We also identified the ORs in the dominant model by using only the homozygote genotype as the reference, but the ORs in many combinations were omitted from the calculation because of small number of references available, especially for APOA1 (data not shown). We also performed the same analysis after excluding the regions not present in HWE, i.e., Shizuoka and Okazaki; excluding these regions from the analysis gave similar results as earlier (data not shown in Table). Increased ORs were also observed in every combination with CETP rs3764261, and the combination was statistically significant (p<0.0033). The combinations of 3 SNPs (APOA1 rs12718465, CETP rs3764261, and SR-B1 rs5888) and 4 SNPs (adding ABCA1 rs2740483) had higher ORs (4.94, 2.13–11.43; and 4.52, 1.43–14.31, respectively) than any other combinations, too. The association in the combination of the SNPs other than with CETP was not apparent, except the combinations of 2 SNPs (ABCA1 rs2740483 and SR-B1 rs5888; 2.89, 1.53–5.39).

thumbnail
Table 5. Odds ratios (ORs) and 95% confidence intervals (CIs) for high HDL-cholesterol with the combination of selected single nucleotide polymorphisms (SNPs) and their interactions.

https://doi.org/10.1371/journal.pone.0082046.t005

A dose-response effect on gene-gene interaction was also assessed. The ORs increased according to the number of alleles, and the dose-responsive trend was statistically significant (p<0.001; Table 6). We also conducted the same analysis after excluding the regions not present in HWE, i.e., Shizuoka and Okazaki. Excluding these regions from the analysis gave similar results as earlier, and the OR with 8 alleles was 6.01 (2.19–16.48) with significant dose-responsive trend (p<0.001; data not shown in Table).

thumbnail
Table 6. Odds ratios (ORs) and 95% confidence intervals (CIs) for high HDL-cholesterol with the number of alleles of selected single nucleotide polymorphisms.

https://doi.org/10.1371/journal.pone.0082046.t006

Discussion

The present study investigated the interaction between gene polymorphisms in ABCA1, APOA1, SR-B1, and CETP for serum HDL-C in a cross-sectional study. We found a multiplier effect with CETP rs3764261 and the 3 other SNPs, with an increasing trend according to the number of alleles for serum HDL-C. To our knowledge, this is the second study to observe a gene-gene interaction for HDL-C [21] and the first study of its kind among the Asian population.

Environmental Factors

The present study identified an association between environmental factors and HDL-C levels among the Japanese study population before investigating the gene-gene interactions. We observed a negative relationship between smoking and HDL-C levels. Previous studies were in concordance with this relationship and showed a potential influence on CETP activity [10], [28]. The relationship between alcohol consumption and increased HDL-C levels has been reported previously [6]. Alcohol consumption increases the expression of ABCA1 [29] and the concentration of APOA1 levels [6] in peripheral blood and decreases the CETP activity [30]. The present study also revealed a positive relationship between alcohol consumption and HDL-C levels. Several interventional studies have reported a positive association between habitual exercise and HDL-C levels [7]. Exercise increases HDL-C levels in a population that already has normal levels of HDL-C but causes no or limited increase in subjects with already low HDL-C levels [31], [32]. Exercise also decreases the volume of fatty tissue and influences HDL-C levels [32], [33]. We observed increased HDL-C levels with daily activity.

Several reports have revealed a positive association between HDL-C levels and the consumption of fish and vegetable oil containing high concentrations of n-3 polyunsaturated fatty acids [9], [34]. An animal study found that the consumption of corn oil results in increased transport rates of HDL protein, as well as increased binding of HDL to liver membranes [35]. However, the results analyzing the association between fish consumption and increased HDL-C levels are inconsistent [36]. The present study also revealed no association between fish consumption and increased HDL-C levels. The present questionnaire did not include information about the types of fish species that were consumed, although different fish species, such as fatty fish, may have a different effect on cholesterol levels.

Genetic Factors

Previous studies have shown that SNPs in the RCT system are responsible for varying levels of HDL-C [11][17]. The results of the present study also showed a significant relationship between the studied SNPs and serum HDL-C levels, except for LCAT, which had a lower MAF. The present results showed the highest coefficient for CETP rs3764261, which has also been observed in previous genome-wide screenings in the Japanese population and in the US EMR-linked biobanks [14], [21], [37]. A study that evaluated the Metabochip in African Americans reported an association between HDL-C and rs12740374 and rs17231520 in CETP [38]. We also observed a modest correlation with ABCA1 rs2740483, and APOA1 rs12718465. These SNPs are thought to play a partial role in influencing HDL-C levels because the combination with CETP rs3764261 shows a multiplier effect. The increasing trend with the number of alleles for serum HDL-C also supports the role of these SNPs, suggesting a dose-dependence on these alleles. The absence of an apparent interaction between these SNPs suggests that they may independently influence HDL-C levels. The previous study on gene-gene interaction using the US EMR-linked biobanks also adjusted for environmental factors such as BMI and smoking status [21]. The present study independently identified the effects of these SNPs after adjusting for a larger number of potential environmental factors in the general linear model.

The genotype frequency of CETP rs3764261 was not in HWE. Geographical differences in the allele frequency of the CETP rs3764261 polymorphism was observed in 2 regions (Okazaki and Shizuoka), and the genotype frequency was not in HWE in these regions. The deviation from HWE in these selected regions may potentially indicate population stratification rather than genotyping errors, because genotype frequencies of CETP polymorphisms in the 6 other regions were in HWE, and all SNPs were genotyped in the same laboratory at the same time. Okazaki and Shizuoka are located at the mainland of Japan, and their subjects were recruited from health checkup examinees of general population. Some selection bias may have influenced the genotype frequency of CETP polymorphisms, but its detail is unclear. We included the region in the analysis to control for regional CETP rs3764261 heterogeneity.

This study has several limitations, which need to be discussed. First, because this study was cross-sectional, the effects of diseases and medications on the results were not accounted for in our study. Therefore, we established strict exclusion criteria to remove these effects. Second, the lifestyle information obtained through the questionnaire includes a potential misclassification. To minimize this, we re-checked the questionnaire by employing trained interviewers. Third, geographic variation in SNPs in the Japanese population, including our study regions, has been previously reported [24], [39], [40]. Therefore, we included the region in question in our analysis to control for regional heterogeneity. Fourth, we did not analyze the fractions of HDL-C levels. Atheroprotective and non-atheroprotective HDL particles have different effects on coronary heart disease risk; therefore, fractional analysis of HDL-C levels is warranted in future studies. Fifth, bioinformaics on these genes and SNPs, and the function of SNPs are partially used to compute the allele dosage. We selected candidate SNPs of the genes in the RCT system by referring to previous reports that partially used bioinformatics and/or SNP functions, in order to compute allele dosage. Then, we selected one SNP from each gene in the RCT system after estimating their effects. Sixth, the potential of an alpha error needs to be considered, because the number of the present subjects was relatively large and multiple comparisons were made. Therefore, we applied Bonferroni correction. On the other hand, reducing an alpha error for null association increases a beta error for those associations that are not null. Several researchers mentioned that no adjustments are needed for multiple comparisons, although extended debate did not achieve consensus [41], [42]. Some results of smaller ORs in the present study become statistical significance, if Bonferroni correction was not applied, which did not influence the conclusion. Further biological relevance will be useful to support the present results.

In conclusion, the present study confirmed that smoking and drinking habits, daily activity, and polymorphisms in the CETP gene are associated with HDL-C levels. We found a multiplier effect from the CETP rs3764261 polymorphism in combination with ABCA1, APOA1, and SR-B1 polymorphisms and a dose-dependency according to the number of alleles, although no gene-gene interaction was apparent.

Acknowledgments

The authors thank Kyota Ashikawa, Tomomi Aoi and other members of the Laboratory for Genotyping Development, Center for Genomic Medicine, RIKEN.

Author Contributions

Conceived and designed the experiments: HN TT. Performed the experiments: MK. Analyzed the data: AN TT. Contributed reagents/materials/analysis tools: KK EM KW NH YN TCT SS KO HU EO SH HM. Wrote the paper: AN NH TT. Summarized the whole cohort: HT.

References

  1. 1. Mackay J, Mensah G (2004) Atlas of Heart Disease and Stroke. Geneva, World Health Organization.
  2. 2. Kubo M, Kiyohara Y, Kato I, Tanizaki Y, Arima H, et al. (2003) Trends in the incidence, mortality, and survival rate of cardiovascular disease in a Japanese community: the Hisayama study. Stroke 34: 2349–2354.
  3. 3. Wilson PW, D'Agostino RB, Levy D, Belanger AM, Silbershatz H, et al. (1998) Prediction of coronary heart disease using risk factor categories. Circulation 97: 1837–1847.
  4. 4. D'Agostino RB, Russell MW, Huse DM, Ellison RC, Silbershatz H, et al. (2000) Primary and subsequent coronary risk appraisal: new results from the Framingham study. Am Heart J 139: 272–281.
  5. 5. Voight BF, Peloso GM, Orho-Melander M, Frikke-Schmidt R, Barbalic M, et al. (2012) Plasma HDL cholesterol and risk of myocardial infarction: a mendelian randomisation study. Lancet 380: 572–580.
  6. 6. Rimm EB, Williams P, Fosher K, Criqui M, Stampfer MJ (1999) Moderate alcohol intake and lower risk of coronary heart disease: meta-analysis of effects on lipids and haemostatic factors. BMJ 319: 1523–1528.
  7. 7. Cauley JA, Kriska AM, LaPorte RE, Sandler RB, Pambianco G (1987) A two year randomized exercise trial in older women: effects on HDL-cholesterol. Atherosclerosis 66: 247–258.
  8. 8. Weggemans RM, Zock PL, Katan MB (2001) Dietary cholesterol from eggs increases the ratio of total cholesterol to high-density lipoprotein cholesterol in humans: a meta-analysis. Am J Clin Nutr 73: 885–891.
  9. 9. Nakamura T, Takebe K, Tando Y, Arai Y, Yamada N, et al. (1995) Serum fatty acid composition in normal Japanese and its relationship with dietary fish and vegetable oil contents and blood lipid levels. Ann Nutr Metab 39: 261–270.
  10. 10. Criqui MH, Wallace RB, Heiss G, Mishkel M, Schonfeld G, et al. (1980) Cigarette smoking and plasma high-density lipoprotein cholesterol. The Lipid Research Clinics Program Prevalence Study. Circulation 62: IV70–76.
  11. 11. Jensen MK, Pai JK, Mukamal KJ, Overvad K, Rimm EB (2007) Common genetic variation in the ATP-binding cassette transporter A1, plasma lipids, and risk of coronary heart disease. Atherosclerosis 195: e172–e180.
  12. 12. Rejeb J, Omezzine A, Rebhi L, Naffeti I, Kchok K, et al. (2008) Association of the cholesteryl ester transfer protein Taq1 B2B2 genotype with higher high-density lipoprotein-cholesterol concentrations and lower risk of coronary artery disease in a Tunisian population. Archives of Cardiovascular Diseases 101: 629–636.
  13. 13. van Acker BAC, Botma GJ, Zwinderman AH, Kuivenhoven JA, Dallinga T, et al. (2008) High HDL cholesterol does not protect against coronary artery disease when associated with combined cholesteryl ester transfer protein and hepatic lipase gene variants. Atherosclerosis 200: 161–167.
  14. 14. Hiura Y, Shen CS, Kokubo Y, Okamura T, Morisaki T, et al. (2009) Identification of genetic markers associated with high-density lipoprotein-cholesterol by genome-wide screening in a Japanese population: the Suita study. Circ J 73: 1119–1126.
  15. 15. Roberts CG, Shen H, Mitchell BD, Damcott CM, Shuldiner AR, et al. (2007) Variants in scavenger receptor class B type I gene are associated with HDL cholesterol levels in younger women. Hum Hered 64: 107–113.
  16. 16. Boes E, Coassin S, Kollerits B, Heid IM, Kronenberg F (2009) Genetic-epidemiological evidence on genes associated with HDL cholesterol levels: A systematic in-depth review. Experimental Gerontology 44: 136–160.
  17. 17. Brown CM, Rea TJ, Hamon SC, Hixson JE, Boerwinkle E, et al. (2006) The contribution of individual and pairwise combinations of SNPs in the APOA1 and APOC3 genes to interindividual HDL-C variability. J Mol Med (Berl) 84: 561–572.
  18. 18. Junyent M, Tucker KL, Smith CE, Garcia Rios A, Mattei J, et al. (2009) The effects of ABCG5/G8 polymorphisms on plasma HDL cholesterol concentrations depend on smoking habit in the Boston Puerto Rican Health Study. J Lipid Res 50: 565–573.
  19. 19. Marques-Vidal P, Bochud M, Paccaud F, Waterworth D, Bergmann S, et al. (2010) No interaction between alcohol consumption and HDL-related genes on HDL cholesterol levels. Atherosclerosis 211: 551–557.
  20. 20. Zhou Y, Yin R, Deng Y, Li Y, Wu J (2008) Interactions between alcohol intake and the polymorphism of rs708272 on serum high-density lipoprotein cholesterol levels in the Guangxi Hei Yi Zhuang population. Alcohol 42: 583–591.
  21. 21. Turner SD, Berg RL, Linneman JG, Peissig PL, Crawford DC, et al. (2011) Knowledge-driven multi-locus analysis reveals gene-gene interactions influencing HDL cholesterol level in two independent EMR-linked biobanks. PLoS One 6: e19586.
  22. 22. Tsompanidi EM, Brinkmeier MS, Fotiadou EH, Giakoumi SM, Kypreos KE (2010) HDL biogenesis and functions: Role of HDL quality and quantity in atherosclerosis. Atherosclerosis 208: 3–9.
  23. 23. Hamajima N (2007) Group JMS (2007) The Japan Multi-Institutional Collaborative Cohort Study (J-MICC Study) to detect gene-environment interactions for cancer. Asian Pac J Cancer Prev 8: 317–323.
  24. 24. Wakai K, Hamajima N, Okada R, Naito M, Morita E, et al. (2011) Profile of Participants and Genotype Distributions of 108 Polymorphisms in a Cross-sectional Study of Associations of Genotypes With Lifestyle and Clinical Factors: A Project in the Japan Multi-institutional Collaborative Cohort (J-MICC) Study. J Epidemiol 21: 223–235.
  25. 25. Naito M, Eguchi H, Okada R, Ishida A, Nishio K, et al. (2008) Controls for monitoring the deterioration of stored blood samples in the Japan Multi-Institutional Collaborative Cohort Study (J-MICC Study). Nagoya J Med Sci 70: 107–115.
  26. 26. Tokudome S, Goto C, Imaeda N, Tokudome Y, Ikeda M, et al. (2004) Development of a data-based short food frequency questionnaire for assessing nutrient intake by middle-aged Japanese. Asian Pac J Cancer Prev 5 (1) 40–43.
  27. 27. Tokudome Y, Goto C, Imaeda N, Hasegawa T, Kato R, et al. (2005) Relative validity of a short food frequency questionnaire for assessing nutrient intake versusthree-day weighed diet records in middle-aged Japanese. J Epidemiol 15 (4) 135–145.
  28. 28. Mero N, Van Tol A, Scheek LM, Van Gent T, Labeur C, et al. (1998) Decreased postprandial high density lipoprotein cholesterol and apolipoproteins A-I and E in normolipidemic smoking men: relations with lipid transfer proteins and LCAT activities. J Lipid Res 39: 1493–1502.
  29. 29. Hoang A, Tefft C, Duffy SJ, Formosa M, Henstridge DC, et al. (2008) ABCA1 expression in humans is associated with physical activity and alcohol consumption. Atherosclerosis 197: 197–203.
  30. 30. Hannuksela M, Marcel YL, Kesaniemi YA, Savolainen MJ (1992) Reduction in the concentration and activity of plasma cholesteryl ester transfer protein by alcohol. J Lipid Res 33: 737–744.
  31. 31. Williams PT, Stefanick ML, Vranizan KM, Wood PD (1994) The effects of weight loss by exercise or by dieting on plasma high-density lipoprotein (HDL) levels in men with low, intermediate, and normal-to-high HDL at baseline. Metabolism 43: 917–924.
  32. 32. Zmuda JM, Yurgalevitch SM, Flynn MM, Bausserman LL, Saratelli A, et al. (1998) Exercise training has little effect on HDL levels and metabolism in men with initially low HDL cholesterol. Atherosclerosis 137: 215–221.
  33. 33. Fonong T, Toth MJ, Ades PA, Katzel LI, Calles-Escandon J, et al. (1996) Relationship between physical activity and HDL-cholesterol in healthy older men and women: a cross-sectional and exercise intervention study. Atherosclerosis 127: 177–183.
  34. 34. Smith KM, Barraj LM, Kantor M, Sahyoun NR (2009) Relationship between fish intake, n-3 fatty acids, mercury and risk markers of CHD (National Health and Nutrition Examination Survey 1999–2002). Public Health Nutr 12: 1261–1269.
  35. 35. van Tol A, Terpstra AHM, van den Berg P, Beynen AC (1999) Dietary corn oil versus olive oil enhances HDL protein turnover and lowers HDL cholesterol levels in hamsters. Atherosclerosis 147: 87–94.
  36. 36. Eslick GD, Howe PR, Smith C, Priest R, Bensoussan A (2009) Benefits of fish oil supplementation in hyperlipidemia: a systematic review and meta-analysis. Int J Cardiol 136: 4–16.
  37. 37. Takeuchi F, Isono M, Katsuya T, Yokota M, Yamamoto K, et al. (2012) Association of genetic variants influencing lipid levels with coronary artery disease in Japanese individuals. PLoS One 7: e46385.
  38. 38. Buyske S, Wu Y, Carty CL, Cheng I, Assimes TL, et al. (2012) Evaluation of the metabochip genotyping array in African Americans and implications for fine mapping of GWAS-identified loci: the PAGE study. PLoS One 7: e35651.
  39. 39. Yamaguchi-Kabata Y, Nakazono K, Takahashi A, Saito S, Hosono N, et al. (2008) Japanese population structure, based on SNP genotypes from 7003 individuals compared to other ethnic groups: effects on population-based association studies. Am J Hum Genet 83: 445–456.
  40. 40. Nishiyama T, Kishino H, Suzuki S, Ando R, Niimura H, et al. (2012) Detailed analysis of Japanese population substructure with a focus on the southwest islands of Japan. PLoS One 7: e35000.
  41. 41. Rhthman KJ (1990) No adjustments are needed for multiple comparisons. Epidemiology 1: 43–46.
  42. 42. Greenland S (2008) Multiple comparisons and association selection in general epidemiology. Int J Epidemiol 37: 430–434.