Next Article in Journal
UK Women’s Views of the Concepts of Personalised Breast Cancer Risk Assessment and Risk-Stratified Breast Screening: A Qualitative Interview Study
Next Article in Special Issue
The Overall Quality of Life and Oncological Outcomes Following Radical Hysterectomy in Cervical Cancer Survivors Results from a Large Long-Term Single-Institution Study
Previous Article in Journal
Dual-Wavelength Fluorescence Monitoring of Photodynamic Therapy: From Analytical Models to Clinical Studies
Previous Article in Special Issue
Socioeconomic Environment and Survival in Patients with Digestive Cancers: A French Population-Based Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Excess Mortality by Multimorbidity, Socioeconomic, and Healthcare Factors, amongst Patients Diagnosed with Diffuse Large B-Cell or Follicular Lymphoma in England

1
Inequalities in Cancer Outcomes Network, Department of Non-Communicable Disease Epidemiology, London School of Hygiene and Tropical Medicine, London WC1E 7HT, UK
2
MRC Clinical Trials Unit, Institute of Clinical Trials and Methodology, University College London, London WC1V 6LJ, UK
3
Noncommunicable Disease and Cancer Epidemiology Group, Instituto de Investigación Biosanitaria de Granada, Ibs.GRANADA, Andalusian School of Public Health, 18012 Granada, Spain
4
Centro de Investigación Biomédica en Red de Epidemiología y Salud Pública (CIBER of Epidemiology and Public Health, CIBERESP), 28029 Madrid, Spain
5
Epidemiology of Childhood and Adolescent Cancers Team, Research Centre in Epidemiology and Biostatistics (CRESS), Inserm UMR 1153, Université de Paris, 94801 Villejuif, France
6
School of Mathematics, University of Nairobi, Nairobi 30197-00100, Kenya
*
Author to whom correspondence should be addressed.
Cancers 2021, 13(22), 5805; https://doi.org/10.3390/cancers13225805
Submission received: 1 October 2021 / Revised: 10 November 2021 / Accepted: 16 November 2021 / Published: 19 November 2021

Abstract

:

Simple Summary

Diffuse large B-cell (DLBCL) and follicular lymphoma (FL) account for most non-Hodgkin lymphoma diagnoses: around 35% and 20% in England, respectively. Despite the vast contrast in survival between the subtypes, similar socioeconomic inequalities in survival have persisted over the past two decades, possibly due to the presence of comorbidities. The aim of our study was to assess the association between socioeconomic status and survival from DLBCL or FL accounting for patient and health system characteristics. We found that, for both DLBCL and FL, the most deprived patients had a higher excess mortality hazard compared to the least deprived, regardless of the comorbidity status. Our results show the need for the current framework of the National Health Service to improve the survival of DLBCL and FL patients in the most deprived areas of England, and further consideration is needed for patient-tailored management plans amongst patients with comorbidities or multimorbidities.

Abstract

(1) Background: Socioeconomic inequalities of survival in patients with lymphoma persist, which may be explained by patients’ comorbidities. We aimed to assess the association between comorbidities and the survival of patients diagnosed with diffuse large B-cell (DLBCL) or follicular lymphoma (FL) in England accounting for other socio-demographic characteristics. (2) Methods: Population-based cancer registry data were linked to Hospital Episode Statistics. We used a flexible multilevel excess hazard model to estimate excess mortality and net survival by patient’s comorbidity status, adjusted for sociodemographic, economic, and healthcare factors, and accounting for the patient’s area of residence. We used the latent normal joint modelling multiple imputation approach for missing data. (3) Results: Overall, 15,516 and 29,898 patients were diagnosed with FL and DLBCL in England between 2005 and 2013, respectively. Amongst DLBCL and FL patients, respectively, those in the most deprived areas showed 1.22 (95% confidence interval (CI): 1.18–1.27) and 1.45 (95% CI: 1.30–1.62) times higher excess mortality hazard compared to those in the least deprived areas, adjusted for comorbidity status, age at diagnosis, sex, ethnicity, and route to diagnosis. (4) Conclusions: Deprivation is consistently associated with poorer survival among patients diagnosed with DLBCL or FL, after adjusting for co/multimorbidities. Comorbidities and multimorbidities need to be considered when planning public health interventions targeting haematological malignancies in England.

1. Introduction

Non-Hodgkin lymphoma (NHL) is a heterogeneous group of malignancies, and is currently the 6th most commonly diagnosed cancer in England; in 2014, approximately 32 males and 23 females per 100,000 person years were diagnosed [1]. The heterogeneity in morphology leads to variation in survival probability; for instance, 5-year survival of follicular lymphoma (FL) (86.3%) is higher than diffuse large B-cell lymphoma (DLBCL) (54.8%) [2].
The healthcare system in England aims to offer equitable access to care for all patients. However, variability in health outcomes amongst patients with similar cancers and sociodemographic characteristics still occur; [2,3,4] convincing reasons for this variability remain a topic of interest. In 2001, the National Health Service (NHS) Cancer Plan [5] recognised, and aimed to reduce, the disparities in survival. Since implementation, there is no evidence that the Plan has had an impact on the inequalities [6,7]. The deprivation gap in survival is still apparent, despite the Plan and successive policies [5,8,9,10], illustrating the incomplete understanding of the mechanisms underlying these inequalities and raising the concern that these policies have missed the relevant targets.
Patients’ comorbidity status may impact timely diagnosis, possibly leading to treatment with more adverse effects [11]; comorbidities are, on average, more prevalent and severe amongst more deprived patients [12]. However, recent evidence indicates that comorbidity explains little of the differential cancer survival between socioeconomic groups [13,14,15]. Variations in healthcare access, such as location of residence, could partly explain the inequalities [16,17,18,19,20].
Since population-based cancer registries rarely hold reliable information on the cause of death, cancer-specific mortality estimates can be estimated with relative survival methods. These methods compare the mortality hazard (i.e., excess mortality hazard) observed in a population of cancer patients to the mortality hazard observed in the general population with identical demographic characteristics. In this context, the survival estimate derived from the excess mortality hazard is termed net survival (or cancer survival), which is interpreted as the survival where death is due directly, or indirectly, to the cancer studied, and death from other causes has been removed [21].
Overall, the association between comorbidity and cancer survival in patients with DLBCL and FL, accounting for other socio-demographic characteristics and the area of residence, remains unclear. We aim to describe the association between comorbidities and cancer survival amongst DLBCL or FL patients, while accounting for sociodemographic and economic factors, hypothesizing that the presence of comorbidities is associated with poorer survival.

2. Methods

2.1. Study Design, Participants, and Data Sources

We developed a population-based multilevel cohort study of adult patients diagnosed with DLBCL or FL between 1 January 2005 and 31 December 2013 in England. Patients were followed up until death or the end of the study on the 31 December 2015, whichever occurred first.
DLBCL and FL were defined according to the 10th revision of the International Statistical Classification of Diseases and Related Problems (ICD-10 codes C82.0–C85.9) [22]. Morphology (cell type) and topography (tumour site) were defined using renewed updates of the ICD for Oncology (ICD-O); ICD-O-3 [23] was used for diagnoses up to 2010, and ICD-O-3.1 [24] for diagnoses after 2011. Information on patients with DLBCL or FL was collected from the linkage of English cancer registry data, the Cancer Analysis System [25] (CAS), and Hospital Episode Statistics [26] (HES) data sets within the National Cancer Registry and Analysis Service (NCRAS). These datasets contained detailed information on patient and tumour characteristics (see details below).

2.2. Outcome, Exposure, and Patients’ Sociodemographic Characteristics

