Long-Term Exposure to Low-Level PM2.5 and Mortality: Investigation of Heterogeneity by Harmonizing Analyses in Large Cohort Studies in Canada, United States, and Europe

BACKGROUND

Studies across the globe generally reported increased mortality risks associated with particulate matter with aerodynamic diameter ≤ 2.5 μ m (PM 2.5 ) exposure with large heterogeneity in the magnitude of reported associations and the shape of concentration-response functions (CRFs). We aimed to evaluate the impact of key study design factors (including confounders, applied exposure model, population age, and outcome definition) on PM 2.5 effect estimates by harmonizing analyses on three previously published large studies in Canada [Mortality-Air Pollution Associations in Low Exposure Environments (MAPLE), 1991-2016], the United States (Medicare, 2000-2016), and Europe [Effects of Low-Level Air Pollution: A Study in Europe (ELAPSE), 2000-2016] as much as possible.

METHODS

We harmonized the study populations to individuals 65 + years of age, applied the same satellite-derived PM 2.5 exposure estimates, and selected the same sets of potential confounders and the same outcome. We evaluated whether differences in previously published effect estimates across cohorts were reduced after harmonization among these factors. Additional analyses were conducted to assess the influence of key design features on estimated risks, including adjusted covariates and exposure assessment method. A combined CRF was assessed with meta-analysis based on the extended shape-constrained health impact function (eSCHIF).

RESULTS

More than 81 million participants were included, contributing 692 million person-years of follow-up. Hazard ratios and 95% confidence intervals (CIs) for all-cause mortality associated with a 5 - μ g / m 3 increase in PM 2.5 were 1.039 (1.032, 1.046) in MAPLE, 1.025 (1.021, 1.029) in Medicare, and 1.041 (1.014, 1.069) in ELAPSE. Applying a harmonized analytical approach marginally reduced difference in the observed associations across the three studies. Magnitude of the association was affected by the adjusted covariates, exposure assessment methodology, age of the population, and marginally by outcome definition. Shape of the CRFs differed across cohorts but generally showed associations down to the lowest observed PM 2.5 levels. A common CRF suggested a monotonically increased risk down to the lowest exposure level. https://doi.org/10.1289/EHP12141.

Inequality in the Distribution of Air Pollution Attributable Mortality Within Canadian Cities

Recent studies have identified inequality in the distribution of air pollution attributable health impacts, but to our knowledge this has not been examined in Canadian cities. We evaluated the extent and sources of inequality in air pollution attributable mortality at the census tract (CT) level in seven of Canada's largest cities. We first regressed fine particulate matter (PM2.5) and nitrogen dioxide (NO2) attributable mortality against the neighborhood (CT) level prevalence of age 65 and older, low income, low educational attainment, and identification as an Indigenous (First Nations, Métis, Inuit) or Black person, accounting for spatial autocorrelation. We next examined the distribution of baseline mortality rates, PM2.5 and NO2 concentrations, and attributable mortality by neighborhood (CT) level prevalence of these characteristics, calculating the concentration index, Atkinson index, and Gini coefficient. Finally, we conducted a counterfactual analysis of the impact of reducing baseline mortality rates and air pollution concentrations on inequality in air pollution attributable mortality. Regression results indicated that CTs with a higher prevalence of low income and Indigenous identity had significantly higher air pollution attributable mortality. Concentration index, Atkinson index, and Gini coefficient values revealed different degrees of inequality among the cities. Counterfactual analysis indicated that inequality in air pollution attributable mortality tended to be driven more by baseline mortality inequalities than exposure inequalities. Reducing inequality in air pollution attributable mortality requires reducing disparities in both baseline mortality and air pollution exposure.

Mortality inequalities of Black adults in Canada

Background

Mortality rates in Canada have been shown to vary by population group (e.g., Indigenous peoples, immigrants) and social economic status (e.g., income levels). Mortality patterns for some groups, including Black individuals, are not as well known. The objective of this study was to assess cause-specific mortality for Black adults living in Canada.

Methods

Mortality inequalities between Black and White cohort members were estimated by sex using Cox proportional hazard models, based on data from the 2001, 2006 and 2011 Canadian Census Health and Environment Cohorts (CanCHECs). The CanCHEC cycles were combined and followed for mortality between Census Day and December 31, 2016 or 2019, resulting in a follow-up period of 15.6, 13.6 or 8.6 years, depending on the CanCHEC cycle.

Results

Ischemic heart disease mortality was the leading cause of death among adult Black males (12.9%) and females (9.8%), as it is for adult White males (16.4%) and females (12.4%). Despite reduced risk of all-cause mortality among Black males and females, compared with White cohort members, there was notable increased risk for some cause-specific mortality. For instance, in the age-adjusted model, among the 25 causes of death examined, Black males had an increased risk of dying from four causes (HIV/AIDS, prostate cancer, diabetes mellitus and cerebrovascular disease), compared with White males. Similarly, Black females were at an increased risk for 6 causes of death (HIV/AIDS, stomach cancer, corpus uteri cancer, lymphomas and multiple myeloma, diabetes mellitus, and endocrine disorders) out of the 27 causes of death examined. These relative increased risks persisted for most causes of death after adjustment for differences in important social determinants of health.

Interpretation

Results showed substantial variability in the risk of dying by cause of death between Black and White cohort members. An important step in reducing health inequities is the routine identification and surveillance of different health outcomes by population groups. This study helps fill that information gap.

How low can you go? Air pollution affects mortality at very low levels

The World Health Organization (WHO) recently released new guidelines for outdoor fine particulate air pollution (PM_2.5) recommending an annual average concentration of 5 μg/m³. Yet, our understanding of the concentration-response relationship between outdoor PM_2.5 and mortality in this range of near-background concentrations remains incomplete. To address this uncertainty, we conducted a population-based cohort study of 7.1 million adults in one of the world's lowest exposure environments. Our findings reveal a supralinear concentration-response relationship between outdoor PM_2.5 and mortality at very low (<5 μg/m³) concentrations. Our updated global concentration-response function incorporating this new information suggests an additional 1.5 million deaths globally attributable to outdoor PM_2.5 annually compared to previous estimates. The global health benefits of meeting the new WHO guideline for outdoor PM_2.5 are greater than previously assumed and indicate a need for continued reductions in outdoor air pollution around the world.

Mortality-Air Pollution Associations in Low Exposure Environments (MAPLE): Phase 2

