The Prevalence of Chronic Obstructive Pulmonary Disease (COPD) and the Heterogeneity of Risk Factors in the Canadian Population: Results from the Canadian Obstructive Lung Disease (COLD) Study.

PURPOSE: To determine the spirometric-based prevalence of COPD across different regions in Canada and to evaluate the site heterogeneity of risk factors. PATIENTS AND METHODS: In this cross-sectional, population-based study, random samples of non-institutionalized adults aged ≥40 years were generated by random digit dialling. Participants answered an interviewer-administered questionnaire and performed spirometry before and after bronchodilator administration. COPD was defined as post-bronchodilator FEV1/FVC <0.70 (fixed ratio, FR) and as FEV1/FVC <5th percentile (lower limits of normal, LLN). Separate logistic regression models were used to compute the risk (adjusted odds ratio, aOR) for COPD. I2 and Tau2 analyses were used to evaluate heterogeneity.
RESULTS: Out of 5176 (95%) participants, 4893 (47% male with mean age 56.6 years (95% confidence interval, 56.0-57.2)) had spirometry that satisfied ATS criteria. The population prevalence of COPD was 16.2% (95% CI, 14.5-17.8) by FR and 11.2% (95% CI, 9.7-12.6) by LLN. Male predominance in prevalence was shown by FR but not by LLN criteria. Patient characteristics associated with an increased risk of COPD included: age (OR 1.56; 95% CI 1.33-1.84); history of physician-diagnosed asthma (OR 3.30; 95% CI 2.42-4.49); and childhood hospitalization for respiratory illness (OR 1.81; 95% CI 1.17-2.80). In terms of smoking-related risk factors, current smoking status had the highest odds ratio (OR 3.49; 95% CI 2.55-4.80). Variance in prevalence among sites was significantly reduced by adjusting for risk factors in Tau2 analyses. Higher odds of exposure for each risk factor was found in more severe COPD, suggesting that a higher risk could be linked to the development of severe disease.
CONCLUSION: This study reports the population prevalence of COPD in nine urban cities which collectively represent the majority of the Canadian population and demonstrates that heterogeneity in prevalence among sites is substantially explained by variation in associated risk factors for COPD.

Introduction

The burden of COPD in Canada and worldwide is substantial and increasing.1 The Global Burden of Disease study estimated that the projection of global increase in COPD from 1990 to 2020 will rise from sixth to third in terms of morbidity.2 Chronic respiratory diseases collectively account for 4.7% of global disability-adjusted life year, with COPD comprising two-third of this total.3 In addition, COPD has a significant impact on healthcare resource utilization and workplace productivity loss. Previous results of Canadian economic analysis from 2003 demonstrated that the annual direct cost of COPD-related primary and secondary care visits, treatments, and laboratory tests was at almost $2000 CAD per patient, with unscheduled care visits accounting for almost 60% of the total direct cost of COPD.4

Spirometry is the most widely available, reproducible test of lung function and is essential to the diagnosis of COPD and reduces underdiagnosis.5 Population-based data suggested that the prevalence of COPD confirmed by spirometry is 2.5 times higher than the rate for self-reported diagnosis and fourfold higher than the previous estimates based on annual community surveys.6,7

Insight regarding the impact of variations in risk factors on COPD is important in interpreting prevalence of this disease. The relationship between tobacco smoking and the development of COPD is well established. The rate of daily or occasional smoking in Canada has decreased from 25.9% to 19.9% between 2001 and 2011,8 yet the prevalence of COPD continues to increase. A decline in smoking may not have an immediate and parallel impact on the prevalence in COPD. As COPD is the result of long-term cumulative exposure to noxious gas and particles such as in tobacco smoking, a time lag between a fall in the prevalence of risk factors and prevalence in COPD is to be expected. Other reasons for a rising COPD prevalence include improved diagnosis and an aging population as the incidence of COPD peaks in older adults.

There is little information on the heterogeneity of risk factors within a country and its impact on the variation in COPD prevalence across different regions.6 Understanding the impact of variations in risk factors on COPD is important in interpreting the prevalence of this disease. Our working hypothesis was that much of the apparent heterogeneity among large urban centers in the same country, with similar health-care delivery, might be explained by differences in demographics (eg, age) and other risk factors such as asthma. It is also unclear whether the association between risk factors and COPD varies with disease severity.

To address these issues, the Canadian Obstructive Lung Disease (COLD) Study, a pan-Canadian initiative was formed to measure COPD across the nation with similar sampling strategy and standardized methods as that of the International Burden of Obstructive Lung Disease (BOLD)9,10 and the PLATINO initiatives11 with the aim to provide accurate population-based estimates of COPD prevalence and risk factors. The study recruited adults aged ≥40 years through random digit dialling from nine Canadian urban city sites and collected demographic and clinical information, including spirometry before and after bronchodilation, between 2006 and 2011.6

The objective of this cross-sectional multisite, nationwide, population-based study was to determine the estimated population prevalence of COPD as defined by two widely used spirometric definitions (FR and LLN), in nine Canadian cities; to evaluate the heterogeneity of risk factors associated with the disease across sites; and to assess whether the variation in prevalence for COPD between site can be reduced by adjustment for common risk factors.

Methods

Study Population and Design

Data for this cross-sectional, multisite, nationwide, population-based study were collected between 2006 and 2011 in 9 sites across Canada. The overall sampling strategy and baseline study protocol of the COLD study were the same as that used in the BOLD initiative, the full details of which have been published previously.6,9,10 Random samples of non-institutionalized adults aged ≥40 years in nine Canadian urban sites in six provinces (Vancouver, Montreal, Toronto, Halifax, Calgary, Ottawa, Kingston, Quebec City, Saskatoon) were drawn as in the BOLD protocol,9 with each sampling frame consisting of a population of at least 150,000 individuals of all ages. Using landline telephone digit dialling to identify eligible individuals, random samples were created by BCStats (Victoria, BC, Canada) using official data from Statistics Canada (Survey and Analysis Section, Victoria, BC, Canada) and recruitment was conducted by NRG Research Group (Vancouver, Canada) at each respective site. Eligible individuals were invited to attend a visit to complete an interviewer-administered respiratory questionnaire and to perform pre- and post-bronchodilator spirometry by trained research staff.

Study Questionnaire and Spirometric Measurements

Trained and certified technicians administered questionnaires in English or French to participants. The standardized questionnaire included demographic details, information on respiratory symptoms and diseases, smoking history (tobacco and marijuana) and other risk factors for COPD, medication use, health-care utilization, activity limitation, health-related quality of life, and occupational and biomass exposure. Lung function data were obtained in the seated position with the EasyOne spirometer (ndd Medical Technologies Inc., Andover, MA, USA) before and 15 minutes after administration of two puffs (200 µg) of salbutamol via a metered-dose inhaler with a spacer. Pulmonary function quality assurance with over reading was conducted, in which all spirograms were reviewed and graded using American Thoracic Society (ATS) and European Respiratory Society (ERS) guidelines with prompt feedback to the technician at each site.12 All sites used the same scoring criteria and protocol. The quality of all test sessions was reviewed centrally in the COLD study pulmonary function reference centre in Vancouver. Only spirometric data that fulfilled the ATS/ERS acceptability and repeatability criteria were used for analysis (at least two acceptable and reproducible tests for both forced expiratory volume in 1 second [FEV1] and forced vital capacity [FVC]).13

Definitions of Airflow Limitation

Chronic airflow limitation was defined, a priori, as post-bronchodilator FEV1/FVC <0.70 or fixed ratio (FR). The severity of airflow limitation was categorized according to the Global Initiative for Chronic Obstructive Lung Disease (GOLD) stages, grades 1–4.14 Our study also computed the prevalence of COPD as defined by the lower limits of normal (LLN) of FEV1/FVC, an alternative spirometric definition for airflow limitation.9,11,14–16

Consent

All participants gave written informed consent, and the study was approved by all respective University and Institutional Ethical Review Boards: UBC/PHC Research Ethics Board, P05-006 (Vancouver); Biomedical-C Research Ethics Board, BMC-06-002 (Montreal); UHN REB, 06–0421-B (Toronto); Capital Health Research Ethics Board, CDHA-RS/2007-255 (Halifax); Conjoint Health Research Ethics Board, ID21258 (Calgary); DMED-1240-09 (Kingston); 2009519–01H (Ottawa); Bio-REB09-162 (Saskatoon); CER20459 (Quebec City).

Statistical Analysis

All data analyses were performed using statistical software (SAS version 9.4, SAS Institute, Cary, NC, USA), and all tests were two-tailed in nature; p ≤ 0.05 was considered significant. Univariate analysis of variables across sites was performed using one-way analysis of variance for continuous variables or the Χ2-test for dichotomous variables.

The prevalence data in this study were weighted to the sex and age profile of the Canadian population. Estimated population prevalence was shown as % (95% CI). Sample prevalence was converted to site-specific population prevalence by adjusting for age and sex (weighted) to that of the local site population. The details for each site population were obtained from the Canadian population (Census 2006). The population estimates were shown by site and sex, and COPD was defined spirometrically in two ways: 1) by the Global Initiative Obstructive Lung Disease GOLD definition,14 that is FEV1/FVC<0.7 (as well as different cut-offs of lung function: FEV1/FVC <0.7 and FEV1% predicted ≥80% (GOLD grade 1); FEV1/FVC <0.7 and FEV1% predicted < 80% (GOLD grades 2 to 4); and 2) by less than the lower limit of normal (LLN) of FEV1/FVC (FEV1/FVC <5th percentile [predicted FEV1/FVC minus 1.645 × standard error of estimate, SEE]).17,18

Separate logistic regression models were constructed to evaluate risk factors independently for COPD which included age, sex, years of education, smoking status, history of physician-diagnosed asthma, and childhood hospitalizations for respiratory illnesses and exposure to dusty jobs. A separate random effects meta-analysis model was constructed for each risk factor to estimate odds ratios and to assess heterogeneity by sites and by sex with use of the I2 measures.19 Individual sites were weighted according to the precision of the estimate at that site and these weights were reflected in the forest plots. The odds ratios shown for each risk factor were adjusted for age and school years, with additional adjustment as follows: 1) for asthma: pack-years, childhood hospitalization, and dusty job exposure; 2) for childhood hospitalization: pack-years and dusty jobs; 3) for pack-years: asthma, childhood hospitalization, and dusty job exposure; 4) for ever smoking: asthma, childhood hospitalization, and dusty job exposure; and 5) for current smoking status: asthma, childhood hospitalization, and dusty job exposure. Further details of weightage are provided in the . The forest plot analyses used the metan procedure in Stata V.10.0 (Stata Corp).

To show whether the variation in prevalence for COPD between site is substantially reduced by adjustment for risk factors in this study, we conducted an additional random-effects meta-analysis analysis19 to assess the change in the variation of the prevalence for LLN COPD among the nine sites (heterogeneity) by adjustment for factors such as age, smoking, asthma and childhood hospitalization using two logistic regressions with random intercept (Glimmix procedure in SAS). The first model did not include any risk factors; the second multiple model included the set of risk factors; and the mean difference for the tau-squared (Tau2) and 95% CI was estimated using 100 bootstrap samples. Additional details of this method are in the .

Results

Participant Demographics and Risk Factors

Table 1 shows the characteristics of the participants from the nine sites recruited between 2006 and 2011. A total of 5176 participants completed the study, comprising 47% men with mean age of 56.6 years (Table 1), of which 4893 (95%) individuals had spirometry results that satisfied ATS/ERS criteria. The mean study response rate was 78% (ranged 63% to 88% across the 9 sites, ). The proportion of individuals aged 40 years and older in the Canadian population increased slightly from 49.34% to 51.60% between the 2006 and 2016 Canadian Census (). The proportions of non-Hispanic white and non-white ethnic groups in this study were 88.5% and 11.5%, respectively ().

Table 1

Characteristics of the Participants in the Nine Sites of the Canadian Obstructive Lung Disease Study

Variables	All Sites	Vancouver	Calgary	Saskatoon	Toronto	Kingston	Ottawa	Montreal	Quebec City	Halifax	Overall P-value
Total, n(%)	5176	856	600	507	548	550	483	557	595	480
Male sex, n(%)	2292(47.2)	362(48.2)	230(48.3)	262(47.5)	236(47.3)	254(47.0)	243(47.1)	237(46.9)	290(46.2)	168(46.3)	<0.0001
Age, years	56.6(0.3)	57.0(0.5)	55.2(0.5)	57.4(0.6)	56.4(0.6)	58.2(0.6)	56.4(0.6)	57.0(0.7)	57.6(0.6)	56.5(0.6)	<0.0001
Race ethnicity, n(%)											<0.0001
Non-Hispanic white	4741(88.5)	666(77.8)	565(93.4)	496(97.6)	449(80.6)	533(96.8)	465(96.5)	509(91.8)	590(99.1)	468(98.0)
Black	60(2.8)	3(0.3)	4(0.8)	1(0.2)	25(4.8)	2(0.4)	1(0.2)	18(3.2)	2(0.4)	4(0.8)
Asian	232(5.2)	134(15.5)	19(4.0)	3(0.6)	50(9.8)	4(0.7)	10(1.7)	10(1.6)	0(0.0)	2(0.3)
Hispanic	35(1.3)	5(0.5)	1(0.2)	1(0.3)	7(1.3)	3(0.6)	3(0.9)	12(2.1)	2(0.4)	1(0.1)
Other	106(2.2)	47(5.9)	11(1.6)	6(1.3)	17(3.5)	8(1.4)	4(0.7)	8(1.3)	1(0.1)	4(0.8)
BMI, kg/m2	27.7(0.1)	26.7(0.2)	28.6(0.3)	28.2(0.3)	28.0(0.2)	28.5(0.3)	28.2(0.3)	27.0(0.3)	27.8(0.2)	29.1(0.3)	<0.0001
Education, years	15.6(0.1)	15.4(0.1)	15.2(0.2)	14.8(0.2)	15.5(0.2)	15.7(0.2)	16.3(0.2)	15.8(0.2)	14.8(0.2)	15.3(0.2)	<0.0001
Pack years of tobacco exposure	24.2(0.7)	24.6(1.9)	21.2(1.1)	25.8(1.7)	23.5(1.4)	23.0(1.3)	25.2(1.4)	25.9(1.3)	22.1(1.1)	21.9(1.1)	0.0445
Current smoker, n(%)	701(15.2)	117(13.5)	69(13.6)	73(14.4)	75(13.7)	67(13.1)	50(11.1)	118(20.0)	74(12.5)	58(12.5)	<0.0001
Ever Smoker, n(%)	2746(53.5)	447(53.1)	276(46.8)	250(46.1)	281(51.0)	288(52.3)	245(49.3)	323(58.5)	363(58.7)	273(56.9)	<0.0001
History of TB, n(%)	71(1.0)	26(3.1)	9(1.7)	9(1.9)	2(0.5)	5(0.7)	3(0.8)	6(1.2)	2(0.4)	9(1.5)	0.0001
History of childhood hospitalization, n(%)	294(4.6)	71(8.2)	40(7.1)	34(6.5)	20(3.9)	37(7.3)	20(4.2)	26(4.1)	26(4.6)	20(4.5)	0.0009
History of asthma, n(%)	825(16.3)	122(13.9)	91(16.2)	77(15.4)	79(13.8)	85(15.8)	93(20.0)	103(18.9)	78(14.1)	97(19.4)	0.0102
≥1 year in dusty job, n(%)	1474(28.0)	257(30.6)	177(32.4)	199(39.5)	149(28.0)	157(29.6)	119(25.2)	152(27.5)	122(21.7)	142(31.0)	<0.0001
Number of years in dusty job	11.5(0.4)	10.2(0.7)	12.5(1.1)	15.7(0.9)	12.1(0.9)	11.6(0.8)	10.3(1.0)	9.7(0.8)	14.9(1.1)	14.0(1.0)	<0.0001