The outcome of the study was the time to death, or censoring, among DLBCL and FL patients 5 years after cancer diagnosis. Net survival was deduced after estimating the excess mortality hazard. Hence, we used England life tables stratified by deprivation, sex, age, and calendar year (2005–2013) to account for the overall mortality rate from the background population [27]. As follow up of patients ended in 2015 and life tables were available until 2013, we assumed that the expected mortality rates plateaued for 2014 and 2015.
Comorbidity status was the main exposure. We defined comorbidity as the existence of other chronic medical disorders, in addition to cancer, the primary disease of interest, which are causally unrelated to the primary disease [28,29]. Records from HES were used to identify patients’ comorbidity status based on a computational algorithm published elsewhere [30]. The algorithm searches for the presence of comorbidities retrospectively and defines a time window of 6 to 24 months prior to cancer diagnosis where comorbidities are recorded to avoid bias due to the presence of comorbidities related to cancer (i.e., cardiological comorbidities due to DLBCL or FL cancer treatment). Patient comorbidity status was adapted from the original Charlson comorbidity index [31] (CCI). We used the Royal College of Surgeons (RCS) modified Charlson Score (Table A1) [32]. The score removes patients with a previous malignancy to avoid bias, does not assign different weights to comorbidities, and categorises comorbidities as: no comorbidities, one comorbidity, and two or more comorbidities (multimorbidity).
Socio-demographic and economic characteristics were collected from the HES dataset. Age was specified at time of diagnosis. Sex is recorded as male or female. Ethnicity was recorded as white or other. Area-level deprivation, classified into one of five quintiles, was determined by the Index of Multiple Deprivation [33] (IMD), which was based on the Lower Super Output Area [34] (LSOA) residence of the patient at the time of cancer diagnosis. LSOA is a geographical location with a median of 1500 inhabitants. We also include the information regarding patients’ diagnosis path (route to diagnosis), a UK-specific programme, classified as: accident and emergency room diagnosis, general practitioner referral (routine and urgent referrals where the patient was not referred under two-week-wait), two-week-wait (urgent GP referral with a suspicion of cancer), and secondary care diagnosis (other outpatient and inpatient elective routes) [35].

2.3. Statistical Analysis

We tabulated the sociodemographic characteristics by DLBCL and FL. To estimate the excess mortality hazard, we used a multilevel excess hazard regression model (EHM) with a cubic B-spline with two knots placed at 1 and 3 years after diagnosis for the baseline hazard λ 0 ( t ) . We accounted for the hierarchical structure of the data via the inclusion of a random effect [36]. The statistical contribution of the random effect to the overall goodness of fit of the model was tested using a likelihood ratio test statistic with a Chi-square mixture distribution [37]. From the estimated excess hazard, we could deduce the net survival via the classical relationship between hazard and survival [38]. Net survival is the survival associated with the cancer under study, after eliminating the other causes of death.
In the EHM we included the following variables: age, sex, comorbidities (categorical, 3 categories), deprivation (categorical, 5 categories), lymphoma subtype, ethnicity, and route of cancer diagnosis. We included the non-linear effect of age using a regression spline (defined using a truncated power basis) with one knot located at 70 years of age. Furthermore, we assumed a time-dependent effect of age at diagnosis, represented by the interaction between B-spline function of time and age. The parameter estimates for the variables were interpreted conditionally on the random effect, i.e., they have a cluster-specific interpretation, where a cluster refers to a given LSOA. From the model we derived the excess mortality hazard ratios (EMHR) and their respective 95% confidence intervals (CI) for all the categorical variables, and the variance of the random effect for the LSOA. Empirical Bayes estimates of the random effect were used to explore the between-LSOA variability in the excess mortality hazard from DLBCL or FL. The random effect was tested for using a likelihood ratio test, with the reference distribution being a mixture of chi-squared distributions with 0 and 1 degrees of freedom, to account for the well-known boundary problem for random effect variances [39,40].

2.4. Missing Data Analysis

We explored the missing data mechanism for the two variables with missing data (ethnicity (FL 24.9%, DLBCL: 22.7%) and route (FL: 7.8%, DLBCL 5.0%)). Due to clustered data and partially observed categorical variables, we used the latent normal joint modelling multiple imputation approach, under a missing at random assumption (MAR) [41]. The imputation model included all fully and partially observed variables, vital status indicator, and the Nelson–Aalen estimate of the cumulative overall hazard, and accounted for clustering of patients within lower super output areas. We generated 10 imputed datasets. The multilevel EHM was fitted to each of these datasets, and results combined using Rubin’s rules [42,43]. Overall tests for the effects of age after multiple imputation were carried out using the F-based procedure for the test of multiple parameters after multiple imputation [41].
We used R software (version 4.1.2, R Development Core Team, 2020, R: A language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, Austria) for all data analyses; the mexhaz [36] package was used for excess hazard modelling and the jomo [44] package for multiple imputation.

3. Results

Overall, 15,516 (34.2%) patients were diagnosed with FL and 29,898 (65.8%) diagnosed with DLBCL in England between 2005 and 2013 (Table 1). The prevalence of at least one comorbidity was higher amongst DLBCL (10.7%) compared to FL (7.5%). The average age was lower amongst FL compared to DLBCL, 63.9 compared to 67.4 years, respectively. The prevalence of DLBCL was higher amongst deprived areas (16.0%) than FL (14.4%). ‘White’ was the most prevalent ethnicity for both FL (94.9%) and DLBCL (94.1%). GP referral was the most common route to diagnosis amongst FL (44.0%), whereas amongst DLBCL, A&E was the most common (33.8%).
In the multivariable analysis (Table 2), amongst DLBCL, and after multiple imputation, patients with comorbidity and multimorbidity showed 23% and 40% increased excess mortality compared to patients without comorbidity (i.e., EMHR: 1.23; 95% CI: 1.14–1.32, and EMHR: 1.40; CI: 1.01–1.94, respectively). Patients living in the most deprived areas had 1.22 (95% CI: 1.18–1.27) times higher excess mortality than those living in the least deprived areas. Patients diagnosed through A&E had nearly three times a higher excess mortality compared to GP referral (i.e., EMHR: 2.75; 95% CI: 2.54–2.98). Females had a significantly lower excess mortality compared to males (i.e., EMHR 0.93; 95% CI: 0.90–0.96). There was, however, no evidence of a difference in excess mortality by ethnicity (Table 2). Using a likelihood ratio test (a mixture of chi-square distributions) there was strong evidence (p < 0.001) that including the random effect improved the fit of the model.
In the multivariable analysis (Table 3), amongst FL, patients with comorbidity and multimorbidity showed 1.52 and 2.19 times the excess mortality compared to patients without comorbidity (i.e., EMHR: 1.52; 95% CI: 1.25–1.84, and EMHR: 2.19; CI: 1.45–3.31, respectively). Patients living in the most deprived areas had 1.45 (95% CI: 1.30–1.62) times higher excess mortality than those living in the least deprived areas. Patients diagnosed through A&E had nearly three times a higher excess mortality compared to GP referral (i.e., EMHR: 3.32; 95% CI: 2.49–4.43). Females had a significantly lower excess mortality compared to males (i.e., EMHR 0.89; 95% CI: 0.81–0.97). There was, however, no evidence of a difference in excess mortality by ethnicity (Table 3). Using a likelihood ratio test (a mixture of chi-square distributions), there was strong evidence (p < 0.001) that including the random effect improved the fit of the model.
Figure 1 and Figure 2 show the EMHR for patients with DLBCL and FL, respectively, according to age at diagnosis at different time since diagnosis (Figure 1A and Figure 2A), and according to time since diagnosis for different age at diagnosis (Figure 1B and Figure 2B). The excess mortality hazard for DLBCL and FL patients for different values of age at diagnosis is shown in the Appendix A (Figure A1 and Figure A2, respectively). These plots were obtained from the three-dimensional plots of EMHR, as shown in the Appendix A (Figure A3 and Figure A4, respectively). For DLBCL (Figure 1), the EMHR was higher for older patients whatever the follow-up time (Figure 1A). For those of older or younger ages, in comparison to 70-year-olds, the EMHR was markedly different immediately after, or at 5 years since, diagnosis, but was most similar around 18 months after diagnosis (Figure 1B).
For FL (Figure 2), the non-linear effect of age was almost similar whatever the time since diagnosis; being older was associated with a higher excess mortality hazard (Figure 2A). For those of older or younger ages, in comparison to 70-year-olds, the EMHR was markedly different immediately after or at 5 years since diagnosis but was most similar around 18 months after diagnosis (Figure 2B).
Figure 3 and Figure 4 show the net survival probability as predicted from the regression model amongst patients with DLBCL and FL, respectively. Amongst DLBCL patients (Figure 3), those living in more deprived areas experienced approximately 7% lower 5-year survival compared to patients in the least deprived areas (e.g., 5-year net survival, amongst those without comorbidities, was 56% for the least deprived compared to 49% for the most deprived). Amongst FL patients (Figure 4), those living in more deprived areas experienced approximately 4% lower 5-year survival compared to those living in the least deprived areas (e.g., 5-year net survival, amongst those without comorbidities, was 86% for the least deprived compared to 82% for the most deprived). For DLBCL only (Figure 3), the deprivation gap in survival was apparent from approximately 6 months after diagnosis, regardless of the comorbidity status.
In the Appendix A, we graphically illustrate the empirical Bayes (EB) estimates of the LSOA random effect for the excess mortality hazard from DLBCL and FL (Figure A5 and Figure A6, respectively). A positive EB estimate indicated a higher excess mortality hazard for a patient from that LSOA in comparison with a patient who has similar observed characteristics but from an LSOA with either a less positive, or negative EB estimate. The EB estimates were grouped by deprivation level, to which the LSOA contributed. For both DLBCL and FL (Figure A5 and Figure A6, respectively), the results show there were no outliers and approximately equal distribution of the EB estimates for each deprivation level.

