Socioeconomic Status and Longitudinal Lung Function of Healthy Mexican Children

Abstract

Introduction

Our aim was to estimate the longitudinal effect of Socioeconomic status (SES) on lung function growth of Mexican children and adolescents.

Materials and Methods

A cohort of Mexican children in third grade of primary school was followed with spirometry twice a year for 6 years through secondary school. Multilevel mixed-effects lineal models were fitted for the spirometric variables of 2,641 respiratory-healthy Mexican children. Monthly family income (in 2002 U.S. dollars [USD]) and parents’ years completed at school were used as proxies of SES.

Results

Individuals with higher SES tended to have greater height for age, and smaller sitting height/standing height and crude lung function. For each 1-year increase of parents’ schooling, Forced expiratory volume in 1 sec (FEV₁) and Forced vital capacity (FVC) increased 8.5 (0.4%) and 10.6 mL (0.4%), respectively (p <0.05) when models were adjusted for gender. Impact of education on lung function was reduced drastically or abolished on adjusting by anthropometric variables and ozone.

Conclusions

Higher parental schooling and higher monthly family income were associated with higher lung function in healthy Mexican children, with the majority of the effect likely due to the increase in height-for-age.

About the Authors:

David Martínez-Briseño

Affiliation: Epidemiology and Social Science in Health, National Institute of Respiratory Diseases (INER), Mexico City, Mexico

Rosario Fernández-Plata

Affiliation: Epidemiology and Social Science in Health, National Institute of Respiratory Diseases (INER), Mexico City, Mexico

Laura Gochicoa-Rangel

Affiliation: Epidemiology and Social Science in Health, National Institute of Respiratory Diseases (INER), Mexico City, Mexico

Luis Torre-Bouscoulet

Affiliation: Epidemiology and Social Science in Health, National Institute of Respiratory Diseases (INER), Mexico City, Mexico

Rosalba Rojas-Martínez

Affiliation: Population Health Research Center, National Institute of Public Health (INSP), Mexico City, Mexico

Laura Mendoza-Alvarado

Affiliation: Population Health Research Center, National Institute of Public Health (INSP), Mexico City, Mexico

Cecilia García-Sancho

Affiliation: Epidemiology and Social Science in Health, National Institute of Respiratory Diseases (INER), Mexico City, Mexico

Rogelio Pérez-Padilla

* E-mail: [email protected]

Affiliation: Epidemiology and Social Science in Health, National Institute of Respiratory Diseases (INER), Mexico City, Mexico

Introduction

Lung function, an important component of the evaluation of children with respiratory problems, is influenced by gender, height, and age, but also by prenatal exposures, genetic factors, ethnicity, obesity, altitude of residence, tobacco smoking, air pollution, nutrition, Socioeconomic status (SES), and lung disease [1]. Data from a variety of studies consistently show that low SES is associated with lower mean levels of pulmonary function [2–5]. However, data of children from developing countries is scarce, especially if studied in longitudinal fashion [6]. The pattern of increase in lung function may differ if obtained from longitudinal or cross-sectional studies [7] because, in the latter, the effects of age on lung function are confounded by secular time and by the presence of the multiple birth cohorts assembled in a cross-sectional study. A longitudinal study may better describe the growth spurt in adolescents [8] and the time sequence of events, such as the impact of general state of health and nutrition and SES on lung growth. The main objective of the present study was to evaluate the effect of SES on the lung function of healthy Mexican children from a cohort assembled in Mexico City.

Materials and Methods

Population and study design

The Metropolitan Study to Evaluate the Chronic Effects of Pollution in School-age Children (EMPECE) was undertaken in Mexico City beginning on April 23, 1996, with children in third grade of primary school [9]. The protocol was approved by the Ethics Committee of the Mexican National Institute of Respiratory Diseases (INER). All parents provided written informed consent for the study subjects.

We selected 10 fixed-site air-monitoring stations in Mexico City and randomly selected 39 elementary schools from among those located within a 2-km radius of the stations. The study cohort consisted of students at the selected schools who were 8 years of age at the beginning of the study, who had not been diagnosed as having asthma, and whose parents had signed a letter of informed consent. A substantial number of children entered or left the cohort during the course of the study. At baseline, a questionnaire was completed by the parents of 1,819 children, and a spirometric test was administered to each child (phase 1). The following evaluations occurred every 6 months during the spring and autumn seasons of each year until the end of the children’s primary school education in 1999. Children remaining in the same schools studied during secondary school were followed for 3 additional years until 2002; thus, we collected information for a total of 12 evaluations [9–11].

Spirometry testing

Spirometry tests were conducted using identical computerized dry-rolling seal spirometers (922 Spirometer; SensorMedics, USA) that were calibrated each morning prior to data collection with a 3-L syringe (SensorMedics). We recorded only the expiratory part of Forced expiratory maneuvers (FEM) and analyzed Forced expiratory volume in 1 sec (FEV₁), Forced vital capacity (FVC), and their ratio (FEV₁/FVC). Tests were performed at the school during morning and early afternoon hours. As many as eight FEM were conducted for each child to obtain three acceptable ones according to 1994 American Thoracic Society (ATS) criteria [12]. Additional details on the spirometry methodology, including sustained quality control along the study, were described in a previous report [10].

Respiratory-healthy children

Children with self-report of asthma, tobacco smoking, chronic respiratory symptoms (cough, wheezing, phlegm, dyspnea), or children in the >95% percentile of Body mass index (BMI) for age according to growth charts from the Centers for Disease Control and Prevention (CDC) [13], were excluded from the analysis because those conditions would tend to reduce lung function and confound a relationship with SES. Children <8 years of age were scarce and were also excluded, as performance of spirometry is less reliable than in older children. Although some of these diseases, symptoms, and exposures could also be associated with poverty, we preferred to analyze the impact of SES on lung function in a healthy population from the respiratory point of view, but also included statistical models with data from the whole cohort.

Socioeconomic status

Socioeconomic status was based on monthly family income (in 2002, log-transformed U.S. dollars [USD]) assessed at each evaluation and by parents’ schooling (years completed at school), both reported by the parents and categorized in tertiles for analysis. Parents’ schooling was built taking the average of the highest completed school years of both parents. Income changed over time, while parents’ schooling was constant.

Ozone measurement

The impact of several air pollutants on the lung function of the cohort was previously reported (9). As an indicator of air pollution, we included ambient O₃ from 10 government air-monitoring stations, assigning for each child data from the station closest to their school. We calculated 8-h ozone means (parts per billion [ppb] between 10 a.m. and 6 p.m.), whereas long-term exposure for each day of the study period was estimated as the average over the previous 6 months of the daily O₃ 8-h mean [9].

Statistical analysis

We estimated tertiles of schooling and compared spirometric and anthropometric variables by Analysis of variance (ANOVA). We then fitted multilevel mixed-effects linear models adjusted for age and gender to determine the association of SES (monthly family income or parents’ years at school) with spirometric variables (log-transformed FEV₁ and FVC) along time with and without height, ozone (average over the previous 6 months of the daily O₃ 8-h mean, as a general indicator of air pollution), interaction between SES variables and ozone and secondhand smoke (the presence of any person smoking inside the home). We also fitted multilevel mixed-effects linear models for spirometric variables expressed as Z-score [14] including SES variables, second hand smoke and ozone. Sitting height was measured only during the 4th evaluation; therefore, we could include sitting height/standing height as additional confounder and also as a consequence of SES only for observations during this phase. We also estimated the height-for-age Z-score [13] and the sitting height/standing height ratio (only available in the 4th evaluation) in order to compare groups (by tertiles), and we established a relationship between these two variables with spirometric and other anthropometric measurements.

The analysis was conducted using the software STATA v.13 software (Stata Corp., College Station, TX, USA).

Results

Of the 3,177 children finally included in the cohort, 536 presented at least one of the exclusion criteria as follows: 28 were <8 years of age; 67 reported having asthma; 190 were smokers, and 251 had obesity. Thus, we collected a total of 14,165 measurements from 2,641 children who were considered respiratory-healthy and who were 8–17 years of age. Mean age at study inclusion was 9.2 years of age (Standard deviation [SD], 1.1 years) for girls and 9.5 years (SD, 1.2 years) for boys. Observations per individual ranged in number from 1‒12 (median, 4 observations, and interquartile range, 4 observations). Mean follow-up duration was 2.5 years (SD, 1.9 years). Study participants’ characteristics by follow-up phase can be observed in (S1 and S2 Tables) [14].

Tables 1 and 2 describe socioeconomic variables and ozone exposure by follow-up phase and gender. Response rate (RR) of the monthly family income variable was 84.5% along the follow-up period. Table 3 depicts the characteristics of the study population in phase 4, the only phase including measurements of sitting height/standing height. Table 4 presents results by tertiles of parents’ schooling, while S3 Table separates data by tertiles of monthly family income. The lowest tertile had a reported mean monthly family income of $119.4 USD and 7.7 years at school, the middle tertile, $360.6 USD and 11.8 years at school, and the highest tertile, $949.3 USD and 16.0 years at school.

[Figure omitted. See PDF.]

Table 1. Socioeconomics status (SES) and ozone exposure of the boys studied (means and Standard deviation [SD]).

https://doi.org/10.1371/journal.pone.0136935.t001

[Figure omitted. See PDF.]

Table 2. Socioeconomics status (SES) and ozone exposure of the girls studied (means and Standard deviation [SD]).

https://doi.org/10.1371/journal.pone.0136935.t002

[Figure omitted. See PDF.]

Table 3. Characteristics of the studied population in phase 4 of the study.

SD = Standard deviation; BMI = Body mass index; FEV₁ = Forced expiratory volume at 1 sec; FVC = Forced vital capacity; FEV₁/FVC = ratio of FEV₁ and FVC.

https://doi.org/10.1371/journal.pone.0136935.t003

[Figure omitted. See PDF.]

Table 4. Characteristics of the children studied by tertile of parents’ schooling.

https://doi.org/10.1371/journal.pone.0136935.t004

Children of families with higher parental education in general were taller for age, had higher crude spirometric variables and, during the 4th evaluation, also had lower sitting height/standing height. Similarly, poorer children had smaller body and lung size, but in addition, smaller legs in relation to trunk size (Table 4).

Mixed regression models predicting log FEV₁ (mL) and log FVC (mL) showed that the higher the monthly family income, the higher the lung function (models 1‒4). Parents’ schooling also was associated with higher lung function (model 1). When monthly income increased by 10%, FEV₁ increased 0.008‒0.03% and FVC, 0.02‒0.04% (models 1‒4). However, only the co-efficient in model 1 was statistically significant for both FEV₁ and FVC. A 1-year increase in parents’ schooling was associated with an increase of 0.37% for FEV₁ and 0.41% for FVC (Table 5). Utilizing model 1 of Table 5 and average income ($300.40 USD), schooling (11.5 years), FEV₁ (2,322 mL), and FVC (2,605 mL), an increment of $100 USD in monthly income increased FEV₁ by 1.98 mL (0.09%) and FVC by 2.88 mL (0.11%), whereas a 1-year increment in parents’ schooling increased FEV₁ by 8.5 mL (0.37%) and FVC by 10.6 mL (0.41%).

[Figure omitted. See PDF.]

Table 5. Longitudinal models for spirometric variables and Socioeconomic status (SES), both genders taken together.

https://doi.org/10.1371/journal.pone.0136935.t005

If anthropometry is added to the model (model 2, Table 5), the impact of SES on lung function became non-significant. Secondhand smoke and ozone levels also exerted adverse effects on lung function, but had no additional impact on the relationship between SES and lung function (model 4). Higher ozone levels were observed in areas with students with higher SES.

Mixed regression models by gender or including healthy and non-healthy children were also performed, observing co-efficients for SES variables were similar to those obtained by models with healthy children and including both genders (S4–S7 Tables). Finally, cross-sectional models obtained from the first observation of each individual were compared with the more adequate longitudinal mixed-regression models illustrated in Table 5, observing similar effects for parents’ schooling, but not for monthly family income compared with longitudinal models. Furthermore, interaction between SES variables and ozone were not statistically significant. Longitudinal models had a lower Akaike information criterion (AIC) than cross-sectional models (S8 Table), again showing better fit of the observations.

S9 Table shows the effect of SES over spirometric variables expressed as Z-score and confirmed the positive effect of monthly family income on FEV₁ (P<0.1) and FVC (P<0.01). However parents’ schooling co-efficient was not statistically significant.

Fig 1 depicts spirometric variables (L), sitting height/standing height ratio (%, from the 4th evaluation), and monthly family income (2009 USD) of the cohort vs. height-for-age Z-score. As height-for-age increases, spirometric variables and income also increase, but sitting height/standing height decreases. Fig 2 shows crude FEV₁ and FVC in L vs. age. Children in the uppermost tertile of monthly family income had a slight higher lung function over time than children in the first two tertiles.

[Figure omitted. See PDF.]

Fig 1. Relationship between height-for-age Z-score and spirometric, sitting height/standing height, and income.

(A) Lung function (FVC and FEV₁), (B) sitting height/standing height, and (C) monthly family income as a function of height for age as Z-score (horizontal axis). Line and shadows represent linear or quadratic regression and 95% Confidence intervals (95% CI) of the regression models. Symbols are means of terciles for height-for-age as described in the text, showing good fit of the models.

https://doi.org/10.1371/journal.pone.0136935.g001

[Figure omitted. See PDF.]

Fig 2. Forced expiratory volume at 1 sec (FEV₁) and Forced vital capacity (FVC) expected as a function of age by tertiles of monthly family income.

(A) FEV₁ and (B) FVC as a function of age (horizontal axis). Lung function increases with growth and development in alinear fashion vs. age, but also varies with Socioeconomic status (SES) (shown as terciles in different lines), especially at the end of the observation.

https://doi.org/10.1371/journal.pone.0136935.g002

Discussion

Our study shows, in a cohort of school age children residing 2,240 m above sea level, a lower FEV₁ and FVC in families with lower monthly income or lower parents’ schooling. Socioeconomic disparities are associated with differences in children’s lung function. The majority of studies to date, many of these cross-sectional, demonstrate a positive correlation between SES (income, education, occupation, and others, or a combination of these in an index) and lung function [15], as our results describe.

The detailed mechanisms of why SES adversely impacts lung function are unknown, but are likely multifactorial and beyond the common explanation of poor nutrition [3]. SES adversely impacts growth and development; therefore, several anthropometric measurements, such as height-for-age and sitting height/standing height [6, 16–19], are important predictors of body size and overall lung size and function, and are likely key intermediate variables between SES and lung function. In fact, in our study, the impact of SES on lung function disappeared on adjusting regression models by anthropometry, observing an expected adverse effect on lung function and air pollutants (exemplified by ozone levels) and passive smoking, supporting the idea that the majority of the effect of SES on lung function is associated with changes in height and sitting height/standing height. In other words, poverty leads to smaller body size, lung size, and crude lung function, but not to height- and age-adjusted lung function. Changes in height-adjusted lung function, for example, when function is expressed as percentage predicted by gender, age, and height, or as Z- score of gender-, age-, and height-predicted values, were mild and inconsistent in different tertiles of income or education.

Compared with children in developed countries, including Mexican-American children residing close to sea level, gender-, age-, and height-adjusted lung function of children from Mexico City have been higher [14] suggesting that an adverse environmental impact on height is more marked than on the lung, especially in individuals residing at moderate or high altitudes. However, crude lung function in Mexican-American children is clearly higher than in children of the same age residing in Mexico City, as also shown in the groups with higher SES within Mexico City; improving nutrition during growth and development could raise the maximal lung function attained and decrease the risk of adult-age lung diseases such as Chronic obstructive pulmonary disease (COPD) [20, 21]. Another source of the adverse impact of low SES on lung function derives from air pollution, which tends to be worse in poorer areas. For example, solid fuel use, very much a pollutant, is linked with poverty; however, in our study, the use of solid fuels is very uncommon and ozone levels were higher in more highly advantaged areas. In addition, previous studies have found a differential impact of air pollutants depending on SES, but not consistently [22]. For example benzene, NO₂, and NOx [21], and Ozone can be modified by maternal educational level or SES, respectively, but particulate matter and NO₂ have a small decrease in lung function, even when SES was taken into account in the models. [23].

Several study limitations should be acknowledged. Income has disadvantages on measuring SES, because it is often under-reported, although income may change over time, reflecting better modifications in family economy during follow-up. Schooling, on the other hand, does not change over time, even with changes in income or poverty, but is able to be measured with greater reliability. More complex indicators of SES have been proposed, with composite indices including existent in-home appliances, but these are not necessarily better than traditional family income or education.

Longitudinal follow-up of the cohort ended when the children were about 15 years of age, during adolescence, but before final height and lung function were reached. In addition, as can be observed in Tables 1 and 2, follow-up decreased during secondary school due to the children’s dropping out of school or changing to a different school, reducing information at the beginning of adolescence. We also lacked contributions of genetic ancestry, because biological samples were not obtained and questions on race or ethnic origin are usually avoided in Mexico, as the majority of the population is considered Mexican mestizo and differs from the so-called Hispanic population in the U.S., with a variety of contributions of African ancestry depending on country-of-origin. On the other hand, the study includes a considerable number of spirometric observations conducted with good quality control, as well as measurements of air pollution and a heterogeneous socioeconomic status of the individuals studied, allowing for relevant analysis.

Conclusions

Lower socioeconomic status, estimated as reported income or parents’ education, is associated with a reduction in lung function and height-for-age and higher sitting height/standing height. According to our data, the majority of the impact of SES on lung function was accomplished by means of a reduction in height. Further studies are needed in order to explore the potential factors underling this association and the impact of residing at a moderate altitude.

Supporting Information

[Figure omitted. See PDF.]

S1 Table. Main characteristics of boys studied (means and Standard deviation [SD]).

Population characteristics.

https://doi.org/10.1371/journal.pone.0136935.s001

(DOC)

S2 Table. Main characteristics of girls studied (means and Standard deviation [SD]).

Population characteristics.

https://doi.org/10.1371/journal.pone.0136935.s002

(DOC)

S3 Table. Anthropometric and spirometric variables by tertiles of monthly family income.

Population characteristics.

https://doi.org/10.1371/journal.pone.0136935.s003

(DOC)

S4 Table. Longitudinal models for Socioeconomic status (SES) and lung function in girls.

https://doi.org/10.1371/journal.pone.0136935.s004

(DOC)

S5 Table. Longitudinal models for Socioeconomic status (SES) and lung function in boys.

https://doi.org/10.1371/journal.pone.0136935.s005

(DOC)

S6 Table. Longitudinal models for Socioeconomic status (SES) and lung function including unhealthy girls.

https://doi.org/10.1371/journal.pone.0136935.s006

(DOC)

S7 Table. Longitudinal models for Socioeconomic status (SES) and lung function including unhealthy boys.

https://doi.org/10.1371/journal.pone.0136935.s007

(DOC)

S8 Table. Cross-sectional models for spirometric variables and Socioeconomic status (SES), both genders.

https://doi.org/10.1371/journal.pone.0136935.s008

(DOC)

S9 Table. Longitudinal models for spirometric variables expressed as Z-score and Socioeconomic status (SES), both genders taken together.

https://doi.org/10.1371/journal.pone.0136935.s009

(DOC)

S1 Dataset. Compressed Dataset in Excel (XLS) format.

https://doi.org/10.1371/journal.pone.0136935.s010

(XLS)

Acknowledgments

We are grateful to all of the participants in the EMPECE study and to their parents and teachers.

Author Contributions

Conceived and designed the experiments: RPP RRM LTB LGR. Performed the experiments: RRM LMA RPP. Analyzed the data: DMB RFP CGS. Contributed reagents/materials/analysis tools: RPP LTB LGR DMB RFP CGS. Wrote the paper: RPP LTB LGR DMB RFP CGS RRM LMA.

Citation: Martínez-Briseño D, Fernández-Plata R, Gochicoa-Rangel L, Torre-Bouscoulet L, Rojas-Martínez R, Mendoza-Alvarado L, et al. (2015) Socioeconomic Status and Longitudinal Lung Function of Healthy Mexican Children. PLoS ONE 10(9): e0136935. https://doi.org/10.1371/journal.pone.0136935

References

1. Lum S, Stocks J. Forced expiratory manoeuvres. Eur Respir Mon. 2010;47:46–65.

2. Steinberg M, Becklake MR. Socio-environmental factors and lung function. S Afr Med J. 1986;70(5):270–4. pmid:3529452.

3. Hegewald MJ, Crapo RO. Socioeconomic status and lung function. Chest. 2007;132(5):1608–14. pmid:17998360.

4. Demissie K, Ernst P, Hanley JA, Locher U, Menzies D, Becklake MR. Socioeconomic status and lung function among primary school children in Canada. Am J Respir Crit Care Med. 1996;153(2):719–23. pmid:8564123.

5. Raju PS, Prasad KV, Ramana YV, Balakrishna N, Murthy KJ. Influence of socioeconomic status on lung function and prediction equations in Indian children. Pediatr Pulmonol. 2005;39(6):528–36. pmid:15789442.

6. Menezes AM, Dumith SC, Perez-Padilla R, Noal RB, Wehrmeister FC, Martinez-Mesa J, et al. Socioeconomic trajectory from birth to adolescence and lung function: prospective birth cohort study. BMC public health. 2011;11:596. pmid:21794105; PubMed Central PMCID: PMC3155502.

7. Lebowitz MD. Age, period, and cohort effects. Influences on differences between cross-sectional and longitudinal pulmonary function results. Am J Respir Crit Care Med. 1996;154(6 Pt 2):S273–7. Epub 1996/12/01. pmid:8970400.

8. Ware JH, Weiss S. Statistical issues in longitudinal research on respiratory health. Am J Respir Crit Care Med. 1996;154(6 Pt 2):S212–6. Epub 1996/12/01. pmid:8970390.

9. Rojas-Martinez R, Perez-Padilla R, Olaiz-Fernandez G, Mendoza-Alvarado L, Moreno-Macias H, Fortoul T, et al. Lung function growth in children with long-term exposure to air pollutants in Mexico City. Am J Respir Crit Care Med. 2007;176(4):377–84. Epub 2007/04/21. pmid:17446338.

10. Perez-Padilla R, Regalado-Pineda J, Mendoza L, Rojas R, Torres V, Borja-Aburto V, et al. Spirometric variability in a longitudinal study of school-age children. Chest. 2003;123(4):1090–5. pmid:12684298

11. Pérez‐Padilla R, Regalado‐Pineda J, Rojas M, Catalán M, Mendoza L, Rojas R, et al. Spirometric function in children of Mexico City compared to Mexican‐American children. Pediatric pulmonology. 2003;35(3):177–83. pmid:12567385

12. Standardization of Spirometry, 1994 Update. American Thoracic Society. Am J Respir Crit Care Med. 1995;152(3):1107–36. pmid:7663792.

13. Kuczmarski RJ, Ogden CL, Guo SS, Grummer-Strawn LM, Flegal KM, Mei Z, et al. 2000 CDC Growth Charts for the United States: methods and development. Vital and health statistics Series 11, Data from the national health survey. 2002;(246):1–190. Epub 2002/06/05. pmid:12043359.

14. Martinez-Briseno D, Fernandez-Plata R, Gochicoa-Rangel L, Torre-Bouscoulet L, Rojas-Martinez R, Mendoza L, et al. Longitudinal lung function growth of Mexican children compared with international studies. PLoS One. 2013;8(10):e77403. pmid:24143231; PubMed Central PMCID: PMC3797091.

15. Krieger N. A glossary for social epidemiology. Journal of epidemiology and community health. 2001;55(10):693–700. pmid:11553651; PubMed Central PMCID: PMC1731785.

16. Bogin B, Varela-Silva MI. Leg length, proportion, health and beauty: a review. Anthropologischer Anzeiger; Bericht uber die biologisch-anthropologische Literatur. 2009;67(4):439–59. Epub 2010/05/06. pmid:20440962.

17. Bogin B, Varela-Silva MI. Leg length, body proportion, and health: a review with a note on beauty. International journal of environmental research and public health. 2010;7(3):1047–75. Epub 2010/07/10. pmid:20617018; PubMed Central PMCID: PMC2872302.

18. Malina RM, Selby HA, Buschang PH, Aronson WL. Growth status of schoolchildren in a rural Zapotec community in the Valley of Oaxaca, Mexico, in 1968 and 1978. Annals of human biology. 1980;7(4):367–74. Epub 1980/07/01. pmid:7436350.

19. Malina RM, Reyes ME, Little BB. Socioeconomic variation in the growth status of urban school children 6–13 years in Oaxaca, Mexico, in 1972 and 2000. American journal of human biology: the official journal of the Human Biology Council. 2009;21(6):805–16. Epub 2009/04/23. pmid:19384866.

20. Stocks J, Hislop A, Sonnappa S. Early lung development: lifelong effect on respiratory health and disease. The Lancet Respiratory medicine. 2013;1(9):728–42. pmid:24429276.

21. Stocks J, Sonnappa S. Early life influences on the development of chronic obstructive pulmonary disease. Therapeutic advances in respiratory disease. 2013;7(3):161–73. pmid:23439689; PubMed Central PMCID: PMC4107852.

22. Sacks JD, Stanek LW, Luben TJ, Johns DO, Buckley BJ, Brown JS, et al. Particulate matter-induced health effects: who is susceptible? Environ Health Perspect. 2011;119(4):446–54. pmid:20961824; PubMed Central PMCID: PMC3080924.

23. Gehring U, Gruzieva O, Agius RM, Beelen R, Custovic A, Cyrys J, et al. Air pollution exposure and lung function in children: the ESCAPE project. Environ Health Perspect. 2013;121(11–12):1357–64. pmid:24076757; PubMed Central PMCID: PMC3855518.