|
|
||||||||
Lipids, Lipoproteins, and Cardiovascular Risk Factors |
1 Unitat de Lípids i Epidemiologia Cardiovascular, Institut Municipal dInvestigació Mèdica (IMIM), Barcelona, Spain; 2 Universitat Autònoma de Barcelona, Barcelona, Spain; 3 Universitat Pompeu Fabra, Barcelona, Spain.
aAddress correspondence to this author at: Unitat de Lípids i Epidemiologia Cardiovascular, Institut Municipal dInvestigació Mèdica, c/ Dr Aiguader 88, Barcelona 08003, Spain. E-mail msenti{at}imim.es.
| Abstract |
|---|
|
|
|---|
Methods: A nested case-control design was used to analyze 3 polymorphisms of the ESR2 gene and their associated haplotypes in 710 myocardial infarction cases from the REGICOR (Registre Gironí del Corazón) study and 2379 controls randomly selected in a representative population of a Spanish cross-sectional study.
Results: The rs1271572 T allele was significantly more common in patients who developed MI (P < 0.001). No association was observed for rs1256049 or rs4986938. Assuming a dominant model of inheritance, the association, as determined by logistic multivariate regression after adjustment for conventional cardiac risk factors, remained statistically significant in men [odds ratio (OR) 1.65, 95% CI 1.18–2.30; P = 0.003) but not in women (P = 0.754). A very common haplotype encompassing the rs1271572 variant was also associated with the risk of MI in the overall population (OR 1.41, 95% CI 1.06–1.87; P = 0.020) and in men (OR 1.57, 95% CI 1.12–2.21; P = 0.009).
Conclusions: The rs1271572 SNP T variant was associated with increased risk of MI in a Spanish population, and this association was found to be limited to men only. Sex differences in the genetic risk merit further investigation.
| Introduction |
|---|
|
|
|---|
Extensive studies of the role of the ER
and ERβ receptors in cardiovascular disease have shown that both receptors are expressed in humans in vascular endothelial, smooth muscle, and myocardial cells (2). ERβ is expressed more highly in vascular smooth muscle cells in women than men, whereas ER
is present in equal quantities (3). In the human heart, ERβ is upregulated by pressure overload, and ERβ mRNA levels have been observed to be increased in hearts of patients with aortic stenosis, with the greatest increases seen in women (4). ER
has also been associated with the development of coronary artery disease in women (5).
Single nucleotide polymorphisms (SNPs) in both ERs have been associated with cardiovascular disease. ESR12 (estrogen receptor 1; formerly ESR) polymorphisms have been associated with coronary artery disease (6)(7)(8), myocardial infarction (9), and blood pressure (10). Polymorphisms in ESR2 [estrogen receptor 2 (ERβ)] have been associated with cardiovascular disease (11) in a North American population and left ventricular mass and left ventricular wall thickness in women with hypertension (12), blood pressure (13)(14), and premature coronary artery disease (15).
In animal models, ERβ knockout mice are characterized by right and left ventricular hypertrophy (16) and systemic hypertension (17). In females, ERβ has been shown to attenuate the hypertrophic response to pressure overload (18), play a role in the protection against ischemia/reperfusion injury in the heart (19) and against left ventricular hypertrophy (20), attenuate clinical and biochemical manifestations of heart failure, and maintain normal ventricular repolarization and ventricular automaticity after MI (21). In males, ERβ mRNA expression is induced after vascular injury in blood vessels (22), and ERβ seems to regulate the small artery tone more strongly than in females (23).
ESR2 was discovered in 1996 by Mosselman et al. (24). ERβ, the nuclear hormone receptor encoded by ESR2, is a 54.2-kDa protein comprising 485 amino acids. It has 86% structural homology with ER
in the DNA-binding domain and 58% in the ligand-binding domain (25). The gene is located on chromosome 14q22–24 and has 8 transcribed exons and a 5'-untranslated exon which exists in 2 isoforms (26). In the absence of hormone, ERβ is sequestered in a multiprotein inhibitory complex in the nuclei. The binding of a ligand causes activation and conformational change of the protein and homodimerization through phosphorylation. The dimers bind with high affinity to estrogen response elements (EREs), located in the regulatory regions of target genes, leading to regulatory functions. Various polymorphisms have been described in the nucleotide sequence of the ESR2 gene; one of them resulted in a nonsynonymous change from leucine to valine, but this polymorphism was not observed in our population (27)(28).
The aim of this study was to investigate the relationship between various SNPs in the human ESR2 gene and the risk of MI through a nested-case control study in a Spanish population.
| Materials and Methods |
|---|
|
|
|---|
We randomly selected 2379 control subjects [640 men and 1739 women, age 46.2 (12.9) years] from a representative population of a cross-sectional study designed to establish the prevalence of the major cardiovascular risk factors in the province of Gerona, Spain (29)(30). Control subjects were judged free of angina or MI by history, physical examination, electrocardiography, and routine laboratory data.
Written consent was obtained from all subjects, and the local ethics committee approved the study. Using a standardized questionnaire, we collected details on socioeconomic and demographic characteristics, in addition to information on smoking, hypertension, diabetes, and current medication use, from all subjects. Presence of a dyslipidemic state was determined by a self-reported history of increased lipid concentrations and/or use of measures controlling dyslipidemia (diet or lipid-lowering drug therapy). Presence of diabetes was determined by a self-reported history of diagnosis of diabetes and/or measures controlling diabetes (diet, oral antidiabetic medication, or injected insulin).
Hypertension was defined according to WHO criteria, and 542 participants were taking antihypertensive drug therapy. Blood pressure measurements were taken after a 5 min rest and 20 min later, and the value used was the mean of the 2 determinations. Standardized hypertension questionnaires were used as described (30).
Body mass index (BMI) was determined as weight divided by squared height (kg/m2). The mean BMI of the participants in the study was 27.2 kg/m2.
screening of esr2 gene
The promoter region and 9 exons encoding the ESR2 RNA were amplified from the genomic DNA of 23 participants. The PCR products were analyzed by direct sequencing (dideoxynucleotide-sequencing method ABI Prism BigDye Terminator 3.0; Applied Biosystems) and confirmed by sequencing of the second strand. The DNA samples of the participants were then analyzed by TaqMan assay (ABI Prism 7900HT; Applied Biosystems). Three interesting SNPs were identified in the ERβ gene (NCBI GeneID 2100):
G transition in exon 6, using the forward and reverse primers 5'-TGGCAGCCAAGCATCAACAT-3' and 5'-TGCTGCTGCACCACAGATTA-3' and the FAMAAGTGCGGCTCTT-TAMRA and VIC-CCAAGTACGGCTCTT-TAMRA probes.
A transition in the nontranslated region downstream of exon 9, using the forward and reverse primers 5'-TCGTCCTTGCCCTTGAGCCTAAAT-3' and 5'-ACGCTGCATTCAAATGTGCCCT-3' and the FAM-ACGCTTCAGCCTGTGA-TAMRA and VIC-CACGCTTCAGCTTGTGA-TAMRA probes.
T transition in the promoter region, using the forward and reverse primers 5'-ATT TGCCAGCGACACACTCT-3' and 5'-AGGCCTTTCGCGTTAGATCA-3' and the FAM-ATTGTGAGACCCCCC-TAMRA and VIC-TTGTGAGAACCCCC-TAMRA probes.
statistical analysis
We assessed deviation from Hardy-Weinberg equilibrium using a
2 test with 1 degree of freedom to compare the observed and expected genotype frequencies among the subjects. We used a
2 or Fisher exact test, as appropriate, to compare categorical variables between groups, and Student t test to compare continuous variables between groups. We estimated odds ratios (ORs) and 95% CIs for the effect of ESR2 polymorphisms on MI risk using logistic regression analyses adjusted for the effects of other cardiovascular risk factors. Logistic regression analysis was also used to test for interactions. Bonferroni correction was used to adjust for multiple comparison testing.
We quantified pairwise linkage disequilibrium (LD) between the polymorphisms using the Shi standardized coefficient D' (|D'|). We examined haplotype frequencies and haplotype-based analysis for confounding variables using a Bayesian algorithm implemented in the R software, available at http://www.R-project.org. The dominant model (where heterozygotes and homozygotes for a particular haplotype are assumed to have equivalent effects) was fitted by the function "haplo.glm" in the R package "haplo.stats." To deal with haplotype phase ambiguity, this function uses a method that performs an iterative 2-step EM (expectation maximization), with the posterior probabilities of pairs of haplotypes per subject used as weights to update the regression coefficients, and the regression coefficients used to update the posterior probabilities (31). P values <0.05 were considered statistically significant. We used the SPSS statistical software package version 15.0 for statistical analysis.
| Results |
|---|
|
|
|---|
|
esr2 genotypes, haplotypes, and allele frequencies
To search for sequence variants of the ESR2 gene, we sequenced the promoter, exons, and intron–exon junctions of the gene in 23 unrelated individuals. Of the 7 SNPs identified, 3 were selected and genotyped in our study population because they seemed interesting for their specific frequencies in this group of patients or their putative biological function, or because they were previously described in the literature.
The first SNP, rs1271572, is located in the promoter region, 839 base pairs from the first exon, and corresponds to a G-to-T transition. The second SNP, rs1256049, is situated in exon 6 and corresponds to a synonymous G-to-A transition. The third SNP, rs4986938, is located 39 base pairs downstream from the last exon 9, and corresponds to a G-to-A transition. Genotype frequencies did not differ significantly from those predicted by Hardy-Weinberg equilibrium for rs1271572 (P = 0.160) or for rs1256049 (P = 0.480), but they differed in the case of rs4986938 (P = 0.028).
Genotype frequencies did not differ between controls and MI subjects for the rs1256049 and rs4986938 SNPs, whereas for rs1271572, a higher percentage of MI patients carried the T allele compared with controls (Table 2
), indicating an association between this polymorphism and MI (P < 0.001). Frequencies of the rs1271572 SNP differed significantly between men and women (P = 0.004), but were not found to differ between controls and MI subjects (data not shown). The 3 SNPs were found to be in strong linkage disequilibrium. In particular, the rs1271572 polymorphism was in strong linkage disequilibrium with rs1256049 (|D' | = 0.959) and rs4986938 (|D'| = 0.828) and rs1256049 with rs4986938 (|D'| = 0.926). Five common haplotypes (frequency >2%) and 3 uncommon haplotypes were detected (Table 3
).
|
|
association between esr2 RS1271572 and mi
The association between the presence of T allele of the rs1271572 SNP and MI remained statistically significant in a multivariate analysis that adjusted for sex, age, BMI, diabetes, smoking, hypertension, and HDL cholesterol (P = 0.003) (Table 4
). When the data were stratified by sex, this association persisted in men (P = 0.003) but not in women (P = 0.754). The comparison between the 3 genotypes revealed that the most relevant difference of frequencies was between the heterozygote carriers and the GG carriers. Assessment of the other 2 SNPs did not reveal an association with MI.
|
association between esr2 haplotype and mi
The haplotype analysis revealed an association between MI and the very common haplotype TGG (defined as T allele for rs1271572, G allele for rs1256049, and G allele for rs4986938), which had an estimated frequency of 38.2% in the population. The association was statistically significant in both the overall population (P = 0.020) and in men considered separately (P = 0.009), with a relative risk of 1.41 (95% CI 1.06–1.87) and 1.57 (95% CI 1.12–2.21), respectively, for each copy of haplotype carried (Table 5
). The most common haplotype, with a frequency of 38.8%, corresponds to GGA.
|
| Discussion |
|---|
|
|
|---|
, there is accumulating evidence that ERβ may also play an important role in the development of cardiovascular disease. Recent studies have shown that ERβ knockout mice suffer from systemic hypertension and heart ventricular hypertrophy (32) but that 17-β-estradiol, through ERβ-mediated mechanisms, protects the murine heart against left ventricular hypertrophy (33). There is also growing evidence suggesting a role for ERs in regulation of vascular healing and proliferation after injury as well as ER-mediated regulation of endothelial-dependent vasodilator reactivity responses (34). In this study, we examined and observed an association between MI and ESR2 rs1271572 SNP in the overall population, an association that seems to involve men but not women. This association did not change after adjustment for clinically relevant cardiovascular risk factors, suggesting that rs1271572 SNP is an independent risk factor for MI risk.
This SNP is located in the promoter region of ESR2 but does not belong to any consensus regulatory sequence. It is possible that the translation to RNA might be affected if the SNP belongs to a site not yet described, such as transcription factor binding sequences. The observed association of the SNP with MI may also be related to a change in the structure of either the DNA or RNA molecules, e.g., a change in the ability of the DNA to wrap itself into a nucleosome around histones, thus affecting the access of regulatory proteins to their binding sites (35). A third possibility is that the SNP could be in linkage disequilibrium with other functional modifications.
The statistical significance of these results persisted after correction for multiple testing by the Bonferroni method. Therefore, it is extremely unlikely that these results are false positives, although the possibility cannot be totally excluded. The use of the Bonferroni method to adjust for multiple testing is controversial in genetics, because it assumes that there is total independence between different polymorphisms. Many polymorphisms are known not to be independent because their variation in the population can be structured into haplotypes, transmitted as blocks of alleles in linkage disequilibrium with each other. Haplotype analysis reduces the dimensionality of association tests and increases the statistical power (36). The present study revealed an increased risk of MI in men carrying the haplotype named TGG. The risk was found also in the overall population but did not appear in women, which could be due to the small number of female cases available in the study. The TGG haplotype represents a unique change in the first SNP, rs1271572, from the major G allele to the minor T allele. The other SNPs involved are represented in the haplotype by their major allele. Therefore, the rs1271572 SNP could be a genetic risk factor for MI.
A few hypotheses may explain the differences in risk between men and women that are conferred by this haplotype. First, ERβ is more highly expressed in women than in men; thus, depending on the number of copies of the haplotype, the increased quantity of normal ERβ protein produced in women could compensate for the defect. Another possibility is that not only the expression and regulation of ERβ, but also its activity (as determined by a complex process that is modulated by estrogens and antiestrogens), could be influenced by the differences in concentrations of circulating hormones between the sexes. Interactions with other genes or nongenetic factors correlated with sex might also explain this phenomenon (37). ERβ may also mediate sex differences in ischemia-reperfusion injury, as ERβ knockout female mice display significantly less functional recovery (and more necrosis) than wild-type female control mice; thus ERβ may have a greater cardioprotective role in females (38).
A study on the same set of 3 polymorphisms in the ESR2 gene was recently reported by Rexrode et al. (11). In that study, an association between rs1271572 and MI was also described, but the association was seen only in women, whereas in our study it was found only in men. The contradictory findings between the Rexrode study and ours might be due to differences between populations. Patients of the Rexrode study were mainly whites from North America, and the frequencies described in the HapMap project for the rs1271572 polymorphism are very different in African and Asian populations. Therefore, it is highly possible that even a small proportion of people from ethnic backgrounds other than white could lead to some differences in observed frequencies. Rexrode et al. also described an association between the rs1256049 SNP and MI in women (11), but this association was not observed in our population. The differences between these 2 studies on the same subject suggest that a more detailed examination of the ESR2 gene will be important.
Techniques available today permit analysis of the entire ESR2 gene using SNP tagging, genome-wide associations, copy number variants, microsatellite analysis, and many other approaches. It would be useful to investigate whether the SNPs considered in this study are in linkage disequilibrium with other functional variants in the ESR2 gene or other genes in close proximity on the chromosome, such as KCNH5 [potassium voltage-gated channel, subfamily H (eag-related), member 5], and other functionally-related genes such as ER
. Functional studies should be performed to analyze the levels of RNA and protein expression and their activity.
A recent genome-wide association study with various cardiovascular disease outcomes in the context of the Framingham Heart Study (39) did not identify ESR2 as a candidate gene. Instead, these investigators identified 4 SNPs located on the same chromosome, 14. Three of the SNPs were associated with atrial fibrillation and 1 with major cardiovascular disease. In particular, the gene GJA4 (gap junction protein,
4, 37 kDa) was described as associated with cardiovascular and coronary heart disease outcomes. In this regard, a genome-wide association study in our population would be necessary to confirm and support those results.
In summary, the present study is the largest reported investigation of the relationship of ESR2 genetic variants and MI. The findings of this study suggest that a variant of ESR2 may contribute to cardiovascular disease; this knowledge may be helpful in elucidating the complex mechanisms of cardiovascular diseases and their genetic components. Functional studies, as well as the analysis of nearby polymorphisms that could explain the phenotype observed, may help in understanding the exact contribution of the variant described to the susceptibility to cardiovascular disease.
| Acknowledgments |
|---|
Financial Disclosures: None declared.
| Footnotes |
|---|
2 Human genes: ESR1, estrogen receptor 1 (formerly ESR); ESR2, estrogen receptor 2 (ER β); KCNH5, potassium voltage-gated channel, subfamily H (eag-related), member 5; GJA4, gap junction protein,
4, 37 kDa (formerly CX37). ![]()
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |