Next Article in Journal
Integrative Metabolomics Reveals Deep Tissue and Systemic Metabolic Remodeling in Glioblastoma
Next Article in Special Issue
Excess Mortality by Multimorbidity, Socioeconomic, and Healthcare Factors, amongst Patients Diagnosed with Diffuse Large B-Cell or Follicular Lymphoma in England
Previous Article in Journal
Detection of Pancreatic Ductal Adenocarcinoma by Ex Vivo Magnetic Levitation of Plasma Protein-Coated Nanoparticles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Socioeconomic Environment and Survival in Patients with Digestive Cancers: A French Population-Based Study

1
‘ANTICIPE’ U1086 INSERM-UCN, Normandie University UNICAEN, Centre François Baclesse, 14000 Caen, France
2
Service de Biostatistique–Bioinformatique, Pôle Santé Publique, Hospices Civils de Lyon, 69000 Lyon, France
3
Laboratoire de Biométrie et Biologie Évolutive, Équipe Biostatistique-Santé, University of Lyon 1, CNRS, UMR 5558, 69100 Villeurbanne, France
4
Digestive Cancer Registry of Burgundy, Dijon University Hospital, INSERM UMR 1231, University of Burgundy, 21079 Dijon, France
5
French Network of Cancer Registries, 31000 Toulouse, France
6
Digestive Tumors Registry of Finistère, EA SPURBO 7479, CHRU Morvan, 29200 Brest, France
7
Digestive Cancer Registry of Calvados, Caen University Hospital, ‘ANTICIPE’ U1086 INSERM-UCN, Normandie University UNICAEN, Centre François Baclesse, 14000 Caen, France
8
Research Department, Caen University Hospital, ‘ANTICIPE’ U1086 INSERM-UCN, Normandie University UNICAEN, Centre François Baclesse, 14000 Caen, France
9
Department of Digestive Surgery, University Hospital of Caen, 14000 Caen, France
10
Loire-Atlantique/Vendée Cancer Registry, 44000 Nantes, France
11
CERPOP, Université de Toulouse, Inserm, UPS, 31000 Toulouse, France
*
Author to whom correspondence should be addressed.
FRANCIM Group are listed in Appendix A.
Cancers 2021, 13(20), 5156; https://doi.org/10.3390/cancers13205156
Submission received: 22 July 2021 / Revised: 7 October 2021 / Accepted: 8 October 2021 / Published: 14 October 2021

Abstract

:

Simple Summary

Studies investigating the social gradient in digestive cancer survival are scarce, and the statistical methods used do not always consider important assumptions in survival analysis for adequate assessment. Using an ecological index (European Deprivation Index), we found a negative impact of social environment in digestive cancers net survival (especially for esophagus, stomach, bile ducts among females; colon and rectum for both sexes) and provided insight into how this social gradient in cancer survival builds up, and at what time of follow-up it appears. These results can guide clinical practice/public health actions to address social inequalities in survival by targeting digestive cancers with the greatest impact and identifying key follow-up periods to implement actions.

Abstract

Social inequalities are an important prognostic factor in cancer survival, but little is known regarding digestive cancers specifically. We aimed to provide in-depth analysis of the contextual social disparities in net survival of patients with digestive cancer in France, using population-based data and relevant modeling. Digestive cancers (n = 54,507) diagnosed between 2006–2009, collected through the French network of cancer registries, were included (end of follow-up 30 June 2013). Social environment was assessed by the European Deprivation Index. Multidimensional penalized splines were used to model excess mortality hazard. We found that net survival was significantly worse for individuals living in a more deprived environment as compared to those living in a less deprived one for esophageal, liver, pancreatic, colon and rectal cancers, and for stomach and bile duct cancers among females. Excess mortality hazard was up to 57% higher among females living in the most deprived areas (vs. least deprived) at 1 year of follow-up for bile duct cancer, and up to 21% higher among males living in the most deprived areas (vs. least deprived) regarding colon cancer. To conclude, we provide a better understanding of how the (contextual) social gradient in survival is constructed, offering new perspectives for tackling social inequalities in digestive cancer survival.

1. Introduction

The role of social environment in health has been established for many years and concerns a wide variety of diseases, both chronic and acute, including digestive pathologies [1]. Research on social determinants of health rely on measures of social situation at the individual level (through indicators such as level of education, income, employment status etc.), or on contextual indexes that encompass more distal factors from the individuals (e.g., percentage of people below the poverty line, unemployed or low-skilled in a neighbor-hood, accessibility to services, medical premises, social network, etc.), providing a more comprehensive and integrated measure of their socio-economic environment and living conditions. Currently, beyond reporting, studies on the social determination of health are warranted to identify and understand the underlying mechanisms in an attempt to guide programs and practices aimed at tackling social inequalities [2]. To that end, and including in the field of cancer, studies based on unbiased general population data covering the entire social spectrum of patients would be helpful. Concerning the survival of patients with cancer in France, in agreement with the international literature based on either individual or contextual measures of social situation [3,4,5], a previous study showed that the prognosis was worse for the most disadvantaged populations than for the least disadvantaged ones, for most cancer sites, with a marked difference for some digestive ones, such as colorectal cancer and liver cancer in males or bile duct cancer and esophageal cancer in females [6].
However, in this previous study, as in the international literature on that subject, some elements that might help understand this social gradient were not fully explored. First, the key moments in the construction of these inequalities may differ according to the site, the mode of diagnosis, the availability of screening, the kind of treatment and other prognostic factors. Second, the gradient may not be expressed in the same way at all ages. Third, the pathway of social inequalities also depends on national contextual elements (e.g., public health policies) concerning the organization of primary prevention, screening and care. To date, much of the research on this topic has not analyzed net survival; thus, it has not been possible to distinguish between mortality due to cancer and that due to other comorbidities [7,8]. Moreover, as in the previous French study [6], among the studies based on the concept of cancer net survival, most used non-parametric analyses. Consequently, they did not account for baseline hazard flexibility and the putative time-dependent and non-linear effect of variables (i.e., social environment in our case) or interaction with age, which could be a limitation in cancer survival analysis [9]. It is possible that inequalities are built throughout the follow-up and add up through the different steps of cancer management (therapeutic choices, medical follow-up, treatment compliance, management of side-effects or relapses, etc.), leading to an increase in the social gradient of cancer survival over time. Conversely, it is possible that specific factors linked to the beginning of the cancer management induce social inequalities in cancer survival, which are no longer present thereafter, leading to a reduction in the social gradient of survival over time [10,11]. Furthermore, age-related factors could increase or reduce social inequalities in survival [12]. For instance, specific and close monitoring of patients in oncogeriatric departments could reduce the social gradient in this population. Conversely, the isolation or lack of autonomy of the elderly may make it worse.
The objective of this study was to provide in-depth analysis of the social disparities in survival at the contextual level in patients with digestive cancer in France for each cancer site, through flexible excess mortality hazard models using multidimensional penalized splines [13,14] and drawing on cancer registry population-based data.

2. Materials and Methods

2.1. Population and Data

The study population, which comprised 32,837 males and 21,670 females with diagnosed digestive cancer, was derived from the population-based data of three specialized digestive and 13 general French cancer registries belonging to the French Network of Cancer Registries (FRANCIM). All digestive cancer cases diagnosed and registered between 1 January 2006, and 31 December 2009 in patients over 15 years old were included, except for the Gironde and Lille area cancer registries for which cases were available only for 2008 and 2009, and for the Haute–Vienne cancer registry for which cases were available only for 2009. Cases were followed-up until the date of death or 30 June 2013 (except for loss to follow-up, which accounted for about 2% of all registered cases/cancers combined [6]). The study was approved by the Consultative Committee for the Processing of Health Research Data (CCTIRS) and the French Data Protection Authority (CNIL, authorization n° 913013).
Cancers were classified according to the topographical and morphological codes of the 3rd edition of the International Classification of Diseases for Oncology (ICD-O-3). Distribution of sites according to sex is shown in Table 1. Survival time corresponded to the difference between the date of death (or last information on vital status) and the date of cancer diagnosis. In about 0.1% of the cases for which the date of death was the same as the date of cancer diagnosis (survival time = 0), a survival time of 0.5 day was assigned. Other available data were age at diagnosis and social environment, assessed by the French version of the European Deprivation Index (EDI). For several years, there has been a wide consensus that socioeconomic status cannot be summed up by a single indicator. At individual level, socioeconomic status is usually explored in three fields: income, education and/or socio-professional category. The collection of individual socio-economic data consistently comes up against the problem of their absence from medical files or from medico-administrative databases (as in the French cancer registries), together with the issue of legal protection. Moreover, individual data may be exposed to a non-response bias in questionnaire surveys. In addition, assessment at individual level cannot account for contextual elements related to factors such as the place of residence (green spaces, criminality, equipment for physical activity, supply of consumer goods). Aggregated composite indexes have been constructed to allow the measurement of the social environment in large unbiased samples. Using a weighted combination of census data, they make it possible to integrate contextual elements and to assess the socio-economic environment at different area levels. For these reasons, we chose to use the EDI in this study, which is a European trans-cultural aggregate measure of relative poverty constructed with data from the European Union Statistics on Income and Living Conditions survey (EU-SILC) of EUROSTAT and declined in several European countries using census data [15]. It is currently available for seven European countries. The French version is available at the “IRIS” (Ilots Regroupés pour l’Information Statistique) level, which are the smallest geographic areas for which census data are available. There are more than 15.000 IRIS in France. The EDI was a continuous variable ranging from −16 to 56 (median value: −0.9) at the national level; the higher the index, the greater the deprivation in the IRIS. All the French IRIS are distributed into national quintiles of deprivation according to the EDI value nationwide, quintile 1 being the least deprived and quintile 5 the most deprived.
For each cancer case in the study population, the patient’s address at the time of diagnosis was geolocalized using Geographic Information Systems (ArcGIS 10.2, ESRI, Redlands, CA, USA) in order to be allocated to an IRIS with its corresponding EDI value and national quintile.