4. Discussion

We found strong evidence of a higher excess mortality amongst DLBCL and FL patients diagnosed with comorbidities compared to patients without comorbidities after adjusting for age, deprivation level, ethnicity, and route to diagnosis and accounting for the patient’s area of residence; we also found a noticeable deprivation gap in cancer survival.
Differences in access to treatments, or risk of adverse effects, may explain some of the disparities in survival among DLBCL patients. Immunotherapy (rituximab) for the treatment of aggressive lymphomas (e.g., DLBCL) is known to be effective for those of an advanced age [45,46,47]. Rituximab is often used in combination with doxorubicin, an increase in dosage of which is associated with an increased incidence of adverse effects (cardiotoxicity), such as congestive heart failure [48]. Guidelines based on National Institute for Health and Care Excellence (NICE) recommend that patients at risk of cardiotoxicity, or low tolerance of intensive therapy, consider a less-intensive treatment regimen [49,50,51]. This less-intensive treatment allocation may partly explain the comorbidity inequalities in survival from DLBCL. For patients with FL, the standard management is ‘watch-and-wait’; thus, in the absence of a treatment, the comorbidity inequalities in survival may be largely explained by the presence of a comorbidity itself rather than being explained by the effect of comorbidity on treatment. Novel treatment strategies are being explored to ascertain the survival and quality of life benefits in comparison to current standards of care [52]; however, access to these treatments will depend on (i) the specialist centre available to the patient, and (ii) the patients’ influence on the decision of a treatment allocation. While the results of clinical trials provide insights into the efficacy of a treatment, they may lack external validity when the treatment is administered within the healthcare system. Pragmatic trials could be developed to understand the real-world benefit of NHL treatments, and how they are delivered, to patients of older ages or with underlying health conditions who would often not be eligible for clinical trials [53]. In addition, real-world data are becoming increasingly more available and would be of great utility to evaluate the effectiveness of treatment on the whole population [54,55].
For FL patients, we showed that the excess mortality hazard among older patients compared to younger patients is highest after just 4 years since diagnosis (Figure 2B). Since we accounted for background population mortality, and adjusted for comorbidity, the higher excess hazard could be because of histological transformation from lower to higher grades of FL. Studies suggest the risk of histological transformation increases by 3% per year since cancer onset [56]. Thus, the increased excess hazard amongst older patients may be because histological transformation, which complicates the treatment and management of FL. The effects of time-varying variables (e.g., progression/relapse, treatment response, transplant, etc.) on the excess mortality hazard would be strong predictive factors, but this information was not available from the data. This represents an interesting topic for further investigation, where we could either use an extension of the flexible parametric models for time-dependent variables, or a landmark approach [57].
The importance of understanding the association of comorbid conditions with cancer patients’ outcomes has been well documented [58]. To our knowledge, this is the first study of England cancer registry data investigating survival by comorbidity status among DLBCL and FL patients. Our results are consistent with previous findings from a Danish study, which showed that the hazard of death increased with severity of comorbidity status [59]; however, the study did not account for missing data and the association with comorbidities was potentially overestimated. The EHR associated with comorbidity decreased after accounting for missing data. The deprivation gap in survival persists even after accounting for prognostic factors such as comorbidity [59,60,61]. Smith et al. [3] reported no deprivation gap in survival; however, their study may have lacked power, and their study used the relative survival ratio, which can be biased over longer-term follow up [62].
Consistent with previous studies [4], survival after GP referral (non-emergency) diagnosis is significantly better compared with A&E. However, our study also finds that patients diagnosed through TWW, who would be expected to have worse symptoms and survival, showed no evidence of a difference in survival compared with GP referral. There are two possible reasons for the absence of a difference in the associations. Firstly, GPs could advocate for a prompt referral even though the patient is not on the TWW pathway, resulting in patients with similar access to healthcare facilities. Secondly, on the other hand, patients referred through the TWW pathway have more severe symptoms and are expected to have a higher excess hazard. Our results show no difference in the excess mortality, indicating that the TWW pathway prevents patients with more severe symptoms from having a higher excess hazard. This suggests that the performance of TWW pathway is at least as beneficial to a patient’s survival as GP referral. Other studies have suggested ways to improve outcomes for patients diagnosed with comorbidities, which include: novel treatment strategies [63], the inclusion of elderly patients in clinical trials [64,65], and the investigation of dose allocation amongst those with higher comorbidity scores [66]. However, further factors associated with the interactions between comorbidities and health care systems leading to poorer survival among DLBCL and FL cancer patients need to be studied.
The strengths of this study are that, firstly, we used a large population-based sample size obtained from cancer registry databases linked to HES, which encompasses all patients in England with a diagnosis of DLBCL and FL between 2005 and 2013. HES data encapsulate a national coverage of comorbidities diagnosed during hospital admission and may have missed comorbidities diagnosed during primary care (e.g., diabetes diagnosed during a GP consultation). However, the addition of information provided from comorbidity records captured during primary care may not improve the prediction of cancer patient survival beyond what is captured in HES data [67]. For example, information on comorbidities, such as diabetes, diagnosed outside of hospital admission are likely to have a minimal impact on the prediction of survival beyond information captured in HES. Secondly, we used the Royal College of Surgeons’ adaptation [32] of the Charlson comorbidity score, which provides a more valid measure of the patient’s comorbidity status, because it was developed within the England population healthcare data setting. Thirdly, we used a latent normal joint modelling multiple imputation to treat missing data in ethnicity and diagnostic route. This approach allows the imputation of a mix of variable types, while accounting for multilevel structures arising from the clustering of patients within LSOAs [41,68,69]. We assumed that missing data on partially observed variables were missing at random, given the observed variables; further analysis could explore the violation of this assumption and impute under a not-missing-at-random assumption.
This study has its limitations. Firstly, individual-level socioeconomic measures are recommended in addition to area-level measures [70]. Information on individual-level socioeconomic status was unavailable, but using area-level measures captures the multidimensional composition of a patient’s deprivation level, in addition to the contextual level [33,71]. Furthermore, using area-level measures, there is greater consistency in the measurement of deprivation between time periods, because deprivation scores have a high concordance amongst updates [33]. Secondly, due to data availability, we did not include tumour stage, which may have partly explained the socioeconomic inequalities in survival. However, even though reliable estimates can be obtained after the multiple imputation of partially observed variables with high proportions of missing data [72], the inclusion of tumour stage may not have provided further information for the prediction of survival beyond that of diagnostic route, because late cancer stage is strongly associated with delayed diagnostic route [73].
Survival at 1 and 5 years since diagnosis of DLBCL and FL in England trails that of other European countries [74]; however, restricting estimates to those surviving at least 1 year after diagnosis (conditional survival) shows a comparable 5-year survival [75]. This indicates that long-term survival differences are largely explained by the increased short-term mortality. Understanding long-term survival from FL is more complex due to the histological transformation of indolent lymphomas, which would require an adaptation of the treatment, support, and management from healthcare facilities. This adaptation could be compounded by the patient’s susceptibility to cardiotoxic treatments. Further studies could focus on the mechanisms and inequalities of short-term mortality, the long-term survival of patients with transformed lymphomas, and the survival of patients at risk of cardiotoxicity.

