Hostname: page-component-7bb8b95d7b-pwrkn Total loading time: 0 Render date: 2024-10-06T09:44:56.008Z Has data issue: false hasContentIssue false

Longitudinal trajectory analysis of antipsychotic response in patients with schizophrenia: 6-week, randomised, open-label, multicentre clinical trial

Published online by Cambridge University Press:  22 October 2020

Minhan Dai
Affiliation:
Mental Health Center and Psychiatric Laboratory, West China Hospital of Sichuan University, China; and West China Brain Research Center, West China Hospital of Sichuan University, China
Yulu Wu
Affiliation:
Mental Health Center and Psychiatric Laboratory, West China Hospital of Sichuan University, China; and West China Brain Research Center, West China Hospital of Sichuan University, China
Yiguo Tang
Affiliation:
Mental Health Center and Psychiatric Laboratory, West China Hospital of Sichuan University, China; and West China Brain Research Center, West China Hospital of Sichuan University, China
Weihua Yue
Affiliation:
Peking University Sixth Hospital (Institute of Mental Health), China; and National Clinical Research Center for Mental Disorders and Key Laboratory of Mental Health, Ministry of Health (Peking University), China
Hao Yan
Affiliation:
Peking University Sixth Hospital (Institute of Mental Health), China; and National Clinical Research Center for Mental Disorders and Key Laboratory of Mental Health, Ministry of Health (Peking University), China
Yamin Zhang
Affiliation:
Mental Health Center and Psychiatric Laboratory, West China Hospital of Sichuan University, China; West China Brain Research Center, West China Hospital of Sichuan University, China
Liwen Tan
Affiliation:
Second Xiangya Hospital, Central South University, China
Wei Deng
Affiliation:
Mental Health Center and Psychiatric Laboratory, West China Hospital of Sichuan University, China; and West China Brain Research Center, West China Hospital of Sichuan University, China
Qi Chen
Affiliation:
Second Xiangya Hospital, Central South University, China
Guigang Yang
Affiliation:
Beijing Anding Hospital, Institute for Brain Disorders, Capital Medical University, China
Tianlan Lu
Affiliation:
Peking University Sixth Hospital (Institute of Mental Health), China; and National Clinical Research Center for Mental Disorders and Key Laboratory of Mental Health, Ministry of Health (Peking University), China
Lifang Wang
Affiliation:
Peking University Sixth Hospital (Institute of Mental Health), China; and National Clinical Research Center for Mental Disorders and Key Laboratory of Mental Health, Ministry of Health (Peking University), China
Fude Yang
Affiliation:
Beijing HuiLongGuan Hospital, China
Fuquan Zhang
Affiliation:
Wuxi Mental Health Center, Nanjing Medical University, China
Jianli Yang
Affiliation:
Institute of Mental Health, Tianjin Anding Hospital, China; and Tianjin Medical University General Hospital, Tianjin Medical University, China
Keqing Li
Affiliation:
Hebei Mental Health Center, China
Luxian Lv
Affiliation:
Second Affiliated Hospital of Xinxiang Medical University, China
Qingrong Tan
Affiliation:
Department of Psychiatry, Xijing Hospital, Fourth Military Medical University, China
Hongyan Zhang
Affiliation:
Wuxi Mental Health Center, Nanjing Medical University, China
Xin Ma
Affiliation:
Beijing Anding Hospital, Institute for Brain Disorders, Capital Medical University, China
Lingjiang Li
Affiliation:
Second Xiangya Hospital, Central South University, China
Chuanyue Wang
Affiliation:
Beijing Anding Hospital, Institute for Brain Disorders, Capital Medical University, Beijing, China
Xiaohong Ma
Affiliation:
Mental Health Center and Psychiatric Laboratory, West China Hospital of Sichuan University, China; and West China Brain Research Center, West China Hospital of Sichuan University, China
Dai Zhang
Affiliation:
Peking University Sixth Hospital (Institute of Mental Health), China; and National Clinical Research Center for Mental Disorders and Key Laboratory of Mental Health, Ministry of Health (Peking University), China
Hao Yu
Affiliation:
Department of Psychiatry, Jining Medical University, China
Liansheng Zhao
Affiliation:
Mental Health Center and Psychiatric Laboratory, West China Hospital of Sichuan University, China; and West China Brain Research Center, West China Hospital of Sichuan University, China
Hongyan Ren
Affiliation:
Mental Health Center and Psychiatric Laboratory, West China Hospital of Sichuan University, China; West China Brain Research Center, West China Hospital of Sichuan University, China
Yingcheng Wang
Affiliation:
Mental Health Center and Psychiatric Laboratory, West China Hospital of Sichuan University, China; and West China Brain Research Center, West China Hospital of Sichuan University, China
Xun Hu
Affiliation:
West China Brain Research Center, West China Hospital of Sichuan University, China; and Biobank, West China Hospital of Sichuan University, China
Guangya Zhang
Affiliation:
Suzhou Psychiatric Hospital, The Affiliated Guangji Hospital of Soochow University, China
Xiaodong Du
Affiliation:
Suzhou Psychiatric Hospital, The Affiliated Guangji Hospital of Soochow University, China
Qiang Wang*
Affiliation:
Mental Health Center and Psychiatric Laboratory, West China Hospital of Sichuan University, China; and West China Brain Research Center, West China Hospital of Sichuan University, China
Tao Li
Affiliation:
Mental Health Center and Psychiatric Laboratory, West China Hospital of Sichuan University, China; and West China Brain Research Center, West China Hospital of Sichuan University, China
*
Correspondence: Qiang Wang. Email: wangqiang130@scu.edu.cn
Rights & Permissions [Opens in a new window]

Abstract

Background

Understanding the patterns of treatment response is critical for the treatment of patients with schizophrenia; one way to achieve this is through using a longitudinal dynamic process study design.

Aims