2.2. Statistical Analyses

All analyses were computed separately for each cancer site.
To model cancer-specific mortality in the absence of available data on the cause of death in the FRANCIM registries, analyses were conducted with the excess mortality framework [16]. Thus, at given values of time (t), age at diagnosis (a) and EDI, the observed mortality hazard h of an individual is as follows:
h t , a , EDI , z   = h E t , a , EDI   +   h P a + t , z
where h E is the excess mortality hazard (EMH), i.e., the mortality directly or indirectly due to cancer, and hp is the expected mortality (hp is the all-cause mortality hazard of the general French population at age a + t, given the demographic characteristics z of that individual). Here, z is composed of the variables sex, year of death and the residence Département (which is the main territorial and administrative division in France). The expected mortality hp was provided by French life tables, produced by the National Institute of Statistics and Economic Studies (Institut National de la Statistique et des Etudes Economiques, INSEE).
The EMH was modeled using multidimensional penalized splines, which allows to model flexible baseline hazard, non-linear and non-proportional (i.e., time-dependent) effects of covariates as well as interactions [13,14]. This novel statistical model offers flexibility by using regression splines while limiting overfitting issues thanks to penalization. Four models based on penalized splines were adjusted and the best one was selected according to the corrected Akaike information criterion (AIC) [17]:
M0: l o g h E t , a   = tensor t , a
M1: l o g h E t , a , E D I   = tensor t , a   +   s EDI
M1b: l o g h E t , a , E D I   = tensor t , a   +   s EDI   + tint t , EDI
M2: l o g h E t , a , E D I   = tensor t , a , EDI
The keywords tensor, s, and tint respectively stand for a penalized tensor product spline, a one-dimensional penalized spline, and a penalized tensor product spline only containing interaction terms. Restricted cubic splines were used as one-dimensional splines or as marginal splines in a tensor product spline. We used 6, 5, and 5 knots for time, age, and EDI, respectively. The locations of these knots correspond to the percentiles of the distribution of time, age and EDI among deceased patients (this choice being justified by previous work [14]).
Smoothing parameters were estimated by optimizing the laplace approximate marginal likelihood (LAML) criterion and regression parameters by maximizing the penalized likelihood of the survival model. If M0 was selected, this meant that the effect of EDI on the EMH was considered as non-significant. If M1 was selected, the effect of EDI on the EMH was considered as significant and steady over time since diagnosis and identical, regardless of age at diagnosis. If M1b was selected, the effect of EDI was considered as significant and time-dependent but not age-dependent. If M2 was selected, the effect of EDI was considered as significant and age-dependent (or time- and age-dependent). The potential non-linearity of the effect of EDI (included as a continuous variable) was considered in all four models. The adequacy of the selected model was checked by comparing the net survival curves predicted by the model and those derived from a non-parametric method (Pohar-Perme) [7], using R software (R Core Team, Vienna, Austria, version 3.5.1) and the ‘relsurv’ (2.2.3) package.
Net survival probabilities and the EMH predicted by the selected model were then computed and plotted as a function of time since diagnosis, according to five key values for deprivation, defined as the median value of EDI in each quintile of the national distribution: mQ1 (least deprived, EDI = −4.2), mQ2 (EDI = −2.4), mQ3 (EDI = −0.9), mQ4 (EDI = 0.8), mQ5 (most deprived, EDI = 5.1). To represent the social gradient of cancer survival, the excess hazard ratio (EHR) of mQ5, mQ4, mQ3 and mQ2 versus mQ1 was computed. This was performed for several times of follow-up if the effect of EDI was found to be time-dependent, i.e., if M1b or M2 was selected.
Net survival methods assume that the death rate in the patient population is higher than the all-causes death rate in the background population. This is a reasonable assumption for cancers (especially digestive cancers), which is why such methods are relevant and commonly used in cancer studies. In addition, if this assumption would have been false, we would have encountered model convergence problems [7], which was not the case.
Since missing data for EDI accounted for less than 1%, we performed complete case analyses.
French life tables provided by INSEE are not stratified on deprivation, although background mortality in the general population might substantially differ according to socio-economic position; thus, social gradient in net survival for patients with cancer may be due at least partly to socially determined comorbidities. Therefore, as in previous studies [5,6,7,8,9,10,11,12,13,14,15,16,17,18], we conducted sensitivity analyses using two sets of simulated deprivation-specific French life tables. The simulations were based on the following: a) the mortality rate ratios by quintiles of the income domain score of the Index of Multiple Deprivation [19] provided by the deprivation-specific England life tables [20], England having large mortality inequalities as in France [21]; and b) the mortality rate ratios by quintiles of net income per consumption unit (individual level) provided by The Permanent Demographic Sample (Echantillon Démographique Permanent, EDP), a large-scale socio-demographic panel established in France [22]. Thus, in both scenarios, we applied the social gradient in mortality observed in the corresponding general population to the original French life tables. Since the external sources used for the simulations provided extreme social gradients in background mortality, our sensitivity analyses were conducted under “extreme correction” of the potential bias.
All the models were fitted using R software (3.5.1) with the “survPen” package (1.0.1) [23].

3. Results

