Serotonin re-uptake transporter gene polymorphisms are associated with imatinib-induced diarrhoea in chronic myeloid leukaemia patients

Tyrosine kinase inhibitors (TKIs), the treatment of choice for chronic myeloid leukaemia (CML), can cause lower gastrointestinal (GI) toxicity which is manifested as diarrhoea. The mechanisms are not fully understood. The enteroendocrine signalling compound, serotonin (5-HT), is important for regulating peristaltic motion, fluid secretion and visceral hypersensitivity in the GI tract, and has been implicated in diseases such as irritable bowel syndrome. In this study, we have evaluated whether TKI-induced diarrhoea may be related to variation in the serotonin re-uptake transporter (SERT) gene. CML patients with and without diarrhoea on the SPIRIT2 trial (imatinib, n = 319; and dasatinib, n = 297) were genotyped for the promoter 5-HTTLPR, intron 2 VNTR and rs25531 polymorphisms by PCR-based methods. Diarrhoea was more prevalent in imatinib, than in dasatinib treated patients (P = 0.015), which when stratified by gender was seen to be driven by female patients (P = 0.036). Logistic regression analysis revealed that age, and the dominant HTTLPR with the rs25531 single nucleotide polymorphism (SNP) model, explained the occurrence of diarrhoea in ~10% of imatinib-treated female CML patients. These data suggest SERT polymorphisms influence imatinib-induced diarrhoea but not that of dasatinib.


Results patient characteristics.
Of the patients recruited to the SPIRIT2 trial, we excluded patients from this study for the following reasons: 1) no consent or sample given; 2) mention of constipation; 3) co-morbidities involving the GI tract; 4) documented co-medications that would influence diarrhoea or serotonin levels including selective serotonin re-uptake inhibitors; and 5) no PCR product due to poor DNA. This left 319 imatinib and 297 dasatinib-treated patients who were finally included in the study ( Supplementary Fig. 1). There were no significant differences between imatinib and dasatinib patient groups regarding age, BMI or platelet count (Table 1).