This study aims to explore the response trajectory of antipsychotics and compare the treatment responses of seven different antipsychotics over 6 weeks in patients with schizoprenia (trial registration: Chinese Clinical Trials Registry Identifier: ChiCTR-TRC-10000934).

Method

Data were collected from a multicentre, randomised open-label clinical trial. Patients were evaluated with the Positive and Negative Syndrome Scale (PANSS) at baseline and follow-up at weeks 2, 4 and 6. Trajectory groups were classified by the method of k-means cluster modelling for longitudinal data. Trajectory analyses were also employed for the seven antipsychotic groups.

Results

The early treatment response trajectories were classified into a high-trajectory group of better responders and a low-trajectory group of worse responders. The results of trajectory analysis showed differences compared with the classification method characterised by a 50% reduction in PANSS scores at week 6. A total of 349 patients were inconsistently grouped by the two methods, with a significant difference in the composition ratio of treatment response groups using these two methods (χ2 = 43.37, P < 0.001). There was no differential contribution of high- and low trajectories to different drugs (χ2 = 12.52, P = 0.051); olanzapine and risperidone, which had a larger proportion in the >50% reduction at week 6, performed better than aripiprazole, quetiapine, ziprasidone and perphenazine.

Conclusions

The trajectory analysis of treatment response to schizophrenia revealed two distinct trajectories. Comparing the treatment responses to different antipsychotics through longitudinal analysis may offer a new perspective for evaluating antipsychotics.

Type
Papers
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Author(s) 2020. Published by Cambridge University Press on behalf of Royal College of Psychiatrists

Background

Antipsychotic drugs are currently the mainstay of treatment for schizophrenia. A meta-analysis of the time course of treatment response reported an early response profile for antipsychotic treatment.Reference Sherwood, Thornton and Honer1 However, previous studies have shown that the treatment response to schizophrenia is heterogeneous. Although most patients’ symptoms improve after treatment, approximately a third of patients respond poorly and are considered to be drug resistant.Reference Elkis and Buckley2,Reference Lehman, Lieberman, Dixon, McGlashan, Miller and Perkins3 Some studies have even found that the drugs available to treat schizophrenia are only effective in approximately 50% of patients.Reference Lieberman, Stroup, McEvoy, Swartz, Rosenheck and Perkins4 A poor response to treatment may lead to worse community functioning,Reference Iasevoli, Giordano, Balletta, Latte, Formato and Prinzivalli5 symptom exacerbation,Reference Liu-Seifert, Adams and Kinon6 relapse,Reference Ayuso-Gutierrez and del Rio Vega7 reduced patient adherence,Reference Perkins8 increased risk of admission to hospital, and heavy social and economic burden.Reference Kennedy, Altar, Taylor, Degtiar and Hornberger9 In particular, in the early stages of treatment with antipsychotic drugs, symptom improvement is not stable; therefore, longitudinal profiles of the early efficacy of different antipsychotics can provide novel insight into their different characteristics.

Importance of understanding response patterns to medication and use of machine learning

Understanding antipsychotic response patterns is critical for successful treatment. The identification of individual response trajectories and potential risk factors associated with poor response trajectories can be a guide for personalised treatments and effective interventions. Previous studies have shown that trajectory analysis may be a useful strategy to reduce heterogeneity in treatment for schizophrenia and may provide novel insight into clinically meaningful patient subgroups.Reference Honer, Jones, Thornton, Barr, Procyshyn and Vila-Rodriguez10 Some studies with small sample sizes have used parametric analyses and discerned four or more trajectories, but some trajectories merely account for 2% of the sample or even less.Reference Stauffer, Case, Kollack-Walker, Ascher-Svanuma, Ball and Kapur11Reference Levine, Rabinowitz, Case and Ascher-Svanum16

When using model-based procedures, the tasks are clustered by data without a priori hypotheses. Thus, this approach would provide a better way to avoid issues concerning model selection. K-means for longitudinal data (Kml) is an approach related to the partitional classification/cluster analysis of machine learning, which may be an appropriate measure to address model selection issues.Reference Genolini and Falissard17 The k-means algorithm uses variance distance or dissimilarity measures to identify and classify trajectories, which, as a non-parametric classification method, does not require any assumption regarding normality or parametric assumptions within clusters or the shape of the trajectory and can help to fit the data more closely.Reference Genolini and Falissard17,Reference Genolini and Falissard18 Various examples of the implementation of this method, including estimating developmental trajectories of hyperactivity/inattention symptoms, can be found in several studies.Reference Pingault, Tremblay, Vitaro, Japel, Boivin and Cote19Reference Abdennour, Reggio, Le Naour, Liu, Poitou and Aron-Wisnewsky21 However, to our knowledge, few studies have applied such an objective data-driven Kml approach to identify homogeneous trajectories of treatment response in patients with schizophrenia. Therefore, our current study is the first to evaluate antipsychotic drugs from a longitudinal perspective. Moreover, our study has a larger sample size than previous studies.Reference Stauffer, Case, Kollack-Walker, Ascher-Svanuma, Ball and Kapur11Reference Leucht, Davis, Engel, Kane and Wagenpfeil13,Reference Levine, Rabinowitz, Faries, Lawson and Ascher-Svanum15,Reference Levine, Rabinowitz, Case and Ascher-Svanum16

Aims

The aims of our study are (a) to explore the response trajectory of antipsychotic drugs with the unsupervised machine learning Kml method and to describe the subgroup features and baseline characteristics of patients with similar treatment responses and (b) to compare the longitudinal treatment responses for seven different antipsychotics.

Method

Participants

