Skip to main content

Genome- and epigenome-wide association study of hypertriglyceridemic waist in Mexican American families

Abstract

Background

There is growing interest in the hypertriglyceridemic waist (HTGW) phenotype, defined as high waist circumference (≥95 cm in males and ≥80 cm in females) combined with high serum triglyceride concentration (≥2.0 mmol/L in males and ≥1.5 mmol/L in females) as a marker of type 2 diabetes (T2D) and cardiovascular disease. However, the prevalence of this phenotype in high-risk populations, its association with T2D, and the genetic or epigenetic influences on HTGW are not well explored. Using data from large, extended families of Mexican Americans (a high-risk minority population in the USA) we aimed to: (1) estimate the prevalence of this phenotype, (2) test its association with T2D and related traits, and (3) dissect out the genetic and epigenetic associations with this phenotype using genome-wide and epigenome-wide studies, respectively.

Results

Data for this study was from 850 Mexican American participants (representing 39 families) recruited under the ongoing San Antonio Family Heart Study, 26 % of these individuals had HTGW. This phenotype was significantly heritable (h 2 r = 0.52, p = 1.1 × 10−5) and independently associated with T2D as well as fasting glucose levels and insulin resistance. We conducted genome-wide association analyses using 759,809 single nucleotide polymorphisms (SNPs) and epigenome-wide association analyses using 457,331 CpG sites. There was no evidence of any SNP associated with HTGW at the genome-wide level but two CpG sites (cg00574958 and cg17058475) in CPT1A and one CpG site (cg06500161) in ABCG1 were significantly associated with HTGW and remained significant after adjusting for the closely related components of metabolic syndrome. CPT1A holds a cardinal position in the metabolism of long-chain fatty acids while ABCG1 plays a role in triglyceride metabolism.

Conclusions

Our results reemphasize the value of HTGW as a marker of T2D. This phenotype shows association with DNA methylation within CPT1A and ABCG1, genes involved in fatty acid and triglyceride metabolism. Our results underscore the importance of epigenetics in a clinically informative phenotype.

Background

As the global epidemic of type 2 diabetes (T2D) continues to expand, better, accurate, and cost-efficient ways to stratify risk of individuals are urgently required. In that vein, it is now fairly well established that hypertriglyceridemic waist (HTGW) is a promising and simple measurement that is clinically feasible, yet accurate for risk stratification of T2D [13]. HTGW is a critical link between obesity (specifically, central obesity) and T2D. The rationale for using a combination of waist circumference and fasting triglyceride levels is to partially overcome the inability of waist circumference to discriminate between subcutaneous and visceral adiposity; raised serum triglycerides in the presence of increased waist circumference are more likely to indicate visceral rather than subcutaneous fat accumulation [4]. The last few years have seen a dramatic increase in studies around the world that have reported the use and importance of HTGW in prediction of T2D and cardiovascular risk [512].

Considering the strong genetic basis of T2D, waist circumference, and serum triglyceride levels, it is likely that there is also a strong genetic contribution to HTGW. Several studies have shown the possible involvement of HNF1A variants, epistatic interactions with genes involved in the very low density lipoprotein pathway and variation within ADIPOQ to be associated with HTGW [1316]. Generally, however, replicative literature support for these observations from large, population based studies has been lacking. Moreover, given the possibility of a significant environmental influence (e.g., through diet and physical activity) on both waist circumference and triglyceride levels [11, 17, 18], it is conceivable that the HTGW phenotype may be epigenetically modifiable. Although recent elegant studies have shown promising results for the association of DNA methylation separately with waist circumference [1921] and triglyceride levels [2224], to our knowledge, currently, there are no studies on the possible epigenetic role in the inter-individual variability in HTGW.

In this study, we investigated both the genetic and epigenetic basis of HTGW in Mexican Americans, a minority population in the USA who are at high risk for both obesity and T2D [2527]. To enhance our power to detect genetic and epigenetic effects, we conducted this study in pedigreed individuals representing large, extended families, recruited in the San Antonio Family Heart Study [28, 29]. Using genome- and epigenome-wide approaches, we show here an association between HTGW and T2D in Mexican Americans and identify three CpG sites in CPT1A and ABCG1 that are significantly and independently associated with the HTGW phenotype.

Results

Study participants

We studied 850 pedigreed Mexican Americans (representing 39 extended families), who were predominantly middle aged (mean age 46.75 years, standard deviation 14.54 years) and 63 % female. At the time of assessment 21 % of individuals had T2D (fasting glucose ≥7 mmol/L), another 17 % had impaired fasting glucose (fasting serum glucose between 5.55 mmol/L and 7 mmol/L), 56 % were obese (body mass index ≥30 kg/m2), and 32 % had hypertension (systolic blood pressure >140 mmHg and/or diastolic blood pressure >90 mmHg, or already receiving antihypertensive treatment). A total of 14, 16, and 24 % of participants were receiving lipid-lowering, anti-diabetic and anti-hypertension medications, respectively. The prevalence of central obesity (waist circumference ≥90 cm for males or ≥85 cm for females) was very high (88 %) in our sample but the prevalence of hypertriglyceridemia (≥2.0 mmol/L for males and ≥1.5 mmol/L in females) was comparatively lower (28 %). There were a total of 223 (26 %) participants with HTGW, with a higher prevalence in females (31 %) compared to males (18 %). To ensure that the potential association of genetic and epigenetic variants with HTGW is not confounded by the presence of comorbidities, we used robust statistical models that accounted for the simultaneous presence of comorbidities (systolic and diastolic blood pressure and presence of diabetes and obesity) in our analyses.

There was a significant phenotypic correlation between waist circumference and triglycerides in the study participants (ρ P = 0.19, p = 6.8 × 10−7). However, results obtained using bivariate trait analyses showed that this phenotypic correlation was primarily because of shared environmental influences (ρ E = 0.32, p = 0.0011) rather than due to shared genetic influences (ρ G = 0.07, p = 0.5898). We therefore conducted an investigation into the potential genetic and epigenetic basis of HTGW.

Association of HTGW with T2D

We first tested the hypothesis that HTGW is an independent determinant of T2D in our sample of Mexican American families. For this, we tested the association of HTGW with four T2D-related phenotypic traits (T2D, fasting blood glucose (FBG), fasting plasma insulin (FPI), and homeostatic model of assessment-insulin resistance (HOMA-IR)). After adjusting for clinically relevant confounders, we found that HTGW was significantly and independently associated with all four traits related to T2D (Table 1). The polygenic regression coefficients, when transformed to odds ratios, indicate that study participants with HTGW have 3.16 times higher odds (as compared to those without HTGW) of T2D (95 % CI 2.00–4.93).

Table 1 Independent association of HTGW with T2D-related traits in San Antonio Family Heart Study participants

We next evaluated whether the presence of both high waist circumference and hypertriglyceridemia is more strongly associated with T2D-related traits than either of these factors alone. For this, we used interactive polygenic regression models, the results of which are shown in Table 1. For all four traits studied, the regression coefficients for the combination of high waist and hypertriglyceridemia were statistically significantly associated with the trait in question. With the exception of T2D, the regression coefficients for the combination of high waist circumference and hypertriglyceridemia were substantially larger as compared to those for high waist circumference or hypertriglyceridemia alone. These results indicate that in the context of T2D-related phenotypes, HTGW is likely to be strongly associated with the underlying pathology of T2D and may therefore be an important T2D risk assessment tool.

Heritability of HTGW

We next estimated the heritability of HTGW using polygenic regression that accounted for familial relationships. Even after accounting for age and sex (and their first and second order interactions), use of lipid lowering drugs, use of anti-hypertensive medications, systolic blood pressure, diastolic blood pressure, and presence of diabetes and obesity, we found that the heritability of HTGW was 0.52 (SE 0.12, p = 1.1 × 10−5).