5. Conclusions

After accounting for sociodemographic factors, healthcare factors, socioeconomic deprivation, and the patient’s area of residence, comorbidities were consistently associated with poorer survival and an increased excess mortality amongst patients with DLBCL or FL in England. Furthermore, survival inequalities between socioeconomic levels in patients with DLBCL or FL persisted after accounting for the presence of comorbidities and multimorbidities. These results show the need for the current framework of the National Health Service to improve the survival of DLBCL and FL patients in the most deprived areas of England, and further consideration is needed for patient-tailored management plans amongst patients with comorbidities or multimorbidities.

Author Contributions

M.J.S., E.N.N., A.B. (Aurélien Belot) and B.R. contributed to the conception of the study and designed the study. E.N.N., A.B. (Aurélien Belot), B.R., S.G. and M.Q. provided advice on statistical methods. M.J.S. conducted the analyses of the data and prepared the draft of the manuscript, tables, and figures. E.N.N., A.B. (Aurélien Belot), B.R. and M.A.L.F. supervised the study and provided comments on the manuscript draft. E.N.N., B.R., M.Q., M.A.L.F., S.G., S.B.M., A.B. (Aurélien Belot) and A.B. (Audrey Bonaventure) provided comments on the final draft of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Cancer Research UK (grant number C7923/A18525) and the Cancer Research UK Studentship in Cancer Survival (grant number C7923/A20987). The analyses, and the writing of the manuscript were solely the responsibility of the authors. The findings and conclusions in this manuscript are those of the authors and do not necessarily represent the views of Cancer Research UK.

Institutional Review Board Statement

We obtained the statutory approvals required for this research from the Confidentiality Advisory Group (CAG) of the Health Research Authority (HRA): PIAG 1–05(c) 2007. Ethical approval was obtained from the Research Ethics Committee (REC) of the Health Research Authority (HRA): 07/MRE01/52.

Informed Consent Statement

Informed consent from participants was waived by the ethics committee. We used anonymised National Cancer Registry and Hospital Episode Statistics data. All methods were carried out in accordance with relevant guidelines and regulations.

Data Availability Statement

The data that support the findings of this study are available via application to the Public Health England Office for Data Release, but restrictions apply to the availability of these data.

Acknowledgments

We would like to thank Adrian Turculet, Data Manager of the LSHTM Inequalities in Cancer Outcomes Network, for his support and assistance with the data linkage.

Conflicts of Interest

The authors declare no potential conflicts of interest. The authors declare no support from any organisations for the submitted work. The design of the study.

Appendix A

Figure A1. Excess mortality hazard (i.e., white males, least deprived, no comorbidities, diagnosed through general practitioner referral within an average LSOA (random effect of zero)) over time since diagnosis, for different ages, amongst those diagnosed with diffuse large B-cell lymphoma (n = 29,898) in England during 2005–2013.
Figure A1. Excess mortality hazard (i.e., white males, least deprived, no comorbidities, diagnosed through general practitioner referral within an average LSOA (random effect of zero)) over time since diagnosis, for different ages, amongst those diagnosed with diffuse large B-cell lymphoma (n = 29,898) in England during 2005–2013.
Cancers 13 05805 g0a1
Figure A2. Excess mortality hazard (i.e., white males, least deprived, no comorbidities, diagnosed through general practitioner referral within an average LSOA (random effect of zero)) over time since diagnosis, for different ages, amongst those diagnosed with follicular lymphoma (n = 15,516) in England during 2005–2013.
Figure A2. Excess mortality hazard (i.e., white males, least deprived, no comorbidities, diagnosed through general practitioner referral within an average LSOA (random effect of zero)) over time since diagnosis, for different ages, amongst those diagnosed with follicular lymphoma (n = 15,516) in England during 2005–2013.
Cancers 13 05805 g0a2
Figure A3. Excess mortality hazard ratio according to age at diagnosis and time since diagnosis for patients diagnosed with diffuse large B-cell lymphoma (n = 29,898) in England between 2005 and 2013.
Figure A3. Excess mortality hazard ratio according to age at diagnosis and time since diagnosis for patients diagnosed with diffuse large B-cell lymphoma (n = 29,898) in England between 2005 and 2013.
Cancers 13 05805 g0a3
Figure A4. Excess mortality hazard ratio according to age at diagnosis and time since diagnosis for patients diagnosed with follicular lymphoma (n = 15,516) in England between 2005 and 2013.
Figure A4. Excess mortality hazard ratio according to age at diagnosis and time since diagnosis for patients diagnosed with follicular lymphoma (n = 15,516) in England between 2005 and 2013.
Cancers 13 05805 g0a4
Figure A5. Empirical Bayes estimates of the random effect of LSOA from the excess mortality hazard model for patients diagnosed with diffuse large B-cell lymphoma (n = 29,898) in England during 2005–2013.
Figure A5. Empirical Bayes estimates of the random effect of LSOA from the excess mortality hazard model for patients diagnosed with diffuse large B-cell lymphoma (n = 29,898) in England during 2005–2013.
Cancers 13 05805 g0a5
Figure A6. Empirical Bayes estimates of the random effect of LSOA from the excess mortality hazard model for patients diagnosed with follicular lymphoma (n = 15,516) in England during 2005–2013.
Figure A6. Empirical Bayes estimates of the random effect of LSOA from the excess mortality hazard model for patients diagnosed with follicular lymphoma (n = 15,516) in England during 2005–2013.
Cancers 13 05805 g0a6
Table A1. Comorbidities and their diagnostic ICD-10 codes.
Table A1. Comorbidities and their diagnostic ICD-10 codes.
ComorbidityICD-10
Myocardial infarctionI21.x, I22.x, I25.2
Congestive heart failureI11.0, I13.0, I13.2, I25.5, I42.0, I42.5–I42.9, I43.x, I50.x, P29.0
Peripheral vascular diseaseI70.x, I71.x, I73.1, I73.8, I73.9, I77.1, I79.0, I79.2, K55.1, K55.8, K55.9, Z95.8, Z95.9
Cerebrovascular diseaseG45.x, G46.x, H34.0, I60.x–I69.x
DementiaF00.x–F03.x, F05.1, G30.x, G31.1
Chronic obstructive pulmonary diseaseI27.9, J40.x–J47.x, J60.x–J67.x, J68.4, J70.1, J70.3
Rheumatic diseaseM05.x, M06.x, M31.5, M32.x–M34.x, M35.1, M35.3, M36.0
Liver diseaseB18.x, K70.0–K70.3, K70.9,
K71.3–K71.5, K71.7, K73.x, K74.x,
K76.0, K76.2–K76.4, K76.8, K76.9,
Z94.4, K71.1, K72.1, K72.9, K76.5, K76.6, K76.7, I85.0, I85.9, I86.4, I98.2, K70.4,
Diabetes without chronic complicationE10.0, E10.1, E10.6, E10.8, E10.9, E11.0, E11.1, E11.6, E11.8, E11.9, E12.0, E12.1, E12.6, E12.8, E12.9, E13.0, E13.1, E13.6, E13.8, E13.9, E14.0, E14.1, E14.6, E14.8, E14.9
Diabetes with chronic complicationE10.7, E11.2–E11.5, E11.7, E12.2–E12.5, E12.7, E13.2–E13.5, E13.7, E14.2–E14.5, E14.7
Hemiplegia or paraplegiaG04.1, G11.4, G80.1, G80.2, G81.x, G82.x, G83.0–G83.4, G83.9
Renal diseaseI12.0, I13.1, N03.2–N03.7, N05.2–N05.7, N18.x, N19.x, N25.0, Z49.0–Z49.2, Z94.0, Z99.2
AIDS/HIVB20.x–B22.x, B24.x
ICD-10: International Classification of Diseases, 10th Revision. Diabetes with/without chronic complication is combined in the RCS Charlson Comorbidity Score.

