Socioeconomic status and COPD among low- and middle-income countries.

BACKGROUND: Socioeconomic status (SES) is a strong social determinant of health. There remains a limited understanding of the association between SES and COPD prevalence among low- and middle-income countries where the majority of COPD-related morbidity and mortality occurs. We examined the association between SES and COPD prevalence using data collected in Argentina, Bangladesh, Chile, Peru, and Uruguay.
METHODS: We compiled lung function, demographic, and SES data from three population-based studies for 11,042 participants aged 35-95 years. We used multivariable alternating logistic regressions to study the association between COPD prevalence and SES indicators adjusted for age, sex, self-reported daily smoking, and biomass fuel smoke exposure. Principal component analysis was performed on monthly household income, household size, and education to create a composite SES index.
RESULTS: Overall COPD prevalence was 9.2%, ranging from 1.7% to 15.4% across sites. The adjusted odds ratio of having COPD was lower for people who completed secondary school (odds ratio [OR] =0.73, 95% CI 0.55-0.98) and lower with higher monthly household income (OR =0.96 per category, 95% CI 0.93-0.99). When combining SES factors into a composite index, we found that the odds of having COPD was greater with lower SES (interquartile OR =1.23, 95% CI 1.05-1.43) even after controlling for subject-specific factors and environmental exposures.
CONCLUSION: In this analysis of multiple population-based studies, lower education, lower household income, and lower composite SES index were associated with COPD. Since household income may be underestimated in population studies, adding household size and education into a composite index may provide a better surrogate for SES.

Background

COPD is a preventable and treatable disease defined by airflow obstruction that is not fully reversible.1,2 The leading causes of COPD worldwide include tobacco smoking and exposure to household air pollution from biomass fuels, although this may change with increasing exposure to tobacco smoking particularly among low- and middle-income countries (LMICs).3 Estimates show that COPD affects 328 million people globally, though population-based data examining the prevalence, risk factors, and clinical outcomes of COPD are sparse outside high-income countries (HIC).4,5 A few population-based studies have been conducted over the last decade that found a wide variation in COPD prevalence among populations.5–7

Socioeconomic status (SES) is one of the most powerful determinants of health across a range of chronic diseases.8,9 SES is defined as an individual’s social and economic standing and serves as proxy for social or economic position or rank in a social group. More than a measure of income, SES encompasses several other measures including education, occupation, housing, assets, and participation in social organization.10,11 While lower SES was found to be associated with greater COPD morbidity and mortality in HIC, few studies have examined the role of low SES in the prevalence of COPD among LMICs.12–14 An understanding of the role of SES and COPD in LMICs will potentially inform public interventions beyond harm reduction, ie, tobacco cessation and decreased biomass fuel smoke exposure, toward those aimed at disparities in SES.

Factors limiting a comprehensive assessment of the relationship between SES and COPD in LMICs include the high cost of spirometry and limited data from which to analyze associations. Here, we describe the relationship between SES and COPD prevalence among the ten low- and middle-income settings in Argentina, Bangladesh, Chile, Peru, and Uruguay, which encompass a diversity of geographies, ethnicity, and variations in altitude and degree of urbanization.

Methods

Study setting

