Genetic mapping of a new heart rate QTL on chromosome 8 of spontaneously hypertensive rats

Background Tachycardia is commonly observed in hypertensive patients, predominantly mediated by regulatory mechanisms integrated within the autonomic nervous system. The genetic loci and genes associated with increased heart rate in hypertension, however, have not yet been identified. Methods An F2 intercross of Spontaneously Hypertensive Rats (SHR) × Brown Norway (BN) linkage analysis of quantitative trait loci mapping was utilized to identify candidate genes associated with an increased heart rate in arterial hypertension. Results Basal heart rate in SHR was higher compared to that of normotensive BN rats (365 ± 3 vs. 314 ± 6 bpm, p < 0.05 for SHR and BN, respectively). A total genome scan identified one quantitative trait locus in a 6.78 cM interval on rat chromosome 8 (8q22–q24) that was responsible for elevated heart rate. This interval contained 241 genes, of which 65 are known genes. Conclusion Our data suggest that an influential genetic region located on the rat chromosome 8 contributes to the regulation of heart rate. Candidate genes that have previously been associated with tachycardia and/or hypertension were found within this QTL, strengthening our hypothesis that these genes are, potentially, associated with the increase in heart rate in a hypertension rat model.


Background
Heart rate is primarily determined by regulatory mechanisms that are integrated within the autonomic nervous system. In arterial hypertension, an increased heart rate is commonly observed as a compensatory response to decreased stroke volume (due to either a diastolic or systolic myocardial dysfunction) [1,2]. Although an increased basal heart rate, accompanying myocardial dysfunction, might represent a secondary phenomenon in established hypertension, an early increase in heart rate may contribute directly to the pathogenesis of the disease [3][4][5][6]. In fact, tachycardia in arterial hypertension seems to be a heritable trait and independent of reflex responses [7]. Few studies, however, have focused on the identification of the genetic basis of tachycardia related to hypertension. A total genome scan to identify quantitative trait loci (QTL) provides a powerful tool to determine the chromosomal location of alleles and polygenes that contribute to complex diseases, such as autoimmune disorders [8,9], cardiovascular diseases [10,11], cancer [12,13], infectious diseases [14] and drug addiction [15,16].
Even though QTLs for complex diseases have been identified, little progress has been made towards the identification of the genes related to complex traits such as increased heart rate associated with hypertension. Recently, Jaworski et al. [17] reported a QTL in rat chromosome 2, centered around the D2Rat61/62 markers, that is related to the air-puff stimulus-induced bradycardia response. In agreement with Jaworski, Alemayehu et al. [18] described a heart rate QTL involved in bradycardia, induced by salt-load diet, on rat chromosome 2. We previously mapped five QTL's (two on chromosome 2, and one each on chromosomes 4, 8 and 16) involved in blood pressure variation in an F2 population of an intercross between Brown Norway (BN) and Spontaneously hypertensive rats (SHR) [19]. In this study, we used linkage analysis for QTL mapping and identification of candidate genes associated with basal levels of heart rate in the progenies of an F2 rat intercross between the BN and SHR strains.

Parental rat strain crosses and F2 phenotyping
To construct a detailed linkage map of all chromosomes and determine the genetic traits for tachycardia associated with hypertension we analyzed an intercross (F2) progeny (N = 188 rats) of the BN and the SHR strains, and normotensive (BN, N = 4) and hypertensive (SHR, N = 23) rats as parental strains. All animals were genotyped for 179 polymorphic markers distributed along the rat genome, according to a standard protocol [20]. Basal heart rate levels of the F2 progeny were obtained from a previously described phenotype database [19], and data from SHR and BN strains (Medical School, University of São Paulo) were evaluated in male 12-14 week old animals. Briefly, an arterial polyethylene catheter (PE-10, 0.28 mm ID, 0.61 mm OD, Biocorp Australia, Huntingdale, Victoria, Australia) was implanted under anesthesia (Pentobarbitone sodium, 30 mg/kg, i.p.) into the abdominal aorta through the left femoral artery at least 24 hours before measurements. Blood pressure was measured in nonanesthetized, freely moving rats, on a beat-to-beat basis at a 1.0 KHz of sample rate (DataQ Instruments, Inc. Ohio, USA) by connecting the arterial cannula to a pressure transducer (Statham P 23 Dd, Hato Rey, Puerto Rico). The pulse pressure signal from the transducer was fed to an amplifier (GPA-4 model 2, Stemtech, Inc., Wood Dale, IL, USA) and then to a 10-bit analog-to-digital converter (DataQ Instruments, Inc., Akron, OH, USA). Basal heart rate was derived from pulsatile blood pressure. All animal experimental procedures were followed in accordance to the guidelines for animal care and use of laboratory of the University of São Paulo, Brazil and the protocol approved by designated ethics committee.

