- 1Children’s Health and Environment Program, Child Health Research Centre, The University of Queensland, Brisbane, QLD, Australia
- 2Superfund Research Centre, Louisiana State University, Baton Rouge, LA, United States
- 3National Centre for Epidemiology and Population Health, Australian National University, Canberra, ACT, Australia
- 4School of Medicine and Dentistry, Griffith University, Gold Coast, QLD, Australia
Introduction: The association between air pollution and adverse health outcomes has been extensively documented, with oxidative stress widely considered a contributing factor. However, the precise underlying mechanism(s) remains unclear. Recent studies suggest that environmentally persistent free radicals (EPFRs) may provide the missing connection between air pollution and its detrimental health effects. Nevertheless, the indoor environment has received limited attention in EPFR research. Therefore, in this study, we measured EPFRs in house dust samples from two locations in Australia and examined household characteristics associated with their presence.
Methods: Household characteristics and behaviours that influence indoor air quality were collected from an online questionnaire; 24-h indoor and outdoor air quality were measured with a TSI DustTrak™ DRX Aerosol monitor 8,533; annual indoor and outdoor air quality were matched to two validated, satellite-based land-use regression models; and dust samples were collected from house vacuums. Dust samples were analyzed using nano electron paramagnetic resonance (EPR) to determine spin concentrations, g-factor, and delta H peak-to-peak (Hp-p). Key variables were identified using Lasso-penalized regression models, followed by unpenalized linear regression and post-selection inference to estimate coefficients and assess the robustness of the findings.
Results: Our analysis revealed that factors such as extractor fan usage during cooking, exposure to traffic-related air pollution and ambient PM2.5 levels, indoor combustion activities, seasonal variation, housing construction type, ventilation, and cleaning practices were significantly associated with EPFR concentrations in Australian homes. Notably, consistent use of extractor fans during cooking was strongly and consistently associated with lower EPFR concentrations in house dust across both study locations.
Discussion: Our research provided insight into the potential impact of household characteristics on EPFR concentrations, which can potentially lead to adverse health effects. Future research should link our research findings on factors affecting indoor EPFRs to their potential health effects.
1 Introduction
The association between air pollution and adverse health outcomes has been long established, but the causal mechanism(s) behind this association are still uncertain (1–4). Recent studies propose that an oxidant component of particle air pollution, environmentally persistent free radicals (EPFRs), may serve as an inducer of oxidative stress (OS), potentially bridging the gap between air pollution and its detrimental health impacts (5–14). EPFRs are formed during combustion processes and are typically associated with fine particulate matter (PM2.5; particles with a mass median aerodynamic diameter of ≤2.5 μm). These radicals are commonly produced from sources such as traffic-related air pollution (TRAP), residential combustion activities, cooking, industrial burning, cigarette smoking/vaping, and wildfires (6, 10–12, 14–17). EPFRs are particularly concerning as they are free radicals that can persist in the environment and biological systems for prolonged periods (8, 10–12, 14). Typically, free radicals have short lifetimes—only a few picoseconds—due to the extreme instability of unpaired electron/s in their outermost orbit (8, 11); however, EPFRs are unusually stable and long-lasting. This stability arises from the continuous cycling of their valence electrons between paired and unpaired states, facilitated by electron transfer from organic particles to the surfaces of redox-active transition metals present in particulate matter (PM) (6, 7). As a result, these stabilized EFPRs can persist for extended periods, lasting up to 5 years in certain environments, and are notably resistant to decomposition (8, 12, 14).
Little is known about the presence of EPFRs in the household environment. Studies that have explored household environments mostly focused on airborne EPFRs produced from solid fuel combustion used for cooking and heating, demonstrating household practices can be important sources for residential EPFRs (13, 14, 18–20). In addition, motor vehicle emissions, a major source of ambient EPFRs, have been extensively reviewed (16, 21, 22). Household practices such as air-conditioning and window opening can influence TRAP penetration into indoor environments (23). Consequently, it is believed that EPFRs generated from TRAP can accumulate in homes as house dust through the infiltration of ambient air. Emerging evidence shows the potential for adverse health effects of household EPFRs (10), particularly for children, due to their increased sensitivity to exposures and continuing respiratory development (24). Sly et al. reported that 89 out of 90 house dust samples in Brisbane, Australia, had detectable EPFRs; high concentrations (≥ 6 × 1017 spins/g) were associated with an increased risk of wheezing in young children compared to low EPFR concentrations (< 4 × 1017 spins/g) (10). The health concerns from EPFR exposure highlight the importance of understanding the factors that influence the presence of EPFRs in indoor environments, particularly in residential homes. We hypothesize that household characteristics and behaviors influence EPFR concentration in house dust. This study aimed to determine household characteristics associated with EPFRs in Australian homes.
2 Materials and methods
2.1 Population
The data used in the present study were collected in two cohort studies in Australia: the Early Life Lung Function (ELLF) study and the Barwon Infant Study (BIS).
The ELLF study is a longitudinal birth cohort from Brisbane, Australia, enrolled during pregnancy and followed until age 7. Brisbane is located in Southeast Queensland, and this city has a subtropical climate with average temperatures ranging from 21°C to 30.4°C during summer and 10.5°C to 23.4°C in winter (25). The current study represents an extension of the ELLF project. Participants who had previously submitted at least one house dust sample and lived approximately 100 km from the research center were invited to participate in home visits following their year 7 follow-up. In total, 50 participants were invited, 25 responded, and 24 agreed to participate. Home visits were conducted between February 2021 and July 2022, with each home visited twice—once in summer and once in winter. A total of 23 participants completed both seasonal visits, while one participant completed only the winter visit and was subsequently lost to follow-up. This resulted in 47 individual observations: 23 during summer visits and 24 during winter visits.
The Barwon Infant Study (BIS) is a pre-birth cohort comprising 1,074 infants, recruited using an unselected sampling frame to investigate the early-life origins of immune dysregulation in the modern environment (26). The study is based in the Barwon region of Victoria, Australia, which encompasses urban, rural, and coastal areas. Compared to Brisbane, the Barwon region has a cooler climate, with temperatures ranging from 5.2°C and 14.8°C in winter and 11.9°C to 25.0°C in summer (25). Participant characteristics broadly reflect those of the general Australian population, although there is a lower representation of families from non-English-speaking backgrounds. Further details are available in the cohort profile (26). Participants were followed up at multiple time points, including a home visit at 9 months of age for a sub-sample (n = 228), during which a household assessment was conducted and a dust sample was collected.
2.2 Air quality
2.2.1 Direct measurements
During ELLF home visits, indoor and outdoor air quality were measured simultaneously. Particulate matter (PM) was measured with a TSI DustTrak™ DRX Aerosol monitor 8,533, while NO2 and CO were measured using a Dräger X-am 5,000 with a XXs NO2 LC and XXS CO LC sensor. The monitors were set up inside the home in a common area, with a second set placed outside in an area protected from rain. Where possible, sampling occurred over a 24-h period. Data were collected in 5-min intervals and averaged over the monitoring period.
2.2.2 Modeled air pollution exposure
The BIS did not capture direct air quality measurements. These records were matched to two validated, satellite-based land-use regression models (27, 28) built for the Australian continent from spatial predictors, including land use and satellite information. The residential addresses of each child in both BIS and ELLF were geocoded and matched to an annual estimate of PM2.5 and NO2 to reflect long-term air pollution exposure.
2.3 Household characteristics survey
The ELLF cohort had previously completed an annual survey, which captured information on housing characteristics. We used expert opinion and previous research to determine the household characteristics and behaviors likely to influence indoor air quality and, thus, the presence of EPFRs. These relationships are depicted in a directed acyclic graph (DAG) (see Supplementary Figure S1). DAG is important in epidemiological research as it illustrates the assumed causal pathway between an exposure and an outcome and helps identify potential biases (29). A survey was developed for the ELLF home visits that matched questions asked in the BIS survey, with some additional questions identified by the DAG. A few examples of household characteristics collected were house age, type of heating system (e.g., gas, electric, wood), type of cooling system (e.g., air conditioning, fan) cooktop type, and the use of an extractor fan during cooking. Complete household characteristics are reported in Table 1. Participants completed the questionnaire online up to 1 week prior to each home visit. Responses were validated by study staff during each visit. In the BIS, the questionnaire was developed to allow pooling with international cohorts and included validated instruments where possible (26). A home visit review was conducted by a trained researcher, including a researcher-facilitated survey, and the research team validated responses.