The participants were recruited from a multicentre, randomised open-label clinical trial that was designed to compare the treatment effectiveness of seven antipsychotic drugs in patients with first-onset psychosis and in those who relapsed. We recruited 3030 patients with schizophrenia from 32 different clinical settings from 6 July 2010 to 30 November 2011, all of whom were enrolled in the China Antipsychotic Pharmacogenomics Consortium. Details of the inclusion and exclusion criteria are provided in the Supplementary File available at https://doi.org/10.1192/bjo.2020.105 (all supplementary figures and tables mentioned can be found in the Supplementary File). Some studies have already been published using these data.Reference Wang, Man Wu, Yue, Yan, Zhang and Tan22,Reference Zhang, Wang, Reynolds, Yue, Deng and Yan23

Consensus diagnoses were made by at least two experienced psychiatrists based on unstructured interviews with the patients and patients’ families, as well as the review of patients’ medical records. According to the study protocol, participants were randomly assigned (1:1:1:1:1:½:½) to five atypical antipsychotics (risperidone, olanzapine, quetiapine, aripiprazole and ziprasidone) and two typical antipsychotics (perphenazine and haloperidol) using a random allocation sequence generated by a computer. The random allocation was unmasked only to patients and their psychiatrists but not to the researchers following baseline assessments.

The first 2 weeks of treatment were used for drug titration, and the next 4 weeks were the maintenance period. The upper limit of the dose of antipsychotic drugs had to be higher than the lower limit of any target dose range, as recommended by the International Consensus Study of Antipsychotic Dosing.Reference Gardner, Murphy, O'Donnell, Centorrino and Baldessarini24 Other psychotropic drugs were not allowed to be used as concomitant drugs, except for short-acting benzodiazepines for insomnia as well as lorazepam for agitation and psychotic anxiety. The anticholinergic drug hyoscine (up to 6 mg per day) could be prescribed for extrapyramidal symptoms and the beta-blocker propranolol (up to 80 mg per day) could be prescribed for akathisia.

The study was approved by the institutional review board at each site and was conducted by Good Clinical Practice guidelines and the Helsinki Declaration. Written informed consent to participate was obtained from both the patients and their legal guardians. The trial was registered with the Chinese Clinical Trials number: ChiCTR-TRC-10000934.

Clinical measurements

Demographic and clinical characteristics were collected by trained researchers at baseline. Illness severity was evaluated every 2 weeks using the Positive and Negative Syndrome Scale (PANSS).Reference Kay, Fiszbein and Opler25 The measurements of clinical manifestations included the total score and the positive symptom subscale, negative symptom subscale and general psychopathology subscale scores.