Allelic frequencies.
The rs25531 SNP and 5-HTTLPR promoter indel polymorphism (rs4795541) did not deviate from the Hardy-Weinberg Equilibrium (HWE) within the imatinib and dasatinib-treated CML cohorts ( Table 2). The tri-allelic STin2 VNTR polymorphism however showed significant deviation from HWE in both imatinib (χ 2 = 24.1537; P = 0.0002) and dasatinib (χ 2 = 19.0302; P = 0.0019) arms. Our observed minor allele frequency (MAF) for rs25531 was 7%, which is consistent with the reported Caucasian frequency of 9% on the HaploReg v4.1 website 23 . A study performed in 2622 European-American subjects, using the National Institute of Mental Health repository, found that the 5-HTTLPR short allele frequency was 43% 24 , which is identical to our findings. The frequencies of the STin2 VNTR 12, 10 and 9 alleles in our cohort were 62%, ~36% and ~2.5% respectively. Our frequency of 12 and 10 alleles are consistent with various addiction and disorder studies within a predominantly Caucasian background, but we have detected twice the frequency of 9 allele carriers than these published studies [25][26][27][28] . Our major deviation from the HWE was also replicated by these groups. We did not identify any carriers of the 7-repeat length STin2 VNTR in our patient cohort.
Logistic regression analysis incorporating all patients (n = 616) showed drug (P = 0.015), gender (P = 0.002), age (P = 0.049), BMI (P = 0.034) and the dominant HTTLPR model (P = 0.049) were significant covariates for the dependent variable of diarrhoea (Table 3). After retaining these covariates in a forward conditional regression analysis, only gender and BMI were retained in the model explaining diarrhoea as an adverse event in 3.4% of the patient cohort.
On further separation of the patients by the drug received, however, none of the tested covariates were significant for the dasatinib cohort. The significant differences were driven by the imatinib-treated patient cohort, and included gender (P = 0.007), 5-HTTLPR genotype (P = 0.020 and P = 0.021), especially of the dominant allele model (P = 0.012 and P = 0.0004), with or without the rs25531 SNP combined, respectively, and the STin2 VNTR (P = 0.043) polymorphism. After including these covariates in a forward conditional regression analysis, gender enhanced the 5-HTTLPR dominant model to explain diarrhoea in 7.8% of the imatinib-treated patients. By further stratifying the imatinib cohort by gender, only the 5-HTTLPR dominant model explained ~6% of male diarrhoea cases, but age and the 5-HTTLPR plus the rs25531 dominant model explained nearly 10% of diarrhoea cases for females (Table 4).
www.nature.com/scientificreports www.nature.com/scientificreports/ To ascertain which alleles influenced imatinib-induced diarrhoea incidence, we performed Chi-square analyses. The 'long' allele of the 5-HTTLPR genotype alone correlated with a greater incidence of diarrhoea (P = 0.001) ( Fig. 2A). The association was still significant when 5-HTTLPR was combined with the rs25531 SNP, and displayed a modest gene dose-dependent trend (P = 0.036) ( Fig. 2B) with the dominant 'long' allele driving this significance for the incidence of diarrhoea (P = 0.013) (Supplementary Table 3). When separated by gender, significant trends were not favoured by any one gender for the 5-HTTLPR dominant genotype (male: P = 0.012; female: P = 0.045) (Supplementary Table 3). Significance was lost when considering the 5-HTTLPR with rs25531 SNP combined, but a trend was observed favouring females with the dominant genotype (P = 0.066).
In patients with imatinib-induced diarrhoea, carriers of the 09 alleles for the STin2 VNTR were more likely to experience diarrhoea: model B (P = 0.032) (Fig. 2C) and recessive model A (P = 0.028) (Fig. 2D). When stratifying by gender this trend appeared to be influenced by females, though not significantly (Supplementary Table 3). However, these data must be viewed with caution owing to the low patient numbers of 09 allele carriers available for analysis: 9/186 (4.8%) male and 3/133 female (2.3%).
SeRt genotype does not correlate with diarrhoea grade. Diarrhoea toxicity grade was not significantly different between drug arms (P = 0.120; 95% CI = −0.020, 0.173) ( Table 1). Diarrhoea experienced by all TKI-treated patients were largely grade 1 and 2 (75% and 14% of patients for imatinib, and 71% and 19% for dasatinib, respectively) ( Table 1). Only three patients had grade 3 diarrhoea and there were no reports of grade 4 diarrhoea. Diarrhoea toxicity grade was not significantly different between drug arms or gender (Table 1). Chi-square analysis with genotype did not show any significant correlations with diarrhoea toxicity grade for either the imatinib or dasatinib arms (Supplementary Table 3).