This study used pooled data from multiple, population-based studies sponsored by the United States National Institutes of Health/National Heart, Lung, and Blood Institute and UnitedHealth Chronic Disease Initiative (http://www.nhlbi.nih.gov/about/org/globalhealth/centers). Specifically, data were compiled from the Pulmonary Risk in South America (PRISA) study conducted by the Institute for Clinical Effectiveness and Health Policy (IECS) in two sites in Argentina (Marcos Paz and Bariloche), one in Chile (Temuco), and one in Uruguay (Canelones), the CRONICAS Cohort Study in Peru conducted by CRONICAS Centre of Excellence for Chronic Diseases at Universidad Peruana Cayetano Heredia, and a longitudinal study in Bangladesh conducted by the Centre for Control of Chronic Diseases at the International Centre for Diarrheal Disease Research, Bangladesh (icddr,b), totaling ten settings of diverse geography and SES in five LMICs. Both PRISA and CRONICAS studies are prospective longitudinal studies with multiple years of follow-up that started in 2010.15,16 The icddr,b conducted a longitudinal study from 2011 to 2012.17 The overall cohort is a mixture of rural, semiurban, and urban populations.15–18

Study design

Both PRISA and CRONICAS used age- and sex-stratified random sampling, whereas the Bangladesh study used simple random sampling of available census data at each site.15–18 All three studies aimed to enroll between 500 and 1,800 participants per site. The PRISA study included participants aged 45–75 years, whereas the CRONICAS and Bangladesh studies enrolled participants ≥35 years and ≥40 years of age, respectively.15–18 We limited the age range to 35–95 years because of the lack of reference equations in older adults.19 All studies required permanent residency in the sampling areas and ability to provide informed consent.15–18 Common exclusion criteria for the PRISA and CRONICAS studies included anyone who had active pulmonary tuberculosis and was pregnant.15,16 Only one person per household was enrolled in all studies.15–17 Informed consent was obtained from all participants, and all field workers completed confidentiality training.15–18 All studies obtained approval from their corresponding local and international internal review boards.15–18

Spirometry

Technicians performed spirometry with EasyOne (PRISA and Bangladesh) and Easy-On-PC spirometers (CRONICAS), both manufactured by ndd Medical Technologies (Zurich, Switzerland).15–18 All studies followed standard guidelines for the conduct of spirometry.15–18 Post-bronchodilator spirometry was conducted 15 minutes after administration of inhaled salbutamol in Bangladesh (5 mg via nebulizer) and CRONICAS (200 μg via spacer) or albuterol in PRISA (200 μg via spacer).16–18 PRISA and CRONICAS conducted pre- and post-bronchodilator spirometry on all participants, whereas Bangladesh only conducted post-bronchodilator spirometry among those who had evidence of obstruction, i.e., a post-bronchodilator forced expiratory volume in one second (FEV1) over forced vital capacity (FVC) ratio less than 0.70.17 Methods for collecting other data for these studies are described in previous publications.16–18

Definitions

We defined COPD as having a post-bronchodilator FEV1/FVC below the lower limit of normal (5th percentile) of the Global Lung Function Initiative (GLI) reference values for mixed ethnic populations20 and conducted sensitivity analyses using Caucasian reference values. Severe COPD was classified as additionally having a post-bronchodilator FEV1 <50% predicted.20 For comparability with other studies, we also used the FEV1/FVC <0.70 definition for COPD. Three main indicators of SES that were commonly available across studies were selected for this analysis: household size, education, and monthly household income. We categorized monthly income into six categories of increasing order at each site; education as below secondary school or secondary school and above. Smoking was classified as daily smoker (≥1 cigarettes/day) or not, while biomass fuel smoke exposure was categorized as using biomass as a primary source of fuel or not. In sensitivity analyses, we categorized smoking as current, former, and never smoking and also calculated pack-years of smoking.

Biostatistical methods

This study had two main analytical objectives: to evaluate the association between SES and COPD prevalence and to evaluate the association between SES and pre-bronchodilator FEV1/FVC Z-scores. We limited our analysis to pre-bronchodilator FEV1/FVC values because the Bangladesh study only conducted post-bronchodilator spirometry among participants with a pre-bronchodilator FEV1/FVC of <0.7. We used multiple variable logistic regressions to evaluate the association between each of the SES indicators and COPD adjusted for the following a priori selected variables: age, sex, and self-reported biomass fuel smoke exposure and daily smoking. Principal component analysis was then performed on monthly household income, household size, and education to create a composite index based on the first principal component (proportion of variance 0.91), which was then included in a multivariable logistic regression model adjusted for the following a priori selected variables: age, sex, and self-reported biomass fuel smoke exposure and daily smoking.21,22 Since there were multiple sites in these analyses leading to a nested design of site and individual within site, we used alternating logistic regressions to account for correlation.23 We scaled the effect of lower composite SES index to the interquartile odds ratio (OR), ie, we calculated the ratio of odds of COPD between the values in the 25th and 75th percentiles of the composite SES index.

We used linear regression to evaluate the association between pre-bronchodilator FEV1/FVC Z-scores and each SES indicator and SES composite index, adjusted for the following a priori selected variables: age, sex, and self-reported biomass fuel smoke exposure and daily smoking. We used linear mixed-effects models with a random intercept by site to model heterogeneity by site24 and scaled the effect of lower composite SES index to the interquartile difference, ie, we calculated absolute difference in pre-bronchodilator FEV1/FVC Z-scores between 75th and 25th percentiles of the composite SES index.

For all the abovementioned analyses, we conducted sensitivity analyses by replacing daily smoking with categories of smoking history, and also with pack-years of smoking. All analyses were performed using R (www.r-project.org).

Results

Participant characteristics

A total of 11,042 participants met inclusion criteria from all ten sites: Bariloche (n=1,105), Canelones (n=938), Dhaka (n=1,878), Lima (n=997), Marcos Paz (n=1,242), Matlab (n=1,846), rural Puno (n=505), Temuco (n=1,063), Tumbes (n=958), and urban Puno (n=510). Pre-bronchodilator lung function measurements were available for 11,042 individuals. Average age was 56.1±10.4 years, with males accounting for 45% of the sample. There was incomplete information for household size, monthly household income, education, or daily smoking in 1,522 participants. There were no differences in prevalence of COPD (P=0.34) or sex (P=0.08) between participants with and without complete data adjusted for site. On average, participants with incomplete data were 0.9 years older (P<0.01) than those with complete data. The interquartile difference for the composite SES index was 2.90 (25th and 75th percentiles were −1.85 and 1.05, respectively). We summarized other sociodemographic and lifestyle characteristics by site (Table 1).

Table 1

Sociodemographic characteristics and COPD prevalence stratified by site

Variables	Bariloche	Canelones	Dhaka	Lima	Marcos Paz	Matlab	Rural Puno	Temuco	Tumbes	Urban Puno
COPD positive, n (%)	82 (7)	96 (10)	182 (10)	56 (6)	175 (14)	284 (15)	43 (9)	52 (5)	16 (2)	25 (5)
SES composite index, mean (SD)	0.60 (1.34)	−0.63 (1.37)	1.02 (1.48)	−0.09 (1.06)	−0.09 (1.15)	−0.55 (1.44)	−1.67 (0.82)	0.81 (1.89)	−0.73 (0.93)	−0.04 (1.27)
SES composite index, interquartile difference (25th and 75th percentile)	2.08 (−0.53, 1.55)	1.09 (−1.93, −0.85)	2.90 (0.50, 2.95)	1.09 (−0.94, 0.14)	1.09 (−1.03, 0.06)	1.99 (−1.59, 0.51)	1.09 (−2.03, −0.95)	3.98 (−0.40, 3.59)	1.09 (−1.04, 0.05)	1.99 (−1.40, 0.59)
Age in years, median (IQR)	56.7 (50.5, 63.0)	59.2 (52.1, 66.2)	50.0 (43.5, 56.5)	54.7 (45.6, 63.7)	58.4 (51.6, 65.3)	53.0 (45.0, 61.0)	54.9 (45.6, 64.3)	58.8 (51.6, 66.0)	54.8 (44.9, 64.6)	54.7 (44.9, 64.4)
Number of males (%)	431 (39)	375 (40)	853 (45)	491 (49)	498 (40)	848 (51)	239 (47)	484 (46)	480 (50)	250 (49)
Height in m, median (IQR)	1.61 (1.54, 1.68)	1.62 (1.55, 1.68)	1.55 (1.48, 1.62)	1.54 (1.48, 1.60)	1.61 (1.54, 1.68)	1.53 (1.47, 1.60)	1.55 (1.49, 1.61)	1.60 (1.53, 1.67)	1.58 (1.52, 1.65)	1.57 (1.50, 1.64)
BMI in kg/m2, mean (SD)	28.9 (5.6)	29.6 (6.2)	24.5 (4.8)	28.4 (4.5)	30.1 (5.8)	20.6 (3.6)	25.2 (3.7)	29.1 (4.9)	28.4 (4.8)	27.9 (4.4)
BMI ≥30 kg/m2, n (%)	417 (38)	398 (42)	234 (13)	324 (33)	577 (47)	29 (2)	54 (10.7)	389 (37)	306 (32)	136 (27)
Daily smokers, n (%)	263 (24)	205 (22)	119 (6)	32 (3)	284 (23)	165 (9)	1 (0)	161 (15)	53 (6)	11 (2)
Pack-years of smoking, mean (SD)	19.7 (7.9, 31.4)	27.7 (12.8, 42.6)	15.6 (4.5, 26.7)	0.2 (0.0, 1.0)	28.7 (16.8, 40.6)	15.6 (5.9, 25.4)	0.1 (0.0, 0.3)	8.5 (1.9, 15.2)	0.6 (0.0, 2.9)	0.1 (0.0, 0.8)
Biomass as primary source of fuel, n (%)	254 (23)	43 (5)	66 (4)	24 (3)	6 (1)	1,805 (98)	481 (95)	233 (22)	163 (17)	10 (2)
Household size, median (IQR)	9 (8, 10)	9 (8, 10)	5 (4, 6)	5 (3.5, 6.5)	9 (7.5, 10.5)	5 (3.5, 6.5)	3 (1.5, 4.5)	9 (7.5, 10.5)	4 (3, 5)	4 (3, 5)
Education ≥ secondary, n (%)	548 (50)	469 (50)	815 (64)	456 (46)	426 (34)	346 (19)	169 (34)	752 (71)	347 (36)	329 (65)
Monthly household income categories, n (%)
Lowest	24 (3)	86 (10)	89 (5)	2 (0)	15 (3)	360 (20)	148 (40)	68 (9)	11 (1)	8 (2)
Second	86 (9)	368 (41)	91 (5)	179 (19)	101 (18)	384 (21)	176 (47)	123 (17)	439 (47)	110 (26)
Third	237 (25)	292 (32)	315 (17)	432 (46)	224 (40)	541 (29)	41 (11)	115 (16)	354 (38)	136 (31)
Fourth	239 (25)	43 (5)	374 (20)	200 (21)	126 (23)	273 (15)	7 (2)	45 (6)	83 (9)	96 (22)
Fifth	184 (20)	12 (1)	359 (19)	83 (9)	66 (12)	174 (9)	0 (0)	22 (3)	30 (3)	45 (10)
Highest	170 (18)	101 (11)	646 (34)	52 (5)	28 (5)	114 (6)	3 (1)	357 (49)	20 (2)	37 (9)

Abbreviations: SES, socioeconomic status; IQR, interquartile range; BMI, body mass index.

COPD prevalence

Overall prevalence of COPD was 9.2% with a range of 1.7% (Tumbes) to 15.4% (Matlab). The prevalence of COPD by site and sex is shown in Figure 1. COPD prevalence was higher among participants who were identified as daily smokers vs nondaily smokers and nonsmokers (18.5% vs 7.9%), used biomass as a primary source of fuel vs not (12.8% vs 7.8%), and males vs females (12.4% vs 6.5%). Severe COPD was 1.5% with a range of 0.0% (urban and rural Puno) to 2.8% (Matlab). When the fixed cutoff (post-bronchodilator FEV1/FVC <0.7) was used, COPD prevalence was 10.9% with a range of 3.5% (Tumbes) to 18.0% (Marcos Paz).

Figure 1

Prevalence of COPD per city using the lower limit of normal (LLN) method of diagnosis with Global Lung Function Initiative (GLI) mixed ethnic population, presented as overall and stratified by sex.

Lung function

Average pre-bronchodilator FEV1, FVC, and FEV1/FVC values were 2.34±0.75 L, 3.08±0.94 L, and 0.76±0.09, respectively. The distributions of pre-bronchodilator FEV1/FVC Z-score per site are shown in Figure 2. Bariloche, Canelones, Marcos Paz, and Temuco had markedly more outliers below the zero line than did other sites. Urban Puno had the highest average FEV1 and FVC at 2.80±0.83 L and 3.70±1.07 L, respectively, while Dhaka had the lowest at 1.90±0.57 L and 2.49±0.71 L. Tumbes had the highest mean pre-bronchodilator FEV1/FVC at 0.79 and 0.81. In Figure 3, we show that altitude of residence was associated with a higher pre-bronchodilator FEV1 (P=0.04) or FVC Z-score (P=0.01). On the other hand, there was not a clear relationship between pre-bronchodilator FEV1/FVC Z-score and altitude of residence (P=0.33).

Figure 2

Pre-bronchodilator FEV1/FVC Z-score using Global Lung Function Initiative (GLI) mixed ethnic population, stratified by age and sex.

Notes: Age is stratified into four categories. Males are represented by red box plots and females by green, in ten panels each representing an individual site.

Figure 3

Mean pre-bronchodilator FEV1 and FEV1/FVC Z-score using Global Lung Function Initiative (GLI) mixed ethnic population, by elevation in meters and site.

Notes: Elevation is in meters on the x-axis with Z-score on the y-axis. City names are beside each corresponding point. The lines represent the fitted values of a regression line of Z-score by elevation.

Determinants of COPD

Analysis using multivariable alternating logistic regression showed that the odds of having COPD was positively associated with use of biomass as a primary fuel source for cooking, daily smoking, not having completed secondary education, lower monthly household income, and older age (Table 2). We found a positive association between the odds of COPD and a lower SES composite score, with an interquartile OR of 1.23 (95% CI 1.05–1.43). In sensitivity analyses, the positive association between low SES and COPD remained significant (interquartile OR 1.29, 95% CI 1.12–1.50) when we used the GLI Caucasian reference population.

Table 2

Unadjusted and adjusted alternating logistic regression models for COPD using the lower limit of normal (LLN) method of diagnosis with Global Lung Function Initiative (GLI) mixed ethnic population, coefficients reported as odds ratios

Variables	Unadjusted		Adjusted model (individual SES indicators)		Adjusted models (SES composite index)
	OR (95% CI)	P-value	OR (95% CI)	P-value	OR (95% CI)	P-value
Monthly household income (change per category)	0.93 (0.87, 1.00)	0.04	0.96 (0.93, 0.99)	0.01
Median household size (household size <4 as reference)	1.10 (0.93, 1.32)	0.26	1.08 (0.90, 1.30)	0.42
Secondary school or higher (below secondary as reference)	0.71 (0.51, 0.98)	0.04	0.73 (0.55, 0.98)	0.04
Composite SES index (interquartile difference)	1.23 (1.01, 1.50)	0.04			1.23 (1.05, 1.43)	0.01
Age (per decade increase)	1.48 (1.35, 1.62)	<0.001	1.41 (1.31, 1.51)	<0.001	1.43 (1.33, 1.55)	<0.001
Female (male as reference)	0.46 (0.29, 0.72)	<0.001	0.51 (0.34, 0.78)	<0.001	0.54 (0.36, 0.80)	0.002
Daily smokers (nondaily and nonsmokers as reference)	2.66 (2.10, 3.37)	<0.001	2.25 (1.74, 2.92)	<0.001	2.25 (1.76, 2.89)	<0.001
Biomass as primary source of fuel	1.73 (1.44, 2.09)	<0.001	1.52 (1.23, 1.89)	<0.001	1.59 (1.31, 1.94)	<0.001

Abbreviations: SES, socioeconomic status; OR, odds ratio.

In further sensitivity analyses, we analyzed differences using pack-years and categories of smoking history (current, former, or never smoker) instead of daily smoking or not, while adjusting for age, sex, biomass fuel smoke exposure, and site. The interquartile OR for SES was 1.20 (95% CI 1.01–1.41) when using pack-years and 1.18 (95% CI 1.03–1.34) when using categories of smoking history. These results were very similar to using daily smoking vs not (interquartile OR for SES of 1.23, 95% CI 1.05–1.43).

We then plotted the overall sample adjusted ORs (left panel) and site-specific heterogeneity in adjusted ORs (right panel) for all SES indicators (Figure 4). Adjusted ORs of COPD by site were all in the same direction for secondary education and in the majority of sites for monthly household income (except for Canelones, Temuco, and Tumbes) and for the interquartile difference in composite SES index (except for Temuco, Tumbes, and Urban Puno).

Figure 4

Odds ratio of COPD by socioeconomic status composite index (SES index), secondary education or higher (below secondary as reference), greater than or equal to median household size of 4 (below as reference), and monthly household income (change per category).

Notes: Adjusted odds ratios are displayed for the overall sample (left) and for individual sites (right). In the left panel, the square represents the adjusted odds ratio and the segments represent 95% CIs. In the right panel, the triangles represent site-specific adjusted odds ratios. For the composite SES index, we calculate the interquartile odds ratio of COPD, ie, the ratio of odds of COPD between the 25th vs 75th percentile of the composite SES index.

Abbreviation: SES, socioeconomic status.

Determinants of pre-bronchodilator FEV1/FVC values

Lower pre-bronchodilator FEV1/FVC Z-scores were associated with individuals who were daily smokers, used biomass as a primary fuel source for cooking, males, the elderly, lower than secondary education, and those who had a lower monthly household income. Household size was not statistically significant (Table 3). A lower SES composite index was associated with lower pre-bronchodilator FEV1/FVC Z-scores. Specifically, a decrease equal to the inter-quartile range was associated with a 0.14 lower (95% CI 0.09–0.19) pre-bronchodilator FEV1/FVC Z-score. When using the Caucasian reference instead of that of a mixed ethnic population, the overall interquartile difference was decreased by a similar amount (0.13, 95% CI 0.08–0.18). We plotted the overall sample-adjusted differences in FEV1/FVC Z-score (left panel) and site-specific heterogeneity (right panel) for all SES indicators (Figure 5). Adjusted differences in FEV1/FVC Z-score by site were all in the same direction for secondary education, for monthly household income, and for the interquartile difference in composite SES index.

Table 3

Unadjusted and adjusted random effects linear regression models for pre-bronchodilator FEV1/FVC Z-score using Global Lung Function Initiative (GLI) mixed ethnic population

Variables	Unadjusted		Adjusted model (individual SES indicators)		Adjusted models (SES composite index)
	Difference (95% CI)	P-value	Difference (95% CI)	P-value	Difference (95% CI)	P-value
Monthly household income (change per category)	0.05 (0.03, 0.07)	<0.001	0.04 (0.02, 0.06)	<0.001
Median household size (household size <4 as reference)	0.04 (−0.01, 0.09)	0.13	0.03 (−0.03, 0.08)	0.32
Secondary school or higher (below secondary as reference)b	0.12 (0.17, 0.77)	<0.001	0.10 (0.05, 0.15)	<0.001
Composite SES index (interquartile difference)	−0.15 (−0.16, −0.13)	<0.001			−0.14 (−0.09, −0.19)	<0.001
Age (per decade increase)	−0.08 (−0.11, −0.06)	<0.001	−0.07 (−0.09, −0.05)	<0.001	−0.08 (−0.10, −0.06)	<0.001
Female (male as reference)	0.19 (0.14, 0.23)	<0.001	0.18 (0.13, 0.23)	<0.001	0.17 (0.12, 0.21)	<0.001
Daily smokers (nondaily and nonsmokers as reference)	−0.40 (−0.47, −0.33)	<0.001	−0.41 (−0.49, −0.33)	<0.001	−0.40 (−0.48, −0.33)	<0.001
Biomass as primary source of fuel	−0.22 (−0.30, −0.13)	<0.001	−0.17 (−0.26, −0.08)	<0.001	−0.18 (−0.27, −0.09)	<0.001

Abbreviation: SES, socioeconomic status.

Figure 5

Absolute difference in pre-bronchodilator FEV1/FVC Z-score by socioeconomic status composite index (SES index), secondary education or higher (below secondary as reference), greater than or equal to household size of 4 (below as reference), and monthly household income (change per category).

Notes: We display adjusted differences in pre-bronchodilator FEV1/FVC Z-score for the overall sample (left) and for individual sites (right). In the left panel, the square represents the adjusted odds ratio and the segments represent 95% CIs. In the right panel, the triangles represent site-specific adjusted odds ratios. For the composite SES index, we calculate the interquartile adjusted difference in pre-bronchodilator FEV1/FVC Z-score, ie, the difference between the 75th and 25th percentiles of the composite SES index.

Abbreviation: SES, socioeconomic status.

Discussion

In this analysis of multicenter population-based studies, we describe the association between SES and COPD among ten sites in Argentina, Bangladesh, Chile, Peru, and Uruguay. This study is among the largest to examine the association between SES and COPD in LMIC settings, pooling data from multiple, population-based studies. We found that lower education, lower household income, and a lower composite SES index were all associated with higher odds of having COPD and a lower pre-bronchodilator FEV1/FVC ratio. As household monthly income may be underestimated in study populations, here we describe a composite index that may provide a better surrogate for SES in resource-poor settings. Although the association between SES and COPD has been well described in HIC, there is value in population-based studies among LMICs to better understand these associations across income settings, identify populations at risk, and examine SES-attributable risk factors such as biomass fuel smoke exposure.8,12,25,26 Moreover, while some aspects of SES cannot be modified, eliminating SES disparities could not only decrease future COPD prevalence but may also have a greater impact on improving COPD outcomes than would current medications.8,27

Our study found similar rates of COPD prevalence (10.9% with a range from 3.5% to 18.0% among settings when using the fixed ratio) compared with other multiregion, population-based studies.1,5 The PLATINO study, conducted in five urban Latin American cities, had an overall prevalence of 14.5% with a range from 7.8% to 19.7%, while the BOLD multicenter study found an overall prevalence of 10% with a range from 8.5% to 22.2%.1,5 Our study also found that SES was inversely related to COPD, consistent with North American and European studies that used a range of measures (education, income, and deprivation scores) to evaluate SES among high-income populations. Specifically, we found a 25% increase in the odds of COPD between low- and high-SES quartiles, a modest increase compared with a previous review examining COPD and SES across income settings.8 This effect was mostly consistent across sites. Additionally, we found a positive relationship between composite SES index and pre-bronchodilator FEV1/FVC Z-score (ie, participants with higher index score had a greater FEV1/FVC for their age and sex), and these effects were consistent across all settings (Figure 5). This measure adds strength to the association between SES and COPD since FEV1/FVC Z-score is a continuous metric that does not need to be dichotomized such as lower limit of normal cutoffs.28 While previous studies have shown a similar relationship, this analysis is unique in that we use a composite index combining multiple SES predictors after controlling for important environmental exposures in a large population-based sample of individuals from LMICs.10 Since SES is a complex measure of a combination of variables, single variables such as education and household income may not have the discriminatory power to make accurate predictions.

Our analysis also has some potential shortcomings. First, we were unable to adjust for occupation, which could influence the association between SES and COPD. While income level and certain occupations correspond, some occupations may pose a risk to getting COPD that is not adjusted for when excluded. Second, monthly household income was collected as six categories, which may not be perfectly standardized and comparable among countries. Third, the GLI reference population used to obtain Z-scores in COPD diagnosis may not perfectly represent the populations used in this study. However, models were run using Caucasian values that showed almost identical results. While this does point to a robust method of diagnosis, further research is needed to investigate variation in results when using other GLI reference groups. Fourth, we used self-report of daily biomass fuel smoke exposure and did not measure personal or household pollutant exposures in all participating households. However, in previous research, we have found that self-report of daily biomass fuel smoke exposure was a good proxy for high vs low levels of exposure.29,30

The large sample size, across ten diverse sites including ranges of altitude and degree of urbanization, is a fundamental strength of this study. We used the same definitions for SES variables, used similar spirometer devices from the same manufacturer (ndd Medical Technologies), all studies were population-based with similar methods of screening, and data collections occurred during concurrent time periods.16–18 Furthermore, we used several SES indicators to establish a composite, more robust SES index. Combining household size, education, and household monthly income into a composite SES index likely accounted for income variability across sites. We controlled for daily smoking and biomass fuel smoke exposure, showing that SES is associated with having COPD above and beyond important environmental exposures. Finally, the study utilized lower limit of normal to diagnose COPD at each site as opposed to a fixed FEV1/FVC ratio utilized in previous studies. Utilization of lower limit of normal likely prevents overdiagnosis of COPD among the elderly where natural lung decline may be confounded with obstructive lung disease.1,31

Through analysis of a large pooled database from multiple LMICs, we found that SES was associated with COPD and lower lung function. Specifically, this study shows that participants of lower SES levels were more likely to have COPD and lower FEV1/FVC values for their age, using the lower limit of normal, which is currently recommended by international respiratory societies. The implications of the link between SES and disease are broader beyond COPD.32 For example, recently biochemical and genetic pathways have been identified showing that SES stressors result in DNA methylation and molecular genetic mechanisms, which in turn result in a host of inflammatory and chronic diseases.33,34 Although SES consists of a range of environmental exposures, this evidence suggests that the association with disease may be mediated by social and psychological stress. Our findings validate prior assumptions and findings that point to a relationship between SES and COPD and support efforts toward LMIC disease burden and priority setting in the global health agendas. Government policies aimed at disparities in SES among high-risk groups, such as campaigns to use cleaner fuels and guidelines for occupational safety, can occur at systems-based level. At an individual level, microenterprise and social entrepreneurship may serve as additional public health interventions beyond tobacco cessation. While global inequality has greatly reduced over the past 30 years, inequity within many LMICs is on the rise.35 SES will continue to be a risk factor, which must be addressed in any public health effort to decrease the burden of COPD.

Conclusion

In this analysis of multiple population-based studies, we found that lower SES was associated with higher odds of having COPD. Since income may be underestimated in study populations, composite scores including household size and education may be better proxies of SES for our analyses. Overall, the relationship between COPD and SES across five LMICs with varying disease burden was highlighted to guide future investigations and implement interventions to reduce the burden of COPD.