Casey Olives

An Integrated Risk Function for Estimating the Global Burden of Disease Attributable to Ambient Fine Particulate Matter Exposure

Background: Estimating the burden of disease attributable to long-term exposure to fine particulate matter (PM2.5) in ambient air requires knowledge of both the shape and magnitude of the relative risk (RR) function. However, adequate direct evidence to identify the shape of the mortality RR functions at the high ambient concentrations observed in many places in the world is lacking. Objective: We developed RR functions over the entire global exposure range for causes of mortality in adults: ischemic heart disease (IHD), cerebrovascular disease (stroke), chronic obstructive pulmonary disease (COPD), and lung cancer (LC). We also developed RR functions for the incidence of acute lower respiratory infection (ALRI) that can be used to estimate mortality and lost-years of healthy life in children < 5 years of age. Methods: We fit an integrated exposure–response (IER) model by integrating available RR information from studies of ambient air pollution (AAP), second hand tobacco smoke, household solid cooking fuel, and active smoking (AS). AS exposures were converted to estimated annual PM2.5 exposure equivalents using inhaled doses of particle mass. We derived population attributable fractions (PAFs) for every country based on estimated worldwide ambient PM2.5 concentrations. Results: The IER model was a superior predictor of RR compared with seven other forms previously used in burden assessments. The percent PAF attributable to AAP exposure varied among countries from 2 to 41 for IHD, 1 to 43 for stroke, < 1 to 21 for COPD, < 1 to 25 for LC, and < 1 to 38 for ALRI. Conclusions: We developed a fine particulate mass–based RR model that covered the global range of exposure by integrating RR information from different combustion types that generate emissions of particulate matter. The model can be updated as new RR information becomes available. Citation: Burnett RT, Pope CA III, Ezzati M, Olives C, Lim SS, Mehta S, Shin HH, Singh G, Hubbell B, Brauer M, Anderson HR, Smith KR, Balmes JR, Bruce NG, Kan H, Laden F, Prüss-Ustün A, Turner MC, Gapstur SM, Diver WR, Cohen A. 2014. An integrated risk function for estimating the global burden of disease attributable to ambient fine particulate matter exposure. Environ Health Perspect 122:397–403; http://dx.doi.org/10.1289/ehp.1307049

Prevalence, Awareness, Treatment, and Control of Hypertension in United States Counties, 2001–2009

Hypertension is an important and modifiable risk factor for cardiovascular disease and mortality. Over the last decade, national-levels of controlled hypertension have increased, but little information on hypertension prevalence and trends in hypertension treatment and control exists at the county-level. We estimate trends in prevalence, awareness, treatment, and control of hypertension in US counties using data from the National Health and Nutrition Examination Survey (NHANES) in five two-year waves from 1999–2008 including 26,349 adults aged 30 years and older and from the Behavioral Risk Factor Surveillance System (BRFSS) from 1997–2009 including 1,283,722 adults aged 30 years and older. Hypertension was defined as systolic blood pressure (BP) of at least 140 mm Hg, self-reported use of antihypertensive treatment, or both. Hypertension control was defined as systolic BP less than 140 mm Hg. The median prevalence of total hypertension in 2009 was estimated at 37.6% (range: 26.5 to 54.4%) in men and 40.1% (range: 28.5 to 57.9%) in women. Within-state differences in the county prevalence of uncontrolled hypertension were as high as 7.8 percentage points in 2009. Awareness, treatment, and control was highest in the southeastern US, and increased between 2001 and 2009 on average. The median county-level control in men was 57.7% (range: 43.4 to 65.9%) and in women was 57.1% (range: 43.0 to 65.46%) in 2009, with highest rates in white men and black women. While control of hypertension is on the rise, prevalence of total hypertension continues to increase in the US. Concurrent increases in treatment and control of hypertension are promising, but efforts to decrease the prevalence of hypertension are needed.

Association between air pollution and coronary artery calcification within six metropolitan areas in the USA (the Multi-Ethnic Study of Atherosclerosis and Air Pollution): a longitudinal cohort study

