Abstract
Introduction:
Nitrogen dioxide (NO2) is known to be a trigger for asthma exacerbation. However, little is known about the role of seasonal variation in indoor and outdoor NO2 levels in childhood asthma in a mixed rural-urban setting of North America.
Methods:
This prospective cohort study, as a feasibility study, included 62 families with children (5-17 years) that had diagnosed persistent asthma residing in Olmsted County, Minnesota. Indoor and outdoor NO2 concentrations were measured using passive air samples over 2 weeks in winter and 2 weeks in summer. We assessed seasonal variation in NO2 levels in urban and rural residential areas and the association with asthma control status collected from participants’ asthma diaries during the study period.
Results:
Outdoor NO2 levels were lower (median: 2.4 parts per billion (ppb) in summer, 3.9 ppb in winter) than the Environmental Protection Agency (EPA) annual standard (53 ppb). In winter, a higher level of outdoor NO2 was significantly associated with urban residential living area (P = .014) and lower socioeconomic status (SES) (P = .027). For both seasons, indoor NO2 was significantly higher (P < .05) in rural versus urban areas and in homes with gas versus electric stoves (P < .05). Asthma control status was not associated with level of indoor or outdoor NO2 in this cohort.
Conclusions:
NO2 levels were low in this mixed rural-urban community and not associated with asthma control status in this small feasibility study. Further research with a larger sample size is warranted for defining a lower threshold of NO2 concentration with health effect on asthma in mixed rural-urban settings.
Introduction
Evidence indicates that exposure to nitrogen dioxide (NO2) is significantly associated with asthma incidence and poorer health outcomes in the United States (U.S.) and elsewhere.1-8 Prospective studies examining the association between NO2 and childhood asthma have been conducted largely in urban communities where traffic volume, level of NO2, and asthma prevalence are higher than in rural or mixed rural-urban areas.3,4,9 Indoor and outdoor NO2 exposure levels were assessed among children with asthma living in urban and suburban Connecticut and Massachusetts, and significant associations were found between increased levels of indoor NO2 (above 6 parts per billion (ppb)) and asthma symptoms and severity, at levels well below the U.S. Environmental Protection Agency (EPA) outdoor standards (53 ppb). 10 These observations indicate an NO2 concentration below the EPA standard might still pose a significant health effect to people with asthma, which is yet unknown whether this is the case for non-urban settings.
Olmsted County, Minnesota (MN), was rated as having excellent air quality by the American Lung Association in 2014 based on the level of ambient ozone. 11 However, the prevalence of current asthma in school age children is 13%, which exceeds the national average (8%).12,13 Given the high prevalence of asthma in this mixed rural-urban area, the reported health effect of low NO2 concentration on asthma, and a relatively long and cold winter season in our study setting, it is worthy of exploring factors associated with residential NO2 concentrations (eg, indoors from gas stove or oven, unvented gas fireplaces; and outdoors from traffic and combustion from electric power plants) and the health effect of exposure to NO2. 14 Measuring and assessing the health effect of indoor and outdoor NO2 concentrations in a mixed rural-urban setting in a Midwest region may help us better define a lower threshold of NO2 concentrations with health effect.15,16 Further, given increasing traffic volume and population within Olmsted County, associated with the multi-year $5.6 billion economic development initiative, it would be worthy of assessing the health effect of indoor and outdoor NO2 concentrations. 17
Our aim was to test if it is feasible to measure residential levels of outdoor and indoor NO2, describe residential characteristics related to spatial and seasonal variations in NO2, and explore associations of residential NO2 exposure with respiratory outcomes among a group of pediatric patients with persistent asthma.
Methods
Study Setting and Population
Olmsted County, MN (area: 1700 km2, population: 156,277 in 2018) contains both urban and rural areas as defined by the United States Census (16% rural population and 91% rural area), 18 and it is demographically similar to the surrounding Midwest region, with the exception that a large proportion of the population is employed in the health care field.19-21 Olmsted County currently has a single monitoring site which measures PM2.5 and ozone concentration but not NO2. Previous studies showed that the prevalence of asthma among children and adolescents in Olmsted County, MN, in 2000 was higher than national average (13%vs 8% at the national level) and asthma was the third highest health care expenditure in children and adolescents in Olmsted County, MN.12,13,22
Study Design
This is a single site prospective cohort study conducted at the Mayo Clinic with a collaboration with the Center for Perinatal, Pediatric, and Environmental Epidemiology, Yale School of Public Health. Informed consent was obtained from Mayo Clinic pediatric primary care patients with persistent asthma residing in Olmsted County (n = 62). We assessed their indoor and outdoor NO2 levels for a total of 4 weeks (2 in winter and 2 in summer) using passive monitors (Palmes tube),23,24 (Figure 1) and collected asthma outcomes during the same period. This study was approved by the Institutional Review Board of Mayo Clinic and the Human Investigation Committee of Yale University.

