Influence of serotonin transporter gene polymorphism (5-HTTLPR polymorphism) on the relation between brain 5-HT transporter binding and heart rate corrected cardiac repolarization interval.

OBJECTIVE: Serotonin transporter gene polymorphism (5-HTTLPR polymorphism) predicts the degree of structural and functional connectivity in the brain, and less consistently the degree of vulnerability for anxiety and depressive disorders. It is less known how 5-HTTLPR polymorphism influences on the coupling between brain and neuronal cardiovascular control. The present study demonstrates the impact of 5-HTTLPR polymorphism on the relations between heart rate (HR) corrected cardiac repolarization interval (QTc interval) and the brain 5-HTT binding.
MATERIAL AND METHODS: Thirty healthy young adults (fifteen monozygotic twin pairs) (mean age 26±1.3 years, 16 females) were imagined with single-photon emission computed tomography (SPECT) using iodine-123 labeled 2β-carbomethoxy-3β-(4-iodophenyl) nortropane (nor-β-CIT). Continuous ECG recording was obtained from each participant at supine rest. Signal averaged QTc interval on continuous ECG was calculated and compared with the brain imaging results.
RESULTS: In the two groups [l homozygotes (n = 16, 10 females), s carriers (n = 14, 8 female)] HR and the length of QTc interval were not influenced by 5-HTTLPR polymorphism. There were no significant relations between HR and 5-HTT binding in the brain. There were significant associations between QTc interval and nor-β-CIT binding in the brain in l homozygotes, but not in s carriers (correlations for QTc interval and nor-β-CIT binding of striatum, thalamus and right temporal region were -0.8--0.9, (p<0.0005), respectively).
CONCLUSION: The finding of longer QTc interval with less 5-HTT binding availability in major serotonergic binding sites in l homozygotes, but not in s carriers, implicate to differentiated control of QTc interval by 5-HTTLPR polymorphism.

Introduction

Serotonin transporter (5-HTT) is coded by a single gene (SLC6A4), which is located in human chromosome 17q11.2 [1]. After the discovery of polymorphism of the promoter region of the 5-HTT gene (5-HTTLPR polymorphism), with longer allele (with 16 repeats) having higher basal and induced transcription rates than shorter (with 14 repeats), the evidence is still inconclusive about the direct impact of 5-HTTLPR polymorphism on the 5-HTT binding in the human brain. Both significant associations and findings with no association have been reported [2] [3] [4] [5] [6] [7] [8].

Instead of direct gene effect on the brain 5-HTT binding, certain patterns of coupling of central nervous functions and structures are more clearly determined by 5-HTTLPR polymorphism. For example, there may be higher structural covariance between amygdala and anterior cingulate in individuals with both long alleles (l homozygotes) than in s allele carriers, and also there may be higher baseline amygdala activity with exaggerated responses to stressful images in s carriers than in l homozygotes [9] [10]. It has been hypothesized that 5-HTTLPR polymorphism causes these structural and functional variations by influencing on the timing and duration of 5-HTT gene expression [8].

In the present study young healthy adults were examined to assess the impact of 5-HTTLPR polymorphism on the functional integration in terms of relations between the brain 5-HTT binding and cardiovascular function, i.e. heart rate (HR), and heart rate corrected QT interval (QTc interval). We were particularly interested on the relations between QTc interval and brain nor-β-CIT binding, as we noticed the association between QTc interval and striatum nor-β-CIT binding in partially the same, but larger sample, but this finding remained without satisfactorial explanation [11]. The present study reanalyzes this interrelation with genotype data, which with functional connectivity context brings essentially new light to understand the relation. QTc interval is also particularly interesting compared with other measurements from ECG, because QTc interval is partially controlled by autonomic nervous system [12] [13].As s allele in 5-HTTLPR polymorphism has been related to less intra brain functional connectivity, we hypothesised similarly divergent pattern of integration between QTc interval and the brain 5-HTT binding by 5-HTTLPR polymorphism.

Study Participants and Methods

Study participants

The sample of thirty healthy Caucasian young adults from fifteen MZ pairs (mean age 26.2±1.3, range 21–28 years; 16 women), discordant or concordant for body mass (mean BMI 26±4, range 19–32 kg/m2), were recruited from the FinnTwin16 study born from 1975 to 1979. They were previously studied to explore metabolic and neurobiological correlates of obesity [14]. None of the study participants had a history of psychiatric disorder, substance abuse or eating disorder based on a structured psychiatric interview. Somatic health was checked by clinical examination (K.P), and by laboratory tests. Six women in three pairs used oral contraceptives, two persons were smokers, and seven reported weekly use of alcohol [14]. The study was approved by the ethical board of the Research Ethics Committee of Hospital district of Northern Savo and written informed consent was obtained from all participants.

Genotyping

DNA was extracted from peripherally drawn blood samples from each study participants using standard methods. The SLC6A4 promoter region containing the long (16)/short (14) (l/s) polymorphism was PCR-amplified using the following primers: forward 5′-CGC TCC TGC ATC CCC CAT TA-3′ and reverse 5′-GGG ATG CGG GGG AAT ACT GGT-3′, which produced 297/253 bp (l/s) product. The genotype was analyzed in 3% MetaPhor (R) agarose (FMC BioProducts, Rockland, Maine) gel electrophoresis.

Imaging procedure and data analysis

Single-photon emission tomography (SPECT), using iodine-123 labelled 2β-carbomethoxy-3β -(4-iodophenyl) nortropane (nor-β-CIT), was used for visualization and quantification of 5-HTT binding in the brain. Nor-β-CIT binds to monoamine transporters, both to 5-HTT and dopamine transporter (DAT) proteins, with ten-fold and 30% higher affinity to 5-HTT than DAT sites in vitro and in vivo, respectively. [15] [16]. As monoamine transporter density varies regionally in the brain, nor-β-CIT binding reflects predominantly DAT binding in striatum, and 5-HTT binding in midbrain, thalamus and cortical regions [17] [18] [19] [16]. Reported repeatability for high 5-HTT regions was good with 4% test-to-retest mean variation [20].

Each twin pair was studied on the same day before noon. Mean activity of injected nor-β-CIT CIT (supplied by MAP Medical Technologies Oy, Tikkakoski, Finland) was 195±5 MBq. Serial SPECT scans (with position lasers head position control) were performed at 5 min, 6 and 24 hours after the tracer injection using dedicated Siemens MultiSPECT 3 gamma camera equipped with fan-beam collimators (Siemens Medical Systems; Hoffman Estates, Ill., USA). The SPECT scans were decay-corrected and reconstructed with Butterworth-filtered back projection in a 128×128 matrix with a pixel size of 3×3 mm, and were attenuation-corrected with a Chang's algorithm. The imaging resolution was 8–9 mm. The SPECT slices were consecutively summarised to the slice thickness of 6 mm and re-aligned using a semi-automatic brain quantification program of Siemens and the Talairach coordinates. The slices were rotated and re-aligned so that transaxial (x-direction), sagittal (y-direction) and coronal (z-direction) ones were at right angles to each other. Region of interest placement was based on a semi-automatic brain quantification program of Siemens. The lower threshold of 60% of the maximum count was used to reduce the volume averaging and partial volume errors. Equilibrium analysis method was applied to estimate specific binding ratio (target-cerebellum)/cerebellum [21]. Peak equilibrium time was assumed to be 24 hours in the striatum and 6 hours in the midbrain and thalamus, temporal and midfrontal regions [22]. Analyses were done by experienced analyser who was blinded to other information about the study participants.

ECG recording

After the brain imaging, continuous ECG recordings took place in a quiet dimly lit room during five minutes supine rest. AD -converted recordings were transferred into PC computer (digitalization with time resolution 200 Hz/channel) with WinCPRS program (Absolut Aliens, Turku, Finland). The Q wave was determined correspondingly as local minimum just before the R peak. The T-wave apex was determined as a local maximum between the S-peak and the P-wave. End of the T wave was determined as the crossing point of the isoelectric line and the tangent of the T-wave positioned at the half amplitude of the T-wave on the right side of the apex. The QT interval was determined as a time elapsed from the onset of Q wave to the end of the T wave, and steady state of 30–60 seconds of ECG recording was used for each sample of signal averaged QTc interval. For heart rate correction of QT interval Bazzet's, Friedricia's and Karjalainen approach were tested, and finally Karjalainen approach chosen [11]. (With Karjalanen approach QTc interval was independent on HR, which was not the case with Bazett's of Friedricia's corrections [23] [24].)

Statistical Analysis

Results are presented as means ± SD. Mann-Whitney two-tailed tests were used for group comparisons of the data, when indicated. Pearson and Spearman correlations were used, when indicated. Using Bonferroni correction P<0.008 was taken to be statistical significant with a two tailed test. To identify those variables which are likely associated to QTc interval variation we employed linear regression analysis and multiple regression model, cluster-corrected to take into account the sampling of twins within twin pairs and thus provide correct standard errors and p-values [25]. SPSS 19.0 (SPSS Inc, Chicago, U.S.A.) and Stata (StataCorp, Texas, U.S.A.) programs were used for statistical analysis of the data.

Results

Automatic detection of the QTc interval failed for three subjects (two l homozygotes and one s allele carrier). Basic characteristics of the study participants and nor-β-CIT binding values are presented in Tables 1 and 2. Means of body mass index, HR and QTc interval (n = 27) were not significantly different between l homozygotes and s carriers (Table 1). Mean thalamus 5-HTT binding was slightly lower in l homozygotes than in s carriers (Table 2).

Table 1

Basic characteristics in the two groups by 5-HTTLPR polymorphism.

	l homozygotes (n = 16)		s carriers (n = 14)
	mean	SD	mean	SD
Age (years)	26	1.5	26	1.1
BMI (kg/m2)	24.9	3.5	26.6	3.6
HR (bpm)	67	11	62	17
QTc (ms)	354	15	345	17

QTc interval was available in thirteen l homozygotes and twelve s carriers. BMI = body mass index. None of the differences between the two genotype groups were statistically significant (p>0.05).

Table 2

Brain nor-β-CIT binding in the two groups by 5-HTTLPR polymorphism.

		Striatum	Thalamus	Midbrain	Tem (R)	Tem (L)	Midfront
l homozygotes (16)	Mean	2.65	1.09*	1.29	0.27	0.26	0.34
	SD	0.33	0.19	0.12	0.06	0.05	0.1
s carriers (14)	Mean	2.75	1.25*	1.36	0.28	0.27	0.31
	SD	0.28	0.17	0.18	0.09	0.11	0.15

 = Thalamus 5-HTT binding was slightly lower in l homozygotes compared to s carriers (p = 0.006). Tem (L/R) = left/right temporal region.

QTc interval was related to regional 5-HTT binding in l homozygotes, as presented in Table 3. There were no significant relations between 5-HTT binding and QTc interval in s carriers (Table 3). In l homozygotes gender, body mass index and QTc interval accounted for 48% and 61% of the variability in thalamic and the right temporal 5-HTT binding; for QTc interval; p = 0.001 and 0.0005, respectively (n = 14).

Table 3

Correlations between heart rate correlated QTc interval and brain nor-β-CIT binding in l homozygotes and s carriers.

		Striatum	Thalamus	Midbrain	Tem (R)	Tem (L)	Midfront
l homozygotes	r	−0.84**	−0.82**	−0.13	−0.87**	−0.70	−0.44
	P -value	<0.0005	<0.0005	0.67	<0.0005	0.01	0.11
s carriers	r	−0.56	−0.24	−0.13	0.16	−0.25	0.04
	P -value	0.04	0.44	0.66	0.60	0.40	0.90

Correlation is significant at the level of 0.0005 (**) (two-tailed). QTc interval was available in fourteen l homozygotes and in thirteen s carriers, respectively.

Discussion

The main finding of the present study suggests impact of 5-HTTLPR polymorphism on the relation between QTc interval and the brain 5-HTT binding; while QTc interval was significantly related to thalamus and temporal 5-HTT binding in l homozygotes, no association was found in s carriers.

There is some evidence for associations between cardiac repolarization and the brain function, linking QTc interval to brain and efferent autonomic physiology. In a large series of patients with prior brain infarction, the location of right or left insula for brain infarction was associated with abnormal cardiac repolarization [26]. Stellatum blockade in healthy individuals was associated with the length of QTc interval [13].

To our knowledge, the association between the brain 5-HTT binding and QTc interval is a novel finding. With a limitation of small sample size this relation is based on high repeatability of both QTc interval and 5-HTT binding in the thalamus [27] [20]. The association is most likely related to autonomic control of cardiovascular system, rather than other factors which regulate QTc interval. Parasympathetic blockade with atropine induced attenuation of QT shortening in HIS bundle paced dogs [28]. In humans electric stimulation of vagus nerve via auricular nerve induced similarly QT and HR shortening [12]. One study showed the effect of QTc prolongation along with significant HR increase with systemic sympathetic stimulation [29]. More detailed study showed prolongation of QTc interval after right stellatum blocade, and the opposite effect, after the left side blockade [13].

The relation between the brain 5-HTT binding and QTc interval was dictated by 5-HTT genotype in the present study. The influence of s allele was dominantly inhibited gene transcription in vitro [3]. Imaging studies have not been able to show direct gene effect from 5-HTTLPR to the brain 5-HTT-binding [8]. Still, there is robust evidence for the significance of 5-HTTLPR on the brain structure and function: increased amygdala activity, increased amygdala stress responsiveness, reduced gray matter volume in perigenual anterior cingulate cortex, diminished corticolimbic structural co-variation, and increased amygdala response for fearful stimuli has been reported in s carriers in comparison with l homozygotes [30]. We hypothesised that there might be heart-brain connection in this genotype dependent coupling, and our findings may widen the picture of 5-HTTLPR dependent functional integrity to encompass cardiac-brain axis. The influence of s allele was absence of cardiac-brain coupling.

Limitations

This study is a preliminary evaluation with a small sample size (fourteen male and sixteen female participations) and generalizations from our findings should be taken with care.

The study participants were originally selected to evaluate metabolic neurobiological correlates of obesity, which is a limitation in principle, but adjustment for weight did not change the results.

Using bi-allelic gene characterization in this study may have influenced on the findings, but not likely towards the alternative hypothesis. The reported frequency of LGLA or LGLG genotype, obtained by tri-allelic gene characterization, is 10% in Finnish population, and it cannot be excluded that there was phenotypes with this genotype in the study participants [31].

We did not have superior accurate [11] C-DASB for 5-HTT imaging. Still, nor-β-CIT can be regarded as regionally specific tracer for 5-HTT imaging, and is has high reproducibility in high density regions of 5-HTTs, in the thalamus and in the midbrain [20].

Conclusions

Our results suggest that in healthy young adults, who were dichotomized by 5-HTTLPR, there was significant relation between QTc interval and the brain 5-HTT binding only in l homozygotes. This may reflect divergent constitutional functional integration of l homozyogtes and s allele carriers by 5-HTTLPR.