Genetic marker map
The 188 animals from the F2 progeny were previously genotyped for 179 genetic markers distributed along the whole genome [19]. Although no new genotypes were incorporated for this study, we have recently re-accessed the genetic distances of each genotyped marker [21] through sequential consultation of different internet databases [22][23][24].

Quantitative trait loci (QTL) mapping
To assess the contribution of a putative QTL to the heart rate trait, we used the composite interval mapping method [25][26][27][28] in the F2 population (described above) using the Windows QTL Cartographer [29]. We then determined the significance of putative QTL's from the interval mapping scan based on a likelihood ratio (LR) threshold and considered 11.5 as suggestive linkage.

Criteria for candidate gene identification
Genes presented on the mapped heart rate QTL were listed by consulting the public database of the rat genome [30]. The selection of candidate genes proceeded according to the following criteria. Firstly, the known genes presented in the mapped QTL were listed; secondly, the known genes were submitted to gene ontology analysis (described below) for refinement of selection; and finally, genes were considered as candidate only if presented in a heart rate related phenotype, confirmed by citations on the NCBI PubMed browser [31].
Annotation of functional classes was carried out using information provided by the Gene Ontology database [32], accessed by NCBI's public repository [33]. The known genes found in the QTL were classified in specific function, process, or component classes through the GO terms assignation. The genes belonging to the same class were grouped and listed, and the candidate genes were elected from the most frequent gene ontology terms or according to gene ontology terms that are directly related to the heart rate trait.

Statistical analysis
Kolmogorov-Smirnov test (with Lilliefors' correction) was used to test data for normality of the estimated underlying population using SigmaStat 3.11 (Systat Software, Inc., San Jose, CA, USA). Normality was assumed for all parametric tests and regression procedures. All phenotype measurements (blood pressure and heart rate) were com-pared between each parental strain (SHR × BN) using the Student t test. Correlation between blood pressure and heart rate values was performed calculating the Pearson r correlation coefficient. Statistical significance was establish according to the P value (p < 0.05) for a two-tailed distribution. All data are reported as mean ± SEM.

Hemodynamic correlations
As expected, inbred parental strains had different levels of basal blood pressure and heart rate ( Table 1). Systolic (188 ± 2 vs. 127 ± 5 mm Hg), diastolic (138 ± 2 vs. 91 ± 4 mm Hg), and mean (162 ± 2 vs. 109 ± 5 mm Hg) blood pressure measurements were significantly increased in SHR compared to the normotensive BN strain. Moreover, the SHR strain presented an increase in heart rate levels (365 ± 3 vs. 314 ± 6 bpm, p < 0.0001) compared to BN. Although heart rate significantly correlated (p = 0.0069 for systolic blood pressure) with blood pressure levels ( Table 2), the R squared value suggested that these variables are poorly correlated only (0.0335 for systolic blood pressure). Thus, we did not adjust heart rate values for systolic or diastolic blood pressure traits before QTL mapping analysis.

QTL mapping
The linkage signals for the composite interval mapping based on 179 markers along the 21 rat chromosomes (Figure 1) revealed the presence of a putative QTL responsible for an increase in heart rate in chromosome 8 ( Figure 1). Figure 2 shows, in detail, the 6.78 cM heart rate QTL mapped on rat chromosome 8 (8q22-q24), flanked by APOA02 and R830 markers. In addition, the likelihood ratio statistics for the higher linkage signal of the putative QTL in chromosome 8 was 18.67. The estimated effects on heart rate are -27.3 bpm for the additive and 5.3 bpm for the dominant terms. Table 3 summarizes the allelic effects on heart rate levels for the APOA02, R1106, and R830 markers in homozygotes for either BN (BN/BN) and SHR (SHR/SHR), and heterozygotes (BN/SHR) in the F2 intercross (SHR × BN) population. Although no significant difference was observed in heart rate level for the APOA02 (p = 0.7614) and the R830 (p = 0.8963) markers, an allelic effect was observed for the R1106 marker (p = 0.0306). SHR homozygotes presented significantly lower heart rate values compared to BN homozygotes (377 ± 6 vs. 370 ± 3 vs. 354 ± 8 bpm, for BN/BN, BN/SHR, and SHR/SHR, respectively).

Discussion and conclusion
A number of reports studying hypertensive patients have demonstrated an increased basal heart rate accompanying myocardial dysfunction [1][2][3][4][5][6]. Furthermore, a positive association between heart rate and all-cause mortality or cardiovascular mortality has been observed [5,6]. As such, it is important to elucidate the genetic factors that lead to increased heart rate levels. Based on the observation of an increased basal heart rate in the SHR strain, and using a powerful genomic exploratory approach, we mapped a new QTL on the rat chromosome 8 that is associated with basal levels of heart rate in SHR and identified several genes that could potentially be responsible for heart rate modulation.
We, herein, show that the parental SHR strain displayed increased heart rate values compared to normotensive BN rats. Moreover, in the F2 intercross animals, the allelic Total genome scan analysis has contributed to the discovery of several QTLs associated with different phenotypes in complex diseases. However, although relative success has been achieved in the identification of monogenic diseases [34,35], limited progress has been made in the identification and confirmation of candidate genes in complex diseases. Different chromosomal regions, associated with heart rate-phenotype traits, have been mapped in both patients and experimental animals [7,17,18,36,37]. Kreutz et al. [36] mapped a QTL on rat chromosome 3 that was associated with both basal heart rate and saltinduced changes in heart rate in an F2 intercross (SHRSP × WKY) population. Jaworski et al. [17] identified independent chromosomal regions related to bradycardia in response to mild stress (airpuff startle). Alemayehu et al. [18] identified one putative heart rate locus on rat chro-mosome 2 and confirmed its presence in different congenic strains after the substitution of the mapped interval from the SHR for the allelic region of a normotensive strain. Interestingly, Kren et al [38] reported a blood pressure QTL on chromosome 8 which maps nearby the heart rate-related QTL reported in the present study. However, a congenic SHR-Lx strain harbouring this particular region showed only changes in basal blood pressure and cardiac mass with no heart rate effects. Therefore, it is less likely that the same gene(s) are influencing the effects on heart rate here described. In a mouse model, Sugiyama et al. [37] conducted a QTL mapping analysis in an F2 progeny (CBA/CaJ × BALB/cJ) and identified 3 loci associated with heart rate, two on the mouse chromosome 2, and one on chromosome 15. We, herein, mapped a new QTL located in chromosome 8 and associated with a negative effect on heart rate in the SHR. Kahn and collaborators [39] found the rat chromosome 2 to be involved in heart rate, but were not able to confirm their findings in a congenic strain after the substitution of the mapped interval from the SHR for the allelic region of a normotensive BN strain. With regard to the rat chromosome 8, we previously mapped a putative QTL [19], flanked by R19 and R850 markers, associated with salt-induced changes in blood pressure that is located relatively close to the QTL, demonstrated herein, that is associated with heart rate trait.
In humans, loci associated with the autonomic nervous system's control of the heart rate have also been mapped [7]. Singh et al. [7] investigated 725 patients from the Framingham Heart Study and mapped 2 QTLs (chromo-  somes 2 and 15) associated with the frequency domain of heart rate variability indexes. Although the authors [7] have not mapped a QTL associated specifically with basal heart rate, they proposed that the region contains genes that share function for both basal heart rate and heart rate variability traits. Interestingly, the QTL mapped on chromosome 15 in humans is a syntenic region compared with the interval we mapped for basal heart rate traits in the rat chromosome 8.
In the present investigation, interval mapping revealed 17 genes that are potentially implicated in basal heart rate in the SHR. Among these candidate genes, just 4 families of genes were identified that maintained the association with the heart rate phenotype after the use of certain selection criteria; these genes were Scn2b (Sodium channel type II, beta polypeptide); Htr3a and Htr3b (5-Hydrox-ytryptamine receptors type 3); Hcn4 (Hyperpolarizationactivated, cyclic nucleotide-gated K + type 4); Chrna3, Chrna5, and Chrnb4 (Neuronal nicotinic acetylcholine receptor, alpha polypeptide 3 and 5, and beta polypeptide 4). Interestingly, association studies and genetic engineered mice data support the involvement of some of the heart rate candidate genes in arterial hypertension.
Knockout mice for the hyperpolarization-activated, cyclic nucleotide -gated K + 4 (HCN4) gene have cardiac cells with "mature" pacemaker potentials, demonstrating that this gene that encodes a hyperpolarization-activated channel is essential for the proper generation of pacemaker potentials in the sinoatrial node [40]. In patients, a missense D553N mutation in the HCN4 gene was associated with sinus node dysfunction, QT prolongation in the electrocardiogram, polymorphic ventricular tachycardia, and Likelihood ratio (LR) statistics for basal heart rate trait for all chromosomes in 188 rats from an F2 intercross (SHR × Brown Norway) population Figure 1 Likelihood ratio (LR) statistics for basal heart rate trait for all chromosomes in 188 rats from an F2 intercross (SHR × Brown Norway) population. Each panel represents the chromosome number with the genetic length of the chromosome (cM) on the abscissa, and the horizontal dashed line denotes an LR score of 11.5 as the threshold for suggestive linkage.
syncope [41]. In addition, the presence of the S672R mutation in the coding region of the HCN4 gene, mimics a mild vagal stimulation in the sinoatrial node, slowing the heart rate by decreasing the inward diastolic current [42]. These data suggest that loss of function of the HCN4 gene is associated with sinus nodal dysfunction.
The neuronal nicotinic acetylcholine receptor (Chrn) subtypes α3, α5, and β4 are members of a ligand-gated ion channel family that affect sodium and potassium transients. Knockout mice strains lacking the genes encoding the subtypes α5 and β4 of the nicotinic receptors showed cardiac autonomic dysfunction, affecting the heart rate response to vagal electrical stimulation. In humans, although intragenic polymorphisms on the nicotinic receptor subunits have been identified [43][44][45], these gene mutations have not yet been associated with cardiovascular phenotypes.
One potential study limitation is the number of genetic markers used in the genome-mapping experiment.
Although a large number of F2 analyses have been reported, genes responsible for SHR hypertension and/or heart rate related phenotypes have not been yet identified. One of the possible explanations for this fact is the small number of markers genotyped in these studies. We used 179 informative markers in our mapping experiment using an intercross between the SHR and BN strains. Although genetic distances and locations of these markers were updated before conducting our analysis, it is possible that QTLs important for the phenotype studied remain unidentified.
Our data suggest that genetic markers on the rat chromosome 8 contribute to heart rate control in experimental hypertension. The recognition of genetic determinants of heart rate may help gain insights into the pathophysiology of the autonomic nervous system and offer clues as to its modulation in hypertension. Further studies will be necessary to identify the specific gene(s) in chromosome 8 that are involved in hypertension-related tachychardia and characterize the molecular mechanisms underlying the genetic influence on heart rate in hypertensive rats.
Rat chromosome 8 heart rate QTL in 188 rats from an F2 intercross (SHR × Brown Norway) population Figure 2 Rat chromosome 8 heart rate QTL in 188 rats from an F2 intercross (SHR × Brown Norway) population. The horizontal dashed line denotes the likelihood ratio (LR) of 11.5 as the threshold for suggestive linkage. See Table 3 for details of markers delimiting the QTL and the effects on heart rate, according to genotype.