Photograph of an NO2 monitor (Palmes tube) clipped to blue plastic weather protector jar and suspended from a tree branch in the backyard of a home. Unexposed end of the Palmes tube within the protector jar is sealed (green cap with white sample label), and screen lid is replaced and screwed on tightly.
Study Cohort
All eligible study cohort participants were Mayo Clinic pediatric primary care patients aged 5 to 17 years who (1) had at least one clinic visit to Mayo Clinic Rochester between 8/1/2015 and 1/31/2019, (2) had a documented physician diagnosis of persistent asthma from clinical notes, 25 (3) currently using controller medication(s) (eg, inhaled corticosteroid), and (4) resided at an address in Olmsted County, MN (n = 947). An invitation letter was sent to these families who, upon indicating interest in the study, were screened by phone interview to check inclusion and exclusion criteria. Exclusion criteria included those who did not have a protocol-acceptable setting for the monitor at least 10 feet away from their residence (for outdoor NO2), non-English speaking, and where the child with asthma was not in the home a minimum of 5 days per week, or had other chronic lung conditions (eg, cystic fibrosis).
Measurement of Indoor and Outdoor NO2 and Collection of Asthma Outcomes
We followed the protocol and methods for measuring indoor and outdoor NO2 previously established for environmental studies in an urban setting by the Yale University study team. 10 A study materials kit was mailed to each family who agreed to participate. Each participant received an initial kit for the winter season followed by a second kit for the summer measurement period. Each kit included 2 NO2 samplers (one for indoor placement and one for outdoors (Figure 1),23,24 a stand for hanging the indoor NO2 sampler, a small plastic open-ended jar to serve as a weather shield for the outdoor monitor, a set of installation instructions, and an asthma diary for recording asthma symptoms. A study coordinator placed a follow-up phone call to families to facilitate proper placement of the passive air samplers and conduct a structured interview collecting baseline information (eg, presence of gas stove). Following each 2-week monitoring period (2 weeks in winter and 2 weeks in summer, for a total of 4 weeks per residence to capture seasonal variability), the study coordinator contacted the families by phone to arrange for removal and return of the samplers. A 2-week integrated average NO2 (ppb) concentration was determined for each period, indoors and outdoors, for each residence.
Asthma Control Status
The parent or guardian of the participant was instructed to record the following information on the asthma diary: (1) asthma day and night symptoms (eg, wheezing, persistent cough, chest tightness, and limited activities), and (2) asthma rescue medication use (eg, albuterol). At the time of each NO2 monitor take-down interview they were asked to refer to the asthma diary to answer questions. With the collected data from participants, an asthma control status, based on the 2007 National Asthma Education and Prevention Program (NAEPP) guideline, was constructed for data analysis. 26 The score of asthma control status consisted of 3 components collected: symptom steps (day and night) and rescue medication steps. We defined the score of asthma control status as well-controlled (Score 1: ≤2 daytime asthma symptoms/week, ≤1 nighttime asthma symptom/month, or ≤2 rescue medication use/week), not well-controlled (Score 2: 2-6 daytime asthma symptoms/week, ≥2 nighttime asthma symptoms/month, or 2-6 rescue medication use/week), and poorly controlled (Score 3: 7 daytime asthma symptoms/week, ≥2 nighttime asthma symptoms/week, or 7 rescue medication use/week). If the patient has a “Score 3” for any category (daytime symptoms, nighttime symptoms, or rescue medication), then they are considered to have “poorly controlled” asthma. If a “Score 2” is the highest score among 3 categories, they are considered to have “not well-controlled” asthma). 27
Other Covariates
As part of the structured study interview conducted by phone, we collected information on second-hand smoke exposure, type and usage of cooking and heating system, presence of electronic air cleaner, presence of gas fireplace or logs, and mother’s educational level. In addition, we assessed the socioeconomic status (SES) of the study subjects using the HOUSES index.28 -33 Participants’ addresses were geocoded to link real property data to each housing unit (eg, living unit for apartment complex), and each property item corresponding to the individual’s address was standardized into a z-score within a county and converted to a quartile such that a higher HOUSES score indicated higher SES. Rural-urban classification of the residence was defined by the U.S. Census Bureau. 34 Average annual daily traffic (AADT) volume data was obtained from Minnesota Geospatial Commons. 35 The GIS package ArcMap (ESRI, Redlands, CA) was used to geocode the residential addresses.
Statistical Analysis
Descriptive statistics include counts (percentage) for categorical variables and medians (25th-75th percentile) for continuous variables (Table 1). Spearman’s correlation coefficients were calculated to measure the degree of correlation between indoor and outdoor NO2 measurements in each season (winter and summer). Univariate analyses were conducted using Kruskal Wallis tests to examine associations between each categorical baseline characteristic (eg, HOUSES) and level of NO2 (indoor and outdoor measurements separately). Univariate association between indoor/outdoor NO2 and asthma outcomes (binary variables; none vs any symptoms or use of rescue medication) were tested using logistic regression models stratified by season and the odds ratios (ORs) and 95% confidence intervals were presented. To test overall association between NO2 and asthma outcomes while considering dependency in data collected in 2 seasons for a given subject, a generalized estimating equation for logistic regression model was used (referred to as “Combined”). As exploratory analyses, we stratified by sex for assessing the association between NO2 levels and asthma outcomes, and also repeated analyses using dichotomized NO2 level (<6 ppb or ≥6 ppb) as an earlier study demonstrated significant associations between increased levels of indoor NO2 and asthma symptoms and severity beginning at levels above 6 ppb. 10 All analyses were performed using SAS statistical software package (version 9.4M6; SAS Institute, Inc, Cary, NC) and a P-value less than .05 was considered to be statistically significant.
Basic Characteristics of Study Subjects.
Greatest parental educational attainment by self-report (1 missing); bHOUSES: Individual-level socioeconomic status measure (8 missing); cAt least one episode for daytime asthma symptoms, nighttime asthma symptoms, and rescue medication usage; dPer National Asthma Education and Prevention Program EPR-3 (See the section 2.5 in detail).
Results
Basic Characteristics
Table 1 summarizes baseline characteristics of the study subjects. The median age was 11 (25th – 75th percentile: 8-14) years with 29 (47%) female, 54 (87%) non-Hispanic White, and 43 (69%) self-reporting any type of allergies. Thirty-two (52%) children had parents whose educational level was greater than a college education, and 26 (48%) families came from the highest SES (HOUSES Q4 in the community), with 9 (15%) families living in rural areas. All children were on asthma controller medications with 58 (94%) experiencing asthma symptoms in the 12 months prior to the enrollment date. None of the participants reported second-hand smoke exposure, and 18 (29%) had a gas stove at home. Overall poorer asthma outcomes (eg, more symptoms, asthma rescue medications) were reported in winter than in summer.
Spatial and Seasonal Variations in NO2 Related to Residential Characteristics
A total of 238 passive air samplers were returned with 2 unanalyzable samples excluded from NO2 analysis. Median level of NO2 in winter was 3.4 ppb (25th – 75th percentile: 2.5-5.9) for indoor and 3.9 ppb (2.8-5.0) for outdoor (P < .001), and 2.4 ppb (1.8-4.9) for indoor and 2.0 ppb (1.7-2.7) for outdoor in summer (P < .001), respectively. Table 2 summarizes seasonal variations in NO2 related to covariates measured with overall higher correlation between indoor and outdoor NO2 in summer than in winter, albeit statistically not significant (Spearman Correlation coefficient: .14 in summer vs −.03 in winter). Median level of NO2 in winter was 3.4 ppb (25th – 75th percentile: 2.5-5.9) for indoor and 3.9 ppb (2.8-5.0) for outdoor (P < .001), and in summer its was 2.4 ppb (1.8-4.9) for indoor and 2.0 ppb (1.7-2.7) for outdoor (P < .001).
Association of Level of Indoor and Outdoor NO2 With Covariates by Season (unit: ppb).
NO2 measured in ppb; Median (25th percentile, 75th percentile) of NO2 level is reported for each subgroup and subgroups are compared with Kruskal-Wallis test; Indoor features were compared only for Indoor NO2 levels.
Spearman correlation coefficient was used to assess correlation between indoor and outdoor NO2 levels in a given season.
P-value comparing between levels of indoor and outdoor NO2 in each season.
Lower SES as defined by HOUSES index was significantly associated with higher median outdoor NO2 in winter (5.2vs 4.4vs 4.3vs 3.2 for Q1 (lowest SES), Q2, Q3, and Q4, respectively; P = .027); this association was not found for indoor NO2 (P = .54). Levels of outdoor NO2 in winter (median outdoor NO2 (25th – 75th percentile)) in urban areas were significantly higher than in rural areas (4.0 (3.0-5.2) for urban vs 2.7 (2.3-3.7) for rural; P = .014) (Table 2 and Figure 2). This difference between urban and rural areas was not observed in the summer (2.0vs 1.8; P = .084). In winter, the median indoor NO2 level in the homes with gas stoves was higher than those without (7.2 (5.0-8.2) vs 2.9 (2.2-4.0); P < .001). The median level of indoor NO2 was higher in rural compared to urban homes for both seasons (6.7 (4.2-8.5) (rural) vs 3.2 (2.4-5.0) (urban) in winter (P = .022) and 5.7 (3.0-6.3) versus 2.3 (1.7-4.0) in summer (P = .030)). Presence of a gas fireplace in winter was not associated with increased indoor NO2 (P = .13).