Table 1 shows descriptive statistics by sex and cancer site as well as distribution of the study population into the national quintiles of deprivation and population net survival 1 month, 1 year and 5 years after cancer diagnosis provided by the best model selected by the AIC (see methods). Median age ranged between 66–77 years old across the cancer sites. As expected, 5-year cancer net survival probabilities were low for pancreas (males: 8.07%; females: 6.69%), liver (males: 14.61%; females: 14.22%), esophagus (males: 14.65%; females: 15.41%), bile ducts (males: 19.18%; females: 15.44%) and stomach (males: 23.7%; females: 27.69%) and higher for small intestines (males: 54.07%; females: 51.34%), rectum (males: 59.69%; females: 60.34%) and colon (males: 60.48%; females: 59.9%). Distribution of patients into the five national quintiles of EDI was around 20% for males, and it was a bit more heterogeneous among females, with less than 15% of patients in Q1 (least deprived) for esophagus or stomach, and 27.4% of patients in Q5 (most deprived) for liver cancer (resulting probably from a social gradient of incidence for these cancers).
As described in the Section 2, different models of the EMH were tested for each site and sex to assess whether net survival was influenced by EDI, and if so (M1, M1b or M2 model selected), whether this influence varied over time since diagnosis (M1b) and according to age at diagnosis (M2). As summarized in Table 2, net survival varied significantly according to EDI for all cancer sites but not for small intestine in both sexes (M0), nor for stomach and bile ducts in males (M0). It was dependent on time since diagnosis (M1b) of pancreas in males and for stomach, colon and bile ducts in females. This effect was not dependent on age at diagnosis for any site (no M2 selected).
Figure 1 shows the prediction of net survival by the selected model for each cancer site in the first five years after diagnosis for males (Figure 1a) and females (Figure 1b) according to medians of EDI national quintiles, when the selected model included an effect of EDI on net survival. Since the EDI effect was never dependent on age, we chose to represent net survival at 70 years of age at diagnosis.
Figure 1 shows how net survival regularly decreased with increasing deprivation score for each cancer site, reflecting the social gradient of survival. The gradient was particularly marked for colon, rectum, esophagus as well as bile ducts for females, with wider gaps between the five curves representing the five levels of EDI as compared to the other digestive cancers (pancreas, liver, as well as stomach among females).
Table 3 provides the estimates of excess mortality hazard ratios (EHR) between median values of national deprivation quintiles using mQ1 (least deprived) as reference for all cancer sites for which the EDI effect was significant. Detailed estimates are given at 1 month, 1 year and 5 years after cancer diagnosis for sites for which this effect was time-dependent.
Table 2 and Table 3 show that the effect of EDI was moderate for most digestive cancers studied. For most sites for which a significant effect of EDI on excess mortality was highlighted, the effect of EDI did not depend on time since diagnosis (4/5 sites for males: esophagus, colon, rectum and liver; and 4/7 for females: esophagus, rectum, liver and pancreas). For these sites, the prognosis progressively worsened with deprivation (Figure 1 and Table 3), with a regular and gradual worsening of the prognosis between the different quintiles of deprivation from the least deprived (first quintile) to the most deprived (fifth quintile). The pejorative effect of deprivation on net survival was particularly marked with esophageal (EHRmQ5 vs. mQ1: 1.44, 95% confidence interval (CI): 1.13–1.83), stomach (significant at 1 month of follow-up and not thereafter: EHRmQ5 vs. mQ1 (1 month): 1.48, 95% CI: 1.09–2.03) and bile ducts (significant at 1 year of follow-up only: EHRmQ5 vs. mQ1 (1 year): 1.57, 95% CI: 1.21–2.02) cancers among females, and with colon and rectal cancers in both males (colon, EHRmQ5 vs. mQ1: 1.21, 95% CI: 1.08–1.35; rectum, EHRmQ5 vs. mQ1: 1.2, 95% CI: 1.07–1.34) and females (colon, EHRmQ5 vs. mQ1 (1 month): 1.2, 95% CI: 1.01–1.43 and EHRmQ5 vs. mQ1 (1 year): 1.23, 95% CI: 1.06–1.43; rectum, EHRmQ5 vs. mQ1: 1.23, 95% CI: 1.09–1.39).
Table 3 and Figure 2 show the pattern of EDI effect changing over time for stomach, colon, and bile duct cancers in females and for pancreatic cancer in males. For stomach and colon among females and pancreas among males, the effect of deprivation on excess mortality was maximal during the first months or years after diagnosis, and then diminished over time, or even tended to reverse in the years furthest from diagnosis, although never significantly. However, for these three sites, the effect of EDI at the beginning of the follow-up was borderline significant, with the lower bound of the 95% CI close to 1.00. The effect of EDI was specific for bile duct cancer among females, with an increasing social gradient according to time since diagnosis. The effect was substantial as the EHR reached a value of 2.1 [1.1–3.9] from 3.8 years after diagnosis.
To assess the impact of potential bias associated with the lack of French life tables stratified on deprivation, two sets of simulated French deprivation-specific life tables were used. In these sensitivity analyses (Table 4), similar results were found overall, except for colon, rectal and esophageal cancers in males and colon cancer in females, for which no significant effect of EDI on net survival was highlighted by the selected model in the sensitivity analyses (Model M0 selected) contrary to main analyses (Model M1 or M1b selected). However, trends remained for these sites when M1 was tested in the sensitivity analyses, as was also the case in non-parametric analyses using the simulated life tables and Pohar-Perme method. In pancreas for males and bile ducts for females, the effect of EDI on net survival moved from time-dependent in the main analyses (Model M1b) to steady over time in the sensitivity analyses (Model M1). In all the other cases, results were steady from main analyses to sensitivity analyses (e.g., for liver cancer among females: EHRmQ5 vs. mQ1 (main analyses): 1.18, 95% CI: 1.02–1.35; EHRmQ5 vs. mQ1 (sensitivity analyses using England life tables): 1.17, 95% CI: 1.01–1.34; EHRmQ5 vs. mQ1 (sensitivity analyses using EDP): 1.16, 95% CI: 1.01–1.34).

4. Discussion

