- Open Access
Diagnosis of chronic obstructive pulmonary disease in lung cancer screening Computed Tomography scans: independent contribution of emphysema, air trapping and bronchial wall thickening
- Onno M Mets1,
- Michael Schmidt2, 3,
- Constantinus F Buckens4,
- Martijn J Gondrie1,
- Ivana Isgum5,
- Matthijs Oudkerk6,
- Rozemarijn Vliegenthart6,
- Harry J de Koning7,
- Carlijn M van der Aalst7,
- Mathias Prokop8,
- Jan-Willem J Lammers9,
- Pieter Zanen9,
- Firdaus A Mohamed Hoesein1,
- Willem PThM Mali1,
- Bram van Ginneken3, 5,
- Eva M van Rikxoort†3 and
- Pim A de Jong†1Email author
© Mets et al.; licensee BioMed Central Ltd. 2013
Received: 1 March 2013
Accepted: 24 May 2013
Published: 27 May 2013
Beyond lung cancer, screening CT contains additional information on other smoking related diseases (e.g. chronic obstructive pulmonary disease, COPD). Since pulmonary function testing is not regularly incorporated in lung cancer screening, imaging biomarkers for COPD are likely to provide important surrogate measures for disease evaluation. Therefore, this study aims to determine the independent diagnostic value of CT emphysema, CT air trapping and CT bronchial wall thickness for COPD in low-dose screening CT scans.
Prebronchodilator spirometry and volumetric inspiratory and expiratory chest CT were obtained on the same day in 1140 male lung cancer screening participants. Emphysema, air trapping and bronchial wall thickness were automatically quantified in the CT scans. Logistic regression analysis was performed to derivate a model to diagnose COPD. The model was internally validated using bootstrapping techniques.
Each of the three CT biomarkers independently contributed diagnostic value for COPD, additional to age, body mass index, smoking history and smoking status. The diagnostic model that included all three CT biomarkers had a sensitivity and specificity of 73.2% and 88.%, respectively. The positive and negative predictive value were 80.2% and 84.2%, respectively. Of all participants, 82.8% was assigned the correct status. The C-statistic was 0.87, and the Net Reclassification Index compared to a model without any CT biomarkers was 44.4%. However, the added value of the expiratory CT data was limited, with an increase in Net Reclassification Index of 4.5% compared to a model with only inspiratory CT data.
Quantitatively assessed CT emphysema, air trapping and bronchial wall thickness each contain independent diagnostic information for COPD, and these imaging biomarkers might prove useful in the absence of lung function testing and may influence lung cancer screening strategy. Inspiratory CT biomarkers alone may be sufficient to identify patients with COPD in lung cancer screening setting.
Computed Tomography (CT)-based lung cancer screening has gained much interest after the National Lung Screening Trial (NLST) reported a 20% lung cancer mortality reduction in the CT arm compared to the chest radiography arm of the screening trial . Although implementation and cost-effectiveness is currently debated, screening has already started after the release of guidelines by two major organizations in the U.S. [2, 3]. Lung cancer screening CT scans enable the evaluation of other smoking-induced diseases besides lung cancer, such as chronic obstructive pulmonary disease (COPD), which accounts for significant morbidity and mortality and is also associated with an increased risk for lung cancer. Additional diagnosis of COPD might potentially be useful in the optimization of benefits and cost-effectiveness of CT-based screening in heavy smokers, and may also aid in the identification of the optimal target population for lung cancer screening  and personalization of lung cancer screening intervals. However, the opportunity of additional disease evaluation is not yet widely acknowledged .
Since lung function testing is not regularly incorporated in lung cancer screening, quantitative imaging biomarkers for COPD are likely to provide important surrogate measurements for disease evaluation. Moreover, computerized quantitative analysis is diagnostically superior to visual scores by human observers . It has recently been shown that automated quantitative analysis of screening CT images may be used to additionally identify COPD among screening participants . In that study, CT emphysema and CT air trapping, together with age, body mass index, smoking status and packyears smoked, provided good accuracy in the identification of COPD. To quantify the extent of air trapping, expiratory CT was added to the lung cancer screening protocol .
Software technology now allows for fully automatic analysis of bronchial wall thickness in inspiratory CT [8, 9]. Bronchial dimensions of the larger airways on CT have been shown to correlate well with those of the small airways measured by pathology , which may suggest that both reflect small airways disease. Another option is that large airway wall thickening is a separate marker of airway disease and visualizes the chronic bronchitis component of COPD [11, 12]. Therefore, quantification of bronchial wall thickness may be able to either replace the CT air trapping measurements (and hence the need for additional expiratory imaging), or may add independent diagnostic value for COPD and therewith improve the diagnostic performance of CT. The objective of this study is to determine the independent diagnostic value of CT emphysema, CT air trapping and CT bronchial wall thickness in the identification of COPD in low-dose lung cancer screening CT.
Materials and methods
This study was performed as a side-study of the Dutch and Belgian lung cancer screening trial (NELSON-trial, ISRCTN 63545820) . The study was approved by both the Dutch Ministry of Health and the local ethical review board of the University Medical Center Utrecht. Each participant provided written informed consent.
The selection of the study population has been described in detail elsewhere . In summary, we included 1140 male subjects with a paired inspiratory and expiratory CT obtained on the same day between July 2007 and September 2008. The included study population was a subsample of the total lung cancer screening population. The subsample did not differ significantly from the total population of subjects screened for lung cancer, as previously reported . Participants included in the lung cancer screening trial were current and former (<10 year) heavy smokers with a smoking history of at least 16.5 packyears, and were at baseline between 50 and 75 years of age. All participants weighted less than 140 kg and were physically fit enough to undergo surgery. Participants provided self-reported presence of respiratory symptoms. Participants were asked to indicate whether cough, sputum production, dyspnea and wheezing were present for at least three months a year. A symptomatic status was assigned when presence of at least 1 of the symptoms was reported. Otherwise, a participant was assigned an asymptomatic status.
Pulmonary function testing
Prebronchodilator spirometry was performed using ZAN equipment (ZAN Messgeräte GmbH, Germany), according to European Respiratory Society and American Thoracic Society guidelines [14, 15]. Forced expiratory volume in one second (FEV1) was expressed as percentage predicted. Forced vital capacity was expressed as percentage. We defined COPD as FEV1/FVC below 70%, with a classification into mild (FEV1 ≥ 80%), moderate (80% > FEV1 ≥ 50%) and severe (50% > FEV1 ≥ 30%).
All subjects received paired inspiratory and expiratory CT, obtained with 16 × 0.75 mm collimation (Brilliance 16P; Philips Medical Systems). Data was acquired according to the low-dose protocol of the lung cancer screening trial; 120 kVp (≤80 kg) or 140 kVp (>80 kg) both at 30 mAs for the inspiratory CT, and 90 kVP (≤80 kg) or 120 kVp (>80 kg) both at 20 mAs for the expiratory CT. The estimated effective dose was at maximum 1.4 millisievert (mSv) for the inspiratory acquisition, and at maximum 0.65 mSv for the expiratory acquisition. CT images were reconstructed from lung bases to apices with slice thickness of 1 mm at 0.7 mm increments, using a smooth reconstruction algorithm (B-filter, Philips).
Quantitative analysis of emphysema and air trapping
Quantitative analysis of bronchial dimensions
To check the bronchial segmentation process in each case, the 3D bronchial tree segmentation, a random selection of about 90 bronchial cross-sections and the plotted regression line were presented to an independent observer. After carefully reviewing the segmentation results cases with major errors in the bronchial tree segmentation (N = 22) or where segmentation failed (N = 6) were excluded. No manual corrections were performed.
To obtain normal distribution, transformation using the natural logarithm was applied on the CT emphysema measure (logIN−950). The other variables did not have to be transformed. Missing values were imputed using multiple imputation technique; quantitative bronchial wall measurements were missing in 28 cases and respiratory symptoms were missing in 55 cases. There were no other variables with missing values.
Various multivariate models were developed using logistic regression analysis with COPD as outcome variable. The modeling started with a baseline model including only demographic variables (i.e. age, body mass index, smoking history and smoking status) as associated variables. Subsequently, quantitative CT emphysema, CT air trapping, and CT bronchial wall thickness were added to the analysis in a stepwise fashion. This resulted in a total of eight different models: 1) demographic variables only; 2) demographic variables + CT emphysema; 3) demographic variables + CT air trapping; 4) demographic variables + CT bronchial wall thickness; 5) demographic variables + CT bronchial wall thickness + CT air trapping; 6) demographic variables + CT emphysema + CT air trapping; 7) demographic variables + CT emphysema + CT bronchial wall thickness; 8) demographic variables + CT emphysema + CT bronchial wall thickness + CT air trapping.
Internal validation was performed for each model by assessing the overoptimism, using bootstrap resampling with 500 iterations [26, 27], and subsequent shrinkage of the initial coefficients by the estimated degree of overoptimism . Calibration (i.e. the agreement between the predicted and observed values) of the shrunk models was assessed visually using the calibration plots. Discrimination (i.e. the ability of a model to differentiate between subjects with and without COPD) was assessed by calculation of the area under the receiver operating characteristic (ROC) curve (i.e. C-statistic). The point of optimal accuracy was defined as the highest number of true positive plus true negative cases. For this cut-off value, sensitivity and specificity was calculated for each model. Using the Net Reclassification Index (NRI) , all eight models were compared to each other to determine the added value of each quantitative marker. The NRI assesses the added usefulness of a marker by showing the increase in correctly categorized cases (i.e. the shift from false negative to true positive and from false positive to true negative). In addition to the primary analysis in the total study population, we also applied each model separately in the subgroups of symptomatic and asymptomatic subjects.
A p-value lower than 0.05 was considered to represent a statistically significant difference. Analysis were performed using R-statistical program (v.2.13.1) using regression modeling strategies v.3.3-0 and ROCR v.1.0-4 . Values are presented as mean ± SD, unless indicated otherwise.
Study population characteristics
Mean ± SD or Number (%)
N = 1140
62.5 ± 5.2
Body mass index, kg/m2
27.1 ± 3.6
38 (28 – 49)
Mild disease ‡
Moderate disease §
Severe disease ll
Quantitative CT analysis
Segmented lung volumes were 6.79 ± 1.13 liter in inspiratory CT and 3.71 ± 0.90 liter in expiratory CT. The quantified amount of CT emphysema showed a median (P25 – P75) value of 0.75 (0.40 – 1.46)%. For CT air trapping this was 0.84 (0.80 – 0.88). On average 1280 ± 403 (range 260 – 3536) bronchial cross-sections throughout the lung were obtained per subject. The Pi10 for all participants was 2.41 ± 0.51 mm.
Diagnostic performance of the multivariate models
Performance measures for the various multivariate models to identify COPD
Baseline model *
+ CT air trapping
+ CT emphysema
+ CT air trapping
+ CT emphysema
+ CT air trapping
+ CT emphysema
+ CT emphysema
+ CT air trapping
When a combination of two CT biomarkers was added to the baseline model, diagnostic performance improved even further. Combinations of emphysema and air trapping, emphysema and bronchial wall thickening, or air trapping and bronchial wall thickening showed a C-statistic of 0.834, 0.861 and 0.801, respectively. The NRI compared to the baseline model was 32.4%, 39.9% and 24.8%, respectively.
Discrimination of the various multivariate models in the identification of COPD in subgroups of asymptomatic and symptomatic subjects
C-index (95% CI)
Baseline model *
0.674 (0.625 -0.722)
0.634 (0.589 - 0.679)
0.739 (0.695 - 0.783)
0.764 (0.725 - 0.803)
+ CT air trapping
0.737 (0.693 - 0.780)
0.794 (0.759 - 0.829)
+ CT emphysema
0.753 (0.707 – 0.800)
0.806 (0.771 - 0.841)
0.771 (0.730 – 0.813)
0.821 (0.788 – 0.855)
+ CT air trapping
+ CT emphysema
0.782 (0.740 - 0.824)
0.872 (0.844 - 0.899)
+ CT air trapping
+ CT emphysema
0.828 (0.790 – 0.866)
0.886 (0.859 - 0.912)
+ CT emphysema
0.832 (0.795 – 0.869)
0.905 (0.881 - 0.929)
+ CT air trapping
This study shows that CT emphysema, CT air trapping and CT bronchial wall thickness, automatically assessed on low-dose lung cancer screening CT images, all have independent diagnostic value for COPD. The best performance measures were shown by the model that included all three quantitative CT biomarkers. This full model yielded a very good diagnostic performance and assigned the correct status the vast majority of the screening participants. Using the full model almost three-quarter of the COPD subjects in our population were automatically identified, with a low false-positive rate. However, the additional diagnostic value of expiratory CT was limited, and performance measures remained sufficient when expiratory CT was not included.
Our results lead to several considerations. Firstly, although the concept of additional COPD detection has yet received limited attention and early diagnosis of COPD is subject to debate, the concept of identifying COPD as early as possible using screening CT is feasible. The actual extent of benefits from early diagnosis and intensified smoking cessation-strategies is yet unknown, however, it has been shown that early smoking cessation lowers disease progression and decreases COPD-related morbidity and mortality [31, 32]. Moreover, in the setting of lung cancer screening there is also a second rationale to diagnose COPD; having COPD or emphysema increases the risk of developing lung cancer and therefore may aid in risk assessment and personalization of the screening strategy . Since lung function testing is generally not included in lung cancer screening, other measurements might be needed to allow for COPD evaluation. Imaging biomarkers as presented in this study may thus prove important surrogate measures in lung cancer screening setting. Notwithstanding, it should be emphasized that the approach of using automated identification of COPD in screening CT cannot yet be widely implemented at this moment as it needs to be externally validated and quantitative CT analyses need further standardization to allow inter-technique and inter-site comparison .
Secondly, our results show that inspiratory CT bronchial wall thickness might usefully replace the expiratory CT air trapping measurements. This is specifically interesting given that expiratory acquisitions are generally not performed in other screening trials, and might not be performed when lung cancer screening is actually implemented. The present study shows that the sensitivity and specificity of the diagnostic model remained acceptable when expiratory air trapping was not included into the model, implying that inspiratory CT biomarkers alone are valuable in identifying patients with COPD in early stages in routine lung cancer screening setting. Although addition of the expiratory CT data did reclassify a further 4.5% of subjects, this effect is limited for diagnostic purposes and may not justify obtaining an expiratory CT scan in each screening participant.
Our study has several potential limitations. Firstly, the generalizability of our study may be limited given that our study population comprises only males. Few women are included in the lung cancer screening trial, since there are few women in the Dutch population with accumulated long-term smoking histories . Also, differences in CT protocols (e.g. CT equipment, dose settings, reconstruction algorithm) may influence absolute values obtained from quantitative image analysis. Although we believe it is highly unlikely that such differences could lead to contradicting findings when applied in another study, close attention should always be paid to technical parameters when comparing quantitative results. Secondly, while the official definition of COPD is based on post-bronchodilator spirometry, spirometry in the lung cancer screening trial was obtained without administration of a bronchodilator due to time restrictions. This may lead to some misclassification due to asthma. However, we do not think this is an important confounder given that the prevalence of asthma in a population of heavy smokers is generally very low. In the Netherlands asthma is present in about 2% of all males between 50 and 75 years of age .
In the absence of lung function testing quantitative CT biomarkers provide useful surrogate measures for COPD evaluation, which may in the future influence screening strategy and enhance screening benefits and cost-effectiveness. The present study shows that quantitatively assessed CT emphysema, CT air trapping and CT bronchial wall thickness each contain independent diagnostic information for COPD, and represent independent imaging biomarkers. As the added diagnostic value of expiratory CT air trapping is limited, inspiratory CT biomarkers alone may be sufficient to identify subjects with COPD in lung cancer screening setting.
The NELSON-trial was supported by the Netherlands Organisation for Health Research and Development; Dutch Cancer Society Koningin Wilhelmina Fonds; Stichting Centraal Fonds Reserves van Voormalig Vrijwillige Ziekenfondsverzekeringen; Siemens Germany; Roche Diagnostics; Rotterdam Oncologic Thoracic Steering Committee; and the G.Ph. Verhagen Trust, Flemish League Against Cancer, Foundation Against Cancer, and Erasmus Trust Fund.
- Aberle DR, Adams AM, Berg CD, et al: Reduced lung-cancer mortality with low-dose computed tomographic screening. N Engl J Med. 2011, 365 (5): 395-409.PubMedView ArticleGoogle Scholar
- Wood DE, Eapen GA, Ettinger DS, et al: Lung cancer screening. J Natl Compr Canc Netw. 2012, 10 (2): 240-265.PubMedGoogle Scholar
- American Lung Association: Providing Guidance for Lung Cancer Screening: The American Lung Association Interim Report on Lung Cancer Screening. 2012, Available from: [http://www.lung.org/finding-cures/research-news/new-screening-guidelines/lung-cancer-screening.pdf]. Accessed February 2013Google Scholar
- Young RP, Hopkins RJ: Lung cancer risk prediction to select smokers for screening CT. Cancer Prev Res. 2012, 5 (4): 697-698. 10.1158/1940-6207.CAPR-11-0531.View ArticleGoogle Scholar
- Mets OM, de Jong PA, Prokop M: Computed tomographic screening for lung cancer: an opportunity to evaluate other diseases. JAMA. 2012, 308 (14): 1433-1434. 10.1001/jama.2012.12656.PubMedView ArticleGoogle Scholar
- Mets OM, Smit EJ, Mohamed Hoesein FA, et al: Visual versus automated evaluation of chest computed tomography for the presence of chronic obstructive pulmonary disease. PLoS One. 2012, 7 (7): e42227-10.1371/journal.pone.0042227.PubMedPubMed CentralView ArticleGoogle Scholar
- Mets OM, Buckens CF, Zanen P, et al: Identification of chronic obstructive pulmonary disease in lung cancer screening computed tomographic scans. JAMA. 2011, 306 (16): 1775-1781. 10.1001/jama.2011.1531.PubMedView ArticleGoogle Scholar
- van Ginneken B, Baggerman W, van Rikxoort EM: Robust segmentation and anatomical labeling of the airway tree from thoracic CT scans. Med Image Comput Comput Assist Interv. 2008, 11 (Pt 1): 219-226.PubMedGoogle Scholar
- Schmidt M, Kuhnigk JM, Krass S, Owsijewitsch M, Hoop BJ, Peitgen H: Reproducibility of airway wall thickness measurements. Proceedings of the SPIE, Medical Imaging. 2010, 7624: 76240-1-76240-10.View ArticleGoogle Scholar
- Nakano Y, Muro S, Sakai H, et al: Computed tomographic measurements of airway dimensions and emphysema in smokers. Correlation with lung function. Am J Respir Crit Care Med. 2000, 162 (3 Pt 1): 1102-1108.PubMedView ArticleGoogle Scholar
- Kim V, Han MK, Vance GB, et al: The chronic bronchitic phenotype of chronic obstructive pulmonary disease: An analysis of the COPDgene study. Chest. 2011, 140 (3): 626-633. 10.1378/chest.10-2948.PubMedPubMed CentralView ArticleGoogle Scholar
- Mair G, Maclay J, Miller JJ, et al: Airway dimensions in COPD: relationships with clinical variables. Respir Med. 2010, 104 (11): 1683-1690. 10.1016/j.rmed.2010.04.021.PubMedView ArticleGoogle Scholar
- van Iersel CA, de Koning HJ, Draisma G, et al: Risk-based selection from the general population in a screening trial: selection criteria, recruitment and power for the Dutch-Belgian randomised lung cancer multi-slice CT screening trial (NELSON). Int J Cancer. 2007, 120 (4): 868-874. 10.1002/ijc.22134.PubMedView ArticleGoogle Scholar
- Miller MR, Crapo R, Hankinson J, et al: General considerations for lung function testing. Eur Respir J. 2005, 26: 153-161. 10.1183/09031936.05.00034505.PubMedView ArticleGoogle Scholar
- Quanjer PH, Tammeling GJ, Cotes JE, Pedersen OF, Peslin R, Yernault JC: Lung volumes and forced ventilatory flows. Report Working Party Standardization of Lung Function Tests, European Community for Steel and Coal. Official Statement of the European Respiratory Society. Eur Respir J Suppl. 1993, 16: 5-40.PubMedView ArticleGoogle Scholar
- van Rikxoort EM, de Hoop B, Viergever MA, Prokop M, van Ginneken B: Automatic lung segmentation from thoracic computed tomography scans using a hybrid approach with error detection. Med Phys. 2009, 36: 2934-2947. 10.1118/1.3147146.PubMedView ArticleGoogle Scholar
- Schilham AM, van Ginneken B, Gietema H, Prokop M: Local noise weighted filtering for emphysema scoring of low-dose CT images. IEEE Trans Med Imaging. 2006, 25: 451-463.PubMedView ArticleGoogle Scholar
- Gevenois PA, de Maertelaer V, De Vuyst P, Zanen J, Yernault JC: Comparison of computed density and macroscopic morphometry in pulmonary emphysema. Am J Respir Crit Care Med. 1995, 152 (2): 653-657. 10.1164/ajrccm.152.2.7633722.PubMedView ArticleGoogle Scholar
- O'Donnell RA, Peebles C, Ward JA, et al: Relationship between peripheral airway dysfunction, airway obstruction, and neutrophilic inflammation in COPD. Thorax. 2004, 59 (10): 837-842. 10.1136/thx.2003.019349.PubMedPubMed CentralView ArticleGoogle Scholar
- Mets OM, Murphy K, Zanen P, et al: The Relationship between Lung Function Impairment and Quantitative Computed Tomography in Chronic Obstructive Pulmonary Disease. Eur Radiol. 2012, 22 (1): 120-128. 10.1007/s00330-011-2237-9.PubMedPubMed CentralView ArticleGoogle Scholar
- Mets OM, Zanen P, Lammers JW, et al: Early Identification of Small Airways Disease on Lung Cancer Screening CT: Comparison of Current Air Trapping Measures. Lung. 2012, 190 (6): 629-633. 10.1007/s00408-012-9422-8.PubMedView ArticleGoogle Scholar
- Mets OM, Isgum I, Mol CP, et al: Variation in quantitative CT air trapping in heavy smokers on repeat CT examinations. Eur Radiol. 2012, 22 (12): 2710-2717. 10.1007/s00330-012-2526-y.PubMedPubMed CentralView ArticleGoogle Scholar
- Mets OM, Willemink MJ, de Kort FP, et al: The effect of iterative reconstruction on computed tomography assessment of emphysema, air trapping and airway dimensions. Eur Radiol. 2012, 22 (10): 2103-2109. 10.1007/s00330-012-2489-z.PubMedPubMed CentralView ArticleGoogle Scholar
- Nakano Y, Wong JC, de Jong PA, et al: The prediction of small airway dimensions using computed tomography. Am J Respir Crit Care Med. 2005, 171 (2): 142-146. 10.1164/rccm.200407-874OC.PubMedView ArticleGoogle Scholar
- Patel BD, Coxson HO, Pillai SG, et al: Airway wall thickening and emphysema show independent familial aggregation in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2008, 178 (5): 500-505. 10.1164/rccm.200801-059OC.PubMedView ArticleGoogle Scholar
- Harrell F: Regression modeling strategies with Applications to Linear Models, Logistic Regression, and Survival Analysis. 2001, New York: SpringerGoogle Scholar
- Efron B, Tibshirani RJ: An Introduction to the Bootstrap. 1993, London: Chapman & HallView ArticleGoogle Scholar
- Steyerberg EW: Clinical Prediction Models: A practical approach to development, validation and updating. 2009, New York: SpringerView ArticleGoogle Scholar
- Pencina MJ, D'Agostino RB, D'Agostino RB, Vasan RS: Evaluating the added predictive ability of a new marker: from area under the ROC curve to reclassification and beyond. Stat Med. 2008, 27 (2): 157-172. 10.1002/sim.2929.PubMedView ArticleGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing, Available from: [http://www.R-project.org/]. Accessed February 2013
- Soriano JB, Zielinski J, Price D: Screening for and early detection of chronic obstructive pulmonary disease. Lancet. 2009, 374 (9691): 721-732. 10.1016/S0140-6736(09)61290-3.PubMedView ArticleGoogle Scholar
- Decramer M, Sibille Y, Bush A, et al: The European Union conference on chronic respiratory disease: purpose and conclusions. Eur Respir J. 2011, 37 (4): 738-742. 10.1183/09031936.00020211.PubMedView ArticleGoogle Scholar
- Mets OM, de Jong PA, van Ginneken B, Gietema HA, Lammers JWJ: Quantitative Computed Tomography in COPD: Possibilities and limitations. Lung. 2012, 190 (2): 133-145. 10.1007/s00408-011-9353-9.PubMedPubMed CentralView ArticleGoogle Scholar
- Nationaal Kompas Volksgezondheid [National Public Health Compass]: Available from: [http://www.nationaalkompas.nl/gezondheid-en-ziekte/ziekten-en-aandoeningen/ademhalingswegen/astma/omvang]. Accessed February 2013
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.