Abstract
Aim:
The aim of this study was to investigate the association of four variants of the hepatic lipase (HL [or LIPC]) gene with various lipid parameters among South Indian subjects with normal glucose tolerance (NGT).
Subjects and Methods:
In total, 747 NGT subjects were randomly selected from the Chennai Urban Rural Epidemiological Study (CURES). Serum triglycerides, serum cholesterol, and high-density lipoprotein cholesterol (HDL-C) were measured using a Hitachi-912 autoanalyzer (Roche Diagnostics GmbH, Mannheim, Germany). Genotyping of HL gene variants was done by the polymerase chain reaction–restriction fragment length polymorphism method, and 20% of samples were sequenced to validate the genotypes obtained. Haplotype analysis was also carried out.
Results:
The TT genotype of the rs1800588 C/T (C-480T) polymorphism was significantly associated with hypertriglyceridemia, with an adjusted odds ratio of 2.58 (95% confidence interval 1.38–4.85, P=0.003), whereas those with the CC genotype of the rs6074 A/C (Thr479Thr) had significantly lower HDL-C levels (41.3±9.8 mg/dL) compared with the AA genotype (43.6±10.2 mg/dL, P=0.02). Haplotype analysis showed the TGC haplotype was significantly associated with low HDL-C levels.
Conclusions:
Among South Indian subjects without diabetes, the rs1800588 C/T (C-480T) and rs6074 C/A (Thr479Thr) variants of the HL gene are associated with hypertriglyceridemia and low HDL-C, respectively. The TGC haplotype was significantly associated with low HDL-C.
Introduction
The prevalence of T2D and CHD has doubled in India in the past three decades. 3,4 It is well established that Asian Indians have an increased susceptibility to T2D and premature CHD compared with Europeans. This has been explained by a higher frequency of hyperinsulinemia, 4 insulin resistance, 5 dyslipidemia with low levels of HDL-C, 6 and increased visceral fat despite a lower body mass index, and this unique phenotype is referred to as the “Asian Indian phenotype.” 7,8 However, it is not clear whether the Asian Indian phenotype represents a unique genetic susceptibility or results from lifestyle-related factors. Support for genetic susceptibility comes from the observation that hyperinsulinemia was demonstrated in newborn Asian Indian children but not in their United Kingdom counterparts 9 and for some differences in genetic association with T2D in Asian Indians. 10,11
One of the most striking features of the Asian Indian phenotype is the presence of dyslipidemia, with hypertriglyceridemia and low HDL-C being the most prominent features. 6 Indeed, a low level of HDL-C is seen in 60–70% of adults in urban Indians 12 and even in urban adolescent children. 13,14 However, few studies have looked at the possible causes of the low level of HDL-C and hypertriglyceridemia. It is likely that lifestyle factors (e.g., high carbohydrate diet 15 or physical inactivity 16 ) contribute to the dyslipidemia in our population. However, it is also possible that genetic factors play a role in this. Association of genetic variants with low HDL levels has been studied earlier in this population, with respect to the lipoprotein lipase gene. 17 Hepatic lipase (HL) is similar to lipoprotein lipase and plays a key role in the metabolism of pro-atherogenic and anti-atherogenic lipoproteins, affecting their plasma levels as well as their physicochemical properties. HL is a glycoprotein that catalyzes the hydrolysis of lipoprotein triacylglycerols and phospholipids. The majority of HL is synthesized and secreted by the liver and is bound to heparin sulfate proteoglycans on the surface of sinusoidal endothelial cells and external surfaces of microvilli of parenchymal cells in the space of Disse, 18 promoting the uptake of HDL and apolipoprotein-B-containing remnant particles. 19 It is estimated that 30–45% of variability in HL activity is genetically determined. 20 Several independent studies have documented the association of some of the prominent single nucleotide polymorphisms (SNPs) in the HL gene. 21 –25 Although there are some studies on the variants in coding regions of the HL gene in Turkish 26 and Japanese 27 population, little is known about the association of variants of the HL gene with dyslipidemia in South Asians.
The aim of the present study was to investigate the association of four variants in the HL gene—namely, rs1800588 (C-480T) (C/T), rs6074 (Thr479Thr) (C/A), rs690 (Val155Val) (G/T) and rs6083 (Ser193Asn) (G/A)—with dyslipidemia in an Asian Indian (urban South Indian) population.
Research Design and Methods
Subjects and study design
The study subjects were selected from the urban component of the Chennai Urban Rural Epidemiology Study (CURES), an epidemiology study conducted on a representative population (>20 years) of Chennai (formerly Madras), the fourth largest city in India. 28 In brief, 26,001 adult subjects (>20 years of age) were recruited in Phase 1 of CURES using a systematic random sampling method covering the whole of Chennai City. This identified 1,529 self-reported “diabetes subjects.” In Phase 2, all self-reported or “known diabetes subjects” were invited to our center for detailed studies, of which 1,382 responded (response rate, 90.4%). In Phase 3, every 10th subject from Phase 1 (n=2,600), excluding those with self-reported diabetes, were invited to undergo an oral glucose tolerance test (OGTT). For the present study, 1,000 unrelated subjects with normal glucose tolerance (NGT), defined as those who had a fasting plasma glucose level of <5.6 mmol/L (100 mg/dL) and 2-h plasma glucose value of 7.8 mmol/L (140 mg/dL), were selected from Phase 3 of CURES. The ratio of males versus females in the present study was 466:534. Those who were on lipid-lowering drugs such as statins, fibrates, and niacin were excluded from the study (n=134). On the basis of the National Cholesterol Education Program Adult Treatment Panel III guidelines, 29 the study population was divided into those with normal HDL-C (≥45 mg/dL among males; ≥50 mg/dL among females) and low HDL-C (≤45 mg/dL among males; ≤50 mg/dL among females), normal triglycerides (≤150 mg/dL) and high triglycerides (≥150 mg/dL), normal LDL-C (≤110 mg/dL) and high LDL-C (≥110 mg/dL), and normal total cholesterol (≤200 mg/dL) and high total cholesterol (≥200 mg/dL).
Written informed consent was obtained from each study participant, and the studies were approved by the Institutional Ethics Committee of the Madras Diabetes Research Foundation.
Phenotype measurements
Anthropometric measurements including weight, height, and waist measurements were obtained using standardized techniques. The body mass index (BMI) was calculated using the formula of weight (kg)/height squared (m2). Fasting blood sample was taken for estimation of levels of plasma glucose and serum lipids. Biochemical assays were carried out using a Hitachi-912 autoanalyzer (Roche Diagnostics GmbH, Mannheim, Germany) using kits supplied by Roche Diagnostics. Glycated hemoglobin was estimated by the high-performance liquid chromatography method using the Variant™ machine (Bio-Rad, Hercules, CA). Serum insulin concentration was estimated using kits from Dako (Glostrup, Denmark).
Levels of serum cholesterol (cholesterol oxidase–peroxidase–amidopyrine method; Roche Diagnostics), serum triglycerides (glycerol phosphate oxidase–peroxidase–amidopyrine method; Roche Diagnostics), and HDL-C (direct method, polyethylene glycol–pretreated enzymes; Roche Diagnostics) were measured using the Hitachi-912 autoanalyzer. Insulin resistance was calculated using the homeostasis assessment model using the formula fasting insulin (IU/mL)×fasting glucose (mmol/L)/22.5. 30
Genetic analysis
EDTA-anticoagulated venous blood samples were collected from all study subjects, and the genomic DNA was isolated from whole blood by proteinase K digestion followed by phenol–chloroform extraction. 31 The concentration and purity was estimated with a NanoDrop™ 1000 spectrophotometer (Thermo Scientific, Wilmington, DE).
We amplified the regions of the HL gene using the GeneAmp® PCR [polymerase chain reaction] system 9700 (Applied Biosystems, Foster City, CA), using the following primers: forward 5′ GGC ATC TTT GCT TCT TCG TC 3′, reverse 5′ CTG GCT CAG GAA AGT GGT GT′; forward 5′ TTCATCCAGGCAGCTCTTCT 3′, reverse 5′ CCAATCTTGTGCGTTCCAC 3′; forward 5′ CTTTCCCATTAGGGCTGGAT 3′, reverse 5′ TCCTATGGGCTGTTTGATGC 3′; and forward 5′ TGC TGT GTT TGC TTC CTG TT 3′, reverse 5′ AAG GCA GCC ATT CCA GAT AA3′ (Sigma, Bangalore, India). Restriction fragment length polymorphism was carried out using NlaIII (rs1800588 C/T), Xmn1 (rs690 G/T), Tsp509I (rs6083 G/A), and ScrfI (rs6074 A/C) restriction enzymes (New England Biolabs, Inc., Ipswich, MA). Agarose gel electrophoresis was used to check the amplification of the PCR and the restriction enzyme–digested products. The percentage of agarose in the gel varied, depending on the expected size(s) of the fragment(s). Ethidium bromide was included in the gel matrix to enable fluorescent visualization of the DNA fragments under ultraviolet light. To ensure that the genotyping was of adequate quality, we performed random duplicates in 20% of the samples. The assays were performed by a technician who was masked to the phenotype, and there was 99% concordance in the genotyping. Furthermore, a few variants were confirmed by direct sequencing with an ABI 310 genetic analyzer (Applied Biosystems).
Statistical analysis
Statistical Package for Social Sciences (SPSS) Windows (version 15.0) (SPSS, Inc., Chicago, IL) was used for statistical analysis. Allele frequencies were estimated by gene counting. Agreement with Hardy–Weinberg expectations was tested using a χ2 goodness-of-fit test. The χ2 test or Fisher's exact t test as appropriate was used to compare the proportions of genotypes or alleles. One-way analysis of variance was used to compare groups for continuous variables. Logistic regression analysis was done using HDL and triglycerides (dichotomized) as the dependent variable and the genotypes (risk factor) and age, gender, and BMI (covariates adjusted in the analysis) as independent variables. Significant P values obtained were then corrected for multiple testing.
Haplotype analysis was performed using Haploview software (available at
Results
The clinical and biochemical characteristics of the study subjects are shown in Table 1. Among the 1,000 NGT study subjects the ratio of male:female NGT subjects in the present study was 466:534. Table 2 shows the distribution of genotype frequencies among the study subjects stratified based on serum triglyceride levels. We found the rs1800588 C-480T variant of the HL gene to be significantly associated with hypertriglyceridemia. All the genotype frequencies of the four variants in the study subjects were in Hardy–Weinberg equilibrium. The frequency of the TT genotype was 11.3% among the hypertriglyceridemic group, compared with 5% among the normotriglyceridemic group (P=0.005). The minor allele frequency of the T allele was also significantly higher among the hypertriglyceridemic group (29.2%) compared with the normotriglyceridemic group (23.0%) (P=0.02). Logistic regression analysis showed that the unadjusted odds ratio (OR) for the TT genotype was 2.51 (95% confidence interval [CI] 1.34–4.68, P=0.004). This association persisted even after adjusting for age, gender, and BMI, with an OR of 2.58 (95% CI 1.38–4.85, P=0.003). None of the other variants showed any significant difference in the distribution of genotype frequency between the two groups. We also performed the power calculation to evaluate whether our sample size had sufficient power to detect the observed difference in the proportion of the CC and CT genotypes in the two groups (controls and cases) and found that the power was 0.75.
Data are mean±SD values.
BMI, body mass index; HDL-C, high-density lipoprotein-cholesterol; LDL-C, low-density lipoprotein-cholesterol.
Adjusted for age, gender, and body mass index.
P value indicates a significant difference.
CI, confidence interval; MAF, minor allele frequency; OR, odds ratio.
Testing the association of the HL variants under various genetic models
Various genetic models were tested for the variant that was significantly associated with serum triglycerides (Table 3).
Adjusted for age, gender, and body mass index.
P value indicates a significant difference.
CI, confidence interval; OR, odds ratio; TG, triglyceride.
The association of the rs1800588 (C-480T) variant of the HL gene with triglycerides was assessed based on the additive genetic model. We obtained an OR of 2.73 (95% CI 1.44–5.18, P=0.002) for the TT genotype, after adjusting for age, sex, and BMI. Under the dominant model, the OR was 1.23 (95% CI 0.86–1.75, P=0.24) for the CT+TT genotype. Under the recessive model, we observed an adjusted OR of 2.69 (95% CI 1.44–5.02, P=0.002) for the TT genotype. All these observations were done after adjusting for age, sex, and BMI.
Table 4 shows the distribution of genotype frequencies among the study subjects stratified based on HDL-C levels according to the modified Adult Treatment Panel III guidelines 29 as mentioned earlier. The rs6074 (Thr479Thr) (A/C) variant of the HL gene was significantly associated with low HDL-C levels. The logistic regression analysis showed that the unadjusted OR for the CC genotype was 1.47 (95% CI 0.97–2.21, P=0.06); after adjusting for age, sex, and BMI, the OR was 1.25 (95% CI 0.81–1.95, P=0.30). Similarly, the rs6083 (Ser193Asn) (G/A) variant of the HL gene was significantly associated with normal HDL-C levels among the NGT subjects. We observed a significantly higher GA frequency (47.5%) among the normal HDL-C group, compared with the low HDL-C group (40.7%) (P=0.05). However, the OR significance was lost after adjusting for age, gender, and BMI: OR=0.72 (95% CI 0.51–1.02, P=0.06).
Adjusted for age, gender, and body mass index.
P value indicates a significant difference.
CI, confidence interval; MAF, minor allele frequency; OR, odds ratio.
Because HDL-C cut points are different in males and females, separate analysis was done in males and females. Table 5 shows that among the males, the GA genotype of the rs6083 (Ser193Asn) (G/A) variant was also significantly associated with normal HDL-C levels. The frequency of the GA genotype was significantly associated with normal HDL-C levels (54.1%) compared with low HDL-C levels (38.3%) (P=0.05). However, the significance was lost after adjusting for age and BMI. There was no significant association of HL gene variants with HDL-C among females. Table 6 shows that among males, none of the four HL variants showed any significant association with serum triglycerides; however, in the case of female NGT subjects, in the variant C-514T, the frequency of the TT genotype was significantly higher among the high triglyceride group (11.9%) compared with the normal triglyceride group (4.1%) (P=0.001). The minor allele frequency of the T allele was also significantly higher among the high triglyceride group (39.7%) compared with the normal triglyceride group (23.3%) (P<0.001).
P value indicates a significant difference.
HDL, high-density lipoprotein; MAF, minor allele frequency.
P value indicates a significant difference.
MAF, minor allele frequency; TG, triglyceride.
In Table 7, a linear regression analysis was done to estimate the effect size β (B) of the four HL gene variants with various quantitative traits among subjects without diabetes. In the case of the −514 C/T variant, we observed a significant effect size (B) of 0.06 with an SE of 0.21 (P=0.003). Also, in the case of the rs6074 (Thr479Thr) variant, a significant B of 0.79 was observed (P=0.01). The other variants did not show any significant association with any of the quantitative traits among the NGT subjects.
Adjusted for age, sex, and body mass index.
P value indicates a significant difference.
HDL-C, high-density lipoprotein-cholesterol; LDL-C, low-density lipoprotein-cholesterol; TG, triglyceride.
Table 8 shows the association between clinical and biochemical parameters with the four variants of the HL gene. In the case of the C-480T variant, the TT genotype was significantly associated with higher triglyceride levels compared with the other genotypes (P=0.002). In the case of the rs6074 (Thr479Thr) (A/C) variant, the CC genotype was significantly associated with lower HDL-C levels (41.3±9.8 mg/dL), compared with the AA genotype (43.6±10.2 mg/dL) (P=0.02). None of the parameters showed any significant association with rs690 (Val155Val) (G/T) variants in the HL gene among the study subjects. The significant associations observed were validated by performing multiple testing.
Data are mean±SD values.
0 is the homozygous normal genotype, 1 is heterozygous, and 2 is the homozygous variant genotype.
P adjusted for age and sex.
P value indicates a significant difference.
HDL-C, high-density lipoprotein-cholesterol; LDL-C, low-density lipoprotein-cholesterol.
Association of the haplotype with HDL-C and serum triglyceride levels
Figure 1 shows the LD pattern observed in the HL gene with respect to the four variants studied among NGT subjects. The LD value between rs690 (G/T) and rs6074 (A/C) was weak with r 2 =0.02, the LD between rs1800588 (C-482T) and rs690 (G/T) was very weak with r 2 =0.00, and the LD value between rs6083 (G/A) and rs6074 (A/C) was moderate with r 2 =0.13.

Linkage disequilibrium pattern of the hepatic lipase gene (four loci). Pairwise linkage disequilibrium coefficients D′×100 are shown in each cell.
Table 9 shows the haplotypes obtained with Haploview version 4.2 software. The haplotypes associated were tested for multiple correction (α/number of possible haplotypes), where α=0.05. Owing to lower LD values (r 2 =0.00) between the rs1800588 (C/T) (−C514T) variant and other variants, the rs1800588 (C/T) (−C514T) variant was not included in the haplotype block. So the other three SNPs were included in the haplotype block. Therefore, in combination with three SNPs, rs6074 (Thr479Thr) (C/A), rs690 (Val155Val) (G/T), and rs6083 (Ser193Asn) (G/A) in the coding region of the HL gene, eight haplotype combinations were obtained, of which the GAA and TGC haplotypes were taken for further analysis with various disease phenotypes. Because eight tests were performed corresponding to the eight haplotypes satisfying the selection criterion, a multiple correction was done using Bonferroni's correction. The true association of the haplotype was tested using permutation analysis in Haploview software, wherein 10,000 permutations were computed, to obtain the differences in haplotype frequencies between cases and controls. We observed that the difference in the proportion in the haplotypes GAA and TGC between the cases and controls was truly significant.
P value after multiple corrections (P value × number of tests/number of haplotypes).
Permutation P values were calculated based on 10,000 permutation replications.
P value indicates a significant difference.
Case, low high-density lipoprotein level; control, normal high-density lipoprotein level; NS, not significant.
For HDL-C levels, the frequency of the GAA haplotype was significantly higher among the normal HDL group (0.268) compared with the low HDL group (0.257) (P=0.006). The frequency of the TGC haplotype was observed to be significantly higher among the low HDL group (0.150) compared with the normal HDL group (0.096) (P=0.002).
In the case of serum triglycerides, the haplotype frequencies were not significantly different between the cases and controls (data not shown).
Association of TGC haplotype with quantitative traits among the subjects without diabetes
Table 10 shows the association of the TGC haplotype with various biochemical parameters. The TGC haplotype was significantly associated with hypertriglyceridemia: the subjects with the TGC haplotype showed significantly higher serum triglyceride levels (2.03±0.22 mg/dL) compared with the other haplotypes (1.99±0.2 mg/dL) (P=0.05). We also observed significantly lower HDL-C levels among the TGC carriers (41.4±9.8 mg/dL) compared with the other haplotypes (43.7±10.0 mg/dL) (P=0.003).
Data are mean±SD values.
P value after multiple corrections (P value × number of tests/number of haplotypes).
P value indicates a significant difference.
HDL-C, high-density lipoprotein; LDL, low-density lipoprotein.
Discussion
In this case-control study, we aimed to investigate the association of the HL gene with dyslipidemia. We genotyped four common variants: rs1800588 (C/T) (also denominated as −480 C/T or −514 C/T) located in the promoter region; rs6083 (Ser193Asn) (G/A) variant located in exon 4; rs690 (Val155Val), located in exon 3; and rs6074 (Thr479Thr) (A/C) variant, located in exon 9 of the HL gene.
We found that the T allele of the variant −480 C/T located in the promoter region of the HL gene was associated with higher triglyceride levels. Also, the TGC haplotype was significantly associated with low HDL-C levels among the study subjects. Studies among a Finnish population 33 also found the −480 C/T polymorphism to be associated with insulin resistance and high triglyceride levels. In our study the TT genotype was found to confer a 2.5 times higher risk of hypertriglyceridemia.
The proximal promoter region of the HL gene is an important binding site for upstream stimulatory factor-1, and the −480 C/T substitution might induce changes in the E-box motif from CACGTG to CATGTG. This disruption in the binding site is probably responsible for lowering the HL gene activity. 34 Because HL catalyzes the hydrolysis of TG from intermediate-density lipoprotein and LDL, 35 such defects in the HL gene due to the promoter variant −480 C/T may lead to hypertriglyceridemia. Studies on an interactive effect between HL (also known as LIPC [lipase, hepatic]) gene variants and apolipoprotein E2 and that included C-514T,G-250A variants showed strong association of the TT genotype of the −514 C/T variant with hypertriglyceridemia among healthy Canadian adults 36 and similar association with serum triglycerides among an Iranian population, 37 strengthening the finding among NGT subjects in the present study. Although a few studies have shown quite a significant association with HDL-C in other ethnic populations, 38 –40 the present study did not find any association with HDL-C.
The rs6074 (Thr479Thr) (A/C) variant located in exon 9 of the HL gene was significantly associated with low HDL-C levels among the study subjects; however, correction for age, gender, and BMI abolished the association. Moreover, the CC genotype was observed to be significantly associated with low HDL-C levels (41.3±9.8 mg/dL) even after correcting for age and sex. An earlier study on two white populations failed to show an association with HDL-C. 41 HL probably has anti-atherogenic functions involving the uptake of HDL-C by the liver as part of the reverse cholesterol transport pathway. 42 Although a synonymous variation, this might facilitate biophysical change in the conformation and energetic aspects of the gene, resulting in enhancement of the gene function. High HL activity is associated with an increase in small dense LDL particles and increased risk of vascular events. 43
With reference to the rs6083 (Ser193Asn) (G/A) variant, located in exon 4 of the HL gene, we observed significantly lower serum triglycerides levels among the subjects with the AA genotype compared with the GG genotype. Exon 4 of the HL gene is the region that binds to the lipoprotein substrate. 44 This points to the potential role played by this variant in enhancing the activity of the HL enzyme in the hydrolysis of the triglycerides of intermediate-density lipoproteins. Studies conducted among white populations 41,45 showed no significant association of this variant with HDL-C levels.
The haplotype analysis representing three variants in the HL gene showed the GAA haplotype as a protective haplotype, conferring lower risk toward low HDL-C levels among the study subjects, and the TGC haplotype proved to be a high-risk haplotype, conferring high risk toward low HDL-C level. The benefit of a haplotype-based analysis is that it captures all of the variation across a region, which aids in improving the ability to detect an association. Hence we can delineate the association of the HL gene variants with major dyslipidemic phenotypes among our study subjects without diabetes. Apart from the −480 C/T variant, there are not many studies that support the association of rs690 (Val155Val) (G/T), rs6083 (Ser193Asn) (G/A), and rs6074 (Thr479Thr) (C/A) variants with HDL-C and triglycerides. By identifying whole chromosomal regions, haplotypes have improved power and reproducibility in elucidation of disease–gene associations.
To circumvent the problem of population stratification, we performed a case-control study at six unlinked marker loci believed to be unrelated to the disease under study, but known to have allelic diversity among different populations. 46 The allele frequency difference was not statistically significant at any of the loci, indicating that the findings in this study were unlikely to be an artifact of population substructuring.
To our knowledge this is the first study of the association of four variants of the HL gene with HDL-C and serum triglycerides in Asian Indians. We conclude that in this ethnic group, the rs1800588 (C/T) (C-480T) variant is associated with hypertriglyceridemia, that the rs6074 (Thr479Thr) (A/C) variant is associated with low HDL-C levels, and that the haplotype TGC confers susceptibility to low HDL-C levels.
Footnotes
Acknowledgments
K.A.A. was supported by the Council of Scientific and Industrial Research. The Chennai Wellingdon Corporate Foundation supported the CURES field studies. This is the 109th publication from the CURES (CURES-109).
Author Disclosure Statement
No competing financial interests exist.