Our results strongly suggest that there is a pejorative influence of deprivation (at the contextual level) on the prognosis of digestive cancers: esophageal, liver, pancreatic, colon and rectal cancers for both sexes, and stomach and bile duct cancers among females. The extent of the influence of social environment on net survival was greater for females than for males. In females, the maximum reached for bile duct cancer with an excess mortality hazard increased by 57% at one year of follow-up in the most deprived quintile compared to the least deprived quintile; in males, the maximum reached for colon cancer with an excess mortality hazard increased by 21% in the most deprived quintile compared to the least deprived quintile.
The use of net survival and flexible modeling of excess mortality due to cancer allowed us to show that the influence of deprivation on the excess mortality was similar in all age groups, that it could be time-dependent for some cancers, and that there was a progressive gradient across the social scale for all digestive cancer sites. The models showed that the social gradient of survival was observable from the first months or years after diagnosis for almost all digestive cancer sites, and that it remained throughout the patient’s care for most of them.
Social environment had a stronger effect on cancer survival in females. Except for esophageal and liver cancer, it is unlikely that this difference was due to differences in the biological or histological nature of the cancers. Similarly, as social environment was assessed in an aggregated manner using a geographical approach, it is unlikely that it was assessed differently for males and females. Therefore, these differences between males and females are likely due to the way in which cancers are diagnosed, managed and treated, as well as to a putative social determinism of participation in screening that is stronger in females than in males, particularly for colon cancer where these differences were marked. Unfortunately, due to the lack of data on the stage of extension at diagnosis or screening practice in our dataset, this hypothesis could not be tested.
Colon and rectal cancers are the cancers in which the impact of social environment on survival has been most studied, particularly in England. Our finding of an excess mortality risk greater than 20% for most deprived people as compared to least is consistent with published studies reporting social disparities in survival at the expense of the most deprived, whether it be colon cancer [4,24,25], rectal cancer [26,27] or colorectal cancer [18,28,29,30,31]. For colon cancer in females, our results suggest that social inequalities accumulate almost exclusively in the first months after diagnosis. This confirms data obtained with different models in England, Ireland and Spain, some of which explained social inequalities in survival mainly by the stage of extension at the time of diagnosis of the disease and treatment [24,27,30,32,33]. Similar results have been reported for rectal cancer with a high frequency of patients presenting in an emergency setting [27] and for both colon and rectal localizations combined [30]. However, other studies suggested that this gradient may develop at a distance from diagnosis, as suggested by the meta-analysis of Malietzis [34], which pointed out the relationship between social status and adjuvant chemotherapy modalities, and the study of Lyratzopoulos [26], which clearly showed that, before release, therapeutic innovations aggravate social inequalities in survival. Unfortunately, we could not investigate such a relationship because those data were unavailable.
Concerning liver cancer, our results show a significant effect of EDI on survival but with a smaller impact than for other digestive localizations, especially in males with an excess mortality risk of around 10% for the most deprived as compared to the least deprived. A pejorative and significant effect of social deprivation has been found in other studies conducted in the United States (SEER Program) [35,36], Korea [37], Australia and Canada [3,38], often with a stronger effect than in our study. The availability of powerful data from the CONCORD 3 study allowed Shao [39] to demonstrate a correlation between the Human Development Index, a composite measure of health, education, and economy, and liver cancer survival.
Similar studies are rarer for other digestive cancer sites. For pancreatic cancer as for other localizations, our results show that the influence of deprivation was maximal in the months following diagnosis, and that it was more marked in males than in females. These results are consistent with previous work that associated deprivation with less surgical resection [40,41,42,43,44].
We evidenced the social determination of survival even for esophageal cancer, which is characterized by a strong social determination of incidence [45] and short overall survival. In esophageal cancer, we found not only a stronger effect but also a clearer social gradient in females than in males, where there seemed to be a plateau effect for the most deprived groups. Our findings are consistent with previous research from Korea [37], the United States [46] and China [47].
Concerning stomach cancer, we did not find a significant effect in males, although the number of cases was equivalent to that of most of the other cancer sites. In females, the results are difficult to interpret with a strong but never significant reversal of the direction of the association between deprivation and excess mortality over the years after diagnosis. Our results are not consistent with those found in Korea in males over 60 years of age but are consistent with those found in England in 2008 [48,49]. Until now, we have been unable to explain the results concerning gastric cancer, in which there was a pattern of risk reversal over time for females. This issue therefore requires further investigation.
Apart from gastric cancer, the localizations for which there was no significant effect of social environment on survival were those with the lowest number of cases in our study. Concerning bile duct cancer in males, there was an excess of mortality that increased with deprivation but never significantly. In females, however, bile ducts were the localization where social deprivation had the greatest impact on survival, with a significant excess mortality risk of more than 50% one year after diagnosis for the most deprived females. A recent US study [42] using SEER data 2007–2015 reported the influence of social status on the surgical resection rate and survival, whereby the more privileged received wider resection.
We did not find any effect of social environment on survival for the small intestine in either sex. These results are consistent with those of the only other study we know of in small bowel cancer that did not evidence any effect of social deprivation on survival, even though the study included more than 5000 patients [50]. Similar results have been found regarding mortality inequalities among females [51].
In our study, a social gradient of survival was found for almost all digestive cancer sites, regardless of their prognosis, the availability of screening, the conditions of their diagnosis or their therapeutic management. The putative mechanisms underlying this phenomenon are numerous. Our results also suggest that there is a social gradient across the social spectrum, with a regular continuum from the most advantaged to the most disadvantaged, each social class having fewer chances of survival than the one immediately above. Similarly, we recently demonstrated how such a social gradient of survival is strong enough to create a social gradient of mortality, including for cancers such as colorectal cancer with the lowest incidence in the most deprived [52]. These findings rely on contextual/environmental social situation only since information at the individual level was not available in our data. Considering both levels and using multilevel analysis would have been more accurate and should be considered for future studies. Nevertheless, aggregated environmental indexes of deprivation have been recognized to be good proxies of the social situation at the individual level [53]. In addition, previous studies have shown that social environment itself may play a role in health related outcomes, especially cancer survival and incidence [54,55]. Our results therefore confirm these previous findings and underline the interest of also investigating the social context in which individuals live, in order to better understand the social determinants of cancer survival.
Our original statistical modeling methods revealing interactions over time showed that the social gradient of survival was not formed exclusively at a distance from diagnosis in any type of digestive cancer. For most sites, the absence of variation in excess mortality over time suggests that the construction of social inequalities occurs throughout the medical course of the disease, thus highlighting the role of the organization of care. However, for several sites, these inequalities are most likely to develop during the first few months following diagnosis. This phenomenon was particularly marked for colorectal cancer, thus highlighting the importance of access to screening in the development of social inequalities in survival [24,30].
Our study has several strengths. First, most studies that have examined this topic classically analyze crude survival with the Cox model. Studies similar to ours that model net survival [3,18,30,35,56] are free of gender- and age-related co-morbidities and can thus model excess mortality directly due to disease. Second, compared to the non-parametric evaluations of net survival, our flexible method allowed an in-depth population-based analysis and may have contributed to uncovering potential underlying mechanisms such as non-proportional and time-dependent effects.
The study also has limitations. First, the analysis was limited by the lack of data on cancer extension and modalities of treatment, which are the most important cancer prognostic factors, often related to social situation themselves. Unfortunately, such parameters are not routinely collected by the French cancer registries (which conversely present the advantage of providing exhaustive and high quality data with large coverage of the French population). A perspective to continue and complete this work would be to conduct a “high resolution” study with collection of various clinical and biological parameters, based on a smaller sample. Nevertheless, we think that our study provides a first highlight of the problem of social inequalities in digestive cancers survival in France and paves the way for future research. Second, in the absence of a mortality table of the general population as a function of the level of social deprivation, models such as ours do not allow socially determined causes of death to be considered and cannot be used to unequivocally attribute differences to a social differentiation in disease management, no matter how sophisticated they are. However, the sensitivity analyses, albeit conducted under a high hypothesis of error and revealing the presence of the expected overestimation bias, did not question the overall findings of the study. EHR reported in this study are probably overestimated due to that limitation, but sensitivity analyses suggest that the error might not be important enough to substantially modify the value of the EHR and to contradict our findings. Although these were only sensitivity analyses, to us, this method, which has already been used and approved in previous published studies [5,8,57], was the best way to confirm our findings. Nonetheless, we are aware that this bias will be correctly accounted for only when French life tables stratified by deprivation for the general population become available. Third, we may have introduced uncertainty because of multiple testing, and by using a two steps procedure to estimate excess hazard rates (i.e., choice of model using AIC and then prediction of excess hazard/net survival from this model) which is a common problem in studies on cancer net survival based on flexible model. However, use of multidimensional penalized splines considerably reduced the number of steps in the model building strategy by reducing the number of candidate models, thus limiting the number of tests and the uncertainty as compared to previous strategies. In the end, the extent of model uncertainty due to such an AIC strategy in a penalized setting, based on only four candidate models, is probably slight.
Considering these limitations, estimates provided in this study must be interpreted with caution, keeping in mind that the social gradient might be slightly overestimated.
From a methodological point of view, when required variables are available, the most relevant studies to analyze the underlying mechanisms of construction of social inequalities in health are those based on a mediation model, since they aim to quantify the relative contribution of the different paths of construction of these inequalities. The mediation analysis conducted by Frederiksen [29] is particularly interesting. After accounting for the potential effects of the stage of extension at diagnosis and the mode of treatment, it raises the possibility of a direct effect (i.e., not mediated by differences in cancer diagnosis or management) of the social environment on an individual’s ability to survive cancer. The highlighting of such a direct effect can provide support for the hypothesis of allostatic load [58], based on the exhaustion of stress defense mechanisms. The assumption is that the lower the financial, social or cultural capital, the greater the individual’s load. When repeatedly placed in conditions of difficulty to meet essential needs, this permanent state of demand (“social stress”) overstretches an individual’s ability to adapt, particularly their ability to cope in such a way that the necessary balance in life is ensured. When external demand exceeds an individual’s capacity to adapt, the mobilization of mechanisms to maintain the balance can become deleterious from a health perspective. Such a hypothesis has been confirmed in the field of cardiovascular and neurological diseases but has received less attention in the field of cancerology. Beyond the repeated observation of social inequalities in survival for patients with cancer in different national contexts, this hypothesis is one of the avenues that requires further exploration. In addition, the dataset supporting our findings is ancient and it would be interesting to reproduce this study with more recent data (when available), which will also enable us to compare periods and investigate evolutions in the social gradient in digestive cancers survival in France.

5. Conclusions

These research perspectives would be interesting to explore for digestive cancers, especially since our results confirm the existence of a systematic and unidirectional social gradient in survival for the majority of digestive cancer sites in France, and since it has been shown that the social gradient in survival contributes more strongly than the social gradient in incidence to the overall social gradient in cancer-related mortality [52].

Author Contributions

Conceptualization, G.L. and O.D.; methodology, L.T., M.F., L.R., G.L. and O.D.; software, L.T., M.F. and L.R.; validation, L.T., M.F., A.-M.B., M.R., V.B., M.C., V.J., O.D., L.R., A.A., FRANCIM Group, F.M. and G.L.; formal analysis, L.T.; investigation, L.T., M.F., A.-M.B., M.R., V.B., M.C., V.J., O.D., L.R., A.A., FRANCIM Group, F.M. and G.L.; resources, G.L. and O.D.; data curation, A.-M.B., M.R., M.C., V.B., V.J., O.D., FRANCIM Group, F.M. and G.L.; writing—original draft preparation, G.L. and L.T.; writing—review and editing, L.T., M.F., A.-M.B., M.R., V.B., M.C., V.J., O.D., L.R., A.A., FRANCIM Group, F.M. and G.L.; visualization, L.T., G.L., O.D. and V.B.; supervision, G.L. and O.D.; project administration, G.L.; funding acquisition, G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the ‘Institut de Recherche en Santé Publique’ (IReSP) (no grant number), the Institut National du Cancer (INCa) (no grant number), the ‘Ligue Nationale Contre le Cancer’ (no grant number) and Santé Publique France (SpF) (no grant number).

Institutional Review Board Statement

The study was approved by the Consultative Committee for the Processing of Health Research Data (CCTIRS) and the French Data Protection Authority (CNIL, authorization n° 913013). Study was performed in accordance with the Declaration of Helsinki.

Informed Consent Statement

Not applicable since patients are not identifiable (anonymous data).

Data Availability Statement

Data and code are available upon reasonable request.

Acknowledgments

The authors thank the staff of each member registry of the FRANCIM network who participated in the collection of data and MapInMed platform. The authors would also like to thank the working group on sensitivity analyses: Aurélien Belot, Bernard Rachet, Camille Maringe, Libby Ellis, Laurent Roche, Nadine Bossard, Zoé Uhry, Valérie Jooste, Pascale Grosclaude, Véronique Bouvier, Sébastien Lamy, Michael Mounie, Anne-Sophie Woronoff, Ludivine Launay, Ophélie Merville.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A