BACKGROUND Long-term exposure to fine particulate matter less than 2.5 μm in diameter (PM2.5) and traffic-related air pollutant concentrations are associated with cardiovascular risk. The disease process underlying these associations remains uncertain. We aim to assess association between long-term exposure to ambient air pollution and progression of coronary artery calcium and common carotid artery intima-media thickness. METHODS In this prospective 10-year cohort study, we repeatedly measured coronary artery calcium by CT in 6795 participants aged 45-84 years enrolled in the Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air) in six metropolitan areas in the USA. Repeated scans were done for nearly all participants between 2002 and 2005, for a subset of participants between 2005 and 2007, and for half of all participants between 2010 and 2012. Common carotid artery intima-media thickness was measured by ultrasound in all participants at baseline and in 2010-12 for 3459 participants. Residence-specific spatio-temporal pollution concentration models, incorporating community-specific measurements, agency monitoring data, and geographical predictors, estimated concentrations of PM2.5 and nitrogen oxides (NOX) between 1999 and 2012. The primary aim was to examine the association between both progression of coronary artery calcium and mean carotid artery intima-media thickness and long-term exposure to ambient air pollutant concentrations (PM2.5, NOX, and black carbon) between examinations and within the six metropolitan areas, adjusting for baseline age, sex, ethnicity, socioeconomic characteristics, cardiovascular risk factors, site, and CT scanner technology. FINDINGS In this population, coronary calcium increased on average by 24 Agatston units per year (SD 58), and intima-media thickness by 12 μm per year (10), before adjusting for risk factors or air pollutant exposures. Participant-specific pollutant concentrations averaged over the years 2000-10 ranged from 9.2-22.6 μg PM2.5/m(3) and 7.2-139.2 parts per billion (ppb) NOX. For each 5 μg PM2.5/m(3) increase, coronary calcium progressed by 4.1 Agatston units per year (95% CI 1.4-6.8) and for each 40 ppb NOX coronary calcium progressed by 4.8 Agatston units per year (0.9-8.7). Pollutant exposures were not associated with intima-media thickness change. The estimate for the effect of a 5 μg/m(3) higher long-term exposure to PM2.5 in intima-media thickness was -0.9 μm per year (95% CI -3.0 to 1.3). For 40 ppb higher NOX, the estimate was 0.2 μm per year (-1.9 to 2.4). INTERPRETATION Increased concentrations of PM2.5 and traffic-related air pollution within metropolitan areas, in ranges commonly encountered worldwide, are associated with progression in coronary calcification, consistent with acceleration of atherosclerosis. This study supports the case for global efforts of pollution reduction in prevention of cardiovascular diseases. FUNDING US Environmental Protection Agency and US National Institutes of Health.

A unified spatiotemporal modeling approach for predicting concentrations of multiple air pollutants in the multi-ethnic study of atherosclerosis and air pollution.

Background: Cohort studies of the relationship between air pollution exposure and chronic health effects require predictions of exposure over long periods of time. Objectives: We developed a unified modeling approach for predicting fine particulate matter, nitrogen dioxide, oxides of nitrogen, and black carbon (as measured by light absorption coefficient) in six U.S. metropolitan regions from 1999 through early 2012 as part of the Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air). Methods: We obtained monitoring data from regulatory networks and supplemented those data with study-specific measurements collected from MESA Air community locations and participants’ homes. In each region, we applied a spatiotemporal model that included a long-term spatial mean, time trends with spatially varying coefficients, and a spatiotemporal residual. The mean structure was derived from a large set of geographic covariates that was reduced using partial least-squares regression. We estimated time trends from observed time series and used spatial smoothing methods to borrow strength between observations. Results: Prediction accuracy was high for most models, with cross-validation R2 (R2CV) > 0.80 at regulatory and fixed sites for most regions and pollutants. At home sites, overall R2CV ranged from 0.45 to 0.92, and temporally adjusted R2CV ranged from 0.23 to 0.92. Conclusions: This novel spatiotemporal modeling approach provides accurate fine-scale predictions in multiple regions for four pollutants. We have generated participant-specific predictions for MESA Air to investigate health effects of long-term air pollution exposures. These successes highlight modeling advances that can be adopted more widely in modern cohort studies. Citation: Keller JP, Olives C, Kim SY, Sheppard L, Sampson PD, Szpiro AA, Oron AP, Lindström J, Vedal S, Kaufman JD. 2015. A unified spatiotemporal modeling approach for predicting concentrations of multiple air pollutants in the Multi-Ethnic Study of Atherosclerosis and Air Pollution. Environ Health Perspect 123:301–309; http://dx.doi.org/10.1289/ehp.1408145

Risk Factors for Long-Term Coronary Artery Calcium Progression in the Multi-Ethnic Study of Atherosclerosis

