Measuring long-term exposure to wildfire PM2.5 in California: Time-varying inequities in environmental burden.

Wildfires have become more frequent and intense due to climate change and outdoor wildfire fine particulate matter (PM2.5) concentrations differ from relatively smoothly varying total PM2.5. Thus, we introduced a conceptual model for computing long-term wildfire PM2.5 and assessed disproportionate exposures among marginalized communities. We used monitoring data and statistical techniques to characterize annual wildfire PM2.5 exposure based on intermittent and extreme daily wildfire PM2.5 concentrations in California census tracts (2006 to 2020). Metrics included: 1) weeks with wildfire PM2.5 < 5 µg/m3; 2) days with non-zero wildfire PM2.5; 3) mean wildfire PM2.5 during peak exposure week; 4) smoke waves (≥2 consecutive days with <15 µg/m3 wildfire PM2.5); and 5) mean annual wildfire PM2.5 concentration. We classified tracts by their racial/ethnic composition and CalEnviroScreen (CES) score, an environmental and social vulnerability composite measure. We examined associations of CES and racial/ethnic composition with the wildfire PM2.5 metrics using mixed-effects models. Averaged 2006 to 2020, we detected little difference in exposure by CES score or racial/ethnic composition, except for non-Hispanic American Indian and Alaska Native populations, where a 1-SD increase was associated with higher exposure for 4/5 metrics. CES or racial/ethnic × year interaction term models revealed exposure disparities in some years. Compared to their California-wide representation, the exposed populations of non-Hispanic American Indian and Alaska Native (1.68×, 95% CI: 1.01 to 2.81), white (1.13×, 95% CI: 0.99 to 1.32), and multiracial (1.06×, 95% CI: 0.97 to 1.23) people were over-represented from 2006 to 2020. In conclusion, during our study period in California, we detected disproportionate long-term wildfire PM2.5 exposure for several racial/ethnic groups.

Air Pollution and Mortality at the Intersection of Race and Social Class.

BACKGROUND: Black Americans are exposed to higher annual levels of air pollution containing fine particulate matter (particles with an aerodynamic diameter of ≤2.5 µm [PM2.5]) than White Americans and may be more susceptible to its health effects. Low-income Americans may also be more susceptible to PM2.5 pollution than high-income Americans. Because information is lacking on exposure-response curves for PM2.5 exposure and mortality among marginalized subpopulations categorized according to both race and socioeconomic position, the Environmental Protection Agency lacks important evidence to inform its regulatory rulemaking for PM2.5 standards. METHODS: We analyzed 623 million person-years of Medicare data from 73 million persons 65 years of age or older from 2000 through 2016 to estimate associations between annual PM2.5 exposure and mortality in subpopulations defined simultaneously by racial identity (Black vs. White) and income level (Medicaid eligible vs. ineligible). RESULTS: Lower PM2.5 exposure was associated with lower mortality in the full population, but marginalized subpopulations appeared to benefit more as PM2.5 levels decreased. For example, the hazard ratio associated with decreasing PM2.5 from 12 µg per cubic meter to 8 µg per cubic meter for the White higher-income subpopulation was 0.963 (95% confidence interval [CI], 0.955 to 0.970), whereas equivalent hazard ratios for marginalized subpopulations were lower: 0.931 (95% CI, 0.909 to 0.953) for the Black higher-income subpopulation, 0.940 (95% CI, 0.931 to 0.948) for the White low-income subpopulation, and 0.939 (95% CI, 0.921 to 0.957) for the Black low-income subpopulation. CONCLUSIONS: Higher-income Black persons, low-income White persons, and low-income Black persons may benefit more from lower PM2.5 levels than higher-income White persons. These findings underscore the importance of considering racial identity and income together when assessing health inequities. (Funded by the National Institutes of Health and the Alfred P. Sloan Foundation.).

Meta-analysis of breast cancer risk for individuals with PALB2 pathogenic variants.