Notes: Continuous data are shown as mean (SD) and categorical data are shown as the number of frequency, n (%); Data are weighted to the local population. Chi-square test for categorical variables and Kruskal–Wallis test for continuous variables; the overall p-value is the comparison between all sites.

Abbreviations: BMI, body mass index; TB, tuberculosis.

As shown in Table 1, the frequency of current smoking was lowest in Ottawa (11.1%) and highest in Montreal (20.0%,) with the heaviest tobacco exposure indicated by the number of pack-years (25.9) for Montreal. Childhood hospitalization for respiratory illness was most common in Vancouver and Kingston, 8.2% and 7.3%, respectively. Saskatoon had the highest percentage (39.5%) of participants with exposure to dusty job, and for the greatest number of years, 15.7 (details by sex are shown in and ).

Estimated Population Prevalence by Site

The nine sampling cities are identified on a map of Canada, collectively representing where the majority of the Canadian population resides (Figure 1). Based on airflow limitation defined by the LLN ratio, the overall population prevalence of COPD was highest in Montreal, 14.7% (95% CI, 11.4–18.0), and lowest in Quebec City, 8.3% (95% CI, 5.8–10.8).

Figure 1

Geographical locations of the nine sampling urban cities in Canada. Overall population prevalence of COPD expressed as percentages. Airflow limitation was defined by FEV1/FVC <5th percentile (LLN). These sites collectively represent where the majority of the Canadian population resides.

The estimated population prevalence for the presence of airflow limitation and its severity based on GOLD grades are shown for the nine study sites and separately for men and women in Table 2. In the overall estimated population prevalence, regardless of the definition (LLN or FR) used to define airflow limitation, the highest prevalence was in the same three locations: Montreal (LLN: 14.7%, FR: 19.0%), Kingston (LLN: 13.9%, FR: 19.9%) and Vancouver (LLN: 13.4%, FR: 19.3%).

Table 2

Estimated Population Prevalence of COPD Shown as Number (%) According to Different Definitions in the Canadian Obstructive Lung Disease (COLD) Cohort

	All Sites	Vancouver	Calgary	Saskatoon	Toronto	Kingston	Ottawa	Montreal	Quebec City	Halifax