References

  1. Smittenaar, C.R.; Petersen, K.A.; Stewart, K.; Moitt, N. Cancer Incidence and Mortality Projections in the UK until 2035. Br. J. Cancer 2016, 115, 1147–1155. [Google Scholar] [CrossRef] [PubMed]
  2. Smith, A.; Crouch, S.; Lax, S.; Li, J.; Painter, D.; Howell, D.; Patmore, R.; Jack, A.; Roman, E. Lymphoma Incidence, Survival and Prevalence 2004–2014: Sub-Type Analyses from the UK’s Haematological Malignancy Research Network. Br. J. Cancer 2015, 112, 1575–1584. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Smith, A.; Crouch, S.; Howell, D.; Burton, C.; Patmore, R.; Roman, E. Impact of Age and Socioeconomic Status on Treatment and Survival from Aggressive Lymphoma: A UK Population-Based Study of Diffuse Large B-Cell Lymphoma. Cancer Epidemiol. 2015, 39, 1103–1112. [Google Scholar] [CrossRef] [Green Version]
  4. Kane, E.; Howell, D.; Smith, A.; Crouch, S.; Burton, C.; Roman, E.; Patmore, R. Emergency Admission and Survival from Aggressive Non-Hodgkin Lymphoma: A Report from the UK’s Population-Based Haematological Malignancy Research Network. Eur. J. Cancer 2017, 78, 53–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. The NHS Cancer Plan: A Plan for Investment, A Plan for Reform. Available online: https://www.thh.nhs.uk/documents/_Departments/Cancer/NHSCancerPlan.pdf (accessed on 24 March 2021).
  6. Exarchakou, A.; Rachet, B.; Belot, A.; Maringe, C.; Coleman, M.P. Impact of National Cancer Policies on Cancer Survival Trends and Socioeconomic Inequalities in England, 1996–2013: Population Based Study. BMJ 2018, 360, k764. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Maringe, C.; Li, R.; Mangtani, P.; Coleman, M.P.; Rachet, B. Cancer Survival Differences between South Asians and Non-South Asians of England in 1986–2004, Accounting for Age at Diagnosis and Deprivation. Br. J. Cancer 2015, 113, 173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Improving Outcomes: A Strategy for Cancer. Available online: https://www.gov.uk/government/publications/the-national-cancer-strategy (accessed on 24 March 2021).
  9. Improving Outcomes in Haematological Cancers: The Manual. Available online: https://www.nice.org.uk/guidance/NG47/documents/supporting-evidence (accessed on 29 March 2020).
  10. Haematological Cancers: Improving Outcomes. Available online: https://www.nice.org.uk/guidance/ng47 (accessed on 29 March 2020).
  11. Renzi, C.; Lyratzopoulos, G.; Hamilton, W.; Maringe, C.; Rachet, B. Contrasting Effects of Comorbidities on Emergency Colon Cancer Diagnosis: A Longitudinal Data-Linkage Study in England. BMC Health Serv. Res. 2019, 19, 311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Fowler, H.; Belot, A.; Ellis, L.; Maringe, C.; Luque-Fernandez, M.A.; Njagi, E.N.; Navani, N.; Sarfati, D.; Rachet, B. Comorbidity Prevalence among Cancer Patients: A Population-Based Cohort Study of Four Cancers. BMC Cancer 2020, 20, 2. [Google Scholar] [CrossRef] [PubMed]
  13. Li, R.; Daniel, R.; Rachet, B. How Much Do Tumor Stage and Treatment Explain Socioeconomic Inequalities in Breast Cancer Survival? Applying Causal Mediation Analysis to Population-Based Data. Eur. J. Epidemiol. 2016, 31, 603–611. [Google Scholar] [CrossRef] [Green Version]
  14. Fowler, H.; Belot, A.; Njagi, E.N.; Luque-Fernandez, M.A.; Maringe, C.; Quaresma, M.; Kajiwara, M.; Rachet, B. Persistent Inequalities in 90-Day Colon Cancer Mortality: An English Cohort Study. Br. J. Cancer 2017, 117, 1396. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Belot, A.; Fowler, H.; Njagi, E.N.; Luque-Fernandez, M.-A.; Maringe, C.; Magadi, W.; Exarchakou, A.; Quaresma, M.; Turculet, A.; Peake, M.D.; et al. Association between Age, Deprivation and Specific Comorbid Conditions and the Receipt of Major Surgery in Patients with Non-Small Cell Lung Cancer in England: A Population-Based Study. Thorax 2019, 74, 51–59. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Rachet, B.; Ellis, L.; Maringe, C.; Chu, T.; Nur, U.; Quaresma, M.; Shah, A.; Walters, S.; Woods, L.; Forman, D.; et al. Socioeconomic Inequalities in Cancer Survival in England after the NHS Cancer Plan. Br. J. Cancer 2010, 103, 446–453. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Comber, H.; De Camargo Cancela, M.; Haase, T.; Johnson, H.; Sharp, L.; Pratschke, J. Affluence and Private Health Insurance Influence Treatment and Survival in Non-Hodgkin’s Lymphoma. PLoS ONE 2016, 11, e0168684. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Woods, L.M.; Rachet, B.; Coleman, M.P. Origins of Socio-Economic Inequalities in Cancer Survival: A Review. Ann. Oncol. 2005, 17, 5–19. [Google Scholar] [CrossRef] [PubMed]
  19. Quaglia, A.; Vercelli, M.; Lillini, R.; Mugno, E.; Coebergh, J.W.; Quinn, M.; Martinez-Garcia, C.; Capocaccia, R.; Micheli, A. Socio-Economic Factors and Health Care System Characteristics Related to Cancer Survival in the Elderly: A Population-Based Analysis in 16 European Countries (ELDCARE Project). Crit. Rev. Oncol. Hematol. 2005, 54, 117–128. [Google Scholar] [CrossRef] [PubMed]
  20. Afshar, N.; English, D.R.; Milne, R.L. Rural–Urban Residence and Cancer Survival in High-Income Countries: A Systematic Review. Cancer 2019, 125, 2172–2184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Belot, A.; Pohar-Perme, M. Social Disparities in Cancer Survival: Methodological Considerations. In Social Environment and Cancer in Europe: Towards an Evidence-Based Public Health Policy; Launoy, G., Zadnik, V., Coleman, M.P., Eds.; Springer International Publishing: Cham, Switzerland, 2021; pp. 39–54. ISBN 978-3-030-69329-9. [Google Scholar]
  22. International Agency for Research on Cancer International Classification of Diseases for Oncology. Available online: http://codes.iarc.fr/ (accessed on 4 October 2019).
  23. Fritz, A.; Percy, C.; Jack, A.; Shanmugaratnam, K.; Sobin, L.H.; Parkin, D.M.; Whelan, S.L. International Classification of Diseases for Oncology, 3rd ed.; World Health Organisation: Geneva, Switzerland, 2000. [Google Scholar]
  24. Campo, E.; Swerdlow, S.H.; Harris, N.L.; Pileri, S.; Stein, H.; Jaffe, E.S. The 2008 WHO Classification of Lymphoid Neoplasms and beyond: Evolving Concepts and Practical Applications. Blood 2011, 117, 5019–5032. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. gov.uk. National Cancer Registry and Analysis Service. Available online: https://www.gov.uk/guidance/national-cancer-registration-and-analysis-service-ncras (accessed on 4 October 2019).
  26. NHS Digital Hospital Episode Statistics. Available online: https://digital.nhs.uk/data-and-information/data-tools-and-services/data-services/hospital-episode-statistics (accessed on 4 October 2019).
  27. Rachet, B.; Maringe, C.; Woods, L.M.; Ellis, L.; Spika, D.; Allemani, C. Multivariable Flexible Modelling for Estimating Complete, Smoothed Life Tables for Sub-National Populations. BMC Public Health 2015, 15, 1240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Porta, M. A Dictionary of Epidemiology; Oxford University Press: Oxford, UK, 2014; ISBN 9780195314496. [Google Scholar]
  29. Feinstein, A.R. The Pre-Therapeutic Classification of Co-Morbidity in Chronic Disease. J. Chronic. Dis. 1970, 23, 455–468. [Google Scholar] [CrossRef]
  30. Maringe, C.; Fowler, H.; Rachet, B.; Luque-Fernandez, M.A. Reproducibility, Reliability and Validity of Population-Based Administrative Health Data for the Assessment of Cancer Non-Related Comorbidities. PLoS ONE 2017, 12, e0172814. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Charlson, M.E.; Pompei, P.; Ales, K.L.; MacKenzie, C.R. A New Method of Classifying Prognostic Comorbidity in Longitudinal Studies: Development and Validation. J. Chronic. Dis. 1987, 40, 373–383. [Google Scholar] [CrossRef]
  32. Armitage, J.N.; van der Meulen, J.H. Identifying Co-Morbidity in Surgical Patients Using Administrative Data with the Royal College of Surgeons Charlson Score. Br. J. Surg. 2010, 97, 772–781. [Google Scholar] [CrossRef]
  33. gov.uk. Indices of Multiple Deprivation. Available online: https://www.gov.uk/government/statistics/english-indices-of-deprivation-2015 (accessed on 4 October 2019).
  34. National Health Service: Data Dictionary Lower Super Output Area. Available online: https://www.datadictionary.nhs.uk/data_dictionary (accessed on 4 October 2019).
  35. Elliss-Brookes, L.; McPhail, S.; Ives, A.; Greenslade, M.; Shelton, J.; Hiom, S.; Richards, M. Routes to Diagnosis for Cancer—Determining the Patient Journey Using Multiple Routine Data Sets. Br. J. Cancer 2012, 107, 1220–1226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Charvat, H.; Remontet, L.; Bossard, N.; Roche, L.; Dejardin, O.; Rachet, B.; Launoy, G.; Belot, A. A Multilevel Excess Hazard Model to Estimate Net Survival on Hierarchical Data Allowing for Non-Linear and Non-Proportional Effects of Covariates. Stat. Med. 2016, 35, 3066–3084. [Google Scholar] [CrossRef]
  37. Verbeke, G.; Molenberghs, G. Linear Mixed Models for Longitudinal Data, 1st ed.; Springer: New York, NY, USA, 2000. [Google Scholar]
  38. Belot, A.; Ndiaye, A.; Luque-Fernandez, M.-A.; Kipourou, D.-K.; Maringe, C.; Rubio, F.J.; Rachet, B. Summarizing and Communicating on Survival Data According to the Audience: A Tutorial on Different Measures Illustrated with Population-Based Cancer Registry Data. Clin. Epidemiol. 2019, 11, 53–65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Molenberghs, G.; Verbeke, G. Models for Discrete Longitudinal Data, 1st ed.; Springer: New York, NY, USA, 2005. [Google Scholar]
  40. Agresti, A. Categorical Data Analysis, 2nd ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2002. [Google Scholar]
  41. Carpenter, J.R.; Kenward, M.G. Multiple Imputation and Its Application, 1st ed.; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2013. [Google Scholar]
  42. Little, R.J.; Rubin, D.B. Statistical Analysis with Missing Data; John Wiley & Sons, Inc.: New York, NY, USA, 1987. [Google Scholar]
  43. Rubin, D.B. Multiple Imputation for Nonresponse in Surveys; Wiley: New York, NY, USA, 1987; ISBN 978-0-471-65574-9. [Google Scholar]
  44. Quartagno, M.; Carpenter, J.R. Jomo: A Package for Multilevel Joint Modeling Multiple Imputation. Stat. Med. 2016. [Google Scholar] [CrossRef]
  45. Coiffier, B.; Lepage, E.; Brière, J.; Herbrecht, R.; Tilly, H.; Bouabdallah, R.; Morel, P.; Van Den Neste, E.; Salles, G.; Gaulard, P.; et al. CHOP Chemotherapy plus Rituximab Compared with CHOP Alone in Elderly Patients with Diffuse Large-B-Cell Lymphoma. N. Engl. J. Med. 2002, 346, 235–242. [Google Scholar] [CrossRef] [PubMed]
  46. Coiffier, B. Rituximab in Combination with CHOP Improves Survival in Elderly Patients with Aggressive Non-Hodgkin’s Lymphoma. Semin. Oncol. 2002, 29, 18–22. [Google Scholar] [CrossRef] [PubMed]
  47. Delarue, R.; Tilly, H.; Mounier, N.; Petrella, T.; Salles, G.; Thieblemont, C.; Bologna, S.; Ghesquieres, H.; Hacini, M.; Fruchart, C.; et al. Dose-Dense Rituximab-CHOP Compared with Standard Rituximab-CHOP in Elderly Patients with Diffuse Large B-Cell Lymphoma (the LNH03-6B Study): A Randomised Phase 3 Trial. Lancet Oncol. 2013, 14, 525–533. [Google Scholar] [CrossRef]
  48. McGowan, J.V.; Chung, R.; Maulik, A.; Piotrowska, I.; Walker, J.M.; Yellon, D.M. Anthracycline Chemotherapy and Cardiotoxicity. Cardiovasc. Drugs Ther. 2017, 31, 63–75. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. National Institute for Health and Care Excellence. Non-Hodgkin’s Lymphoma: Diagnosis and Management; National Institute for Health and Care Excellence: London, UK, 2016. [Google Scholar]
  50. Tilly, H.; Gomes da Silva, M.; Vitolo, U.; Jack, A.; Meignan, M.; Lopez-Guillermo, A.; Walewski, J.; André, M.; Johnson, P.W.; Pfreundschuh, M.; et al. Diffuse Large B-Cell Lymphoma (DLBCL): ESMO Clinical Practice Guidelines for Diagnosis, Treatment and Follow-Up. Ann. Oncol. 2015, 26, v116–v125. [Google Scholar] [CrossRef] [PubMed]
  51. Bröckelmann, P.J.; McMullen, S.; Wilson, J.B.; Mueller, K.; Goring, S.; Stamatoullas, A.; Zagadailov, E.; Gautam, A.; Huebner, D.; Dalal, M.; et al. Patient and Physician Preferences for First-Line Treatment of Classical Hodgkin Lymphoma in Germany, France and the United Kingdom. Br. J. Haematol. 2019, 184, 202–214. [Google Scholar] [CrossRef] [Green Version]
  52. Solimando, A.G.; Ribatti, D.; Vacca, A.; Einsele, H. Targeting B-Cell Non Hodgkin Lymphoma: New and Old Tricks. Leuk. Res. 2016, 42, 93–104. [Google Scholar] [CrossRef] [PubMed]
  53. Merali, Z.; Wilson, J.R. Explanatory Versus Pragmatic Trials: An Essential Concept in Study Design and Interpretation. Clin. Spine Surg. 2017, 30, 404–406. [Google Scholar] [CrossRef]
  54. Ghesquières, H.; Rossi, C.; Cherblanc, F.; Le Guyader-Peyrou, S.; Bijou, F.; Sujobert, P.; Fabbro-Peray, P.; Bernier, A.; Belot, A.; Chartier, L.; et al. A French Multicentric Prospective Prognostic Cohort with Epidemiological, Clinical, Biological and Treatment Information to Improve Knowledge on Lymphoma Patients: Study Protocol of the “REal World DAta in LYmphoma and Survival in Adults” (REALYSA) Cohort. BMC Public Health 2021, 21, 432. [Google Scholar] [CrossRef]
  55. El-Galaly, T.C.; Cheah, C.Y.; Villa, D. Real World Data as a Key Element in Precision Medicine for Lymphoid Malignancies: Potentials and Pitfalls. Br. J. Haematol. 2019, 186, 409–419. [Google Scholar] [CrossRef] [PubMed]
  56. Lossos, I.S.; Gascoyne, R.D. Transformation of Follicular Lymphoma. Best Pract. Res. Clin. Haematol. 2011, 24, 147–163. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Putter, H.; van Houwelingen, H.C. Understanding Landmarking and Its Relation with Time-Dependent Cox Regression. Stat. Biosci. 2017, 9, 489–503. [Google Scholar] [CrossRef] [Green Version]
  58. Sogaard, M.; Thomsen, R.W.; Bossen, K.S.; Sorensen, H.T.; Norgaard, M. The Impact of Comorbidity on Cancer Survival: A Review. Clin. Epidemiol. 2013, 5, 3–29. [Google Scholar] [CrossRef] [Green Version]
  59. Frederiksen, B.L.; Dalton, S.O.; Osler, M.; Steding-Jessen, M.; de Nully Brown, P. Socioeconomic Position, Treatment, and Survival of Non-Hodgkin Lymphoma in Denmark--a Nationwide Study. Br. J. Cancer 2012, 106, 988–995. [Google Scholar] [CrossRef]
  60. Rachet, B.; Mitry, E.; Shah, A.; Cooper, N.; Coleman, M.P. Survival from Non-Hodgkin Lymphoma in England and Wales up to 2001. Br. J. Cancer 2008, 99, S104–S106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Bray, C.; Morrison, D.S.; McKay, P. Socio-Economic Deprivation and Survival of Non-Hodgkin Lymphoma in Scotland. Leuk. Lymphoma 2008, 49, 917–923. [Google Scholar] [CrossRef] [PubMed]
  62. Pohar Perme, M.; Estève, J.; Rachet, B. Analysing Population-Based Cancer Survival—Settling the Controversies. BMC Cancer 2016, 16, 933. [Google Scholar] [CrossRef] [Green Version]
  63. Kobayashi, Y.; Miura, K.; Hojo, A.; Hatta, Y.; Tanaka, T.; Kurita, D.; Iriyama, N.; Kobayashi, S.; Takeuchi, J. Charlson Comorbidity Index Is an Independent Prognostic Factor among Elderly Patients with Diffuse Large B-Cell Lymphoma. J. Cancer Res. Clin. Oncol. 2011, 137, 1079–1084. [Google Scholar] [CrossRef]
  64. Saygin, C.; Jia, X.; Hill, B.; Dean, R.; Pohlman, B.; Smith, M.R.; Jagadeesh, D. Impact of Comorbidities on Outcomes of Elderly Patients with Diffuse Large B-cell Lymphoma. Am. J. Hematol. 2017, 92, 989–996. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Chihara, D.; Westin, J.R.; Oki, Y.; Ahmed, M.A.; Do, B.; Fayad, L.E.; Hagemeister, F.B.; Romaguera, J.E.; Fanale, M.A.; Lee, H.J.; et al. Management Strategies and Outcomes for Very Elderly Patients with Diffuse Large B-Cell Lymphoma. Cancer 2016, 122, 3145–3151. [Google Scholar] [CrossRef] [Green Version]
  66. Janssen-Heijnen, M.L.; van Spronsen, D.J.; Lemmens, V.E.; Houterman, S.; Verheij, K.D.; Coebergh, J.W. A Population-Based Study of Severity of Comorbidity among Patients with Non-Hodgkin’s Lymphoma: Prognostic Impact Independent of International Prognostic Index. Br. J. Haematol. 2005, 129, 597–606. [Google Scholar] [CrossRef] [PubMed]
  67. Crooks, C.J.; West, J.; Card, T.R. A Comparison of the Recording of Comorbidity in Primary and Secondary Care by Using the Charlson Index to Predict Short-Term and Long-Term Survival in a Routine Linked Data Cohort. BMJ Open 2015, 5, e007974. [Google Scholar] [CrossRef] [PubMed]
  68. Carpenter, J.; Goldstein, H.; Kenward, M. REALCOM-IMPUTE Software for Multilevel Multiple Imputation with Mixed Response Types. J. Stat. Softw. 2011, 45, 1–20. [Google Scholar] [CrossRef]
  69. Quartagno, M.; Carpenter, J.R. Multiple Imputation for Discrete Data: Evaluation of the Joint Latent Normal Model. Biom. J. Biom. Z. 2019, 61, 1003–1019. [Google Scholar] [CrossRef] [PubMed]
  70. Ingleby, F.C.; Belot, A.; Atherton, I.; Baker, M.; Elliss-Brookes, L.; Woods, L.M. Assessment of the Concordance between Individual-Level and Area-Level Measures of Socio-Economic Deprivation in a Cancer Patient Cohort in England and Wales. BMJ Open 2020, 10, e041714. [Google Scholar] [CrossRef]
  71. Belot, A.; Remontet, L.; Rachet, B.; Dejardin, O.; Charvat, H.; Bara, S.; Guizard, A.-V.; Roche, L.; Launoy, G.; Bossard, N. Describing the Association between Socioeconomic Inequalities and Cancer Survival: Methodological Guidelines and Illustration with Population-Based Data. Clin. Epidemiol. 2018, 10, 561–573. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Madley-Dowd, P.; Hughes, R.; Tilling, K.; Heron, J. The Proportion of Missing Data Should Not Be Used to Guide Decisions on Multiple Imputation. J. Clin. Epidemiol. 2019, 110, 63–73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Carnerio, I.; Fry, A.; Ironmonger, L.; Connor, K.; Ormiston-Smith, N.; Henson, K.; Johnson, S.; McPhail, S.; Elliss-Brookes, L. Variation in the Routes to Cancer Diagnosis and Stage for Ten Cancer Sites. Available online: https://www.researchgate.net/publication/311602679_Variation_in_the_Routes_to_Cancer_Diagnosis_and_Stage_for_Ten_Cancer_Sites (accessed on 15 August 2021).
  74. Allemani, C.; Matsuda, T.; Di Carlo, V.; Harewood, R.; Matz, M.; Nikšić, M.; Bonaventure, A.; Valkov, M.; Johnson, C.J.; Estève, J.; et al. Global Surveillance of Trends in Cancer Survival 2000–14 (CONCORD-3): Analysis of Individual Records for 37,513,025 Patients Diagnosed with One of 18 Cancers from 322 Population-Based Registries in 71 Countries. Lancet 2018, 391, 1023–1075. [Google Scholar] [CrossRef] [Green Version]
  75. Thomson, C.S.; Forman, D. Cancer Survival in England and the Influence of Early Diagnosis: What Can We Learn from Recent EUROCARE Results? Br. J. Cancer 2009, 101, S102–S109. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Excess mortality hazard ratios according to (A) age at diagnosis at different time since diagnosis (3 months, 1 and 5 years), and (B) time since diagnosis for different age groups, amongst patients diagnosed with diffuse large B-cell lymphoma (n = 29,898) in England during 2005–2013.