Multigene panel testing now allows efficient testing of many cancer susceptibility genes leading to a larger number of mutation carriers being identified. They need to be counseled about their cancer risk conferred by the specific gene mutation. An important cancer susceptibility gene is PALB2. Multiple studies reported risk estimates for breast cancer (BC) conferred by pathogenic variants in PALB2. Due to the diverse modalities of reported risk estimates (age-specific risk, odds ratio, relative risk, and standardized incidence ratio) and effect sizes, a meta-analysis combining these estimates is necessary to accurately counsel patients with this mutation. However, this is not trivial due to heterogeneity of studies in terms of study design and risk measure. We utilized a recently proposed Bayesian random-effects meta-analysis method that can synthesize estimates from such heterogeneous studies. We applied this method to combine estimates from 12 studies on BC risk for carriers of pathogenic PALB2 mutations. The estimated overall (meta-analysis-based) risk of BC is 12.80% (6.11%-22.59%) by age 50 and 48.47% (36.05%-61.74%) by age 80. Pathogenic mutations in PALB2 makes women more susceptible to BC. Our risk estimates can help clinically manage patients carrying pathogenic variants in PALB2.

Comparing traditional and causal inference methodologies for evaluating impacts of long-term air pollution exposure on hospitalization with Alzheimer's disease and related dementias.

Alzheimer's disease and related dementias (ADRD) present a growing public health burden in the United States. One actionable risk factor for ADRD is air pollution: multiple studies have found associations between air pollution and exacerbation of ADRD. Our study builds on previous studies by applying modern statistical causal inference methodologies-generalized propensity score (GPS) weighting and matching-on a large, longitudinal dataset. We follow 50 million Medicare enrollees to investigate impacts of three air pollutants-fine particular matter (PM${}_{2.5}$), nitrogen dioxide (NO${}_2$), and summer ozone (O${}_3$)-on elderly patients' rate of first hospitalization with ADRD diagnosis. Similar to previous studies using traditional statistical models, our results found increased hospitalization risks due to increased PM${}_{2.5}$ and NO${}_2$ exposure, with less conclusive results for O${}_3$. In particular, our GPS weighting analysis finds IQR increases in PM${}_{2.5}$, NO${}_2$, or O${}_3$ exposure results in hazard ratios of 1.108 (95% CI: 1.097-1.119), 1.058 (1.049-1.067), or 1.045 (1.036-1.054), respectively. GPS matching results are similar for PM${}_{2.5}$ and NO${}_2$ with attenuated effects for O${}_3$. Our results strengthen arguments that long-term PM${}_{2.5}$ and NO${}_2$ exposure increases risk of hospitalization with ADRD diagnosis. Additionally, we highlight strengths and limitations of causal inference methodologies in observational studies with continuous treatments. Keywords: Alzheimer's disease and related dementias, air pollution, Medicare, causal inference, generalized propensity score.

Assessing the causal effects of a stochastic intervention in time series data: are heat alerts effective in preventing deaths and hospitalizations?

The methodological development of this article is motivated by the need to address the following scientific question: does the issuance of heat alerts prevent adverse health effects? Our goal is to address this question within a causal inference framework in the context of time series data. A key challenge is that causal inference methods require the overlap assumption to hold: each unit (i.e., a day) must have a positive probability of receiving the treatment (i.e., issuing a heat alert on that day). In our motivating example, the overlap assumption is often violated: the probability of issuing a heat alert on a cooler day is near zero. To overcome this challenge, we propose a stochastic intervention for time series data which is implemented via an incremental time-varying propensity score (ItvPS). The ItvPS intervention is executed by multiplying the probability of issuing a heat alert on day $t$-conditional on past information up to day $t$-by an odds ratio $\delta_t$. First, we introduce a new class of causal estimands, which relies on the ItvPS intervention. We provide theoretical results to show that these causal estimands can be identified and estimated under a weaker version of the overlap assumption. Second, we propose nonparametric estimators based on the ItvPS and derive an upper bound for the variances of these estimators. Third, we extend this framework to multisite time series using a spatial meta-analysis approach. Fourth, we show that the proposed estimators perform well in terms of bias and root mean squared error via simulations. Finally, we apply our proposed approach to estimate the causal effects of increasing the probability of issuing heat alerts on each warm-season day in reducing deaths and hospitalizations among Medicare enrollees in 2837 US counties.

Rejoinder to the discussion on "Bayesian meta-analysis of penetrance for cancer risk".

The five discussions of our paper provide several modeling alternatives, extensions, and generalizations that can potentially guide future research in meta-analysis. In this rejoinder, we briefly summarize and comment on some of those points.

Bayesian meta-analysis of penetrance for cancer risk.