Side-effects were evaluated using the Barnes Akathisia Scale,Reference Barnes26 Abnormal Involuntary Movement ScaleReference Guy27 and Extrapyramidal Symptom Rating Scale.Reference Simpson and Angus28 All raters had received systematic training, and high interrater reliability was achieved. Investigators, staff and patients were masked to the treatment assignments. Where no empirical data for dose conversion were available, we assumed the daily defined dose (i.e. the average maintenance dose per day calculated from the dose recommendations in each drug's product information according to the World Health Organization) to be equivalent.Reference Patel, Arista, Taylor and Barnes29 So we collected the drug doses of all participants at the end of the second week and converted them into risperidone doses by daily defined dose for comparison.

The per cent change in the PANSS score from baseline was calculated by subtracting each value from baseline, dividing by the baseline score, subtracting 30 (the lowest possible PANSS total score), and finally multiplying by 100.Reference Leucht, Davis, Engel, Kissling and Kane30 As the subscale score of some participants at baseline was equal to the lowest score of the subscale, the reduction rate in the PANSS subscale was done by subtracting each value from baseline and dividing by the baseline score.

Data values that were 1.5 × interquartile range (IQR) higher than the third quartile or 1.5 × IQR lower than the first quartile were considered outliers. These values were replaced by values at the 0.01 and 0.99 quantiles. Multiple imputation using the R package MICE was applied to impute missing data.Reference Buuren31

Statistical analysis

The chi-squared test or Fisher's exact test was used to compare the differences in the distribution of categorical variables among the seven drug groups. Group differences regarding outcome criteria were additionally tested using the chi-squared test or Fisher's exact test for count data; z-test/t-test was applied for continuous variables; the Wilcoxon rank-sum test was applied for continuous data if the variables did not conform to a normal distribution. After running Fisher's exact test to compare the seven groups, post hoc pairwise comparisons adjusted by the false discovery rate (FDR) were performed to identify the differences between groups.

Treatment response trajectories

To identify similar trajectories of treatment response (reduction rate) during the 6 weeks, k-means for clustering the longitudinal data was used with the R Package Kml.Reference Genolini and Falissard17,Reference Genolini and Falissard18 Kml was applied with default settings, allowing k-means to run for two to six clusters 20 times each. This non-parametric method classified participants into different trajectories, i.e. homogeneous subgroups with similar treatment response patterns. Briefly, each patient was first assigned arbitrarily to one initial trajectory. Second, the centre (k-mean) of each trajectory was calculated, and each participant was reassigned to the closest trajectory. The operation was repeated until convergence was achieved. The process from assignment to convergence was then repeated (20 times in this study) to ensure that the solution was not dependent on the initial assignment. Finally, the best solution was determined by a criterion. The Calinski–Harabasz index was applied in this study. The Calinski–Harabasz index describes the closeness through the intraclass dispersion matrix and the separation of the interclass dispersion matrix. The larger the Calinski–Harabasz index is, the better the clustering result is.

All the statistical analyses described above were performed in R 3.6.1.

Results

Participant characteristics

Of the 3030 patients with schizophrenia, 20 were excluded from further evaluation because of the failure to meet the inclusion criteria or to participate in the study. In total, 2630 participants finished the study. The study sample and the reasons for dropping out are shown in Fig. 1. Participants who dropped out (n = 380, 12.62%) tended to be younger and drug-naive and to have higher educational degrees, a shorter duration of illness and milder symptoms. However, there was no difference in the type of medication between the patients who completed the study and those who dropped out (the details are shown in Supplementary Table 1). The baseline demographic characteristics of the patients treated with the seven antipsychotic drugs are shown in Supplementary Table 2. The pattern of missing data is shown in Supplementary Fig. 1, and multiple imputation was applied to impute missing data. There was no significant difference in the data after imputation (P > 0.05) (Supplementary Fig. 2).

Fig. 1 The trial profile.

Trajectories of treatment response

In the trajectory analyses of treatment response, values of fit criteria improved steadily from a six-trajectory solution to a two-trajectory solution. The Calinski–Harabasz index of two-trajectory solution was the largest (see Supplementary Table 3). According to the fit criterion, we yielded a two-trajectory solution as the best model (Fig. 2).

Fig. 2 Kml trajectory analysis of treatment response (all patients).

(a) The optimal number of clusters to separate patients into groups with homogeneous treatment response over time. The x-axis represents the number of runs for two to six clusters. The y-axis represents the Calinski–Harabasz index. The curve marked with the number 2 represents the two-trajectory solution and explains the data best as it has the highest Calinski–Harabasz index. (b) The treatment trajectory of all patients. The y-axis represents the reduction rate of PANSS scale. Light green line (B) corresponds to a high trajectory or better treatment response (48.9% of patients) and dark green line (A) corresponds to the low trajectory or worse treatment response (51.1% of patients). The thin lines (black) represent individual patient profiles.

In this pattern, 1471 (48.9% of the patients who entered the study) participants were assigned to the high-trajectory group, and they responded better to treatment. The average rates of change in PANSS total score in the high-trajectory group were 35.84%, 60.51%, and 73.58% at weeks 2, 4 and 6, respectively. The low-trajectory group included 1539 (51.1% of the patients who entered the study) patients. The average rates of change in the PANSS total score were 15.16%, 28.10%, and 31.99% at weeks 2, 4 and 6, respectively. The treatment response of the low-trajectory group was also poorer in terms of positive symptoms, negative symptoms and general symptoms (see details in Supplementary Table 4 and Fig. 3). There was a significant difference in the reduction rate of the PANSS total score and subscale scores between patients in the high- and low-trajectory groups (Table 1, Supplementary Table 4 and Figs 3 and 4). In addition, we performed the same analysis on patients in a first episode, and the results were similar (Supplementary Fig. 5).

Table 1 Comparisons of demographics and baseline characteristics of patients separated into high- and low trajectories (all patients)

PANSS, Positive and Negative Syndrome Scale.

a. The drug doses was converted into risperidone doses by daily defined dose.

The baseline characteristics and outcomes of the high- and low-trajectory groups are shown in Table 1. Participants in the high-trajectory group had a later age of onset, a shorter duration of disease and a lower dose of drugs. There were also differences in the distribution of educational degrees between the two groups. There was no significant difference in other baseline characteristics between the two groups. Among the participants in a first-episode, there was no difference in education between the trajectory groups, and the remaining results were very similar (Supplementary Table 5).

Comparisons between trajectory analysis and the dichotomous thresholds method

Previous studies have suggested that at least a 50% reduction in the PANSS score from baseline is the cut-off value for defining treatment response.Reference Leucht, Kane, Kissling, Hamann, Etschel and Engel32 Based on that, we compared the difference between the trajectory analysis and simple dichotomous thresholds (as indentified by the responder) to evaluate the treatment response. A total of 349 patients were inconsistently grouped by the two methods. Forty-seven patients were in the high-trajectory group with the reduction rate in the PANSS score in the sixth week at less than 50%; 302 patients were in the low-trajectory group, but the reduction rate in the sixth week was more than 50% (Supplementary Fig. 6). The average rates of change in the PANSS total score of these four groups is shown in Supplementary Fig. 7.

In addition, there was a significant difference in the composition ratio of treatment response groups using these two methods (χ2 = 43.37, P < 0.001). The percentage of patients with a >50% reduction in the PANSS score in the sixth week (n = 1726, 57.34%) was significantly higher than that of patients in the high-trajectory group (n = 1471, 48.87%) (see Fig. 3(a)).

Fig. 3 Comparison of findings using the trajectory analysis and dichomtomous threshold methods. (a) The proportion of patients using trajectory analysis and dichotomous thresholds methods for a good response (dark green) versus a poor response (light green). (b) The proportion of patients in a first-episode and patients who relapsed by trajectories (high or low). Data for patients treated with atypical antipsychotics and typical antipsychotics are shown separately. (c) The proportion of patients taking each of the seven antipsychotic drugs by trajectories (high or low). (d) The proportion of patients taking each of seven antipsychotic drugs by dichotomous thresholds methods. Reduction rate at week 6 >50% or <50%.

Percentage figures are indicated within bars.

We also compared the baseline characteristics of the treatment response groups by dichotomous threshold methods. The results showed that there were significant differences in the age at first onset, the duration of illness and the level of education between the two groups (P < 0.05), which was consistent with the results of the Kml grouping (Supplementary Table 6).

Comparison of antipsychotic drugs

Longitudinal trajectory analyses were run separately for each antipsychotic drug group. Furthermore, the analyses were performed on those participants in a first episode and participants who had relapsed in each antipsychotic drug group. These results showed that all drug groups could accept a two-trajectory solution as the best model, and the profiles of response to seven antipsychotics were similar (see Supplementary Figs. 8–15).

In the groups of participants with typical antipsychotic drugs, the proportion was different between the first-episode group and the relapse group (χ2 = 10.51, P = 0.001). More participants were classified into the high-trajectory group in the first-onset group (n = 85, 59.9%), whereas more patients were classified into the low-trajectory group among the relapse group (n = 154, 43.80%). This pattern was not found in the atypical antipsychotic drug groups (χ2 = 1.33, P = 0.25) (Supplementary Table 7 and Fig. 3(b)).

Comparing the treatment response trajectories of the seven drugs separately, there was no differential contribution of the composition of the high-trajectory and low-trajectory groups to the different drugs (χ2 = 12.52, P = 0.051). The results of the analyses are shown in Fig. 3(c) and Supplementary Table 7.

We also compared the difference in the composition ratio of seven drugs grouped by a cut-off value of 50% for the reduction rate at week 6, and the results were significantly different. Olanzapine and risperidone, which had a higher proportion of >50% reduction in the sixth week, performed better than aripiprazole, quetiapine, ziprasidone and perphenazine. The difference was still found after FDR correction (P FDR < 0.05) (Supplementary Table 8).

Discussion

Treatment response trajectories

In this study, we performed unsupervised machine learning Kml to classify homogeneous clusters of 6-week response trajectories in patients with schizophrenia. The patients in the high-trajectory group (48.9%) responded well to the treatment, as evidenced by the approximately 70% reduction in the PANSS total score in the sixth week. The level of improvement was considered to be the degree of ‘very much improved’ at a global clinical level.Reference Leucht, Kane, Etschel, Kissling, Hamann and Engel33 Furthermore, the other group of patients (51.1%) showed an approximately 30% improvement in their symptoms across the course of the 6-week trial, which was only marginally discernible.Reference Leucht, Kane, Etschel, Kissling, Hamann and Engel33

The treatment response trajectory in our study is different from the results of previous trajectory analyses.Reference Stauffer, Case, Kollack-Walker, Ascher-Svanuma, Ball and Kapur11,Reference Levine, Rabinowitz, Faries, Lawson and Ascher-Svanum15,Reference Levine, Rabinowitz, Case and Ascher-Svanum16 Most previous studies used the total score of the scale rather than the reduction rate to analyse the trajectory of treatment response, which may be one of the reasons for the difference in results. However, clinical trials often used the change from baseline scores to evaluate the effects of different interventions, and ‘response’ is defined as a scale score reduction rate that exceeds a specific value.Reference Leucht, Davis, Engel, Kane and Wagenpfeil13 Thus, the use of the reduction rate of the total PANSS score in our study may be more appropriate. Moreover, our study used the Kml method, which does not require any normality or parametric assumptions within clusters.Reference Genolini and Falissard18 This is important in the present case because of the nature of the variable that was clustered, i.e. the reduction rate of the total PANSS score. A previous trajectory analysis used the Kml technique and they found that over 70% of the patients with schizophrenia treated with clozapine and chlorpromazine were assigned to the good treatment response trajectory, whereas this figure was approximately 50% in the refractory group of patients.Reference Honer, Jones, Thornton, Barr, Procyshyn and Vila-Rodriguez10 However, the small size of the study may have resulted in an overestimation of the proportion of patients expected to have a good treatment response trajectory.

Whether baseline characteristics affect the treatment response has shown contradictory results in previous studies.Reference Carbon and Correll34 Some have been able to come to such a conclusion, finding that gender,Reference Stauffer, Case, Kollack-Walker, Ascher-Svanuma, Ball and Kapur11,Reference Levine and Rabinowitz35 age,Reference Stauffer, Case, Kollack-Walker, Ascher-Svanuma, Ball and Kapur11,Reference Stentebjerg-Olesen, Jeppesen, Pagsberg, Fink-Jensen, Kapoor and Chekuri36 age at onsetReference Levine and Rabinowitz35,Reference Schennach-Wolff, Jäger, Mayr, Meyer, Kühn and Klingberg37,Reference Smart, Kępińska, Murray and MacCabe38 and illness durationReference Rabinowitz, Werbeloff, Caers, Mandel and Kapur39 affect the treatment response. Our study found that participants in the high-trajectory group had a later age at onset and a shorter duration of disease, which is similar to previous studies. The dose of the drug was lower in the high-trajectory group, which may be related to the better response to treatment in this group. The results obtained by trajectory analysis were similar to those obtained by the dichotomous thresholds method, which also increases the credibility of our results. However, the identified trajectory groups may have differed at baseline in other ways that were not assessed in this analysis. For example, genetic testing and imaging examination might help to further differentiate response trajectories within this large, pooled patient population.

Comparing treatment trajectories among different antipsychotics

In regard to the trajectories of different antipsychotic drugs, all drug groups could be classified into two trajectories. In the results of the longitudinal effect of typical and atypical antipsychotics, typical antipsychotics performed better in the first-onset group than the relapse group, whereas atypical antipsychotics performed equally in the two groups. Previous studies also found that antipsychotic treatment responses were reduced or delayed in the face of relapse.Reference Takeuchi, Siu, Remington, Fervaha, Zipursky and Foussias40 However, the smaller sample size of typical antipsychotic drugs may also have resulted in the difference.

Our study also proposes a new perspective to evaluate antipsychotic drugs. To date, there have been limited longitudinal studies on antipsychotic drugs. An analysis of two 6-week trials showed that antipsychotic treatment with haloperidol or olanzapine was more likely than placebo to be associated with a trajectory of ‘dramatic response’.Reference Marques, Arenovich, Agid, Sajeev, Muthén and Chen12 Previously, some trajectory analyses suggested possible differences in treatment responses among different antipsychotic drugs. One study showed that patients taking olanzapine had a higher proportion in the trajectory of patients with the most improvement than patients in the other drug groups.Reference Levine, Rabinowitz, Faries, Lawson and Ascher-Svanum15 Another analysis observed that patients in trajectories with clear worsening of symptoms were more likely to have been treated with quetiapine or risperidone and less likely to have been treated with olanzapine or aripiprazole.Reference Stauffer, Case, Kollack-Walker, Ascher-Svanuma, Ball and Kapur11 However, fewer kinds of drugs were included in these studies, and there was no analysis or comparison of the specific trajectories of each drug. The effect of drugs on a high trajectory might be more stable. Our results suggested that there was no difference in the distribution of seven antipsychotic drugs between the high- and low-trajectory groups by comparing the longitudinal drug effect.

Interestingly, by comparing the treatment response of the seven antipsychotics by dichotomous thresholds methods (comparing the composition ratio of patients whose reduction rate in the PANSS score was more than and less than 50% at the sixth week), the results showed that individuals taking olanzapine and risperidone responded better to treatment, which is significantly different from the results derived from the comparison using trajectory analysis. End-point threshold-based cut-off values, though convenient, may be misleading in terms of the underlying biology. For example, patients with 49% remission and those showing a 51% response would be divided into two groups, even if they belonged to the same treatment response trajectory. The end-point cut-off value does not take the instability of symptoms and the difference in the onset time of different patients into account. It is difficult to distinguish between transient and permanent results. Compared with dichotomous threshold methods, longitudinal analysis allows changes in treatment responses over time, takes into account all available data from each time point, and identifies treatment response groups consisting of patients whose response patterns are sufficiently similar to one another and sufficiently different from patients in other groups in an outcome measure over time.Reference Stauffer, Case, Kollack-Walker, Ascher-Svanuma, Ball and Kapur11,Reference Marques, Arenovich, Agid, Sajeev, Muthén and Chen12 The two methods for evaluating treatment response have a different emphasis. However, our study follow-up consisted of only 6 weeks. These approaches need to be applied in further research with a longer follow-up to validate the current conclusions.

Strengths and limitations

The current study also has several limitations. First, the trajectory analysis approach is exploratory rather than definitive and requires replication. Second, we focused on the treatment response trajectory only during the first 6 weeks of treatment. Given that schizophrenia is a chronic disease, a longer-term study is suggested as a key point of future studies to validate the current conclusions. Third, patients with a history of treatment resistance were excluded from the trial, which may have obscured the extent of treatment response and potential trajectories with extreme scores. Finally, the dose of each drug was constant. In the real world, the dose of drugs of individual patients varies with the change of symptoms. The constant dose in our study may not reflect the best effect of drugs, so it is necessary to study the treatment response further in a real-world sample. However, our study is the first to compare seven different antipsychotic drugs through their longitudinal effects, and it includes a large sample size that can better represent the population with schizophrenia. In addition, the use of the Kml method does not rely on prespecified hypotheses regarding the final number of clusters or the most reliable thresholds and contains information from several time points. Moreover, trajectory analysis uses PANSS per cent reduction rather than PANSS total scores to better evaluates treatment response.

Implications

In conclusion, the current study identified that patients with schizophrenia could be divided into two trajectories: a high-trajectory group with a better treatment response and a low-trajectory group with a poorer treatment response. Longitudinal comparisons of treatment response trajectories and dichotomous threshold methods of drug treatment response were significantly different. Therefore, comparing the treatment responses of different antipsychotics through longitudinal analysis may offer a new perspective to evaluate antipsychotic drugs.

Supplementary material

Supplementary material is available online at https://doi.org/10.1192/bjo.2020.105.

Data availability

The data are not publicly available due to their containing information that could compromise the privacy of research participants.

Acknowledgements

We express our appreciation to the China Antipsychotic Pharmacogenomics Consortium and all the participants who participated in this study for their cooperation and assistance.

Author contributions

Xin.M., F.Y., Li.L, C.W., W.Y., D.Z. and Tao.L. designed, supervised the study and revised the article. M.D., Yu.W. and Q.W. did the statistical analyses and prepared the tables and figures. M.D., Yu.W., Q.W., Y.T., and Tao.L. wrote and revised the manuscript. W.Y., Y.Z., W.D., Xia.M., L.Z., H.R., Xin.M., Hao.Yu, X.H., G.Z., X.D. and Yi.W. provided further data interpretation and revised the manuscript. Hao.Yan, L.T., G.Y., Ti.L., L.W., F.Z., J.Y., K.L, Lux.L., Q.T., H.Z., and Q.C. were responsible for patient recruitment and revised the manuscript. All authors read and approved the final manuscript.

Funding

This study was funded in part by grants (81771446, 81630030 and 81920108018) from the National Nature Science Foundation of China; grants 2016YFC0904300 and 2016YFC1307003 from the National Key Research and Development Program of the Ministry of Science and Technology of China; grants SZYJTD201715 from Introduction Project of Suzhou Clinical Expert Team; grants SZZX201509 from Suzhou Key Medical Center for Psychiatric Diseases; grants SS201752 from The Scientific and Technological Program of Suzhou; the Department of Science and Technology of Sichuan provincial government (2019YFS0153) and grants ZY2016103, ZY2016203 and ZYGD20004 from the 1.3.5 Project for Disciplines of Excellence of the West China Hospital, Sichuan University.

Declaration of interest

None.

ICMJE forms are in the supplementary material, available online at http://doi.org/10.1192/bjo.2020.105.

Footnotes

The online version of this article has been updated since original publication. A notice detailing the changes has also been published at https://doi.org/10.1192/bjo.2021.995.

References

Sherwood, M, Thornton, AE, Honer, WG. A meta-analysis of profile and time-course of symptom change in acute schizophrenia treated with atypical antipsychotics. Int J Neuropsychopharmacol 2006; 9: 357–66.CrossRefGoogle ScholarPubMed
Elkis, H, Buckley, PF. Treatment-resistant schizophrenia. Psychiatr Clin North Am 2016; 39: 239–65.Google ScholarPubMed
Lehman, AF, Lieberman, JA, Dixon, LB, McGlashan, TH, Miller, AL, Perkins, DO, et al. Practice guideline for the treatment of patients with schizophrenia, second edition. Am J Psychiatry 2004; 161 (2 suppl): 156.Google ScholarPubMed
Lieberman, JA, Stroup, TS, McEvoy, JP, Swartz, MS, Rosenheck, RA, Perkins, DO, et al. Effectiveness of antipsychotic drugs in patients with chronic schizophrenia. New Engl J Med 2005; 353: 1209–23.CrossRefGoogle ScholarPubMed
Iasevoli, F, Giordano, S, Balletta, R, Latte, G, Formato, MV, Prinzivalli, E, et al. Treatment resistant schizophrenia is associated with the worst community functioning among severely-ill highly-disabling psychiatric conditions and is the most relevant predictor of poorer achievements in functional milestones. Prog Neuropsychopharmacol Biol Psychiatry 2016; 65: 3448.CrossRefGoogle ScholarPubMed
Liu-Seifert, H, Adams, DH, Kinon, BJ. Discontinuation of treatment of schizophrenic patients is driven by poor symptom response: a pooled post-hoc analysis of four atypical antipsychotic drugs. BMC Med 2005; 3: 21.CrossRefGoogle ScholarPubMed
Ayuso-Gutierrez, JL, del Rio Vega, JM. Factors influencing relapse in the long-term course of schizophrenia. Schizophr Res 1997; 28: 199206.CrossRefGoogle ScholarPubMed
Perkins, DO. Predictors of noncompliance in patients with schizophrenia. J Clin Psychiatry 2002; 63: 1121–8.CrossRefGoogle ScholarPubMed
Kennedy, JL, Altar, CA, Taylor, DL, Degtiar, I, Hornberger, JC. The social and economic burden of treatment-resistant schizophrenia: a systematic literature review. Int Clin Psychopharmacol 2014; 29: 6376.CrossRefGoogle ScholarPubMed
Honer, WG, Jones, AA, Thornton, AE, Barr, AM, Procyshyn, RM, Vila-Rodriguez, F. Response trajectories to clozapine in a secondary analysis of pivotal trials support using treatment response to subtype schizophrenia. Can J Psychiatry 2015; 60: S1925.Google Scholar
Stauffer, V, Case, M, Kollack-Walker, S, Ascher-Svanuma, H, Ball, T, Kapur, S, et al. Trajectories of response to treatment with atypical antipsychotic medication in patients with schizophrenia pooled from 6 double-blind, randomized clinical trials. Schizophr Res 2011; 130: 11–9.CrossRefGoogle ScholarPubMed
Marques, TR, Arenovich, T, Agid, O, Sajeev, G, Muthén, B, Chen, L, et al. The different trajectories of antipsychotic response: antipsychotics versus placebo. Psychol Med 2011; 41: 1481–8.CrossRefGoogle ScholarPubMed
Leucht, S, Davis, JM, Engel, RR, Kane, JM, Wagenpfeil, S. Defining ‘response’ in antipsychotic drug trials: recommendations for the use of scale-derived cutoffs. Neuropsychopharmacology 2007; 32: 1903–10.CrossRefGoogle ScholarPubMed
Saeed, F, Ofer, A, George, F, Remington, G. Using treatment response to subtype schizophrenia: proposal for a new paradigm in classification. Schizophr Bull 2013; 39: 1169–72.Google Scholar
Levine, SZ, Rabinowitz, J, Faries, D, Lawson, AH, Ascher-Svanum, H. Treatment response trajectories and antipsychotic medications: examination of up to 18 months of treatment in the CATIE chronic schizophrenia trial. Schizophr Res 2012; 137: 141–6.CrossRefGoogle ScholarPubMed
Levine, SZ, Rabinowitz, J, Case, M, Ascher-Svanum, H. Treatment response trajectories and their antecedents in recent-onset psychosis: a 2-year prospective study. J Clin Psychopharmacol 2010; 30: 446–9.CrossRefGoogle ScholarPubMed
Genolini, C, Falissard, B. KmL: K-means for longitudinal data. Comput Stat 2010; 25: 317–28.CrossRefGoogle Scholar
Genolini, C, Falissard, B. KmL: a package to cluster longitudinal data. Comput Methods Programs Biomed 2011; 104: e112–21.CrossRefGoogle ScholarPubMed
Pingault, JB, Tremblay, RE, Vitaro, F, Japel, C, Boivin, M, Cote, SM. Early nonparental care and social behavior in elementary school: support for a social group adaptation hypothesis. Child Dev 2015; 86: 1469–88.CrossRefGoogle ScholarPubMed
Pingault, JB, Cote, SM, Galera, C, Genolini, C, Falissard, B, Vitaro, F, et al. Childhood trajectories of inattention, hyperactivity and oppositional behaviors and prediction of substance abuse/dependence: a 15-year longitudinal population-based study. Mol Psychiatry 2013; 18: 806–12.CrossRefGoogle ScholarPubMed
Abdennour, M, Reggio, S, Le Naour, G, Liu, Y, Poitou, C, Aron-Wisnewsky, J, et al. Association of adipose tissue and liver fibrosis with tissue stiffness in morbid obesity: links with diabetes and BMI loss after gastric bypass. J Clin Endocrinol Metab 2014; 99: 898907.CrossRefGoogle ScholarPubMed
Wang, Q, Man Wu, H, Yue, W, Yan, H, Zhang, Y, Tan, L. Effect of damaging rare mutations in synapse-related gene sets on response to short-term antipsychotic medication in chinese patients with schizophrenia: a randomized clinical trial. JAMA Psychiatry 2018; 75: 1261–9.Google ScholarPubMed
Zhang, Y, Wang, Q, Reynolds, GP, Yue, W, Deng, W, Yan, H. Metabolic effects of 7 antipsychotics on patients with schizophrenia: a short-term, randomized, open-label, multicenter, pharmacologic trial. J Clin Psychiatry 2020; 81(3): 19m12785.CrossRefGoogle ScholarPubMed
Gardner, DM, Murphy, AL, O'Donnell, H, Centorrino, F, Baldessarini, RJ. International consensus study of antipsychotic dosing. Am J Psychiatry 2010; 167: 686–93.CrossRefGoogle ScholarPubMed
Kay, SR, Fiszbein, A, Opler, LA. The positive and negative syndrome scale (PANSS) for schizophrenia. Schizophr Bull 1987; 13: 261–76.CrossRefGoogle Scholar
Barnes, TR. A rating scale for drug-induced akathisia. Br J Psychiatry 1989; 154: 672–6.CrossRefGoogle ScholarPubMed
Guy, W. ECDEU Assessment Manual for Psychopharmacology. US Department of Health, Education, and Welfare, Public Health Service, 1976.Google Scholar
Simpson, G, Angus, J. A rating scale for extrapyramidal side effects. Acta Psychiatr Scand 1970; 45: 11–9.CrossRefGoogle Scholar
Patel, M, Arista, I, Taylor, M, Barnes, T. How to compare doses of different antipsychotics: a systematic review of methods. Schizophr Res 2013; 149: 141–8.CrossRefGoogle ScholarPubMed
Leucht, S, Davis, JM, Engel, RR, Kissling, W, Kane, JM. Definitions of response and remission in schizophrenia: recommendations for their use and their presentation&nbsp. Acta Psychiatr Scand 2008; 119 (s438): 714.Google Scholar
Buuren, S. Groothuis-Oudshoorn K. mice: multivariate imputation by chained equations in R. J Stat Softw 2011; 45: 168.Google Scholar
Leucht, S, Kane, JM, Kissling, W, Hamann, J, Etschel, E, Engel, RR. What does the PANSS mean? Schizophr Res 2005; 79: 231–8.CrossRefGoogle ScholarPubMed
Leucht, S, Kane, JM, Etschel, E, Kissling, W, Hamann, J, Engel, RR. Linking the PANSS, BPRS, and CGI: clinical implications. Neuropsychopharmacology 2006; 31: 2318–25.CrossRefGoogle ScholarPubMed
Carbon, M, Correll, CU. Clinical predictors of therapeutic response to antipsychotics in schizophrenia. Dialogues Clin Neurosci 2014; 16: 505–24.Google Scholar
Levine, SZ, Rabinowitz, J. Trajectories and antecedents of treatment response over time in early-episode psychosis. Schizophr Bull 2010; 36: 624–32.CrossRefGoogle ScholarPubMed
Stentebjerg-Olesen, M, Jeppesen, P, Pagsberg, AK, Fink-Jensen, A, Kapoor, S, Chekuri, R. Early nonresponse determined by the clinical global impressions scale predicts poorer outcomes in youth with schizophrenia spectrum disorders naturalistically treated with second-generation antipsychotics. J Child Adolesc Psychopharmacol 2013; 23: 665–75.CrossRefGoogle ScholarPubMed
Schennach-Wolff, R, Jäger, M, Mayr, A, Meyer, S, Kühn, K-U, Klingberg, S. Predictors of response and remission in the acute treatment of first-episode schizophrenia patients–is it all about early response? Eur Neuropsychopharmacol 2011; 21: 370–8.CrossRefGoogle ScholarPubMed
Smart, SE, Kępińska, AP, Murray, RM, MacCabe, JH. Predictors of treatment resistant schizophrenia: a systematic review of prospective observational studies. Psychol Med [Epub ahead of print] 29 Aug 2019. Available from: https://doi.org/10.1017/S0033291719002083.Google ScholarPubMed
Rabinowitz, J, Werbeloff, N, Caers, I, Mandel, FS, Kapur, S. Determinants of antipsychotic response in schizophrenia: implications for practice and future clinical trials. J Clin Psychiatry 2014; 75: e308–16.CrossRefGoogle ScholarPubMed
Takeuchi, H, Siu, C, Remington, G, Fervaha, G, Zipursky, RB, Foussias, G. Does relapse contribute to treatment resistance? Antipsychotic response in first- vs. second-episode schizophrenia. Neuropsychopharmacology 2019; 44: 1036–42.CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1 The trial profile.

Figure 1

Fig. 2 Kml trajectory analysis of treatment response (all patients).(a) The optimal number of clusters to separate patients into groups with homogeneous treatment response over time. The x-axis represents the number of runs for two to six clusters. The y-axis represents the Calinski–Harabasz index. The curve marked with the number 2 represents the two-trajectory solution and explains the data best as it has the highest Calinski–Harabasz index. (b) The treatment trajectory of all patients. The y-axis represents the reduction rate of PANSS scale. Light green line (B) corresponds to a high trajectory or better treatment response (48.9% of patients) and dark green line (A) corresponds to the low trajectory or worse treatment response (51.1% of patients). The thin lines (black) represent individual patient profiles.

Figure 2

Table 1 Comparisons of demographics and baseline characteristics of patients separated into high- and low trajectories (all patients)

Figure 3

Fig. 3 Comparison of findings using the trajectory analysis and dichomtomous threshold methods. (a) The proportion of patients using trajectory analysis and dichotomous thresholds methods for a good response (dark green) versus a poor response (light green). (b) The proportion of patients in a first-episode and patients who relapsed by trajectories (high or low). Data for patients treated with atypical antipsychotics and typical antipsychotics are shown separately. (c) The proportion of patients taking each of the seven antipsychotic drugs by trajectories (high or low). (d) The proportion of patients taking each of seven antipsychotic drugs by dichotomous thresholds methods. Reduction rate at week 6 >50% or <50%.Percentage figures are indicated within bars.

Supplementary material: PDF

Dai et al. supplementary material

Dai et al. supplementary material 1

Download Dai et al. supplementary material(PDF)
PDF 1.3 MB
Supplementary material: File

Dai et al. supplementary material

Dai et al. supplementary material 2

Download Dai et al. supplementary material(File)
File 2.6 MB
Submit a response

eLetters

No eLetters have been published for this article.