Men
N	2096	344	209	226	221	229	226	223	258	160
Airflow obstruction (LLN ratio)	241(10.9)	44(14.5)	23(10.3)	31(11.5)	18(7.9)	36(14.4)	19(7.2)	34(16.0)	19(6.6)	17(8.6)
Airflow obstruction (fixed ratio < 70%)	428(18.3)	68(22.0)	46(19.3)	61(21.4)	41(16.8)	63(23.9)	37(12.9)	47(21.9)	35(12.3)	30(15.0)
GOLD 1	248(10.1)	40(12.7)	22(9.5)	38(13.8)	25(9.6)	45(16.6)	21(7.0)	23(11.7)	18(6.7)	16(7.7)
GOLD 2	157(6.7)	28(9.3)	22(9.1)	20(6.7)	14(6.1)	17(7.0)	14(5.5)	20(7.9)	13(4.3)	9(5.0)
GOLD 3–4	23(1.4)	0(0.0)	2(0.7)	3(0.9)	2(1.0)	1(0.3)	2(0.5)	4(2.4)	4(1.3)	5(2.3)
Women
N	2797	483	360	232	300	277	230	318	299	298
Airflow obstruction (LLN ratio)	316(11.4)	58(12.5)	42(10.8)	26(9.8)	29(10.1)	41(13.4)	24(10.6)	42(13.6)	26(9.6)	28(9.8)
Airflow obstruction (fixed ratio < 70%)	413(14.4)	77(16.8)	58(14.7)	43(16.1)	37(12.9)	50(16.4)	32(13.4)	53(16.7)	30(11.2)	33(11.0)
GOLD 1	219(7.4)	45(9.6)	29(7.4)	24(9.6)	16(5.5)	34(11.3)	19(8.3)	28(9.7)	10(4.2)	14(4.4)
GOLD 2	159(5.9)	24(5.5)	25(6.1)	17(6.0)	17(6.4)	13(4.3)	13(5.1)	21(6.0)	15(5.4)	14(4.6)
GOLD 3–4	35(1.0)	8(1.8)	4(1.1)	2(0.6)	4(1.1)	3(0.8)	0(0.0)	4(1.0)	5(1.6)	5(2.0)
Overall
N	4893	827	569	458	521	506	456	541	557	458
Airflow obstruction (LLN ratio)	557(11.2)	102(13.4)	65(10.6)	57(10.6)	47(9.1)	77(13.9)	43(9.0)	76(14.7)	45(8.3)	45(9.3)
Airflow obstruction (fixed ratio < 70%)	841(16.2)	145(19.3)	104(16.9)	104(18.5)	78(14.7)	113(19.9)	69(13.2)	100(19.0)	65(11.7)	63(12.9)
GOLD 1	467(8.6)	85(11.1)	51(8.4)	62(11.5)	41(7.4)	79(13.7)	40(7.7)	51(10.6)	28(5.3)	30(5.9)
GOLD 2	316(6.4)	52(7.3)	47(7.5)	37(6.3)	31(6.3)	30(5.6)	27(5.3)	41(6.8)	28(4.9)	23(4.8)
GOLD 3–4	58(1.2)	8(0.9)	6(0.9)	5(0.7)	6(1.0)	4(0.6)	2(0.2)	8(1.6)	9(1.4)	10(2.2)

Notes: Data are shown as n and % indicating number and population prevalence, respectively. Site-specific prevalence estimates are weighted to the local site population, sex and age profiles obtained from the Canadian population census 2006. GOLD Grades are defined as GOLD 1 = FEV1/FVC<0.7 and FEV1% predicted ≥80%; GOLD 2 = FEV1/FVC<0.7 and FEV1% predicted <80% ≥50%; GOLD 3–4 = FEV1/FVC<0.7 and FEV1% predicted <50%.

Table 2 shows that the overall population prevalence of COPD in men was 18.3% (n = 428) by FR and 10.9% (n = 241) by LLN ratio, while in women, the overall prevalence was 14.4% (n = 413) and 11.4% (n =316), respectively. Overall, across all nine Canadian cities, the prevalence was 16.2% (95% CI, 14.5–17.8) by FR and 11.2% (9.7–12.6) by LLN, with a difference of 5% between the two spirometric definitions. The majority of the population with airflow obstruction is classified as GOLD grade 1 (8.6%) and grade 2 (6.4%), indicating mild and moderate airflow limitation, respectively.

Heterogeneity of Risk Factors by Site and by Sex

The participants’ characteristics of age, history of asthma, history of childhood hospitalization, smoking history, smoking status, and cumulative exposure in pack-years and their association with COPD as defined by airflow limitation below the LLN were summarized as odds ratio adjusted for potential confounding covariates in forest plots (Figures 2–7). The forest plots show: 1) risk association by site for men and women separately and the overall pooled estimate; 2) the site heterogeneity denoted by I2 and its significance on the overall risk association with COPD in the population. The overall adjusted odds ratio [aOR, 95% CI] for age was 1.56 (1.33, 1.84): 1.64 (1.30, 2.08) in men; 1.49 (1.18, 1.88) in women (Figure 2). A history of physician-diagnosed asthma (Figure 3) and childhood hospitalization for respiratory illness (Figure 4) were risk factors, with overall OR 3.30 (2.42, 4.49) and 1.81 (1.17, 2.80), respectively. In particular, men with a history of asthma were at an even higher risk for COPD, OR 4.28 (2.18, 8.37), compared to women, OR 2.86 (1.98, 4.14). With regards to smoking as a risk factor for COPD, our data compared several aspects, including ever vs never-smoking history, current vs former smoking status, and number of pack-years. Current smoking was associated with the highest odds ratio for COPD (OR 3.49; 95% CI 2.55, 4.80) (Figure 5), followed by ever-smoking (OR 2.25; 95% CI 1.81, 2.78) (Figure 6), and total pack-years smoked (OR 1.32; 95% CI 1.22, 1.43) (Figure 7). There was significant heterogeneity among the nine sampling sites in age and pack-years of smoking history, I-squared = 40.8%, p = 0.037 and 48.5%, p = 0.011, respectively, but not for the other risk factors.

Figure 2

Age as a risk factor for COPD, showing significant heterogeneity in risk among the nine sampling sites. Results are shown for men and women and for the whole cohort. The x-axis is the log-scaled, but the values labelled under the line are the actual odds ratio. COPD was defined by airflow limitation (FEV1/FVC < LLN); ES (95% CI) is the adjusted odds ratio and 95% confidence interval, aOR (95% CI), for every 10-year increase in age, adjusted for school years, pack-years, asthma, childhood hospitalization, and dusty job exposure.

Figure 3

Asthma as a risk factor for COPD across nine sites. Results are shown for men and women and for the whole cohort. COPD was defined by airflow limitation (FEV1/FVC < LLN); ES (95% CI) is the adjusted odds ratio and 95% confidence interval, aOR (95% CI), adjusted for age, school years, pack-years, and dusty job exposure. There is no site heterogeneity demonstrated.

Figure 4

Childhood hospitalization as a risk factor for COPD across nine sites. Results are shown for men and women and for the whole cohort. COPD was defined by airflow limitation (FEV1/FVC < LLN); ES (95% CI) is the adjusted odds ratio and 95% confidence interval, aOR (95% CI), adjusted for age, school years, pack-years, and dusty job exposure There is no site heterogeneity demonstrated.

Figure 5

Current vs former smoking as a risk factor for COPD across nine sites. Results are shown for men and women and for the whole cohort. COPD was defined by airflow limitation (FEV1/FVC < LLN); ES (95% CI) is the adjusted odds ratio and 95% confidence interval, aOR (95% CI), adjusted for age, school years, asthma, childhood hospitalization, and dusty job exposure. There is no significant site heterogeneity demonstrated.

Figure 6

Ever-smoking as a risk factor for COPD across nine sites. Results are shown for men and women and for the whole cohort. COPD was defined by airflow limitation (FEV1/FVC < LLN); ES (95% CI) is the adjusted odds ratio and 95% confidence interval, aOR (95% CI), adjusted for age, school years, asthma, childhood hospitalization, and dusty job exposure. There is no significant site heterogeneity demonstrated.

Figure 7

Smoking in pack-years as a risk factor for COPD, showing significant site heterogeneity in risk among the nine sampling sites. Results are shown for men and women and for the whole cohort. COPD was defined by airflow limitation (FEV1/FVC < LLN); ES (95% CI) is the adjusted odds ratio and 95% confidence interval, aOR (95% CI) for every 10-pack-year increase, adjusted for age, school years, asthma, childhood hospitalization, and dusty job exposure.

The impact of these risk factors is further magnified with increasing spirometric severity as demonstrated by the higher odds ratio in moderate-severe COPD (Table 3), particularly in those with a history of asthma, aOR 5.22 (3.54, 7.70).

Table 3

Overall Adjusted Odds Ratio (aOR) for the Risk Factors for COPD of Different Severity: Mild (GOLD Grade 1) and Moderate-Severe (GOLD Grade 2–4)

	Mild COPD (GOLD Grade 1)aOR	Moderate-to-Severe COPD (GOLD Grade 2–4)aOR
Age, years	1.99 (1.80, 2.20)	2.22 (1.87, 2.63)*
Asthma (Y vs N)	2.69 (2.01, 3.59)	5.22 (3.54, 7.70)*
Childhood hospitalization (Y vs N)	1.89 (1.21, 2.94)	3.07 (2.02, 4.66)
Ever-smoker status (Y vs N)	1.65 (1.33, 2.05)	3.05 (2.24, 4.16)
Lifetime tobacco exposure, pack-years#	1.18 (1.10, 1.28)	1.36 (1.25, 1.48)*
Smoking status (current vs former)	2.34 (1.67, 3.29)	4.09 (2.78, 6.01)

