1. Introduction
Chronic heart failure (CHF) is the end result of various insults to the myocardium, the most common of which is ischemic heart disease [
1]. The prevalence of CHF is high in older individuals, and it is a major cause of morbidity, mortality, hospitalizations, and disability [
2]. Typically 1%–2% of the adult population in developed countries suffers from CHF, with the prevalence rising to ≥10% among persons aged ≥70 years [
3]. The pathogenesis of CHF and genetic factors underlying the predisposition to it remain largely unknown. Identification of the genes associated with CHF is important to improve the understanding of the etiology of this multifactorial disorder.
Adenosine is a naturally occurring endogenous purine nucleoside that is composed of an adenine molecule attached to a ribose sugar moiety (ribofuranose) via a β-N9-glycosidic bond (6-amino-9-β-
d-ribofuranosyl-9-H-purine) [
4]. Levels of adenosine increase in response to physiological and pathological stimuli. The resultant shifts in adenosine levels impact all major aspects of cardiovascular function, including the rate and strength of the heartbeat, conduction of the cardiac impulse, autonomic control of the heart, coronary perfusion, cardiovascular growth or remodeling, and cardiac and vascular resistance to injurious insult [
5].
Adenosine regulates cell function via four structurally related G-protein-coupled receptors: A
1 (ADORA1), A
2 (ADORA2), and A
3 (ADORA3). In addition, ADORA2 was subsequently sub-divided according to the existence of distinct high-(ADORA2A)
vs. low-affinity (ADORA2B) binding sites. Layland
et al. [
6] revealed that the use of adenosine for stress testing and induction of systemic (and coronary) hyperemia is primarily related to the activation of
ADORA2A and the resultant increase in myocardial blood flow. Furthermore, deletion of the
ADORA2A gene has been shown to reduce the heart rate [
7], and in a head-to-head comparison in gene-targeted mice for all four individual adenosine receptors, researchers observed that mice deficient in the gene encoding
ADORA2B are particularly prone to myocardial ischemia [
8,
9,
10,
11,
12,
13]. Overall,
ADORA2A and
ADORA2B play an important role in heart activity.
The activities of ADORA2A and ADORA2B could be affected by gene polymorphisms in vivo. It is therefore reasonable to hypothesize that single-nucleotide polymorphisms (SNPs) in ADORA2A and ADORA2B can influence the activity of the adenosine A2 receptors, and further affect the susceptibility and progression of CHF. To the best of our knowledge, this hypothesis has not been investigated previously. We therefore investigated whether four selected ADORA2A SNPs (rs2236625, rs2236624, rs4822489, and rs5751876) and one ADORA2B SNP (rs7208480) are associated with CHF risk or severity in a north Chinese population.
3. Discussion
This is the first report on the impact of ADORA2A and ADORA2B polymorphisms on CHF susceptibility and progression in Chinese patients. After adjustment for age, sex, and traditional cardiovascular risk factors, no significant differences were found between CHF patients and controls. In addition, the possible haplotype associated with CHF risk was also not found by the haplotype analysis of the ADORA2A gene. However, the rs4822489 SNP in ADORA2A was significantly associated with NYHA class in CHF patients; this SNP may therefore affect the progression of CHF.
Adenosine is a potent coronary vasodilator that is produced in cardiac myocytes via the catabolism of adenine nucleotides [
14]. It binds with four evolutionarily well-conserved receptors that are expressed ubiquitously [
15,
16]. Activation of
ADORA2A and
ADORA2B produces potent vasodilation of most vascular beds including the coronary circulation, resulting in an increase in myocardial blood flow [
17]. Furthermore,
ADORA2A indirectly alters cardiac contractility by modulating the ADORA1 antiadrenergic effect [
18]. The gene encoding
ADORA2A is located on chromosome 22q11.23. The rs4822489 SNP, for which the minor allele frequency is 0.467 in the Chinese Han population, is located in the intron areas of
ADORA2A, thus causing no change of encoding protein in theory. However, an intronic SNP might alter the splicing of primary transcripts or gene expression [
19]. Howe
et al. [
20] revealed that one specific subset that exerted distinct effects on a multifunctional intron retention reporter was likely preferred for the regulation of endogenous intron retention events. That is, intron retention could affect alternative splicing of pre-mRNA by influencing the exonic splicing silencers. Through alternative splicing, a gene could produce different mature mRNAs and then translated proteins with different functions, and in extreme cases with opposite functions [
20,
21,
22]. Therefore, we hypothesized that the intronic SNP rs4822489 could affect the expression of
ADORA2A by influencing the alternative splicing of its mRNA, leading to alteration of the number and function of
ADORA2A, which might in turn weaken or eliminate the cardioprotective effects of that receptor, including expanding the coronary artery and increasing cardiac contractility and output, to ultimately aggravate any heart damage and accelerate the progression of heart failure. Another possible explanation of the association between the intronic rs4822489 and NYHA class could be that this SNP is in linkage disequilibrium with the “causal” genetic variant.
It has been shown that the rs4822489 variant is associated with other diseases. Bashira
et al. [
23] explored the association between variants of
ADORA2A and proliferative diabetic retinopathy (PDR) in a cohort of patients with type 1 diabetes. They found that carriers homozygous for the T allele of SNP rs4822489 were less likely to progress to PDR during follow-up than were carriers of any G allele. Hider
et al. [
24] examined the role of
ADORA2A polymorphisms on the outcome of methotrexate (MTX) treatment in rheumatoid arthritis (RA), but they did not identify any association between rs4822489 and adverse events in RA patients for this SNP because of the genotyping success rate of <85%. Three other
ADORA2A SNPs were detected in the present study: rs2236624, rs2236625, and rs5751876. rs2236624 is an intronic SNP and is in linkage disequilibrium with 15 currently known SNPs, 4 of which are within exons but are synonymous coding SNPs; 5 are intronic, 6 lie in the 3' untranslated region, 1 is a nonsynonymous SNP, and another is a frame-shift polymorphism [
25]. Christine
et al. [
26] found that in both healthy people and patients with autism spectrum disorder, those with the rs2236624 CC genotype were more susceptible to anxiety and autism. Bashira
et al. [
23] concluded that type 1 diabetes patients with the rs2236624TT genotype were not susceptible to PDR. Hider
et al. [
24] considered that the rs2236624 T-allele was significantly associated with the adverse events of MTX treatment in RA. The rs5751876 SNP is a synonymous substitution located in a tyrosine codon of exon 2, which is reportedly associated with neuropsychiatric phenotypes [
27,
28,
29] and age at onset in patients with Huntington’s disease [
30,
31]. In addition, Xiaoyan
et al. [
32] found that the frequency of heterozygotes with the rs5751876 genotype was significantly lower in subjects with a high degree of myopia. The rs5751876 SNP is also reportedly associated with anxiety [
29,
33,
34]. However, there have been no relevant research reports for rs2236625.
ADORA2B is located on chromosome 21p12; the rs7208480 (intergenic) SNP of this gene was selected in this study to explore its effects on CHF. Desai
et al. [
35] concluded that rs7208480 was significantly associated with pulmonary hypertension (PH) and sickle cell disease. PH is common in patients with CHF and is associated with more severe symptoms and worse outcomes [
36]. It is also a strong and independent predictor of mortality among patients with heart failure and provides incremental and clinically relevant prognostic information independently of known predictors of outcomes [
37,
38,
39]. However, the present study found no association between the rs7208480 polymorphism and susceptibility to CHF and progression thereof.
In this study, we also found that the prevalence of the rs2236625 CC genotype and C allele were significantly higher in smokers among CHF patients, which suggested that this SNP might have a certain effect on the cardiovascular function of smokers. It is reported that adenosine could increase myocardial contractility, which might be due to increased coronary blood flow through activation of ADORA2A receptors (Gregg’s phenomenon) [
40,
41]. Smoking is known to be associated with the impairment of coronary endothelial function [
42,
43,
44]. Therefore, we hypothesized that rs2236625 of ADORA2A might be related to the influence of smoking on coronary blood flow. The hypothesis should be validation by further studies with larger patient samples.
Exploration of the HapMap database revealed that the minor allele frequencies of rs7208480 and rs5751879 in the Chinese Han population were both <0.05. The experimental results of this study show that subjects with the rs7208480TT genotype—the distribution of which in the control group deviated from the HWE—constituted only 0.023 in the entire study population. This may be due to the characteristics of genetic drift, since the subjects were from the same geographical area.
This study has produced interesting outcomes and the first evidence of a relationship between
ADORA2A polymorphism and the severity of CHF in a Chinese Han population; however, several limitations of the study need to be considered. First, the sample was relatively small, and so the reported results need to be validated in a larger sample; Second, the study subjects were drawn from a northern Chinese population, and similar studies in other populations were needed to confirm our findings; Third, we selected one SNP (rs7208480) that was found to be associated with PH by a literature search, without studying all the SNPs of
ADORA2B [
35]; Finally, the functional relevance of the identified polymorphisms in CHF remains to be determined.
4. Experimental Section
4.1. Patients and Controls
Chinese CHF patients (age > 18 years) were recruited in 2013 and 2014 from the First Affiliated Hospital, College of Medicine, Xi’an Jiaotong University. Patients were classified as having heart failure by experienced cardiologists using the Guidelines for Diagnosis and Treatment of Heart Failure in China (2013) and the cause of the disease was determined according to clinical assessment and echocardiography. The following patient exclusion criteria were applied: age ≤ 18 years, presence of severe hepatic or renal insufficiency, tumors or malignant disease, acute attack of CHF, severe acute infection, metabolic disorders, and acute myocardial infarction without revascularization within the previous 2 weeks. The sample collection was completed according to strict criteria by two visiting staff, and was reviewed by the Chief Physician of Cardiology. Eventually, 300 patients were enrolled.
The healthy subjects, from the Medical Examination Center of the same hospital (age > 18 years) and without known personal or family history of cardiovascular disease, were invited to participate in this study as controls. To identify symptom-free subjects and exclude suspected cardiovascular disease subjects, a detailed interview was conducted in the framework of a physical examination by expert physicians. Four hundred controls were finally included, comparable for sex and ethnicity.
To obtain data on demographic and clinical characteristics, such as gender, age, and traditional cardiovascular risk factors, interviewer-administered health-risk questionnaires were performed. Hypertensive patients were defined according to the Guidelines for Prevention and Treatment of Hypertension in China (2013). Someone taking antihypertensive drugs was also considered to have hypertension. In addition, dyslipidemia and diabetes were defined according to the Guidelines on Prevention and Treatment of Blood Lipid Abnormality in Chinese Adults (2007) and the Guidelines for Prevention and Treatment of Diabetes in China (2013), respectively. Information regarding smoking status was self-reported. Individuals were designated as “smokers” if they reported that they had smoked an average of at least five cigarettes daily for at least the past 12 months [
45].
The present study was performed with strict protocol under the Ethics Committee of the First Affiliated Hospital, College of Medicine, Xi’an Jiaotong University. A statement of informed consent was obtained from all participants after the procedure had been fully explained to them.
4.2. Genotyping Assay
The blood samples were collected into tubes containing ethylenediaminetetraacetic acid. After centrifugation, the samples were stored at −80 °C until analysis. The standard phenol-chloroform extraction method was used to extract genomic DNA from the whole blood. The concentration of DNA was measured by spectrometry (DU530 UV/VIS spectrophotometer, Beckman Instruments, Fullerton, CA, USA). Three tag-SNPs (
i.e., rs2236625, rs2236624, and rs4822489), which captured most of the known common variations of ADORA2A according to the Chinese Han population data from HapMap (available online:
http://www.hapmap.org), and one SNP (rs5751876) of ADORA2A yielded by a literature search were selected for the present study. In addition, we learned that rs7208480 of
ADORA2B was associated with PH, which was common in CHF patients, so this SNP was also studied in this study. Sequenom MassARRAY Assay Design 3.0 software (Sequenom Inc., San Diego, CA, USA) was used to design a Multiplexed SNP MassEXTEND assay [
46,
47,
48]. SNP genotyping was performed using the Sequenom MassARRAY RS1000 system according to the standard protocol recommended by the manufacturer [
48]. The corresponding primers used for each SNP in the present study are listed in
Table 8. Sequenom Typer 4.0 software (Sequenom Inc., San Diego, CA, USA) was used to perform data management and analyses [
46,
47,
48,
49]. The association between haplotypes and CHF risk was analyzed using the Haploview soft version 4.2 (Mark Daly’s laboratory, Massachusetts Institute of Technology/Harvard Broad Institute, Cambridge, MA, USA).
Table 8.
Primers used for this study.
Table 8.
Primers used for this study.
SNP_ID | Allele | 1st-PCRP | 2nd-PCRP | UEP-SEQ |
---|
rs2236625 | C/T | ACGTTGGATGTGGTTGGAGGAAATTAACGC | ACGTTGGATGATCTTGGTTAAGGCTACGAC | CTTTAGAGAAAAATCATGTTTTCC |
rs2236624 | T/C | ACGTTGGATGTCTGGAGGGTGGCTTTACTG | ACGTTGGATGGACAGGGTATGGAGTACAAG | ACAAGAGGCCTAGATCC |
rs4822489 | T/G | ACGTTGGATGCAGTATGGAAATCCCTGGTC | ACGTTGGATGACACGGGACTTTCTTTGCAG | GAATACTCCCCTTGTGGGTTCCC |
rs5751876 | T/C | ACGTTGGATGTTGGGCACTCCCTCCACTCA | ACGTTGGATGAGGAGTGTGGGCCAACGGCA | GCGGAGGCCCAATGGCTA |
rs7208480 | C/T | ACGTTGGATGGCAATGCAGAAGCACATTCG | ACGTTGGATGGTCAGTACTACTGTTAGTGG | GTGGGAGTAGCAATCTTGAA |
4.3. Statistical Analysis
SPSS 18.0 for Windows (PASW Statistics, SPSS, Chicago, IL, USA) was used for the statistical analyses. Deviation of genotype distributions from Hardy-Weinberg equilibrium (HWE) was determined using the χ2-test, which was also used to assess differences in genotype and allele frequencies between study groups. Logistic regression analyses with dominant, recessive, and additive genetic models were used to assess the associations between ADORA2A and ADORA2B SNPs and CHF.
The association between haplotypes and CHF risk was analyzed using the Haploview soft version 4.2 (Mark Daly’s laboratory, Massachusetts Institute of Technology/Harvard Broad Institute, Cambridge, MA, USA) and haplotype Stats software obtained from Catalan Institute of Oncology (Barcelona, Spain, available online:
http://bioinfo.iconcolgia.net).
Differences in parameters related to the severity of CHF, such as ejection fraction (EF), left-ventricular end-diastolic diameter (LVEDD), and left-ventricular end-systolic diameter (LVESD), were evaluated among patients, stratified according to the three different genotypes by analysis of variance, and significant differences among pairs of means were tested using Dunnett’s test. Differences in genotype distribution and allele frequency for the significant SNPs according to functional New York Heart Association (NYHA) class were also analyzed by χ2-test. Variables (traditional cardiovascular risk factors) with a retention p value of <0.10 in the univariate analysis were entered into a multivariable ordinal regression model to evaluate the association between SNPs and NYHA class. The cutoff for statistical significance was set at p < 0.05 (two-sided).