Genome-wide association analysis of HTGW

Results of the association between each of the 728,175 single nucleotide polymorphisms (SNPs) and liability of HTGW are shown in Fig. 1a, with details provided in Additional file 1: Table S1. The Q-Q plot for these analyses is shown in Additional file 2: Figure S1. There was no evidence of a systematic inflation of significance values (genomic inflation factor = 1.0251). After correcting for clinically relevant confounders as well as the top four principal components that reflected underlying population admixture, we did not find any statistically significant associations between genetic variants and liability of HTGW after correcting for multiple comparisons by controlling for a false discovery rate (FDR) of 5 %. We observed that even though clusters of SNP-HTGW association were visible on chromosomes 2 (intergenic), 4 (flanking the 3’ untranslated region of the LOC391698 locus), 13 (flanking the 5’ untranslated region of RFC3), and 20 (flanking the 3’ untranslated region of the LOC441940 locus), none of these associations reached genome-wide significance.

Fig. 1
figure 1

Manhattan plots showing the genome-wide association of DNA sequence variants and epigenome-wide association of DNA methylation with liability of HTGW. a Genome-wide association study. Genome-wide significance is indicated by the red horizontal line and chromosomal locations are color-coded. Results are from polygenic regression models that accounted for age, age2, sex, age × sex interaction, and age2 × sex interaction, the top four principal components quantifying ancestry-based population admixture, and use of anti-lipid, anti-hypertensive and anti-diabetic medications. b Epigenome-wide association study. Epigenome-wide significance is indicated by the red horizontal line. Significantly associated CpG sites are labeled and chromosomal locations are color-coded. Results are from polygenic regression models adjusted for age, age2, sex, age × sex interaction, and age2 × sex interaction, Illumina Sentrix® ID and Sentrix® position (to account for batch effects), estimated cell counts, and use of anti-lipid, anti-hypertensive and anti-diabetic medications. Q-Q plots corresponding to Panel A and B are shown in Additional file 2: Figures S1 and S2, respectively

Epigenome-wide analysis of HTGW

Of the 485,577 interrogated loci on the array, association analyses using polygenic regression models could be run without convergence failure on a total of 458,716 (94.5 %) autosomal CpG sites. Of these, 1385 (0.30 %) probes had detection p values >0.01 in >5 % of the study participants. We therefore excluded these probes from further analyses. Thus, our study included analysis on a total of 457,331 CpG sites. Full results of the association analyses between DNA methylation at each CpG site and liability of HTGW are shown graphically in Fig. 1b and provided in detail in Additional file 3: Table S2. The genomic inflation factor for these analyses was 0.6152 indicating a lack of genomic inflation of significance values (Q-Q plot for these analyses is shown in Additional file 2: Figure S2). After correcting for the clinical confounders as well as for multiple testing (FDR < 0.05), we found three CpG sites (cg00574958, cg17058475, and cg06500161, Fig. 1b and Table 2) to be significantly associated with HTGW. The top two sites were in the 5’ UTR of CPT1A while the third site was in ABCG1. It is also worth mentioning that the CpG site cg19693031 in the 3’ UTR of TXNIP was only marginally non-significant (FDR-corrected p = 0.0563, Additional file 3: Table S2). We previously found no evidence for the influence of genetic variation (±50 kb, assessed by SNP microarrays) on any of the CpG sites whose methylation levels were associated with HTGW [30]. Similarly, we had previously used whole-genome sequencing of a subset of our cohort (n = 197) to test for associations between probe-based SNPs and DNA methylation levels [30]. In this cohort subset, the mutant allele for rs78442314 within the CPT1A probe for cg00574958 was not present, and while the mutant allele for rs9982016 within the ABCG1 probe for cg06500161 was present in 2.79 % of the cohort subset, it was not significantly associated with DNA methylation levels detected by the probe (p = 0.1895). Together, this suggests that the significant associations seen between DNA methylation levels and HTGW are not driven by known sequence variation.

Table 2 Characteristics of the CpG sites significantly associated with HTGW

Notably, the median β values (Table 2) at both CpG sites near CPT1A were only slightly smaller in participants with HTGW as compared to those without HTGW (1.0 % for cg00574958 and 1.1 % for cg17058475). Conversely, the median methylation levels were 1.6 % higher in individuals with HTGW as compared to those without for the CpG site near the ABCG1. These results indicate that seemingly small changes in DNA methylation around the CPT1A and ABCG1 loci are associated with different clinical risk profiles. To determine whether the observed associations were driven by a small subset of individuals with extreme values of β, we studied their distribution in the study participants. We observed that there was no marked skew in the distribution of methylation at any of the three significantly associated sites (Fig. 2a, c, and e). A comparison of the methylation values categorized on the basis of presence or absence of HTGW also suggested that the observed associations were unlikely to be a result of skewed observations (Fig. 2b, d. and f). It was noteworthy that the DNA methylation levels at all of these sites were significantly heritable (Table 2; based on data from [30]). We also observed that none of the three probes were previously reported to be cross-reactive by Chen et al. [31].

Fig. 2
figure 2

Distribution of DNA methylation scores for the significantly associated CpG sites in the study participants. Panels a, c, and e show histograms based on all study participants while panels b, d, and f show box plots for the corresponding CpG site in those with (red boxes) or without (yellow boxes) HTGW

To determine the specificity of association of the three CpG sites with HTGW, we further adjusted for presence of type 2 diabetes, obesity, and systolic and diastolic blood pressures. We found that there was a negligible loss in the strength of association of the three sites with HTGW after accounting for the confounding comorbidities (Table 2). Using an interactive model, we tested for association between HTGW, high waist circumference (with normal triglyceride levels), or high triglyceride levels (with normal waist circumference), and DNA methylation at each of the sites. For each CpG site, we found the strongest association to be with the HTGW phenotype. We also examined whether the strategy to include medication use as a covariate in the regression models was confounding the associations. For this, we ran two sets of models for each of these sites, one ignoring concurrent medication use and the second restricting the analyses to only those individuals who were not receiving any medication. The results of these analyses (see Additional file 4: Table S3) showed that all associations remained significant irrespective of the strategy used to account for medication use. We had previously found these sites to be associated with T2D [30], and next considered whether this is mediated through their association with HTGW. For this, we used Sobel’s method of estimating mediation [32] and found that the mediation parameter was statistically significant for each of these three sites. Together, these results show that not only are these three CpG sites strongly and independently associated with HTGW, but it is likely that the association of these sites with T2D is partly but significantly explained through HTGW.

We further considered whether the three significantly associated CpG sites contain redundant statistical information. For this, we ran a series of four models (Table 3), in which each significant CpG site was added to the polygenic regression model in a forward stepwise manner. For each step, we estimated the Kullback–Leibler R 2 (K-L R 2) as a measure of information content. We observed that the addition of each CpG site was associated with a statistically significant improvement in the K-L R 2 statistic. The final model that included the clinical covariates along with all three CpG sites accounted for 19.67 % of the inter-individual variability in HTGW (Table 3). Considering that the base model shown in Table 3 accounted for 10.15 % of HTGW variability, our results indicate that the top three significantly associated CpG sites accounted for an additional 9.52 % of the inter-individual variability in HTGW. Together, these results show that the three sites were significantly, specifically, and independently associated with the HTGW phenotype.

Table 3 Improvement in Kullback-Leibler R 2 by including CpG sites significantly associated with HTGWa

Validation of HumanMethylation450 BeadChip array by pyrosequencing