INTRODUCTION

Mortality is associated with long-term exposure to fine particulate matter (particulate matter ≤2.5 μm in aerodynamic diameter; PM2.5), although the magnitude and form of these associations remain poorly understood at lower concentrations. Knowledge gaps include the shape of concentration-response curves and the lowest levels of exposure at which increased risks are evident and the occurrence and extent of associations with specific causes of death. Here, we applied improved estimates of exposure to ambient PM2.5 to national population-based cohorts in Canada, including a stacked cohort of 7.1 million people who responded to census year 1991, 1996, or 2001. The characterization of the shape of the concentration-response relationship for nonaccidental mortality and several specific causes of death at low levels of exposure was the focus of the Mortality-Air Pollution Associations in Low Exposure Environments (MAPLE) Phase 1 report. In the Phase 1 report we reported that associations between outdoor PM2.5 concentrations and nonaccidental mortality were attenuated with the addition of ozone (O3) or a measure of gaseous pollutant oxidant capacity (Ox), which was estimated from O3 and nitrogen dioxide (NO2) concentrations. This was motivated by our interests in understanding both the effects air pollutant mixtures may have on mortality and also the role of O3 as a copollutant that shares common sources and precursor emissions with those of PM2.5. In this Phase 2 report, we further explore the sensitivity of these associations with O3 and Ox, evaluate sensitivity to other factors, such as regional variation, and present ambient PM2.5 concentration-response relationships for specific causes of death.

METHODS

PM2.5 concentrations were estimated at 1 km2 spatial resolution across North America using remote sensing of aerosol optical depth (AOD) combined with chemical transport model (GEOS-Chem) simulations of the AOD:surface PM2.5 mass concentration relationship, land use information, and ground monitoring. These estimates were informed and further refined with collocated measurements of PM2.5 and AOD, including targeted measurements in areas of low PM2.5 concentrations collected at five locations across Canada. Ground measurements of PM2.5 and total suspended particulate matter (TSP) mass concentrations from 1981 to 1999 were used to backcast remote-sensing-based estimates over that same time period, resulting in modeled annual surfaces from 1981 to 2016.

Annual exposures to PM2.5 were then estimated for subjects in several national population-based Canadian cohorts using residential histories derived from annual postal code entries in income tax files. These cohorts included three census-based cohorts: the 1991 Canadian Census Health and Environment Cohort (CanCHEC; 2.5 million respondents), the 1996 CanCHEC (3 million respondents), the 2001 CanCHEC (3 million respondents), and a Stacked CanCHEC where duplicate records of respondents were excluded (Stacked CanCHEC; 7.1 million respondents). The Canadian Community Health Survey (CCHS) mortality cohort (mCCHS), derived from several pooled cycles of the CCHS (540,900 respondents), included additional individual information about health behaviors. Follow-up periods were completed to the end of 2016 for all cohorts. Cox proportional hazard ratios (HRs) were estimated for nonaccidental and other major causes of death using a 10-year moving average exposure and 1-year lag. All models were stratified by age, sex, immigrant status, and where appropriate, census year or survey cycle. Models were further adjusted for income adequacy quintile, visible minority status, Indigenous identity, educational attainment, labor-force status, marital status, occupation, and ecological covariates of community size, airshed, urban form, and four dimensions of the Canadian Marginalization Index (Can-Marg; instability, deprivation, dependency, and ethnic concentration). The mCCHS analyses were also adjusted for individual-level measures of smoking, alcohol consumption, fruit and vegetable consumption, body mass index (BMI), and exercise behavior.

In addition to linear models, the shape of the concentration-response function was investigated using restricted cubic splines (RCS). The number of knots were selected by minimizing the Bayesian Information Criterion (BIC). Two additional models were used to examine the association between nonaccidental mortality and PM2.5. The first is the standard threshold model defined by a transformation of concentration equaling zero if the concentration was less than a specific threshold value and concentration minus the threshold value for concentrations above the threshold. The second additional model was an extension of the Shape Constrained Health Impact Function (SCHIF), the eSCHIF, which converts RCS predictions into functions potentially more suitable for use in health impact assessments. Given the RCS parameter estimates and their covariance matrix, 1,000 realizations of the RCS were simulated at concentrations from the minimum to the maximum concentration, by increments of 0.1 μg/m3. An eSCHIF was then fit to each of these RCS realizations. Thus, 1,000 eSCHIF predictions and uncertainty intervals were determined at each concentration within the total range.

Sensitivity analyses were conducted to examine associations between PM2.5 and mortality when in the presence of, or stratified by tertile of, O3 or Ox. Additionally, associations between PM2.5 and mortality were assessed for sensitivity to lower concentration thresholds, where person-years below a threshold value were assigned the mean exposure within that group. We also examined the sensitivity of the shape of the nonaccidental mortality-PM2.5 association to removal of person-years at or above 12 μg/m3 (the current U.S. National Ambient Air Quality Standard) and 10 μg/m3 (the current Canadian and former [2005] World Health Organization [WHO] guideline, and current WHO Interim Target-4). Finally, differences in the shapes of PM2.5-mortality associations were assessed across broad geographic regions (airsheds) within Canada.

RESULTS

The refined PM2.5 exposure estimates demonstrated improved performance relative to estimates applied previously and in the MAPLE Phase 1 report, with slightly reduced errors, including at lower ranges of concentrations (e.g., for PM2.5 <10 μg/m3).

Positive associations between outdoor PM2.5 concentrations and nonaccidental mortality were consistently observed in all cohorts. In the Stacked CanCHEC analyses (1.3 million deaths), each 10-μg/m3 increase in outdoor PM2.5 concentration corresponded to an HR of 1.084 (95% confidence interval [CI]: 1.073 to 1.096) for nonaccidental mortality. For an interquartile range (IQR) increase in PM2.5 mass concentration of 4.16 μg/m3 and for a mean annual nonaccidental death rate of 92.8 per 10,000 persons (over the 1991-2016 period for cohort participants ages 25-90), this HR corresponds to an additional 31.62 deaths per 100,000 people, which is equivalent to an additional 7,848 deaths per year in Canada, based on the 2016 population. In RCS models, mean HR predictions increased from the minimum concentration of 2.5 μg/m3 to 4.5 μg/m3, flattened from 4.5 μg/m3 to 8.0 μg/m3, then increased for concentrations above 8.0 μg/m3. The threshold model results reflected this pattern with -2 log-likelihood values being equal at 2.5 μg/m3 and 8.0 μg/m3. However, mean threshold model predictions monotonically increased over the concentration range with the lower 95% CI equal to one from 2.5 μg/m3 to 8.0 μg/m3. The RCS model was a superior predictor compared with any of the threshold models, including the linear model.