Background Coronary artery calcium (CAC) detected by noncontrast cardiac computed tomography scanning is a measure of coronary atherosclerosis burden. Increasing CAC levels have been strongly associated with increased coronary events. Prior studies of cardiovascular disease risk factors and CAC progression have been limited by short follow-up or restricted to patients with advanced disease. Methods and Results We examined cardiovascular disease risk factors and CAC progression in a prospective multiethnic cohort study. CAC was measured 1 to 4 times (mean 2.5 scans) over 10 years in 6810 adults without preexisting cardiovascular disease. Mean CAC progression was 23.9 Agatston units/year. An innovative application of mixed-effects models investigated associations between cardiovascular disease risk factors and CAC progression. This approach adjusted for time-varying factors, was flexible with respect to follow-up time and number of observations per participant, and allowed simultaneous control of factors associated with both baseline CAC and CAC progression. Models included age, sex, study site, scanner type, and race/ethnicity. Associations were observed between CAC progression and age (14.2 Agatston units/year per 10 years [95% CI 13.0 to 15.5]), male sex (17.8 Agatston units/year [95% CI 15.3 to 20.3]), hypertension (13.8 Agatston units/year [95% CI 11.2 to 16.5]), diabetes (31.3 Agatston units/year [95% CI 27.4 to 35.3]), and other factors. Conclusions CAC progression analyzed over 10 years of follow-up, with a novel analytical approach, demonstrated strong relationships with risk factors for incident cardiovascular events. Longitudinal CAC progression analyzed in this framework can be used to evaluate novel cardiovascular risk factors.

Bayes-LQAS: Classifying the Prevalence of Global Acute Malnutrition

Lot Quality Assurance Sampling (LQAS) applications in health have generally relied on frequentist interpretations for statistical validity. Yet health professionals often seek statements about the probability distribution of unknown parameters to answer questions of interest. The frequentist paradigm does not pretend to yield such information, although a Bayesian formulation might. This is the source of an error made in a recent paper published in this journal. Many applications lend themselves to a Bayesian treatment, and would benefit from such considerations in their design. We discuss Bayes-LQAS (B-LQAS), which allows for incorporation of prior information into the LQAS classification procedure, and thus shows how to correct the aforementioned error. Further, we pay special attention to the formulation of Bayes Operating Characteristic Curves and the use of prior information to improve survey designs. As a motivating example, we discuss the classification of Global Acute Malnutrition prevalence and draw parallels between the Bayes and classical classifications schemes. We also illustrate the impact of informative and non-informative priors on the survey design. Results indicate that using a Bayesian approach allows the incorporation of expert information and/or historical data and is thus potentially a valuable tool for making accurate and precise classifications.

Multiple Category-Lot Quality Assurance Sampling: A New Classification System with Application to Schistosomiasis Control

Background Originally a binary classifier, Lot Quality Assurance Sampling (LQAS) has proven to be a useful tool for classification of the prevalence of Schistosoma mansoni into multiple categories (≤10%, >10 and <50%, ≥50%), and semi-curtailed sampling has been shown to effectively reduce the number of observations needed to reach a decision. To date the statistical underpinnings for Multiple Category-LQAS (MC-LQAS) have not received full treatment. We explore the analytical properties of MC-LQAS, and validate its use for the classification of S. mansoni prevalence in multiple settings in East Africa. Methodology We outline MC-LQAS design principles and formulae for operating characteristic curves. In addition, we derive the average sample number for MC-LQAS when utilizing semi-curtailed sampling and introduce curtailed sampling in this setting. We also assess the performance of MC-LQAS designs with maximum sample sizes of n = 15 and n = 25 via a weighted kappa-statistic using S. mansoni data collected in 388 schools from four studies in East Africa. Principle Findings Overall performance of MC-LQAS classification was high (kappa-statistic of 0.87). In three of the studies, the kappa-statistic for a design with n = 15 was greater than 0.75. In the fourth study, where these designs performed poorly (kappa-statistic less than 0.50), the majority of observations fell in regions where potential error is known to be high. Employment of semi-curtailed and curtailed sampling further reduced the sample size by as many as 0.5 and 3.5 observations per school, respectively, without increasing classification error. Conclusion/Significance This work provides the needed analytics to understand the properties of MC-LQAS for assessing the prevalance of S. mansoni and shows that in most settings a sample size of 15 children provides a reliable classification of schools.

