Introduction
Individuals with first episode psychosis (FEP) have a reduced life expectancy of 15–20 years (Kallio et al., Reference Kallio, Korkeila, Malmberg, Gunn, Rautava, Korhonen and Kytö2022; Plana-Ripoll et al., Reference Plana-Ripoll, Pedersen, Agerbo, Holtz, Erlangsen, Canudas-Romo, Andersen, Charlson, Christensen, Erskine, Ferrari, Iburg, Momen, Mortensen, Nordentoft, Santomauro, Scott, Whiteford, Weye and Laursen2019; Wahlbeck et al., Reference Wahlbeck, Westman, Nordentoft, Gissler and Laursen2011) with high mortality risks even decades after illness onset. A meta-analysis reported that nearly 13% of individuals with psychosis deceased within 10 years of diagnosis, and this mortality rate continues to rise over time (Chesney, Goodwin, & Fazel, Reference Chesney, Goodwin and Fazel2014). Two nationwide studies from Denmark and Finland (Lomholt et al., Reference Lomholt, Andersen, Sejrsgaard-Jacobsen, Øzdemir, Graff, Schjerning, Jensen, Straszek, Licht, Grøntved and Nielsen2019; Tanskanen, Tiihonen, & Taipale, Reference Tanskanen, Tiihonen and Taipale2018) suggested that the substantial mortality gap between individuals with schizophrenia and the general population remained over decades. Understanding the multifactorial contributors of excess mortality in individuals with psychosis, both unnatural and natural causes would be crucial for interventions aiming to narrow this gap (Laursen, Nordentoft, & Mortensen, Reference Laursen, Nordentoft and Mortensen2014; Nordentoft et al., Reference Nordentoft, Wahlbeck, Hällgren, Westman, Osby, Alinaghizadeh, Gissler and Laursen2013).
Many studies have investigated excess mortality in individuals with FEP, but few have had longer-term follow-ups. One study (Doyle et al., Reference Doyle, O’Keeffe, Hannigan, Kinsella, Behan, Kelly, Sheridan, Madigan, Lawlor and Clarke2019) in Ireland, tracking FEP patients for over 20 years revealed a persistently high mortality gap, particularly among young males, with mortality risk being highest during the first 8-year follow-up. Another study (Yuen et al., Reference Yuen, Harrigan, Mackinnon, Harris, Yuen, Henry, Jackson, Herrman and McGorry2014) focusing on the FEP population aged 14–30 years over a similar period, found that a significant portion of excess mortality was initially attributed to non-suicide unnatural causes, followed by an increase in natural deaths over time. These temporal variation in causes of death and mortality rates over time highlights the need for longer-term studies to provide a comprehensive and life-course perspective. Some longer-term follow-up studies (Cheng et al., Reference Cheng, Chang, Tsai, Li, Tsai, Bai, Lin, Su, Chen and Chen2023; Crump, Winkleby, Sundquist, & Sundquist, Reference Crump, Winkleby, Sundquist and Sundquist2013; Dickerson et al., Reference Dickerson, Khan, Origoni, Rowe, Katsafanas, Harvin, Yang and Yolken2024; Starzer et al., Reference Starzer, Hansen, Hjorthøj, Speyer, Albert and Nordentoft2023) focusing on individuals with schizophrenia-spectrum disorder (SSD), including a 20-year follow-up study of individuals with SSD in Denmark, which demonstrates a rising proportion of deaths attributable to physical illnesses, particularly cardiovascular and metabolic diseases (Starzer et al., Reference Starzer, Hansen, Hjorthøj, Speyer, Albert and Nordentoft2023). Studies of individuals with narrower diagnostic categories provide a more homogenous and focused understanding. However, given the longitudinal diagnostic instability of psychosis (Bromet et al., Reference Bromet, Kotov, Fochtmann, Carlson, Tanenberg-Karant, Ruggero and Chang2011; Peralta et al., Reference Peralta, Janda, de Jalón, Moreno-Izco, Sánchez-Torres, Cuesta and Peralta2024) and the provision of specialized services to individuals with FEP in many parts of the world, exploring the longitudinal mortality pattern of FEPs in general may provide a more comprehensive understanding.
Excess mortality patterns in individuals with FEP are likely to be shaped by sociodemographic, clinical, and behavioral factors. Sociodemographic determinants, including unemployment (Koutsouleris et al., Reference Koutsouleris, Kahn, Chekroud, Leucht, Falkai, Wobrock, Derks, Fleischhacker and Hasan2016), living alone (Starzer et al., Reference Starzer, Hansen, Hjorthøj, Speyer, Albert and Nordentoft2023), male gender (Lieslehto et al., Reference Lieslehto, Tiihonen, Lähteenvuo, Leucht, Correll, Mittendorfer-Rutz, Tanskanen and Taipale2024), and ethnic minority status (Morgan et al., Reference Morgan, Fearon, Lappin, Heslin, Donoghue, Lomas, Reininghaus, Onyejiaka, Croudace, Jones, Murray, Doody and Dazzan2017) have been suggested to be associated with increased mortality risk. Clinical factors such as hospitalization at illness onset (Lieslehto et al., Reference Lieslehto, Tiihonen, Lähteenvuo, Leucht, Correll, Mittendorfer-Rutz, Tanskanen and Taipale2024), prolonged duration of untreated psychosis (DUP), poor premorbid functioning (Stephens, Richard, & McHugh, Reference Stephens, Richard and McHugh1999), and low cognitive functioning (Doyle et al., Reference Doyle, O’Keeffe, Hannigan, Kinsella, Behan, Kelly, Sheridan, Madigan, Lawlor and Clarke2019) have all been demonstrated to be risk factors of mortality. Mental health comorbidities, particularly suicidality (Simon et al., Reference Simon, Stewart, Yarborough, Lynch, Coleman, Beck, Operskalski, Penfold and Hunkeler2018) and substance use (Correll et al., Reference Correll, Solmi, Croatto, Schneider, Rohani-Montez, Fairley, Smith, Bitter, Gorwood, Taipale and Tiihonen2022; Lieslehto et al., Reference Lieslehto, Tiihonen, Lähteenvuo, Leucht, Correll, Mittendorfer-Rutz, Tanskanen and Taipale2024), further amplify the risk of mortality. Another important clinical factor is the use of antipsychotics, with its relationship with mortality remains a subject of ongoing debate. Continuous use of antipsychotics has been shown to be associated with decreased all-cause mortality (Taipale et al., Reference Taipale, Tanskanen, Mehtälä, Vattulainen, Correll and Tiihonen2020; Tiihonen et al., Reference Tiihonen, Suokas, Suvisaari, Haukka and Korhonen2012) with a recent review highlighting the protective effects of any antipsychotics against all-cause and natural-cause mortality (Correll et al., Reference Correll, Solmi, Croatto, Schneider, Rohani-Montez, Fairley, Smith, Bitter, Gorwood, Taipale and Tiihonen2022). However, a systematic review has suggested that long-term antipsychotic medication may increase mortality risk, primarily due to adverse effects on physical health (Weinmann, Read, & Aderhold, Reference Weinmann, Read and Aderhold2009). Antipsychotic medications have been implicated in weight gain (Dayabandara et al., Reference Dayabandara, Hanwella, Ratnatunga, Seneviratne, Suraweera and de Silva2017), metabolic syndrome (De Hert et al., Reference De Hert, Detraux, van Winkel, Yu and Correll2011), and an increased risk of cardiovascular disease (De Hert et al., Reference De Hert, Detraux, van Winkel, Yu and Correll2011), which are major contributors to natural deaths. The long-term influence of antipsychotics on survival outcomes is further complicated by findings that symptom remission without antipsychotic medication and recovery during the follow-up period can lower mortality risk (Harrow, Jobe, Faull, & Yang, Reference Harrow, Jobe, Faull and Yang2017). Therefore, understanding the nuanced impact of antipsychotic treatment on long-term mortality remains critical for optimizing treatment strategies for individuals with FEP over time.
Apart from these sociodemographic and clinical determinants, healthcare and social welfare systems might be important contextual factors that contribute to the varied mortality rates. Therefore, studies from different regions other than Western countries may provide valuable information to understand the longitudinal excess mortality pattern in individuals with psychosis. Hong Kong has a population of 7.5 million with a life expectancy of 85.6 in 2024, the third-highest life expectancy in the world (Life Expectancy by Country 2025, 2025). Most of the mental health services in Hong Kong, particularly for people with severe mental illness, are provided by the public health care system. The current study aims to explore long-term mortality patterns and associated factors among individuals with FEP over 20 years in Hong Kong. The specific objectives are: (1) to characterize changing patterns of mortality of individuals with FEP over 20 years; (2) to evaluate the associations between baseline sociodemographic and clinical characteristics and subsequent mortality; and (3) to assess the impact of antipsychotic medication use during the first 10 years following FEP onset on subsequent mortality outcomes. By identifying risk factors of long-term mortality, this research seeks to inform clinical and public health strategies to improve long-term survival and outcomes in individuals with FEP.
Methods
Study design and sample
This is a 25-year follow-up study of FEPs, which was originally established for a territory-wide historical control comparison study exploring outcomes of early intervention service (EIS) for FEPs in Hong Kong (Chen, Reference Chen2004; Chen et al., Reference Chen, Tang, Hui, Chiu, Lam, Law, Yew, Wong, Chung, Tso, Chan, Yip, Hung and Honer2011; Wong et al., Reference Wong, Luo, Tang, Moore, Gur, Suen, Hui, Lee, Chang, Yan, Chui, Poon, Lo, Cheung, Kan, Chen and Chan2024). Participants were those with FEP and received the EIS (the Easy Assessment Service for Young People with Psychosis service [EASY]) for the first time between 1 July 2001 to 30 June 2003 with an age range of 15–25 (N = 700), and age, gender, and diagnosis matched FEPs received standard care services (SCS) for the first time between 1st July 1998 and 30th June 2001 (N = 700). The inclusion and exclusion criteria of this cohort were comprehensively reported in the previous study (Chan et al., Reference Chan, So, Hui, Chang, Lee, Chung, Tso, Hung, Yip, Dunn and Chen2015) and summarized in the Supplemental methods. All clinical information was obtained from the medical records and electronic health records from all the public psychiatric units provided by the Hospital Authority of Hong Kong, including the Clinical Management System (CMS) and Clinical Data Analysis and Reporting System (CDARS). The CMS is a real-time, front-end system for recording patient information and clinical activities, while CDARS is a centralized, back-end system for retrospective analysis, providing anonymous data on hospitalizations, medication use, investigation results, and other healthcare service utilization, covering approximately 90% of secondary and tertiary care services (Sek et al., Reference Sek, Cheung, Choy, Wong, Tong, Fung, Fung and Ho2007). The reliability of clinical information from CMS and CDARS has been verified (Zhou, Tang, Chan, & Luo, Reference Zhou, Tang, Chan and Luo2025). Participants were excluded from the study if their identification information was invalid or unmatched between the CMS and CDARS datasets. A total of 1,389 participants (99.2%) with FEPs were included (Figure 1).
Flowchart of the study (participants’ collection and sample size of included participants for each step).
All procedures conducted in this study adhered to the Declaration of Helsinki (1975, as amended in 2008) and were approved by the Institutional Review Board of the University of Hong Kong/Hospital Authority Hong Kong West Cluster (HKU/HA HKW IRB) (IRB reference number: UW 22-613). The requirement for informed consent was waived by the institutional review boards and ethics committees. This study followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guideline (von Elm et al., Reference von Elm, Altman, Egger, Pocock, Gøtzsche and Vandenbroucke2007).
Data collection
Mortality information
All mortality information, including cause and date of death, till the end date of this study of 23 November 2023, was obtained from the CMS and verified with the coroner’s court report. Deaths that were due to accidents, suicide, and actions with unclear intentions were considered as unnatural deaths, whereas deaths due to documented medical conditions were considered as natural deaths (Harrow et al., Reference Harrow, Jobe, Faull and Yang2017). The sex- and age-standardized death rates of all-cause mortality of the general population in Hong Kong from 1998 to 2022 were obtained from the Census and Statistics Department of the Government of the Hong Kong Special Administrative Region (The Mortality Trend in Hong Kong, 1991 to 2022, 2024) for the calculation of the standardized mortality ratios.
Baseline and long-term sociodemographic and clinical information
Baseline sociodemographic and clinical information was systematically extracted by trained researchers from the CMS and written medical records. Sociodemographic variables included sex and educational attainment. Clinical information included diagnosis, age at illness onset, DUP, social functioning, admission and discharge dates of psychiatric hospitalization, and symptom severity at onset, which included positive, negative, and affective symptoms assessed with the Clinical Global Impressions–Schizophrenia (CGI-SCH) scale (Haro et al., Reference Haro, Kamath, Ochoa, Novick, Rele, Fargas, Rodríguez, Rele, Orta, Kharbeng, Araya, Gervin, Alonso, Mavreas, Lavrentzou, Liontos, Gregor and Jones2003). DUP was defined as the duration (in days) between the first manifestation of psychotic symptoms and the initiation of effective psychiatric treatment, as determined by clinicians. Participants who were hospitalized within the first month of presentation to the psychiatric services were considered hospitalized at baseline. Social functioning was evaluated using the Social and Occupational Functioning Assessment Scale (SOFAS) (Goldman, Skodol, & Lave, Reference Goldman, Skodol and Lave1992). Biweekly consensus meetings were held for quality assurance during data collection. An experienced clinician and two researchers reviewed 12 medical records to validate the data collection process. Satisfactory validity and inter-rater reliability of key factors, including DUP, functioning, and symptom severity, were reported in a previous study (Chan et al., Reference Chan, Chan, Pang, Yan, Hui, Chang, Lee and Chen2018) and briefly summarized in the Supplement methods.
Clinical information over the first 10 years of those who survived was extracted from the CDARS by trained researchers. These include dose, types, and duration of antipsychotic medications prescribed, psychiatric and non-psychiatric hospitalization durations, and physical comorbidities. The proportion of psychiatric and non-psychiatric hospitalization duration was calculated as the total duration of psychiatric or non-psychiatric hospitalization divided by the whole follow-up period within the first 10 years. The Charlson Comorbidity Index (Charlson, Carrozzino, Guidi, & Patierno, Reference Charlson, Carrozzino, Guidi and Patierno2022) was used to assess the physical comorbidities of the participants. Polypharmacy is operationally defined as the concurrent prescription of ≥2 antipsychotics for 28 days or longer. Mean square successive difference (MSSD) (von Neumann, Kent, Bellinson, & Hart, Reference von Neumann, Kent, Bellinson and Hart1941), and mean of monthly defined daily dose (DDD) (Nosè et al., Reference Nosè, Tansella, Thornicroft, Schene, Becker, Veronese, Leese, Koeter, Angermeyer and Barbui2008) of antipsychotics and monthly treatment duration of antipsychotics were calculated. MSSD estimates variability by calculating the mean of the squared differences between consecutive data points. Details on the definitions and calculations of clinical factors for the first 10 years are summarized in Supplementary Table 1.
Statistical analysis
For the objective of the study, standard mortality ratios (SMRs) of all-cause were calculated by dividing the observed number of deaths in each gender and age group of the study sample by the expected number of deaths based on the population mortality rates in the Hong Kong census survey (Tsai, Hardy, & Wen, Reference Tsai, Hardy and Wen1992). As the age group was categorized in 5-year increments in the Hong Kong census, the consecutive 5-year average SMRs with 95% confidence intervals (CIs) were calculated and plotted for this study from 1998 to 2022. The causes of mortality over the follow-up period were described. The Aalen–Johansen estimator, a standard nonparametric method for estimating cumulative incidence functions in competing risks, was used to compare trends in natural and unnatural deaths over the study period (Aalen & Johansen, Reference Aalen and Johansen1978; Andersen, Borgan, Gill, & Keiding, Reference Andersen, Borgan, Gill and Keiding2012). As a sensitivity analysis, SMR trends and causes of mortality over the follow-up period were examined in the subgroup of individuals diagnosed with SSD at baseline.
To examine the second objective of the study, multivariable Cox proportional hazards regression models were employed to explore associations between sociodemographic and clinical characteristics at baseline and all-cause, natural, and unnatural mortality, with Bonferroni correction applied for multiple comparisons. Group differences in baseline characteristics between deceased and surviving individuals over the entire follow-up period were assessed, with continuous variables analyzed using the Monte Carlo simulation for t-test, to account for skewed distributions and small sample sizes (Sokal, Reference Sokal1992).
To study the impact of antipsychotic medication usage during the first 10 years on later mortality in our objective three, the analysis excluded participants who deceased within the first 10 years. Group differences in clinical data, including antipsychotic medication, were analyzed between individuals who later deceased and those who survived until the study’s endpoint. Continuous variables were again assessed using the Monte Carlo simulation for t-test. To address potential non-linear relationships among dependent variables and validate the significance of antipsychotic medication patterns, distance correlation tests (Gábor, Maria, & Nail, Reference Gábor, Maria and Nail2007) were employed to examine associations with psychiatric and non-psychiatric hospitalizations, as well as comorbidities, over the first 10 years. Multivariate Cox regression was conducted in two steps: (1) assessing the associations of monthly antipsychotic treatment duration and DDD with all-cause, natural, and unnatural mortality, adjusting for baseline education, gender, and psychiatric diagnosis, (2) reanalyzing these associations with adjustment for identified baseline risk factors. As a sensitivity analysis, a multivariate Cox regression was applied to examine the association between antipsychotic medication patterns and mortality within the SSD subgroup. Bonferroni correction was applied to Cox regression analyses for multiple comparisons across all-cause, unnatural, and natural mortality outcomes. A two-sided P-value <0.05 was considered statistically significant. All statistical analyses were performed using R, version 4.3.2 (The R Foundation).
Results
The study sample consisted of 1389 FEPs at baseline, with 967 (69.6%) participants having a diagnosis of SSD and a median observation period of 22 years. The baseline median age was 21.6 (interquartile range, IQR = 5.42) and 51.4% of the sample were male. To preliminarily validate our antipsychotic medication pattern metrics (mean and MSSD of monthly DDD and treatment duration, and proportion of polypharmacy duration), we evaluated their associations with hospitalization and comorbidity outcomes during the first 10 years (Supplementary Figure 1), highlighting links between antipsychotic patterns and psychiatric hospitalizations.
Mortality patterns over 20 years
Among the 137 (9.86%) deceased individuals during the follow-up period, 102 of them were due to unnatural causes and 35 were due to natural causes, with details of cause of death shown in Supplementary Table 2. The overall SMR was 6.56 (95% CI, 5.50–7.71). The highest all-cause SMR was 11.82 (95% CI, 7.73–15.92) in the first 5 years, which dropped continuously to 3.33 (95% CI, 2.00–4.66) in the last five years (22 years on average after the initial onset) (Figure 2). The cumulative incidence rate of natural deaths increased after 10 years, while that of unnatural deaths rose sharply during the first decade before stabilizing (Figure 3). Similar SMR trends were observed (Supplementary Figure 2) among the individuals with SSD, yielding an overall SMR of 6.98 (95% CI, 5.82–8.25) with similar patterns of SMRs over time. The cumulative incidence patterns for natural and unnatural deaths of SSD were reported in Supplementary Figure 3.
The trend for standardized mortality rates over the years.
(a) This figure indicates the trend of the 5-year average standardized mortality ratio (SMR) of all-cause mortality over twenty-five years. Since the whole cohort started in 1998, the first 5-year average SMR was calculated by using the general Hong Kong population in the matched age-gender specific groups in 2002, and the following 5-year average SMRs were calculated using the mortality rates of the general Hong Kong population in the matched specific age-gender groups in 2007, 2012, 2017, and 2022. (b) This table shows the exact SMR values for each calculated period.
The cumulative incidence of natural and unnatural deaths over the follow-up period.
Baseline factors associated with long-term mortality
Compared with those survived at the end of the study, individuals with FEP deceased during the follow-up period were more likely to be male (alive vs death: 50.16% vs 65.69%, p < 0.01), with a high proportion of low educational level (alive vs death: 73.65% vs 84.56%, p = 0.01), having hospitalization in the first month (alive vs death: 67.41% vs 79.56%, p < 0.01), and low SOFAS scores (alive vs death: 84.83% vs 96.35%, p < 0.01) at baseline (Supplementary Table 3). Multivariate Cox regression analyses found male participants (hazard ratio, HR: 1.72, 95% CI, 1.20–2.45) and low SOFAS scores at baseline (HR: 3.26, 95% CI, 1.29–8.23) were associated with increased risks of all-cause mortality; male participants (HR: 1.91, 95% CI, 1.26–2.90), participants with a low educational level (HR: 2.07, 95% CI, 1.17–3.67), and hospitalization in the first month had increased risks of unnatural mortality after Bonferroni correction (Table 1).
Multivariate Cox regression analysis of baseline predictors of mortality (N = 1389)
Note: Bold indicates significant associations after Bonferroni correction, *denotes P < 0.017 (the P-value threshold after Bonferroni correction).
Abbreviations: CGI, Clinical Global Impressions Scale; DUP, Duration of untreated psychosis.
a Educational levels at baseline were categorized based on years of education: low educational level (<12 years) and high educational level (≥12 years).
b The Social and Occupational Functioning Assessment Scale (SOFAS) category at Month 1 was classified based on SOFAS scores: low scores (0–55) = 1, high scores (>55) = 2 (assessed using the SOFAS scale).
c Diagnoses at baseline were categorized as follows: “Schizophrenia spectrum disorder” included schizophrenia and schizoaffective disorder; “Other psychotic disorders” encompassed acute and transient psychotic disorder (ATPD) and psychosis not otherwise specified (NOS); “Affective psychosis” included bipolar affective disorder (BAD)/mania with psychotic symptoms and severe depression with psychotic symptoms.
d Physical disease history at baseline was determined from hospital records predating the baseline, excluding psychiatric diagnoses (ICD-9 codes 219–319).
Association of antipsychotic medication use patterns of the first 10 years and mortality after first 10 years
After excluding those who had deceased within the first 10 years, 1,310 participants were eligible for examining the associations between long-term antipsychotic medication and the risk of later mortality (Figure 1). A total of 74 participants deceased during the follow-up period beyond the first 10 years. Among participants who survived the first 10 years, those who later deceased had a higher average monthly antipsychotics duration (alive vs deaths: 19.88 vs 21.18, p = 0.01, Cohen’s d = 0.29) during the first 10 years, and those deceased from natural causes had a significantly higher average monthly DDD of antipsychotics during the first 10 years (unnatural deaths vs natural deaths: 0.44 vs 0.79, p = 0.04, Cohen’s d = 0.57) than those who deceased from unnatural causes (Supplementary Table 4). Multivariate Cox regression analyses revealed that higher MSSD of monthly antipsychotic treatment duration in the first 10 years was associated with an increased risk of later all-cause mortality (HR 1.01, 95% CI, 1.00–1.03, Bonferroni-corrected; Table 2). These findings persisted after additional adjustment for baseline risk factors, including sex, education level, hospitalization status, SOFAS score categories, and baseline psychiatric diagnosis (Supplementary Table 5). Similar patterns were observed among participants with SSD (Supplementary Table 6). In this subsample, greater MSSD of monthly antipsychotic duration was linked to increased all-cause mortality (HR 1.01, 95% CI, 1.00–1.03) after Bonferroni correction.
Multivariate Cox regression analysis of the impact of antipsychotic patterns of the first 10 years of follow-ups on the mortality for the following period
Note: A total of 1310 participants were included in this analysis; models were adjusted with gender, educational levels, and diagnosis at baseline. Bold indicates significant associations after Bonferroni correction, * denotes P < 0.017 (the P-value threshold after Bonferroni correction.
Abbreviations: DDD, defined daily dose; MSSD, Mean square successive difference.
a Polypharmacy was defined as the concurrent use of two or more antipsychotics for more than 28 days. The proportion of polypharmacy duration represents the total time a patient spent on multiple antipsychotics divided by their total duration of antipsychotic treatment.
Discussion
This study investigated the mortality patterns of 1,389 individuals with FEP followed for up to 25 years (median of 22 years), with about 70% of the study population having an SSD diagnosis at baseline. Throughout the follow-up period, 137 participants (9.86%) deceased, giving an overall SMR of 6.56. The cumulative incidence rate of natural deaths increased after the first decade of the illness, while that of unnatural deaths rose sharply during the first decade and stabilized in the second decade of the illness. Males with lower educational levels, poorer social functioning, and being hospitalized at baseline were associated with increased risk of mortality over 20 years. Greater fluctuations in monthly antipsychotic medication duration during the first 10 years were linked to an increased risk of subsequent all-cause mortality.
The all-cause long-term mortality rate observed in our study was 9.86%, consistent with findings from prior longitudinal studies, which reported mortality rates ranging between 10% and 14% among individuals with psychotic disorders (Crump et al., Reference Crump, Winkleby, Sundquist and Sundquist2013; Dickerson et al., Reference Dickerson, Khan, Origoni, Rowe, Katsafanas, Harvin, Yang and Yolken2024; Hor & Taylor, Reference Hor and Taylor2010; Starzer et al., Reference Starzer, Hansen, Hjorthøj, Speyer, Albert and Nordentoft2023). However, the overall SMR of all-cause mortality (6.56) and even the lowest SMR of all-cause mortality during the last 5 years of follow-up period (3.33) of this study population, as well as the subsample of SSD (overall all-cause SMR = 6.98, lowest all-cause SMR = 3.62), were higher than the median SMR of all-cause mortality of schizophrenia reported in a systematic review (2.56) (Saha, Chant, & McGrath, Reference Saha, Chant and McGrath2007). The relatively young study population at the start of the study might be a possible reason. The high life expectancy of Hong Kong might also explain this bigger excess mortality gap found in the study, reflecting more prominent health inequalities in people with psychotic disorders. Similar to the 20-year follow-up of OPUS study (Starzer et al., Reference Starzer, Hansen, Hjorthøj, Speyer, Albert and Nordentoft2023), our study found that suicide is the leading cause of death in the FEP population (70.8%), particularly during the first 10 years of illness onset. The significantly higher risk of suicide death during the early stage of the illness reported in this study, as well as previous studies (Brown, Reference Brown1997; Melle et al., Reference Melle, Olav Johannesen, Haahr, Ten Velden Hegelstad, Joa, Langeveld, Larsen, Ilner Opjordsmoen, Qin, Ivar Røssberg, Rishovd Rund, Simonsen, Vaglum, McGlashan and Friis2017), highlights the crucial need to establish and implement effective interventions to reduce suicide during the early stage of the psychotic illness. Even though EIS for psychosis has been reported to be effective in reducing suicide death (Chai et al., Reference Chai, Tang, Ma, Luo and Chan2024; Chan et al., Reference Chan, Chan, Pang, Yan, Hui, Chang, Lee and Chen2018) and many countries and regions have already implemented the EIS (Correll et al., Reference Correll, Solmi, Veronese, Bortolato, Rosson, Santonastaso, Thapa-Chhetri, Fornaro, Gallicchio, Collantoni, Pigato, Favaro, Monaco, Kohler, Vancampfort, Ward, Gaughran, Carvalho and Stubbs2017), the excess mortality due to suicide of psychosis is still prominent. Therefore, establishing further effective intervention strategies to reduce suicide within the EIS context for people with FEP is still crucial. Identifying risk populations, possibly with algorithm-based approaches (Lieslehto et al., Reference Lieslehto, Tiihonen, Lähteenvuo, Leucht, Correll, Mittendorfer-Rutz, Tanskanen and Taipale2024; Wong et al., Reference Wong, Luo, Tang, Moore, Gur, Suen, Hui, Lee, Chang, Yan, Chui, Poon, Lo, Cheung, Kan, Chen and Chan2024), to provide targeted monitoring and intervention might be one possible strategy.
The current study revealed a change in the mortality pattern over the follow-up period. There was a gradual decline of SMR of all-cause mortality in the first 15 years of illness onset, with a plateau of incidence rate of unnatural cause of death (mostly suicide), and then a sharp reduction of SMR of all-cause mortality with an increase of natural death between 15 and 25 years of follow-up. This mortality pattern was also seen in previous long-term follow-up studies in Australia (Yuen et al., Reference Yuen, Harrigan, Mackinnon, Harris, Yuen, Henry, Jackson, Herrman and McGorry2014) and Denmark (Starzer et al., Reference Starzer, Hansen, Hjorthøj, Speyer, Albert and Nordentoft2023). About one quarter of the natural deaths in this study were due to cardiovascular conditions. This highlighted the importance of physical health care among individuals with FEP to reduce risk factors related to cardiovascular conditions, particularly after the initial few years of illness (Correll et al., Reference Correll, Solmi, Veronese, Bortolato, Rosson, Santonastaso, Thapa-Chhetri, Fornaro, Gallicchio, Collantoni, Pigato, Favaro, Monaco, Kohler, Vancampfort, Ward, Gaughran, Carvalho and Stubbs2017). One of which is preventing the development of metabolic syndrome during the early stage of psychosis in the young population. Using a well-validated risk calculator for the development of metabolic syndrome in the young FEP population (Perry et al., Reference Perry, Osimo, Upthegrove, Mallikarjun, Yorke, Stochl, Perez, Zammit, Howes, Jones and Khandaker2021; Tse et al., Reference Tse, Khandaker, Zhou, Luo, Yan, Siu, Poon, Lee, Zhang, Upthegrove, Osimo, Perry and Chan2024) may be an important strategy to reduce the risk of long-term mortality. Another important reason for the excess mortality and morbidities due to physical health has been suggested to be the disparities in physical health access, utilization, and provision (Lawrence & Kisely, Reference Lawrence and Kisely2010; Nasrallah et al., Reference Nasrallah, Meyer, Goff, McEvoy, Davis, Stroup and Lieberman2006; Osborn, King, & Nazareth, Reference Osborn, King and Nazareth2003). Therefore, integrating physical health care within the mental health service provision to allow streamlined service delivery should be considered.
Our study identified male gender, low educational attainment, poor social functioning, and hospitalization at baseline as significant factors associated with long-term all-cause and unnatural mortality in individuals with FEP. Given that most of the unnatural causes of mortality in this study were suicide (95%), and the same set of factors were found to be significantly different between the natural and unnatural deaths, indicating these factors are likely to be risk factors of suicide in this population. These factors are similar to findings of previous studies (Björkenstam et al., Reference Björkenstam, Björkenstam, Hjern, Bodén and Reutfors2014; Hawton et al., Reference Hawton, Sutton, Haw, Sinclair and Deeks2005; Robinson et al., Reference Robinson, Harris, Harrigan, Henry, Farrelly, Prosser, Schwartz, Jackson and McGorry2010), including the OPUS 20-year follow-up study (Starzer et al., Reference Starzer, Hansen, Hjorthøj, Speyer, Albert and Nordentoft2023). In fact, male gender, low educational attainment, and poorer baseline social functioning have generally been considered as risk factors for poorer outcomes of people with psychosis (Santesteban-Echarri et al., Reference Santesteban-Echarri, Paino, Rice, González-Blanch, McGorry, Gleeson and Alvarez-Jimenez2017). Hospitalization at baseline likely reflects higher illness severity at initial presentation. A machine learning-based model in the Swedish and Finnish population has identified hospitalization as a key risk factor of mortality in first-episode psychosis (FEP) (Lieslehto et al., Reference Lieslehto, Tiihonen, Lähteenvuo, Leucht, Correll, Mittendorfer-Rutz, Tanskanen and Taipale2024). Substance use disorders (SUD) at baseline were not identified as a significant risk factor for later mortality during the follow-up period in this study. This may be attributed to the low prevalence of SUD in this FEP population, which reflects the low general SUD in the Hong Kong population (Lau, Kim, & Tsui, Reference Lau, Kim and Tsui2005).
Some risk factors that contribute to the excess mortality and morbidity due to physical illness in the psychosis population may be modifiable, including unhealthy lifestyle and the use of antipsychotics (Correll et al., Reference Correll, Solmi, Veronese, Bortolato, Rosson, Santonastaso, Thapa-Chhetri, Fornaro, Gallicchio, Collantoni, Pigato, Favaro, Monaco, Kohler, Vancampfort, Ward, Gaughran, Carvalho and Stubbs2017; Saha et al., Reference Saha, Chant and McGrath2007). Our study found that greater fluctuations of antipsychotic usage in the first 10 years were associated with higher subsequent mortality risks, with half of them due to physical causes. Inconsistent use of antipsychotics may reflect a course of recurring illness, which will increase vulnerability to adverse events such as suicide and accidents, and poor functioning, which then contributes to an unhealthy lifestyle (Haddad, Brain, & Scott, Reference Haddad, Brain and Scott2014; Tiihonen et al., Reference Tiihonen, Lönnqvist, Wahlbeck, Klaukka, Niskanen, Tanskanen and Haukka2009). Furthermore, fluctuations in antipsychotics intensify cardiovascular risks, including metabolic syndrome, QT-interval abnormalities, and diabetes-related complications, contributing to natural mortality (De Hert et al., Reference De Hert, Detraux, van Winkel, Yu and Correll2011; Henderson et al., Reference Henderson, Vincenzi, Andrea, Ulloa and Copeland2015; Ray et al., Reference Ray, Chung, Murray, Hall and Stein2009). Additionally, irregular treatment elevates respiratory and infectious deaths through immune dysfunction and smoking-related health deterioration (Miller et al., Reference Miller, Buckley, Seabolt, Mellor and Kirkpatrick2011). Though the reason for unstable use of antipsychotics could not be specified, our cross-sectional analysis on the relationship between antipsychotic patterns and psychiatric hospitalization over the first 10 years revealed that fluctuation of the duration of antipsychotic usage was significantly related to both psychiatric and nonpsychiatric hospitalization. These results highlighted the importance of consistency of antipsychotic usage and thus illness stability in relation to long-term mortality.
Several limitations of the study should be considered in interpreting our findings. First, the relatively small sample size might limit the power to identify specific risk factors for some causes of death, such as suicide. Second, interpreting longitudinal SMR trends may be limited by this cohort’s fixed design, which lacks new entrants over time, and this may cause survivor bias that could attenuate SMR estimates as the remaining population reflects stronger survivors. Third, the associations observed between baseline risk factors, antipsychotic treatment, and mortality may be confounded by unmeasured individual biopsychosocial variables, such as genetic predispositions, childhood trauma, and family history of self-harm (Hor & Taylor, Reference Hor and Taylor2010). Meanwhile, all the information of the study was obtained from the clinical records of the public mental health services. Information on participants who had emigrated or stopped receiving care from the public health services during the follow-up period will not be available, leading to bias. The quality of the information obtained depends on the quality of clinical documentation. Furthermore, only baseline diagnoses were considered without examining diagnostic transitions during follow-up. This may underestimate or obscure the distinct mortality risks of types of psychotic disorders. Finally, the cohort was drawn from a specific geographic and healthcare context, which may limit the generalizability of the results.
Conclusions
This study highlights the persistent excess mortality among individuals with FEP despite advances in healthcare, revealing a distinct temporal shift: natural deaths increase significantly after 10 years, while unnatural deaths peak during the first decade and stabilize thereafter. Male gender, low educational attainment, poor social functioning, and hospitalization at baseline were significantly associated factors of mortality. High fluctuations in antipsychotic medication in the first 10 years increased subsequent mortality risk among FEP. Results underscore the need for individualized and integrated mental and physical health care services to reduce long-term mortality in patients with FEP.
Supplementary material
The supplementary material for this article can be found at http://doi.org/10.1017/S0033291725102286.
Acknowledgements
The authors thank all the psychiatric units in Hong Kong for logistic support and the patients who participated in the study.
Funding statement
This work was supported by the General Research Fund of Hong Kong (Sherry Kit Wa Chan, grant reference numbers 17106021).
Competing interests
The funding institutions of this study had no further role in the study design, the collection, analysis, and interpretation of data, the writing of this paper, or the decision to submit for publication. The authors report no financial relationships with commercial interests.
Open access