In the mCCHS cohort analyses inclusion of behavioral covariates did not substantially change the results for both linear and nonlinear models. We examined the sensitivity of the shape of the nonaccidental mortality-PM2.5 association to removal of person-years at or above the current U.S. and Canadian standards of 12 μg/m3 and 10 μg/m3, respectively. In the full cohort and in both restricted cohorts, a steep increase was observed from the minimum concentration of 2.5 μg/m3 to 5 μg/m3. For the full cohort and the <12 μg/m3 cohort the relationship flattened over the 5 to 9 μg/m3 range and then increased above 9 μg/m3. A similar increase was observed for the <10 μg/m3 cohort followed by a clear decline in the magnitude of predictions over the 5 to 9 μg/m3 range and an increase above 9 μg/m3. Together these results suggest that a positive association exists for concentrations >9 μg/m3 with indications of adverse effects on mortality at concentrations as low as 2.5 μg/m3.

Among the other causes of death examined, PM2.5 exposures were consistently associated with an increased hazard of mortality due to ischemic heart disease, respiratory disease, cardiovascular disease, and diabetes across all cohorts. Associations were observed in the Stacked CanCHEC but not in all other cohorts for cerebrovascular disease, pneumonia, and chronic obstructive pulmonary disease (COPD) mortality. No significant associations were observed between mortality and exposure to PM2.5 for heart failure, lung cancer, and kidney failure.

In sensitivity analyses, the addition of O3 and Ox attenuated associations between PM2.5 and mortality. When analyses were stratified by tertiles of copollutants, associations between PM2.5 and mortality were only observed in the highest tertile of O3 or Ox. Across broad regions of Canada, linear HR estimates and the shape of the eSCHIF varied substantially, possibly reflecting underlying differences in air pollutant mixtures not characterized by PM2.5 mass concentrations or the included gaseous pollutants. Sensitivity analyses to assess regional variation in population characteristics and access to healthcare indicated that the observed regional differences in concentration-mortality relationships, specifically the flattening of the concentration-mortality relationship over the 5 to 9 μg/m3 range, was not likely related to variation in the makeup of the cohort or its access to healthcare, lending support to the potential role of spatially varying air pollutant mixtures not sufficiently characterized by PM2.5 mass concentrations.

CONCLUSIONS

In several large, national Canadian cohorts, including a cohort of 7.1 million unique census respondents, associations were observed between exposure to PM2.5 with nonaccidental mortality and several specific causes of death. Associations with nonaccidental mortality were observed using the eSCHIF methodology at concentrations as low as 2.5 μg/m3, and there was no clear evidence in the observed data of a lower threshold, below which PM2.5 was not associated with nonaccidental mortality.

Long-term exposure to wildfires and cancer incidence in Canada: a population-based observational cohort study

BACKGROUND

Wildfires emit many carcinogenic pollutants that contaminate air, water, terrestrial, and indoor environments. However, little is known about the relationship between exposure to wildfires and cancer risk. We aimed to assess the associations between residential exposure to wildfires and the incidence of several cancer outcomes (lung cancer, brain cancer, non-Hodgkin lymphoma, multiple myeloma, and leukaemia) in Canada.

METHODS

We did a population-based observational cohort study of participants in the 1996 Canadian Census Health and Environment Cohort. The 1996 Canadian Census Health and Environment Cohort is a nationally representative sample of Canadian adults, followed up for cancer incidence and mortality from 1996 to 2015. For this analysis, we excluded participants who lived in major Canadian cities (with a population size greater than 1·5 million people), recent immigrants, and individuals younger than 25 years or 90 years of age or older at baseline. Exposures to wildfires were assigned on the basis of area burned within a 20 km or 50 km radius of residential locations and updated for annual residential mobility. Multivariable Cox proportional hazards models were used to estimate associations between exposure to wildfires and specific cancers associated with carcinogenic compounds released by wildfires, including lung and brain cancer, non-Hodgkin lymphoma, multiple myeloma, and leukaemia, adjusted for many personal and neighbourhood-level covariates.

FINDINGS

Our analyses included more than 2 million people followed up for a median of 20 years, for a total of 34 million person-years. Wildfire exposure was associated with slightly increased incidence of lung cancer and brain tumours. For example, cohort members exposed to a wildfire within 50 km of residential locations in the past 10 years had a 4·9% relatively higher incidence (adjusted hazard ratio [HR] 1·049, 95% CI 1·028-1·071) of lung cancer than unexposed populations, and a 10% relatively higher incidence (adjusted HR 1·100, 1·026-1·179) of brain tumours. Similar associations were observed for the 20 km buffer size. Wildfires were not associated with haematological cancers in this study, and concentration-response trends were not readily apparent when area burned was modelled as a continuous variable.

INTERPRETATION

Long-term exposure to wildfires might increase the risk of lung cancer and brain tumours. Further work is needed to develop long-term estimates of wildfire exposures that capture the complex mixture of environmental pollutants released during these events.

FUNDING

Canadian Institute for Health Research and Fonds de recherche du Quebec.

Neighbourhood walkability and mortality: Findings from a 15-year follow-up of a nationally representative cohort of Canadian adults in urban areas

BACKGROUND

Using a nationally representative cohort of Canadian adults, we assessed associations between neighbourhood walkability and cause-specific mortality and investigated whether they differed by socioeconomic status.

METHODS

The study population was drawn from the 2001 Canadian Census Health and Environment Cohort, which contains individual-level data from a random sample of 20% of Canadian households mandated to complete the long-form census. We included those aged ≥ 25 years at baseline who lived in urban and suburban areas. The national death registry was used to ascertain annual vital status. Linkages to annual income tax data provided place of residence. The Canadian Active Living Environments, a national index that summarizes walkability across Canadian neighbourhoods, was assigned to individuals' residential history. The Cox proportional hazards model was used to assess associations between walkability and cause-specific mortality.

RESULTS

A total of 1.8 million participants (52.5% female) accrued 27.3 million person-years and 265 710 deaths during the 15-year follow-up. The adjusted hazard ratio (HR) for living in a highly walkable neighbourhood relative to living in the least walkable neighbourhoods was associated with a 9% (HR: 0.91 [0.88, 0.95]) and 3% (HR: 0.97 [0.94, 0.99]) reduced risk of cardiovascular and all non-accidental mortality, respectively. The strongest benefits of walkability were found among individuals within the lowest education and household income categories, and who lived in the most deprived neighbourhoods. There were no significant associations (most [class 5] versus least [class 1] walkable HR: 0.84 [0.61-1.16]) seen for accidental traffic mortality.

CONCLUSIONS

Canadian adults who live in walkable neighbourhoods have lower rates of cardiovascular and non-accidental mortality, with the greatest benefits seen in those from the lowest socioeconomic groups.

The association between ambient air pollution concentrations and psychological distress

BACKGROUND

A growing number of epidemiological studies have linked air pollution exposure to psychological conditions. Laboratory studies indicate that air pollutants can activate the neuroendocrine stress axis and modulate stress hormone levels, which could contribute to the development or exacerbation of psychological distress. The present study examined the spatial associations between air pollutants (fine particulate matter [PM_2.5], nitrogen dioxide [NO₂] and ground-level ozone [O₃]) and psychological distress among subjects in the most populous provinces in Canada.

DATA AND METHODS

Subjects were sampled from the Canadian Community Health Survey in three regions (Quebec in 2005 [n=25,800], British Columbia and Alberta in 2005 [n=23,000], and Ontario in 2011 [n=36,000]), and were assigned estimates of annual exposure to three ambient air pollutants (PM_2.5, NO₂ and O₃) for the same years. Individual psychological distress was assessed using the Kessler Psychological Distress Scale (K10), based on anxiety and depressive symptoms in the past month. Regression models (both ordinary least squares and simultaneous autoregressive models) were applied to estimate associations between K10 distress scores and each air pollutant, after adjusting for individual (demographic, socioeconomic and behavioural) and neighbourhood covariates..

RESULTS

Psychological distress was positively associated with PM_2.5 and NO₂ in all three regions, and with O₃ in Quebec. However, after further adjusting for individual and neighbourhood covariates, the associations between distress and air pollution remained statistically significant only in Quebec.

DISCUSSION

Some evidence for positive associations between psychological distress and ambient air pollution after adjusting for spatial autocorrelation was found.

Disease assimilation: The mortality impacts of fine particulate matter on immigrants to Canada

BACKGROUND

Immigrants make up 20% of the Canadian population; however, little is known about the mortality impacts of fine particulate matter (PM_2.5) air pollution on immigrants compared with non-immigrants, or about how impacts may change with duration in Canada.

DATA AND METHODS

This study used the 2001 Canadian Census Health and Environment Cohort, a longitudinal cohort of 3.5 million individuals, of which 764,000 were classified as immigrants (foreign-born). Postal codes from annual income tax files were used to account for mobility among respondents and to assign annual PM_2.5 concentrations from 1998 to 2016. Exposures were estimated as a three-year moving average prior to the follow-up year. Cox survival models were used to determine hazard ratios (HRs) for cause-specific mortality, comparing the Canadian and foreign-born populations, with further stratification by year of immigration grouped into 10-year cohorts.

RESULTS

Differences in urban-rural settlement patterns resulted in greater exposure to PM_2.5 for immigrants compared with non-immigrants (mean = 9.3 vs. 7.5 μg/m³), with higher exposures among more recent immigrants. In fully adjusted models, immigrants had higher HRs per 10 μg/m³ increase in PM_2.5 concentration compared with Canadian-born individuals for cardiovascular mortality (HR [95% confidence interval] = 1.22 [1.12 to 1.34] vs. 1.12 [1.07 to 1.18]) and cerebrovascular mortality (HR = 1.25 [1.03 to 1.52] vs. 1.03 [0.93 to 1.15]), respectively. However, tests for differences between the two groups were not significant when Cochran's Q test was used. No significant associations were found for respiratory outcomes, except for lung cancer in non-immigrants (HR = 1.10 [1.02 to 1.18]). When stratified by year of immigration, differences in HRs across varied by cause of death.

DISCUSSION

In Canada, PM_2.5 is an equal-opportunity risk factor, with immigrants experiencing similar if not higher mortality risks compared with non-immigrants for cardiovascular-related causes of death. Some notable differences also existed with cerebrovascular and lung cancer deaths. Continued reductions in air pollution, particularly in urban areas, will improve the health of the Canadian population as a whole.

Residential greenness and indicators of stress and mental well-being in a Canadian national-level survey

BACKGROUND

Residential proximity to greenness in urban areas has been shown to confer a number of health benefits, including improved mental health. We investigated whether greenness was associated with self-reported stress, distress, and mental health among adult participants of multiple cycles of a national Canadian health survey, and whether these associations varied by sex, age, income, and neighbourhood characteristics.

METHODS

Our study population included 397,900 participants of the Canadian Community Health Survey, 18 years of age or older, who lived in census metropolitan areas between 2000 and 2015. We used the Normalized Difference Vegetation Index (NDVI) to characterize participants' exposure to greenness within 250 m, 500 m, and 1 km buffers from a representative location of their postal code. Health outcomes included: self-reported perceptions of life stress, psychological distress, and self-rated mental health. We used multiple regression models, adjusted for relevant individual and neighbourhood-level variables to estimate associations (and 95% confidence intervals) between each outcome and exposure to greenness.

FINDINGS

In models with all participants, we observed 6% lower odds of poor self-rated mental health per increase in the interquartile range (i.e., 0.12) of NDVI within a 500 m buffer. Across the three outcomes, we found substantial heterogeneity in effect size across categories of sex, age, and community-level indicators of deprivation and urban form. For example, each incremental increase in greenness exposure was associated with a reduction of 0.61 (95% CI: 0.81 to -0.51) on the K10 psychological distress score among those living in the active core of cities, and with an increase of 0.07 (95% CI: 0.03-0.12) on this score among those living in the most suburban areas.

CONCLUSIONS

Our results indicate that the potential benefits of residential greenness on mental health vary across personal and neighbourhood-level characteristics and are sensitive to how the outcome is measured. Additional research is needed to understand which features of greenness are most relevant to different sub-groups of the population to maximize these health benefits.

Self-rated stress, distress, mental health, and health as modifiers of the association between long-term exposure to ambient pollutants and mortality

BACKGROUND

Individual and neighbourhood-scale socioeconomic characteristics modify associations between exposure to air pollution and mortality. The role of stress, which may integrate effects of social and environmental exposures on health, is unknown. We examined whether an individual's perspective on their own well-being, as assessed using self-rated measures of stress and health, modifies the pollutant-mortality relationship.

METHODS

The Canadian Community Health Survey (CCHS)-mortality cohort includes respondents from surveys administered between 2001 and 2012 linked to vital statistics and postal codes from 1981 until 2016. Annual fine particulate matter (PM_2.5), nitrogen dioxide (NO₂), and ozone (O₃) exposure estimates were attached to a sample of cohort members aged 30-89 years (n = 398,300 respondents/3,848,400 person-years). We examined whether self-rated stress, distress, mental health, and general health modified associations between long-term exposure to each pollutant (three-year moving average with one-year lag) and non-accidental mortality using Cox survival models, adjusted for individual- (i.e. socioeconomic and behavioural) and neighbourhood-scale covariates.

RESULTS

In fully-adjusted models, the relationship between exposure to pollutants and mortality was stronger among those with poor self-rated mental health, including a significant difference for NO₂ (hazard ratio (HR) = 1.15, 95% CI 1.06-1.25 per IQR) compared to those with very good/excellent mental health (HR = 1.05, 95% CI 1.01-1.08; Cochran's Q = 4.01; p < 0.05). Poor self-rated health was similarly associated with higher pollutant-associated HRs, but only in unadjusted models. Stress and distress did not modify pollutant-mortality associations.

CONCLUSIONS

Poor self-rated mental and general health were associated with increased mortality attributed to exposure to ambient pollutants.

The association between walkable neighbourhoods and physical activity across the lifespan

BACKGROUND

Walkability is positively associated with physical activity in adults. Walkability is more consistently associated with walking for transportation than recreational walking. The purpose of this study is to examine how the association between walkable neighbourhoods and physical activity varies by age and type of physical activity using a new Canadian walkability database.

DATA AND METHODS

The 2016 Canadian Active Living Environments (Can-ALE) database was attached to two cross-sectional health surveys: the Canadian Health Measures Survey (CHMS; 2009 to 2015) and the Canadian Community Health Survey (CCHS; 2015 to 2016). Physical activity was measured in the CHMS using the Actical accelerometer (n = 10,987; ages 3 to 79). Unorganized physical activity outside of school among children aged 3 to 11 was reported by parents in the CHMS (n = 4,030), and physical activity data by type (recreational, transportation-based, school-based, and household and occupational) was self-reported by respondents in the CCHS (n = 105,876; ages 12 and older).

RESULTS

Walkability was positively associated with accelerometer-measured moderate-to-vigorous physical activity in youth (p < 0.05), younger adults (p < 0.0001) and older adults (p < 0.05), while walkability was negatively associated with light physical activity in youth (ages 12 to 17) and older adults (ages 60 to 79) (p < 0.05). Walkability was positively associated with self-reported transportation-based physical activity in youth (p < 0.001) and adults of all ages (p < 0.0001). Walkability was negatively associated with parent-reported unorganized physical activity of children aged 5 to 11, and children living in the most walkable neighbourhoods accumulated 10 minutes of physical activity less-on average-than those living in the least walkable neighbourhoods.

DISCUSSION

The results of this study are consistent with previous studies indicating that walkability is more strongly associated with physical activity in adults than in children and that walkability is associated with transportation-based physical activity. Walkability is one of many built environment factors that may influence physical activity. More research is needed to identify and understand the built environment factors associated with physical activity in children and with recreational or leisure-time physical activity.

An examination of the associations between walkable neighbourhoods and obesity and self-rated health in Canadians

BACKGROUND

Two-thirds of Canadian adults and one-third of Canadian children and youth are overweight or obese. There is increased interest in identifying features of the built environment-such as walkability-that facilitate lifestyle habits associated with reduced obesity and improved health. The purpose of this study is to examine how the associations between walkability and both obesity and self-rated health vary by age in Canadians using a new walkability dataset.

DATA AND METHODS

The 2016 Canadian Active Living Environments (Can-ALE) database was attached to Canadian Health Measures Survey (CHMS; 2009 to 2015) data. Moderate-to-vigorous physical activity (MVPA), light physical activity (LPA) and step counts were measured in the CHMS using the Actical accelerometer (n = 10,852; ages 3 to 79). Body mass index (BMI) and waist circumference were measured in a mobile clinic. Self-rated general and mental health were assessed using a questionnaire.

RESULTS

The percentage of adults aged 40 to 59 classified as overweight or obese was 28 percentage points lower in the most walkable Can-ALE category than in the least walkable category (49.1% vs. 77.5%, p < 0.0125). There was a significant downward linear trend in measured BMI and waist circumference across Can-ALE categories (from least to most walkable) for adults aged 18 to 59, but not for children and youth or older adults aged 60 to 79. MVPA was a significant mediating factor in the association between the Can-ALE index and BMI in adults aged 40 to 79 (and in the waist circumference of respondents aged 40 to 59). Young adults (aged 18 to 39) were more likely than older adults (aged 60 to 79) to report very good or excellent general health as walkability increased.

DISCUSSION

Using a new and freely-available Canadian walkability index, this study observed a positive association between walkability and both measured obesity and self-rated general health in adults. Walkability is one of many built environment characteristics that should be considered when trying to understand the relative contribution of the built environment to a person's weight and overall health.

Cohort profile: The Canadian Census Health and Environment Cohorts (CanCHECs)

The Canadian Census Health and Environment Cohorts (CanCHECs) are population-based linked datasets of the household population at the time of census collection. The CanCHECs combine data from respondents to the long-form census or the National Household Survey between 1991 and 2011 with administrative health data (e.g., mortality, cancer incidence, hospitalizations, emergency ambulatory care) and annual mailing address postal codes. The CanCHEC datasets are rich national data resources that can be used to measure and examine health inequalities across socioeconomic and ethnocultural dimensions for different periods and locations. These datasets can also be used to examine the effects of exposure to environmental factors on human health. Because of their large size, the CanCHECs are an excellent resource for examining rare health outcomes and small population groups. They are ideally suited for environmental health research because of their geographic coverage across all regions of Canada, their long follow-up periods and their linkage to annual postal code history.

Mortality-Air Pollution Associations in Low-Exposure Environments (MAPLE): Phase 1

INTRODUCTION

Fine particulate matter (particulate matter ≤2.5 μm in aerodynamic diameter, or PM2.5) is associated with mortality, but the lower range of relevant concentrations is unknown. Novel satellite-derived estimates of outdoor PM2.5 concentrations were applied to several large population-based cohorts, and the shape of the relationship with nonaccidental mortality was characterized, with emphasis on the low concentrations (<12 μg/m3) observed throughout Canada.

METHODS

Annual satellite-derived estimates of outdoor PM2.5 concentrations were developed at 1-km2 spatial resolution across Canada for 2000-2016 and backcasted to 1981 using remote sensing, chemical transport models, and ground monitoring data. Targeted ground-based measurements were conducted to measure the relationship between columnar aerosol optical depth (AOD) and ground-level PM2.5. Both existing and targeted ground-based measurements were analyzed to develop improved exposure data sets for subsequent epidemiological analyses.

Residential histories derived from annual tax records were used to estimate PM2.5 exposures for subjects whose ages ranged from 25 to 90 years. About 8.5 million were from three Canadian Census Health and Environment Cohort (CanCHEC) analytic files and another 540,900 were Canadian Community Health Survey (CCHS) participants. Mortality was linked through the year 2016. Hazard ratios (HR) were estimated with Cox Proportional Hazard models using a 3-year moving average exposure with a 1-year lag, with the year of follow-up as the time axis. All models were stratified by 5-year age groups, sex, and immigrant status. Covariates were based on directed acyclical graphs (DAG), and included contextual variables (airshed, community size, neighborhood dependence, neighborhood deprivation, ethnic concentration, neighborhood instability, and urban form). A second model was examined including the DAG-based covariates as well as all subject-level risk factors (income, education, marital status, indigenous identity, employment status, occupational class, and visible minority status) available in each cohort. Additional subject-level behavioral covariates (fruit and vegetable consumption, leisure exercise frequency, alcohol consumption, smoking, and body mass index [BMI]) were included in the CCHS analysis.

Sensitivity analyses evaluated adjustment for covariates and gaseous copollutants (nitrogen dioxide [NO2] and ozone [O3]), as well as exposure time windows and spatial scales. Estimates were evaluated across strata of age, sex, and immigrant status. The shape of the PM2.5-mortality association was examined by first fitting restricted cubic splines (RCS) with a large number of knots and then fitting the shape-constrained health impact function (SCHIF) to the RCS predictions and their standard errors (SE). This method provides graphical results indicating the RCS predictions, as a nonparametric means of characterizing the concentration-response relationship in detail and the resulting mean SCHIF and accompanying uncertainty as a parametric summary.

Sensitivity analyses were conducted in the CCHS cohort to evaluate the potential influence of unmeasured covariates on air pollution risk estimates. Specifically, survival models with all available risk factors were fit and compared with models that omitted covariates not available in the CanCHEC cohorts. In addition, the PM2.5 risk estimate in the CanCHEC cohort was indirectly adjusted for multiple individual-level risk factors by estimating the association between PM2.5 and these covariates within the CCHS.

RESULTS

Satellite-derived PM2.5 estimates were low and highly correlated with ground monitors. HR estimates (per 10-μg/m3 increase in PM2.5) were similar for the 1991 (1.041, 95% confidence interval [CI]: 1.016-1.066) and 1996 (1.041, 1.024-1.059) CanCHEC cohorts with a larger estimate observed for the 2001 cohort (1.084, 1.060-1.108). The pooled cohort HR estimate was 1.053 (1.041-1.065). In the CCHS an analogous model indicated a HR of 1.13 (95% CI: 1.06-1.21), which was reduced slightly with the addition of behavioral covariates (1.11, 1.04-1.18). In each of the CanCHEC cohorts, the RCS increased rapidly over lower concentrations, slightly declining between the 25th and 75th percentiles and then increasing beyond the 75th percentile. The steepness of the increase in the RCS over lower concentrations diminished as the cohort start date increased. The SCHIFs displayed a supralinear association in each of the three CanCHEC cohorts and in the CCHS cohort.

In sensitivity analyses conducted with the 2001 CanCHEC, longer moving averages (1, 3, and 8 years) and smaller spatial scales (1 km2 vs. 10 km2) of exposure assignment resulted in larger associations between PM2.5 and mortality. In both the CCHS and CanCHEC analyses, the relationship between nonaccidental mortality and PM2.5 was attenuated when O3 or a weighted measure of oxidant gases was included in models. In the CCHS analysis, but not in CanCHEC, PM2.5 HRs were also attenuated by the inclusion of NO2. Application of the indirect adjustment and comparisons within the CCHS analysis suggests that missing data on behavioral risk factors for mortality had little impact on the magnitude of PM2.5-mortality associations. While immigrants displayed improved overall survival compared with those born in Canada, their sensitivity to PM2.5 was similar to or larger than that for nonimmigrants, with differences between immigrants and nonimmigrants decreasing in the more recent cohorts.

CONCLUSIONS

In several large population-based cohorts exposed to low levels of air pollution, consistent associations were observed between PM2.5 and nonaccidental mortality for concentrations as low as 5 μg/m3. This relationship was supralinear with no apparent threshold or sublinear association.

Low concentrations of fine particle air pollution and mortality in the Canadian Community Health Survey cohort

BACKGROUND

Approximately 2.9 million deaths are attributed to ambient fine particle air pollution around the world each year (PM2.5). In general, cohort studies of mortality and outdoor PM2.5 concentrations have limited information on individuals exposed to low levels of PM2.5 as well as covariates such as smoking behaviours, alcohol consumption, and diet which may confound relationships with mortality. This study provides an updated and extended analysis of the Canadian Community Health Survey-Mortality cohort: a population-based cohort with detailed PM2.5 exposure data and information on a number of important individual-level behavioural risk factors. We also used this rich dataset to provide insight into the shape of the concentration-response curve for mortality at low levels of PM2.5.

METHODS

Respondents to the Canadian Community Health Survey from 2000 to 2012 were linked by postal code history from 1981 to 2016 to high resolution PM2.5 exposure estimates, and mortality incidence to 2016. Cox proportional hazard models were used to estimate the relationship between non-accidental mortality and ambient PM2.5 concentrations (measured as a three-year average with a one-year lag) adjusted for socio-economic, behavioural, and time-varying contextual covariates.

RESULTS

In total, 50,700 deaths from non-accidental causes occurred in the cohort over the follow-up period. Annual average ambient PM2.5 concentrations were low (i.e. 5.9 μg/m3, s.d. 2.0) and each 10 μg/m3 increase in exposure was associated with an increase in non-accidental mortality (HR = 1.11; 95% CI 1.04-1.18). Adjustment for behavioural covariates did not materially change this relationship. We estimated a supra-linear concentration-response curve extending to concentrations below 2 μg/m3 using a shape constrained health impact function. Mortality risks associated with exposure to PM2.5 were increased for males, those under age 65, and non-immigrants. Hazard ratios for PM2.5 and mortality were attenuated when gaseous pollutants were included in models.

CONCLUSIONS

Outdoor PM2.5 concentrations were associated with non-accidental mortality and adjusting for individual-level behavioural covariates did not materially change this relationship. The concentration-response curve was supra-linear with increased mortality risks extending to low outdoor PM2.5 concentrations.

Examining the Shape of the Association between Low Levels of Fine Particulate Matter and Mortality across Three Cycles of the Canadian Census Health and Environment Cohort

BACKGROUND

Ambient fine particulate air pollution with aerodynamic diameter ≤2.5 μm (PM2.5) is an important contributor to the global burden of disease. Information on the shape of the concentration-response relationship at low concentrations is critical for estimating this burden, setting air quality standards, and in benefits assessments.

OBJECTIVES

We examined the concentration-response relationship between PM2.5 and nonaccidental mortality in three Canadian Census Health and Environment Cohorts (CanCHECs) based on the 1991, 1996, and 2001 census cycles linked to mobility and mortality data.

METHODS

Census respondents were linked with death records through 2016, resulting in 8.5 million adults, 150 million years of follow-up, and 1.5 million deaths. Using annual mailing address, we assigned time-varying contextual variables and 3-y moving-average ambient PM2.5 at a 1×1 km spatial resolution from 1988 to 2015. We ran Cox proportional hazards models for PM2.5 adjusted for eight subject-level indicators of socioeconomic status, seven contextual covariates, ozone, nitrogen dioxide, and combined oxidative potential. We used three statistical methods to examine the shape of the concentration-response relationship between PM2.5 and nonaccidental mortality.

RESULTS

The mean 3-y annual average estimate of PM2.5 exposure ranged from 6.7 to 8.0 μg/m3 over the three cohorts. We estimated a hazard ratio (HR) of 1.053 [95% confidence interval (CI): 1.041, 1.065] per 10-μg/m3 change in PM2.5 after pooling the three cohort-specific hazard ratios, with some variation between cohorts (1.041 for the 1991 and 1996 cohorts and 1.084 for the 2001 cohort). We observed a supralinear association in all three cohorts. The lower bound of the 95% CIs exceeded unity for all concentrations in the 1991 cohort, for concentrations above 2 μg/m3 in the 1996 cohort, and above 5 μg/m3 in the 2001 cohort.

DISCUSSION

In a very large population-based cohort with up to 25 y of follow-up, PM2.5 was associated with nonaccidental mortality at concentrations as low as 5 μg/m3. https://doi.org/10.1289/EHP5204.

Evaluation of a method to indirectly adjust for unmeasured covariates in the association between fine particulate matter and mortality

BACKGROUND

Indirect adjustment via partitioned regression is a promising technique to control for unmeasured confounding in large epidemiological studies. The method uses a representative ancillary dataset to estimate the association between variables missing in a primary dataset with the complete set of variables of the ancillary dataset to produce an adjusted risk estimate for the variable in question. The objective of this paper is threefold: 1) evaluate the method for non-linear survival models, 2) formalize an empirical process to evaluate the suitability of the required ancillary matching dataset, and 3) test modifications to the method to incorporate time-varying exposure data, and proportional weighting of datasets.