Multi-gene panel testing allows many cancer susceptibility genes to be tested quickly at a lower cost making such testing accessible to a broader population. Thus, more patients carrying pathogenic germline mutations in various cancer-susceptibility genes are being identified. This creates a great opportunity, as well as an urgent need, to counsel these patients about appropriate risk-reducing management strategies. Counseling hinges on accurate estimates of age-specific risks of developing various cancers associated with mutations in a specific gene, ie, penetrance estimation. We propose a meta-analysis approach based on a Bayesian hierarchical random-effects model to obtain penetrance estimates by integrating studies reporting different types of risk measures (eg, penetrance, relative risk, odds ratio) while accounting for the associated uncertainties. After estimating posterior distributions of the parameters via a Markov chain Monte Carlo algorithm, we estimate penetrance and credible intervals. We investigate the proposed method and compare with an existing approach via simulations based on studies reporting risks for two moderate-risk breast cancer susceptibility genes, ATM and PALB2. Our proposed method is far superior in terms of coverage probability of credible intervals and mean square error of estimates. Finally, we apply our method to estimate the penetrance of breast cancer among carriers of pathogenic mutations in the ATM gene.

Statistical methods for Mendelian models with multiple genes and cancers.

Risk evaluation to identify individuals who are at greater risk of cancer as a result of heritable pathogenic variants is a valuable component of individualized clinical management. Using principles of Mendelian genetics, Bayesian probability theory, and variant-specific knowledge, Mendelian models derive the probability of carrying a pathogenic variant and developing cancer in the future, based on family history. Existing Mendelian models are widely employed, but are generally limited to specific genes and syndromes. However, the upsurge of multigene panel germline testing has spurred the discovery of many new gene-cancer associations that are not presently accounted for in these models. We have developed PanelPRO, a flexible, efficient Mendelian risk prediction framework that can incorporate an arbitrary number of genes and cancers, overcoming the computational challenges that arise because of the increased model complexity. We implement an 11-gene, 11-cancer model, the largest Mendelian model created thus far, based on this framework. Using simulations and a clinical cohort with germline panel testing data, we evaluate model performance, validate the reverse-compatibility of our approach with existing Mendelian models, and illustrate its usage. Our implementation is freely available for research use in the PanelPRO R package.

A pedigree-based prediction model identifies carriers of deleterious de novo mutations in families with Li-Fraumeni syndrome.

De novo mutations (DNMs) are increasingly recognized as rare disease causal factors. Identifying DNM carriers will allow researchers to study the likely distinct molecular mechanisms of DNMs. We developed Famdenovo to predict DNM status (DNM or familial mutation [FM]) of deleterious autosomal dominant germline mutations for any syndrome. We introduce Famdenovo.TP53 for Li-Fraumeni syndrome (LFS) and analyze 324 LFS family pedigrees from four US cohorts: a validation set of 186 pedigrees and a discovery set of 138 pedigrees. The concordance index for Famdenovo.TP53 prediction was 0.95 (95% CI: [0.92, 0.98]). Forty individuals (95% CI: [30, 50]) were predicted as DNM carriers, increasing the total number from 42 to 82. We compared clinical and biological features of FM versus DNM carriers: (1) cancer and mutation spectra along with parental ages were similarly distributed; (2) ascertainment criteria like early-onset breast cancer (age 20-35 yr) provides a condition for an unbiased estimate of the DNM rate: 48% (23 DNMs vs. 25 FMs); and (3) hotspot mutation R248W was not observed in DNMs, although it was as prevalent as hotspot mutation R248Q in FMs. Furthermore, we introduce Famdenovo.BRCA for hereditary breast and ovarian cancer syndrome and apply it to a small set of family data from the Cancer Genetics Network. In summary, we introduce a novel statistical approach to systematically evaluate deleterious DNMs in inherited cancer syndromes. Our approach may serve as a foundation for future studies evaluating how new deleterious mutations can be established in the germline, such as those in TP53.

A Literature Review of the Effects of Air Pollution on COVID-19 Health Outcomes Worldwide: Statistical Challenges and Data Visualization.

Several peer-reviewed papers and reviews have examined the relationship between exposure to air pollution and COVID-19 spread and severity. However, many of the existing reviews on this topic do not extensively present the statistical challenges associated with this field, do not provide comprehensive guidelines for future researchers, and review only the results of a relatively small number of papers. We reviewed 139 papers, 127 of which reported a statistically significant positive association between air pollution and adverse COVID-19 health outcomes. Here, we summarize the evidence, describe the statistical challenges, and make recommendations for future research. To summarize the 139 papers with data from geographical locations around the world, we also present anopen-source data visualization tool that summarizes these studies and allows the research community to contribute evidence as new research papers are published.