We validated the measurements made by the HumanMethylation450 BeadChip array by pyrosequencing. We have previously demonstrated that for the cg06500161 CpG site within ABCG1, there was a significant concordance between the results obtained by HumanMethylation450 BeadChip array and by pyrosequencing (Spearman’s ρ = 0.42, p < 0.0001) [30]. Here, we demonstrate that a significant correlation also exists between the results of these two techniques for the most significantly associated CpG site in CPT1A, cg00574958. We observed a statistically significant correlation between the levels of methylation measured by the microarray and by pyrosequencing (Spearman’s ρ = 0.47, p < 0.0001; Fig. 3a). In general, the microarray reported inflated levels of DNA methylation compared to that measured by pyrosequencing. Therefore, we formally tested the agreement between these two methods by using Bland-Altman plot analysis (Fig. 3b). We found that on average, pyrosequencing reported methylation levels about 6.2 % less than those of the microarray. However, the limits of agreement between these two methods ranged from −0.102 to −0.022 and only 27 (~3.2 %) samples lay outside these limits (Fig. 3b, dots colored red). Pitman’s test of unequal variances was also not significant (p = 0.144). Thus, our results indicate that there was a moderate correlation and a very good agreement between the results obtained by microarray and pyrosequencing techniques. Consistently, we also noticed that the median DNA methylation measured by pyrosequencing was lower in individuals with HTGW as compared to those without HTGW (Fig. 3c). Further, when we ran polygenic regression models and adjusted for the same covariates as mentioned in Model 1, Table 2, we found that the regression coefficient for pyrosequencing-based methylation at cg00574958 was −0.28 (p = 4.4 × 10−7). Thus, our results using the HumanMethylation450 BeadChip array were reproducible and in agreement with those derived using pyrosequencing.

Fig. 3
figure 3

Agreement between the results of HumanMethylation450 BeadChip array and pyrosequencing for the cg00574958 CpG site. a Correlation scatter plot. b Bland-Altman plot of the difference in the measurements of the two methods (ordinate) versus the mean (abscissa). Limits of agreement (LAG) are shown pictorially using dashed horizontal lines. Samples that fall outside the LAG are colored red. c Distribution of the DNA methylation as measured by pyrosequencing in individuals with (red box) and without (yellow box) HTGW

Discussion

Our study reports several novel findings: (1) the prevalence of HTGW in our sample of Mexican Americans is ~26 %; (2) HTGW is a significant independent predictor of T2D as well as fasting glucose, insulin levels and insulin resistance; (3) HTGW is significantly heritable; (4) DNA methylation within the 5’UTR region of CPT1A and in the gene body of ABCG1 is significantly associated with HTGW, independent of related comorbidities; (5) methylation at the top three significantly associated CpG sites (in CPT1A and ABCG1) together account for 9.52 % of the variability of HTGW, which is unlikely to be due to surrounding sequence variants. We also observed that seemingly small changes in DNA methylation levels (<2 % at all three identified CpG sites) are associated with HTGW. This finding implies that DNA methylation levels may be regulated with considerable plasticity and that small departures from such plasticity may have implications in disease pathogenesis through altered gene expression. It is important to consider these observations in the light of existing literature, study limitations, and clinical implications.

Using a previously established definition of HTGW [5, 3336], our estimate of the prevalence of HTGW is within the range of 12.1–35 % prevalence previously reported by various studies across the world [6, 8, 3741], and is comparable to that reported from a nationally representative US sample [42]. However, within our cohort, the prevalence of high waist circumference exceeded the prevalence of hypertriglyceridemia, which raises the possibility that in our study population, there was higher likelihood of accumulation of subcutaneous rather than visceral adipose tissue. Interestingly, such a finding has been confirmed with computerized tomography studies in Mexican American women [43]. This hypothesis needs to be tested and confirmed in future studies.

HTGW remains a largely understudied phenotype and neither its genetic nor epigenetic basis is known. Our study provides two initial observations in this regard. First, it reconfirms the association between HTGW and T2D-related traits [3, 7, 8, 41, 44, 45], which has previously been observed in other populations, and we now replicate in Mexican American families. Second, estimation of heritability is often considered an essential first step towards understanding the genetic basis of a phenotypic trait, and we find that HTGW is indeed significantly heritable (h 2 r = 0.52). These results advocate a need for the HTGW phenotype to be intensively investigated from a genetic perspective. It should be noted that since HTGW captures two components of the metabolic syndrome (central obesity and dyslipidemia), this phenotype is likely to be associated with coexisting conditions like hypertension, hyperglycemia, insulin resistance, and type 2 diabetes in epidemiological studies. In this study, we used robust statistical models to account for the presence of these comorbidities, and we suggest that future studies on the HTGW phenotype also be carefully designed so as to understand the true associations, not confounded by comorbidities.

We did not observe any genome-wide statistically significant associations between SNPs and HTGW. To our knowledge, this is the first genome-wide association study examining HTGW. Genome-wide association analyses have identified at least 36 loci associated with triglyceride levels, including prominent replicable associations in well implicated lipid genes LPL, APOA1, LIPC, and CETP [4649]. Conversely, only a few loci have been implicated in waist circumference via genome-wide association studies, namely MCR4, TFAP2B, MSRA, NRXN3, and MAP3K1, although in most cases, these loci have not been replicated [5054]. Our genome-wide association study of HTGW did not establish any of the previously identified loci for either triglyceride levels or waist circumference, and our most significant loci were on chromosomes 2, 4, 11, 13, and 20. While these loci may represent novel determinants of HTGW, it should be remembered that given the small effect sizes that specific genotypes typically inflict upon complex phenotypes, it is likely that our cohort is underpowered to detect such effects.

Although we did not find any evidence for association between sequence variation and HTGW in our cohort, we identified three CpG sites whose DNA methylation levels were associated with HTGW, and these associations were largely independent of other metabolic syndrome phenotypes or DNA sequence variation. Our most striking observation is the highly significant association between the CPT1A-related CpG sites and HTGW, which was independent of related phenotypes (blood pressure, obesity and T2D). Carnitine palmitoyltransferase acts at the outer mitochondrial membrane to escort long-chain fatty acids into the mitochondria for β oxidation [55]. The CPT1A gene that regulates the expression of this enzyme is 60 kb long and contains 20 exons, and is located on chromosome 11 [56, 57]. It is regulated by PPARA, a nuclear transcription factor that plays a critical role in the oxidation of fatty acids [58, 59]. Considering this molecular nexus, the role of CPT1A in triglyceride metabolism is now well established. However, the role of this gene in waist girth is less obvious. Two population based studies [24, 60] have shown association of sequence variants in CPT1A with various indices of obesity including waist circumference, but the exact mechanism by which CPT1A might contribute to altered adiposity remains unknown. We cannot comment on the precise mechanisms involved, but this opens up a possible and interesting thread of research for future studies. Our results are in line with the emerging evidence that methylation in the promoter region of CPT1A and specifically at the cg00574958 site is associated with several aspects of triglyceride biochemistry including associations with lipoprotein sub-fractions, hypertriglyceridemia, and the effect of lipid lowering drugs [22, 23, 61].

Our results indicate that the other gene that might play a role in HTGW by way of altered methylation levels is the ATP binding cassette G1 (ABCG1) gene. There is increasing evidence to show that this gene plays an important role in triglyceride metabolism [6265]. For example, a recent epigenome-wide study [63] has found a strong association between the same CpG site (cg06500161) that we observed and plasma triglycerides. Although the exact mechanism for the implication of this gene in triglyceride metabolism is unclear, it has been posited that the gene modulates bioavailability of plasma lipoprotein lipase and thus induces lipid accumulation in a triglyceride-rich environment [66]. Moreover, another recent epigenome-wide association study identified methylation at ABCG1 to be significantly associated with waist circumference [67]. It is of interest that we previously observed a strong association of the cg06500161 CpG site with several type 2 diabetes-related traits in this Mexican American cohort [30], and others also found it to be associated with insulin resistance [68]. Thus, these findings raise the possibility that methylation at ABCG1 may provide critical insights into the reported [69] association of central obesity (HTGW) and type 2 diabetes.

Interestingly, we found a marginally significant association (at the level of the epigenome) between the CpG site cg19693031, in TXNIP, and HTGW. The TXNIP protein regulates intra- and extra-cellular reduction-oxidation cycles [70, 71] and genetic variations in this gene have previously been shown to be associated with hypertriglyceridemia in individuals with T2D [72]. Therefore, TXNIP also appears to be an attractive candidate involved in the pathogenesis of T2D mediated by HTGW. Future replications and mechanistic studies are required to definitively support the associations observed in this study.

It should be noted that two of the three significantly associated CpG sites contained a SNP within the probe sequence, and it is therefore possible that these SNPs may drive the methylation-HTGW association. However, in a subset of 197 Mexican Americans for whom deep sequencing data was available, we previously found that the SNP was either not present in our cohort subset (for the rs78442314/cg00574958-containing probe), or that it does not significantly influence DNA methylation levels (for the rs9982016/cg06500161-containing probe) [30]. Further, we had observed that association between a SNP and methylation score was much more likely if the SNP was at the CpG site rather than elsewhere within the probe. Together these observations suggest that the CpG sites significantly associated with HTGW in our cohort are unlikely to have been driven by the SNPs contained within the probes and highlights the need for careful consideration of probe-exclusion criteria to avoid the potential loss of important biological associations.

Some limitations of our study need to be considered before generalizing these results. First, our study cannot directly demonstrate the directionality of the methylation ➔ gene expression ➔ phenotype association. However, it has been shown that methylation at cg00574958 can influence expression of CPT1A and that methylation within ABCG1 also influences expression of the gene [22, 63]. Second, DNA methylation is regulated both genetically and environmentally. Since there is a strong environmental correlation between waist circumference and triglycerides, methylation at CPT1A and ABCG1 is associated with HTGW, and HTGW is a significant predictor of type 2 diabetes related traits, our study points towards possible environmental links with known modifiers of methylation such as diet, physical activity, and pharmacological interventions [73, 74]. These possibilities need to be explored further in the future, and may provide stronger justification for putative preventive interventions.

Conclusions

Prevalence of HTGW and its importance as a marker of T2D has not been clearly demonstrated in Mexican Americans—a high-risk, minority population in the USA. In addition to characterizing HTGW in our study cohort, our results also raise the possibility for an epigenetic basis of HTGW. Our epigenome-wide and genome-wide association studies show that HTGW may be mediated by epigenetic factors, and may also be influenced by the environment. We have identified two CpG sites (cg00574958 and cg17058475) in the 5’ UTR of CPT1A and one CpG site (cg06500161) in the body of ABCG1, which are associated with HTGW. These genes are involved in β oxidation of long-chain fatty acids and triglyceride storage, respectively. Our results highlight the role of epigenetics in HTGW, which is an important marker of T2D and cardiovascular disease.

Methods

Study participants

Participants in this study were from the San Antonio Family Heart Study [28, 29], an ongoing prospective evaluation of Mexican American families living in San Antonio. A total of 850 participants from 39 families were included in the analysis. The recruitment and ascertainment of participants and their phenotyping and genotyping have been extensively described elsewhere [29, 75]. This data was collected in the third wave of ascertainment (2002–2006). Peripheral blood samples were collected following an overnight fast and extensive anthropometric phenotyping was conducted.

Ethics, consent, and permissions

Written consent was obtained for all individuals in this study. This study was approved by the Institutional Review Board at The University of Texas Health Science Center at San Antonio.

Phenotypes

Our main phenotype of interest was HTGW, which does not presently have a uniform definition across studies [5, 6, 39, 44, 45, 76]. In our study, we defined HTGW as high waist circumference (≥90 cm in males and ≥85 cm in females) combined with high serum triglyceride concentration (≥2.0 mmol/L in males and ≥1.5 mmol/L in females) [5, 3336]. Other phenotypes included in this study were: age, sex, systolic and diastolic blood pressure, use of anti-lipid, anti-hypertensive and anti-diabetic medications, presence of type 2 diabetes (diagnosed using the ADA criteria [77], fasting glucose ≥7 mmol/L), triglycerides, waist circumference, and presence of obesity (body mass index ≥30 kg/m2). Methods used to measure these and other phenotypes have been described previously [29, 75].

Genotyping

Study participants were previously genotyped for approximately one million single nucleotide polymorphism markers using several Illumina genotyping arrays, including the HumanHap550v3, HumanExon510Sv1, Human1Mv1, and Human1M-Duov3. The Infinium Whole-Genome Genotyping Assay was employed according to manufacturers’ instructions. Details of the data cleaning and imputing steps for this genotypic data have been described previously [20]. Out of a total of 995,320 SNPS genotyped, we incorporated 759,809 SNPs into the association analyses, which had ≥97 % call rate, a minor allele frequency ≥5 % and a Hardy-Weinberg significance value ≥0.001.

DNA methylation assays and data preprocessing

DNA samples (500 ng) obtained from peripheral blood cells were bisulfite-converted using the EZ-96 DNA Methylation™ Kit (Zymo Research, Irvine, CA) and were subjected to methylation profiling using the Illumina Infinium HumanMethylation450 BeadChip assay (Illumina, San Diego, CA). The array interrogated 485,577 CpG sites across the genome and incorporated both Infinium I and Infinium II bead types. The Illumina iScan was used to scan the BeadChips and raw data was imported into GenomeStudio (V2011.1) to extract image intensities, following background subtraction and normalization to inbuilt controls on the arrays.

Methylation at each CpG site was quantified on a scale from zero (representing fully unmethylated) to one (representing fully methylated) as a methylation score (β). Probes that were located on the sex chromosomes (n = 11,648), that were non-CpG loci (n = 2,994) or that analyzed SNPs (n = 65) were excluded. To correct for differences due to Infinium I and Infinium II probe designs, we used the BMIQ method implemented in the R-based software, BMIQ [78].

We included only those probes for which heritability analyses could be successfully completed without convergence failures. For 12,154 (2.5 %) probes SOLAR was unable to achieve convergence, leaving a total of 458,716 CpG sites available for analysis. Of these, 1385 probes had detection p values >0.01 in >5 % of the samples and therefore we excluded these probes leaving a total of 457,331 CpG sites that were finally included in this study. To minimize loss of informative associations, we did not exclude SNP-containing or cross-reactive probes but rather investigated whether the significantly associated CpG sites could have been confounded due to these characteristics.

Pyrosequencing

For validation of the Illumina microarray data, we performed pyrosequencing on our most significant association (cg00574958 in the CPT1A). For each sample, 500 ng of genomic DNA was bisulfite converted, PCR-amplified, and subjected to pyrosequencing with the PyroMark96 MD (Qiagen, Valencia, CA). Percent DNA methylation was determined using PyroMark CpG software 1.0.11.14. The PCR was carried out at 95 °C for 5 min, followed by 40 cycles of 95 °C for 1 min, 56.6 °C for 1 min and 72 °C for 1 min, and a final extension at 72 °C for 7 min. Pyrosequencing was carried out according to the manufacturers’ instructions, and PCR and sequence primers (designed using the Pyromark Assay Design 2.0 software) were designated as

  • Forward primer: GTTTTTGGTATTGAGGTAAAATTAA

  • Reverse primer (biotinylated): AACCTTTCCAAATTCTTTAAAAC

  • Sequence primer: TTTTTGGTATTGAGGTAAAATTAAT