Members of the French Network of Cancer Registries (FRANCIM): Françoise GALATEAU-SALLE (Registre multicentrique du mésothéliome à vocation nationale (Mesonat)), Anne-Marie BOUVIER (Registre Bourguignon des cancers digestifs), Simona BARA (Registre des cancers de la Manche), Clarisse JOACHIM-CONTARET (Registre des cancers de la Martinique), Bénédicte LAPOTRE-LEDOUX (Registre général des cancers de la Somme), Claire SCHVARTZ (Registre des cancers de la thyroïde Marne Ardennes), Sandrine PLOUVIER (Registre des cancers de Lille et de sa region), Guy LAUNOY (Registre des cancers digestifs du Calvados), Karima HAMMAS (Registre des cancers du Haut-Rhin), Sandrine DABAKUYO (Registre des cancers du sein et des cancers gynécologiques de Côte d’Or), Pascale GROSCLAUDE (Registre des cancers généraux du Tarn), Xavier TROUSSARD (Registre des hémopathies malignes de Basse-Normandie), Marc MAYNADIE (Registre des hémopathies malignes de Côte d’Or), Alain MONNEREAU (Registre des hémopathies malignes de la Gironde), Brigitte TRETARRE (Registre général des tumeurs de l’Hérault), Florence MOLINIE and Anne COWPPLI-BONY (Registre des tumeurs de Loire-Atlantique/Vendée), Anne-Sophie WORONOFF (Registre des tumeurs du Doubs et du Territoire de Belfort), Isabelle BALDI (Registre des tumeurs primitives du système nerveux en Gironde), Jean-Baptiste NOUSBBAUM (Registre Finistérien des tumeurs digestives), Gaëlle COUREAU (Registre général des cancers de la Gironde), Jacqueline DELOUMEAUX (Registre général des cancers de la Guadeloupe), Marc COLONNA (Registre général des cancers de l’Isère), Michel VELTEN (Registre général des cancers du Bas-Rhin), Tania D’ALMEIDA (Registre général des cancers en région Limousin), Anne-Valérie GUIZARD (Registre général des tumeurs du Calvados), Jacqueline CLAVEL (Registre national des hémopathies malignes de l’enfant (RNHME)), Brigitte LACOUR (Registre national des tumeurs solides de l’enfant (RNTSE)), Françoise BORSON-CHAZOT (Registre Rhône Alpin des cancers thyroïdiens), Gautier DEFOSSEZ (Registre des cancers de Poitou-Charentes), Sylvie LAUMOND (Cancers généraux-Nouvelle Calédonie), Emmanuel CHIRPAZ (Registre des Cancers de la Réunion), Laure-Manuella DESROZIERS-IMOUNGA (Registre des cancers généraux de Guyane); [email protected]

References

  1. Marmot, M.; Allen, J.; Bell, R.; Bloomer, E.; Goldblatt, P.; Consortium for the European Review of Social Determinants of Health and the Health Divide. WHO European review of social determinants of health and the health divide. Lancet 2012, 380, 1011–1029. [Google Scholar] [CrossRef]
  2. Whitehead, M. A typology of actions to tackle social inequalities in health. J. Epidemiol. Commun. Health 2007, 61, 473–478. [Google Scholar] [CrossRef]
  3. Stanbury, J.F.; Baade, P.D.; Yu, Y.; Yu, X.Q. Cancer survival in New South Wales, Australia: Socioeconomic disparities remain despite overall improvements. BMC Cancer 2016, 16, 1–9. [Google Scholar] [CrossRef] [Green Version]
  4. Kwak, M.; Kim, C. Disparities by Age, Sex, Tumor Stage, Diagnosis Path, and Area-level Socioeconomic Status in Survival Time for Major Cancers: Results from the Busan Cancer Registry. J. Korean Med. Sci. 2017, 32, 1974–1983. [Google Scholar] [CrossRef] [PubMed]
  5. Nakaya, T.; Miyashiro, I.; Ito, Y.; Nakayama, T.; Ioka, A.; Tsukuma, H.; Rachet, B. Socioeconomic inequalities in cancer survival: A population-based study of adult patients diagnosed in Osaka, Japan, during the period 1993–2004. Acta Oncol. 2014, 53, 1423–1433. [Google Scholar] [CrossRef]
  6. Tron, L.; Belot, A.; Fauvernier, M.; Remontet, L.; Bossard, N.; Launay, L.; Bryere, J.; Monnereau, A.; Dejardin, O.; Launoy, G.; et al. Socioeconomic environment and disparities in cancer survival for 19 solid tumor sites: An analysis of the French Network of Cancer Registries (FRANCIM) data. Int. J. Cancer 2019, 144, 1262–1274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Perme, M.P.; Stare, J.; Estève, J. On Estimation in Relative Survival. Biometrics 2011, 68, 113–120. [Google Scholar] [CrossRef]
  8. Schaffar, R.; Rachet, B.; Belot, A.; Woods, L.M. Estimation of net survival for cancer patients: Relative survival setting more robust to some assumption violations than cause-specific setting, a sensitivity analysis on empirical data. Eur. J. Cancer 2017, 72, 78–83. [Google Scholar] [CrossRef] [Green Version]
  9. Danieli, C.; Remontet, L.; Bossard, N.; Roche, L.; Belot, A. Estimating net survival: The importance of allowing for informative censoring. Stat. Med. 2012, 31, 775–786. [Google Scholar] [CrossRef]
  10. 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, 561–573. [Google Scholar] [CrossRef] [Green Version]
  11. Rachet, B.; Woods, L.M.; Mitry, E.; Riga, M.; Cooper, N.; Quinn, M.J.; Steward, J.; Brenner, H.; Esteve, J.; Sullivan, R.; et al. Cancer survival in England and Wales at the end of the 20th century. Br. J. Cancer 2008, 99, S2–S10. [Google Scholar] [CrossRef] [PubMed]
  12. Shafique, K.; Morrison, D. Socio-Economic Inequalities in Survival of Patients with Prostate Cancer: Role of Age and Gleason Grade at Diagnosis. PLoS ONE 2013, 8, e56184. [Google Scholar] [CrossRef] [PubMed]
  13. Fauvernier, M.; Roche, L.; Uhry, Z.; Tron, L.; Bossard, N.; Remontet, L.; The Challenges in the Estimation of Net Survival Working Survival Group. Multi-dimensional penalized hazard model with continuous covariates: Applications for studying trends and social inequalities in cancer survival. J. R. Stat. Soc. Ser. C. 2019, 68, 1233–1257. [Google Scholar] [CrossRef]
  14. Remontet, L.; Uhry, Z.; Bossard, N.; Iwaz, J.; Belot, A.; Danieli, C.; Charvat, H.; Roche, L.; The CENSUR Working Survival Group. Flexible and structured survival model for a simultaneous estimation of non-linear and non-proportional effects and complex interactions between continuous variables: Performance of this multidimensional penalized spline approach in net survival trend analysis. Stat. Methods Med. Res. 2018, 28, 2368–2384. [Google Scholar] [CrossRef]
  15. Guillaume, E.; Pornet, C.; Dejardin, O.; Launay, L.; Lillini, R.; Vercelli, M.; Marí-Dell’Olmo, M.; Fontelo, A.F.; Borrell, C.; Ribeiro, A.I.; et al. Development of a cross-cultural deprivation index in five European countries. J. Epidemiol. Commun. Health 2016, 70, 493–499. [Google Scholar] [CrossRef] [Green Version]
  16. Estève, J.; Benhamou, E.; Croasdale, M.; Raymond, L. Relative survival and the estimation of net survival: Elements for further discussion. Stat. Med. 1990, 9, 529–538. [Google Scholar] [CrossRef]
  17. Wood, S. Generalized Additive Models: An Introduction with R, 2nd ed.; Chapman & Hall/CRC: London, UK, 2017. [Google Scholar]
  18. Antunes, L.; Mendonça, D.; Bento, M.J.; Rachet, B. No inequalities in survival from colorectal cancer by education and socioeconomic deprivation—A population-based study in the North Region of Portugal, 2000–2002. BMC Cancer 2016, 16, 1–12. [Google Scholar] [CrossRef] [Green Version]
  19. Jordan, H.; Roderick, P.; Martin, D. The Index of Multiple Deprivation 2000 and accessibility effects on health. J. Epidemiol. Commun. Health 2004, 58, 250–257. [Google Scholar] [CrossRef] [PubMed]
  20. Di Carlo, V.R.B.; Bannon, F.; Woods, L.M.; Maringe, C.; Bonaventure, A.; Coleman, M.P.; Allemani, C. Life Tables for the CONCORD Programme. Available online: http://csg.lshtm.ac.uk/life-tables (accessed on 8 January 2020).
  21. Mackenbach, J.P.; Bopp, M.; Deboosere, P.; Kovacs, K.; Leinsalu, M.; Martikainen, P.; Menvielle, G.; Regidor, E.; de Gelder, R. Determinants of the magnitude of socioeconomic inequalities in mortality: A study of 17 European countries. Health Place 2017, 47, 44–53. [Google Scholar] [CrossRef]
  22. Blanpain, N. L’espérance de Vie par Niveau de Vie—Méthode et Principaux Résultats; Documents de Travail n° F1801; Insee: Paris, France, 2018; pp. 1–47. [Google Scholar]
  23. Fauvernier, M.; Remontet, L.; Uhry, Z.; Bossard, N.; Roche, L. survPen: An R package for hazard and excess hazard modelling with multidimensional penalized splines. J. Open Source Softw. 2019, 4, 1434. [Google Scholar] [CrossRef]
  24. 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–1404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Donnelly, C.; Hart, N.; McCrorie, A.D.; Donnelly, M.; Anderson, L.; Ranaghan, L.; Gavin, A. Predictors of an early death in patients diagnosed with colon cancer: A retrospective case–control study in the UK. BMJ Open 2019, 9, e026057. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Lyratzopoulos, G.; Barbiere, J.M.; Rachet, B.; Baum, M.; Thompson, M.R.; Coleman, M. Changes over time in socioeconomic inequalities in breast and rectal cancer survival in England and Wales during a 32-year period (1973–2004): The potential role of health care. Ann. Oncol. 2011, 22, 1661–1666. [Google Scholar] [CrossRef] [PubMed]
  27. Comber, H.; Sharp, L.; Cancela, M.D.C.; Haase, T.; Johnson, H.; Pratschke, J. Causes and outcomes of emergency presentation of rectal cancer. Int. J. Cancer 2016, 139, 1031–1039. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Jang, B.; Chang, J.H. Socioeconomic status and survival outcomes in elderly cancer patients: A national health insurance service-elderly sample cohort study. Cancer Med. 2019, 8, 3604–3613. [Google Scholar] [CrossRef] [Green Version]
  29. Frederiksen, B.L.; Osler, M.; Harling, H.; Ladelund, S.; Jørgensen, T. Do patient characteristics, disease, or treatment explain social inequality in survival from colorectal cancer? Soc. Sci. Med. 2009, 69, 1107–1115. [Google Scholar] [CrossRef]
  30. Møller, H.; Sandin, F.; Robinson, D.; Bray, F.; Klint, Å.; Linklater, K.M.; Lambert, P.; Påhlman, L.; Holmberg, L.; Morris, E. Colorectal cancer survival in socioeconomic groups in England: Variation is mainly in the short term after diagnosis. Eur. J. Cancer 2012, 48, 46–53. [Google Scholar] [CrossRef]
  31. Godbole, C.; Bhangu, A.; Bowley, D.M.; Subramanian, T.; Kamarajah, S.K.; Karandikar, S. Impact of socioeconomic deprivation on short-term outcomes and long-term overall survival after colorectal resection for cancer. Int. J. Color. Dis. 2019, 34, 2101–2109. [Google Scholar] [CrossRef]
  32. Askari, A.; Nachiappan, S.; Currie, A.; Latchford, A.; Stebbing, J.; Bottle, A.; Athanasiou, T.; Faiz, O. The relationship between ethnicity, social deprivation and late presentation of colorectal cancer. Cancer Epidemiol. 2017, 47, 88–93. [Google Scholar] [CrossRef]
  33. Nuñez, O.; Rodríguez Barranco, M.; Fernández-Navarro, P.; Redondo Sanchez, D.; Luque Fernández, M.; Pollán Santamaría, M.; Sánchez, M.J. Deprivation gap in colorectal cancer survival attributable to stage at diagnosis: A population-based study in Spain. Cancer Epidemiol. 2020, 68, 101794. [Google Scholar] [CrossRef]
  34. Malietzis, G.; Mughal, A.; Currie, A.C.; Anyamene, N.; Kennedy, R.H.; Athanasiou, T.; Jenkins, J.T. Factors Implicated for Delay of Adjuvant Chemotherapy in Colorectal Cancer: A Meta-analysis of Observational Studies. Ann. Surg. Oncol. 2015, 22, 3793–3802. [Google Scholar] [CrossRef] [PubMed]
  35. Wang, S.; Sun, H.; Xie, Z.; Li, J.; Hong, G.; Li, D.; Mallampati, S.; Zhou, X.; Zhou, C.; Zhang, H.; et al. Improved survival of patients with hepatocellular carcinoma and disparities by age, race, and socioeconomic status by decade, 1983–2012. Oncotarget 2016, 7, 59820–59833. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Abdel-Rahman, O. Treatment choices and outcomes of non-metastatic hepatocellular carcinoma patients in relationship to neighborhood socioeconomic status: A population-based study. Int. J. Clin. Oncol. 2020, 25, 861–866. [Google Scholar] [CrossRef]
  37. Lee, H.-E.; Zaitsu, M.; Kim, E.-A.; Kawachi, I. Occupational Class and Cancer Survival in Korean Men: Follow-Up Study of Nation-Wide Working Population. Int. J. Environ. Res. Public Health 2020, 17, 303. [Google Scholar] [CrossRef] [Green Version]
  38. Jembere, N.; Campitelli, M.A.; Sherman, M.; Feld, J.J.; Lou, W.; Peacock, S.; Yoshida, E.; Krahn, M.D.; Earle, C.; Thein, H.-H. Influence of Socioeconomic Status on Survival of Hepatocellular Carcinoma in the Ontario Population; A Population-Based Study, 1990–2009. PLoS ONE 2012, 7. [Google Scholar] [CrossRef] [Green Version]
  39. Shao, S.-Y.; Hu, Q.; Wang, M.; Zhao, X.-Y.; Wu, W.-T.; Huang, J.-M.; Liang, T.-B. Impact of national Human Development Index on liver cancer outcomes: Transition from 2008 to 2018. World J. Gastroenterol. 2019, 25, 4749–4763. [Google Scholar] [CrossRef]
  40. Thobie, A.; Mulliri, A.; Dolet, N.; Eid, Y.; Bouvier, V.; Launoy, G.; Alves, A.; Dejardin, O. Socioeconomic status impacts survival and access to resection in pancreatic adenocarcinoma: A high-resolution population-based cancer registry study. Surg. Oncol. 2018, 27, 759–766. [Google Scholar] [CrossRef] [PubMed]
  41. Engberg, H.; Steding-Jessen, M.; Øster, I.; Jensen, J.W.; Fristrup, C.W.; Møller, H. Regional and socio-economic variation in survival after a pancreatic cancer diagnosis in Denmark. Dan. Med. J. 2020, 67, A08190438. [Google Scholar]
  42. Swords, D.S.; Mulvihill, S.J.; Brooke, B.S.; Firpo, M.A.; Scaife, C.L. Size and Importance of Socioeconomic Status-Based Disparities in Use of Surgery in Nonadvanced Stage Gastrointestinal Cancers. Ann. Surg. Oncol. 2019, 27, 333–341. [Google Scholar] [CrossRef]
  43. Azap, R.A.; Paredes, A.Z.; Diaz, A.; Hyer, J.M.; Pawlik, T.M. The association of neighborhood social vulnerability with surgical textbook outcomes among patients undergoing hepatopancreatic surgery. Surgery 2020, 168, 868–875. [Google Scholar] [CrossRef]
  44. Van Roest, M.H.; van der Aa, M.A.; van der Geest, L.G.; de Jong, K.P. The Impact of Socioeconomic Status, Surgical Resection and Type of Hospital on Survival in Patients with Pancreatic Cancer. A Population-Based Study in The Netherlands. PLoS ONE 2016, 11, e0166449. [Google Scholar] [CrossRef]
  45. Bryere, J.; Dejardin, O.; Launay, L.; Colonna, M.; Grosclaude, P.; Launoy, G. Socioeconomic status and site-specific cancer incidence, a Bayesian approach in a French Cancer Registries Network study. Eur. J. Cancer Prev. 2018, 27, 391–398. [Google Scholar] [CrossRef]
  46. Erhunmwunsee, L.; Gulack, B.C.; Rushing, C.; Niedzwiecki, D.; Berry, M.F.; Hartwig, M.G. Socioeconomic Status, Not Race, Is Associated With Reduced Survival in Esophagectomy Patients. Ann. Thorac. Surg. 2017, 104, 234–244. [Google Scholar] [CrossRef] [Green Version]
  47. Kou, K.; Baade, P.D.; Guo, X.; Gatton, M.; Cramb, S.; Lu, Z.; Fu, Z.; Chu, J.; Xu, A.; Sun, J. Area socioeconomic status is independently associated with esophageal cancer mortality in Shandong, China. Sci. Rep. 2019, 9, 6388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Leigh, Y.; Seagroatt, V.; Goldacre, M.; McCulloch, P. Impact of socio-economic deprivation on death rates after surgery for upper gastrointestinal tract cancer. Br. J. Cancer 2006, 95, 940–943. [Google Scholar] [CrossRef] [PubMed]
  49. Gossage, J.A.; Forshaw, M.J.; Khan, A.A.; Mak, V.; Moller, H.; Mason, R.C. The effect of economic deprivation on oesophageal and gastric cancer in a UK cancer network. Int. J. Clin. Pr. 2009, 63, 859–864. [Google Scholar] [CrossRef]
  50. Shack, L.G.; Wood, H.E.; Kang, J.Y.; Brewster, D.; Quinn, M.J.; Maxwell, J.D.; Majeed, A. Small intestinal cancer in England & Wales and Scotland: Time trends in incidence, mortality and survival. Aliment. Pharmacol. Ther. 2006, 23, 1297–1306. [Google Scholar] [CrossRef]
  51. Sundquist, J.; Li, X.; Sundquist, K. Neighborhood deprivation and mortality in individuals with cancer: A multilevel analysis from Sweden. Eur. J. Cancer. Prev. 2012, 21, 387–394. [Google Scholar] [CrossRef]
  52. Bryere, J.; Tron, L.; Menvielle, G.; Launoy, G.; (Francim), T.F.N.O.C.R.; Galateau-Salle, F. The respective parts of incidence and lethality in socioeconomic differences in cancer mortality. An analysis of the French network Cancer registries (FRANCIM) data. Int. J. Equity Health 2019, 18, 1–11. [Google Scholar] [CrossRef] [PubMed]
  53. Krieger, N. Overcoming the absence of socioeconomic data in medical records: Validation and application of a census-based methodology. Am. J. Public Health 1992, 82, 703–710. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Hagedoorn, P.; Vandenheede, H.; Vanthomme, K.; Gadeyne, S. Socioeconomic position, population density and site-specific cancer mortality: A multilevel analysis of Belgian adults, 2001–2011. Int. J. Cancer 2017, 142, 23–35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Bryere, J.; Menvielle, G.; Dejardin, O.; Launay, L.; Molinie, F.; Stucker, I.; Luce, D.; Launoy, G. Neighborhood deprivation and risk of head and neck cancer: A multilevel analysis from France. Oral Oncol. 2017, 71, 144–149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Newnham, A.; Quinn, M.J.; Babb, P.; Kang, J.Y.; Majeed, A. Trends in oesophageal and gastric cancer incidence, mortality and survival in England and Wales 1971–1998/1999. Aliment. Pharmaco.l Ther. 2003, 17, 655–664. [Google Scholar] [CrossRef] [PubMed]
  57. Ellis, L.; Coleman, M.; Rachet, B. The impact of life tables adjusted for smoking on the socio-economic difference in net survival for laryngeal and lung cancer. Br. J. Cancer 2014, 111, 195–202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Seeman, M.; Stein Merkin, S.; Karlamangla, A.; Koretz, B.; Seeman, T. Social status and biological dysregulation: The “status syndrome” and allostatic load. Soc. Sci. Med. 2014, 118, 143–151. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Net survival over time since diagnosis by level of deprivation for each cancer site for which EDI effect was significant, predicted by selected model for 70-year-old (a) males and (b) females. EDI: European Deprivation Index; mQi: median value of national deprivation quintile i (see Section 2).
Figure 1. Net survival over time since diagnosis by level of deprivation for each cancer site for which EDI effect was significant, predicted by selected model for 70-year-old (a) males and (b) females. EDI: European Deprivation Index; mQi: median value of national deprivation quintile i (see Section 2).
Cancers 13 05156 g001aCancers 13 05156 g001b
Figure 2. Excess hazard mortality ratio (EHR) between mQ5 (most deprived) and mQ1 (least deprived), over time for cancer sites for which effect of EDI on net survival was time-dependent (M1b selected). EHR: excess mortality hazard ratios; mQi: median value of national deprivation quintile i (see Section 2).
Figure 2. Excess hazard mortality ratio (EHR) between mQ5 (most deprived) and mQ1 (least deprived), over time for cancer sites for which effect of EDI on net survival was time-dependent (M1b selected). EHR: excess mortality hazard ratios; mQi: median value of national deprivation quintile i (see Section 2).
Cancers 13 05156 g002
Table 1. Sociodemographic characteristics, EDI national quintile distribution and populational net survival (computed from the selected model, at 5 years of follow-up) according to cancer site.
Table 1. Sociodemographic characteristics, EDI national quintile distribution and populational net survival (computed from the selected model, at 5 years of follow-up) according to cancer site.
Cancer SiteTopography Codenn DeathsMedian AgeQ1 EDI (%)Q2 EDI (%)Q3 EDI (%)Q4 EDI (%)Q5 EDI (%)5-Year Net Survival [95%CI]
Males
Esophagus C15325028316617.320.220.819.921.814.65 [11.98;17.66]
Stomach C16349327777218.120.619.819.422.123.70 [20.89;26.66]
Small intestine C175122616819.721.520.320.717.854.07 [46.62;60.94]
ColonC1810,11948157218.621.320.920.019.260.48 [57.97;62.9]
RectumC19-C21622029177018.721.422.718.818.459.69 [56.69;62.57]
LiverC22497943086919.220.719.420.220.514.61 [12.52;16.91]
Bile ductsC23–248487107320.017.219.721.321.719.18 [15.01;23.80]
PancreasC25341631556919.320.121.218.720.78.07 [6.06;10.5]
Females
EsophagusC157055997214.819.919.424.421.615.41 [10.37;21.61]
Stomach C16190514077714.919.720.221.723.627.69 [23.09;32.62]
Small intestine cancerC173992017018.520.119.019.023.351.34 [43.78;58.56]
ColonC18866940377517.220.121.020.121.659.9 [57.24;62.43]
RectumC19-C21467920647217.621.019.920.421.160.34 [57.05;63.50]
LiverC2211159567415.118.419.619.527.414.22 [10.73;18.3]
Bile ductsC23–2410118737715.320.318.423.722.315.44 [11.55;20.01]
PancreasC25318729627515.919.220.621.422.86.69 [5.01;8.72]
EDI: European Deprivation Index; Qi EDI (%): proportion of individuals in population study belonging to national deprivation quintile i; 95%CI: 95% confidence interval.
Table 2. Effect of deprivation assessed by EDI on net survival according to cancer site and sex, as assessed by selected flexible model.
Table 2. Effect of deprivation assessed by EDI on net survival according to cancer site and sex, as assessed by selected flexible model.
Cancer SiteSignificant Effect of EDIEffect of EDI Time-Dependent †Effect of EDI Age-Dependent †Model Selected ‡
Males
EsophagusYESNONOM1
StomachNO------M0
Small IntestineNO------M0
ColonYESNONOM1
RectumYESNONOM1
LiverYESNONOM1
Bile ductsNO------M0
PancreasYESYESNOM1b
Females
EsophagusYESNONOM1
StomachYESYESNOM1b
Small IntestineNO------M0
ColonYESYESNOM1b
RectumYESNONOM1
LiverYESNONOM1
Bile ductsYESYESNOM1b
PancreasYESNONOM1
EDI: European Deprivation Index; †: not applicable (---) if EDI effect was not significant; ‡: effect of EDI on excess mortality hazard: M0: not significant, M1: significant, steady over time since diagnosis and identical regardless of age at diagnosis, M1b: significant, time-dependent but not age-dependent.
Table 3. Excess mortality hazard ratio (EHR) estimates for median value of national deprivation quintiles according to sex and cancer site.
Table 3. Excess mortality hazard ratio (EHR) estimates for median value of national deprivation quintiles according to sex and cancer site.
Cancer SiteTime of Follow-UpEHREHR [95% CI]EHR [95% CI]EHR [95% CI]EHR [95% CI]
Cancer Site †Time of Follow-Up ‡mQ1 (ref)mQ2mQ3mQ4mQ5
Males
EsophagusNA11.06 [1.01–1.12]1.11 [1.02–1.21]1.14 [1.03–1.25]1.14 [1.01–1.29]
Colon NA11.06 [1.01–1.11]1.11 [1.03–1.19]1.15 [1.06–1.25]1.21 [1.08–1.35]
RectumNA11.04 [1.01–1.06]1.07 [1.03–1.11]1.1 [1.04–1.17]1.2 [1.07–1.34]
LiverNA11.02 [1–1.05]1.04 [1–1.09]1.06 [1.01–1.13]1.11 [1.02–1.21]
Pancreas1 month11.07 [1.02–1.13]1.11 [1.03–1.21]1.14 [1.03–1.25]1.18 [1.03–1.34]
1 year11.05 [1–1.1]1.07 [0.99–1.16]1.07 [0.98–1.17]1.06 [0.94–1.19]
5 years10.96 [0.87–1.06]0.91 [0.76–1.09]0.84 [0.64–1.09]0.67 [0.41–1.08]
Females
EsophagusNA11.09 [1.01–1.18]1.17 [1.03–1.33]1.25 [1.06–1.48]1.44 [1.13–1.83]
Stomach 1 month11.18 [1.02–1.36]1.33 [1.06–1.66]1.43 [1.1–1.86]1.48 [1.09–2.03]
1 year10.94 [0.83–1.06]0.8 [0.65–0.99]0.73 [0.56–0.95]0.77 [0.57–1.05]
5 years10.89 [0.63–1.25]0.79 [0.44–1.41]0.68 [0.3–1.54]0.49 [0.12–1.97]
Colon 1 month11.09 [1.01–1.17]1.15 [1.03–1.29]1.18 [1.03–1.35]1.2 [1.01–1.43]
1 year11.06 [1–1.11]1.11 [1.01–1.22]1.15 [1.03–1.3]1.23 [1.06–1.43]
5-year10.9 [0.81–1]0.83 [0.68–1]0.75 [0.56–1]0.61 [0.36–1.04]
RectumNA11.04 [1.02–1.07]1.08 [1.03–1.12]1.12 [1.05–1.19]1.23 [1.09–1.39]
Liver NA11.03 [1–1.06]1.06 [1.01–1.11]1.09 [1.01–1.18]1.18 [1.02–1.35]
Bile ducts1 month11.05 [0.95–1.16]1.1 [0.94–1.28]1.16 [0.95–1.4]1.25 [0.95–1.63]
1 year11.09 [0.99–1.21]1.19 [1.02–1.38]1.31 [1.08–1.58]1.57 [1.21–2.02]
5 years11.15 [0.93–1.43]1.31 [0.89–1.93]1.51 [0.86–2.68]2.06 [0.73–5.78]
PancreasNA11.02 [1–1.03]1.03 [1–1.06]1.04 [0.99–1.1]1.08 [0.99–1.19]
CI: confidence interval; EHR: excess mortality hazard ratios; mQi: median value of national deprivation quintile i (see Section 2); NA: not applicable. † Except sites for which no significant effect of EDI on net survival was highlighted (M0 model selected); ‡ NA if no time-dependent effect of EDI was highlighted (M1 model selected).
Table 4. Selected model and excess hazard mortality ratio (EHR) between mQ5 (most deprived) and mQ1 (least deprived) according to analysis (main or sensitivity) for each cancer site and sex.
Table 4. Selected model and excess hazard mortality ratio (EHR) between mQ5 (most deprived) and mQ1 (least deprived) according to analysis (main or sensitivity) for each cancer site and sex.
Cancer SiteMain AnalysesSA (1)SA (2)
Selected Model †Time of Follow UpEHR mQ5/mQ1 [95% CI]Selected Model †Time of Follow UpEHR mQ5/mQ1 [95% CI]Selected Model †Time of Follow UpEHR mQ5/mQ1 [95% CI]
Males
StomachM0NANAM0NANAM0NANA
LiverM1NA1.11 [1.02;1.21]M1NA1.07 [0.99;1.16]M1NA1.07 [0.99;1.15]
EsophagusM1NA1.14 [1.01;1.29]M0NA1.11 [0.99;1.26] ‡M0NA1.1 [0.98;1.25] ‡
PancreasM1b1 month1.18 [1.03;1.34]M1NA1.07 [0.98;1.17]M1NA1.07 [0.98;1.16]
M1b1 year1.06 [0.94;1.19]
M1b5 years0.67 [0.41;1.08]
ColonM1NA1.21 [1.08;1.35]M0NA1.05 [0.95;1.17] ‡M0NA1.03 [0.93;1.13] ‡
RectumM1NA1.2 [1.07;1.34]M0NA1.1 [0.98;1.24] ‡M0NA1.07 [0.94;1.21] ‡
Bile ductsM0NANAM0NANAM0NANA
Females
StomachM1b1 month1.48 [1.09;2.03]M1b1 month1.48 [1.08;2.02]M1b1 month1.47 [1.08;2.01]
M1b1 year0.77 [0.57;1.05]M1b1 year0.76 [0.55;1.03]M1b1 year0.75 [0.55;1.02]
M1b5 years0.49 [0.12;1.97]M1b5 years0.44 [0.1;1.86]M1b5 years0.43 [0.1;1.85]
LiverM1NA1.18 [1.02;1.35]M1NA1.17 [1.01;1.34]M1NA1.16 [1.01;1.34]
EsophagusM1NA1.44 [1.13;1.83]M1NA1.41 [1.12;1.79]M1NA1.41 [1.11;1.78]
PancreasM1NA1.08 [0.99;1.19]M1NA1.08 [0.98;1.18]M1b1 year1.03 [0.92;1.15]
M1b5 years0.72 [0.39;1.3]
ColonM1b1 year1.23 [1.06;1.43]M0NA1.06 [0.97;1.17] ‡M0NA1.05 [0.95;1.15] ‡
M1b5 years0.61 [0.36;1.04]
RectumM1NA1.23 [1.09;1.39]M1NA1.18 [1.05;1.34]M1NA1.17 [1.03;1.32]
Bile ductsM1b1 year1.57 [1.21;2.02]M1NA1.42 [1.14;1.76]M1NA1.41 [1.13;1.75]
M1b5 years2.06 [0.73;5.78]
Males&Females
Small intestineM0NANAM0NANAM0NANA
CI: confidence interval; EHR: excess mortality hazard ratios; mQi: median value of national deprivation quintile i (see Section 2); NA: not applicable; SA: sensitivity analyses: (1) Simulated deprivation-specific French life tables derived from deprivation-specific England life tables; (2) Simulated deprivation-specific French life tables derived from mortality rate ratios by net income per consumption unit provided by The Permanent Demographic Sample (Echantillon Démographique Permanent, EDP). † Effect of EDI on excess mortality hazard; M0: not significant; M1: significant, steady over time since diagnosis and identical regardless of age at diagnosis; M1b: significant, time-dependent but not age-dependent; ‡ EHR estimated from model M1.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tron, L.; Fauvernier, M.; Bouvier, A.-M.; Robaszkiewicz, M.; Bouvier, V.; Cariou, M.; Jooste, V.; Dejardin, O.; Remontet, L.; Alves, A.; et al. Socioeconomic Environment and Survival in Patients with Digestive Cancers: A French Population-Based Study. Cancers 2021, 13, 5156. https://doi.org/10.3390/cancers13205156

AMA Style

Tron L, Fauvernier M, Bouvier A-M, Robaszkiewicz M, Bouvier V, Cariou M, Jooste V, Dejardin O, Remontet L, Alves A, et al. Socioeconomic Environment and Survival in Patients with Digestive Cancers: A French Population-Based Study. Cancers. 2021; 13(20):5156. https://doi.org/10.3390/cancers13205156

Chicago/Turabian Style

Tron, Laure, Mathieu Fauvernier, Anne-Marie Bouvier, Michel Robaszkiewicz, Véronique Bouvier, Mélanie Cariou, Valérie Jooste, Olivier Dejardin, Laurent Remontet, Arnaud Alves, and et al. 2021. "Socioeconomic Environment and Survival in Patients with Digestive Cancers: A French Population-Based Study" Cancers 13, no. 20: 5156. https://doi.org/10.3390/cancers13205156

APA Style

Tron, L., Fauvernier, M., Bouvier, A. -M., Robaszkiewicz, M., Bouvier, V., Cariou, M., Jooste, V., Dejardin, O., Remontet, L., Alves, A., FRANCIM Group, Molinié, F., & Launoy, G. (2021). Socioeconomic Environment and Survival in Patients with Digestive Cancers: A French Population-Based Study. Cancers, 13(20), 5156. https://doi.org/10.3390/cancers13205156

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