Clinical implications of conflicting variant interpretations in the cancer genetics clinic.

PURPOSE: The aim of this study was to describe the clinical impact of commercial laboratories issuing conflicting classifications of genetic variants. METHODS: Results from 2000 patients undergoing a multigene hereditary cancer panel by a single laboratory were analyzed. Clinically significant discrepancies between the laboratory-provided test reports and other major commercial laboratories were identified, including differences between pathogenic/likely pathogenic and variant of uncertain significance (VUS) classifications, via review of ClinVar archives. For patients carrying a VUS, clinical documentation was assessed for evidence of provider awareness of the conflict. RESULTS: Fifty of 975 (5.1%) patients with non-negative results carried a variant with a clinically significant conflict, 19 with a pathogenic/likely pathogenic variant reported in APC or MUTYH, and 31 with a VUS reported in CDKN2A, CHEK2, MLH1, MSH2, MUTYH, RAD51C, or TP53. Only 10 of 28 (36%) patients with a VUS with a clinically significant conflict had a documented discussion by a provider about the conflict. Discrepant counseling strategies were used for different patients with the same variant. Among patients with a CDKN2A variant or a monoallelic MUTYH variant, providers were significantly more likely to make recommendations based on the laboratory-reported classification. CONCLUSION: Our findings highlight the frequency of variant interpretation discrepancies and importance of clinician awareness. Guidance is needed on managing patients with discrepant variants to support accurate risk assessment.

Investigating Use of Low-Cost Sensors to Increase Accuracy and Equity of Real-Time Air Quality Information.

U.S. Environmental Protection Agency (EPA) air quality (AQ) monitors, the "gold standard" for measuring air pollutants, are sparsely positioned across the U.S. Low-cost sensors (LCS) are increasingly being used by the public to fill in the gaps in AQ monitoring; however, LCS are not as accurate as EPA monitors. In this work, we investigate factors impacting the differences between an individual's true (unobserved) exposure to air pollution and the exposure reported by their nearest AQ instrument (which could be either an LCS or an EPA monitor). We use simulations based on California data to explore different combinations of hypothetical LCS placement strategies (e.g., at schools or near major roads), for different numbers of LCS, with varying plausible amounts of LCS device measurement errors. We illustrate how real-time AQ reporting could be improved (or, in some cases, worsened) by using LCS, both for the population overall and for marginalized communities specifically. This work has implications for the integration of LCS into real-time AQ reporting platforms.

Evaluation of Model-Based PM2.5 Estimates for Exposure Assessment during Wildfire Smoke Episodes in the Western U.S.

Investigating the health impacts of wildfire smoke requires data on people's exposure to fine particulate matter (PM2.5) across space and time. In recent years, it has become common to use machine learning models to fill gaps in monitoring data. However, it remains unclear how well these models are able to capture spikes in PM2.5 during and across wildfire events. Here, we evaluate the accuracy of two sets of high-coverage and high-resolution machine learning-derived PM2.5 data sets created by Di et al. and Reid et al. In general, the Reid estimates are more accurate than the Di estimates when compared to independent validation data from mobile smoke monitors deployed by the US Forest Service. However, both models tend to severely under-predict PM2.5 on high-pollution days. Our findings complement other recent studies calling for increased air pollution monitoring in the western US and support the inclusion of wildfire-specific monitoring observations and predictor variables in model-based estimates of PM2.5. Lastly, we call for more rigorous error quantification of machine-learning derived exposure data sets, with special attention to extreme events.

Statistical approaches for meta-analysis of genetic mutation prevalence.

Estimating the prevalence of rare germline genetic mutations in the general population is of interest as it can inform genetic counseling and risk management. Most studies that estimate the prevalence of mutations are performed in high-risk populations, and each study is designed with differing inclusion criteria, resulting in ascertained populations. Quantifying the effects of ascertainment is necessary to estimate the prevalence in the general population. This quantification is difficult as the inclusion criteria is often based on disease status and/or family history. Combining estimates from multiple studies through a meta-analysis is challenging due to the variety of study designs and ascertainment mechanisms as well as the complexity of quantifying the effect of these mechanisms. We provide guidelines on how to quantify the ascertainment mechanism for a wide range of settings and propose a general approach for conducting a meta-analysis in these complex settings by incorporating study-specific ascertainment mechanisms into a joint likelihood function. We implement the proposed likelihood-based approach using both frequentist and Bayesian methodologies. We evaluate these approaches in simulations and show that the methods are robust and produce unbiased estimates of the prevalence. An advantage of the Bayesian approach is that it can easily incorporate uncertainty in ascertainment probability values. We apply our methods to estimate the prevalence of PALB2 mutations in the United States by combining data from multiple studies and obtain a prevalence estimate of around 0.02%.

