Open Access
How to translate text using browser tools
23 April 2024 Updated Analysis of Radon Exposure and Lung Cancer Mortality in the Cohort of Newfoundland Fluorspar Miners (1950–2016)
Paul J. Villeneuve, Howard I. Morrison, Rachel Lane
Author Affiliations +
Abstract

The commercial mining of fluorspar in St. Lawrence Newfoundland began in 1933. Miners who worked underground were exposed to high levels of radon progeny, especially before ventilation was introduced into the mines in 1960. The mean cumulative radon exposure for underground miners in this cohort was 380.9 working level months (WLM). A series of studies of this cohort have characterized the increased risks of lung cancer mortality due to radon. We have extended the follow-up of this cohort an additional 15 years to provide additional insights on the risks of low levels of radon exposure, and the modifying effects of time since exposure, age at first exposure, attained age, duration of exposure, and cigarette smoking. The cohort consisted of 1,735 underground and 315 male surface miners who, combined, accrued 81,650 person-years of follow-up. The mortality experience of the cohort was determined from 1950–2016 through record linkage to Canadian national death data. Individual-level estimates of exposure to radon progeny, in WLMs, were determined for each year of employment. We compared the mortality experience of the underground miners to Newfoundland men using the standardized mortality ratio (SMR). Poisson regression models were fit to estimate excess relative risks (ERR) per 100 WLM. There were 236 lung cancer deaths identified, and of these, 221 occurred among underground workers. The SMR for lung cancer among underground miners compared to Newfoundland men was 2.67 (95% CI: 2.33, 3.04). The ERR per 100 WLM for lung cancer mortality, assuming a 5-year exposure lag, was 0.41 (95% CI: 0.23, 0.59). Attained age and time since exposure were important modifiers to the radon-lung cancer relationship. The joint relationship between smoking and radon on lung cancer risk was sub-additive, however, the smoking data were limited and available for only half of the cohort.

INTRODUCTION

Radon and its radioactive decay products (radon progeny)2 are established carcinogens, and in Canada, residential radon exposure was estimated to have caused 6.9% of incident lung cancers diagnosed in 2015 (1). The overall number of lung cancers in Canada due to occupational exposure to radon has been estimated to be substantially lower (0.8%), and workplace exposure to radon is highly variable across occupational groups (2). Underground miners, particularly for uranium, have substantially higher levels of exposure to radon than the background concentrations that the general population is exposed to (3). Epidemiological studies of underground miners have played a prominent role in the decision to classify radon as a human carcinogen (4).

There are relatively few cohorts of radon-exposed miners in the world, and the Newfoundland fluorspar miners cohort has been included in previous pooled analyses that have explored variations in lung cancer risk by time since exposure, age at first exposure, attained age, duration of exposure, as well as the role of cigarette smoking (4, 5). The cohort of Newfoundland fluorspar miners has several important features including radon exposures that are higher than most other cohorts included in the BEIR VI report (4), and other pooled analyses (5). Another important feature of this cohort is that, unlike other cohorts of radon-exposed miners, the source of exposure to radon was from the ground water that flowed through the mines, and not from the ore. As a result, the fluorspar miners were not exposed to the other recognized lung carcinogens of uranium miners: gamma, thoron, and radioactive dusts.

The mining of fluorspar in St. Lawrence, Newfoundland began in 1933, and continued until the closure of the mines in 1978. Mining operations for fluorspar in this area recently restarted in 2018; however, operations again ceased in 2022. The analyses presented in this paper only considered those miners who worked until 1978. When the mines first opened in the 1930s, they provided employment opportunities to coastal communities devastated by a 1929 tsunami and the challenges of residents of these towns are well-chronicled (6, 7). The fluorspar mines were especially active during the Second World War, and the mined ores were used to produce aluminum. Radon exposures in the series of these underground mines were quite high, particularly before the introduction of mechanical ventilation in 1960. For context, the current recommended occupational levels for radon exposure is 4 working level months (WLM) per year (3). Before 1950, the typical estimated median annual radon exposure among underground miners was above 100 WLM, and the maximum annual exposure was 564.7 WLM.

In St. Lawrence, a local physician first noticed an increased incidence of lung cancer in the 1950s. More formal assessments of the increased risks of lung cancer due to occupational radon exposure in these miners were published by Parsons et al. in 1969 (8) and then Wright and Couves in 1977 (9). We published a series of studies thereafter (1012) with the most recent publication including deaths identified through December 31, 2001 (10). In this paper, we present findings from an updated analysis of this cohort that includes an additional 15 years of mortality follow-up; 67 years of mortality follow-up overall (1950–2016).