Table 1. Descriptive results of EPFR concentration, air pollution, and season from the Early Life Lung Function (ELLF), Brisbane (n = 46) and Barwon Infant Study (BIS), Geelong, Australia (n = 178).
We combined adjacent categories of categorical variables with low counts (≤ 4 in ELLF and ≤ 13 in BIS, corresponding to 10% of the sample size) and excluded binary variables with low counts (e.g., type of home (house vs. apartment): 43/46 in ELLF and 134/136 in BIS). In the ELLF cohort, the variables type of home, type of oven, and if any family member smokes or vapes were not able to be analyzed due to low counts; in the BIS cohort, the variables type of home, major renovations, use of mosquito coils, open-plan kitchen and heating in the child’s room were not analyzed for the same reason.
2.4 Dust collection
The primary outcome of this study was the presence of EPFRs in house dust. Dust samples were collected from household vacuums by the research team during home visits. Samples were frozen at −20°C in the laboratory of the Centre for Children’s Health Research, Brisbane, or Barwon Biomedical Research Laboratory, Geelong, prior to shipping to Louisiana State University for analysis. Dust samples were processed using a sieve shaker on a mesh size (<38 μm), in which radical concentrations were expressed in terms of spins per gram. The dust samples were analyzed using nano electron paramagnetic resonance (EPR) to determine spin concentrations, g-factor, and delta H peak-to-peak (ΔHp-p). (EMX nano, 2016). The correction for the content of oxygen-centered radicals (O-EFPRs) in the samples was based on the analysis of the EPFR spectra parameters, using a linear combination of the 3rd power of the g-tensor shift with respect to the position of the pure oxygen-centered radical (2.0049) and spectral broadening represented by the value of the delta H peak-to-peak. The resulting adjustment factor was then multiplied by the spin-concentration number. EPFRs were reported as the number of spins per gram of dust analyzed (spins/g) and analyzed as continuous variables. Samples with EPFR concentrations below the limit of detection (LOD) were imputed using half the lowest value of EPFR concentration for ELLF (n = 11) and BIS cohort (n = 46), respectively (30). Samples below the LOD for oxygen-centered environmentally persistent free radicals (O-EPFRs) were excluded from the analysis.
2.5 Statistical analysis
We used lasso-penalized regression to select predictors of EPFR and O-EPFR. We repeated 5-fold cross-validation 100 times and chose the lasso tuning parameter lambda (which determines the extent of the penalization) using percentile-lasso to stabilize our estimates against cross-validation variability, using the recommended 75th percentile (31). We used unpenalized linear regression to obtain coefficient estimates and naive confidence intervals for the chosen variables and exact post-selection inference to obtain confidence intervals that account for the uncertainty in the variable selection (32, 33). For the ELLF cohort, in which up to two observations were available per household, we used generalized estimating equations with exchangeable correlation structures to account for the correlation between repeated measures. Although exact post-selection inference was the recommended approach for lasso selective inference in a recent review and simulation study, it does not account for the correlation across repeated measures and can produce confidence intervals that are too wide in high-dimensional data (low signal-to-noise) settings (33, 34).
In contrast, confidence intervals from linear regression on selected variables may be too narrow as they do not account for the uncertainty in the variable selection. We thus report two sets of confidence intervals: from linear regression on the selected variables and from post-selection inference. We used the R package glmnet for lasso regression (35) and the selectiveInference package for exact post-selection inference (32).
The outcome variable, EPFR concentration, was log-transformed, while pollutant variables—including 24-h average indoor and outdoor PM2.5, NO2, and CO, as well as annual ambient PM2.5 and NO2—were log2-transformed. Indoor pollutant concentrations below the detection limit were imputed as half the detection limit.
Indoor-measured PM2.5, CO, and NO2 were included in the ELLF main analysis, whereas modeled ambient PM2.5 and NO2 were included for BIS. Sensitivity analysis was conducted on measured and modeled outdoor PM2.5, CO, and NO2 for ELLF. Random forest regression was used to examine whether our models identified a similar set of variables. Five-fold cross-validation was performed to ensure 100% of the data were used to build the random forest model. Random forest models with optimal hyperparameters were analyzed for variable importance (see Supplementary materials).
3 Results
The EPFRs, air pollution, and household characteristics are shown in Tables 1, 2. The sample size for complete cases in the ELLF cohort was 46 homes, and in the BIS cohort, it was 136 homes. The median (interquartile range, IQR) EPFRs detected in the ELLF and BIS households were 1.65 × 1017 spins/g (IQR = 3.60 × 1017 spins/g) and 5.78 × 1016 spins/g (IQR = 1.21 × 1017 spins/g), respectively (Table 1). Based on the violin plots in Figure 1, the distribution of EPFR concentration in the ELLF cohort is much more spread out than in the BIS cohort, indicating greater variability in the Brisbane region. The EPFR concentration measured in both cohorts was below the low concentration group (< 4 × 1017 spins/g), as reported by Sly et al. (10), which found a protective association against wheezing in young children when compared to the high concentration group (≥ 6 × 1017 spins/g). Seasonal differences regarding EPFR concentration were inconsistent between the two geographical cohorts. In the ELLF cohort, lower median EPFRs were seen in warmer months of summer (7.06 × 1014 spins/g, IQR = 1.83 × 1017 spins/g) in comparison to colder months of winter (3.10 × 1017 spins/g, IQR = 3.13 × 1017 spins/g) (Figure 1). In the BIS cohort, there were slightly greater median EPFRs in summer (6.31 × 1016 spins/g, IQR = 4.53 × 1016 spins/g) when compared to winter (5.74 × 1016, IQR = 1.20 × 1017 spins/g) (Figure 1). However, the violin plots demonstrated longer upper tails in winter for both cohorts, indicating that the highest concentrations of EPFR tend to occur during colder weather.

Table 2. Descriptive results of household characteristics from the Early Life Lung Function (ELLF), Brisbane (n = 46), and Barwon Infant Study (BIS), Geelong, Australia (n = 178).

Figure 1. Distribution of EPFR concentration by season in the ELLF and BIS cohorts. EPFR concentration is presented on a logarithmic scale of base 10.
Annual air pollution exposure was similar for both cohorts. The average annual ambient PM2.5 concentration was 6.13 μg/m3 (SD = 0.32 μg/m3) in the ELLF cohort and 7.25 μg/m3 (SD = 0.69 μg/m3) in the BIS cohort. For annual average NO2, concentrations were 12.25 μg/m3 (SD = 3.97 μg/m3) in ELLF and 5.36 μg/m3 (SD = 2.08 μg/m3) in the BIS cohort. In the ELLF cohort, the monitored 24-h average PM2.5 concentrations were 13.55 μg/m3 indoors and 11.63 μg/m3 outdoors.
3.1 ELLF cohort
Seven household characteristics were associated with EPFR concentrations in house dust (Table 3 and Figure 2). The winter season was associated with the largest increase in EPFR concentration, followed by factors that affect dust re-suspension, such as using dry dusting methods to clean surfaces or major renovations to the home (Figure 2). Factors associated with lower household EPFR concentrations included always using an extractor fan during cooking, the presence of fabric flooring in the living area, more frequent floor cleaning, and more frequent opening of windows or doors (Figure 2). However, except for the winter season variable, post-selection inference revealed considerable uncertainty regarding the direction of these associations.

Figure 2. Household characteristics associated with EPFR concentrations in the ELLF cohort presented with a difference in mean log-EPFR concentration, β (log(spins/g)).
The following household characteristics were not selected by the percentile lasso model: 24-h indoor concentrations of PM2.5, CO, and NO2; type of heating and cooling in the living area and child’s room; use of candles/incense; cooktop type; housing material; garage type; house age; surface cleaning method and frequency; presence of pets; use of mosquito coils; occasional use of an extractor fan during cooking; household size; presence of a fireplace; and neighborhood traffic (see Supplementary Table S1).
3.2 BIS cohort
A total of 15 household characteristics were associated with household EPFR concentrations (Table 4 and Figure 3). Factors associated with increased estimates of EPFR concentrations were annual ambient PM2.5, high and some neighborhood traffic, factors related to indoor combustion practices such as use of candles or incense and having a fireplace, age of house, weatherboard housing material, higher frequency of cleaning in the child’s room, and presence of carpet or rug in the house (Figure 3). Meanwhile, lower EPFR concentrations were associated with the spring season, always using the extractor fan during cooking, homes with enclosed garages, and gas ovens (Figure 3). However, with the exception of two variables (frequency of cleaning in the child’s room and type of oven), there was substantial uncertainty in the direction of the associations when we accounted for the uncertainty in the variable selection.

Figure 3. Household characteristics associated with EPFR concentrations in the BIS cohort presented with the difference in mean log-EPFR concentration, β (log(spins/g)).
Household variables not chosen by the percentile lasso included annual ambient NO2, frequency of opening windows or doors, method and frequency of cleaning in the living room, type of cooktop, occasionally using the extractor fan during cooking, the heating method in the living room, cooling method in the living room and child’s bedroom, smoking or vaping, total number of cigarettes per day, and other seasons (autumn and winter) (see Supplementary Table S2).
3.3 BIS and ELLF cohort
While the household characteristics collection differed slightly between the ELLF and BIS cohorts, two variables were found to be consistently associated with EPFR concentration in both cohorts: the use of an extractor fan and the presence of carpet/rug in the living room (Figure 4). Using an extractor fan during cooking consistently reduced EPFR concentration in both cohorts compared to never using one. However, the presence of a rug/carpet demonstrated an inconsistent direction of association of both cohorts, with the BIS cohort displaying an increase in EPFR concentration, while the ELLF cohort showed a decrease in EPFR concentration (Figure 4).

