Introduction
Bipolar disorder is a severe mental disorder which affects more than 1% of the population and constitutes one of the main causes of disability worldwide (McIntyre et al., Reference McIntyre, Berk, Brietzke, Goldstein, López-Jaramillo, Kessing, Malhi, Nierenberg, Rosenblat, Majeed, Vieta, Vinberg, Young and Mansur2020). People with bipolar disorder have markedly elevated risk of premature mortality (Hayes et al., Reference Hayes, Miles, Walters, King and Osborn2015), with approximately 8–15-years shorter lifespan relative to the general population (Chang et al., Reference Chang, Hayes, Perera, Broadbent, Fernandes, Lee, Hotopf and Stewart2011; Laursen, Reference Laursen2011; Kodesh et al., Reference Kodesh, Goldshtein, Gelkopf, Goren, Chodick and Shalev2012; Ajetunmobi et al., Reference Ajetunmobi, Taylor, Stockton and Wood2013; Crump et al., Reference Crump, Sundquist, Winkleby and Sundquist2013; Kessing et al., Reference Kessing, Vradi and Andersen2015; Pan et al., Reference Pan, Yeh, Chan and Chang2020). Excess deaths associated with bipolar disorder are mainly attributed to natural causes, particularly cardiovascular diseases, respiratory diseases and cancers (Crump et al., Reference Crump, Sundquist, Winkleby and Sundquist2013; Hayes et al., Reference Hayes, Miles, Walters, King and Osborn2015). Accumulating evidence indicates that such mortality gap has persisted (Ajetunmobi et al., Reference Ajetunmobi, Taylor, Stockton and Wood2013; Crump et al., Reference Crump, Sundquist, Winkleby and Sundquist2013) or widened (Hayes et al., Reference Hayes, Marston, Walters, King and Osborn2017; Pan et al., Reference Pan, Yeh, Chan and Chang2017; Hansen et al., Reference Hansen, Laursen, Grøntved, Straszek, Licht and Nielsen2019; Lomholt et al., Reference Lomholt, Andersen, Sejrsgaard-Jacobsen, Ozdemir, Graff, Schjerning, Jensen, Straszek, Licht, Grontved and Nielsen2019) in recent decades despite overall improvement in life expectancy in the general population due to enhanced healthcare. Physical health inequalities experienced by people with bipolar disorder thus represent a serious public health challenge that warrants urgent attention. Comprehensive evaluation of premature mortality patterns associated with bipolar disorder is crucial for developing effective strategies and optimising healthcare service delivery to reduce avoidable deaths in this vulnerable population.
Notably, literature on excess mortality in bipolar disorder were mostly derived from Western countries (Hayes et al., Reference Hayes, Miles, Walters, King and Osborn2015) and their findings may not be generalisable to other areas owing to substantial cross-regional variation in healthcare systems, sociocultural context and population health indices. Until now, very few studies have been conducted in Asia in this respect (Saku et al., Reference Saku, Tokudome, Ikeda, Kono, Makimoto, Uchimura, Mukai and Yoshimura1995; Pan et al., Reference Pan, Yeh, Chan and Chang2017, Reference Pan, Yeh, Chan and Chang2020) and were limited by small sample size and short follow-up duration. In general, previous research was hampered by several important methodological constraints including sampling patients with psychiatric inpatient treatment only with subsequent bias towards greater illness severity (Saku et al., Reference Saku, Tokudome, Ikeda, Kono, Makimoto, Uchimura, Mukai and Yoshimura1995; Ösby et al., Reference Ösby, Brandt, Correia, Ekbom and Sparén2001; Ajetunmobi et al., Reference Ajetunmobi, Taylor, Stockton and Wood2013; Hoang et al., Reference Hoang, Goldacre and Stewart2013; Laursen et al., Reference Laursen, Wahlbeck, Hällgren, Westman, Osby, Alinaghizadeh, Gissler and Nordentoft2013), lack of evaluation for cause-specific mortality (Chang et al., Reference Chang, Hayes, Broadbent, Fernandes, Lee, Hotopf and Stewart2010; Kodesh et al., Reference Kodesh, Goldshtein, Gelkopf, Goren, Chodick and Shalev2012; Ajetunmobi et al., Reference Ajetunmobi, Taylor, Stockton and Wood2013; Medici et al., Reference Medici, Videbech, Gustafsson and Munk-Jørgensen2015; Pan et al., Reference Pan, Yeh, Chan and Chang2017; Hansen et al., Reference Hansen, Laursen, Grøntved, Straszek, Licht and Nielsen2019; Lomholt et al., Reference Lomholt, Andersen, Sejrsgaard-Jacobsen, Ozdemir, Graff, Schjerning, Jensen, Straszek, Licht, Grontved and Nielsen2019) or focus only on one specific natural cause of death for analysis (Saku et al., Reference Saku, Tokudome, Ikeda, Kono, Makimoto, Uchimura, Mukai and Yoshimura1995; Laursen et al., Reference Laursen, Wahlbeck, Hällgren, Westman, Osby, Alinaghizadeh, Gissler and Nordentoft2013; Westman et al., Reference Westman, Hällgren, Wahlbeck, Erlinge, Alfredsson and Ösby2013; Hayes et al., Reference Hayes, Marston, Walters, King and Osborn2017). Alternatively, a growing number of studies have assessed the impact of premature mortality on survival in bipolar disorder using years of life lost (YLLs), which are based on estimating life expectancy at a single fixed age (mostly at birth or 15 years). This approach, however, does not reflect the underlying age-of-onset distribution of the disorder and implies an assumption that the estimated life expectancy is interpreted as the expected lifespan for patients who have had the disorder since the set fixed age (Plana-Ripoll et al., Reference Plana-Ripoll, Canudas-Romo, Weye, Laursen, McGrath and Andersen2020). In fact, recent Danish register-based studies which investigated the longevity gap between people with mental disorders and the general population using life-years lost (LYLs; Erlangsen et al., Reference Erlangsen, Andersen, Toender, Laursen, Nordentoft and Canudas-Romo2017; 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, McGrath and Laursen2019), a novel and more realistic measure taking into consideration different ages of onset of the disorder (Andersen, Reference Andersen2017), revealed that past estimates of reduced life expectancy tend to overestimate mortality difference. Thus far, there has been no study examining LYLs for bipolar disorder.
In this large population-based cohort study, we aimed to comprehensively examine the risk of premature mortality associated with bipolar disorder over 11 years in Hong Kong (HK), a metropolitan city located at the southeastern tip of China with a population of approximately 7.5 million, utilising data retrieved from a territory-wide medical-record database of public healthcare services. Specifically, we adopted two mortality metrics, namely SMR (for all-cause and cause-specific mortality) and LYLs to quantify the magnitude of excess mortality among patients with bipolar disorder compared with the general population. Changes in SMRs across the study period were also assessed to clarify whether the mortality gap improved or worsened over time.
Methods
Data source
Population statistics and information on all registered deaths in HK between 2008 and 2018 (each calendar year from 2008 to 2018, inclusive) were obtained from the Census and Statistics Department. Data of the patient cohort were extracted from the Clinical Data Analysis and Reporting System (CDARS; Hospital Authority, 2003), a territory-wide electronic health-record database developed by the Hospital Authority (HA) which is a statutory body delivering government-subsidised, universal health coverage to all HK residents (approximately 92% being Chinese) by managing all public hospitals, specialist and general outpatient clinics in HK. Detailed description of CDARS has been reported elsewhere (Cheung et al., Reference Cheung, Fung, Wong, Tong, Sek, Greyling, Tse and Fung2007). Briefly, CDARS is an integrated, longitudinal patient electronic record system capturing clinical data across all healthcare settings of HA facilities. The database contains patients’ demographics and clinical information including diagnoses, attendances to outpatient clinics and emergency departments and hospital admissions. Data on dates and causes of death were retrieved from CDARS via internal linkage to regional death registries from the Immigration Department. Clinical data are collected and entered into computerised clinical-management system (CMS) by treating clinicians and other healthcare professionals, and are then transferred to CDARS for audit and research purposes. CDARS generates unique, anonymised patient identifiers to protect privacy and to link all medical records. This database has been used to conduct high-quality population-based studies on various psychiatric disorders including schizophrenia and other psychoses (Chang et al., Reference Chang, Chan, Wong, Hai, Or and Chen2020; Yung et al., Reference Yung, Wong, Chan, Or, Chen and Chang2020, Reference Yung, Wong, Chan, Chen and Chang2021; Chan et al., Reference Chan, Wong, Or, Chen and Chang2021).
Study population
The study period was between 1 January 2008 and 31 December 2018. We identified all individuals who were diagnosed with bipolar disorder for public psychiatric inpatient admissions or outpatient care between 1 January 2002 and 31 December 2018 (computerised CMS for psychiatric services implemented since 1 January 2000), and aged ⩾15 years during the study period as the study population. Diagnosis of bipolar disorder was recorded and verified by the International Classification of Diseases, 10th revision (ICD10 codes: F30 and F31). Final diagnostic ascertainment took into consideration the longitudinal illness course (Chang et al., Reference Chang, Pang, Chung and Chan2009), and those identified patients with subsequent diagnostic change to schizophrenia or schizoaffective disorder (as their most-recently assigned principal diagnosis) before the end of follow-up were excluded. Follow-up of the patient cohort began on the date of first-recorded diagnosis of bipolar disorder. For patients who had been assigned with a diagnosis of bipolar disorder before the study period, their follow-up start-date was defined as 1 January 2008. The cohort was followed forward until the date of death or 31 December 2018, whichever came first. The study was approved by the Institutional Review Board of the University of Hong Kong/Hospital Authority Hong Kong West Cluster. The study data were anonymised and individual patient records were completely unidentifiable during the analysis. Since our study was based on medical-record database, the requirement for informed consent was waived.
Study outcomes
Causes of death were classified according to ICD10 codes (online Supplementary Table S1), and were divided into natural and unnatural causes. Natural causes were categorised into infectious and parasitic diseases (A00–B99), neoplasms (C00–D48), cardiovascular diseases (I00–I99), respiratory diseases (J00–J99), digestive diseases (K00–K93) and genitourinary diseases (N00–N99). An array of specific natural causes was also identified for analyses. As data on specific unnatural causes (V01–Y98) were not available, we treated unnatural deaths as a single category for analyses.
Statistical analysis
Standardised mortality ratios (SMRs) for all-cause and cause-specific deaths were calculated as the primary mortality measures to quantify the relative mortality rate between patients with bipolar disorder and the general population. First, number of observed deaths and person-years of follow-up were computed for each calendar year (2008–2018), sex and age category (by 5-year age bands from 15 to 84 years, and ⩾85 years) for the patient group. The number of person-years for each stratum was multiplied by the corresponding mortality rate in the general population to produce expected number of deaths, indirectly standardising overall mortality ratio by age, sex and calendar year. SMRs were estimated by dividing the observed number of deaths by the expected number of deaths. Crude mortality rates (CMRs) per 100 000 person-years as well as SMRs for all-cause and cause-specific deaths of the patient group over the whole study period, stratified by sex and four broader age groups (15–34, 35–49, 50–64 and ⩾65 years), were then calculated, with 95% confidence intervals (CIs) of SMRs being derived by mid-P exact tests. To assess trends in all-cause, natural-cause and unnatural-cause SMRs over time, joinpoint regression analysis was performed which estimated the optimal number of linear slopes and joinpoints with modified Bayesian information criteria (Kim et al., Reference Kim, Fay, Feuer and Midthune2000). The models incorporated estimated variation for each data point using the standard error of SMRs, and calculated annual percentage change (APC) in SMRs with 95% CIs.
Excess life-years lost (LYLs) for all-cause mortality were also calculated as a complementary mortality measure. Following the method adopted in previous research for LYL estimation (Andersen, Reference Andersen2017; Erlangsen et al., Reference Erlangsen, Andersen, Toender, Laursen, Nordentoft and Canudas-Romo2017; 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, McGrath and Laursen2019; Yung et al., Reference Yung, Wong, Chan, Chen and Chang2021), we first computed average life expectancy at first-recorded diagnosis of bipolar disorder (as a proxy for age of onset), which takes into account varying ages at illness onset. Life expectancies of patients were calculated for every age at diagnosis from 15 years until a set upper-age limit (95 years in the current study). These were then combined into average life expectancy until age 95 years, weighted by the number of patients at each of the particular ages at diagnosis (Erlangsen et al., Reference Erlangsen, Andersen, Toender, Laursen, Nordentoft and Canudas-Romo2017; 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, McGrath and Laursen2019). LYLs of patients referred to the average number of years lost due to death and was computed by subtracting the average life expectancy from 95 years (i.e. remaining 95-year restricted life expectancy). Similarly, in general population, LYLs denoted the average number of years lost before 95 years for a group of people of the same sex who were alive at ages corresponding to the age-of-onset distribution of those patients with bipolar disorder. Excess LYLs for men and women with bipolar disorder were then calculated as the difference in LYLs between patients and the general population. Thus, excess LYLs referred to the average number of years that patients lost in excess of that observed in the general population of the same sex and age (Andersen, Reference Andersen2017; Erlangsen et al., Reference Erlangsen, Andersen, Toender, Laursen, Nordentoft and Canudas-Romo2017). The 95% CIs of excess LYLs were derived from non-parametric bootstrap with 500 iterations. We also decomposed total LYLs for patients into loss attributable to specific causes of death based on a completing risk model to generate cause-specific LYLs (Andersen, Reference Andersen2013). This enables clarification of how specific death causes contributing to premature mortality in bipolar disorder. Excess LYLs were estimated using the R package lillies (Plana-Ripoll et al., Reference Plana-Ripoll, Canudas-Romo, Weye, Laursen, McGrath and Andersen2020).
Results
The study population included 12 556 patients with bipolar disorder (men: 4928; women: 7628) with 106 147 persons-years of follow-up (men: 40 739; women: 65 408), and a total of 1042 deaths, of which 821 (78.8%) had a known cause. The numbers of person-years and cause-specific death counts for each demographic subgroup of patients are listed in online Supplementary Table S2.
All-cause mortality rates and excess LYLs
All-cause SMR for bipolar disorder was significantly increased in the total sample (SMR = 2.60, 95% CI: 2.45–2.76) and in each demographic subgroup relative to the general population (Tables 1 and 2). All-cause SMR for women and for men was 2.82 (95% CI: 2.59–3.06) and 2.38 (95% CI: 2.18–2.60), respectively. It decreased with increasing age and was particularly high in the youngest age-group (SMR = 10.18, 95% CI: 7.94–12.71). Bipolar disorder was associated with shorter remaining life expectancy relative to the general population. Estimated LYLs for men with bipolar disorder was 19.90 years compared with 13.12 years for men in the general population, with a difference of 6.78 (95% CI 6.00–7.84) excess LYLs. Women with bipolar disorder and in the general population had LYLs of 15.86 and 8.51 years, respectively, with a difference of 7.35 (95% CI 6.75–8.06) excess LYLs.
Crude mortality rates, all-cause and cause-specific standardised mortality ratios of patients with bipolar disorder
CI, confidence interval; CMR, crude mortality rate; SMR, standardised mortality ratio (standardised for age, sex and calendar year).
a CMRs are presented as the numbers of deaths per 100 000 person-years.
Crude mortality rates, all-cause and cause-specific standardised mortality ratios of patients with bipolar disorder by age group
CI, confidence interval; CMR, crude mortality rate; COPD, chronic obstructive pulmonary disease, SMR, standardised mortality ratio (standardised for age, sex and calendar year).
a CMRs are presented as the numbers of deaths per 100 000 person-years.
Natural-cause and unnatural-cause mortality rates
SMRs for natural and unnatural causes were significantly increased in the total sample (natural-cause SMR = 1.90, 95% CI: 1.76–2.05; unnatural-cause SMR = 8.63, 95% CI: 7.34–10.03) and in each demographic subgroup of patients with bipolar disorder, except the youngest age-group who showed no significant increase in natural-cause SMR (Tables 1 and 2). Both natural-cause and unnatural-cause SMRs were higher for women than for men, and generally decreased with age, with markedly elevated unnatural-cause SMR in the youngest age-group (16.58, 95% CI: 12.14–21.72).
Natural causes accounted for most of the known-cause deaths in patients. Approximately two-thirds of all known-cause deaths were attributed to respiratory diseases, cancers and cardiovascular diseases (online Supplementary Table S2). Respiratory diseases represented the leading cause of death and accounted for around one-thirds of all known-cause deaths. Cancers and cardiovascular diseases contributed to about 1 in 5 and 1 in 7 known-cause deaths, respectively. Unnatural causes accounted for approximately one-fifths of known-cause deaths. Men had generally higher CMRs for most listed causes of death than women (Table 1). Mortality rate for most natural causes increased with age, while that for unnatural causes was comparatively higher in the youngest and the oldest age-groups (Table 2). Analyses on cause-specific LYLs for bipolar disorder generally revealed consistent results. Respiratory diseases accounted for the largest share of LYLs for patients of both sexes and the degree of contribution increased with age, particularly in men. Cancers and cardiovascular diseases also constituted a major share to LYLs due to natural-cause deaths for both men and women with bipolar disorder, with the proportion of their contributions mildly increased with age. By contrast, unnatural causes accounted for a decreasing share of LYLs with age in patients of both sexes (online Supplementary Fig. S3 for detailed results).
As shown in Table 1, respiratory diseases displayed the highest SMR among natural-cause categories (4.09, 95% CI: 3.63–4.59), with pneumonia (non-aspiration and aspiration) showing the highest SMRs across all listed specific natural causes. Conversely, SMR for chronic obstructive pulmonary disease (COPD) was not elevated. We observed increased SMR for cardiovascular disease (1.33, 95% CI: 1.10–1.59), with subgroup analyses revealing significantly elevated relative mortality rate among women and individuals aged 50–64 years. Women, but not men or the total sample, with bipolar disorder had elevated SMR for cancers (1.25, 95% CI: 1.01–1.51). Among various cancer types, SMR was only significantly increased for colon cancer, specifically among women and age-groups ⩾50 years (Table 1 and 2). Bipolar disorder was also associated with significantly increased SMRs for infectious and parasitic diseases, digestive diseases and genitourinary diseases.
Trends in mortality rates over study period
Joinpoint regression fitted a linear model of all-cause, natural-cause and unnatural-cause SMRs with no joinpoints over time for the total sample, men and women subgroups. As shown in Fig. 1, both the total sample and women exhibited a stable trend of all-cause SMRs over a period of 2008–2018. A trend in all-cause SMRs in men with bipolar disorder indicated a small degree of reduction over study period, with a significant APC in rate of −2.80% (95% CI: −5.24 to −0.29). There was no significant APC in natural-cause and unnatural-cause SMRs over time for the total sample or in either sex (Table 3, Supplementary Figs S1 and S2).
Annual all-cause standardised mortality ratios of patients with bipolar disorder over study period: (a) total sample, (b) men and (c) women.
Note: APC, annual percentage change; CI, confidence interval.
Annual percentage change in all-cause, natural-cause and unnatural-cause standardised mortality ratios of patients with bipolar disorder from 2008 to 2018a
APC, annual percentage change; CI, confidence interval; SMR, standardised mortality ratio.
a Joinpoint regression analyses were performed to calculate annual percentage change of SMRs with 95% CIs over the study period.
Discussion
To our knowledge, this is the first study using LYLs complementary with relative mortality risk measure (SMR) to evaluate the excess mortality associated with bipolar disorder. Our results showed that patients with bipolar disorder had 2.6-fold increased mortality risk compared with the general population. Both men and women with bipolar disorder exhibited a substantially shorter lifespan than the general population, with approximately 7 years of excess LYLs. Natural causes, specifically respiratory diseases, cardiovascular diseases and cancers, accounted for the majority of known-cause deaths, while unnatural cause had markedly elevated SMR, particularly among patients aged <35 years. Generally, the mortality gap persisted over 11 years, albeit slightly improved in men with bipolar disorder.
Our finding on all-cause SMR for bipolar disorder is broadly consistent with the literature which reported 2–3 times greater overall mortality risk than the general population (Chang et al., Reference Chang, Hayes, Broadbent, Fernandes, Lee, Hotopf and Stewart2010; Crump et al., Reference Crump, Sundquist, Winkleby and Sundquist2013; Laursen et al., Reference Laursen, Wahlbeck, Hällgren, Westman, Osby, Alinaghizadeh, Gissler and Nordentoft2013; Westman et al., Reference Westman, Hällgren, Wahlbeck, Erlinge, Alfredsson and Ösby2013), and is in line with a recent meta-analysis reporting all-cause summary SMR of 2.05 for bipolar disorder (Hayes et al., Reference Hayes, Miles, Walters, King and Osborn2015). In particular, our results are similar to the findings of Taiwan studies (the only recently published Asian studies in this respect), which revealed that all-cause SMR for bipolar disorder ranged from 2.39 to 2.65 (Pan et al., Reference Pan, Yeh, Chan and Chang2017, Reference Pan, Yeh, Chan and Chang2020). Of note, although we affirmed that bipolar disorder was associated with considerably shortened life expectancy after diagnosis, our observed magnitude of longevity gap is smaller than most previous estimates based on YLLs, ranging between 10 and 15 years shorter in lifespan than the general population (Chang et al., Reference Chang, Hayes, Perera, Broadbent, Fernandes, Lee, Hotopf and Stewart2011; Ajetunmobi et al., Reference Ajetunmobi, Taylor, Stockton and Wood2013; Laursen et al., Reference Laursen, Wahlbeck, Hällgren, Westman, Osby, Alinaghizadeh, Gissler and Nordentoft2013; Pan et al., Reference Pan, Yeh, Chan and Chang2020). Recent research, however, indicates that LYL method yields more precise, though more conservative, estimates than prior measures of reduced life expectancy by incorporating variations in age of onset of the disorder into life expectancy calculation (Erlangsen et al., Reference Erlangsen, Andersen, Toender, Laursen, Nordentoft and Canudas-Romo2017; 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, McGrath and Laursen2019). This also concurs with the results of a Danish study which demonstrated that difference in remaining life expectancy between bipolar disorder patients and the general population decreased from 12.0 to 8.7 years for men and 10.6 to 8.3 years for women, when the set age of onset for lifespan estimation increased from 25 to 45 years (Kessing et al., Reference Kessing, Vradi and Andersen2015). Taken together, adoption of LYLs enables more accurate quantification of life expectancy reduction for bipolar disorder, which in turn facilitates policy planning, resource allocation and service development for improvement in physical health outcomes.
Until now, very few studies have comprehensively examined cause-specific SMRs for natural deaths in bipolar disorder. Consistent with the literature (Crump et al., Reference Crump, Sundquist, Winkleby and Sundquist2013; Hayes et al., Reference Hayes, Miles, Walters, King and Osborn2015), we found that cardiovascular diseases, respiratory diseases and cancers represented the major contributors to natural-cause mortality in bipolar disorder. This also accords with the results of our cause-specific LYL analyses demonstrating that these three disease categories accounted for the major share of total LYLs due to natural-cause deaths for patients of both sexes. In fact, evidence has shown that bipolar disorder is associated with an increased incidence of cardiovascular diseases (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), but inequitable cardiac care including lower receipt of revascularisation interventions (Wu et al., Reference Wu, Chen, Juang, Fang, Liu, Sun, Kao, Dewey, Prince and Stewart2013; Heiberg et al., Reference Heiberg, Nesvåg, Balteskard, Bramness, Hultman, Naess, Reichborn-Kjennerud, Ystrom, Jacobsen and Hoye2020) and under-prescription of cardioprotective medications (Smith et al., Reference Smith, Martin, MacLean, Langan, Guthrie and Mercer2013; Laursen et al., Reference Laursen, Mortensen, MacCabe, Cohen and Gasse2014). Our result of SMR for cardiovascular diseases was in the lower end of those previously reported (mostly around 2-fold increased mortality risk; Ösby et al., Reference Ösby, Brandt, Correia, Ekbom and Sparén2001; Crump et al., Reference Crump, Sundquist, Winkleby and Sundquist2013; Laursen et al., Reference Laursen, Wahlbeck, Hällgren, Westman, Osby, Alinaghizadeh, Gissler and Nordentoft2013; Westman et al., Reference Westman, Hällgren, Wahlbeck, Erlinge, Alfredsson and Ösby2013; Hayes et al., Reference Hayes, Miles, Walters, King and Osborn2015). We observed that women, but not men, with bipolar disorder had raised cardiovascular-related SMR. This is contrary to most earlier research which found elevated mortality risk in both sexes (Ösby et al., Reference Ösby, Brandt, Correia, Ekbom and Sparén2001; Crump et al., Reference Crump, Sundquist, Winkleby and Sundquist2013; Westman et al., Reference Westman, Hällgren, Wahlbeck, Erlinge, Alfredsson and Ösby2013; Pan et al., Reference Pan, Yeh, Chan and Chang2020), although some revealed slightly higher cardiovascular mortality in women than in men (Ösby et al., Reference Ösby, Brandt, Correia, Ekbom and Sparén2001; Crump et al., Reference Crump, Sundquist, Winkleby and Sundquist2013). A recent study further demonstrated that female patients with bipolar disorder were more likely to die of undiagnosed cardiovascular diseases than male counterparts (Heiberg et al., Reference Heiberg, Jacobsen, Balteskard, Bramness, Naess, Ystrom, Reichborn-Kjennerud, Hultman, Nesvag and Hoye2019). Lifestyle modification, early detection and optimal treatment of comorbid diabetes, dyslipidaemia and hypertension should be implemented to reduce cardiovascular mortality in bipolar disorder. Our findings that respiratory diseases had the highest SMR across all natural-cause categories accord with past studies indicating respiratory diseases as one of the leading causes of death for bipolar disorder with high relative mortality rate (Ösby et al., Reference Ösby, Brandt, Correia, Ekbom and Sparén2001; Crump et al., Reference Crump, Sundquist, Winkleby and Sundquist2013; Hayes et al., Reference Hayes, Miles, Walters, King and Osborn2015). Pneumonia accounted for the majority of our observed respiratory-related deaths for bipolar disorder. This aligns with evidence showing that bipolar disorder patients display higher risk of pneumonia and poorer post-pneumonia outcomes than the general population (Schoepf and Heun, Reference Schoepf and Heun2014; Davydow et al., Reference Davydow, Ribe, Pedersen, Fenger-Gron, Cerimele, Vedsted and Vestergaard2016; Li et al., Reference Li, Hung, Yang, Pan, Liao, Tsai, Chen and Kuo2018). However, we noted that the proportion of respiratory-related deaths attributed to COPD was small and SMR for COPD was not elevated in bipolar disorder. This is at odds with one previous study demonstrating increased COPD-related mortality risk in bipolar disorder patients (Crump et al., Reference Crump, Sundquist, Winkleby and Sundquist2013). One possible explanation is that a significant proportion of our patients with pneumonia as their assigned death cause might have underlying COPD, a risk factor of and one of the most frequent comorbidities with pneumonia which is a major cause of death among COPD patients. As our mortality analysis was based on a single-cause of death per deceased person, this may introduce inaccuracy in death-cause assignment for patients with multi-morbidity, which is nonetheless more common in bipolar disorder than the general population. Alternatively, our finding of non-elevated mortality rate for cancers in the overall sample is in contrast to a meta-analysis reporting modestly increased cancer-related SMR for bipolar disorder (Hayes et al., Reference Hayes, Miles, Walters, King and Osborn2015). Yet, in agreement with two prior register-based studies (Ösby et al., Reference Ösby, Brandt, Correia, Ekbom and Sparén2001; Crump et al., Reference Crump, Sundquist, Winkleby and Sundquist2013), we found that women (but not men) with bipolar disorder had increased cancer-related SMR (1.25), similar to their reported estimates (ranged: 1.2–1.3). Furthermore, bipolar disorder was associated with significantly elevated SMR for colon cancer only (but not other cancer types), specifically among women. This is consistent with the only previous study which also examined SMRs for specific cancer types in bipolar disorder (Crump et al., Reference Crump, Sundquist, Winkleby and Sundquist2013). It should, however, be noted that our cancer-related SMRs might still be underestimated owing to the missing data on patients’ death cause and potential misclassification bias in cause-of-death ascertainment. Non-elevated cancer-related SMRs might also be attributable to a markedly reduced lifespan in bipolar disorder patients, rendering them more likely to die from non-cancer causes, especially among elderly-age group (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, McGrath and Laursen2019).
Patients with bipolar disorder had approximately 8 times greater risk of dying from unnatural causes than the general population. This finding is concordant with the unnatural-cause summary SMR of 7.42 reported by a recent meta-analysis (Hayes et al., Reference Hayes, Miles, Walters, King and Osborn2015). Such mortality rate was particularly elevated among the youngest-age group, with up to 16-fold increased risk for unnatural-cause deaths. Our cause-specific LYLs analyses also revealed convergent findings that the degree of contribution to the overall LYLs by unnatural-cause deaths was greatest at younger ages, and decreased with age. Notably, our lack of information regarding specific unnatural causes precludes us from clarifying mortality risk of suicide which has been consistently shown to predominate unnatural deaths for bipolar disorder. Nonetheless, our results are in line with past studies indicating that suicide risk is highest during the early course of illness (Ösby et al., Reference Ösby, Brandt, Correia, Ekbom and Sparén2001; Crump et al., Reference Crump, Sundquist, Winkleby and Sundquist2013; Medici et al., Reference Medici, Videbech, Gustafsson and Munk-Jørgensen2015), which typically emerges in late adolescence and early adulthood. This underscores the importance of early intervention for suicide prevention in bipolar disorder.
Although our analysis suggested small, albeit statistically significant, reduction in annual all-cause mortality rate over time in men with bipolar disorder, the overall mortality gap remained largely unchanged over 11 years. Both natural-cause and unnatural-cause mortality differences also showed no significant improvement across the study period in either sex or in the overall sample. Our results thus agree with most previous studies (Medici et al., Reference Medici, Videbech, Gustafsson and Munk-Jørgensen2015; Hayes et al., Reference Hayes, Marston, Walters, King and Osborn2017; Pan et al., Reference Pan, Yeh, Chan and Chang2017, Reference Pan, Yeh, Chan and Chang2020; Hansen et al., Reference Hansen, Laursen, Grøntved, Straszek, Licht and Nielsen2019; Lomholt et al., Reference Lomholt, Andersen, Sejrsgaard-Jacobsen, Ozdemir, Graff, Schjerning, Jensen, Straszek, Licht, Grontved and Nielsen2019), and suggest that bipolar disorder patients still experience pronounced physical health disparities and have not benefited equally from enhanced medical care with life expectancy improvement as compared to the general population.
Several study limitations should be noted. First, several demographic variables including educational level, socioeconomic or employment status and a number of premature mortality risk factors such as unhealthy lifestyles and obesity were not adequately recorded in medical database and thus were not included in the analysis. Second, missing data on patients’ death causes may compromise the accuracy in evaluating cause-specific SMRs. Third, as information on specific unnatural causes was not available, we were not able to investigate mortality risk for individual categories of unnatural deaths including suicide and accidents. Fourth, the study data did not contain information denoting the disorder subtypes, precluding us from examining potential differences between bipolar I and II disorders in mortality rate and life expectancy. Fifth, patients’ age of illness onset was defined by age of first-recorded diagnosis of bipolar disorder. However, previous studies showed that there could be significant delays between illness onset and ascertainment of bipolar disorder diagnosis (Dagani et al., Reference Dagani, Signorini, Nielssen, Bani, Pastore, Girolamo and Large2017), especially among patients who initially present with depressive episode (Berk et al., Reference Berk, Dodd, Callay, Berk, Fitzgerald, de Castella, Filia, Filia, Tahtalian, Biffin, Kelin, Smith, Montgomery and Kulkarni2007). Our study may thus overestimate the actual age of illness onset, and then underestimate LYLs for bipolar disorder. Sixth, a proportion of patients may have first-recorded diagnosis before the study period. This would result in survival bias and may contribute to our results of reduced relative mortality risk with age. Seventh, our lack of individual-level mortality data for the general population precludes us from estimating cause-specific excess LYLs for bipolar disorder (Plana-Ripoll et al., Reference Plana-Ripoll, Canudas-Romo, Weye, Laursen, McGrath and Andersen2020). Eighth, patients’ data were retrieved from medical-record database of public healthcare services managed by HA, and patients who were under private psychiatric care were not included in the study. However, HA is the predominant provider of psychiatric services to individuals with severe mental disorders in HK, and hence the risk of selection bias or missing treated cases of bipolar disorder was minimised. Lastly, as HK is a highly urbanised, densely populated city and is categorised by the World Bank as a high-income economy (World Bank, 2020), our findings may not be generalisable to mainland China or other Asian regions.
In conclusion, this large population-based study indicates that bipolar disorder patients have increased mortality risk and markedly reduced life expectancy in a predominantly Chinese population, with excess mortality being primarily attributable to natural causes. Our findings of persistent mortality gap across the study period highlights an urgent need for further research to identify its underlying causes and to implement multi-level, targeted interventions to promote physical health of bipolar disorder patients so as to significantly reduce their risk of preventable physical morbidity and premature mortality.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S2045796021000305
Data
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Acknowledgements
The authors would like to thank the colleagues in the Hospital Authority of Hong Kong for their kind assistance in data extraction for the current investigation.
Financial support
The study was supported by the Hong Kong Research Grants Council (grant number: 17124715). Additional financial support to undertake this research was provided by the State Key Laboratory of Brain & Cognitive Sciences, the University of Hong Kong.
Conflict of interest
The authors declared no conflicts of interest in relation to the subject of this study.
Ethical standards
The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2000.
Open access