Discussion
Our analysis shows that in the SPIRIT2 trial, the incidence of diarrhoea was greater in imatinib-than in dasatinib-treated chronic myeloid leukaemia (CML) patients. Imatinib-induced diarrhoea was influenced by gender and SERT genotype. Carriers of the 5-HTTLPR 'long' allele or the STin2 VNTR 9-repeat allele were more likely to experience diarrhoea. Logistic regression identified that the dominant 'long' allele of the 5-HTTLPR explained ~6% of diarrhoea incidence in men, while age and the dominant 'long' allele of the 5-HTTLPR plus rs25531 SNP explained 10% of cases of diarrhoea in females. Furthermore, females younger than 65 were less likely to experience diarrhoea than those over the age of 65 ( Supplementary Fig. 2). Although this trend was not Missing data 11 (9.7) 4 (7.3) 7 (11.9) 6 (7.6) 4 (8.9) 2 (5.9) Table 1. Patient clinical information stratified by drug and gender. A T-test of independent samples was performed to compare differences between drug groups. Text in bold indicates statistical significance (P < 0.05). BMI = body mass index. Prevalence of diarrhoea is classified as incidence (yes or no) and toxicity grade (1 to 4). (2020) 10:8394 | https://doi.org/10.1038/s41598-020-65350-0 www.nature.com/scientificreports www.nature.com/scientificreports/ significant, it is consistent with the findings with bosutinib where 92% over the age of 65 years experienced diarrhoea 6 . By contrast, dasatinib-induced diarrhoea was not correlated with any covariate tested.
To our knowledge this is the first study to investigate the role of SERT genotype in drug-induced diarrhoea. Our patients largely had grade 1 or 2 diarrhoea and we had no patients with grade 4 diarrhoea. However, grade 1 or 2 diarrhoeas (up to 1500 ml stool per day) can significantly impact a patient's quality of life. The SERT genotype has been investigated in the diarrhoea predominant irritable bowel syndrome (IBS-D), but studies have been contradictory in their assessment of which 5HTTLPR allele, 'long' or 'short' , is responsible for enhancing diarrhoea symptoms 29 . It is important to note that not all studies considered age or gender, or the rs25531 SNP that lies within the 5-HTTLPR region. Thus, many patients may have been incorrectly genotyped as 'long' when a G variant (L G ) would in fact have converted this to a 'short' allele phenotype. Our analyses demonstrated some loss of significance with the incorporation of the rs25531 SNP within the 5-HTTLPR genotype highlighting the subtle difference this may cause in changing phenotype and reported outcomes. Furthermore, a modest allele-dose-dependent pattern, with diarrhoea incidence, was observed with the incorporation of the rs25531 SNP, compared with the 5-HTTLPR indel alone, suggesting that combining the 5-HTTLPR indel and rs25531 SNP displays a more accurate dose dependent phenotype.
Less emphasis has been placed on the STin2 VNTR polymorphism with fewer published studies investigating this polymorphism compared with 5-HTTLPR. The shorter 10-repeat was shown to reduce SERT expression compared to the full length 12-repeat allele 14 . The shorter 9-repeat was not investigated, but it may be inferred by analogy that this allele would cause further reduced expression of SERT 30 . Owing to the low frequency of the 09 allele, previous studies may have been under-powered (~100 patients) to determine any association 19 . Our cohort, however, comprised ~300 patients per drug treatment which allowed us to show a correlation between carriers of the 09 allele and imatinib-induced diarrhoea. With reduced expression of SERT being associated with a shorter www.nature.com/scientificreports www.nature.com/scientificreports/ STin2 VNTR sequence, our data support the hypothesis that a reduction in SERT activity would result in elevated circulating 5-HT concentrations leading to diarrhoea.
The 5-HTTLPR and rs25531 polymorphisms within our CML patients conformed to HWE. The deviation of the STin2 VNTR polymorphism from HWE may be related to its complex nature; we identified patients with 12, 10 and 9 alleles in our cohort of 616 patients but did not identify an individual with the 7-repeat allele, and hence our sample size may not have been sufficient to determine the true population allele frequency of a polymorphism which consists of four, or possibly more, yet to be identified, alleles. Furthermore, we did not identify any significant differences in patient demographics, such as age, gender, BMI and platelet count, between the imatinib and dasatinib arms, and therefore any differences seen in diarrhoea AEs between the TKI arms may be attributed to the effect of the drug.
IBS-D is more prevalent in women than in men 31 , and the 5-HT 3 receptor antagonist alosetron is more effective against IBS-D in women than men 32 . It is therefore of interest that we found trends suggesting that the effect of the SERT genotype may be more important in females with imatinib-induced diarrhoea. Our study however would have benefited from the inclusion of more females to support, or refute, these trends into statistically significant findings.
We and other groups have only tested the three most commonly reported polymorphisms despite there being over 9700 documented SNPs within the SLC6A4/SERT gene 33 . Our a priori decision in not including further SNPs was due to the vast number and lack of functional information regarding all of these SNPs. Of course, the serotonin transporter is just one component in the serotonin pathway, which includes multiple 5-HT receptors (5-HTRs), the enzyme tryptophan hydroxylase 1 (TPH1) that is required to synthesise the supply of serotonin, as well as genes involved in post-receptor signalling cascades. Investigation of these SNPs may be worthwhile in the future, but will require a much larger sample size, and should probably proceed via a genome-wide association approach.
Our investigation was predicated on the hypothesis that TKI-induced diarrhoea was mediated through increased circulating 5-HT levels caused by SERT genotypes responsible for reducing SERT expression or activity. However, the associations we observed: (a) the 'long' allele 5-HTTLPR genotype, with or without rs25531, would reduce circulating 5-HT level, which is counterintuitive to our hypothesis; (b) carriers of the STin2 VNTR 09 allele, which reduces SERT expression, were more likely to experience diarrhoea, which does fit our hypothesis. It is clear that SERT functional activity is dependent on a number of factors, genetic and non-genetic, and is therefore highly complex. It is also possible that circulating 5HT levels may not necessarily correlate with 5HT concentrations in tissues, making the relationship with genotype even more complex. Nevertheless, the reason  for undertaking this study does have a strong biological plausibility, and the findings of the study suggest that further investigations of the role of 5HT in the pathogenesis of diarrhoea caused by TKIs, which is particularly troublesome for patients, is worthwhile. It may also allow the use of 5HT modulating drugs in the treatment of patients so that they can continue their therapy for CML. An interesting finding from our study is that the association was seen with diarrhoea induced by imatinib but not by dasatinib. Although dasatinib and imatinib are both classed as Bcr-abl inhibitors to treat CML, they vary with respect to the off-target tyrosine kinases that they inhibit 34 , and it is therefore possible that kinases within the bowel may be differentially affected by these drugs as off-target effects. An interaction between the 5HT system and the kinase may therefore be required to elicit diarrhoea. Since the mechanisms of TKI-induced diarrhoea have not been fully elucidated, this is an area which merits further investigation.
In conclusion, our study has identified an association between diarrhoea induced by imatinib, but not by dasatinib, and SERT polymorphisms. The strength of our study is that the patients were recruited from a randomised controlled trial (SPIRIT2) where patients will have had standardised assessment of efficacious and adverse events, as well as being randomised into different arms. However, despite the relatively large sample size, further studies will be needed to replicate this finding. It will also be important to understand the mechanistic basis of the role of 5HT in causing diarrhoea induced by TKIs, which should also include assessment of the off-target effects of the TKIs on the bowel kinome and the interaction with the 5HT pathway.
Methods patient recruitment. The SPIRIT2 trial is a phase 3 prospective randomised comparison of imatinib (400 mg) versus dasatinib (100 mg) each daily in newly-diagnosed chronic phase chronic myeloid leukaemia (CML) patients (www.ClinicalTrials.gov; study identifier: NCT01460693). Ethical permission for the trial was granted by the London Research Ethics Committee, UK (REC ref# 07/H0718/90). Written informed consent was obtained from patients prior to their inclusion and the study conducted according to the declaration of Helsinki. The principal results from the study have been presented previously 35 . The patients' ethnicity was not catalogued, but they were assumed largely to be Caucasians as they were recruited within the UK. Each arm of the study recruited 407 patients with the primary outcome measure described as 5-year event free survival. GI adverse events (AEs) were logged as occurrence of diarrhoea (i.e. yes or no) and grade was classified according to the Common Terminology Criteria for Adverse Events (CTCAE) system 36   SERT/SLC6A4 Genotyping. The SERT/SLC6A4 polymorphisms multiplex PCR assay was adapted from Wendland et al. 37 Primers were synthesised by MWG Eurofins, Germany (Supplementary Table 1). Patient genomic DNA (80 ng) was amplified using OneTaq enzyme (Promega, Southampton, UK) with GC Buffer (New England Biolabs, Herts, UK), dNTPs (10 mM final conc.), and equimolar concentrations (50pmol final concentration) of forward and reverse primers for both 5-HTTLPR and VNTR. PCR cycling conditions were: 94 °C for 30 secs; 35 cycles of 94 °C for 30 secs; 64 °C for 30 secs; 68 °C for 30 secs; and 68 °C for 5 mins. To determine presence of the rs25531 SNP within the 5-HTTLPR region, a restriction enzyme digest was performed on the PCR product. Seven microliters of product was digested with MspI in CutSmart buffer (New England Biolabs, Herts, UK) at 37 °C for 90 mins. PCR products were separated by 2.5% 0.5x TBE-buffered agarose gel and visualised with ethidium bromide using a UV trans-illuminator.
Interpretation of the PCR products expected is described in Supplementary Table 1. Ten percent of the samples were duplicated, for assay reliability, and genotypes were reproducible and in complete concordance. SERT genotypes, dominant and recessive models, were classified as shown in Supplementary Table 2.
Statistical analyses. Differences between treatment arms for age, gender, body mass index (BMI), platelet count and diarrhoea incidence were performed using a comparison of means test for independent samples. Correlations between genotype and categorical and ordinal outcomes were performed by Chi-square (χ 2 ) analyses, generating Pearson and Fisher's Exact tests. One-way ANOVA tests were performed with linear covariates. Logistic regression analyses were performed to describe the influence of covariates on diarrhoea. These analyses were performed using the IBM SPSS statistics software v24 (IBM Corp., Armonk, NY, USA).

Data availability
All data generated or analysed during this study are included in this published article (and its Supplementary  Information files).  Table 2. Actual patient numbers are incorporated within the histogram bars.