We were motivated to extend the follow-up of this cohort for several reasons. First, the additional 15 years of follow-up provided us the opportunity to better describe changes in lung cancer risk with an increased number of years since last exposure. Second, the size of the cohort and the number of previously identified lung cancer deaths is relatively small, and the identification of additional outcomes increases the precision of our risk estimates. Third, the increased follow-up provided an opportunity to revisit the risks of lung cancer from radon among those miners employed after 1960 when levels were much lower.

MATERIALS AND METHODS

Study Population

Previous publications have described the study cohort (10, 12). In brief, occupational records were used to create the initial cohort of fluorspar miners, and it included 2,661 miners. Of these, 2,111 miners worked underground at some point, while the remaining 550 worked exclusively on the surface. All individuals included in these analyses were male. Annual working histories were constructed from employment records. We excluded from analysis miners who lacked sufficient personal identifying information needed for record linkage to national mortality data. As detailed in our previous papers (10, 12) most of these were short-term workers employed during World War II. We reconstructed the cohort using archived data stored by the Canadian Nuclear Safety Commission.

Radon Exposure

An annual estimate of cumulative radon exposure was derived for each miner. These exposures incorporated details about mine location, and occupational data on working hours and location from 1933–1960 that were available for each cohort member (13). Radon exposure estimates were based on sampling measures (14), as well as the general mine layout, amount of water inflow, mining method, and general working conditions (13).

Over time, there were important improvements in the working environments including going from three to two eight-hour shifts per day (at the end of World War II) with four hours between shifts. Between shifts, compressed air was used to provide some ventilation. In 1960, before the introduction of mechanical ventilation in the mines, 80 samples of radon were taken at 50 different underground locations in the Director mine (the principle mine). Values ranged from 0.4 to 183 Working Levels (WLs) in unventilated areas, and from below the level of detection to 12 WLs in ventilated areas of the mine. Unfortunately, there are no realistic means to verify the accuracy of radon exposures before 1960. The introduction of mechanical ventilation resulted in large decreases in radon exposure. For example, from 1955–1959, the mean annual exposure to radon for underground workers was 64.8 WLMs, whereas from 1961–1965 it was 1.3 WLMs. From 1986 onwards, radon levels for each worker were measured daily based on the site location in the mine. The mean exposure for underground miners from 1968–1978 was 0.6 WLM. Cumulative exposures were expressed in working level months (WLMs).3

Previous analyses of the cohort incorporated a lag interval of 5 years into the calculation of cumulative exposure (5), and we followed this approach herein. This assumes that the risk of lung cancer mortality is unrelated to radon exposures received in the five years before death.

Cigarette Smoking

Miners provided smoking information by completing a series of smoking surveys. We administered the last smoking survey in 2003 to update previously collected information on smoking behaviors and sought to identify former miners with no smoking data. Smoking surveys were previously administered in 1966, 1970, 1978, and 1993. For the 2003 survey, we identified participants using the personal identifying information (name, birth date, years worked in mine, and place of residence) from past smoking surveys. Internet search engines were used to help interviewers (hired from the St. Lawrence region) to locate former miners. For this current paper, we compiled data from all five smoking surveys to classify miners based on last known smoking status (current, former, never), and according to the average number of cigarettes smoked per day. When the reported number of cigarettes smoked per day varied across surveys, we derived the average of these values of all surveys. Smoking data were available for approximately half the cohort. It is important to also note for miners with smoking data, data were usually only collected from a single smoking survey. Therefore, it was not possible to construct time-dependent measures of smoking behaviors, nor model lung cancer risk in relation to the time since smoking cessation.

Ascertainment of Mortality (1950–2016)

We created a linkage file with personal identifying information using archived data maintained by the Canadian Nuclear Safety Commission. This file combined data from occupational records, smoking surveys, and vital status from previous linkages. This file contained 2,121 observations each of which represented a unique miner. Statistics Canada staff linked these individual-level data to the Canadian Vital Statistics Death Database (CVSD) and the Canadian Mortality Database (CMDB). The CVSD is an administrative survey that collects demographic and cause of death information annually from all provincial and territorial vital statistics registries on all deaths in Canada (15). Death records were sourced from both the Canadian Mortality Database (CMDB) (data from 1950–2011 inclusive) (16), and Canadian Vital Statistics Database (CVSD) that covered the periods from 1930–1949, and 2011 onwards. Linkage assignments employed both probabilistic and deterministic methodologies, using the “matchit” program on the STATA 14.2 platform. Following the identification of deaths from the record linkage, the lead author (PV) performed a manual review of all death links. He had previously reviewed mortality linkages for past analyses of this cohort that were done in 1994 and 2005.

Miners for whom no death link was found were assumed to be alive at the end of follow-up (December 31, 2016). We made an exception to this rule for miners who were born before 1910 with no link. Namely, these miners (n = 55) were assumed deceased but lost to follow-up and we excluded them from analyses. The final cohort analysis file consisted of 2,055 miners, which is slightly less than the 2,070 miners that formed the basis of our previous analyses (10). The ability to link occupational cohorts to national mortality has been previously demonstrated (1719). We obtained date of death, and underlying cause of death from the record linkage. The International Classification of Diseases 9th Revision (ICD-9) determined the underlying cause of death for deaths before 2000, while the ICD-10 classified deaths from 2000 onwards. Lung cancer deaths were identified with the ICD-9 rubric of 162, and ICD-10 codes of C33 and C34.

We adhered to Statistics Canada's rules of disclosure governing the use of files linked to national mortality. We performed all analyses within the Statistics Canada's Research Data Center Network (20). We were bound to adhere to the rule that all frequencies, including deaths, would be rounded to the base unit of five. For this reason, unlike previous analyses of this cohort, all frequencies presented and referred to in this paper have been rounded to base five. However, the risk estimates derived from the regression models were not subject to rounding.

External Comparisons of Mortality to the General Male Newfoundland Population

We applied the method of indirect standardization to compare the patterns of mortality for the underground fluorspar miners to the general male population of Newfoundland from 1950–2016. Analyses were restricted to the underground miners because we wanted to describe the relative differences in mortality between those with above background levels of exposure to radon to the general population. We obtained cause-specific mortality rates for the Newfoundland male population from Statistics Canada, and tabulated these rates by 5-year age-groupings and 10-year calendar periods. These cause-specific rates were multiplied by the corresponding number of person-years of follow-up accrued over the same calendar and age-group stratum among the underground fluorspar miners to obtain the expected number of deaths. For each cause of death, we calculated the standardized mortality ratio (SMR) by dividing the observed by the expected number of deaths, and estimated the accompanying 95% confidence interval using formulae published by Breslow and Day (21).

Internal Cohort Comparisons of Mortality

Internal cohort analyses explored variations in risk across different features of radon exposure. These analyses are more powerful than external comparisons as they provide a means to characterize exposure-response relationships while adjusting for the influence of other individual-level risk factors including cigarette smoking. They also are less prone to bias from the healthy worker effect than comparisons to the general population (22). The standard approach to characterize risks in these cohorts is the calculation of an excess relative risk (ERR) per 100 WLM using Poisson regression. Person-years, cross-classified across a series of categorical variables, were tabulated using a program adapted from that published by Pearce and Checkoway (23). The person-years of follow-up and number of lung cancer deaths were tabulated across the following categorical variables: attained age, age at first exposure, cumulative radon exposure, calendar period, smoking status, number of cigarettes smoked daily, and a binary variable that captured the periods before and after mechanical ventilation was introduced into the mine. All variables were modelled as categorical variables. For cumulative exposure to radon, we used the following categories: 0, > 0–<2.5, 2.5–<5, 5–10, 10–<25, 25–<75, 75–<150, 150–<300, 300–<600, 600–<1,200, 1,200–<1,800, 1,800–<2,500, 2,500–<3,500, and ≥3,500 WLM. We chose these categories so that there were enough lung cancer deaths (n >10) in each category. For each exposure category, we calculated a mean exposure weighted by the person-years of follow-up. The number of deaths in each stratum were Poisson distributed with the variance estimated by the observed number of deaths.

The AMFIT program of the software package EPICURE was used for all Poisson regression modeling (24). We first estimated relative risks for each cumulative WLM category with adjustment for attained age and calendar period. Linear excess relative risk models were fit using the mean WLM generated for each cumulative WLM category. Mathematically, the form of this model was: RR = 1 + βWLM, where β describes an increase in the ERR per unit increase in cumulative WLM. All linear ERR models were adjusted for attained age and calendar period.

Effect Modification

Similar to methodology used in other cohorts of radiation exposure workers (4), we evaluated variations in the ERR/WLM across a series of potential effect modifiers that included: attained age, age at first exposure, exposure rate, time since last exposure, and duration of exposure. This was done by first fitting the linear ERR model: RR = 1 + βWLM. We evaluated heterogeneity in the ERR/WLM across categories of the potential effect modifier by fitting the model RR = 1 + βjWLM where j represents the number of categories of the modifying factor. The ERR/ WLM for each of the j categories was obtained from this model, and we assessed whether these risks were heterogeneous across these categories by deriving a test of significance for the improvement in fit based on the likelihood ratio test statistic with j-1 degrees of freedom. We evaluated effect modifiers one at a time because, given the small number of lung cancer cases and the nature of the cross-classified data, it was not possible to test for effect modification in more complex models.

Joint Effects of Smoking and Radon on Lung Cancer Mortality

We evaluated the associations between smoking, radon exposure and lung cancer following the methodological approach by Lubin et al. (5). We fitted multiplicative and additive main effects models, as well as a full model. The form of the joint association between smoking and WLM was assessed relative to a joint relationship that was either multiplicative or additive. As described by Lubin (25), with a linear ERR form for RR, a multiplicative joint association for lung cancer mortality would be expressed as the following:

e01_59.gif

whereas an additive RR association would be given by the following:

e02_59.gif

We then considered models where both the multiplicative and additive associations are nested. One such model was the geometric mixture model. This is defined as the following:

e03_59.gif

We proceeded to fit a series of models with different values of λ until we identified a model that minimized the deviance from the full model. The value of fi01_59.gif from this iterative approach was used to characterize the joint association between radon and smoking on lung cancer risk.

Funding and Ethics

The Canadian Nuclear Safety Commission funded this study. The Carleton Research Ethics Board provided ethics approval (Clearance no.: 108564). Finally, Statistics Canada reviewed and approved the study as part of their record linkage processes.

RESULTS

Table 1 describes the characteristics of 1,735 underground and 315 surface fluorspar miners employed from 1933–1976. A total of 1,365 deaths were identified in the cohort through December 31, 2016; of these, 235 were from lung cancer (220 among underground miners, and 15 among surface miners). This represents an additional 30 deaths from lung cancer since our last analyses. Smoking status was available for approximately 54% of the cohort, and of these individuals, 85% were ever-smokers.

TABLE 1

Descriptive Characteristics of the Newfoundland Fluorspar Cohort with Mortality Follow-up from 1950–2016

img-AqhJ_59.gif

After accounting for age and calendar period effects, the overall mortality experience of underground miners in the cohort was not different from Newfoundland men (SMR = 1.03; 95% CI: 0.97, 1.09) (Table 2). The SMR for all cancer deaths was statistically significant (SMR = 1.50; 95% CI: 1.36, 1.65), due predominantly to the excess lung cancer deaths at 2.67 (95% CI: 2.33, 3.04), but there was not a statistically significant excess of mortality observed for any other specific cancer site examined. Increased risks of mortality, relative to the provincial population, were also found for accidents, poisonings and violent deaths (SMR = 1.31; 95% CI: 1.02, 1.66). This excess appears to be largely driven by a higher-than-expected number of deaths for falls (SMR = 2.50, 95% = 1.40, 4.13).

TABLE 2

Standardized Mortality Ratios (SMRs) for Newfoundland Underground Fuorspar Miners Relative to the Mortality Experience of Newfoundland Men, 1950–2016

img-AWS_59.gif

Internal cohort analyses found that cumulative radon exposure was strongly associated with lung cancer mortality (Fig. 1). Those in the highest exposure category (>3,500 WLM) had a statistically significant relative risk of lung cancer mortality of 39.4 (95% CI: 19.1, 81.0) relative to surface miners. The ERR among underground miners, in the model adjusted for age and calendar period, per 100 WLM was 0.41 (95% CI: 0.23, 0.59). When analyses were restricted to cumulative exposures less than 100 WLM, the strength of the association was much stronger, specifically, the ERR per 100 WLM was 2.90 (95% CI:0.60, 5.20) (Fig. 2). The ERR per 100 WLM decreased substantially with increasing length of follow-up (Fig. 3).

FIG. 1

Relative risk (RR)of lung cancer mortality from cumulative radon exposure in working level months (WLMs), Newfoundland Fluorspar cohort, 1950–2016. RRs (and 95% CIs) are plotted at the mean WLM for each exposure category and were adjusted for age and calendar period. The solid line was created from fitting model RR = 1 + ERR(WLM) with adjustment for attained age and calendar period. The dashed lines represent the upper and lower 95% confidence intervals for the ERR. The ERR per 100 WLM was 0.41 (95% CI: 0.23, 0.59).

img-z6-1_59.jpg

FIG. 2

Relative risk (RR) of lung cancer mortality from cumulative radon exposure less than 100 working level months (WLMs), Newfoundland Fluorspar cohort, 1950–2016. RRs (and 95% CIs) are plotted at the mean WLM for each exposure category and were adjusted for age and calendar period. The solid line was created from fitting model RR = 1 + ERR(WLM) with adjustment for attained age and calendar period. The dashed lines represent the upper and lower 95% confidence intervals for the ERR. The ERR per 100 WLM was 2.90 (95% CI: 0.60, 5.20).

img-z6-8_59.jpg

FIG. 3

Excess Relative Risk (ERR) per 100 WLM by follow-up period, Newfoundland Fluorspar miners, 1950–2016. ERR/100 WLMs were adjusted for 5-year age groupings, cigarettes smoked daily, and calendar periods.

img-z7-1_59.jpg

Effect Modifiers

We found no evidence of heterogeneity between cumulative radon exposure and lung cancer deaths across age at first exposure categories (P = 0.39) (Table 3). In contrast, the ERR\100 WLM with attained aged decreased substantially with increasing age (P = 0.001). Similarly, the ERR\WLM decreased with increasing time since last exposure categories (P < 0.001) and those last exposed more than 35 years ago had increased excess risk (ERR\100 WLM = –0.24; 95% CI: –0.82, 3.29). The ERR\100 WLM increased with increasing duration of exposure, which is supportive of a pattern of inverse dose response, but this finding was not statistically significant (P = 0.11). While a smaller ERR\100 WLM was observed among those first employed after 1960 compared to those before this year, this difference was not statistically significant (P = 0.07).

TABLE 3

ERR/100 WLM Across Categories of Attained Age, Time Since Last Exposure, Age at First Exposure, Exposure Rate and, Year First Employed Among Underground Newfoundland Fluorspar Miners, 1950–2016

img-ATus_59.gif

We present the risks of lung cancer by time windows of exposure to radon in Table 6. The ERR\100 WLM was much stronger for exposure received more recently. Specifically, the ERR\100 WLM for exposure received 5–14 years ago was 1.84 (95% CI: 0.19, 3.48), while for exposures received 35 years ago the ERR\100 WLM was 0.08 (95% CI: –0.00, 0.17).

TABLE 6

Excess Relative Riska (ERR/100 WLM) of Lung Cancer Mortality by Time Since Exposure Windows, Newfoundland Fluorspar Underground Miners, 1950–2016

img-Aj2_59.gif

The ERR\100 WLMs were similar between ever- and never-smokers (P > 0.50) (Table 4). In contrast, heterogeneity in the ERR\100 WLM was observed with cigarettes smoked daily (P < 0.05). For the latter, higher ERR\100 WLM were observed among those who smoked on average a pack or more per day. The joint effect of cigarette smoking and radon exposure was found to be sub-additive with the best estimate of fi01_59.gif = –0.36 (Table 5).

TABLE 4

Excess Relative Risk (ERR/100 WLM) of Lung Cancer, per WLM, Across Categories of Number of Cigarettes Smoked Daily and Smoking Status, Newfoundland Fluorspar Cohort, 1950–2016

img-AOrc_59.gif

TABLE 5

Relative Risk (RR) of Lung Cancer by Cumulative Radon Exposure (WLM) and Cigarettes Smoked Daily, Newfoundland Fluorspar Miners, 1950–2016

img-Ax31_59.gif

DISCUSSION

This paper is another update, and likely the last, of the risk of lung cancer among Newfoundland fluorspar miners occupationally exposed to high levels of radon exposure. Many of the patterns of risk are consistent with previously published values. However, the extension of the mortality follow-up of this cohort provides increased precision in the risk estimates, particularly given the relative size of the cohort and corresponding number of lung cancer deaths. As with other cohorts of radon-exposed miners, the ERR\100 WLM for lung cancer mortality has decreased since the last set of analyses. However, this change was quite modest, specifically, in our previous analyses the ERR per 100 WLM for all underground miners was 0.43 (95% CI: 0.23, 0.62) (10), while in our current study the value is 0.41 (95% CI: 0.23, 0.59). Our earlier analyses involving mortality through the end of 1990 found a value of 0.70 per 100 WLM (95% C: 0.44, 1.11) (26). While approximately one-third of the cohort was determined to be alive at the end of the follow-up interval, it is our view that further extensions of the follow-up are unlikely to materially change this overall ERR\100 WLM.

A positive association between cumulative radon exposure and lung cancer risk was observed when analyses were restricted to those with less than 100 WLM of lifetime exposure. The ERR\100 WLM was 2.90 (95% CI: 0.60, 5.20). While we did observe a stronger ERR\100 WLM among workers who were first exposed after mechanical ventilation was introduced into the mines (7.3 vs. 0.44; P = 0.067), this estimate was based on only 35 lung cancer deaths. As these miners were younger relative to the other fluorspar miners (as they were employed in a more recent time), this ERR\100 WLM is expected to decrease with increased follow-up. This finding is consistent with other recent work (27).

As with previous analyses of this cohort, attained age was a strong modifier of the ERR\100 WLM. The ERRs per 100 WLM were 5.3, 0.6, 0.2, and 0.1 for those with attained ages of <50, 50–59, 60–69, and 70+, respectively. The decrease in ERR per unit exposure with increasing categories of attained age was similarly reported in the BEIR VI report and were 1, 0.66, 0.32, and 0.17 among those aged 45–54, 55–64, 65–74, and 75 years or older. In the more recent pooled analyses of the underground uranium miners (27), the decreases in the ERR\WLM for those aged <55, 55–64, 65–74, and 75+ years of age were 100, 0.64, 0.22 and 0.17.

A strength of our study was the availability of some smoking data that are often not available for many of the historical cohorts or radiation-exposed miners. Despite this, we acknowledge our smoking data have important limitations. We lacked the ability to characterize changes in smoking behaviors over time in our cohort as we typically did not have smoking data for each miner across multiple smoking surveys. This limitation precluded us from being able to reliably characterize “time since last smoked” for the miners for which we had smoking data. There were also some changes in the wording of the smoking surveys that were done as early as 1966 and as late as 2001 which presented challenges from harmonizing data. Where possible, we used data from multiple surveys to characterize smoking behaviors, particularly for the measure “average cigarettes smoked daily.” Despite these limitations, we are confident that the data collected in the surveys provides a valid mechanism for identifying ever- and never-smokers, and those who were heavier smokers.

We found no substantial difference in the ERR\WLM among ever-smokers compared to never-smokers, and there was no clear pattern in the ERR\WLM based on the average number of cigarettes smoked daily. A recent analysis of German Wismut uranium miners (28) found that the centered ERR/WLM was higher for non/light smokers relative to those who smoked more (0.022 vs. 0.013). In contrast, Kelly-Reif and colleagues found a substantially higher ERR per 100 WLM among smokers ((ERR/100WLM = 1.35, 95% CI: 0.84, 2.15) relative to never-smokers (ERR/100WLM = 0.12; 95% CI: –0.05, 0.49) in a cohort of Czech uranium miners (29). Our analysis of the association between cigarette smoking, radon and lung cancer revealed that the joint model was sub-additive. This is consistent with the work by Lubin et al. (5) who determined that the optimal fit for the geometric mixture model for this fluorspar cohort was λ = – 0:1: The value we found was λ = – 0:36: Additional support for a sub-additive model was found with the analyses of three European case-control studies of miners (30). The BEIR VII report concluded that the data on smoking, radon and lung cancer in miners studies were compatible with a multiplicative effect and not with an additive one, but that the estimated interaction was sub-multiplicative (31). Hornung et al. remarked (32) that with increasing follow-up, associations between smoking, radon and lung cancer have a tendency to move from multiplicative to additive. We had limited power to discriminate between competing models that included radon and smoking due to the small number of lung cancers and the high prevalence of smoking in our cohort, and because smoking data were only available for approximately half of all miners.

As previously mentioned, the source of exposure for these miners differed from that of other radon-exposed cohorts as radon came from water that ran through the mines. Nonoccupational exposure to radon was likely higher for the community relative to other communities based on water samples taken before 1960. Specifically, concentrations of radon in water samples in the Director mine, and Hares Ears shaft were 476 Bq/L and 445 Bq/L, respectively, while the corresponding concentrations in the town well were 532 Bq/L and 341 Bq/L (13). Residential radon measures taken in St. Lawrence in 1970 demonstrated substantial variability from home to home (33), and a more recent publication reported that indoor dwelling concentrations of radon in St. Lawrence Newfoundland were higher than for 16 other cities based on measures taken in the 1970s (34). The only other city with higher indoor concentrations of radon than St. Lawrence was Winnipeg with an arithmetic mean of 146 Bq/m3 (vs. 125 Bq/m3). While nonoccupational exposures to radon may be higher for residents of St. Lawrence Newfoundland when compared to other jurisdictions, we know of no reason why these exposures would be correlated with occupational exposures in the fluorspar cohort. This would be a condition necessary to introduce confounding.

In summary, the present study provides risk estimates like past analyses of this cohort, albeit with improved precision, and an ability to characterize risks with longer latency. Future analyses of this cohort that could be achieved by increasing the follow-up are unlikely to yield risk estimates for lung cancer that would be meaningfully different. Nonetheless, the updated data should be incorporated with other occupational radon-exposed cohorts when future pooled analyses are considered.

ACKNOWLEDGMENTS

The study was funded by the Canadian Nuclear Safety Commission. We would like to thank Dr. Toyib Olaniyan for his assistance with the record linkage of the cohort to Canadian national mortality data. The analysis took place at Statistics Canada's Research Data Centres (RDC) based at the University of Ottawa, and Carleton University. The RDC is supported by funds to the Canadian Research Data Centre Network (CRDCN) from the Social Sciences and Humanities Research Council (SSHRC), the Canadian Institute for Health Research (CIHR), the Canadian Foundation for Innovation (CFI), and Statistics Canada. Although the research and analysis are based on data from Statistics Canada, the opinions expressed do not represent the views of Statistics Canada. Statistics Canada vetted the data for confidentiality but had no role in the design, collection, analysis, and interpretation of the data.

©2024 by Radiation Research Society. All rights of reproduction in any form reserved.

REFERENCES

1.

Gogna P, Narain TA, O'Sullivan DE, Villeneuve PJ, Demers PA, Hystad P, et al., Estimates of the current and future burden of lung cancer attributable to residential radon exposure in Canada. Prev Med 2019; 122:100–08. Google Scholar

2.

2. Burden of occupational cancer in Canada: Major workplace carcinogens and prevention of exposure. Toronto: Ontario Occupational Cancer Research Centre; 2019. Google Scholar

3.

Canada C, Radon Profile. 2021. Google Scholar

4.

4. National Research Council CotBEoIR, Health Effects of exposure to radon: BEIR VI. Washington, D.C.: National Academy Press; 1999. Google Scholar

5.

5. Radon and Lung Cancer Risk: A joint analysis of 11 underground miners studies. Report No. NIH Publication No. 94-3644: U.S. Department of Health and Human Services, Public Health Service, National Institutes of Health; 1994. Google Scholar

6.

Rennie R, The St. Lawrence Fluorspar Mines: A Brief History. Place Heritage Newfoundland and Labrador: Heritage Newfoundland and Labrador. (  https://www.heritage.nf.ca/articles/economy/st-lawrence-mines.phpGoogle Scholar

7.

McIntyre L, The Wake: The Deadly Legacy of a Newfoundland Tsunami. Toronto, Canada: HarperCollins; 2019. Google Scholar

8.

Parsons WD, Devilliers AJ, Bartlett LS, Becklake MR, Lung cancer in a Fluorspar mining community. II. Prevalence of respiratory symptoms and disability. Br J Ind Med 1964; 21:110–6. Google Scholar

9.

Wright ES, Couves CM, Radiation-induced carcinoma of the lung–the St. Lawrence tragedy. J Thorac Cardiovasc Surg 1977; 74:495–8. Google Scholar

10.

Villeneuve PJ, Morrison HI, Lane R, Radon and lung cancer risk: an extension of the mortality follow-up of the Newfoundland fluorspar cohort. Health Phys 2007; 92:157–69. Google Scholar

11.

Morrison HI, Villeneuve PJ, Lubin JH, Schaubel DE, Radon-progeny exposure and lung cancer risk in a cohort of Newfoundland fluorspar miners. Radiat Res 1998; 150:58–65. Google Scholar

12.

Morrison HI, Semenciw RM, Mao Y, Wigle DT, Cancer mortality among a group of fluorspar miners exposed to radon progeny. Am J Epidemiol 1988; 128:1266–75. Google Scholar

13.

13. A retrospective study of radon daughter concentrations in the workplace in the fluorspar mines in St. Lawrence, Newfoundland. Ottawa: Atomic Energy Control Board; 1984. Google Scholar

14.

14. Dust levels in the mines of Newfoundland Ottawa: National Health and Welfare; 1958. Google Scholar

15.

15. Statistics Canada, Vital Statistics - Death Database. Place Statistics Canada: Statistics Canada; 2019. (  https://www23.statcan.gc.ca/imdb/p2SV.pl?Function=getSurvey&Id=1530624Google Scholar

16.

16. Statistics Canada, Microdata linkage program. 2019. (  https://www.statcan.gc.ca/en/record/2019Google Scholar

17.

Aronson KJ, Howe GR, Utility of a surveillance system to detect associations between work and cancer among women in Canada, 1965-1991. J Occup Med 1994; 36:1174–9. Google Scholar

18.

Goldberg MS, Carpenter M, Theriault G, Fair M, The accuracy of ascertaining vital status in a historical cohort study of synthetic textiles workers using computerized record linkage to the Canadian Mortality Data Base. Can J Public Health 1993; 84:201–4. Google Scholar

19.

Schnatter AR, Acquavella JF, Thompson FS, Donaleski D, Theriault G, An analysis of death ascertainment and follow-up through statistics Canada's Mortality Data Base system. Can J Public Health 1990; 81:60–5. Google Scholar

20.

20. Statistics Canada, Research Data Centres. Place Government of Canada: Government of Canada; 2020. (  https://www.statcan.gc.ca/en/microdata/data-centresGoogle Scholar

21.

Breslow NE, Day NE, Statistical methods in cancer research. Volume II–The design and analysis of cohort studies. IARC Sci Publ 1987; 1–406. Google Scholar

22.

Choi BC, Definition, sources, magnitude, effect modifiers, and strategies of reduction of the healthy worker effect. J Occup Med 1992; 34:979–88. Google Scholar

23.

Pearce N, Checkoway H, A simple computer program for generating person-time data in cohort studies involving time-related factors. Am J Epidemiol 1987; 125:1085–91. Google Scholar

24.

Preston DL, Epicure, version 2.0. Place Risk Sciences: Risk Sciences; 2015. Google Scholar

25.

25. Radon and Lung Cancer Risk: A Joint Analysis of 11 Underground Miners Studies. Bethesda, US: US National Institutes of Health; 1994. Google Scholar

26.

26. Radon progeny exposure and lung cancer risk: analyses of a cohort of Newfoundland fluorspar miners. Report No. INFO-0570. Ottawa: Atomic Energy Control Board; 1995. Google Scholar

27.

Richardson DB, Rage E, Demers PA, Do MT, Fenske N, Deffner V, et al., Lung cancer and radon: Pooled analysis of uranium miners hired in 1960 or later. Environ Health Perspect 2022; 130:57010. Google Scholar

28.

Kreuzer M, Sobotzki C, Schnelzer M, Fenske N, Factors modifying the radon-related lung cancer risk at low exposures and exposure rates among German uranium miners. Radiat Res 2018; 189:165–76. Google Scholar

29.

Kelly-Reif K, Sandler DP, Shore D, Schubauer-Berigan M, Troester M, Nylander-French L, et al., Lung and extrathoracic cancer incidence among underground uranium miners exposed to radon progeny in the Pribram region of the Czech Republic: a case-cohort study. Occup Environ Med 2022; 79:102–8. Google Scholar

30.

Leuraud K, Schnelzer M, Tomasek L, Hunter N, Timarche M, Grosche B, et al., Radon, smoking and lung cancer risk: results of a joint analysis of three European case-control studies among uranium miners. Radiat Res 2011; 176:375–87. Google Scholar

31.

31. Health Risks from Exposure to Low Levels of Ionizing Radiation: BEIR VII Phase 2. Washington, DC; 2006. Google Scholar

32.

Hornung RW, Deddens J, Roscoe R, Modifiers of exposure-response estimates for lung cancer among miners exposed to radon progeny. Environ Health Perspect 1995; 103 Suppl 2:49–53. Google Scholar

33.

McGregor RG, Vasudev P, Letourneau EG, McCullough RS, Prantl FA, Taniguchi H, Background concentrations of radon and radon daughters in Canadian homes. Health Phys 1980; 39:285–9. Google Scholar

34.

Chen J, Moir D, An updated assessment of radon exposure in Canada. Radiat Prot Dosimetry 2010; 140:166–70. Google Scholar

Notes

[1] 2 We refer to radon and its radioactive decay products (radon progeny) as “radon” synonymously throughout this article.

[2] 3 The concentration of radon progeny in workplace air is expressed in working levels (WL), where 1 WL is the concentration of radon progeny per liter of air that would result in the ultimate release of 1.3 × 105 MeV of potential alpha-particle energy. Occupational exposure to radon progeny is the product of time in the workplace and the concentration of radon progeny in the workplace air, measured in working level months (WLM), where 1 WLM is equivalent to one working month (170 hours) at a concentration of 1 WL.

Paul J. Villeneuve, Howard I. Morrison, and Rachel Lane "Updated Analysis of Radon Exposure and Lung Cancer Mortality in the Cohort of Newfoundland Fluorspar Miners (1950–2016)," Radiation Research 202(1), 59-69, (23 April 2024). https://doi.org/10.1667/RADE-23-00114.1
Received: 14 June 2023; Accepted: 9 April 2024; Published: 23 April 2024
Back to Top