METHODS

We used the association between fine particle air pollution (PM_2.5) with mortality in the 2001 Canadian Census Health and Environment Cohort (CanCHEC, N = 2.4 million, 10-years follow-up) as our primary dataset, and the 2001 cycle of the Canadian Community Health Survey (CCHS, N = 80,630) as the ancillary matching dataset that contained confounding risk factor information not available in CanCHEC (e.g., smoking). The main evaluation process used a gold-standard approach wherein two variables (education and income) available in both datasets were excluded, indirectly adjusted for, and compared to true models with education and income included to assess the amount of bias correction. An internal validation for objective 1 used only CanCHEC data, whereas an external validation for objective 2 replaced CanCHEC with the CCHS. The two proposed modifications were applied as part of the validation tests, as well as in a final indirect adjustment of four missing risk factor variables (smoking, alcohol use, diet, and exercise) in which adjustment direction and magnitude was compared to models using an equivalent longitudinal cohort with direct adjustment for the same variables.

RESULTS

At baseline (2001) both cohorts had very similar PM_2.5 distributions across population characteristics, although levels for CCHS participants were consistently 1.8-2.0 μg/m³ lower. Applying sample-weighting largely corrected for this discrepancy. The internal validation tests showed minimal downward bias in PM_2.5 mortality hazard ratios of 0.4-0.6% using a static exposure, and 1.7-3% when a time-varying exposure was used. The external validation of the CCHS as the ancillary dataset showed slight upward bias of -0.7 to -1.1% and downward bias of 1.3-2.3% using the static and time-varying approaches respectively.

CONCLUSIONS

The CCHS was found to be fairly well representative of CanCHEC and its use in Canada for indirect adjustment is warranted. Indirect adjustment methods can be used with survival models to correct hazard ratio point estimates and standard errors in models missing key covariates when a representative matching dataset is available. The results of this formal evaluation should encourage other cohorts to assess the suitability of ancillary datasets for the application of the indirect adjustment methodology to address potential residual confounding.

A Comparative Exposure Assessment of Campylobacter in Ontario, Canada

To inform source attribution efforts, a comparative exposure assessment was developed to estimate the relative exposure to Campylobacter, the leading bacterial gastrointestinal disease in Canada, for 13 different transmission routes within Ontario, Canada, during the summer. Exposure was quantified with stochastic models at the population level, which incorporated measures of frequency, quantity ingested, prevalence, and concentration, using data from FoodNet Canada surveillance, the peer-reviewed and gray literature, other Ontario data, and data that were specifically collected for this study. Models were run with @Risk software using Monte Carlo simulations. The mean number of cells of Campylobacter ingested per Ontarian per day during the summer, ranked from highest to lowest is as follows: household pets, chicken, living on a farm, raw milk, visiting a farm, recreational water, beef, drinking water, pork, vegetables, seafood, petting zoos, and fruits. The study results identify knowledge gaps for some transmission routes, and indicate that some transmission routes for Campylobacter are underestimated in the current literature, such as household pets and raw milk. Many data gaps were identified for future data collection consideration, especially for the concentration of Campylobacter in all transmission routes.