Variation in cancer risk among families with genetic susceptibility.

Germline mutations in many genes have been shown to increase the risk of developing cancer. This risk can vary across families who carry mutations in the same gene due to differences in the specific variants, gene-gene interactions, other susceptibility mutations, environmental factors, and behavioral factors. We develop an analytic tool to explore this heterogeneity using family history data. We propose to evaluate the ratio between the number of observed cancer cases in a family and the number of expected cases under a model where risk is assumed to be the same across families. We perform this analysis for both carriers and noncarriers in each family, using carrier probabilities when carrier statuses are unknown, and visualize the results. We first illustrate the approach in simulated data and then apply it to data on colorectal cancer risk in families carrying mutations in Lynch syndrome genes from Creighton University's Hereditary Cancer Center. We show that colorectal cancer risk in carriers can vary widely across families, and that this variation is not matched by a corresponding variation in the noncarriers from the same families. This suggests that the sources of variation in these families are to be found predominantly in variants harbored in the mutated MMR genes considered, or in variants interacting with them.

Penetrance of male breast cancer susceptibility genes: a systematic review.

PURPOSE: Several male breast cancer (MBC) susceptibility genes have been identified, but the MBC risk for individuals with a pathogenic variant in each of these genes (i.e., penetrance) remains unclear. We conducted a systematic review of studies reporting the penetrance of MBC susceptibility genes to better summarize current estimates of penetrance. METHODS: A search query was developed to identify MBC-related papers indexed in PubMed/MEDLINE. A validated natural language processing method was applied to identify papers reporting penetrance estimates. These penetrance studies' bibliographies were reviewed to ensure comprehensiveness. We accessed the potential ascertainment bias for each enrolled study. RESULTS: Fifteen penetrance studies were identified from 12,182 abstracts, covering five purported MBC susceptibility genes: ATM, BRCA1, BRCA2, CHEK2, and PALB2. Cohort (n = 6, 40%) and case-control (n = 5, 33%) studies were the two most common study designs, followed by family-based (n = 3, 20%), and a kin-cohort study (n = 1, 7%). Seven of the 15 studies (47%) adjusted for ascertainment adequately and therefore the MBC risks reported by these seven studies can be considered applicable to the general population. Based on these seven studies, we found pathogenic variants in ATM, BRCA2, CHEK2 c.1100delC, and PALB2 show an increased risk for MBC. The association between BRCA1 and MBC was not statistically significant. CONCLUSION: This work supports the conclusion that pathogenic variants in ATM, BRCA2, CHEK2 c.1100delC, and PALB2 increase the risk of MBC, whereas pathogenic variants in BRCA1 may not be associated with increased MBC risk.

A validation of models for prediction of pathogenic variants in mismatch repair genes.

PURPOSE: Models used to predict the probability of an individual having a pathogenic homozygous or heterozygous variant in a mismatch repair gene, such as MMRpro, are widely used. Recently, MMRpro was updated with new colorectal cancer penetrance estimates. The purpose of this study was to evaluate the predictive performance of MMRpro and other models for individuals with a family history of colorectal cancer. METHODS: We performed a validation study of 4 models, Leiden, MMRpredict, PREMM5, and MMRpro, using 784 members of clinic-based families from the United States. Predicted probabilities were compared with germline testing results and evaluated for discrimination, calibration, and predictive accuracy. We analyzed several strategies to combine models and improve predictive performance. RESULTS: MMRpro with additional tumor information (MMRpro+) and PREMM5 outperformed the other models in discrimination and predictive accuracy. MMRpro+ was the best calibrated with an observed to expected ratio of 0.98 (95% CI = 0.89-1.08). The combination models showed improvement over PREMM5 and performed similar to MMRpro+. CONCLUSION: MMRpro+ and PREMM5 performed well in predicting the probability of having a pathogenic homozygous or heterozygous variant in a mismatch repair gene. They serve as useful clinical decision tools for identifying individuals who would benefit greatly from screening and prevention strategies.

A Bayesian Multi-Outcome Analysis of Fine Particulate Matter and Cardiorespiratory Hospitalizations.