Figure 4. Household characteristics associated with EPFR in both cohorts presented with the difference in mean log-EPFR concentration, β (log(spins/g)). The 95% confidence intervals are GEE for ELLF and linear regression for BIS.
3.4 Sensitivity analysis
We examined the effect of 24-h monitored outdoor PM2.5 and NO2, as well as CAND-modeled annual outdoor PM2.5 and NO2 in the ELLF cohort. Neither the short-term (24-h average) nor the annual ambient concentrations of particulate or gaseous pollutants were significantly associated with EPFRs or O-EPFRs (see Supplementary Tables S3–S6).
Household characteristics associated with O-EPFR concentration are included in Supplementary materials.
Random forest regression models showed that 17 out of the 28 household characteristics were chosen as important variables in predicting EPFR concentrations in the ELLF cohort, while 18 out of the 28 variables were chosen as important variables in the BIS cohort (see Supplementary Figure S2). In the ELLF cohort, all the variables chosen by lasso (Table 3) were also chosen by random forest. In the BIS cohort, several household characteristics, including some neighborhood traffic, always using an extractor fan during cooking, spring seasonal variation, and fabric flooring in other areas of the house, were chosen by lasso but were not chosen by random forest (see Supplementary Figure S3). However, the majority of the household characteristics chosen by lasso were also chosen by random forest (Table 4). To assess the robustness of model estimates, we conducted a non-parametric bootstrap analysis for the ELLF and BIS cohorts (see Supplementary Tables S7, S8).
4 Discussion
EPFRs are a relatively recent discovery, and research on residential exposure remains limited. To the best of our knowledge, the present study is the first to explore household characteristics associated with EPFR concentrations in house dust. In the ELLF study, we identified seven household characteristics or practices associated with EPFR concentrations, and in the BIS cohort, we identified 15 household characteristics. These household characteristics were indoor combustion activities, exposure to TRAP and ambient PM2.5 levels, seasonal variation, housing construction type, ventilation, and cleaning practices. Although we controlled for different variables in the two study regions, always using an extractor fan during cooking was consistently associated with a lower concentration of EPFRs in the home.
The lower EPFR concentrations observed in homes with consistent extractor fan usage during cooking highlight the significant impact of household cooking practices on indoor air quality. Research on cooking fume by-products has found that toxic compounds such as polycyclic aromatic hydrocarbons (PAH), aldehydes, alkanoic acids, heterocyclic aromatic amines, and other harmful products can be released when cooking at higher temperatures (36–38). It is established in the literature that EPFR concentrations can be generated from solid fuel combustion used for cooking (14, 20), but a close co-existence of EPFRs and cooking fume by-products, particularly PAHs, has also been reported (39–41). Regardless of whether EPFR originates from the fuel source used in cooking or from heating the food at high temperatures, the reduction of cooking emissions through regular use of an extractor appears to lower this oxidant component in house dust. It should be noted that the range of the confidence interval obtained indicates considerable uncertainty regarding the concise magnitude and direction of effects. These wide intervals may be due to limited statistical power. Further investigations with larger populations would help to validate and refine these estimates.
A systematic review found that major indoor PM2.5 sources were related to combustion activities, including cooking, smoking, wood stoves, and the use of candles and incense (42). Since studies have established the presence of EPFRs on PM surfaces, particularly PM2.5 generated from combustion activities (6, 7, 11, 12, 14), the use of candles, incense, gas stoves/ovens or fireplaces in homes could potentially increase EPFR concentrations in house dust. Our results from the BIS cohort support this hypothesis, demonstrating higher EPFR concentrations for homes with a fireplace and increased use of candles or incense. However, our results found that gas ovens were associated with lower EPFR concentrations than electric ovens in the BIS cohort. This association was contrary to prior assumptions, but this may be explained by interactions between EPFRs and other pollutants released during gas combustion rather than being solely attributed to PMs. A systematic review reported that indoor gas cooking and heating are found to be major sources of indoor NO2 (42), but the relationship between EPFRs and NO2 is still not yet established. One study reported that EPFRs were positively correlated with vehicle emissions on highways, particularly elemental carbon (EC) and NO2 (43). Another study supports the finding, demonstrating that PM2.5 exposure to NO2 increased EPFRs by 5–8 times, whereas PM2.5 exposure to NO and ozone did not change EPFR signals (44). Conversely, an experimental study demonstrated that EPFR concentrations were reduced from interacting with NO2 as EPFRs have the ability to donate electrons to form stable products (45). Additional research is needed to further understand the relationship between EPFRs and NO2 and to assess the impact of gas ovens on indoor air quality.
Xu et al. measured EPFR concentrations in PM2.5 released from various sources and found that PM2.5 from vehicle exhaust emissions generated the highest amount of EPFRs (22). The concentration of EPFRs was also observed to differ between agricultural areas and residential areas, with lower EPFR signals reported in agricultural areas due to lower vehicle emissions (16). The results of the present study demonstrated that annual ambient PM2.5 and neighborhoods with moderate to high traffic were associated with higher EPFR concentrations in the BIS cohort. Our findings on the importance of TRAP and ambient PM2.5 in house dust further suggest that outdoor infiltration significantly contributes to indoor air pollution.
We found the season to be important in both settings. In ELLF, a cohort in a sub-tropical region, EPFRs were higher in winter compared with summer. In BIS, a cool temperate region, EPFRs were lower in spring than summer after controlling for a different set of variables. A limitation of this study is that we did not have samples collected in all seasons in ELLF – study visits were only conducted in summer and winter. This finding was inconsistent with the literature, as previous studies have shown seasonal differences. For example, studies conducted in China found that EPFR concentrations were higher in colder months, with a study estimating an equivalent EPFR exposure approximating 23–73 cigarettes per day during the cold periods (12, 13, 20, 21, 46). The authors attribute this to the widespread use of high central heating, especially coal combustion, during winter. In addition, studies have identified that EPFR concentration with an adjacent oxygen atom was more prominent in colder weather, suggesting incomplete combustion (12, 20, 46). Coal combustion is not a heat source in Brisbane, so more work is needed to understand the potential sources in this setting. Rainfall may also influence seasonal variation in EPFR concentrations. The concentration of EPFRs associated with airborne PM2.5 in Beijing, China, was lower in the summer months, a result attributed to the clearing of PM2.5 by precipitation (13). Brisbane typically experiences higher rainfall during summer (47), whereas Geelong experiences its highest rainfall in spring and lowest rainfall during summer (48). These seasonal differences in precipitation across locations may play a role in the seasonal variation reported in the current study.
In the present study, housing construction type, garage type, house age, and housing material were associated with EPFRs in house dust. A systematic review reported that ventilation, house age, and attached garages are important household characteristics that predict indoor air pollutants in Canadian homes (49). These household structures likely influence the overall airtightness of a home (50–52). The air exchange rate (AER) is typically used to understand the leakiness of a house by determining the rate of outdoor air replacing indoor air in a specified space, and AER is one of the main factors that determine indoor air quality (42, 53). Measured AER > 1 (exchange per hour) indicates that outdoor sources are more pronounced, and indoor PM concentrations tend to closely resemble outdoor levels due to indoor air being replaced with outdoor air (54). An Australian report found older houses (over 50 years) in Australia tended to have higher AER (above the 90th percentile of 1.3 h−1), signifying higher outdoor air penetration than newer houses (55). Further, houses built with weatherboard cladding (1.0 h−1) have higher mean AER values than brick houses (0.5 h−1) (55). Higher air penetration in timber-constructed houses may be due to an increased likelihood of structural degradation from termites, sunlight, and water damage (56). Furthermore, a study investigating air exchange rates (AERs) reported a mean AER of 0.44 h−1 with windows closed and 1.57 h−1 with windows open, highlighting the greater influence of outdoor air on the indoor environment when ventilation is increased (57). Our findings in the BIS cohort were consistent with those of (55), demonstrating higher EPFR concentrations of EPFRs in older homes and those constructed with weatherboard compared to brick. In the ELLF cohort, we found that higher natural ventilation and increased frequency of windows or doors opened in the ELLF cohort were associated with lower levels of EPFRs in house dust after controlling for ambient PM2.5 concentrations. High AER can either improve or worsen indoor air quality, depending on the location of the households and the levels of outdoor air pollution. The difference in EPFR concentration associated with high AER observed across different household characteristics in Brisbane and Geelong can be attributed to the lower annual outdoor PM2.5 levels in the ELLF cohort compared to the BIS cohort. Furthermore, the type of garage can influence indoor air quality, with previous research demonstrating that attached garages or underground parking garages led to higher concentrations of pollutants (volatile organic compounds, PM, CO, and NO2) inside the houses when compared to non-attached garages (42, 51, 52, 58). Pollutants generated in enclosed and attached garages from the car can accumulate and penetrate inside the home through walls and ceilings (51, 52). Therefore, we hypothesized that enclosed garages would be associated with greater EPFR concentrations due to greater levels of pollutants. However, we found a negative association in the BIS cohort. The reason for the observed negative association remains unclear, but we believe that the interaction between NO₂ and EPFRs, as previously discussed, may contribute.
Finally, measures to maintain cleanliness in the house were found to be important household characteristics of EPFRs in Australian homes. Measures such as removing shoes before entering homes, using doormats, removing carpets, dusting with damp cloths, and regular cleaning can reduce dust concentration and indoor air pollutants (57). We hypothesize that these cleaning measures can potentially lower EPFRs generated from indoor combustion or infiltration of outdoor air pollution. Our results showed that only the ELLF cohort demonstrated higher EPFRs for households that do not practice regular cleaning and use a dry cloth or dusting wand to clean surfaces. The presence of fabric flooring, carpets, or rugs throughout the home was associated with greater EPFR concentrations in the BIS cohort but lower EPFR concentrations in the living area for the ELLF cohort. The inconsistencies in the direction of associations between settings, with the BIS cohort demonstrating higher EPFR concentrations with increased cleaning frequency and fabric flooring throughout the home, may be attributed to ineffective cleaning practices. In addition, our findings demonstrated that houses that had undergone a major renovation were associated with higher EPFR concentrations in the ELLF cohort. Research on increased small PM (aerodynamic diameter of 1.7 μm) from house renovations (59) and EPFRs has not been extensively studied. The current understanding is that EPFRs are generated through combustion (6, 10–12, 14). Therefore, more research is warranted to understand the importance and mechanism of increased PM from major construction and EPFR concentrations.
While the results demonstrated that indoor EPFRs can be generated through household activities, components of TRAP—particularly PM2.5—also significantly influence EPFR concentrations within homes.
Our research provided insight into the potential impact of household characteristics on EPFR concentrations, which can potentially lead to adverse health effects. Future research should link our research findings on EPFRs—affected by factors such as the use of extractor fans during cooking, indoor combustion activities, exposure to TRAP and ambient PM2.5 levels, seasonal variation, housing construction type, ventilation, and cleaning practices—to their potential health effects.
This study has a few limitations that need to be acknowledged. Although the Barwon Infant Study (BIS) provided a valuable comparison cohort to the ELLF study, we were unable to age-match the children’s health outcomes, as dust samples in BIS were collected only at 9 months of age. As a result, our analysis focused solely on investigating household characteristics associated with EPFR concentrations in house dust across the two locations. Another limitation is the relatively small sample size, which constrained the statistical power of the study. Recruitment in both cohorts was limited by factors such as time, equipment availability, research staff, and funding limitations. Despite these challenges, the study also has several notable strengths. First, we anticipated that the small sample size in the ELLF cohort could lead to instability in model selection and wide confidence intervals. To address these concerns, we included data from the larger BIS cohort in Geelong. We also applied multiple approaches to evaluate the robustness of our results: (1) confidence intervals were estimated using both selective inference following lasso variable selection and standard linear models, and (2) we conducted a bootstrap analysis to assess the stability of results under resampling. These complementary methods allowed for a more comprehensive assessment of uncertainty (34). In addition, this prospective study reduced the risk of recall biases and provided more reliable information on exposures. Finally, the air pollution data for the ELLF cohort were obtained from monitors placed directly in participants’ homes, providing a more accurate representation of personal pollutant exposure compared to relying solely on publicly available environmental data.
5 Conclusion
This study found several household characteristics associated with EPFR levels in Australian homes, including the use of extractor fans during cooking, exposure to TRAP and ambient PM2.5 levels, indoor combustion activities, seasonal variations, housing construction type, ventilation, and cleaning practices. Environmentally Persistent Free Radicals (EPFRs) are oxidant components of air pollution and may play a role in human health impacts of air pollution exposure. This study found EPFRs are present in house dust and that key factors influencing EPFR concentrations included the usage of extractor fans while cooking, traffic-related air pollution, ambient PM2.5, indoor combustion activities, season, housing structures, ventilation, and cleaning habits. Notably, consistent use of extractor fans while cooking was associated with lower EPFR concentrations in household dust. These findings underscore the importance of household practices and environmental factors in determining indoor air quality and the persistence of pollutants within the home.
Data availability statement
The datasets presented in this article are not readily available because data is only available on application to cohorts. Requests to access the datasets should be directed to ZWxsZkB1cS5lZHUuYXU=, YmlzQGJhcndvbmhlYWx0aC5vcmcuYXU=.
Ethics statement
The studies involving humans were approved by the Children’s Health Queensland Human Research Ethics Committee (HREC/13/QRCH/156). The studies were conducted in accordance with local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
WL: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft. PD: Formal analysis, Investigation, Writing – review & editing. GL: Formal analysis, Methodology, Writing – review & editing. NL: Conceptualization, Supervision, Writing – review & editing. ZX: Conceptualization, Supervision, Writing – review & editing. SC: Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Writing – review & editing. SL: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing. PS: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing. DV: Conceptualization, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. Data collection was supported by the National Institute of Environmental Health Sciences of the National Institutes of Health under Award Number P42ES013648.
Acknowledgments
The authors would like to thank Dr. Trish Gilholm for her expertise, time, and advice on the modeling of our random forests.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Author disclaimer
The content is solely the authors’ responsibility and does not necessarily reflect the official views of the National Institutes of Health.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2025.1603114/full#supplementary-material
References
1. Kelly, FJ. Oxidative stress: its role in air pollution and adverse health effects. Occup Environ Med. (2003) 60:612–6. doi: 10.1136/oem.60.8.612
2. Kelly, FJ, and Fussell, JC. Linking ambient particulate matter pollution effects with oxidative biology and immune responses. Ann N Y Acad Sci. (2015) 1340:84–94. doi: 10.1111/nyas.12720
3. Mudway, IS, Kelly, FJ, and Holgate, ST. Oxidative stress in air pollution research. Free Radic Biol Med. (2020) 151:2–6. doi: 10.1016/j.freeradbiomed.2020.04.031
4. Risom, L, Møller, P, and Loft, S. Oxidative stress-induced DNA damage by particulate air pollution. Mutat Res-Fund Mol Mech Mutagen. (2005) 592:119–37. doi: 10.1016/j.mrfmmm.2005.06.012
5. Dugas, TR, Lomnicki, S, Cormier, SA, Dellinger, B, and Reams, M. Addressing emerging risks: scientific and regulatory challenges associated with environmentally persistent free radicals. Int J Environ Res Public Health. (2016) 13:1–8. doi: 10.3390/ijerph13060573
6. Gehling, W, and Dellinger, B. Environmentally persistent free radicals and their lifetimes in PM2.5. Environ Sci Technol. (2013) 47:8172–8. doi: 10.1021/es401767m
7. Guan, X, Truong, L, S, ML, R, LT, and S, AC. Developmental hazard of environmentally persistent free radicals and protective effect of TEMPOL in zebrafish model. Toxics. (2021) 9, 1–2. doi: 10.3390/toxics9010012
8. Pan, B, Li, H, Lang, D, and Xing, B. Environmentally persistent free radicals: occurrence, formation mechanisms and implications. Environ Pollut. (2019) 248:320–31. doi: 10.1016/j.envpol.2019.02.032
9. Ruan, X, Sun, Y, Du, W, Tang, Y, Liu, Q, Zhang, Z, et al. Formation, characteristics, and applications of environmentally persistent free radicals in biochars: a review. Bioresour Technol. (2019) 281:457–68. doi: 10.1016/j.biortech.2019.02.105
10. Sly, PD, Cormier, SA, Lomnicki, S, Harding, JN, and Grimwood, K. Environmentally persistent free radicals: linking air pollution and poor respiratory health? Am J Respir Crit Care Med. (2019) 200:1062–3. doi: 10.1164/rccm.201903-0675LE
11. Vejerano, EP, Rao, G, Khachatryan, L, Cormier, SA, and Lomnicki, S. Environmentally persistent free radicals: insights on a new class of pollutants. Environ Sci Technol. (2018) 52:2468–81. doi: 10.1021/acs.est.7b04439
12. Xu, Y, Qin, L, Liu, G, Zheng, M, Li, D, and Yang, L. Assessment of personal exposure to environmentally persistent free radicals in airborne particulate matter. J Hazard Mater. (2021) 409:125014. doi: 10.1016/j.jhazmat.2020.125014
13. Xu, Y, Yang, L, Wang, X, Zheng, M, Li, C, Zhang, A, et al. Risk evaluation of environmentally persistent free radicals in airborne particulate matter and influence of atmospheric factors. Ecotoxicol Environ Saf. (2020) 196:110571. doi: 10.1016/j.ecoenv.2020.110571
14. Zhao, J, Shi, L, Duan, W, Li, H, Yi, P, Tao, W, et al. Emission factors of environmentally persistent free radicals in PM(2.5) from rural residential solid fuels combusted in a traditional stove. Sci Total Environ. (2021) 773:145151. doi: 10.1016/j.scitotenv.2021.145151
15. Hasan, F, Khachatryan, L, and Lomnicki, S. Comparative studies of environmentally persistent free radicals on Total particulate matter collected from electronic and tobacco cigarettes. Environ Sci Technol. (2020) 54:5710–8. doi: 10.1021/acs.est.0c00351
16. Wang, C, Huang, Y, Zhang, Z, and Cai, Z. Levels, spatial distribution, and source identification of airborne environmentally persistent free radicals from tree leaves. Environ Pollut. (2020) 257:113353. doi: 10.1016/j.envpol.2019.113353
17. Chen, Q, Sun, H, Wang, J, Shan, M, Yang, X, Deng, M, et al. Long-life type — the dominant fraction of EPFRs in combustion sources and ambient fine particles in Xi'an. Atmos Environ. (2019) 219:117059. doi: 10.1016/j.atmosenv.2019.117059
18. Jia, H, Li, S, Wu, L, Li, S, Sharma, VK, and Yan, B. Cytotoxic free radicals on air-borne soot particles generated by burning wood or low-maturity coals. Environ Sci Technol. (2020) 54:5608–18. doi: 10.1021/acs.est.9b06395
19. Qian, R, Zhang, S, Peng, C, Zhang, L, Yang, F, Tian, M, et al. Characteristics and potential exposure risks of environmentally persistent free radicals in PM2.5 in the three gorges reservoir area, southwestern China. Chemosphere. (2020) 252:126425. doi: 10.1016/j.chemosphere.2020.126425
20. Zhao, J, Shen, G, Shi, L, Li, H, Lang, D, Zhang, L, et al. Real-world emission characteristics of environmentally persistent free radicals in PM2.5 from residential solid fuel combustion. Environ Sci Technol. (2022) 56:3997–4004. doi: 10.1021/acs.est.1c08449
21. Chen, Q, Sun, H, Mu, Z, Wang, Y, Li, Y, Zhang, L, et al. Characteristics of environmentally persistent free radicals in PM2.5: concentrations, species and sources in Xi'an, northwestern China. Environ Pollut. (2019) 247:18–26. doi: 10.1016/j.envpol.2019.01.015
22. Xu, M, Wu, T, Tang, Y-T, Chen, T, Khachatryan, L, Iyer, PR, et al. Environmentally persistent free radicals in PM2.5: a review. Waste Disposal & Sustainable. Energy. (2019) 1:177–97. doi: 10.1007/s42768-019-00021-z
23. Sbihi, H, Brook, JR, Allen, RW, Curran, JH, Dell, S, Mandhane, P, et al. A new exposure metric for traffic-related air pollution? An analysis of determinants of hopanes in settled indoor house dust. Environ Health. (2013) 12:48. doi: 10.1186/1476-069X-12-48
24. Grant, T, Brigham, EP, and McCormack, MC. Childhood origins of adult lung disease as opportunities for prevention. J Allergy Clin Immunol Pract. (2020) 8:849–58. doi: 10.1016/j.jaip.2020.01.015
25. Bureau of Meteorology. Climate statistics for Australian locations: Australian government (2023). Available online at: http://www.bom.gov.au/climate/data/index.shtml (Accesed September 4, 2023)
26. Vuillermin, P, Saffery, R, Allen, KJ, Carlin, JB, Tang, ML, Ranganathan, S, et al. Cohort profile: the Barwon infant study. Int J Epidemiol. (2015) 44:1148–60. doi: 10.1093/ije/dyv026
27. Knibbs, LD, van Donkelaar, A, Martin, RV, Bechle, MJ, Brauer, M, Cohen, DD, et al. Satellite-based land-use regression for continental-scale long-term ambient PM2.5 exposure assessment in Australia. Environ Sci Technol. (2018) 52:12445–55. doi: 10.1021/acs.est.8b02328
28. Knibbs, LD, Hewson, MG, Bechle, MJ, Marshall, JD, and Barnett, AG. A national satellite-based land-use regression model for air pollution exposure assessment in Australia. Environ Res. (2014) 135:204–11. doi: 10.1016/j.envres.2014.09.011
29. Feeney, T, Hartwig, FP, and Davies, NM. How to use directed acyclic graphs: guide for clinical researchers. BMJ. (2025) 388:e078226. doi: 10.1136/bmj-2023-078226
30. Hornung, RW, and Reed, LD. Estimation of average concentration in the presence of nondetectable values. Appl Occup Environ Hyg. (1990) 5:46–51. doi: 10.1080/1047322X.1990.10389587
31. Roberts, S, and Nowak, G. Stabilizing the lasso against cross-validation variability. Comput Stat Data Anal. (2014) 70:198–211. doi: 10.1016/j.csda.2013.09.008
32. Taylor, J, and Tibshirani, R. Post-selection inference for ℓ(1)-penalized likelihood models. Can J Stat. (2018) 46:41–61. doi: 10.1002/cjs.11313
33. Lee, JD, Sun, DL, Sun, Y, and Taylor, JE. Exact post-selection inference, with application to the lasso. Ann Stat. (2016) 44:907–27. doi: 10.1214/15-AOS1371
34. Kammer, M, Dunkler, D, Michiels, S, and Heinze, G. Evaluating methods for lasso selective inference in biomedical research: a comparative simulation study. BMC Med Res Methodol. (2022) 22:206. doi: 10.1186/s12874-022-01681-y
35. Friedman, JH, Hastie, T, and Tibshirani, R. Regularization paths for generalized linear models via coordinate descent. J Stat Softw. (2010) 33:1–22. doi: 10.18637/jss.v033.i01
36. Mehta, BM. Nutritional and toxicological aspects of the chemical changes of food components and nutrients during heating and cooking. In: PCK Cheung and BM Mehta, editors. Handbook of food chemistry. Berlin, Heidelberg: Springer Berlin Heidelberg (2015). 897–936.
37. Svedahl, S, Svendsen, K, Qvenild, T, Sjaastad, AK, and Hilt, B. Short term exposure to cooking fumes and pulmonary function. J Occup Med Toxicol. (2009) 4:9. doi: 10.1186/1745-6673-4-9
38. Chiang, T-A, Wu, P-F, and Ko, Y-C. Identification of carcinogens in cooking oil fumes. Environ Res. (1999) 81:18–22. doi: 10.1006/enrs.1998.3876
39. Sun, Y, Chen, J, Qin, W, Yu, Q, Xin, K, Ai, J, et al. Gas-PM2.5 partitioning, health risks, and sources of atmospheric PAHs in a northern China city: impact of domestic heating. Environ Pollut. (2022) 313:120156. doi: 10.1016/j.envpol.2022.120156
40. Jia, H, Zhao, S, Shi, Y, Zhu, L, Wang, C, and Sharma, VK. Transformation of polycyclic aromatic hydrocarbons and formation of environmentally persistent free radicals on modified montmorillonite: the role of surface metal ions and polycyclic aromatic hydrocarbon molecular properties. Environ Sci Technol. (2018) 52:5725–33. doi: 10.1021/acs.est.8b00425
41. Lammel, G, Kitanovski, Z, Kukučka, P, Novák, J, Arangio, AM, Codling, GP, et al. Oxygenated and nitrated polycyclic aromatic hydrocarbons in ambient air—levels, phase partitioning, mass size distributions, and inhalation bioaccessibility. Environ Sci Technol. (2020) 54:2615–25. doi: 10.1021/acs.est.9b06820
42. Vardoulakis, S, Giagloglou, E, Steinle, S, Davis, A, Sleeuwenhoek, A, Galea, KS, et al. Indoor exposure to selected air pollutants in the home environment: a systematic review. Int J Environ Res Public Health. (2020) 17:8972. doi: 10.3390/ijerph17238972
43. Qin, Y, Zhang, X, Huang, W, Qin, J, Hu, X, Cao, Y, et al. Measurement report: impact of emission control measures on environmental persistent free radicals and reactive oxygen species – a short-term case study in Beijing. Atmos Chem Phys. (2024) 24:8737–50. doi: 10.5194/acp-24-8737-2024
44. Gehling, WM Jr Environmentally persistent free radicals (EPFRs) in PM₂.₅: their contribution to hydroxyl radical formation and atmospheric transformation: Louisiana State University (2013)
45. Borrowman, CK, Zhou, S, Burrow, TE, and Abbatt, JP. Formation of environmentally persistent free radicals from the heterogeneous reaction of ozone and polycyclic aromatic compounds. Phys Chem Chem Phys. (2016) 18:205–12. doi: 10.1039/c5cp05606c
46. Jia, S-M, Wang, D-Q, Liu, L-Y, Zhang, Z-F, and Ma, W-L. Size-resolved environmentally persistent free radicals in cold region atmosphere: implications for inhalation exposure risk. J Hazard Mater. (2023) 443:130263. doi: 10.1016/j.jhazmat.2022.130263
47. Australian Government Bureau of Meteorology. Climate statistics for Australian locations: summary statistics BRISBANE REGIONAL OFFICE: Australian government Bureau of Meteorology; (2025) Available online at: http://www.bom.gov.au/climate/averages/tables/cw_040214.shtml (Accessed November 26, 2024).
48. Australian Government Bureau of Meteorology. Climate statistics for Australian locations: summary statistics GEELONG SEC: Australian government Bureau of Meteorology; (2025) Available online at: http://www.bom.gov.au/climate/averages/tables/cw_087025.shtml (Accessed November 26, 2024).
49. Héroux, M-E, Clark, N, Ryswyk, KV, Mallick, R, Gilbert, NL, Harrison, I, et al. Predictors of indoor air concentrations in smoking and non-smoking residences. Int J Environ Res Public Health. (2010) 7:3080–99. doi: 10.3390/ijerph7083080
50. Ambrose, M, and Syme, M. House energy efficiency inspections project: CSIRO; (2015). Available online at: https://www.energy.gov.au/sites/default/files/house-energy-efficiency-inspections-project-2015.pdf (Accessed May 22, 2025).
51. Mallach, G, St-Jean, M, MacNeill, M, Aubin, D, Wallace, L, Shin, T, et al. Exhaust ventilation in attached garages improves residential indoor air quality. Indoor Air. (2017) 27:487–99. doi: 10.1111/ina.12321
52. Oh, H-J, and Kim, J. Monitoring air quality and estimation of personal exposure to particulate matter using an indoor model and artificial neural network. Sustain For. (2020) 12:3794. doi: 10.3390/su12093794
53. Guo, H, Morawska, L, He, C, and Gilbert, D. Impact of ventilation scenario on air exchange rates and on indoor particle number concentrations in an air-conditioned classroom. Atmos Environ. (2008) 42:757–68. doi: 10.1016/j.atmosenv.2007.09.070
54. Abt, E, Suh, HH, Allen, G, and Koutrakis, P. Characterization of indoor particle sources: a study conducted in the metropolitan Boston area. Environ Health Perspect. (2000) 108:35–44. doi: 10.1289/ehp.0010835
55. Metropolitan Fire and Emergency Services Board A best practice approach to shelter-in-place for Victoria. East Melbourne, Victoria; (2011)
56. Reardon, C, Davis, T, and Downton, P. Yourhome – lightweight framing. Australian Government; (2020). Available online at: https://www.yourhome.gov.au/materials/lightweight-framing (Accessed May 22, 2025).
57. Roberts, JW, Wallace, LA, Camann, DE, Dickey, P, Gilbert, SG, Lewis, RG, et al. Monitoring and reducing exposure of infants to pollutants in house dust. In: DM Whitacre, editor. Reviews of environmental contamination and toxicology, vol. 201. Boston, MA: Springer US (2009). 1–39.
58. Baeza_Romero, MT, Dudzinska, MR, Amouei Torkmahalleh, M, Barros, N, Coggins, AM, Ruzgar, DG, et al. A review of critical residential buildings parameters and activities when investigating indoor air quality and pollutants. Indoor Air. (2022) 32:e13144. doi: 10.1111/ina.13144
Keywords: environmentally persistent free radicals, EPFRs, indoor air pollution, particulate matter, house dust
Citation: Lee WR, Dangal P, Langan G, Lazarevic N, Xu Z, Cormier SA, Lomnicki S, Sly PD and Vilcins D (2025) Associations between household characteristics and environmentally persistent free radicals in house dust from two Australian locations. Front. Public Health. 13:1603114. doi: 10.3389/fpubh.2025.1603114
Edited by:
Tong Wang, Duke University, United StatesReviewed by:
Xue Wu, University of California, San Francisco, United StatesXiaoyi Zhang, Jacobi Medical Center, United States
Copyright © 2025 Lee, Dangal, Langan, Lazarevic, Xu, Cormier, Lomnicki, Sly and Vilcins. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dwan Vilcins, ZC52aWxjaW5zQHVxLmVkdS5hdQ==