Evidence for the benefits of food chain interventions on E. coli O157:H7/NM prevalence in retail ground beef and human disease incidence: A success story

OBJECTIVES

Human infection with Escherichia coli O157:H7/NM has historically been associated with consumption of undercooked ground beef. The purpose of this paper is to investigate the correlation of the decline in E. coli O157:H7/NM infections in Canada with the introduction of control efforts in ground beef by industry.

METHODS

The human incidence of E. coli O157:H7/NM, prevalence in ground beef and interventions from 1996 to 2014 were analyzed. Pathogen prevalence data were obtained from federal government and industry surveillance and inspection/compliance programs. A survey of the largest ground beef producers in Canada was conducted to identify when interventions were implemented.

RESULTS

The incidence of E. coli O157:H7/NM infections in Canada declined from ≈4 cases/100 000 to ≈1 case/100000 from 2000 to 2010. Verotoxigenic Escherichia coli (VTEC) prevalence in ground beef sold at retail declined from about 30% around the year 2000 to <2% since 2012. Other measures of the prevalence of E. coli, VTEC, and E. coli O157:H7/NM in beef and ground beef also declined. The number and types of interventions implemented in the major beef processing establishments in Canada increased from 1996 to 2016.

CONCLUSION

The observed decline in human illnesses and pathogen levels in relation to retail meats was associated with the introduction of control efforts by industry, federal and provincial/territorial governments, and the general population. Industry-led changes in beef processing along with the introduction of food safety policies, regulations, and public education have led to improved food safety in Canada.

Campylobacter spp. Prevalence and Levels in Raw Milk: A Systematic Review and Meta-Analysis

Campylobacteriosis is the leading bacterial gastrointestinal disease internationally, contributing significantly to the enteric illness burden. Cases have been associated with the consumption of raw milk, a behavior that has garnered attention recently. Estimates of the prevalence and levels of Campylobacter spp. in raw milk are lacking, which hinders risk assessment attempts. This article is a systematic review and meta-analysis of reported prevalence and levels of zoonotic Campylobacter spp. in the raw milk of cows, goats, and sheep in Canada, the United States, Europe, Australia, and New Zealand. The relevant literature was reviewed, and trained reviewers examined the results for inclusion of articles in the meta-analysis. Relevant data (prevalence and/or level of Campylobacter in raw milk, country of origin, animal species, sample source, Campylobacter species identified, etc.) were extracted, and a meta-analysis was performed in Stata v. 12 (Metaprop command). The weighted mean prevalence of Campylobacter spp. in raw milk samples was 1.18%. Subgroup analyses were conducted to examine how prevalence varied by study characteristics, with the highest prevalence values in studies from the United Kingdom (by country, 6.4%), about cows (by animal species, 1.3%), and including samples taken from inline filters (by sample source, 1.75%) and in studies that included species that are not pathogenic to humans (by Campylobacter species, 1.14%). Two articles each included a single Campylobacter level, 0.16 ± 0.3 and approximately 0.047 most probable number per ml. Despite a relatively low prevalence, consumption of raw milk is inherently risky because no treatment has been used to inactivate pathogens. This potential risk further supports maintaining regulations to limit the sales of raw milk.

Use of Admail and a geographic information system to send surveys to target populations

INTRODUCTION

This paper briefly describes the use of Canada Post Unaddressed Admail and a geographic information system (GIS) for survey distribution to a specific target population in a large, sparsely populated geographic area, and the effectiveness of this approach.

METHODS

Surveys were sent as Unaddressed Admail via Canada Post to a target population of people living within 5 km of a wind turbine in southwestern Ontario.

RESULTS

The overall response rate from 8 wind farms (in 8 counties) was 8.1%.

CONCLUSION

This approach has the potential to save time and money, but low response rates are common, distribution is not precise and there is potential for selection bias. Despite these flaws, Unaddressed Admail is worth consideration for delivery of information, study-recruitment materials and surveys to rural, remote and specific target populations.

A Systematic Review and Meta-Analysis of the Campylobacter spp. Prevalence and Concentration in Household Pets and Petting Zoo Animals for Use in Exposure Assessments

Animal contact is a potential transmission route for campylobacteriosis, and both domestic household pet and petting zoo exposures have been identified as potential sources of exposure. Research has typically focussed on the prevalence, concentration, and transmission of zoonoses from farm animals to humans, yet there are gaps in our understanding of these factors among animals in contact with the public who don’t live on or visit farms. This study aims to quantify, through a systematic review and meta-analysis, the prevalence and concentration of Campylobacter carriage in household pets and petting zoo animals. Four databases were accessed for the systematic review (PubMed, CAB direct, ProQuest, and Web of Science) for papers published in English from 1992–2012, and studies were included if they examined the animal population of interest, assessed prevalence or concentration with fecal, hair coat, oral, or urine exposure routes (although only articles that examined fecal routes were found), and if the research was based in Canada, USA, Europe, Australia, and New Zealand. Studies were reviewed for qualitative synthesis and meta-analysis by two reviewers, compiled into a database, and relevant studies were used to create a weighted mean prevalence value. There were insufficient data to run a meta-analysis of concentration values, a noted study limitation. The mean prevalence of Campylobacter in petting zoo animals is 6.5% based on 7 studies, and in household pets the mean is 24.7% based on 34 studies. Our estimated concentration values were: 7.65x10³cfu/g for petting zoo animals, and 2.9x10⁵cfu/g for household pets. These results indicate that Campylobacter prevalence and concentration are lower in petting zoo animals compared with household pets and that both of these animal sources have a lower prevalence compared with farm animals that do not come into contact with the public. There is a lack of studies on Campylobacter in petting zoos and/or fair animals in Canada and abroad. Within this literature, knowledge gaps were identified, and include: a lack of concentration data reported in the literature for Campylobacter spp. in animal feces, a distinction between ill and diarrheic pets in the reported studies, noted differences in shedding and concentrations for various subtypes of Campylobacter, and consistent reporting between studies.