Statistical analysis

To ensure compatibility with the variance component framework and that the observed associations were unaffected by any undetected skew, we used an inverse normalization preprocessing step for the BMIQ-normalized methylation β values to circumvent any distributional aberrations. This step included ranking the raw values, generating cumulative density functions, and determining z-values based on the cumulative densities. All the transformed values were thus distributed as N(0,1) and could be represented on a common, comparable metric of z-values.

Family studies such as the present one have an added advantage that they can shed genetic light on the phenotypes under study. We first estimated the heritability (defined as the proportion of variability explained by genetic similarity among individuals) of HTGW using a polygenic regression model as follows:

$$ l(HTGW)=\mu +\boldsymbol{ba}+{g}_i+{e}_i $$

where l(HTGW) is the liability function of HTGW, μ is the overall mean liability, b is the regression coefficient vector corresponding to the covariate matrix a, g i is the polygenic effect, and e i is the measurement error. Heritability was then estimated as the ratio of variance due to genetic similarity (modeled as g i in the equation above) and the total phenotypic variability, Ω. The covariates used for the estimation of heritability of HTGW were: age, age2, sex, age × sex interaction, age2 × sex interaction, use of anti-lipid, anti-hypertensive and anti-diabetic medications, systolic blood pressure, diastolic blood pressure, and presence of type 2 diabetes and obesity (body mass index ≥30 kg/m2). Statistical significance of heritability (H o: heritability = 0) was tested by constraining the heritability to 0 and comparing the likelihood ratio statistics of the constrained and unconstrained models.

For genome-wide association analyses, we used polygenic regression models that accounted for age, age2, sex, age × sex interaction, and age2 × sex interaction, the top four principal components to quantify ancestry-based population admixture, and use of anti-lipid, anti-hypertensive, and anti-diabetic medications.

Our methylation studies used blood which contains a mixture of cell types. Reinius et al. [79] and Houseman et al. [80] have demonstrated the influence of differential cell proportions on DNA methylation signatures using different array platforms. Jaffe et al. [81] have extended this procedure to the Illumina Infinium HumanMethylation450 array. We estimated the proportion of CD4+ T cells, CD8+ T cells, B cells, natural killer cells, and granulocytes in each sample using the procedure described by Jaffe et al. [81] and adjusted all the polygenic regression models for these estimated cell counts as covariates. To test the association of DNA methylation at each CpG site with the liability function of HTGW, we ran polygenic regression models for each CpG site. In each model, we used age, age2, sex, age × sex interaction, and age2 × sex interaction, Illumina Sentrix® ID, Sentrix® position, estimated cell counts, and use of anti-lipid, anti-hypertensive, and anti-diabetic medications as covariates. Additional covariates were used for specific analyses and are described in the Results section. In particular, where warranted, we accounted for comorbidities (systolic blood pressure, diastolic blood pressure and presence of type 2 diabetes and obesity) that might influence the HTGW phenotype. Statistical significance for association was tested by constraining the regression coefficient to 0 and comparing the likelihood ratio statistics of the constrained and unconstrained model.

For heritability analyses as well as association analyses, we first estimated the genomic inflation factor (λ median) which was defined as the median χ 2 LL/invchi(0.5,1) where invchi(p,d) is the inverse χ 2 function for probability (p) and degrees of freedom (d), and corrected the nominal significance values for the estimated λ median. Additionally, we used the Benjamini-Hochberg procedure of false discovery rate (FDR) control for multiple testing correction. Significance was assessed at a global type I error rate of 0.05. The odds ratio (OR) for the association between HTGW and inverse-normalized methylation score was determined as \( {e}^{-\sqrt{\pi}\beta } \) where β represents the polygenic regression coefficient. The association of DNA methylation with type 2 diabetes that was mediated through HTGW was estimated using Sobel’s parameter [32]. For this, we ran two regression models for each site—the first model contained HTGW as the outcome and methylation (along with other covariates) as a predictor while the second model used T2D as the outcome and HTGW, methylation and other covariates as predictors. The regression coefficients from these two models were then multiplied to derive a quantified measure of mediation. Standard errors for this parameter were estimated as described by Sobel [32].

To parse out the genetic and environmental covariance, we used the methods entailed under bivariate trait analyses [8284]. Under this analytical framework, the phenotypic covariance (ρ 2 P ) is regarded as a function composed of the additive genetic (ρ 2 G ) and environmental (ρ 2 E ) covariances between two traits (denoted below as i and j).

$$ {\rho}_{P^{\left(i,j\right)}}={\rho}_{G^{\left(i,j\right)}}\sqrt{{h_i}^2{h_j}^2}+{\rho}_{E^{\left(i,j\right)}}\sqrt{\left(1-{h_i}^2\right)\left(1-{h_j}^2\right)} $$

These parameters are estimated using an estimation-maximization algorithm by jointly utilizing all available pedigree information with a multivariate normal model for continuous traits and liability threshold model for discrete traits [8587].

Other statistical methods used were Spearman’s correlation scatter plots and Bland-Altman plots to test for the agreement between different methods of measuring DNA methylation. Heritability, association, and bivariate trait analyses were conducted using the SOLAR software package [88], and all other statistical analyses were conducted using the Stata 12.0 (Stata Corp, College Station, TX) package.

References

  1. Onat A, Can G, Yuksel H. Dysfunction of high-density lipoprotein and its apolipoproteins: new mechanisms underlying cardiometabolic risk in the population at large. Turk Kardiyol Dern Ars. 2012;40:368–85.

    Article  PubMed  Google Scholar 

  2. Despres JP, Cartier A, Cote M, Arsenault BJ. The concept of cardiometabolic risk: bridging the fields of diabetology and cardiology. Ann Med. 2008;40:514–23.

    Article  CAS  PubMed  Google Scholar 

  3. Lemieux I, Poirier P, Bergeron J, Almeras N, Lamarche B, Cantin B, et al. Hypertriglyceridemic waist: a useful screening phenotype in preventive cardiology? Can J Cardiol. 2007;23 Suppl B:23B–31.

    Article  PubMed  Google Scholar 

  4. Despres JP, Lemieux I. Abdominal obesity and metabolic syndrome. Nature. 2006;444:881–7.

    Article  CAS  PubMed  Google Scholar 

  5. Wang A, Li Z, Zhou Y, Wang C, Luo Y, Liu X, et al. Hypertriglyceridemic waist phenotype and risk of cardiovascular diseases in China: results from the Kailuan Study. Int J Cardiol. 2014;174:106–9.

    Article  PubMed  Google Scholar 

  6. Gasevic D, Carlsson AC, Lesser IA, Mancini GJ, Lear SA. The association between "hypertriglyceridemic waist" and sub-clinical atherosclerosis in a multiethnic population: a cross-sectional study. Lipids Health Dis. 2014;13:38.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Diaz-Santana MV, Suarez Perez EL, Ortiz Martinez AP, Guzman Serrano M, Perez Cardona CM. Association between the hypertriglyceridemic waist phenotype, prediabetes, and diabetes mellitus among adults in Puerto Rico. J Immigr Minor Health. 2016;18:102-9.

  8. Daniel M, Paquet C, Kelly SJ, Zang G, Rowley KG, McDermott R, et al. Hypertriglyceridemic waist and newly-diagnosed diabetes among remote-dwelling Indigenous Australians. Ann Hum Biol. 2013;40:496–504.

    Article  CAS  PubMed  Google Scholar 

  9. Hobkirk JP, King RF, Gately P, Pemberton P, Smith A, Barth JH, et al. The predictive ability of triglycerides and waist (hypertriglyceridemic waist) in assessing metabolic triad change in obese children and adolescents. Metab Syndr Relat Disord. 2013;11:336–42.

    Article  CAS  PubMed  Google Scholar 

  10. Egeland GM, Cao Z, Young TK. Hypertriglyceridemic-waist phenotype and glucose intolerance among Canadian Inuit: the International Polar Year Inuit Health Survey for Adults 2007–2008. CMAJ. 2011;183:E553–8.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Alavian SM, Motlagh ME, Ardalan G, Motaghian M, Davarpanah AH, Kelishadi R. Hypertriglyceridemic waist phenotype and associated lifestyle factors in a national population of youths: CASPIAN Study. J Trop Pediatr. 2008;54:169–77.

    Article  PubMed  Google Scholar 

  12. Esmaillzadeh A, Mirmiran P, Azadbakht L, Azizi F. Prevalence of the hypertriglyceridemic waist phenotype in Iranian adolescents. Am J Prev Med. 2006;30:52–8.

    Article  PubMed  Google Scholar 

  13. Foucan L, Maimaitiming S, Larifla L, Hedreville S, Deloumeaux J, Joannes MO, et al. Adiponectin gene variants, adiponectin isoforms and cardiometabolic risk in type 2 diabetic patients. J Diabetes Investig. 2014;5:192–8.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Brisson D, St-Pierre J, Santure M, Hudson TJ, Despres JP, Vohl MC, et al. Genetic epistasis in the VLDL catabolic pathway is associated with deleterious variations on triglyceridemia in obese subjects. Int J Obes (Lond). 2007;31:1325–33.

    Article  CAS  Google Scholar 

  15. Pollex RL, Hanley AJ, Zinman B, Harris SB, Hegele RA. Clinical and genetic associations with hypertriglyceridemic waist in a Canadian aboriginal population. Int J Obes (Lond). 2006;30:484–91.

    Article  CAS  Google Scholar 

  16. St-Pierre J, Vohl MC, Despres JP, Gaudet D, Poirier P. Genetic aspects of diabetes and its cardiovascular complications: contribution of genetics to risk assessment and clinical management. Can J Cardiol. 2005;21:199–209.

    CAS  PubMed  Google Scholar 

  17. Liu Y, Wang J, Zhang R, Zhang Y, Xu Q, Zhang J, et al. A good response to oil with medium- and long-chain fatty acids in body fat and blood lipid profiles of male hypertriglyceridemic subjects. Asia Pac J Clin Nutr. 2009;18:351–8.

    CAS  PubMed  Google Scholar 

  18. Xue C, Liu Y, Wang J, Zhang R, Zhang Y, Zhang J, et al. Consumption of medium- and long-chain triacylglycerols decreases body fat and blood triglyceride in Chinese hypertriglyceridemic subjects. Eur J Clin Nutr. 2009;63:879–86.

    Article  CAS  PubMed  Google Scholar 

  19. Luttmer R, Spijkerman AM, Kok RM, Jakobs C, Blom HJ, Serne EH, et al. Metabolic syndrome components are associated with DNA hypomethylation. Obes Res Clin Pract. 2013;7:e106–15.

    Article  PubMed  Google Scholar 

  20. Carless MA, Kulkarni H, Kos MZ, Charlesworth J, Peralta JM, Goring HH, et al. Genetic effects on DNA methylation and its potential relevance for obesity in Mexican Americans. PLoS One. 2013;8:e73950.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Huang RC, Galati JC, Burrows S, Beilin LJ, Li X, Pennell CE, et al. DNA methylation of the IGF2/H19 imprinting control region and adiposity distribution in young adults. Clin Epigenetics. 2012;4:21.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Irvin MR, Zhi D, Joehanes R, Mendelson M, Aslibekyan S, Claas SA, et al. Epigenome-wide association study of fasting blood lipids in the genetics of lipid-lowering drugs and diet network study. Circulation. 2014;130:565–72.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Frazier-Wood AC, Aslibekyan S, Absher DM, Hopkins PN, Sha J, Tsai MY, et al. Methylation at CPT1A locus is associated with lipoprotein subfraction profiles. J Lipid Res. 2014;55:1324–30.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Lemas DJ, Wiener HW, O'Brien DM, Hopkins S, Stanhope KL, Havel PJ, et al. Genetic polymorphisms in carnitine palmitoyltransferase 1A gene are associated with variation in body composition and fasting lipid traits in Yup'ik Eskimos. J Lipid Res. 2012;53:175–84.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Razzouk L, Muntner P. Ethnic, gender, and age-related differences in patients with the metabolic syndrome. Curr Hypertens Rep. 2009;11:127–32.

    Article  PubMed  Google Scholar 

  26. Kurian AK, Cardarelli KM. Racial and ethnic differences in cardiovascular disease risk factors: a systematic review. Ethn Dis. 2007;17:143–52.

    PubMed  Google Scholar 

  27. Cossrow N, Falkner B. Race/ethnic issues in obesity and obesity-related comorbidities. J Clin Endocrinol Metab. 2004;89:2590–4.

    Article  CAS  PubMed  Google Scholar 

  28. MacCluer JW, Stern MP, Almasy L, Atwood LA, Blangero J, Comuzzie AG, et al. Genetics of atherosclerosis risk factors in Mexican Americans. Nutr Rev. 1999;57:S59–65.

    Article  CAS  PubMed  Google Scholar 

  29. Voruganti VS, Lopez-Alvarenga JC, Nath SD, Rainwater DL, Bauer R, Cole SA, et al. Genetics of variation in HOMA-IR and cardiovascular risk factors in Mexican-Americans. J Mol Med (Berl). 2008;86:303–11.

    Article  CAS  Google Scholar 

  30. Kulkarni H, Kos MZ, Neary J, Dyer TD, Kent Jr JW, Goring HH, et al. Novel epigenetic determinants of type 2 diabetes in Mexican-American families. Hum Mol Genet. 2015;24:5330–44.

    Article  CAS  PubMed  Google Scholar 

  31. Chen YA, Lemire M, Choufani S, Butcher DT, Grafodatskaya D, Zanke BW, et al. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics. 2013;8:203–9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Sobel ME. Asymptotic confidence intervals for indirect effects in structural equation models. In: Leinhardt S, editor. Sociological Methodology. Washington, D.C.: American Sociological Association; 1982. p. 290–312.

    Google Scholar 

  33. Arsenault BJ, Lemieux I, Despres JP, Wareham NJ, Kastelein JJ, Khaw KT, et al. The hypertriglyceridemic-waist phenotype and the risk of coronary artery disease: results from the EPIC-Norfolk prospective population study. CMAJ. 2010;182:1427–32.

    Article  PubMed Central  PubMed  Google Scholar 

  34. Chaput JP, McNeil J, Despres JP, Bouchard C, Tremblay A. Short sleep duration as a risk factor for the development of the metabolic syndrome in adults. Prev Med. 2013;57:872–7.

    Article  PubMed  Google Scholar 

  35. Haack RL, Horta BL, Gigante DP, Barros FC, Oliveira I, Silveira VM. The hypertriglyceridemic waist phenotype in young adults from the Southern Region of Brazil. Cad Saude Publica. 2013;29:999–1007.

    Article  PubMed  Google Scholar 

  36. Padwal RS, Pajewski NM, Allison DB, Sharma AM. Using the Edmonton obesity staging system to predict mortality in a population-representative cohort of people with overweight and obesity. CMAJ. 2011;183:E1059–66.

    Article  PubMed Central  PubMed  Google Scholar 

  37. Lemieux I, Almeras N, Mauriege P, Blanchet C, Dewailly E, Bergeron J, et al. Prevalence of 'hypertriglyceridemic waist' in men who participated in the Quebec Health Survey: association with atherogenic and diabetogenic metabolic risk factors. Can J Cardiol. 2002;18:725–32.

    PubMed  Google Scholar 

  38. Solati M, Ghanbarian A, Rahmani M, Sarbazi N, Allahverdian S, Azizi F. Cardiovascular risk factors in males with hypertriglycemic waist (Tehran Lipid and Glucose Study). Int J Obes Relat Metab Disord. 2004;28:706–9.

    Article  CAS  PubMed  Google Scholar 

  39. Czernichow S, Bruckert E, Bertrais S, Galan P, Hercberg S, Oppert JM. Hypertriglyceridemic waist and 7.5-year prospective risk of cardiovascular disease in asymptomatic middle-aged men. Int J Obes (Lond). 2007;31:791–6.

    CAS  Google Scholar 

  40. Li Y, Zhou C, Shao X, Liu X, Guo J, Zhang Y, et al. Hypertriglyceridemic waist phenotype and chronic kidney disease in a chinese population aged 40 years and older. PLoS One. 2014;9:e92322.

    Article  PubMed Central  PubMed  Google Scholar 

  41. Zhang M, Gao Y, Chang H, Wang X, Liu D, Zhu Z, et al. Hypertriglyceridemic-waist phenotype predicts diabetes: a cohort study in Chinese urban adults. BMC Public Health. 2012;12:1081.

    Article  PubMed Central  PubMed  Google Scholar 

  42. Kahn HS, Valdez R. Metabolic risks identified by the combination of enlarged waist and elevated triacylglycerol concentration. Am J Clin Nutr. 2003;78:928–34.

    CAS  PubMed  Google Scholar 

  43. Keller C, Chintapalli K, Lancaster J. Correlation of anthropometry with CT in Mexican-American women. Res Nurs Health. 1999;22:145–53.

    Article  CAS  PubMed  Google Scholar 

  44. He S, Zheng Y, Shu Y, He J, Wang Y, Chen X. Hypertriglyceridemic waist might be an alternative to metabolic syndrome for predicting future diabetes mellitus. PLoS One. 2013;8:e73292.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Carlsson AC, Riserus U, Arnlov J. Hypertriglyceridemic waist phenotype is associated with decreased insulin sensitivity and incident diabetes in elderly men. Obesity (Silver Spring). 2014;22:526–9.

    Article  CAS  Google Scholar 

  46. Willer CJ, Schmidt EM, Sengupta S, Peloso GM, Gustafsson S, Kanoni S, et al. Discovery and refinement of loci associated with lipid levels. Nat Genet. 2013;45:1274–83.

    Article  CAS  PubMed  Google Scholar 

  47. Zhou L, He M, Mo Z, Wu C, Yang H, Yu D, et al. A genome wide association study identifies common variants associated with lipid levels in the Chinese population. PLoS One. 2013;8:e82420.

    Article  PubMed Central  PubMed  Google Scholar 

  48. Weissglas-Volkov D, Aguilar-Salinas CA, Nikkola E, Deere KA, Cruz-Bautista I, Arellano-Campos O, et al. Genomic study in Mexicans identifies a new locus for triglycerides and refines European lipid loci. J Med Genet. 2013;50:298–308.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Asselbergs FW, Guo Y, van Iperen EP, Sivapalaratnam S, Tragante V, Lanktree MB, et al. Large-scale gene-centric meta-analysis across 32 studies identifies multiple lipid loci. Am J Hum Genet. 2012;91:823–38.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Chambers JC, Elliott P, Zabaneh D, Zhang W, Li Y, Froguel P, et al. Common genetic variation near MC4R is associated with waist circumference and insulin resistance. Nat Genet. 2008;40:716–8.

    Article  CAS  PubMed  Google Scholar 

  51. Lindgren CM, Heid IM, Randall JC, Lamina C, Steinthorsdottir V, Qi L, et al. Genome-wide association scan meta-analysis identifies three Loci influencing adiposity and fat distribution. PLoS Genet. 2009;5:e1000508.

    Article  PubMed Central  PubMed  Google Scholar 

  52. Heard-Costa NL, Zillikens MC, Monda KL, Johansson A, Harris TB, Fu M, et al. NRXN3 is a novel locus for waist circumference: a genome-wide association study from the CHARGE Consortium. PLoS Genet. 2009;5:e1000539.

    Article  PubMed Central  PubMed  Google Scholar 

  53. Liu CT, Monda KL, Taylor KC, Lange L, Demerath EW, Palmas W, et al. Genome-wide association of body fat distribution in African ancestry populations suggests new loci. PLoS Genet. 2013;9:e1003681.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Randall JC, Winkler TW, Kutalik Z, Berndt SI, Jackson AU, Monda KL, et al. Sex-stratified genome-wide association studies including 270,000 individuals show sexual dimorphism in genetic loci for anthropometric traits. PLoS Genet. 2013;9:e1003500.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  55. McGarry JD, Brown NF. The mitochondrial carnitine palmitoyltransferase system. From concept to molecular analysis. Eur J Biochem. 1997;244:1–14.

    Article  CAS  PubMed  Google Scholar 

  56. Gobin S, Bonnefont JP, Prip-Buus C, Mugnier C, Ferrec M, Demaugre F, et al. Organization of the human liver carnitine palmitoyltransferase 1 gene (CPT1A) and identification of novel mutations in hypoketotic hypoglycaemia. Hum Genet. 2002;111:179–89.

    Article  CAS  PubMed  Google Scholar 

  57. Bonnefont JP, Djouadi F, Prip-Buus C, Gobin S, Munnich A, Bastin J. Carnitine palmitoyltransferases 1 and 2: biochemical, molecular and medical aspects. Mol Aspects Med. 2004;25:495–520.

    Article  CAS  PubMed  Google Scholar 

  58. AlSaleh A, Sanders TA, O'Dell SD. Effect of interaction between PPARG, PPARA and ADIPOQ gene variants and dietary fatty acids on plasma lipid profile and adiponectin concentration in a large intervention study. Proc Nutr Soc. 2012;71:141–53.

    Article  CAS  PubMed  Google Scholar 

  59. Patterson AD, Slanar O, Krausz KW, Li F, Hofer CC, Perlik F, et al. Human urinary metabolomic profile of PPARalpha induced fatty acid beta-oxidation. J Proteome Res. 2009;8:4293–300.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Robitaille J, Houde A, Lemieux S, Perusse L, Gaudet D, Vohl MC. Variants within the muscle and liver isoforms of the carnitine palmitoyltransferase I (CPT1) gene interact with fat intake to modulate indices of obesity in French-Canadians. J Mol Med (Berl). 2007;85:129–37.

    Article  CAS  Google Scholar 

  61. Gagnon F, Aissi D, Carrie A, Morange PE, Tregouet DA. Robust validation of methylation levels association at CPT1A locus with lipid plasma levels. J Lipid Res. 2014;55:1189–91.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  62. Murphy AJ, Yvan-Charvet L. Adipose modulation of ABCG1 uncovers an intimate link between sphingomyelin and triglyceride storage. Diabetes. 2015;64:689–92.

    Article  CAS  PubMed  Google Scholar 

  63. Pfeiffer L, Wahl S, Pilling LC, Reischl E, Sandling JK, Kunze S, et al. DNA methylation of lipid-related genes affects blood lipid levels. Circ Cardiovasc Genet. 2015;8:334–42.

    Article  CAS  PubMed  Google Scholar 

  64. Frisdal E, Le Lay S, Hooton H, Poupel L, Olivier M, Alili R, et al. Adipocyte ATP-binding cassette G1 promotes triglyceride storage, fat mass growth, and human obesity. Diabetes. 2015;64:840–55.

    Article  CAS  PubMed  Google Scholar 

  65. Guay SP, Brisson D, Lamarche B, Gaudet D, Bouchard L. Epipolymorphisms within lipoprotein genes contribute independently to plasma lipid levels in familial hypercholesterolemia. Epigenetics. 2014;9:718–29.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  66. Olivier M, Tanck MW, Out R, Villard EF, Lammers B, Bouchareychas L, et al. Human ATP-binding cassette G1 controls macrophage lipoprotein lipase bioavailability and promotes foam cell formation. Arterioscler Thromb Vasc Biol. 2012;32:2223–31.

    Article  CAS  PubMed  Google Scholar 

  67. Demerath EW, Guan W, Grove ML, Aslibekyan S, Mendelson M, Zhou YH, et al. Epigenome-wide association study (EWAS) of BMI, BMI change and waist circumference in African American adults identifies multiple replicated loci. Hum Mol Genet. 2015;24:4464–79.

    Article  CAS  PubMed  Google Scholar 

  68. Hidalgo B, Irvin MR, Sha J, Zhi D, Aslibekyan S, Absher D, et al. Epigenome-wide association study of fasting measures of glucose, insulin, and HOMA-IR in the genetics of lipid lowering drugs and diet network study. Diabetes. 2014;63:801–7.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  69. Ding J, Reynolds LM, Zeller T, Muller C, Mstat KL, Nicklas BJ, et al. Alterations of a cellular cholesterol metabolism network is a molecular feature of obesity-related type 2 diabetes and cardiovascular disease. Diabetes. 2015;64:3464-74.

  70. Yoshihara E, Masaki S, Matsuo Y, Chen Z, Tian H, Yodoi J. Thioredoxin/Txnip: redoxisome, as a redox switch for the pathogenesis of diseases. Front Immunol. 2014;4:514.

    Article  PubMed Central  PubMed  Google Scholar 

  71. Sano R, Reed JC. ER stress-induced cell death mechanisms. Biochim Biophys Acta. 1833;2013:3460–70.

    Google Scholar 

  72. van Greevenbroek MM, Vermeulen VM, Feskens EJ, Evelo CT, Kruijshoop M, Hoebee B, et al. Genetic variation in thioredoxin interacting protein (TXNIP) is associated with hypertriglyceridaemia and blood pressure in diabetes mellitus. Diabet Med. 2007;24:498–504.

    Article  PubMed  Google Scholar 

  73. McGuinness D, McGlynn LM, Johnson PC, MacIntyre A, Batty GD, Burns H, et al. Socio-economic status is associated with epigenetic differences in the pSoBid cohort. Int J Epidemiol. 2012;41:151–60.

    Article  PubMed  Google Scholar 

  74. Youngson NA, Morris MJ. What obesity research tells us about epigenetic mechanisms. Philos Trans R Soc Lond B Biol Sci. 2013;368:20110337.

    Article  PubMed Central  PubMed  Google Scholar 

  75. Mitchell BD, Kammerer CM, Blangero J, Mahaney MC, Rainwater DL, Dyke B, et al. Genetic and environmental contributions to cardiovascular risk factors in Mexican Americans. The San Antonio Family Heart Study. Circulation. 1996;94:2159–70.

    Article  CAS  PubMed  Google Scholar 

  76. Gomez-Huelgas R, Bernal-Lopez MR, Villalobos A, Mancera-Romero J, Baca-Osorio AJ, Jansen S, et al. Hypertriglyceridemic waist: an alternative to the metabolic syndrome? Results of the IMAP Study (multidisciplinary intervention in primary care). Int J Obes (Lond). 2011;35:292–9.

    Article  CAS  Google Scholar 

  77. Expert Committee on the Diagnosis and Classification of Diabetes Mellitus. Report of the expert committee on the diagnosis and classification of diabetes mellitus. Diabetes Care. 2003;26(1):S5-20.

  78. Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, et al. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013;29:189–96.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  79. Reinius LE, Acevedo N, Joerink M, Pershagen G, Dahlen SE, Greco D, et al. Differential DNA methylation in purified human blood cells: implications for cell lineage and studies on disease susceptibility. PLoS One. 2012;7:e41361.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  80. Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012;13:86.

    Article  PubMed Central  PubMed  Google Scholar 

  81. Jaffe AE, Irizarry RA. Accounting for cellular heterogeneity is critical in epigenome-wide association studies. Genome Biol. 2014;15:R31.

    Article  PubMed Central  PubMed  Google Scholar 

  82. Melton PE, Rutherford S, Voruganti VS, Goring HH, Laston S, Haack K, et al. Bivariate genetic association of KIAA1797 with heart rate in American Indians: the Strong Heart Family Study. Hum Mol Genet. 2010;19:3662–71.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  83. Diego VP, Goring HH, Cole SA, Almasy L, Dyer TD, Blangero J, et al. Fasting insulin and obesity-related phenotypes are linked to chromosome 2p: the Strong Heart Family Study. Diabetes. 2006;55:1874–8.

    Article  CAS  PubMed  Google Scholar 

  84. Almasy L, Dyer TD, Blangero J. Bivariate quantitative trait linkage analysis: pleiotropy versus co-incident linkages. Genet Epidemiol. 1997;14:953–8.

    Article  CAS  PubMed  Google Scholar 

  85. Williams JT, Begleiter H, Porjesz B, Edenberg HJ, Foroud T, Reich T, et al. Joint multipoint linkage analysis of multivariate qualitative and quantitative traits. II. Alcoholism and event-related potentials. Am J Hum Genet. 1999;65:1148–60.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  86. Williams JT, Van Eerdewegh P, Almasy L, Blangero J. Joint multipoint linkage analysis of multivariate qualitative and quantitative traits. I. Likelihood formulation and simulation results. Am J Hum Genet. 1999;65:1134–47.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  87. Almasy L, Blangero J. Variance component methods for analysis of complex phenotypes. Cold Spring Harb Protoc. 2010;2010:pdb top77.

    Article  PubMed Central  PubMed  Google Scholar 

  88. Almasy L, Blangero J. Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998;62:1198–211.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Financial support for this study was provided by the NIDDK (DK087749). SNP genotyping and SOLAR is supported by NIMH (MH078143, MH078111 MH083824, MH059490). Collection of family data and blood samples in the San Antonio Family Heart Study was supported by a NHLBI program project grant (HL045522). This investigation was conducted in facilities constructed with support from Research Facilities Improvement Program Grant Number C06 RR017515 from the National Center for Research Resources, National Institutes of Health.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Manju Mamtani.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MAC and JB conceived the study. MAC generated the DNA methylation data. JB provided the samples for analysis. MM, HK, TDD, HHHG, JN, SAC, JWK, LA, MCM, and JB contributed to pedigree verification and statistical analysis. AGC, MCM, and JB generated phenotypic data. MAC, SK, DCG, JEC, RD, and JB provided genotyping. MM, HK, and MAC drafted the manuscript with editing from JB. All authors read and approved the final draft.

Additional files

Additional file 1:

This file contains Supplementary Table 1. (XLSX 782 MB)

Additional file 2:

This file contains Supplementary Figs. 1 and 2. (DOCX 208 KB)

Additional file 3:

This file contains Supplementary Table 2. (XLSX 314 MB)

Additional file 4:

This file contains Supplementary Table 3. (XLSX 1982 KB)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mamtani, M., Kulkarni, H., Dyer, T.D. et al. Genome- and epigenome-wide association study of hypertriglyceridemic waist in Mexican American families. Clin Epigenet 8, 6 (2016). https://doi.org/10.1186/s13148-016-0173-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13148-016-0173-x

Keywords