A class of non-linear exposure-response models suitable for health impact assessment applicable to large cohort studies of ambient air pollution

The effectiveness of regulatory actions designed to improve air quality is often assessed by predicting changes in public health resulting from their implementation. Risk of premature mortality from long-term exposure to ambient air pollution is the single most important contributor to such assessments and is estimated from observational studies generally assuming a log-linear, no-threshold association between ambient concentrations and death. There has been only limited assessment of this assumption in part because of a lack of methods to estimate the shape of the exposure-response function in very large study populations. In this paper, we propose a new class of variable coefficient risk functions capable of capturing a variety of potentially non-linear associations which are suitable for health impact assessment. We construct the class by defining transformations of concentration as the product of either a linear or log-linear function of concentration multiplied by a logistic weighting function. These risk functions can be estimated using hazard regression survival models with currently available computer software and can accommodate large population-based cohorts which are increasingly being used for this purpose. We illustrate our modeling approach with two large cohort studies of long-term concentrations of ambient air pollution and mortality: the American Cancer Society Cancer Prevention Study II (CPS II) cohort and the Canadian Census Health and Environment Cohort (CanCHEC). We then estimate the number of deaths attributable to changes in fine particulate matter concentrations over the 2000 to 2010 time period in both Canada and the USA using both linear and non-linear hazard function models.

Historical Prediction Modeling Approach for Estimating Long-Term Concentrations of PM2.5 in Cohort Studies before the 1999 Implementation of Widespread Monitoring.

INTRODUCTION Recent cohort studies have used exposure prediction models to estimate the association between long-term residential concentrations of fine particulate matter (PM2.5) and health. Because these prediction models rely on PM2.5 monitoring data, predictions for times before extensive spatial monitoring present a challenge to understanding long-term exposure effects. The U.S. Environmental Protection Agency (EPA) Federal Reference Method (FRM) network for PM2.5 was established in 1999. OBJECTIVES We evaluated a novel statistical approach to produce high-quality exposure predictions from 1980 through 2010 in the continental United States for epidemiological applications. METHODS We developed spatio-temporal prediction models using geographic predictors and annual average PM2.5 data from 1999 through 2010 from the FRM and the Interagency Monitoring of Protected Visual Environments (IMPROVE) networks. Temporal trends before 1999 were estimated by using a) extrapolation based on PM2.5 data in FRM/IMPROVE, b) PM2.5 sulfate data in the Clean Air Status and Trends Network, and c) visibility data across the Weather Bureau Army Navy network. We validated the models using PM2.5 data collected before 1999 from IMPROVE, California Air Resources Board dichotomous sampler monitoring (CARB dichot), the Childrens Health Study (CHS), and the Inhalable Particulate Network (IPN). RESULTS In our validation using pre-1999 data, the prediction model performed well across three trend estimation approaches when validated using IMPROVE and CHS data (R2 = 0.84-0.91) with lower R2 values in early years. Model performance using CARB dichot and IPN data was worse (R2 = 0.00-0.85) most likely because of fewer monitoring sites and inconsistent sampling methods. CONCLUSIONS Our prediction modeling approach will allow health effects estimation associated with long-term exposures to PM2.5 over extended time periods ≤ 30 years. Citation: Kim SY, Olives C, Sheppard L, Sampson PD, Larson TV, Keller JP, Kaufman JD. 2017. Historical prediction modeling approach for estimating long-term concentrations of PM2.5 in cohort studies before the 1999 implementation of widespread monitoring. Environ Health Perspect 125:38-46; http://dx.doi.org/10.1289/EHP131.

Cluster designs to assess the prevalence of acute malnutrition by lot quality assurance sampling: a validation study by computer simulation

Traditional lot quality assurance sampling (LQAS) methods require simple random sampling to guarantee valid results. However, cluster sampling has been proposed to reduce the number of random starting points. This study uses simulations to examine the classification error of two such designs, a 67×3 (67 clusters of three observations) and a 33×6 (33 clusters of six observations) sampling scheme to assess the prevalence of global acute malnutrition (GAM). Further, we explore the use of a 67×3 sequential sampling scheme for LQAS classification of GAM prevalence. Results indicate that, for independent clusters with moderate intracluster correlation for the GAM outcome, the three sampling designs maintain approximate validity for LQAS analysis. Sequential sampling can substantially reduce the average sample size that is required for data collection. The presence of intercluster correlation can impact dramatically the classification error that is associated with LQAS analysis.