The expected burden of mesothelioma mortality in Great Britain from 2002 to 2050.

The British mesothelioma register contains all deaths from 1968 to 2001 where mesothelioma was mentioned on the death certificate. These data were used to predict the future burden of mesothelioma mortality in Great Britain. Poisson regression analysis was used to model male mesothelioma deaths from 1968 to 2001 as a function of the rise and fall of asbestos exposure during the 20th century, and hence to predict numbers of male deaths in the years 2002-2050. The annual number of mesothelioma deaths in Great Britain has risen increasingly rapidly from 153 deaths in 1968 to 1848 in 2001 and, using our preferred model, is predicted to peak at around 1950 to 2450 deaths per year between 2011 and 2015. Following this peak, the number of deaths is expected to decline rapidly. The eventual death rate will depend on the background level and any residual asbestos exposure. Between 1968 and 2050, there will have been approximately 90,000 deaths from mesothelioma in Great Britain, 65,000 of which will occur after 2001.

Mesothelioma is a formerly rare cancer that principally affects the pleura and the peritoneum (Greenberg and Lloyd Davies, 1974) and is almost always caused by asbestos exposure. The disease is rapidly fatal, most of those affected dying within a year of diagnosis (Peto ). There is a long latent period between first exposure to asbestos and diagnosis of mesothelioma that is seldom less than 15 years and often exceeds 60 years (Bianchi ). In all, 85% of deaths are among men, and the risk is highest in occupations with substantial exposure to asbestos (McElvenny ).

Previous predictions of future numbers of mesothelioma deaths in Britain were based on observed male deaths between 1968 and 1991 (Peto ; Hodgson ) and were derived from a simple birth cohort model in which mesothelioma risk was related independently to age and to date of birth. These analyses suggested that male mesothelioma deaths would peak at about 2700–3300 deaths around 2020. The data have conformed to this model well up to 1991, but from the 1990s show a departure from fit for some of the later cohorts (Figure 1A and B). The slower increase with age for mesothelioma death rates in the 1990s implies that the simple age–cohort model is not reliable. We fitted the age–cohort model as confirmation of these observations, but our main models related current mortality to past asbestos exposure levels and did not assume the same age distribution of mortality in different birth cohorts.

Figure 1

(A) Mesothelioma death rates by 5-year cohort and age group for males born 1878–1937 (note: The earlier cohorts are to the lower right of the chart). (B) Mesothelioma death rates by 5-year cohort and age group for males born 1938–1982 (note: The earlier cohorts are to the upper left of the chart).

MATERIALS AND METHODS

The mesothelioma register is described in detail elsewhere (McElvenny ). Briefly, deaths occurring in England, Scotland and Wales, where mesothelioma was mentioned on the death certificate, are notified annually to the Health and Safety Executive by the Office for National Statistics and the General Register Office for Scotland. Originally, this was done manually as part of the coding of causes of death, but since 1992 (1996 in Scotland), the computerised national mortality database has been scanned. Checks for completeness include searching deaths in England and Wales for alternative spellings of mesothelioma and flagging cancer registrations with a morphology code consistent with mesothelioma for death notification at the National Health Service Central Registers in Edinburgh and Southport.

Details of the analysis, which was based on deaths in men aged 20–89 years between 1968 and 2001, are given in Appendix A1. The model parameters were estimated iteratively, and the deviance (based on observed and expected numbers aggregated over 5-year groups by age, date of birth and year of death) was used to assess overall goodness of fit. Two of the model parameters, k (the power of time defining the increase in risk after exposure) and H (the half-life for clearance of asbestos from the lung), are closely correlated and cannot be independently estimated. The effect of reducing the half-life is to increase the best-fit value of k, but the fit of the model is affected only slightly. Since fitting the model with both H and k is unstable, we used two versions of the model, one with (effectively) no clearance (H=1000 years) and the other with a clearance half-life of 15 years – a value suggested from the modelling of mortality of the Wittenoom workforce (Berry, 1991).

The adequacy of the model was tested by examining deviance residuals and by comparing observed and fitted numbers of mesothelioma deaths. The uncertainty of the estimate was quantified by calculating an approximate 95% confidence (CI) interval for the peak in predicted future number of mesothelioma deaths and the year in which this peak was expected to occur. Our estimate of the predicted peak was scaled to include mesothelioma deaths in men aged 90 years or older, and deaths among women.

Mesothelioma mortality is very low up to 20 years after first asbestos exposure even among heavily exposed workers (Peto ), so death rates up to 2001 give virtually no information on exposures in the previous 20 years. The level and timing of the predicted peak in mortality are virtually unaffected by exposure since 1980, but for longer-term prediction recent exposures become increasingly relevant, and some assumption must be made about asbestos exposure since the 1970s. An assessment of current exposure levels for a recent HSE Regulatory Impact Assessment (HSE, 2002) suggested that population exposure in the 2000 was around 4% of the peak value reached in the 1960s. For the present projections, we assumed a continuing decline in asbestos exposure, from 4% of the peak level in 2000 to 2% by 2010 and 0.75% by 2050.

RESULTS

The annual number of mesothelioma deaths in Great Britain has risen increasingly rapidly by about 12-fold from 1968 (the first year in which ascertainment is believed to be complete) to 2001 (Figure 2). The mesothelioma death rate in males has continued to increase in older age groups, but has decreased among younger men in recent years (Figure 3). In 2001, there were 1579 male deaths, 85% of the total, the majority being at ages 60–79 years, with relatively few aged less than 50 years. Less than 1% of deaths in men have occurred at age 90 or older, but the number is likely to increase substantially over time. Men born around 1940 have suffered higher death rates than any previous or subsequent birth cohort (Figures 1A and B), and they will not reach 90 years of age until 2030.

Figure 2

Mesothelioma deaths by sex and year.

Figure 3

Average annual male mesothelioma death rates per million by age and time period.

Parameter values estimated from the two versions of the Poisson regression model are shown in the Table 1. The two versions of the model produced effectively equivalent fits overall, and showed very similar patterns of residuals. This is illustrated (for the nonclearance model) in Figures 4A–C which show plots of fitted and observed deaths by year of birth, age and year of death. These show close agreement. A two-dimensional array (not shown) of deviance residuals by 5-year age and year of birth groupings also showed an unbiased pattern with no strong clustering of residuals of the same sign. While the fit of the clearance and nonclearance versions of the model to past data were equivalent, their predictions of future mortality were different. Under the nonclearance model, the annual number of mesothelioma deaths in men aged below 90 years is predicted to reach a peak of 1846 deaths (95% CI 1650–2100) in 2013 (95% CI 2011–2015). The clearance model predicts a peak of 1983 deaths (95% CI 1835–2233) in 2011 (95% CI 2009–2013). As more data accumulate on the pattern of reduced mortality related to exposure levels after the mid-1960s peak, it may be possible to estimate the parameters k and H separately.

Table 1

Estimated parameter values from optimal Poisson's regression model

Parameter	Estimates for non clearance model (half life=1000 years)	Estimates for model with clearance half life of 15 years
Power of time since first exposure, k	2.6	4.1
Year of exposure maximum	1967	1962
Diagnostic trend (percentage decrease in missed cases per year), DxT	5	2
Relative exposure potential for age group (years)
0–4	0.00	0.00
5–15	0.03	0.07
16–19	0.21	1.18
20–29	1.00 (baseline)	1.00 (baseline)
30–39	1.24	1.95
40–49	1.11	1.02
50–59	0.00	0.08
60–64	0.00	0.00
65+	0.00	0.00
Change in exposure index (percentage per year) in peak year±X years	% change	% change
−45	29	30
−35	6	8
−25	11	11
−15	9	9
−5	5	10
0 (peak year)	0	0
+5	−14	−17
+15	−39	−4

Figure 4

(A) Observed and fitted deaths by year of birth. (B) Observed and fitted deaths by age. (C) Observed and fitted mesothelioma deaths by year of death, with derived exposure index. (D) Number of observed and fitted mesothelioma deaths for 1955–1980 yearly birth cohorts.

Allowing for deaths occurring in men aged 90 or over, the estimated peak among males is 1857. Assuming a linear relationship between annual numbers of male M and female F mesothelioma deaths (F=22.9+0.136M; R2=0.95), the maximum number of female deaths was estimated to be 276, giving a predicted total of 2133 mesothelioma deaths (95% CI 1950–2450) at the peak of the epidemic. The total number of mesothelioma deaths to 2050 is predicted to be around 90 000, with 65 000 of these occurring from 2002 onwards.

DISCUSSION

The age–cohort model was fitted (residual deviance 121.5 on 72 degrees of freedom) and confirmed the observed departure of the data from that model by providing a significantly worse fit than our preferred model (residual deviance 235.6 on 182 degrees of freedom). The existence of two versions of the model (one with clearance and the other with no clearance) with equivalent overall fits but with different future implications raises the obvious question as to which is closest to the truth.

One way of assessing which model is the most realistic is to match the implied exposure patterns from each of the models with the actual pattern of asbestos imports, taking account of the difference between fibre types. Figure 5A shows the fitted exposure index for the nonclearance model together with the best approximating linear combination of the three fibre-specific import profiles, and Figure 5B shows this information for the lung clearance model (chrysotile had zero weight in both). It is evident that overall approximation is less good for the clearance than for the nonclearance exposure index. Although this argument is not definitive, we interpret the better agreement of imports with the exposure profile associated with the nonclearance model as justifying a modest preference for this model in developing our current projections. This preference is reinforced by the fact that the observed 2001 total number of mesothelioma deaths is a more extreme outlier in relation to the clearance model (1563 observed vs 1463 expected, P=0.0022) than in relation to the nonclearance model (1488 expected, P=0.0143).

Figure 5

(A) Fitted import index for nonclearance model with best approximating weighting of actual import series (amosite+0.7 × crocidolite. (B). Fitted import index for clearance model with best approximating weighting of actual import series (amosite+2.2 × crocidolite).

Our assumption that exposure in 2000 was 4% of the peak level hardly affects the predicted peak in mesothelioma mortality. It does, however, influence the predicted number of cases following the peak. Thus, the total burden of mesothelioma mortality to 2050 remains very uncertain. The estimated power of time since first exposure of 2.6 is in the range expected on the basis of fits of similar models to cohort data, providing a measure of support for the application of the model at that population level. The diagnostic trend, a decrease of 5% per year in the proportion of cases that are undiagnosed, implies that in 1968 about 90% of cases were diagnosed. This parameter improves the fit but has virtually no influence on the projections.

The model was fitted by minimising the total deviance. Recent birth cohorts have fewer mesothelioma deaths, so the fitted model is dominated by earlier cohorts with longer follow-up. The 1920–1924 cohort has the most deaths. Residuals for the most recent cohorts therefore provide an indicator of predictive reliability. Figures 4A and D show observed and fitted deaths by year of birth, with some divergence for the most recent births (Figure 4D). From 1965 onwards there are eight deaths compared to a fitted value of 3.5, a borderline significant excess (P=0.054). However, the fitted model does not include an allowance for a background rate of mesothelioma. If there were, as widely assumed, around 1–2% of male mesothelioma deaths not due to asbestos exposure, this equates to around 50 male deaths annually in Great Britain. On this basis, around four deaths would be expected among those not exposed to asbestos and born from 1965 onwards, bringing observed and fitted into good agreement. As data accumulates in these more recently born groups, it will be important to make explicit allowance in the model for the possibility of a background rate. In terms of projecting the timing and level of peak mortality, however, the numbers generated by this background rate can be ignored.

Predictions of the future number of mesothelioma cases have been attempted in other countries. For example, in Denmark, the 1912 cases of malignant mesothelioma reported to the Danish cancer registry between 1943 and 1993 were used to predict a peak incidence of 93 cases among men born before 1955 to around the year 2015 (Kjaergaard, 2000). In the United States, using mesothelioma incidence data from the Surveillance, Epidemiology and End Results programme, it has been estimated that there will be a peak around the years 2000 to 2004 of approximately 2000 cases and a return to background incidence by 2055 (Price and Ware, 2004). In Australia, the incidence of mesothelioma is expected to peak around 700 cases per year in 2010 (Leigh and Driscoll, 2003). In Sweden, the preventive measures of the mid-1970s can probably not be evaluated with reasonable precision until around 2005 (Jarvholm ). In the Netherlands, it has been predicted that pleural mesothelioma will peak around the year 2028, with up to 900 cases per year (Segura ). In France, the number of deaths is predicted to reach a peak at about 2200 per year some time after 2020 (Ilg ). Based on combined data from Britain, France, Germany, Italy, the Netherlands, Switzerland and Hungary, it has been suggested that the number of men dying from mesothelioma will almost double over the next 20 years, peaking at about 90 000 cases around 2018 (Peto ). These projections – like the earlier British projections (Peto ) – were based on a simple age and birth cohort model. More recent death rates in three of these countries (France, Germany and Italy) are lower than were predicted, and the timing and level of the peak of Europe-wide mesothelioma deaths may prove to be earlier and lower, respectively, than previous projections suggested.

CONCLUSIONS

The data for 2001 showed around 1850 mesothelioma deaths among males and females in Great Britain. Using a statistical modelling approach, mesothelioma mortality in Great Britain is predicted to peak at around 1950–2450 deaths per year some time between 2011 and 2015. Around 90 000 deaths are predicted to occur by 2050, with 65 000 of these occurring from 2002 onwards.