Notes: GOLD Grades are defined as GOLD 1 = FEV1/FVC < 0.7 and FEV1% predicted ≥80%; GOLD 2 = FEV1/FVC < 0.7 and FEV1% predicted ≥50% and <80%; GOLD 3 = FEV1/FVC < 0.7 and FEV1% predicted ≥30% and <50%; GOLD 4 = FEV1/FVC < 0.7 and FEV1% predicted <30%. *Indicates presence of significant site heterogeneity. #Pack years = (number of cigarettes per day/20) × number of years smoked. Separate models were constructed for each risk factor. Age adjusted for school years and ever-smoking history. Asthma adjusted for age, school years, and ever-smoking history. Childhood hospitalization adjusted for age, school years, and ever-smoking history. Ever-smoker status adjusted for age and school years. Lifetime tobacco exposure adjusted for age and school years. Smoking status adjusted for age and school years.

Abbreviation: aOR, overall adjusted odds ratio for men and women.

Table 4 shows the results of the additional heterogeneity analysis to test whether the variation in prevalence for COPD among sites is substantially reduced by adjustment for risk factors in this study. The results of Tau2 analyses further confirmed that there was significant heterogeneity among the nine sites without the influence of risk factors, with COPD as the outcome. When these risk factors were introduced into the model and Tau2 was adjusted for sex, age, smoking status, asthma, and childhood hospitalization, the variance decreased. The mean difference between the two analyses was significant, 0.023 (0.001, 0.037); p < 0.001. The variation among sites is therefore explained by differences in risk factors.

Table 4

Heterogeneity (Tau2) for the Prevalence of COPD (LLN) Among the Nine Sites

	Tau2 without Risk Factors	*Tau2 with Risk Factors	#Mean Difference (95% CI)	p-value
Prevalence of COPD (LLN)	0.053	0.030	0.023 (0.001,0.037)	<0.001$

Notes: *Tau2 is adjusted for sex, age, smoking status, asthma, and childhood hospitalization; #Mean difference for Tau2 (95% CI) with and without risk factors; $Indicates a significant difference of Tau2 with and without risk factors.

Discussion

The present study reported population prevalence for COPD in Canada using two spirometric criteria and provided new data on the heterogeneity of prevalence across the country. We found heterogeneity in age and pack-years of smoking (reflecting exposure burden) across sites, while a history of asthma, childhood hospitalization, and smoking patterns were homogeneous across sites in Canada; and that the variation in prevalence across sites can be reduced by adjustment for these risk factors. We have also shown that the risk factors associated with COPD are qualitatively similar across all grades of severity of COPD but the strength of association increased with severity, suggesting that a higher risk could be linked to the development of severe disease. Taken together, accurate population data on spirometry-based COPD prevalence and a better understanding of the heterogeneity of major risk factors in different regions of the country can help inform public policymakers on target strategies in reducing COPD and its severity in the country.

COPD Prevalence in Canada

Determination of spirometry-based population prevalence is crucial to understanding the impact of COPD in Canada. Previous reports on the prevalence of COPD in Canada were largely based on surveys with self-reported history of chronic bronchitis or emphysema,20 administrative data-based criteria in single or multiple centers,21,22 or criteria that only relied on pre-bronchodilator spirometry.23 Following the methodology in the BOLD study,9 the use of random sampling and a standardized protocol, including the definition of COPD by strict post-bronchodilator criteria, has allowed for more accurate determination of prevalence. Although some data on COPD prevalence on five Canadian sites were previously published,6 this study adds to the literature by providing robust and extensive data on the Canadian population prevalence of COPD documented by the current two widely used clinical measures of post-bronchodilator FEV1/FVC ratio (FR < 0.70 and LLN <5th percentile). Our study emphasizes the importance of the spirometric definition when quoting prevalence of COPD and showed the impact on the interpretation on sex difference in prevalence.

An interesting observation is that more than half of the overall COPD cases were classified as mild (GOLD 1, FEV1/FVC <0.7 and FEV1% predicted >80%) in the population. This pattern is expected from a population-based sample of non-institutionalized men and women. Yet, it is pertinent to health-care planning as this group though “mild” can still experience exacerbations and utilize health-care resources with a significant economic and social impact.24 These patients should be provided education on COPD, and be encouraged on lifestyle optimization and smoking cessation to prevent disease progression. Early intervention is essential as our findings have shown that the impact of smoking-related risk factors is further magnified in more severe disease.

The Canadian prevalence of GOLD grade 2 to 4 COPD (FEV1/FVC <0.7 and FEV1% predicted <80%) is lower (7.5, SD 0.6) compared to that reported in the BOLD study (10.1%, SE 4.8).9 In the PLATINO study, the prevalence of COPD varied significantly among the five Latin American locations, ranging from 7.8% (95% CI, 5.9–9.7) in Mexico City to almost 20% (95% CI, 17.2–22.2) in Montevideo.11 In a systematic review that did not specifically use spirometry to define COPD, the global pooled prevalence of COPD was 9.9% (95% CI, 8.2–11.8) in the population over 40 years old.25

