1. Introduction
Lipid disorders are one of the most important predictors of atherosclerotic coronary artery disease (CAD). Polymorphic variants of key genes involved in lipid metabolism may also significantly influence the risk of CAD, modifying the system’s response. The
USF1 (upstream stimulatory factor 1) gene (1q23.3) is one of the regulatory genes of lipid metabolism [
1]. Its locus was initially identified by linkage mapping studies and it was associated with familial combined hyperlipidemia [
2,
3,
4]. The
USF1 gene encodes the USF1 polypeptide, a subunit of the dimeric transcription factor USF, which belongs to the DNA-binding proteins [
5]. USF regulates expression of numerous genes involved in cellular proliferation, cell cycle, stress and immune response, as well as lipid and carbohydrate metabolism, including apolipoprotein E gene (
APOE), ATP-binding cassette, sub-family A, member 1 gene (
ABCA1), apolipoprotein A-V gene (
APOA5), angiotensinogen gene (
AGT) and many more [
4,
5]. Through the impact on target genes, USF is involved in many processes like the synthesis of fatty acids, β-oxidation, binding and transport of lipids in the bloodstream, and prostanoid metabolism.
In many studies,
USF1 gene polymorphisms have been associated either with familial combined hyperlipidemia or one of its component traits [
2,
4,
6,
7,
8,
9,
10,
11,
12,
13,
14,
15]. They were also analyzed in the context of diabetes mellitus type II, metabolic syndrome [
8,
9,
10,
16] and CAD [
11,
13,
14]. However, there are significant differences in the observed results, especially those concerning the impact of the
USF1 gene polymorphisms on lipid levels and on the phenotype of the disease [
6,
17,
18,
19].
The presence of those substantial differences suggests that additional factors may affect the influence of USF1 genotypes on lipid parameters. Therefore, in the current study we tried to identify the factors that may modify these effects. We analyzed two SNPs of the USF1 gene, namely rs2516839 (c.161043331C>T) and rs3737787 (c.161039733C>T), and their possible associations with premature CAD. In addition to the classical case-control analysis, we also analyzed: (1) the impact of the USF1 genotype variants on serum lipid parameters; (2) possible interactions of both USF1 gene polymorphisms with traditional risk factors of atherosclerosis increasing the risk of CAD and (3) we searched for traditional risk factors influencing the differences in lipid parameters, dependent on USF1 gene alleles.
3. Discussion
The present work is the first study analyzing polymorphisms of the
USF1 gene in the Polish population. A limitation of the present study is the fact that the analyses were performed in the group with a relatively small number of participants. But it has also its advantages, because of ethnical homogeneity of the subjects (only white polish Caucasians, inhabitants of Upper Silesia region were included). The frequencies of alleles of the rs2516839 in our study are similar to those observed in other Caucasian populations including German [
9,
10], Finnish [
6,
13], French-Canadian [
14] and Australian of European origin [
7]. The frequencies of alleles of the rs3737787 polymorphism also did not differ significantly from those observed previously in Caucasians [
10,
14,
15].
There was no effect of analyzed polymorphisms on the risk of CAD in the present study. Because of the marginal frequency of the CA haplotype in our population, obtained results should be interpreted with great caution, however, rare haplotypes should not be ignored in the genetic studies of diseases with complex pathophysiology [
20]. Case-control association studies concerning the role of
USF1 gene polymorphisms on CAD risk are not numerous. In the study of Komulainen
et al. [
13] the effect of
USF1 polymorphisms on CAD risk and overall mortality was showed, but only in females. The effect concerned allelic variants of rs2073658 polymorphism, but not rs2516839 [
13]. In a study on U.S. whites, the frequencies of genotypes and alleles of rs3737787 SNP did not differ between CAD and non-CAD subjects [
11]. There are also few functional studies concerning the influence of
USF1 polymorphisms on atherosclerotic phenotype. In the autopsy series of 700 Finnish middle-aged men,
USF1 polymorphisms were analyzed with quantitative morphometric measurements of coronary atherosclerosis [
17]. The TT and CT genotypes of rs2516839 were associated with the proportion of advanced atherosclerotic plaques, area of calcification of coronary arteries and an increased risk of sudden cardiac death [
17]. Contrary to these results, in the study of Fan
et al. conducted on younger Finnish population, the C allele was associated with a higher mean carotid intima-media thickness, but only in females [
6]. These results are supported by the study carried out on the Chinese population of atherosclerotic stroke patients [
19], where the T allele carriers of rs2516839 SNP had a lower total unstable carotid plaque area than the CC homozygotes.
In the present work, the levels of triglycerides differed between the rs2516839 genotypes. A trend of increasing TG levels in relation to the C allele dose was observed, however, its statistical significance disappeared after the correction for multiple testing. Our findings are in agreement with previous studies, where C allele was associated with increased lipid levels [
6,
13,
14]. Although some studies imply that the T allele is rather related to the presence of lipid abnormalities [
7,
8] whereas other show no effect of rs2516839 SNP on lipid levels [
15,
19]. Similar differences concern the rs3737787 polymorphism [
8,
10,
11,
12,
18]. The presence of these discrepancies results from the fact that additional factors may modulate the functioning of both the
USF1 gene and the phenotype dependent on
USF1 polymorphisms. For the time being, sex-related allelic differences [
6,
13] and
USF1-other genes interactions have been documented [
7,
9].
In our study, a synergistic effect of cigarette smoking and rs2516839 on CAD risk was found. Cigarette smoking increased the risk of CAD especially in the C allele carriers. Our results also indicated that TG levels differentiated particular genotypes but only in smokers, independently of other traditional risk factors.
The molecular mechanism of interaction between the
USF1 gene and cigarette smoking remains unknown, but it may be partly explained by the fact that USF1 is considered as a stress-responsive transcription factor. This issue was comprehensively discussed in the review of Corre and Galibert [
5], summarizing that USF1 has antiproliferative properties and regulates many genes of cell cycle and tumor suppression (
i.e., cyclin B1 gene—
CCNB1, cyclin-dependent kinase 1 gene—
CDK1, breast cancer 2 gene—
BRCA2, tumor protein p53 gene—
TP53, adenomatous polyposis coli gene—
APC and others). Accordingly, a loss of transcriptional activity of
USF1 in many cancer cell lines was observed (discussed in Corre and Galibert [
5]). Interestingly, according to Fan
et al. study [
6], the
USF1 expression was significantly lower also in atherosclerotic plaque specimens compared with the control tissue [
6], and CC homozygotes of rs2516839 SNP had a lower expression of the
USF1 gene in relation to the T allele carriers. In this context, reduced expression of
USF1 in C allele carriers may impair the cellular response to cigarette smoking, and increase the predisposition to atherosclerotic CAD. The question whether the USF transcription factor is involved in the cellular response to smoking (nicotine or any other component of tobacco smoke) remains open, however, it was demonstrated that nicotine enhanced the USF1 translocation from the cytoplasm to the nucleus in odontoblast cells [
21].
There is, however, strong evidence for links between cigarette smoking and phenotype partly dependent on
USF1 gene polymorphisms [
22,
23]. It is well documented that cigarette smoking influences lipid concentrations. Results of a large meta-analysis of Craig
et al. [
22] indicated that cigarette smokers had significantly higher serum levels of TC (3.0%), TG (9.1%), VLDL-C (10.4%) and LDL-C (1.7%), and lower serum levels of HDL-C (−5.7%) compared with non-smokers, and the effect was dose-related. Cigarette smoking also influences lipid metabolism. The activity of lipoprotein lipase in skeletal muscles is reduced and this in turn decreases lipid clearance. Lipoprotein lipase activity is insulin-dependent and insulin resistance among smokers was observed in many studies [
24,
25].
In conclusion, C allele carriers of rs2516839 SNP seem to be particularly at risk of CAD, when exposed to cigarette smoking. Cigarette smoking influences TG levels, which increase with the number of copies of the C allele. The existence of numerous positive feedbacks between the USF1 gene, targets of USF transcription factor, cigarette smoking and lipid abnormalities may intensify the effects associated with each of the factors analyzed individually. In order to understand the relationship between USF1 gene polymorphisms, cigarette smoking and lipid levels on atherosclerotic phenotype, further functional studies are necessary.
4. Experimental Section
4.1. Subjects
Four hundred and seventy Polish Caucasians, inhabitants of Upper Silesia were studied. The first group (CAD) consisted of 235 patients with angiographically proven premature CAD (70 females and 165 males), aged 27–55 years (mean 44.61 ± 5.90). Group 2 (BD) included 235 blood donors matched by age and gender, with no signs of CAD and with negative familial history of CAD as inclusion criteria. CAD subjects were selected from patients admitted to the: (1) 1st Department and Clinic of Cardiology at the Upper Silesian Center of Cardiology in Katowice; (2) 1st Department of Cardiac Surgery at the Upper Silesian Center of Cardiology in Katowice. They were classified for the study by the same cardiologist. Controls were recruited from the Regional Centers of Blood Donation and Blood Treatment in Katowice and Raciborz. Following recommendations of the Polish Centers of Blood Donation and Blood Treatment, the blood samples were obtained only from subjects with systolic blood pressure <140 and diastolic blood pressure <90 on the day of blood collection. Inclusion and exclusion criteria, details of the medical interview, diagnosis and evaluation as well as criteria for CAD, myocardial infarction and traditional risk factors were described previously [
26]. The study protocol was approved by the Ethics Committee of the Medical University of Silesia in Katowice, Poland. All subjects gave written informed consents. The methods used in this study were in accordance with the Helsinki Declaration of 1975 and its further revisions.
4.2. Serum Lipids Measurement
Total serum cholesterol, HDL cholesterol and triglycerides were measured by enzymatic colorimetric methods (Analco, Warsaw, Poland). LDL cholesterol levels were calculated according to the Friedewald formula [
27] in subjects with triglycerides levels below 4.4 mmol/L.
4.3. DNA Extraction and Genotyping
Genomic DNA was extracted from peripheral leukocytes using the MasterPure genomic DNA purification kit (Epicentre Technologies, Madison, WI, USA). The USF1 gene polymorphisms were genotyped using the TaqMan® Pre-designed SNP Genotyping Assay (Applied Biosystems, Foster City, CA, USA). The total volume of 20 µL reaction mix included: 10 µL of TaqMan® Genotyping Master Mix (Cat. # 4371355), 1 µL of probe (TaqMan® Pre-designed SNP Genotyping Assay, Cat. # 4351376: ID; C_1839183_10 for rs2516839 or ID; C_1459759_20 for rs3737787 polymorphism), 1 µL of a DNA template (15 ng/µL) and 8 µL of deionized water. The probe was diluted with the TE buffer (1:1) before the reaction. Polymerase chain reaction amplification was performed according to the manufacturer’s specifications. Genotyping was performed using the 7300 Real-Time PCR System (Applied Biosystems). Correctness of genotyping of both polymorphisms was checked by re-genotyping 10%–15% of the samples. Repeatability of results was 100%.
4.4. Statistical Analysis
Data were analyzed using the Statistica 10.0 software (STATSOFT, Tulsa, OK, USA). Normality of distribution was assessed by the Shapiro-Wilk test and then a comparison of quantitative data was performed by the Mann-Whitney U test (for variables with non-normal distribution) or the student’s t test (for variables with normal distribution). The Hardy-Weinberg equilibrium was tested in all groups by a χ2 test as well as comparisons of genotypes and alleles frequencies between cases and controls. When the number of subjects in the subgroup was lower than 10 the Fisher’s correction was used. Odds ratios (OR) as well as their 95% confidence intervals (CI) were computed using an univariate analysis. Multiple logistic regression analysis was used after adjustment for age, sex and traditional risk factors of CAD (cigarette smoking, lipid levels, BMI value, hypertension, family history of CAD and diabetes mellitus). Risk ratio values (95% CI) were used when the number of individuals in any of the analyzed subgroups was 0. The effective sample size and statistical power of association analyzes were computed using Epi Info™ 7.1.1.0 developed by Centers for Disease Control and Prevention (Atlanta, GA, USA).
Haplotype blocks were defined by the HaploView software [
28] using the Four Gamete Rule method. D′ and r
2 values were used as linkage disequilibrium measures. A haplotype phase was determined on the basis of the population genetic data using the modified PHASE algorithm [
29] used in the SNPator script [
30].
To determine possible interactions between
USF1 genotypes and traditional risk factors of CAD, the 4 × 2 table approach was used. OR values obtained from 4 × 2 table comparisons were used for calculating multiplicative synergy indices (SIMs). SIM is a ratio of the observed effect with the joint exposure to genetic and traditional factors (OR
11) divided by the effect predicted for the joint exposure assuming multiplication of the effects observed in the presence of either a traditional (OR
01) or genetic factor (OR
10). This method was used only in the case when the number of subjects from each compared group (not exposed, exposed only to genetic factor, exposed only to traditional risk factor or exposed to both factors) was >0. The following formula of SIM was used [
31]:
Confidence intervals (at a 95% confidence level) for SIMs were calculated using the script described by Cortina-Borja
et al. [
32].
In univariate analyses statistical significance was accepted at p < 0.05. The p values were adjusted for multiple comparisons using the Bonferroni correction in the case of analyses of gene-traditional risk factors interactions and the effects of genotypes on serum lipid levels. Because each adjustment procedure included one SNP and one tested factor, the significance level was reduced to p < 0.025.
Because the statistics of biological interactions preclude the use of quantitative variables, an effect of USF1 gene variants and risk factors on CAD risk and lipid levels was tested using univariate and multiple logistic regression approach, after adjustment for age, gender and traditional risk factors of CAD.