Abstract
Wildfire smoke is associated with short-term respiratory outcomes including asthma exacerbation in children. As investigations into developmental wildfire smoke exposure on children’s longer-term respiratory health are sparse, we investigated associations between developmental WSE and first use of respiratory medications. Prescription claims from IBM MarketScan Commercial Claims and Encounters database were linked with wildfire smoke plume data from NASA satellites based on Metropolitan Statistical Area (MSA). A retrospective cohort of live infants (2010-2016) born into MSAs in six western states, having prescription insurance, and whose birthdate was estimable from claims data was constructed (N=184,703); of these, gestational age was estimated for 113,154 infants. MSA, gestational age, and birthdate were used to estimate average weekly smoke exposure days (smoke-day) for each developmental period: three trimesters, and two sequential 12-week periods post-birth. Medications treating respiratory tract inflammation were classified using active ingredient and mode of administration into three categories: : ‘upper respiratory’, ‘lower respiratory’, ‘systemic anti-inflammatory’. To evaluate associations between wildfire smoke exposure and medication usage, Cox models associating smoke-days with first observed prescription of each medication category were adjusted for infant sex, birth-season, and birthyear with a random intercept for MSA. Smoke exposure during postnatal periods was associated with earlier first use of upper respiratory medications (1-12 weeks: hazard ratio (HR)=1.094 per 1-day increase in average weekly smoke-day, 95%CI: (1.005,1.191); 13-24 weeks: HR=1.108, 95%CI: (1.016,1.209)); sex-specific HRs varied by post-birth exposure window. Protective associations were observed during gestational windows for both lower respiratory and systemic anti-inflammatory medications; it is possible that these associations may be a consequence of live-birth bias. These findings suggest wildfire smoke exposure during early postnatal developmental periods impact subsequent early life respiratory health.
Introduction
Over the past two decades, frequency and severity of wildfire events has increased, resulting in greater land area being burned each decade (3.3 and 6.8 million average acres per year in 1990s and 2010s, respectively(1)) and substantial economic impact ($347 billion in 2016) that is expected to grow as wildfire frequency increases (2). Wildfire smoke is a complex chemical mixture of both gases and small particles and a major contributor to particulate matter (PM), including the fine fraction of PM (PM2.5)(3) which has been implicated in over 8.8 million premature deaths worldwide (4).
Particulate matter (PM) exposure during wildfires is associated with acute respiratory and other outcomes in children (5); in addition, wildfire exposures may increase exacerbation-of-asthma events as indicated by hospitalizations, emergency department visits, and outpatient visits (6). In adults, controlled exposure studies of short-term exposures to woodsmoke have underscored the role of smoke exposure in respiratory dysfunction, and provide evidence of both systemic and respiratory inflammation (7). In studies of biomass – a primary fuel source for wildfires – burning and woodsmoke, exposure to these pollutants was associated with asthma related symptoms including wheeze in children (8). Epidemiologic studies (9–14) provide further evidence that wildfire, woodsmoke, and biomass smoke exposures are associated with an increased risk of respiratory infection and reduced lung function in children (15,16) and adults (17–19).
While the acute effects of wildfire or biomass burning exposure have been somewhat established, there is little known about the effects of in utero or in early life wildfire particulate matter exposure on longer term respiratory health outcomes. Lung development continues through gestation into the postnatal period, a likely critical window of respiratory susceptibility to air pollution (20). Additionally, PM exposure during gestation (e.g., 21,22) and early childhood (e.g., 23) is linked to sex-specific respiratory health outcomes in humans (21). Additional evidence comes from recent studies of non-human primates where exposure to wildfire smoke in the post-natal period resulted in sex-specific attenuation of host-defense mediators and impaired lung function in adolescence (24).
As similar respiratory effects may also be possible in humans, we consider the research question: Does exposure to wildfire smoke during gestation or in the early postnatal period result in earlier first use of anti-inflammatory respiratory medications in early childhood? By investigating the association of smoke exposure during multiple developmental periods with first use of respiratory medication, we identify critical windows of development where exposure to wildfire smoke may result in increased respiratory vulnerability in the population of young children.
Methods
Population
We constructed a retrospective cohort of infants from the IBM MarketScan® Commercial Claims and Encounters Research Database (MarketScan) born into six western U.S. states. MarketScan is a proprietary deidentified claims database, comprised of data from private U.S.-based insurance companies. This private claims dataset includes hospitalizations, outpatient visits, services, and prescription claims data.
Children born between January 1, 2010 and December 31, 2016, in metropolitan statistical areas (MSAs) in California, Oregon, Washington, Idaho, Montana, or Nevada (Figure 1; Table S1) were eligible. Utilizing a validated algorithm (25,26) and ICD-9 or ICD-10 codes, birthdates were estimated for each live birth (Appendix 1, Supplementary Material). The MSA was the smallest geographic unit available in MarketScan, and therefore the primary spatial unit for assessing wildfire smoke exposure.
Cohort construction. WA, OR, CA, MT, NV, and ID refer to the states of Washington, Oregon, California, Montana, Nevada, and Idaho, respectively. MSA is Metropolitan Statistical Area and CONUS is Continental United States.
Birthdate and MSA were required to estimate exposure to wildfire smoke for all observations.
Each child was required to have continuous enrollment with the insurer with at least 1 week of prescription coverage starting at birth. For each child, follow-up begins after risk period; specifically, follow-up begins at birth for prenatal exposures, and after the post-natal risk period for postnatal exposures. For each outcome, the end of follow-up, measured in weeks, was either date of the outcome, the end date of continuous enrollment, or December 31, 2018, whichever came first. This resulted in two analytic cohorts: the full cohort with 182,387 liveborn children and a sub-cohort with 113,154 liveborn children that also had estimates for gestational age (hereafter, GA sub-cohort). This study (IRB No:20-2817) was evaluated and approved by the Institutional Review Board at the University of North Carolina at Chapel Hill.
Gestational age estimation
Gestational age was used in the GA sub-cohort to estimate with greater precision the average weekly number of smoke exposure days for the portion of each trimester completed. For example, if a child had a GA of 34 weeks, T1 and T2 estimates would be calculated as described above, but T3 would be the average weekly smoke exposure from the 27th week to the end of the 34th week.
As the MarketScan dataset does not explicitly contain variables for gestational age (GA) at delivery, infant records were linked, where possible, to the birthing parent to allow for GA determination. To estimate GA at birth, we used an algorithm based on ICD-9 codes for the number of completed weeks of gestation, as described in Appendix 2 (Supplementary Material) (Adapted from 27); these codes were obtained either from the infants’ or the birthing parents’ records. Observations with GA were not evenly distributed across time due to the shift from ICD-9 to ICD-10 codes in the fall of 2015. Due to marked changes in GA coding between ICD-9 and ICD-10 codes, we chose to estimate GA based on ICD-9 codes only as they were used for the majority of the study period.
Outcome: First use of prescription respiratory medication
Three outcomes were considered: the first prescription use of (1) upper respiratory, (2) lower respiratory, and (3) systemic anti-inflammatory medication. Classifications for each of the three outcomes are briefly summarized in Table S2. As systemic anti-inflammatory medications are often used in acute cases of inflammation or infection that could involve the lower respiratory tract, the upper respiratory tract or other organ system, we chose to analyze them separately from medications used more specifically for the upper or lower respiratory tract. If a child was prescribed medications from two outcome classes on separate or the same claim date, they were assumed to have the event in analysis of both outcomes. Medications were selected using a combination of mode of administration (MSTFMDS variable in MarketScan), therapeutic drug classes (THERCLS), and therapeutic drug groups (THERGRP) as defined in Micromedex RedBook (28) in consultation with a pediatric allergist, and referencing Up-To-Date (29); final classification of medications were validated by two physicians. The outcome event was defined as date of the first fill of prescription respiratory medication in drug claims data.
Exposure: Average weekly wildfire exposure smoke-days
Wildfire smoke exposure data were obtained from the publicly available NASA satellite imagery-based Hazard Mapping System Fire and Smoke product (30,31) and were previously used to generate smoke-day exposure estimates for each ZIP Code Tabulation Area (ZCTA) in the study area (10). These data use visible plumes observed via satellite imaging as an approximation of true wildfire smoke exposure. Using these data, we constructed MSA-level daily smoke-day exposures as follows: First, we used the official crosswalk files released by the United States Office of Housing and Urban Development to define membership within an MSA for each ZCTA in the study area, for each year of the study (32,33). A MSA was defined as having a smoke-day on a given date if at least 25% of the population of the MSA experienced a smoke-day, based on weighting the binary ZCTA smoke-days by year-specific ZCTA population; we refer to this exposure measurement as the “25% threshold” smoke-day.
In sensitivity analyses, we further examined how estimated proportion of the MSA population exposed in developing the exposure metric might impact observed associations. Because each MSA contained a differing number of ZCTAs, we also defined two additional exposure measures for use in this sensitivity analyses: (1) if any ZCTA within that MSA had a smoke-day on a given date, the MSA was defined as having a smoke-day on that date [hereafter referred to as “0% threshold”]; and (3) if 50% of the population of the MSA experienced a smoke-day on a given date, based on weighting the binary ZCTA smoke-days by year-specific ZCTA population, the MSA was defined as having a smoke-day on that date [50% threshold]. Thus, the “0% threshold” and “50% threshold” roughly correspond to average weekly smoke-days experienced by greater than 0% or 50% of the MSA’s population, respectively (Figure 2). All analyses were conducted using each of these exposure thresholds, with the 25% threshold smoke-day exposure utilized in the primary analyses, and the 0% and 50% thresholds applied in sensitivity analyses.
Spatial distribution of total number of smoke-days (2010-2016) at each MSA for the (a) 0%, (b) 25%, and (c) 50% thresholds. Thresholds are defined as the population-weighted percentage of zip codes within an MSA required to experience a smoke-day, in order to assign a smoke-day to the MSA.
Smoke-day exposures were linked to each liveborn child according to residential MSA at birth. For each child, the number of smoke-days were summed in each exposure period and divided by the number of weeks in the period to produce average weekly smoke-day estimates. This transformation allows comparison across periods with differing durations and across individuals with various GAs. Average weekly smoke-day estimates for each period were generated for each of the three exposure thresholds.
In the full cohort, we estimated the average weekly wildfire smoke-days for the following time periods: each trimester of gestation, and two consecutive 12-week periods after birth. First, second, and third trimesters (T1, T2, & T3) were assumed to be 280-197, 196-99, & 98-1 days before the estimated date of birth, respectively. In the full cohort, postnatal exposure periods, P1 and P2, were defined 1-84 days (1-12 weeks) and 85-168 days (13-24 weeks) after birth, respectively.
In the GA sub-cohort, we additionally used GA to more accurately estimate the average weekly smoke-days for the portion of each trimester completed (clarifying example in Supplementary Material). To distinguish between exposures obtained under the assumption of 40-week gestation (full cohort) and those obtained using estimated GA (GA sub-cohort), we use the suffix ‘-40w’ or ‘-GA’, respectively.
Statistical Analyses
For each trimester and postnatal period, mixed effects Cox proportional hazards models were used to estimate the hazard ratio (HR) for the association between the weekly average smoke-days and each prescription outcome. Follow-up time was counted in weeks. The association between average weekly smoke exposure and first use of medications was estimated using the function below:
where λij(t) is the probability of the first filled respiratory prescription at time t for child j in the ith MSA; the MSA into which a child was born is the random intercept, and wi is vector of coefficients for the random intercept. Smoke is the average weekly number of smoke-days for the given trimester or postnatal period. In addition to the two models of postnatal smoke-day exposures, separate models were estimated for the six gestational smoke-day exposure periods (T1-40w, T2-40w, T3-40w, T1-GA, T2-GA, and T3-GA). A proxy for sex-specific hormones during development, sex was dichotomous (male/female).
Birthyear was included to adjust for time trends in both the propensity to prescribe medications and the increasing number of smoke-days. Socioeconomic status was not available in the MarketScan database within which all children were privately insured. The 4-category covariate Birth Season, defined as birth during Spring (March-May), Summer (June-August), Fall (September-November), or Winter (December-February), was included to adjust for potential confounding.
Results for each trimester exposure are presented both for full cohort (all eligible births, assuming a 40-week gestation) and for the GA sub-cohort (the subset of births with estimable gestational age).
Results for each postnatal exposure are presented for the full cohort only. Hazard ratios (HR (95% confidence interval(CI))) are estimated for a one day increase in weekly average smoke-days. Results are presented for all children as well as stratified on sex.
Though this analysis focuses on population-level effects rather than specific sub-populations or etiologies (work-in-progress), we recognize that preterm birth (PTB) that may drive our main findings. As such, we also stratified by PTB (defined as <37wks) within the GA-sub-cohort.
RESULTS
Across the study region of 60 MSAs in six western states (U.S.A.), each MSA had between 31 (Walla Walla, WA) and 24,860 (Los Angeles-Long Beach-Glendale, California) births recorded in MarketScan during 2010-2016, with mean and median of 3,039 and 914 births per MSA, respectively (Table S1). Los Angeles, Portland, and Seattle contributed the MSAs with the greatest number of births. Of 182,387 eligible births with estimable birthdates (based on claims codes summarized in Appendix 1), 89,066 (48.8%) were female, and 113,154 (62.0%) also had an estimable GA (Table 1). The number of births per year ranged from 18,578 to 32,721, and the sex distribution remained similar across years. In the full cohort, Cetirizine hydrochloride (54.3%) and Mometasone furoate (34.8%) were the most frequently filled upper respiratory prescriptions; Albuterol sulfate (90.0%) was the most frequently filled lower respiratory prescription; Methylprednisolone, Prednisolone, or Prednisone (86.3%) was the most frequently filled systemic anti-inflammatory prescriptions (Table S3).
Cohort composition in the full cohort and the gestational age sub-cohort, overall and by sex.
On average, children were exposed to less than one day of smoke per week in each exposure period; the mean trimester exposure ranged from 0.28 to 0.70 weekly smoke-days, and postnatal mean exposure ranged from 0.37 to 0.74, across the thresholds (Table S4). Descriptive statistics for the full cohort, the GA sub-cohort and those without GA were relatively comparable (Table S5).
Postnatal period
In both the first (P1) and second (P2) postnatal periods, the models reflected an association between smoke-day exposure with earlier first use of upper respiratory medication (P1: HR=1.094 per 1-day increase in weekly smoke-day (1.005, 1.191); P2: HR=1.108 (1.016, 1.209)); Figure 3, Table 2). In the P1 period, female children had the larger effect size and stronger association (HR=1.123 (0.994, 1.27)) as compared to males (HR=1.055 (0.940, 1.185)). In contrast, male children had a larger effect size (HR=1.128 (1.004, 1.266)) and stronger association in P2, as compared to female children (HR=1.072 (0.942, 1.22); Figure 3, Table 2). Associations in primary models (i.e., using the 25% threshold) of lower respiratory as well as systemic anti-inflammatory outcomes were largely null (Figure 3, Table 2).
Mixed effects cox model of smoke exposure during post-partum trimester at the 25% exposure threshold, in all eligible children and stratified by sex. P1 and P2 periods correspond to 0-12 and 13-24 weeks, respectively. Models are adjusted for birth-season, birthyear and, where applicable, sex, with a random intercept for MSA. All exposures are the average weekly number of smoke-days in the given post-partum period. RX = prescribed medication claim.
Hazard ratio and confidence intervals for the relationship of average weekly smoke-days and first use of prescription outcome, in each averaging period at the 25% exposure threshold, for all children (green square), and within strata of male children (purple diamonds) and female children (orange circle). The Cox proportional hazards models were adjusted for birth season, birth year, infant sex and included a random intercept for metropolitan statistical area; infant sex covariate was omitted in infant sex-stratified models. rx = prescribed medication claim; T1-T3 refer to the first-third trimesters; P1 and P2 refer to the first and second 12-week post-natal period.
Gestational period
Positive associations with smoke exposure were observed only in the first trimester for the upper respiratory outcome. The 40w- and GA-exposure models produced HRs in the same direction with comparable magnitude, and with some variation in precision (Figure 3, Table 3).
Mixed effects cox model results of smoke exposure during gestational periods at the 25% exposure threshold, in all eligible children and stratified by sex. Models are adjusted for birth season, birthyear and, where applicable, sex, with a random intercept for MSA. The first and second N values correspond to the entire cohort (full cohort) and the cohort portion whose gestational age estimate was available to estimate wildfire smoke-day exposure (GA sub-cohort), respectively. All exposures are the average weekly number of smoke-days in the given period. RX = prescribed medication claim.
Upper respiratory
In the primary model, the impact of smoke-days on the upper respiratory medications outcome varied by gestational period, with increased risk in the first trimester, and null or decreased risk in the second and third trimester (Figure 3, Table 3). In the first trimester, smoke-day exposure was associated with shorter times to first upper respiratory medication, among all (HR=1.082 (0.985, 1.188)), and female children (HR=1.140 (0.997, 1.303)) among the full cohort (T1-40w), but not in male children (HR=1.026 (0.902, 1.167); Figure 3); these first trimester HRs are comparable in the GA sub-cohort (T1-GA), albeit with less precision in the 95%CI. For all other exposure windows, generally null associations were consistently observed in both the full cohort and the GA sub-cohort.
Lower respiratory
Smoke-days had a null or slightly inverse association across all three trimesters in the 40-week exposure modeled in both the full cohort and the GA sub-cohort (Figure 3, Table 3). In contrast to first trimester upper respiratory outcome findings, first trimester smoke-day exposure among the full cohort (T1-40w) was associated with delayed first lower respiratory medications among all (HR=0.977 (0.951, 1.003)) and male children (HR=0.960 (0.927, 0.993)), but not in female children (HR=0.997 (0.956, 1.039); Figure 3).
Systemic Anti-inflammatory
Second trimester smoke-day exposures in the GA sub-cohort (T2-GA) were inversely associated with longer time to the use of systemic anti-inflammatory medication in all (HR=0.946 (0.905, 0.989)) and male (HR=0.921 (0.869, 0.987)) children (Figure 3, Table 3). Similar, but smaller magnitude associations were observed in the first trimester (T1-GA) (All: HR=0.962 (0.927, 0.997); male: HR=0.959 (0.914, 1.006); Figure 3)). Interestingly, a similar protective association was found among female children in the third trimester (T3-GA) (HR=0.935 (0.892, 0.992)), but not among male children. In the full cohort, the exposures across all trimesters and groups produce consistent and similarly precise associations that are closer to the null than their counterparts in the GA sub-cohort.
Sensitivity analyses
In the postnatal period, HRs and 95%CIs for both upper respiratory and systemic anti-inflammatory outcomes are approximately similar across the 0%, 25% and 50% exposure thresholds (Figure S1). For the lower respiratory outcome, however, results in the exposure sensitivity analyses are more variable across thresholds.
For all trimester outcomes, the various thresholds used to generate smoke-day exposure estimates generally had little effect on HR estimates, with some marginal changes in the confidence interval precision (Figure S1). HR did not markedly differ between term births and PTBs.
DISCUSSION
Using private insurance claims data, we conducted a time-to-event analysis to estimate the association between developmental exposure to wildfire smoke and first use of respiratory medication. We found an association between smoke exposure during the first trimester, and both postnatal periods with shorter time to first use of upper respiratory medication; these results were robust across exposure definitions. This positive association between wildfire exposure and upper respiratory medication used was larger for female children in the first trimester and the first 12-week postnatal period, while it was stronger for male children in the 13-24 postnatal week period. These stratified results may suggest sex-specific windows of vulnerability.
Additionally, for both systemic anti-inflammatory medication and lower respiratory medication, protective effects of wildfire exposure were observed during most gestational windows. Most prominently, we observed a shorter time to first lower respiratory medication use among male children and all children during the first trimester; this finding was relatively consistent across cohort specification (i.e., full cohort vs GA sub-cohort). Similarly, all children and male children had shorter times to first systemic anti-inflammatory medication with increased smoke exposure during the first and second trimesters; this was true also of girls in the first trimester. Whether the mechanism for these observed protective effects is an artifact of observational epidemiology or a causal biological finding cannot be ascertained in this dataset.
Our findings suggest a complex epidemiological relationship between gestational exposure to wildfire smoke and time to first symptoms requiring respiratory prescription use. The juxtaposition of below-null HRs during gestation with evidence in the literature suggesting that ambient air pollution is not protective indicates that this analysis may suffer from selection bias (as do most perinatal epidemiological studies (34)). In administrative datasets, these analyses are necessarily limited to live births, and do not capture pregnancy loss. If the in utero response to wildfire smoke exposure is strong, restricting the cohort to live births could lead to attrition of susceptible pregnancies and to an in utero-outcome association which is biased to appear less harmful or even protective. Indeed, non-human primate work shows pregnancy loss is associated with wildfire smoke exposure (35).
The differing exposure pathways in the gestational and postnatal periods could also account for the discrepancy in findings between these two developmental windows. Broadly, possible pathways for air pollutants to exert influence on a fetus include altered placentation or placenta-mediated processes, and direct fetal exposure. Recent work showed that PM can deposit in placental tissue and, to a lesser extent, pass through the placenta to fetal capillaries (36–38). The increased tidal volume in the birthing parent as well as the increased respiratory rate in young children may result in a larger dose of air pollutants (39,40). Should this increased deposition of wildfire smoke pollutants in the birthing parent impair placentation, there may be increased fetal loss as observed in non-human primates and mentioned above (35). In the surviving fetuses, however, our inconsistent or null findings in the gestational period as compared to positive upper respiratory postnatal findings, could indicate that direct exposure to nasal epithelia, as occurs in the postnatal period, produces a more prominent and detectable increase in use of respiratory medications.
Previous studies of PM exposures suggest specific respiratory conditions that could result from wildfire smoke exposure, and may or may not require prescriptions treatment. For example, developmental exposures to PM or biomass burning are associated with early life deficits in lung function (41–43), and increased susceptibility to acute respiratory infection (ARI) or infant bronchiolitis (44,45). While ARI is often treated with prescription medication in young children, subclinical deficits in lung function associated with developmental air pollution exposure may not always be identified and treated with prescription medication unless they lead to respiratory infection (21,43,46–51). It is possible that study cohort members experienced respiratory symptoms that were transient, sub-clinical or well-managed with non-prescription medication, and thus were not captured in this analysis blunting any potential association.
Because lung development occurs in phases starting in early gestation and continuing into adolescence, damage in early life may have lasting negative effects. As PM is a component of wildfire smoke exposure, the mechanisms by which such negative effects occur likely include a combination of oxidative stress and inflammation, which may further induce physiological (52–54) or epigenetic changes in the birthing parent and/or offspring (55,56). Early gestational exposure may lead to abnormal cellular differentiation or branching due to a disruption of essential redox signaling (52), while later exposure may induce structural or functional abnormalities (53,54). In non-human primates, early life exposure to woodsmoke has been linked to reduced lung function metrics (24) which may be due to epigenetic changes to respiratory tissues (57).
To our knowledge, this is the first analysis describing the impact of exposure to wildfire smoke during development on prescription-treated respiratory conditions in human children. There are several additional strengths in this analysis. We were able to represent multiple exposure scenarios by using various thresholds for assigning smoke-days to an MSA. For example, because the definition of smoke-day is based on the appearance of a plume from satellite data, adjacent areas may indeed be exposed to a lower level of smoke that, while not visible, is still impactful; this scenario is captured by the 0% threshold. We also explored refining exposure estimate using GA. Despite losing 38% of the dataset by including claims-based estimates of GA, exaggerated effect estimates observed in the GA sub-cohort suggest that accounting for GA captures important aspects of exposure relative to fetal development that the assumption of 40-week gestation does not. The use of prescription claims data more accurately captures medication usage than diagnosis or prescription delivery captured in electronic medical records. Further the thorough categorization of these medications by two physicians ensures that they are meaningfully classified. Finally, we were the first to leverage two widely used datasets to answer our research question: well-vetted wildfire exposure data and private claims data able to capture prescription usage.
First use of respiratory prescription medication as described in a private insurance claims database serves as an indicator of respiratory health and vulnerability; this approach is not without limitations. The MarketScan Commercial Claims and Encounters database is limited by omission of those who are not on private insurance and those living in rural areas; this could make our results less generalizable to other groups because of socioeconomic factors. While the Marketscan claims database allows the examination of a large population distributed over a large wildfire-impacted region, the availability of individual-level confounders, such as race, SES, and maternal or second-hand smoking, is unfortunately limited. As observed with other air pollutant exposures (58), SES could modify the effect of smoke exposure on respiratory health. Relatedly, this analysis assumes that all participants are exposed to the same proportion of outdoor air pollution (e.g., that they spend the same amount of time outside and all have the same/no indoor air purification systems), which may contribute to the potential for exposure bias differential on SES or other related factors.
Imprecision in the estimates of residence may lead to exposure misclassification, though this impact is expected to be small given the size of a MSA in which exposures are defined as uniform due to the spatial coarseness of the data. MSA at birth provides a relatively coarse estimate of exposure which could be improved upon by using a dataset which includes more granular residential information. We also assumed that the MSA at time of birth was the same during the entirety of gestation; while likely not universally true, the impact of this assumption is expected to be small with potential bias toward the null (59,60). Additionally, the coarseness of gestational age available in the ICD-9 coding system, wherein codes are available at 2-week intervals between 24-37 weeks (inclusive), somewhat limits the accuracy of our exposure estimates in the GA sub-cohort; such exposure misclassification may also induce bias in the effect estimates. Finally, some respiratory outcomes may be lost by claims data, as it is only able to provide information on conditions treated by prescription for those who have access to necessary care.
The modest differences in associations by sex found in this analysis contribute to growing evidence supporting sex-differential effects of wildfire smoke exposure. There are well-documented sex differences in lung development and asthma (61), including stress-induced oxidative response among fetuses (62). Noted sex-specific reductions in lung function in non-human primates following wildfire smoke exposure during infancy (24) may be exacerbated by prenatal PM exposure (21).
Overall, these findings suggest that first trimester and postnatal wildfire smoke exposure is associated with shorter time to the first upper respiratory medication usage, while gestational wildfire smoke exposure across trimesters is associated with longer time to systemic anti-inflammatory respiratory prescription usage. This research supports the growing body of literature indicating that wildfire smoke exposure in early life poses a health risk for pediatric populations, and adds granularity to the current understanding of sex-specific and trimester-specific effects of in utero wildfire smoke exposure.
Clinicians may consider discussing the potential benefits of reducing exposure using evidenced-based measures such as observing air quality index alerts, placing air purifiers in their home, and wearing face masks to reduce particulate exposures. Future work should investigate specific clinical outcomes, such as acute respiratory infection, and utilize finer spatial and temporal granularity to refine estimates of wildfire smoke exposure. This work also highlights the continued need to evaluate measures to protect against wildfire smoke during key developmental stages.
Data Availability
All wildfire smoke satellite data are publicly available through NASA at https://www.ospo.noaa.gov/Products/land/hms.html Office of Housing and Urban Development data is available at https://www.huduser.gov/portal/datasets/usps_crosswalk.html#data Claims data are proprietary and may be obtained form IBM MarketScan (https://www.ibm.com/products/marketscan-research-databases).
https://www.ospo.noaa.gov/Products/land/hms.html
https://www.huduser.gov/portal/datasets/usps_crosswalk.html#data
Declarations
Authors contributions
RD, MER, and IJ conceived of the study question. AGR, CWC, and IJ obtained study funding. RD, CK, YX, and BSS built the dataset and carried out statistical analyses. KR, AGR, and CW-C advised on statistical analyses and cohort development from claims data. RD classified medications as to the outcome and AC and MH reviewed these classifications. RD, CK, and HVJ drafted the manuscript.
Consent for publication
All authors revised the manuscript critically for intellectual content, gave final approval of the version to be published, and agreed to be accountable for all aspects of the work.
Competing interests
The authors declare they have no competing interests to disclose.
Funding
Supported by the Environmental Protection Agency (EPA) Cooperative Agreement (U.S. EPA CR83578501); HVJ was supported in part by an appointment to the ORISE participant research program supported by an interagency agreement between EPA and DOE. CK was supported in part by a training grant from the National Institute of Environmental Health Sciences [T32ES007018]. BSS was supported by the Genetic Epidemiology of Heart, Lung, and Blood Traits Training Grant [T32HL129982].
Disclaimer
The views expressed in this article are those of the authors and do not necessarily represent the views or policies of the U.S. Environmental Protection Agency.
Availability of data and materials
IBM MarketScan® Research Databases, a proprietary claims dataset, on which these analyses are based are available through IBM (https://www.ibm.com/products/marketscan-research-databases). The environmental data from the NOAA Hazards Mapping System that support the conclusions of this article are available in the U.S. EPA Environmental Dataset Gateway (https://www.ospo.noaa.gov/Products/land/hms.html - maps).
Ethics approval
This study (IRB No:20-2817) was evaluated and approved by the Institutional Review Board at the University of North Carolina at Chapel Hill.
Acknowledgements
The authors would like to acknowledge Michele Jonnson-Funk, Virginia Pate, and Tom Luben for their support and guidance in this developing this work.
Footnotes
Revisions include clarifications and minor corrections that reflect feedback from reviewers in previous journal submissions. This preprint may also be deposited at Environmental Health Journal, at the request of journal editors.
Abbreviations
- WSE
- Wildfire smoke exposure
- MSA
- Metropolitan Statistical Area
- PM
- Particulate matter
- PM2.5
- fine fraction of particulate matter
- MarketScan
- IBM MarketScan® Commercial Claims and Encounters Research Database
- ICD
- International Classification of Disease
- ZCTA
- ZIP Code Tabulation Area
- MSTFMDS
- mode of administration
- THERCLS
- therapeutic drug classes
- THERGRP
- therapeutic drug groups
- GA
- Gestational Age
- T1, T2, T3
- 1st, 2nd and 3rd trimesters, respectively
- P1 and P2
- 1st and 2nd consecutive 12-week postnatal period, respectively
- HR
- hazard ratio
- 95% CI
- 95% confidence interval
- PTB
- pre-term birth