Heterogeneity of Risk Factors

Our study indicates that the burden of main risk factors underlying the development of COPD differs among Canadian cities. It is conceivable that the heterogeneity of COPD prevalence within Canada may also result from other unique factors not included in this analysis. For example, geographic variation including altitude and temperature has also been proposed as explanation for differences in COPD prevalence in other countries.11,26 Although sex, age, and smoking contribute to two-thirds of the variability in COPD prevalence, a significant portion of the variability remains unexplained.26

Our study found that a history of asthma or childhood hospitalization for respiratory illnesses is just as important as smoking in the development of spirometrically confirmed COPD. While long-term cigarette smoking behaviour is the most important and well-studied cause of COPD in the Western World,27–29 COPD can occur in never-smokers which suggests the importance of other risk factors in determining an individual’s susceptibility to COPD.2,30 In never-smokers, occupational exposures including organic and inorganic dusts, chemical agents, and fumes are related to the development of COPD. In patients with established COPD, these exposures can worsen the morbidity, quality of life, and exacerbation risk.31,32 Medical events during gestation and birth, and exposures during childhood and adolescence affect lung growth and subsequent function in adult life.33,34 A history of severe childhood respiratory infection and hospitalization is associated with reduced lung function and increased respiratory symptoms in adulthood.35,36 In a prospective observational study over 20 years, subjects with physician-confirmed active asthma were identified as a risk factor for COPD (HR 12.5, 95% CI 6.83–22.8), chronic bronchitis (HR 10.0, 4.90–20.2), and emphysema (HR 17.0, 8.30–34.8).37 It is intriguing that the risk burden of asthma and childhood hospitalization for COPD, unlike the risk burden of exposure to cigarette smoke, did not show regional variation in Canada.

Gender Differences

Recent research on gender-specific prevalence suggested that the prevalence has been rising more rapidly in females compared to males.38,39 The discordance in the sex difference in COPD prevalence defined by the two definitions of fixed ratio and lower limits of normal ratio cut-offs is intriguing. We think this is explained by the fact that men have larger absolute lung volumes than women (ie increased FVC); thus, it is easier to develop a reduced FEV1/FVC ratio with progressive airflow impairment. This issue can be largely mitigated when the FEV1/FVC ratio is adjusted for age, sex and ethnicity and presented as % of predicted value.40

The phenomenon of higher prevalence for men when using fixed ratio and lower or almost similar sex prevalence when using LLN has been shown previously in a paper from the Netherlands where the use of LLN for stage 1 [FEV1/FVC

Finally, our finding that significant site heterogeneity in the risk for COPD due to cumulative tobacco exposure (pack-years) was seen in men but not in women suggests that the site relative proportion of men and women could also be an explanation for some of the variation in prevalence, consistent with our previous observation that site variation in prevalence appeared to be driven by age and sex.6

Limitations

There are several limitations of this study. The presented data may not necessarily reflect rural settings and certain ethnic groups in Canada. For example, the Aboriginal population of Canada has a higher incidence and burden of disease compared to non-Aboriginal Canadian population, likely as a result of disparities in risk factors including smoking.41 A previous systematic review demonstrated that COPD prevalence is generally lower in rural settings.42 Nonetheless, our study obtained samples from cities where the collective majority of the Canadian population resided. The delay in reporting these data, which were collected more than 10 years ago, is a limitation due to the potential influence of changes in age profile and ethnic distribution of the population over the past decade. Comparing the Canadian Census data in 2006 and 2016, the proportion of individuals over age 40 has increased only slightly by 2.3% (). The ethnicity of the Canadian population has shifted towards an increase in non-white group by 7.6% between 2006 and 2016 (). Visible minority Canadians were underrepresented in this study. The gap in health data and research on visible minorities in Canada has been previously documented43 and may be due to various barriers to participation. Finally, some subjects with airflow limitation may have incompletely reversible asthma rather than COPD in this study based on spirometrically defined COPD.

Conclusions

Our study has confirmed that COPD remains a significant health issue in Canada, as shown by the population prevalence based on post-bronchodilator spirometry and defined by two globally popular definitions, thus, providing representative and reliable population estimates of this disease. Sex difference in prevalence is influenced by the choice of spirometric definition for COPD. Significant site heterogeneity in risk factor was demonstrated in the burden of tobacco exposure (increasing pack-years) rather than the pattern of smoking and that the variation in prevalence across sites can be reduced by adjustment for risk factors. The strength of association between COPD and the analysed risk factors was magnified with disease severity. The study emphasizes the importance of understanding regional heterogeneity of risk factors for COPD, as this knowledge is essential for interpretation of the prevalence of COPD and help inform health policy planning for COPD within a country. More studies are needed to better understand the link between variations in risk factors and COPD.