BACKGROUND: Short-term fine particulate matter (PM2.5) exposure is positively associated with acute cardiovascular and respiratory events. Understanding whether this association varies across specific cardiovascular and respiratory conditions has important biologic, clinical, and public health implications. METHODS: We conducted a time-stratified case-crossover study of hospitalizations from 2000 through 2014 among United States Medicare beneficiaries aged 65+. The outcomes were hospitalizations with any of 57 cardiovascular and 32 respiratory discharge diagnoses. We estimated associations with two-day moving average PM2.5 as a piecewise linear term with a knot at PM2.5 = 25 g/m3. We used Multi-Outcome Regression with Tree-structured Shrinkage (MOReTreeS) to identify de novo groups of related diseases such that PM2.5 associations are: (1) similar within outcome groups; but (2) different between outcome groups. We adjusted for temperature, humidity, and individual-level characteristics. We introduce an R package, moretrees. RESULTS: Our dataset included 16,007,293 cardiovascular and 8,690,837 respiratory hospitalizations. Of 57 cardiovascular diseases, 51 were grouped and positively associated with PM2.5. We observed a stronger positive association for heart failure, which formed a separate group. We observed negative associations for groups containing the outcomes other aneurysm and intracranial hemorrhage. Of 32 respiratory outcomes, 31 were grouped and were positively associated with PM2.5. Influenza formed a separate group with a negative association. CONCLUSIONS: We used a new statistical approach, MOReTreeS, to uncover variation in the association between short-term PM2.5 exposure and hospitalizations for cardiovascular and respiratory causes controlling for patient characteristics, time trends, and environmental confounders.

Practical implementation of frailty models in Mendelian risk prediction.

There are numerous statistical models used to identify individuals at high risk of cancer due to inherited mutations. Mendelian models predict future risk of cancer by using family history with estimated cancer penetrances (age- and sex-specific risk of cancer given the genotype of the mutations) and mutation prevalences. However, there is often residual risk heterogeneity across families even after accounting for the mutations in the model, due to environmental or unobserved genetic risk factors. We aim to improve Mendelian risk prediction by incorporating a frailty model that contains a family-specific frailty vector, impacting the cancer hazard function, to account for this heterogeneity. We use a discrete uniform population frailty distribution and implement a marginalized approach that averages each family's risk predictions over the family's frailty distribution. We apply the proposed approach to improve breast cancer prediction in BRCAPRO, a Mendelian model that accounts for inherited mutations in the BRCA1 and BRCA2 genes to predict breast and ovarian cancer. We evaluate the proposed model's performance in simulations and real data from the Cancer Genetics Network and show improvements in model calibration and discrimination. We also discuss alternative approaches for incorporating frailties and their strengths and limitations.

Nationwide Study of Short-term Exposure to Fine Particulate Matter and Cardiovascular Hospitalizations Among Medicaid Enrollees.

BACKGROUND: Fine particulate matter (PM2.5) has been consistently linked to cardiovascular disease (CVD). Although studies have reported modification by income, to our knowledge, no study to date has examined this relationship among adults in Medicaid, which provides health coverage to low-income and/or disabled Americans. METHODS: We estimated the association between short-term PM2.5 exposure (average of PM2.5 on the day of hospitalization and the preceding day) and CVD admissions rates among adult Medicaid enrollees in the continental United States (2000-2012) using a time-stratified case-crossover design. We repeated this analysis at PM2.5 concentrations below the World Health Organization daily guideline of 25 µg/m. We compared the PM2.5-CVD association in the Medicaid ≥65 years old versus non-Medicaid-eligible Medicare enrollees (≥65 years old). RESULTS: Using information on 3,666,657 CVD hospitalizations among Medicaid adults, we observed a 0.9% (95% CI = 0.6%, 1.1%) increase in CVD admission rates per 10 µg/m PM2.5 increase. The association was stronger at low PM2.5 levels (1.3%; 95% CI = 0.9%, 1.6%). Among Medicaid enrollees ≥65 years old, the association was 0.9% (95% CI = 0.6%, 1.3%) vs. 0.8% (95% CI = 0.6%, 0.9%) among non-Medicaid-eligible Medicare enrollees ≥65 years old. CONCLUSION: We found robust evidence of an association between short-term PM2.5 and CVD hospitalizations among the vulnerable subpopulation of adult Medicaid enrollees. Importantly, this association persisted even at PM2.5 levels below the current national standards.