Spatial distribution of residences classified by levels of observed outdoor NO2 in winter, and average annual daily traffic (AADT) on numbered roadways, Olmsted County, MN.
Associations of Asthma Outcomes With NO2 Level
In the univariate analysis, no association was found between NO2 levels (indoor or outdoor) and asthma control status, which was confirmed by repeated measures analysis in this cohort (Table 3). In our exploratory analysis where we further assessed if sex might modify this association, the association of daytime asthma symptoms with higher indoor NO2 levels during the summer was found among females (n = 29) but not males while poor asthma control status in winter was associated with indoor NO2 levels among males but not females (data not shown). Further exploratory analysis was conducted to evaluate the relationship between dichotomized NO2 level and asthma control status; the lack of association remained using dichotomized NO2 level (<6 ppb or ≥6 ppb) (data not shown).
Univariate Odds Ratios for Each Asthma Outcome by Level of NO2 Stratified by Season.
Logistic regression models using each asthma outcome as the dependent variable with reported odds ratio and 95% confidence intervals corresponding to a 1 -unit (ppb) change in NO2 level.
Per National Asthma Education and Prevention Program EPR-3 (see Method section in detail) with Well-Controlled as a reference.
Mixed effect logistic regression model treating subject as random effect.
Discussion
It is feasible to measure residential levels of outdoor and indoor NO2 to describe residential characteristics related to spatial and seasonal variations in NO2. Both indoor and outdoor levels of NO2 were low in this mixed U.S. rural-urban community. Observed spatial and seasonal variation of this air pollutant was largely attributable to urban-rural status, individual level SES, and presence of gas stove in participant’s home. Specifically, children from lower socioeconomic status living in urban areas were exposed to higher outdoor NO2, especially in winter. Indoor NO2 was higher in rural areas as well as in families with a gas stove regardless of seasons.
For seasonal and geographic variability, the median indoor level of NO2 in rural areas was higher than in urban areas in both seasons (6.7vs 3.2 for winter; P = .022, 5.7 vs 2.3 for summer; P = .030), but not outdoor. While the proportion of families who had gas stoves in their homes in this study population was lower than the U.S. average (29%vs 38%), we did find higher indoor NO2 levels in homes with gas stoves in both seasons which is consistent with the literature10,36 and may explain higher indoor levels of NO2 in rural areas of our study setting compared to urban areas as presence of a gas stove was statistically higher in families living in rural areas (vs urban areas; 78%vs 21%; P = .001). Albeit statistically not significant, the correlation between indoor and outdoor NO2 concentrations was higher in the summer (.14) compared to winter (−.03), which may be influenced by the frequency and duration of opening windows during summer. Other unmeasured factors such as emissions from indoor sources and NO2 reactivity might also contribute to differential exposure pattern. 15
Apart from seasonal and geographic factors, residential characteristics played an important role in determining NO2 concentrations. While literature shows inconsistent results on the association between outdoor NO2 and SES measures,37,38 a negative relationship between outdoor NO2 and individual-level SES as defined by the HOUSES index was observed in both seasons in our study. This relationship was not found for indoor NO2. Literature suggests that elevated outdoor NO2 is strongly associated with urban areas with high traffic volume39,40 which is consistent with our data, especially in winter. This observation might be pertinent to asthma epidemiology and contribute to disparities in asthma incidence. For example, we previously reported that children living in homes located in high-traffic areas have a higher risk of developing asthma (hazard ratios, 1.3-1.6 depending on the matching methods) compared with the matched counterparts in this community. 41 Outdoor NO2 levels are known to be higher in winter than in summer, partly attributable to increased fossil fuel use for domestic heating and transportation,42,43 and lower mixing layer height in winter. 44 This may contribute to a disproportionate burden of asthma among people from lower SES, especially those living in urban areas and areas with high traffic volumes. Given seasonality of asthma exacerbation, future studies with a larger sample size need to assess the effect of traffic volume and NO2 concentration and other air quality measures in the context of seasonal and geographic variability on asthma outcomes.
Our exploratory analysis suggested that sex might be an effect modifier for the association between NO2 and asthma outcomes; higher indoor NO2 was associated with increased daytime asthma symptoms among females in summer and associated with overall poorer asthma control status among males in winter, suggesting potentially greater impact of indoor NO2 than outdoor NO2 (data not shown). 45 Sex differences in behaviors by season as well as sex differences in lung growth rates and hormone exposure might contribute to this finding. Interestingly, a review paper highlighting the effects of sex on air pollution epidemiology included studies that showed stronger effects among boys in early life and among girls in later childhood. 46
Our study findings have a few implications on asthma care and research. Our study results may help clinicians to contextualize the potential impact of air quality on health outcomes in mixed rural-urban communities. For example, in rural areas, they may target families with gas stoves to mitigate the potential impact of indoor NO2, while in urban areas, they may want to focus on families with lower SES to mitigate the impact of outdoor NO2, especially for patients with asthma or chronic obstructive pulmonary diseases. To reduce exposure to air pollutants for families who use a gas stove and live in rural areas with low outdoor NO2 levels (eg, 1.8-2.7 ppb), more frequent indoor ventilation or switching from using gas stoves to electric stoves for cooking may be recommended.8,47 Acknowledging the impact of these factors would be helpful when providing guidance to patients with respiratory or cardiovascular diseases on actional interventions that would be most effective at decreasing the potential burden of environmental factors. As our study showed a relatively low level of outdoor NO2 and limited association with asthma control status in children, further research may be warranted for examination of other outdoor triggers (eg, fine particle exposures such as Particulate Matter 2.5 (PM2.5)).
A major strength of this study includes the 62 broad collection sites across Olmsted County (Figure 2). This is an important step toward developing spatial models of outdoor NO2 concentrations which can be applied to diverse health outcomes such as asthma, chronic obstructive pulmonary diseases, cardiovascular diseases, mental disorders, mortality, etc., for large-sample sized population-based studies.48-50 Another strength is using the HOUSES index, an objective individual-level socioeconomic status (SES) measure with precision, which demonstrated an inverse association with outdoor NO2 levels (Table 2). This study also has some limitations. First, our cohort was subject to selection bias by including a disproportionately higher number of families from a higher SES background (all parents or guardians with a college degree or above). As a feasibility study, we made the decision to enroll all those who chose to accept the study invitation, but future investigators may find it prudent to use the HOUSES index to guide enrollment in order to procure a well-balanced sample set regarding SES by monitoring the distribution of SES among enrolled subjects. Another limitation is the short collection period (4 weeks total). However, by collecting NO2 during the heating and non-heating seasons, we were able to account for seasonality by reporting both repeated measures analysis and stratified analysis by season. A few additional factors might have impacted our study results. Families who had gas stoves in our study also tended to have an air cleaner (33%) compared to those without gas stoves (21%), albeit not statistically significant. The inconsistent placement or location of the Palmes tube/monitor during the measurement period (ie, within the living room) with respect to where the air cleaner is located (data not collected) may affect this result. Also, only 12 families (about 20%) reported they had used a gas fireplace more than half of the time during the study period. Lastly, given the small sample size for assessing the association between the level of NO2 and asthma outcomes, we could not control for other confounding effects, such as SES, other airborne contaminants (eg, PM2.5), and living in close proximity to areas with high grass cover, which have been shown to affect asthma incidence or outcomes.5,6,51
Conclusion
In conclusion, NO2 levels were low in this mixed rural-urban community with geographical and seasonal variations. A rural setting and the presence of a gas stove were associated with indoor NO2 concentrations while an urban setting and lower SES (outdoor) with outdoor NO2 concentrations. These findings may help clinicians contextualize the potential impact of air quality on health outcomes in mixed rural-urban communities and guide researchers to define a lower threshold of NO2 concentration with health effect.
Footnotes
Credit Author Statement
All authors meet the criteria for authorship based on the following four requirements: (1) substantial contributions to the conception or design of the work; or the acquisition, analysis, or interpretation of data for the work; (2) drafting the work or revising it critically for important intellectual content; (3) final approval of the version to be published; and (4) agreement to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. Specifically, CW had full access to all the data in the study, takes responsibility for the integrity of the data and the accuracy of the data analysis, and had authority over manuscript preparation and the decision to submit the manuscript for publication. Study concept and design: CW, JG, ER, PW, BL, and YJ. Acquisition, analysis, or interpretation of data: CW, JG, JB, KK, ER, KS, JP (Plano), LM, JP(Porcher), SC, AD, KP, BS, BL, and YJ. Drafting of the manuscript: CW. Critical revision of the manuscript for important intellectual content: all authors.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This was supported by Marcia T. Kreyling Career Development Award, Mayo Foundation. This study was made possible using the resources of the HOUSES program of the Precision Population Science Lab of the Mayo Clinic. The content is solely the responsibility of the authors and does not necessarily represent the official views of the Mayo Foundation or the HOUSES program.
