Fingolimod-induced decrease in heart rate may predict subsequent decreasing degree of lymphocytes
Low serum neurofilament light chain values identify optimal responders to dimethyl fumarate in multiple sclerosis treatment
Paulette Esperanza Walo-Delgado, Susana Sainz de la Maza, … Luisa María Villar
Peripheral proinflammatory Th1/Th17 immune cell shift is linked to disease severity in amyotrophic lateral sclerosis
Mengmeng Jin, Rene Günther, … Tjalf Ziemssen
Selective cognitive dysfunction and physical disability improvement after autologous hematopoietic stem cell transplantation in highly active multiple sclerosis
N. Giedraitiene, R. Kizlaitiene, … G. Kaubrys
Disease-modifying therapy aids cognition in multiple sclerosis
Maria Pia Amato & Lauren B. Krupp
Short-term changes in frequencies of circulating leukocytes associated with narrowband UVB phototherapy in people with clinically isolated syndrome
Stephanie Trend, Anderson P. Jones, … Prue H. Hart
Apparent changes in the epidemiology and severity of multiple sclerosis
Nils Koch-Henriksen & Melinda Magyari
The role of the immune system and the biomarker CD3 + CD4 + CD45RA−CD62L− in the pathophysiology of migraine
Zbyšek Pavelek, Ondřej Souček, … Martin Vališ
Epidemiology and treatment of multiple sclerosis in elderly populations
Caila B. Vaughn, Dejan Jakimovski, … Bianca Weinstock-Guttman
Patient's perspective in clinical practice to assess and predict disability in multiple sclerosis
S. Gil-Perotin, L. Bernad, … L. Solís-Tarazona
Tokunori Ikeda ORCID: orcid.org/0000-0002-5520-31851,2,
Tatsuyuki Kakuma3,
Mari Watari1 &
Yukio Ando1
Scientific Reports volume 8, Article number: 16430 (2018) Cite this article
An Author Correction to this article was published on 28 August 2019
Here, we determined whether degree of decreased heart rate due to fingolimod treatment correlates with decreasing degree of lymphocytes in relapse-remitting multiple sclerosis (RRMS). In total, 30 patients with RRMS were treated with 0.5 mg fingolimod and their heart rate recorded every 30 minutes for 24 hours. Time trends of heart rate were characterised as three individual amplitudes and phase angles from three cosine curves using a mixed-effect model. Spearman's correlation coefficient and regression analysis were used to determine the effect of heart rate information on change in lymphocyte count pre- and post-fingolimod treatment. Moreover, the degree of decreased lymphocytes induced by fingolimod treatment on heart rate was compared between low and high influence groups. Positive correlation between amplitude from the second curve and difference in lymphocyte number (p = 0.006) was observed. Regression analysis was also significant (p = 0.002). Moreover, the second curve derived from the high amplitude group exhibited a greater decrease in lymphocyte number after fingolimod treatment than the low amplitude group (p < 0.001). We suggest that the degree of decreased lymphocytes after fingolimod treatment (main effect) may be predicted by estimating the influence of degree in heart rate (side effect).
Fingolimod is an effective oral medicine for relapse-remitting multiple sclerosis (RRMS)1,2. Fingolimod is a structural analog of sphingosine and mainly acts as a functional antagonist at sphingosine-1-phosphate 1 (S1P1) receptors3,4. Fingolimod is phosphorylated in vivo and binds to S1P1 receptors on leukocytes, especially lymphocytes, which retains the lymphocytes in secondary lymphoid tissue5. Consequently, lymphopenia in peripheral blood is induced and maintained in a steady state 1 month later6,7,8. Moreover, fingolimod has other side effects, e.g., decreased heart rate, which can present because of S1P1 receptor activation on cardiomyocytes9. Therefore, with initial administration of fingolimod, the appearance of excessive decreases in heart rate or abnormal arrhythmia are monitored in patients at hospital admission. Thus, it is possible to obtain heart rate information within hours of fingolimod treatment.
In the present study, we assumed a similar mechanism for decreased heart rate and leukopenia, and determined whether the degree of decrease in heart rate correlates with decreasing number of leukocytes, with a specific focus on lymphocytes.
Statistical model of heart rate
To represent 24 hours of continuous heart rate after fingolimod treatment, we used a combination of three cosine curves. Our results showed correct expression of the statistical heart rate model (Fig. 1a,b). Next, heart rate curves for each patient were estimated using a mixed-effect model with three cosine curves and adjustments for age and sex. Accordingly, the shape of individual curves can be roughly classified into two groups based on whether the second wave is clearly visible or not (Fig. 1c,d).
Model of circadian heart rate rhythm using three cosine curves. Mean values (a) and median values (b) of heart rate were expressed using a combination of three cosine curves, as described in the Methods. Heart rate curves for each patient were estimated using a mixed-effect model adjusted for age and sex. Representative patients shown include those with a visible (c) and not visible (d) second wave. Each dot, dotted line, and solid line represents heart rate at a specific time point, with a connected line from actual heart rate and estimated line from three cosine curves model.
Correlation between heart rate and lymphocytes
For each individual, we obtained information in the form of three amplitudes and phase angles for individual predicted curves. We then examined Spearman's correlation coefficients between amplitude or phase angle and difference in leukocyte or subset (lymphocyte, monocyte, and neutrophil) count pre- and post-fingolimod treatment (Table 1). Strong correlation was detected between differences in leukocyte and lymphocyte (0.74, p < 0.001) or neutrophil (0.77, p < 0.001) count. In addition, moderate correlation was detected between difference in leukocyte and monocyte number (0.47, p = 0.009). Positive correlation was confirmed between amplitude from the second cosine curve and difference in lymphocyte number (r = 0.49, p = 0.006). While there was weak correlation with leukocytes (r = 0.36, p = 0.05), monocytes (r = 0.31, p = 0.09), but not neutrophils (r = 0.08, p = 0.69). In contrast, there was virtually no correlation between other heart rate information and difference in leukocyte or subset counts, except for correlation between phase angle from the first cosine curve and difference in monocyte number (−0.61, p < 0.001). Next, we used regression analysis (Fig. 2 and Table 2) to examine differences in values for leukocytes (p = 0.028, 95% confidence interval [95%CI]: 62.24–754.93) and lymphocytes (p = 0.002, 95%CI: 112.45–426.72) in relation to amplitude from the second curve. To confirm this regression model, we used bootstrap testing (n = 50,000), and corroboratively identified the same tendency (leukocytes, p < 0.001, 95%CI: 332.46–1043.49; lymphocytes, p = 0.004, 95%CI: 78.41–347.30). There was a weak relationship between difference in value for monocytes and amplitude from the second cosine curve (p = 0.12, 95%CI: −4.65–47.25). In this regard, bootstrap testing confirmed a positive relationship between difference in monocytes and amplitude from the second curve (p = 0.045, 95%CI: 1.80–54.20). A relationship with neutrophils was not confirmed by either test (regression, p = 0.39, 95%CI: −131.47–345.99; bootstrap, p = 0.38, 95%CI: −137.16–364.32).
Table 1 Spearman's correlation coefficients.
Linear regression plot. A linear regression plot showing the relationship between difference in number of pre-post leukocytes (a), lymphocytes (b), monocytes (c), and neutrophils (d) and amplitude from the second cosine curve.
Table 2 Results of multivariate regression analysis.
Comparison of difference in lymphocyte number in low and high amplitude groups
Underweight women are reported to be at risk for low lymphocytes10. Accordingly, we examined correlation and regression between body weight and lymphocyte number before fingolimod treatment. We found weak correlation (r = 0.39, p = 0.06) and a strong linear relationship (p = 0.002, 95%CI: 19.04–64.88) between these two variables in females (Fig. 3a) but not males (correlation, r = 0.09, p = 0.92; linear relationship, p = 0.97, 95%CI: −90.55–94.45) (Fig. 3b). Thus, our findings support the assumption that underweight women are in general at latent risk for low lymphocyte number. Finally, to discriminate between difference in lymphocyte number in individuals showing higher and lower effectiveness for fingolimod, we performed decision tree analysis using amplitude from the second cosine curve and a decided cut-off value. As shown in Fig. 3c, the high amplitude group (cut-off value ≥ 4.53) showed a greater decrease in lymphocyte number than the low amplitude group (p < 0.001, 95%CI: −1238.00–−317.58). We also examined body weight in the high amplitude group (n = 6) and found that all the patients were female but not necessarily the lightest in weight.
Linear regression plot showing the relationship between lymphocyte count and body weight (kg) in females (a) and males (b). Comparison in decreased lymphocyte number by classification of low and high amplitude groups. To distinguish patients with low or high sensitivity for fingolimod, a cut-off value was determined by decision tree analysis using the amplitude value from the second curve. Amplitude <4.53 was designated as the low amplitude group and ≥4.53 was the high amplitude group. Boxplot comparing the difference in lymphocyte number between both groups. Wilcox rank sum test was performed (c).
Because baseline lymphocyte count before fingolimod treatment was different in each individual, we next performed analysis of covariance (ANCOVA) and adjusted lymphocyte count number before fingolimod treatment between both groups. Accordingly, lymphocyte count number before fingolimod treatment was related to decreased lymphocyte number, as previously reported10 (Table 3). Moreover, our analysis showed decreased lymphocyte number in the high amplitude group compared with the low amplitude group (p = 0.015, 95%CI: 70.27–506.88) (Table 3).
Table 3 Results of analysis of covariance.
Lymphopenia after fingolimod treatment is the desired outcome for therapy of RRMS, yet it is possible that an excessive decrease of lymphocytes increases the risk of developing infection. Therefore, a marker is needed that can predict an excessive decrease of lymphocytes after fingolimod treatment. Previously, patients with a low lymphocyte number before fingolimod treatment and underweight women were prone to lymphopenia after fingolimod treatment11. In addition, underweight women are at risk for low lymphocytes10, yet the incidence rate of females with MS is higher than males12. Therefore, the risk of excessive lymphopenia must be considered when administering fingolimod to underweight women with RRMS.
Additionally, in this study, we focused on the functional mechanism of fingolimod. Fingolimod is a functional antagonist for four of five S1P receptors (namely S1P1, S1P3, S1P4, and S1P5)13,14,15. S1P1 and S1P3 receptors are expressed on cardiac myocytes, on which fingolimod promotes activation of G-protein-coupled inwardly rectifying potassium (GIRK) channels16,17. Although S1P3 receptors are involved in the fingolimod-induced decrease in heart rate, a novel S1P receptor modulator, siponimod (BAF312) (which influences S1P1 and S1P5 but not S1P3) was previously shown to induce a decrease in heart rate. Subsequently, considerable involvement of S1P1 receptors are recognised in this decreased heart rate18,19. In fact, S1P1 receptor expression levels are higher in cardiac myocytes than S1P3 receptors17. Further, half maximal effective concentration (EC50) in binding of S1P1 receptors and phosphorylated fingolimod is low compared with S1P3 receptors20. In contrast, no clinical adverse effects on cardiac rhythm have been reported with a selective S1P1 receptor modulator, amiselimod (MT-1303)21,22. Although the reason for the discrepancy in these results is not fully understood, amiselimod shows low activation of GIRK channels on cardiac myocytes compared with fingolimod22. Altogether, it is more likely that S1P1 receptors on cardiac myocytes are a key player in the decrease of heart rate, at least regarding fingolimod. Leukocytes are mainly composed of neutrophils, monocytes, and lymphocytes, and these cell types also express S1P receptor subtypes, including S1P1 receptors5,23,24,25,26,27,28,29. Decreased circulating lymphocytes and monocytes due to fingolimod or siponimod treatment involves S1P1 receptors23,25. However, the influence of fingolimod on circulating neutrophils is weak, and S1P4 and S1P1 receptors are the main S1P receptor subtypes in neutrophils27,28,29. Thus, S1P1 receptors are also involved in development of leukopenia via fingolimod.
Fingolimod binds S1P1 receptors, which are internalised through repetitive stimulation by fingolimod30. As a result, cardiac myocytes escape the influence of fingolimod to signal downstream, with subsequent early recovery of heart rate rhythm despite the presence of fingolimod. Without fingolimod treatment, S1P1 receptors on lymphocytes are naturally internalised by high concentration of physiological ligand (i.e., S1P) in blood and efferent lymph30. However, because internalisation of S1P1 receptors is reversible, expression of S1P1 receptors is restored in secondary lymphoid tissue (e.g., lymph nodes), in which S1P concentration is at a low density30,31. When lymphocytes migrate from secondary lymphoid tissue to efferent lymph, they utilise the density gradient of S1P between secondary lymphoid tissue (low concentration) and efferent lymph (high concentration)30. In contrast, with fingolimod treatment, fingolimod concentration is high in secondary lymphoid tissue, which induces S1P1 receptor internalisation on lymphocytes and interferes with lymphocyte migration into efferent lymph. Therefore, fingolimod retains lymphocytes in secondary lymphoid tissue and a decreased lymphocyte count persists in blood, unlike in cardiac myocytes32,33.
Nonetheless, there is a time lag between decreased heart rate and lymphocyte count, regardless of the mechanism induced by fingolimod via S1P1 receptors7. Our hypothesis for this is based on the finding that a decreasing heart rate by fingolimod is reportedly caused by sinus node cells in cardiac myocytes34. While total cardiac myocyte number is estimated at approximately 2–3 × 109 cells35, sinus node cell number is not known. However, the sinus node area is approximately 13.5 ± 2.5 mm in length, 1.2 ± 0.3 mm in height, and 5.6 ± 1.4 mm in width, with a low occupancy in heart tissue36. Therefore, it is highly likely that the number of sinus node cells is small. Furthermore, total cell number of lymphocytes is approximately 2 × 1012 cells in adults26. In addition, there are 400–600 lymph nodes (secondary lymphoid tissue) distributed throughout the body, and lymphocytes circulate within secondary lymphoid tissue, blood and efferent lymph, unlike sinus node cells which are fixed37. Because the body contains many secondary lymphoid tissues, the time for fingolimod distribution may be different in each secondary lymphoid tissue. In addition, the usual state of circulating lymphocytes and presence of larger number of lymphocytes may affect the probability that fingolimod encounters S1P1 receptors on lymphocytes. Therefore, although the decrease in peripheral lymphocyte count starts at an early stage, similar to the influence for sinus node cells (which are smaller in number and exist as a solid organ), it may take time to reach the steady state decrease in lymphocyte count induced by fingolimod treatment. Actually, lymphocyte count in blood gradually decreases in a day-dependent manner to reach a steady state38. In addition, decreased lymphocyte count in blood begins at 4–5 hours after fingolimod administration in mice4. Moreover, deceasing degree of pulse rate in a seated position and lymphocyte count after fingolimod administration shows dose-dependency7,8.
Finally, the mechanisms underlying a decrease in heart rate via fingolimod are assumed to be similar to leukopenia. Therefore, we determined whether the degree of decreased heart rate by fingolimod correlates with later decreasing degree of leukocytes or its subsets (including lymphocytes, monocytes, and neutrophils).
To validate this, we accurately aggregated continuous heart rate information among individuals, and then isolated shapes from certain variables of this data. In this respect, although the spline model can easily collate continuous heart rate information, it cannot appropriately isolate heart rate information as variables. Hence, the circadian rhythm of heart rate after fingolimod treatment was expressed by cosine curves using a mixed-effect model, with amplitude and phase angle for each of the three cosine curves obtained39,40. By this means, continuous heart rate information was summarised as minimal variables formed by three amplitudes and phase angles. In addition, we used a mixed-effect model to account for differences in medication time of fingolimod in each individual (which are evident as differences in behaviour or fluctuation for each heart rate). Usually a mixed-effect model analyses repeated measurement data (such as heart rate in our data), and contains random effects to consider variations in each subject41.
Hence, we confirmed whether any variables are linked to differences in number of leukocytes or subsets pre- and post-fingolimod treatment. Accordingly, we identified a positive relationship between amplitude from the second cosine curve and difference in lymphocyte number. A similar tendency was observed for leukocytes and monocytes but not neutrophils. Each of these relationships corresponded to the degree of S1P1 receptor expression in each subset, as described above. Altogether, these results suggest that amplitude from the second cosine curve is most important. Because the heart rate wave naturally presents as a monomodal wave, we assumed that the influence of fingolimod treatment could be extracted from the second wave in a fingolimod-influenced bimodal wave. Therefore, each amplitude and phase angle from the first and second cosine curves were mainly represented as first and second waves of the bimodal wave, respectively. Moreover, amplitude and phase angle from the third cosine curve acted as modulators of the bimodal wave. Actually, although we employed three cosine curve models, by reference to Akaike's information criterion and expectation, as modulators of the third cosine curve for the bimodal wave, amplitude from the second cosine curve in the two cosine curve model also manifested a similar outcome as the three cosine curve model. Hence, we interpreted the importance of amplitude from the second cosine curve based on the following rationale. When fingolimod strongly affects cardiac myocytes (sinus node cells), the fingolimod-related decrease in heart rate reaches its lowest point within 6 hours7, which is represented by a deep trough between the first and second waves of the bimodal wave. After that, because the body attempts to restore heart rate to its natural state, we inferred that amplitude from the second cosine curve increases depending on the trough depth, with a high amplitude reflecting the strength of fingolimod's effect. Certainly, in our data, the high amplitude group had a value of ≥4.53 in the second cosine curve, which was exhibited by a decreased difference in lymphocyte number pre- and post-fingolimod treatment. Moreover, we used 24, 12, and 8 hours in our model as each respective cycle of the three cosine curves. Accordingly, the second cosine curve has a cycle of 12 hours and half time of 6 hours, which is the lowest point of this cycle. Therefore, our statistical model suggests a decreasing heart rate peak within 6 hours of fingolimod treatment.
Otherwise, we also expected that fingolimod wholly impacted on the circadian rhythm of heart rate, and that the power of this influence is apparent from the first wave of the bimodal wave, namely, lower amplitude from the first cosine curve. However, this was not observed. In this regard, there was a negative relationship between phase angle from the first cosine curve and difference in monocyte number. Although the reason for this is not obvious, earlier appearance of the second wave in the bimodal wave (which indicates a strong effect of fingolimod) affected duration of the first wave. Consequently, a shorter phase angle might occur, and monocytes may be sensitive to this phenomenon.
In our analysis, six patients were treated with interferon (IFN)-β at induction of fingolimod treatment. Because IFN-β decreases lymphocyte number at considerable frequency42, lymphocyte count before fingolimod treatment may affect IFN-β in these patients. However, this situation is a usual occurrence in daily medical practice, and analyses must be performed by applying real-world situations. In addition, because we used ANCOVA as the statistical model for adjusting baseline lymphocyte count, the clinical state before fingolimod treatment was considered.
Above all, although our study has limitations regarding the number of RRMS patients, it is possible that degree of decreased lymphocytes after fingolimod treatment (main effect) may be predicted by estimating the influence of degree in amplitude from the second cosine curve (side effect). If so, by monitoring heart rate for 24 hours after fingolimod treatment, neurologists may be able to use amplitude information to detect an excessive influence of fingolimod on the main effect. Consequently, an earlier return visit may be considered for fingolimod-treated RRMS patients to check the degree of decreased lymphocytes.
S1P modulator drugs for MS that do not influence heart rate will be available in the future. However, this will take time; therefore MS patients should be treated as occasionally having an advantageous side effect without simply considering an adverse event.
In total, 30 patients with RRMS were evaluated. They were administered 0.5 mg fingolimod once daily at Kumamoto University Hospital (Kumamoto, Japan) between 2012 and 2016 (Table 4). Of these patients, fingolimod was taken around 9 am or 1 pm in 24 and 6 patients, respectively. Six patients switched from IFN-β therapy to fingolimod, while the other patients were not administered any disease-modifying drugs (including IFN-β) at induction of fingolimod. During initial treatment of fingolimod, all patients underwent 24 hours of continuous electrocardiogram (ECG) monitoring in hospital, with heart rate recorded every 30 minutes in the resting state. The number of leukocytes and its subsets (i.e., lymphocytes, monocytes, and neutrophils) were measured at a central clinical laboratory of Kumamoto University Hospital. This study was approved by the Institutional Review Board of Kumamoto University (Permit Number: 1279). Informed consent was obtained from all patients for participation in this study, and all experiments were performed in accordance with the Declaration of Helsinki.
Table 4 Characteristics in MS patients.
A retrospective longitudinal observational study was performed at Kumamoto University Hospital. Sample size was determined with consideration of the number of inpatients to Kumamoto University Hospital during the survey period. Shapiro–Wilk's test was performed to identify normal distribution of variables. Heart rate manifested as circadian rhythms, as described previously39. To investigate the relationship between 24 hours of heart rate monitoring data and leukocyte or subset (lymphocyte, monocyte, and neutrophil) count during the treatment period, two analytical steps were performed. Repeated measures of heart rate were treated as a predicting variable, while change in leukocyte or subset count was modelled as the response variable. In the first step, time trend of heart rate data was modelled as a circadian rhythm. Because circadian rhythms statistically express cosine curves40, amplitude and phase angle characterised the time trend of heart rate, and were predicted for each patient based on the mixed-effect model as follows:
$$\begin{array}{rcl}{y}_{ij} & = & {\beta }_{0}+{A}_{i}\,\cos \,(\tfrac{2\pi ({t}_{ij}-{\psi }_{1i})}{P})+{B}_{i}\,\cos \,(\tfrac{4\pi ({t}_{ij}-{\psi }_{2i})}{P})+{C}_{i}\,\cos \,(\tfrac{6\pi ({t}_{ij}-{\psi }_{3i})}{P})\\ & & +\,{\beta }_{1}{x}_{1i}+{\beta }_{2}{x}_{2i}+{\varepsilon }_{ij}\\ & = & {\beta }_{0}+{\alpha }_{i}\,\cos \,(\tfrac{2\pi {t}_{ij}}{P})+{\beta }_{i}\,\sin \,(\tfrac{2\pi {t}_{ij}}{P})+{\gamma }_{i}\,\sin \,(\tfrac{4\pi {t}_{ij}}{P})\\ & & +\,{\delta }_{i}\,\sin \,(\tfrac{4\pi {t}_{ij}}{P})+{\epsilon }_{i}\,\sin \,(\tfrac{6\pi {t}_{ij}}{P})+{\theta }_{i}\,\sin \,(\tfrac{6\pi {t}_{ij}}{P})\\ & & +\,{\beta }_{1}{x}_{1i}+{\beta }_{2}{x}_{2i}+{\varepsilon }_{ij}\\ {\alpha }_{i} & = & {A}_{i}\,\cos \,(\tfrac{2\pi {\psi }_{1i}}{P}),{\beta }_{i}={A}_{i}\,\sin \,(\tfrac{2\pi {\psi }_{1i}}{P}),{A}_{i}=\sqrt{{\alpha }_{i}^{2}+{\beta }_{i}^{2}},{\psi }_{1i}=\tfrac{P}{2\pi }\,\arctan \,(\tfrac{{\beta }_{i}}{{\alpha }_{i}})\\ {\gamma }_{i} & = & {B}_{i}\,\cos \,(\tfrac{4\pi {\psi }_{2i}}{P}),{\delta }_{i}={B}_{i}\,\sin \,(\tfrac{4\pi {\psi }_{2i}}{P}),{B}_{i}=\sqrt{{\gamma }_{i}^{2}+{\delta }_{i}^{2}},{\psi }_{2i}=\tfrac{P}{4\pi }\,\arctan \,(\tfrac{{\delta }_{i}}{{\gamma }_{i}})\\ {\epsilon }_{i} & = & {C}_{i}\,\cos \,(\tfrac{6\pi {\psi }_{3i}}{P}),{\theta }_{i}={C}_{i}\,\sin \,(\tfrac{6\pi {\psi }_{3i}}{P}),{C}_{i}=\sqrt{{\epsilon }_{i}^{2}+{\theta }_{i}^{2}},{\psi }_{3i}=\tfrac{P}{6\pi }\,\arctan \,(\tfrac{{\theta }_{i}}{{\epsilon }_{i}})\end{array}$$
yij: Value of heart rate at a time (hour) in each individual. β0: Common intercept. π: Circle ratio. P: 24 hours. x1i: Age (year) in each individual. x2i: Sex of each individual. Male was used as the base. \({\alpha }_{i},{\beta }_{i},{\gamma }_{i},\,{\delta }_{i},\,{{\epsilon }}_{i},{\theta }_{i}\): Total estimated coefficient value from fixed and random effects in the mixed-effect model. ψ1i: Individual phase angle of the first cosine curve at peak time in each individual. ψ2i: Individual phase angle of the second cosine curve at peak time in each individual. ψ3i: Individual phase angle of the third cosine curve at peak time in each individual. tij: Time point of day (hour) in each individual. Ai: Individual amplitude of first cosine curve. Bi: Individual amplitude of second cosine curve. Ci: Individual amplitude of third cosine curve. εij: Residual error (which followed a normal distribution).
In this model, age and sex affect heart rate; therefore these two variables were used as adjusted factors43. In addition, heart rate >90 bpm was regarded as an outlier (from the perspective of normal heart rate i.e., resting state) and was excluded from the model. In this mixed model, three types of cosine and sine curves showed a random effect. Three amplitudes (Ai, Bi, and Ci) and phase angles (ψ1i, ψ2i, and ψ3i) were calculated in each individual. In some evaluations, mean or median values from all MS patients were used to express the above three cosine curves. In the second step, Spearman's correlation coefficients and a regression model were used to examine the effect of heart rate on change in leukocyte or subset counts. The components of age and sex for heart rate information were included in the mixed-effect model, and were therefore not used in the regression model. Because change in leukocytes or subsets was stable after 1 month of fingolimod treatment7,8, difference in leukocyte and subset count was used. This was calculated by the difference in two mean values from 1 year to the period immediately before fingolimod treatment, and from 1 month to 1 year after fingolimod treatment. In addition, decision tree analysis was used to identify the cut-off value for determining whether amplitude was related to decreased lymphocyte number. The Wilcoxon rank sum test was performed to compare the difference in lymphocyte number in low and high amplitude groups. Difference in lymphocyte number in these two groups was performed by ANCOVA based on leukocyte count before fingolimod treatment. Analyses were performed using R version 3.3.2 (The R Foundation for Statistical Computing, Vienna, Austria) and SAS Version 9.4 (SAS Institute Inc., Cary, North Carolina, USA), with the level of statistical significance set at p < 0.05.
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
An amendment to this paper has been published and can be accessed via a link at the top of the paper.
Kappos, L. et al. A placebo-controlled trial of oral fingolimod in relapsing multiple sclerosis. N Engl J Med 362, 387–401, https://doi.org/10.1056/NEJMoa0909494 (2010).
Cohen, J. A. et al. Oral fingolimod or intramuscular interferon for relapsing multiple sclerosis. N Engl J Med 362, 402–415, https://doi.org/10.1056/NEJMoa0907839 (2010).
Mandala, S. et al. Alteration of lymphocyte trafficking by sphingosine-1-phosphate receptor agonists. Science (New York, N.Y.) 296, 346–349, https://doi.org/10.1126/science.1070238 (2002).
Pham, T. H., Okada, T., Matloubian, M., Lo, C. G. & Cyster, J. G. S1P1 receptor signaling overrides retention mediated by G alpha i-coupled receptors to promote T cell egress. Immunity 28, 122–133, https://doi.org/10.1016/j.immuni.2007.11.017 (2008).
Matloubian, M. et al. Lymphocyte egress from thymus and peripheral lymphoid organs is dependent on S1P receptor 1. Nature 427, 355–360, https://doi.org/10.1038/nature02284 (2004).
Article ADS CAS PubMed Google Scholar
Kowarik, M. C. et al. Differential effects of fingolimod (FTY720) on immune cells in the CSF and blood of patients with MS. Neurology 76, 1214–1221, https://doi.org/10.1212/WNL.0b013e3182143564 (2011).
Kira, J. et al. Fingolimod (FTY720) therapy in Japanese patients with relapsing multiple sclerosis over 12 months: results of a phase 2 observational extension. BMC Neurol 14, 21, https://doi.org/10.1186/1471-2377-14-21 (2014).
Saida, T. et al. Long-term efficacy and safety of fingolimod in Japanese patients with relapsing multiple sclerosis: 3-year results of the phase 2 extension study. BMC Neurol 17, 17, https://doi.org/10.1186/s12883-017-0794-5 (2017).
DiMarco, J. P. et al. First-dose effects of fingolimod: Pooled safety data from three phase 3 studies. Multiple sclerosis and related disorders 3, 629–638, https://doi.org/10.1016/j.msard.2014.05.005 (2014).
Nishida, T. & Sakakibara, H. Association between underweight and low lymphocyte count as an indicator of malnutrition in Japanese women. Journal of women's health (2002) 19, 1377–1383, https://doi.org/10.1089/jwh.2009.1857 (2010).
Warnke, C. et al. Initial lymphocyte count and low BMI may affect fingolimod-induced lymphopenia. Neurology 83, 2153–2157, https://doi.org/10.1212/wnl.0000000000001049 (2014).
Alonso, A. & Hernan, M. A. Temporal trends in the incidence of multiple sclerosis: a systematic review. Neurology 71, 129–135, https://doi.org/10.1212/01.wnl.0000316802.35974.34 (2008).
Alexander, S. P., Mathie, A. & Peters, J. A. Guide to Receptors and Channels (GRAC), 5th edition. Br J Pharmacol 164(Suppl 1), S1–324, https://doi.org/10.1111/j.1476-5381.2011.01649_1.x (2011).
Pyne, S. & Pyne, N. Sphingosine 1-phosphate signalling via the endothelial differentiation gene family of G-protein-coupled receptors. Pharmacology & therapeutics 88, 115–131 (2000).
Brinkmann, V. et al. The immune modulator FTY720 targets sphingosine 1-phosphate receptors. J Biol Chem 277, 21453–21457, https://doi.org/10.1074/jbc.C200176200 (2002).
Zhang, J. et al. Signals from type 1 sphingosine 1-phosphate receptors enhance adult mouse cardiac myocyte survival during hypoxia. Am J Physiol Heart Circ Physiol 293, H3150–3158, https://doi.org/10.1152/ajpheart.00587.2006 (2007).
Means, C. K. & Brown, J. H. Sphingosine-1-phosphate receptor signalling in the heart. Cardiovasc Res 82, 193–200, https://doi.org/10.1093/cvr/cvp086 (2009).
Gergely, P. et al. The selective sphingosine 1-phosphate receptor modulator BAF312 redirects lymphocyte distribution and has species-specific effects on heart rate. Br J Pharmacol 167, 1035–1047, https://doi.org/10.1111/j.1476-5381.2012.02061.x (2012).
Legangneux, E., Gardin, A. & Johns, D. Dose titration of BAF312 attenuates the initial heart rate reducing effect in healthy subjects. Br J Clin Pharmacol 75, 831–841, https://doi.org/10.1111/j.1365-2125.2012.04400.x (2013).
Albert, R. et al. Novel immunomodulator FTY720 is phosphorylated in rats and humans to form a single stereoisomer. Identification, chemical proof, and biological characterization of the biologically active species and its enantiomer. J Med Chem 48, 5373–5377, https://doi.org/10.1021/jm050242f (2005).
Harada, T. et al. Cardiac effects of amiselimod compared with fingolimod and placebo: results of a randomised, parallel-group, phase I study in healthy subjects. Br J Clin Pharmacol 83, 1011–1027, https://doi.org/10.1111/bcp.13203 (2017).
Kappos, L. et al. Safety and efficacy of amiselimod in relapsing multiple sclerosis (MOMENTUM): a randomised, double-blind, placebo-controlled phase 2 trial. Lancet Neurol 15, 1148–1159, https://doi.org/10.1016/s1474-4422(16)30192-2 (2016).
Durafourt, B. A. et al. Differential responses of human microglia and blood-derived myeloid cells to FTY720. J Neuroimmunol 230, 10–16, https://doi.org/10.1016/j.jneuroim.2010.08.006 (2011).
Lewis, N. D. et al. Circulating monocytes are reduced by sphingosine-1-phosphate receptor modulators independently of S1P3. J Immunol 190, 3533–3540, https://doi.org/10.4049/jimmunol.1201810 (2013).
Chun, J. & Hartung, H. P. Mechanism of action of oral fingolimod (FTY720) in multiple sclerosis. Clin Neuropharmacol 33, 91–101, https://doi.org/10.1097/WNF.0b013e3181cbf825 (2010).
Blum, K. S. & Pabst, R. Lymphocyte numbers and subsets in the human blood. Do they mirror the situation in all organs? Immunol Lett 108, 45–51, https://doi.org/10.1016/j.imlet.2006.10.009 (2007).
Gorlino, C. V. et al. Neutrophils exhibit differential requirements for homing molecules in their lymphatic and blood trafficking into draining lymph nodes. J Immunol 193, 1966–1974, https://doi.org/10.4049/jimmunol.1301791 (2014).
Sun, X. J. et al. Sphingosine-1-phosphate and its receptors in anti-neutrophil cytoplasmic antibody-associated vasculitis. Nephrol Dial Transplant 32, 1313–1322, https://doi.org/10.1093/ndt/gfw427 (2017).
Rahaman, M., Costello, R. W., Belmonte, K. E., Gendy, S. S. & Walsh, M. T. Neutrophil sphingosine 1-phosphate and lysophosphatidic acid receptors in pneumonia. Am J Respir Cell Mol Biol 34, 233–241, https://doi.org/10.1165/rcmb.2005-0126OC (2006).
Lo, C. G., Xu, Y., Proia, R. L. & Cyster, J. G. Cyclical modulation of sphingosine-1-phosphate receptor 1 surface expression during lymphocyte recirculation and relationship to lymphoid organ transit. The Journal of experimental medicine 201, 291–301, https://doi.org/10.1084/jem.20041509 (2005).
Schwab, S. R. et al. Lymphocyte sequestration through S1P lyase inhibition and disruption of S1P gradients. Science (New York, N.Y.) 309, 1735–1739, https://doi.org/10.1126/science.1113640 (2005).
Chiba, K. FTY720, a new class of immunomodulator, inhibits lymphocyte egress from secondary lymphoid tissues and thymus by agonistic activity at sphingosine 1-phosphate receptors. Pharmacology & therapeutics 108, 308–319, https://doi.org/10.1016/j.pharmthera.2005.05.002 (2005).
Chiba, K., Matsuyuki, H., Maeda, Y. & Sugahara, K. Role of sphingosine 1-phosphate receptor type 1 in lymphocyte egress from secondary lymphoid tissues and thymus. Cellular & molecular immunology 3, 11–19 (2006).
Meissner, A. & Limmroth, V. Update on the cardiovascular profile of fingolimod in the therapy of relapsing-remitting multiple sclerosis (MS). Multiple sclerosis and related disorders 8, 19–26, https://doi.org/10.1016/j.msard.2016.04.002 (2016).
Tirziu, D., Giordano, F. J. & Simons, M. Cell communications in the heart. Circulation 122, 928–937, https://doi.org/10.1161/circulationaha.108.847731 (2010).
Sanchez-Quintana, D. et al. Sinus node revisited in the era of electroanatomical mapping and catheter ablation. Heart (British Cardiac Society) 91, 189–194, https://doi.org/10.1136/hrt.2003.031542 (2005).
Hoorweg, K. & Cupedo, T. Development of human lymph nodes and Peyer's patches. Seminars in immunology 20, 164–170, https://doi.org/10.1016/j.smim.2008.02.003 (2008).
Kahan, B. D. et al. Pharmacodynamics, pharmacokinetics, and safety of multiple doses of FTY720 in stable renal transplant patients: a multicenter, randomized, placebo-controlled, phase I study. Transplantation 76, 1079–1084, https://doi.org/10.1097/01.tp.0000084822.01372.ac (2003).
Korpelainen, J. T., Sotaniemi, K. A., Huikuri, H. V. & Myllyla, V. V. Circadian rhythm of heart rate variability is reversibly abolished in ischemic stroke. Stroke 28, 2150–2154 (1997).
Tong, Y. L. Parameter estimation in studying circadian rhythms. Biometrics 32, 85–94 (1976).
Saville, B. R. & Herring, A. H. Testing random effects in the linear mixed model using approximate bayes factors. Biometrics 65, 369–376, https://doi.org/10.1111/j.1541-0420.2008.01107.x (2009).
Article MathSciNet PubMed PubMed Central MATH Google Scholar
Pohl, D., Rostasy, K., Gartner, J. & Hanefeld, F. Treatment of early onset multiple sclerosis with subcutaneous interferon beta-1a. Neurology 64, 888–890, https://doi.org/10.1212/01.wnl.0000153570.33845.6a (2005).
Mason, J. W. et al. Electrocardiographic reference ranges derived from 79,743 ambulatory subjects. J Electrocardiol 40, 228–234, https://doi.org/10.1016/j.jelectrocard.2006.09.003 (2007).
This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors. We thank Rachel James, Ph.D., from Edanz Group (www.edanzediting.com/ac) for editing a draft of this manuscript.
Department of Neurology, Graduate School of Medical Sciences, Kumamoto University, Kumamoto, Japan
Tokunori Ikeda, Mari Watari & Yukio Ando
Department of Clinical Investigation, Kumamoto University Hospital, Kumamoto, Japan
Tokunori Ikeda
Department of Biostatistics Center, Kurume University, Kurume, Japan
Tatsuyuki Kakuma
Mari Watari
Yukio Ando
T.I. and T.K. initiated and organised the study. T.I. and M.W. were involved in data collection. T.I. and T.K. analysed data and performed statistical analysis. T.I., T.K. and Y.A. contributed to data discussion. T.I., T.K. and Y.A. wrote the manuscript. All authors approved the final version of the manuscript.
Correspondence to Tokunori Ikeda or Tatsuyuki Kakuma.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Ikeda, T., Kakuma, T., Watari, M. et al. Fingolimod-induced decrease in heart rate may predict subsequent decreasing degree of lymphocytes. Sci Rep 8, 16430 (2018). https://doi.org/10.1038/s41598-018-34797-7
Fingolimod Treatment
Cosine Curve
Heart Rate Information
Group Amplitude
S1P1 Receptors
About Scientific Reports
Guide to referees
Journal highlights
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)