Figure 1. Excess mortality hazard ratios according to (A) age at diagnosis at different time since diagnosis (3 months, 1 and 5 years), and (B) time since diagnosis for different age groups, amongst patients diagnosed with diffuse large B-cell lymphoma (n = 29,898) in England during 2005–2013.
Cancers 13 05805 g001
Figure 2. Excess mortality hazard ratios according to (A) age at diagnosis at different time since diagnosis (3 months, 1 and 5 years), and (B) time since diagnosis for different age groups, amongst patients diagnosed with follicular lymphoma (n = 15,516) in England during 2005–2013.
Figure 2. Excess mortality hazard ratios according to (A) age at diagnosis at different time since diagnosis (3 months, 1 and 5 years), and (B) time since diagnosis for different age groups, amongst patients diagnosed with follicular lymphoma (n = 15,516) in England during 2005–2013.
Cancers 13 05805 g002
Figure 3. Net survival model-based prediction for diffuse large B-cell lymphoma for each comorbidity status by deprivation level (n = 29,898) in England between 2005 and 2013. The values here are presented for 70-year-old white males diagnosed via a general practitioner referral. Values will change for other covariates, but the pattern observed here will remain.
Figure 3. Net survival model-based prediction for diffuse large B-cell lymphoma for each comorbidity status by deprivation level (n = 29,898) in England between 2005 and 2013. The values here are presented for 70-year-old white males diagnosed via a general practitioner referral. Values will change for other covariates, but the pattern observed here will remain.
Cancers 13 05805 g003
Figure 4. Net survival model-based prediction for follicular lymphoma for each comorbidity status by deprivation level (n = 15,516) in England between 2005 and 2013. The values here are presented for 70-year-old white males diagnosed via a general practitioner referral. Values will change for other covariates, but the pattern observed here will remain.
Figure 4. Net survival model-based prediction for follicular lymphoma for each comorbidity status by deprivation level (n = 15,516) in England between 2005 and 2013. The values here are presented for 70-year-old white males diagnosed via a general practitioner referral. Values will change for other covariates, but the pattern observed here will remain.
Cancers 13 05805 g004
Table 1. Distribution of cancer subtypes by patient and healthcare system characteristics for patients (n = 45,414) diagnosed with non-Hodgkin lymphoma in England during the period 2005–2013.
Table 1. Distribution of cancer subtypes by patient and healthcare system characteristics for patients (n = 45,414) diagnosed with non-Hodgkin lymphoma in England during the period 2005–2013.
Patient CharacteristicsSubtype of NHL
FLDLBCL
N = 15,516N = 29,898
Age (mean, SD)63.9 (13.6)67.4 (14.9)
Sex, n (%)
Male7318(47.2%)16,215(54.2%)
Female8198(52.8%)13,683(45.8%)
Deprivation quintiles (Q), n (%)
Least deprived (Q1)3547(22.9%)6340(21.2%)
Q23517(22.7%)6663(22.3%)
Q33294(21.2%)6246(20.9%)
Q42925(18.9%)5863(19.6%)
Most deprived (Q5)2233(14.4%)4786(16.0%)
Comorbidity status, n (%)
No comorbidity14,343(92.4%)26,718(89.4%)
One comorbidity641(4.1%)1570(5.3%)
Multimorbidity532(3.4%)1610(5.4%)
Route of diagnosis, n (%)
GP referral6297(44.0%)8157(28.7%)
A&E1869(13.1%)9617(33.8%)
Secondary care2222(15.5%)3724(13.1%)
TWW3912(27.4%)6918(24.4%)
Missing *1216(7.8%)1482(5.0%)
Ethnicity, n (%)
White11,052(94.9%)21,739(94.1%)
Others600(5.2%)1369(5.9%)
Missing *3864(24.9%)6790(22.7%)
GP: general practitioner referral, A&E: accident and emergency room, TWW: two-week-wait. Complete case analysis: missing ethnicity 23.5%; missing route to diagnosis 5.9%. * Proportions are of the total number of patients.
Table 2. Adjusted excess mortality hazard ratios for age, sex, deprivation, comorbidity, cancer subtype, route of diagnosis, ethnicity, and LSOA as random intercept for (i) complete case analysis, and (ii) after multiple imputation for patients (n = 29,898) diagnosed with diffuse large B-cell lymphoma in England during the period 2005–2013.
Table 2. Adjusted excess mortality hazard ratios for age, sex, deprivation, comorbidity, cancer subtype, route of diagnosis, ethnicity, and LSOA as random intercept for (i) complete case analysis, and (ii) after multiple imputation for patients (n = 29,898) diagnosed with diffuse large B-cell lymphoma in England during the period 2005–2013.
Patient CharacteristicsModel (i): Complete CaseModel (ii): After Imputation
HRCIp–ValueHRCIp–Value
Sex
MaleRefRef RefRef
Female0.930.89–0.980.0030.930.90–0.96<0.001
Ethnicity
WhiteRefRef RefRef
Other0.970.87–1.080.5560.990.91–1.080.809
Deprivation quintiles (Q)
Least deprived Q1RefRef RefRef
Q21.030.96–1.110.3721.000.93–1.080.922
Q31.081.00–1.160.0451.071.00–1.14 0.045
Q41.171.08–1.26<0.0011.131.04–1.230.003
Most deprived Q51.261.16–1.37<0.0011.221.18–1.27<0.001
Comorbidity status
No comorbidityRefRef RefRef
One comorbidity1.261.15–1.38<0.0011.231.14–1.32<0.001
Multimorbidity1.501.38–1.64<0.0011.401.01–1.940.043
Route of diagnosis
GP referralRefRef RefRef
A&E2.752.60–2.91<0.0012.752.54–2.98<0.001
Secondary Care1.431.22–1.67<0.0011.231.11–1.36<0.001
TWW1.331.23–1.45<0.0010.830.56–1.240.362
Random Effect
SD (SE)0.48 (0.08)--0.39 (0.04)--
GP: general practitioner referral. A&E: accident and emergency room. TWW: two-week-wait.
Table 3. Adjusted excess mortality hazard ratios for age, sex, deprivation, comorbidity, cancer subtype, route of diagnosis, ethnicity, and LSOA as random intercept for (i) complete case analysis, and (ii) after multiple imputation for patients (n = 15,516) diagnosed with follicular lymphoma in England during the period 2005–2013.
Table 3. Adjusted excess mortality hazard ratios for age, sex, deprivation, comorbidity, cancer subtype, route of diagnosis, ethnicity, and LSOA as random intercept for (i) complete case analysis, and (ii) after multiple imputation for patients (n = 15,516) diagnosed with follicular lymphoma in England during the period 2005–2013.
CharacteristicsModel (i): Complete CaseModel (ii): After Imputation
HRCIp–ValueHRCIp–Value
Sex
MaleRefRef RefRef
Female0.860.76–0.960.0100.890.81–0.970.009
Ethnicity
WhiteRefRef RefRef
Other0.590.41–0.830.0030.760.60–0.960.019
Deprivation quintiles (Q)
Least deprived Q1RefRef RefRef
Q21.090.91–1.310.3641.100.92–1.320.309
Q31.231.02–1.480.0301.110.96–1.29 0.166
Q41.371.13–1.650.0011.341.06–1.690.015
Most deprived Q51.691.38–2.06<0.0011.451.30–1.62<0.001
Comorbidity status
No comorbidityRefRef RefRef
One comorbidity1.511.19–1.91<0.0011.521.25–1.84<0.001
Multimorbidity2.381.90–3.00<0.0012.191.45–3.31<0.001
Route of diagnosis
GP referralRefRef RefRef
A&E3.182.69–3.76<0.0013.322.49–4.43<0.001
Secondary Care1.270.86–1.900.2331.220.96–1.550.107
TWW1.170.98–1.400.0841.060.63–1.780.830
Random Effect
SD (SE)0.87 (0.14)--0.69 (0.16)--
GP: general practitioner referral. A&E: accident and emergency room. TWW: two-week-wait.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Smith, M.J.; Belot, A.; Quartagno, M.; Luque Fernandez, M.A.; Bonaventure, A.; Gachau, S.; Benitez Majano, S.; Rachet, B.; Njagi, E.N. Excess Mortality by Multimorbidity, Socioeconomic, and Healthcare Factors, amongst Patients Diagnosed with Diffuse Large B-Cell or Follicular Lymphoma in England. Cancers 2021, 13, 5805. https://doi.org/10.3390/cancers13225805

AMA Style

Smith MJ, Belot A, Quartagno M, Luque Fernandez MA, Bonaventure A, Gachau S, Benitez Majano S, Rachet B, Njagi EN. Excess Mortality by Multimorbidity, Socioeconomic, and Healthcare Factors, amongst Patients Diagnosed with Diffuse Large B-Cell or Follicular Lymphoma in England. Cancers. 2021; 13(22):5805. https://doi.org/10.3390/cancers13225805

Chicago/Turabian Style

Smith, Matthew James, Aurélien Belot, Matteo Quartagno, Miguel Angel Luque Fernandez, Audrey Bonaventure, Susan Gachau, Sara Benitez Majano, Bernard Rachet, and Edmund Njeru Njagi. 2021. "Excess Mortality by Multimorbidity, Socioeconomic, and Healthcare Factors, amongst Patients Diagnosed with Diffuse Large B-Cell or Follicular Lymphoma in England" Cancers 13, no. 22: 5805. https://doi.org/10.3390/cancers13225805

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop