Abstract
Lung cancer screening trials provide an opportunity to study the natural history of emphysema by using computed tomography (CT) lung density as a surrogate parameter.
In the Danish Lung Cancer Screening Trial, 2,052 participants were included. At screening rounds, smoking habits were recorded and spirometry was performed. CT lung density was measured as the volume-adjusted 15th percentile density (PD15). A mixed effects model was used with former smoking males with <30 pack-yrs and without airflow obstruction (AFO) at entry as a reference group.
At study entry, 893 (44%) participants had AFO. For the reference group, PD15 was 72.6 g·L−1 with an annual decline of -0.33 g·L−1. Female sex and current smoking increased PD15 at baseline, 17.3 g·L−1 (p<0.001) and 10 g·L−1 (p<0.001), respectively; and both increased the annual decline in PD15 (female: -0.3 g·L−1; current smoking: -0.4 g·L−1). The presence and severity of AFO was a strong predictor of low PD15 at baseline (Global Initiative for Chronic Obstructive Lung Disease (GOLD) I: -1.4 g·L−1; GOLD II: -6.3 g·L−1; GOLD III: -17 g·L−1) and of increased annual decline in PD15 (GOLD I: -0.2 g·L−1; GOLD II: -0.5 g·L−1; GOLD III: -0.5 g·L−1).
Female sex, active smoking and the presence of AFO are associated with accelerated decline in lung density.
The mortality and morbidity of chronic obstructive pulmonary disease (COPD) is increasing worldwide. COPD, in which tobacco smoke is the single biggest risk factor, accounts for about 400,000 deaths annually in industrialised countries [1]. In COPD, loss of lung tissue due to emphysema is a major component. Measurement of lung density by computed tomography (CT) is a surrogate marker for emphysema and reflects both disease severity [2–4] and physiological impairment [5–6]. Cross-sectional studies of the influence of smoking on lung density by CT have reported inconsistent results, with some researchers finding no influence of smoking habit on lung density [7], while others report decreased lung density with increasing pack-yrs [8–9]. This is not surprising given the fact that smoking has opposing effects on lung density. On one hand, current smoking triggers inflammation [10] that increases lung density [11], and on the other hand, smoking results in emphysema and decreased lung density in a substantial minority of smokers.
There is no evidence from longitudinal studies that smoking cessation halts or slows down the progression of emphysema. In recent studies, we have shown that smoking cessation is associated with decrease in lung density [11–12]; nevertheless, the influence on lung density changes of the burden of smoking, continuing smoking and the number of cigarettes smoked is largely unknown. In this regard, lung cancer screening trials with annual CT of the lungs provides a unique opportunity to study the influence of various factors on lung density including smoking habits and the presence and severity of airflow obstruction (AFO). In Denmark, we conducted a lung cancer screening trial [13], in which half of the 4,104 participants were randomised to annual CT of the lungs for 5 yrs between 2005 and 2009. These and similar data from ongoing lung cancer screening trials provide invaluable insight into the pathophysiology and natural history of COPD.
The literature on quantitative CT of emphysema is extensive, consisting primarily of cross-sectional descriptive studies because of the challenge imposed by lung volume changes during repeated scans. Change in lung volume is the most important confounder of the measurement of emphysema by CT resulting in poor reproducibility. To overcome this obstacle, we have applied statistical modelling to standardise lung density to lung volume both derived from the same scan, and in this way reduce measurement variations [14]. Using this method, we have shown a significant progression of emphysema [15] with a clear time trend that is superior to lung function measures [16]. In the current study, we used data from the Danish Lung Cancer Screening Trial (DLCST) and applied statistical modelling to adjust for lung volume changes between scans to study the influence of age, sex, smoking habits and the presence and degree of AFO on the decline in CT lung density.
MATERIAL AND METHODS
Study population
The study population represents those participants randomised to annual low-dose CT in the DLCST [13], which is a 5-yr trial investigating the effect of screening on lung cancer mortality. Individuals volunteered for the trial in response to advertisements in local free newspapers. From October 2004 to March 2006, 4,104 participants were enrolled and randomised to either annual low-dose CT or a control group, who were not offered CT screening. Participants in the DLCST were males and females who were 50–70 yrs of age without lung cancer related symptoms. Other inclusion criteria were: a history of cigarette smoking of at least 20 pack-yrs; current or ex-smoker; ex-smokers had to have quit after the age of 50 yrs and less than 10 yrs ago; forced expiratory volume in 1 s (FEV1) of at least 30% predicted normal at baseline; and provision of a written informed consent. Exclusion criteria were: body weight >130 kg; previous treatment for lung cancer, breast cancer, malignant melanoma or hypernephroma; any other malignancies within the last 5 yrs; tuberculosis within the last 2 yrs; or any serious illness that would shorten life expectancy to <10 yrs. In addition, measurements from individuals who changed their smoking habit during the study were excluded from the current analysis from the date of change of smoking habit.
In the screening group, 2,052 current or ex-smokers aged 50–70 yrs were screened annually for 5 yrs (2005–2009) with low-dose CT. At annual screening rounds, smoking habits were recorded, the carbon monoxide level in exhaled breath was measured and spirometry was performed.
The DLCST is funded by a government grant and was approved by the regional ethics committee. The trial is performed in one institution: Gentofte Hospital in Copenhagen, Denmark. A detailed description of the study design and study population of the DLCST has been previously published [13].
Imaging
All CT scans were performed on the same multi-detector CT scanner with 16 detector rows (Philips Mx 8000; Philips Medical Systems, Eindhoven, the Netherlands). Scans were performed in the supine position at full inspiration with a caudocranial direction including the whole lungs. Participants were instructed by voice recordings to take a deep breath and to hold their breath during scanning (10–12 s), and then when the scan was completed to breathe normally again. A low-dose technique was applied (120 kV and 40 mAs) with the following acquisition parameters: field of view 40 cm; collimation 16×0.75 mm; rotation time 0.5 s; and pitch 1.5. Images were reconstructed with 3 mm slice thickness using a soft reconstruction algorithm (kernel A). The scanner was usually maintained, quality-controlled and calibrated daily for air and at regular intervals for water according to the manufacturer's recommendations.
Image analysis
All CT scans were analysed by the Image Group at the Dept of Computer Science, University of Copenhagen (Copenhagen, Denmark) using in-house developed software designed to segment the lung and calculate various densitometric parameters [17]. The lung segmentation algorithm starts by detecting the trachea in the top slice. A special region growing algorithm segments the trachea down to the main carina. This region is not included in the calculations. Subsequently, the lungs are segmented starting from the main bronchi using a competing region growing algorithm and a threshold of -400 HU for the interface between lung and extrapulmonary tissue. The whole analysis is fully automated.
In a previous publication [11], the validity of the segmentation process was assessed by two physicians independently checking the segmentation of a sample of 200 CT scans. The observers checked whether the right and left lungs were correctly segmented and whether air from surrounding structures such as the oesophagus or bowel were erroneously included. Problems were encountered in 4% of segmentations and were all minor, corresponding to a CT volume <1 mL, except for one case in which 40 mL of bowel air was included in lung segmentation. No error in the separation of left and right lung was encountered in the reviewed scans.
From the segmented area, the total lung volume (TLV) was calculated. A frequency distribution of voxel densities (histogram) was generated from which the 15th percentile density (PD15) was extracted as the density in g·L−1 at which 15% of the voxels have lower densities. Density values were converted into g·L−1 by adding 1,000 to the density value in HU (e.g. PD15 value of -920 HU equals 80 g·L−1).
Lung density is expressed as the volume-adjusted PD15 throughout this paper. Changes in lung volume have substantial influence on lung density, which more than doubles from full inspiration to full expiration; therefore, PD15 was adjusted for the change in inspiration level between scans by physiological modelling using a so called “sponge model” [18]. From a theoretical point of view, the lungs could be considered a sponge-like structure, in which a proportional decrease in lung volume would yield an equally proportional increase in density, as compression would be mass-preserving. As a result, the product of PD15 and TLV is constant and independent of the level of inspiration, and therefore we can adjust PD15 to the predicted total lung capacity (TLC) by multiplying PD15 by TLV and dividing by predicted TLC [14].
Lung function tests
Spirometry was performed annually on all participants according to recommendations by the European Respiratory Society using electronic spirometers (Spirotrac IV; Vitalograph, Buckingham, UK). Results are expressed in absolute values and as per cent of predicted normal values according to European reference equations [19].
Statistics
The influence of age, sex, smoking and AFO on PD15 was analysed in a linear mixed effects model with random intercept and random linear effect of age (slope over time). Former smoking males, with no AFO and 20–30 pack-yrs at entry to the screening study, were chosen as a reference group. A linear decline in PD15 over time was assumed and the time point corresponding to age 60 yrs was used for cross-sectional comparisons. Furthermore, interactions between age (the time variable) and sex, smoking and AFO were included in the model as indications of the influences of these variables on the change of PD15 over time (that is the slopes in fig. 1). Because an interaction between age and another covariate indicates an influence on slope, it is interesting from a clinical point of view, and we decided to include these interactions in our model even if they were not statistically significant. Several other interactions between various covariates were tested as well, but none of them reached statistical significance, and they were not included in the final model.
The influence of various factors on lung density by computed tomography. a) Sex, b) accumulated previous smoking (pack-yrs), c) current smoking status and d) presence and degree of airflow obstruction (AFO). Solid lines represent a reference group of former smoking males without AFO and with 20–30 pack-yrs at entry to the study. aPD15: adjusted 15th percentile density; GOLD: Global Initiative for Chronic Obstructive Lung Disease.
We used a random coefficients model with linear time/age effect, allowing the intercept and slope to vary between subjects. It takes the form Yij = αi + βitij + eij where αi and βi are the individual intercept and slope associated with subject i, tij is the age of subject i at the time of the jth CT scan, and Yij is the corresponding CT lung density (PD15) at this time. The eij are independent random errors associated with Yij and are normally distributed with mean 0. The effects of covariates on the intercept and slope are modelled via αi = α + δ11×X1i + δ12×X2i + … + δ1k×Xki + ai and βi = β + δ21×X1i + δ22×X2i + … + δ2k×Xki + bi where α and β are fixed population effects (intercept and slope), Xm are the m=1, … , k covariates of interest, δ1m and δ2m are the m=1, … , k coefficients associated with the covariates, and ai and bi are the random intercept and slope for subject i. The random effects ai and bi for the same subject may be dependent, but it is assumed that each of them is independently, normally distributed with mean 0 across subjects.
RESULTS
Patient characteristics
The characteristics of the study population of the DLCST are comprehensively described elsewhere [13]. The number of participants in the screening group in the DLCST was 2,052. Of those, 10 participants were excluded from the current analysis, because they withdrew after randomisation but before undergoing a CT scan. The baseline characteristics of the 2,042 participants included are shown in table 1. There were 1,142 (56%) males in the screening group. The mean±sd age was 57.9±4.8 yrs, and males were on average 0.9 yrs older than females (p<0.001). All participants were heavy current 1,549 (76%) or ex-smokers 493 (24%) with a mean smoking history of 36.4±13.4 pack-yrs with males having smoked on average 6 pack-yrs more than females (p<0.001). For current smoking, the median tobacco consumption on entry to the trial was 20 cigarettes per day (range 1–65) with males smoking slightly more than females. Throughout the 5-yr screening period 1,134 (56%) participants remained active smokers, 424 (21%) remained ex-smokers and 484 (24%) changed their smoking habit (403 (20%) restarted and 81 (4%) stopped). At baseline, AFO defined as FEV1/forced vital capacity <0.7 was detected in 892 (44%) participants, of those 574 (64%) were in Global Initiative for Chronic Obstructive Lung Disease (GOLD) stage I; 291 (33%) in GOLD stage II; and 27 (3%) in GOLD stage III. Only one participant did not have valid spirometry data. There were no participants in GOLD stage IV as they were excluded from the study because if the screen detected lung cancer they could not be offered surgery.
Influence of age and sex on CT lung density
A total of 8,063 CT scans were obtained and analysed. Throughout the paper, all subsequent comparisons are made to the reference group mentioned in the statistics section. PD15±se in the reference group was 72.6±0.8 g·L−1 at age 60 yrs with an annual decline of -0.33±0.09 g·L−1 (p<0.001). Female sex was predictive of a higher PD15 at age 60 yrs (17.3±0.6 g·L−1, p<0.001) and was associated with increased annual loss of lung density of -0.3±0.1 g·L−1 (p<0.001) (table 2 and fig. 1a). The sex difference at age 60 yrs from a similar mixed effects model using raw PD15 data without volume adjustment was 12.7±0.8 g·L−1 (p<0.001), showing that the sex difference was not introduced by volume adjustment.
Influence of smoking on CT lung density
Smoking had a significant impact on lung density at age 60 yrs and on the annual decline in lung density. Smoking is represented by two variables in the model, the accumulated previous smoking (pack-yrs) and current smoking status. The number of pack-yrs was a strong predictor of a higher PD15 at age 60 yrs. Compared with 20–30 pack-yrs, individuals with 30–40 pack-yrs had an increased PD15±se level of 2.5±0.7 g·L−1 (p<0.001); and this was even higher in individuals with >40 pack-yrs, 5.5±0.8 g·L−1 (p<0.001) (table 2). There was no significant interaction between pack-yrs and age, thus the number of pack-yrs was not predictive of a steeper decline in PD15 (table 2 and fig. 1b).
In current smokers, the number of cigarettes smoked daily was a strong predictor of a higher PD15 compared with the reference group of ex-smokers and was PD15±se 8.98±0.73 g·L−1 (p<0.001) higher in those smoking 1–10 cigarettes per day; 10.33±0.71 g·L−1 (p<0.001) higher in those smoking 11–20 cigarettes per day; and 12.63±0.73 g·L−1 (p<0.001) higher in those smoking >20 cigarettes per day (table 2).
Continuous smoking during the study was also a strong predictor of accelerated annual decline in PD15. Compared with the reference group of former smoking males, participants smoking 1–10 cigarettes daily had an additional decline in PD15±se of -0.28±0.09 g·L−1 (p<0.01); those smoking 11–20 cigarettes daily had an additional decline of -0.48±0.08 g·L−1 (p<0.001); and those smoking >20 cigarettes a day had an additional decline of -0.42±0.09 g·L−1 (p<0.001) (table 2 and fig. 1c).
Influence of AFO on CT lung density
Severity of AFO according to GOLD stages was a strong predictor of lower PD15 at age 60 yrs: GOLD stage I PD15±se -1.42±0.70 g·L−1 (p=0.042); GOLD stage II -6.30±0.92 g·L−1 (p<0.001); and GOLD stage III -16.99±2.73 g·L−1 (p<0.001) (table 2). Likewise, severity of airway obstruction at baseline was a strong predictor of the magnitude of decline in PD15: GOLD stage I -0.21±0.08 g·L−1 (p=0.01); GOLD stage II -0.48±0.11 g·L−1 (p<0.001); and GOLD stage III -0.50±0.30 g·L−1 (p=0.10) (table 2 and fig. 1d). There was no significant interaction between pack-yrs and the severity of AFO (GOLD stages).
DISCUSSION
Analysis of this large cohort from the DLCST shows that age, sex, smoking habits and the presence of AFO have a significant impact on lung density. The presence of AFO is in addition associated with the largest annual loss in lung density. Furthermore, female sex and the number of cigarettes smoked daily throughout the trial are associated with an increased annual loss in lung density. Given the design of the trial, we have no data on lung density changes in healthy never-smokers; therefore, the reference group comprises ex-smokers without AFO. The highly significant time trend in the reference group probably indicates an age-related decline in lung density, a finding supported by earlier cross-sectional pathological and radiological studies [20, 21].
Based on 8,063 CT scans from 2,042 participants we have constructed a statistical model describing the influence of different factors including smoking habits on lung density. The robustness of the model is indicated by highly significant coefficients for the explanatory variables (most p<0.001), and the model is also fairly comprehensive including all significant interaction terms. Smoking is represented by two variables in the model, that is the accumulated previous smoking (pack-yrs) and current smoking status (0 cigarettes per day, i.e. ex-smoker, 1–10 cigarettes per day, 11–20 cigarettes per day and >20 cigarettes per day), and both these variables were strongly correlated to lung density at age 60 yrs. Three additional covariates were included in the model: age, sex and GOLD stage. The clinical interpretation of the model is somewhat simplified by the fact that only interaction terms involving age were statistically significant as can be seen in the different slopes in figure 1. Possible clinical interpretations of the model for each variable are discussed next.
Lung density in females was significantly higher than in males, and there was a significant interaction with age that is the mean annual loss of lung tissue was higher for females than for males (table 2 and fig. 1a). In patients with COPD, a number of studies reported more severe emphysema in males [22–24] and more rapid progression of emphysema in females [9, 25]. It is generally thought that biological differences make females more susceptible to the harmful effects of smoking than males [26]. In a recent study, Sverzillati et al. [25] reported certain morphological differences in emphysema suggesting a sex-related phenotype characterised by smaller emphysema lesions less prominent in the lung core. In order to rule out the possibility that the sex-related difference in PD15 was an artefact introduced by the volume adjustment, the analysis was repeated on raw PD15 data. The difference was 13 g·L−1, which indeed reflects a large difference between males and females that might be attributed to lower lung volumes in females. Independent of the burden of smoking and the presence of COPD, females have a significantly larger decline in PD15 than males; however, this difference does not seem to influence prognosis as females with COPD have better survival rates than males [27].
While the influence of sex on lung density seems to be straightforward, the influence of smoking is a bit more complex. There seem to be two pathological influences, one of inflammation and some degree of scarring (fibrosis) causing increased lung density, and in a substantial minority of individuals one of emphysema and reduced lung density. For the amount of previous smoking (i.e. pack-yrs) the model indicated a proportional increase in lung density (fig. 1b). The increase per pack-yr was modest (0.2–0.3 g·L−1), but highly significant (p<0.001), and no significant interaction was observed with the other explanatory variables including age. This increase in lung density may be counter-intuitive given the fact that smoking causes emphysema that implies loss of lung tissue. As will be discussed in more details later, low lung density at age 60 yrs was only observed in the presence of obstructive lung disease (GOLD stages) in our model, and this was not related to the previously smoked pack-yrs (i.e. the interaction between pack-yrs and GOLD stages was not significant). People without AFO who had stopped smoking prior to inclusion demonstrated increasing lung density for each pack-yr they had smoked. The limitation of our study was that we did not include subjects with <20 pack-yrs, and therefore, we have no information on lung density for such people. However, based on the fairly linear relationship between pack-yrs and lung density observed in our study population, it is tempting to assume that this relationship continues below 20 pack-yrs indicating that never smokers have less dense lungs than ex-smokers without COPD. Other investigators found that the risk of emphysema progression is higher in current than in ex-smokers; but nevertheless in line with our findings they reported no significant impact of pack-yrs on emphysema progression [28].
A possible explanation for this interesting finding could be the microscopic damage to lung tissue due to the noxious effect of smoking with inflammation and microscopic scarring (fibrosis). There is clear evidence that smoking induces an inflammatory reaction in the lung, which is evident even in young smokers [29]. Inflammation is associated with or followed by a certain degree of microscopic fibrosis. Auerbach et al. [30] in a microscopic study of over 1,800 autopsy lungs reported that the degree of fibrosis increased with increasing amounts of cigarette smoking. More recently, researchers have shown pathological and radiological evidence of interstitial fibrosis in smokers even in the absence of clinical evidence of interstitial lung disease [31, 32]. Another possible explanation for this proportional increase in lung density with the number of pack-yrs regardless of the presence of COPD is irreversible self-perpetuating inflammatory changes in the lung in smokers.
Current smoking increased lung density at age 60 yrs and at the same time interacted with age indicating an aggravated annual loss of lung density (fig. 1c). This finding is supported by Bellomi et al. [28], who found that the percentage of emphysema over a 2-yr period was significantly higher in current than in ex-smokers. The coefficients of the three smoking categories (1–10 cigarettes per day, 11–20 cigarettes per day and >20 cigarettes per day) showed a fairly clear dose–response relationship between the number of cigarettes smoked daily and both lung density at age 60 yrs and the annual loss of lung density. These opposing effects of smoking are not surprising, and are probably explained by smoking induced inflammation with increased movement of plasma and inflammatory cells from the circulation into the lung interstitium in the short-term and in the long-term remodelling and loss of lung tissue. The opposing effects of smoking on lung density may also explain the inconsistent results of previous cross-sectional studies [7–9]. In order to avoid large inter-individual variations, CT-measurements of participants who changed their smoking habits were excluded from the date of change. This might have introduced a selection bias, if individuals with rapidly declining lung density are more likely to quit smoking.
It is interesting that less than half of our participants (44%) had AFO. The presence of airways obstruction (GOLD stage 1–3) was the only explanatory variable that was related to low lung density at age 60 yrs and, independent of smoking, airway obstruction implied an exaggerated annual loss of lung density (fig. 1d). For both lung density at age 60 yrs and the annual loss of lung density, the relation had the character of a dose–response relationship. This is interesting and the clinical interpretation is obvious given that parenchymal destruction with loss of lung tissue (emphysema) is an important pathological mechanism in the development of COPD. There was no significant interaction between GOLD stages and prior or current smoking, which might suggest that this pathological process continued even after patients had stopped smoking, while continued smoking further aggravated the annual loss of lung density in patients with COPD.
Our study had some limitations. Partially excluding data from patients who changed their smoking habit may have introduced a selection bias because they might have stopped smoking because of rapid decline in lung density. However, smoking cessation induces sudden decline in lung density [12]. Due to the instructions during scanning, some participants may have performed a Valsalva manoeuvre, which affects lung perfusion by increasing intrathoracic pressure and decreasing venous blood flow. Variations in lung perfusion might have a potential influence on the reported results [33]. Finally, it would have been interesting to correlate the results on density with respiratory symptoms, but these data are currently not available.
In conclusion, females had a higher lung density and significantly larger annual decline than males. Both previous smoking (pack-yrs) and current smoking, independent of each other, related to an increase in lung density at age 60 yrs, although current smoking also seemed to increase the annual loss of lung density. Reduced lung density at age 60 yrs was only observed in individuals with AFO (GOLD stage 1–3). Furthermore, these individuals with AFO showed an exaggerated annual loss of lung density that, although more tempered, continued after they had stopped smoking.
Footnotes
Support Statement
The Danish Lung Cancer Screening Trial received financial support from the Danish Ministry of Interior and Health. The densitometric measurements reported in this paper were supported by unrestricted grants from the Danish Council for Strategic Research under the Programme Commission for Nanoscience and Technology, Biotechnology and IT (NABIIT) and from AstraZeneca R&D, Sweden.
Statement of Interest
Statements of interest for A. Dirksen, M. de Bruijne and for the study itself can be found at www.erj.ersjournals.com/site/misc/statements.xhtml
- Received November 27, 2011.
- Accepted February 10, 2012.
- ©ERS 2012