Skip to main content

The public health exposome and pregnancy-related mortality in the United States: a high-dimensional computational analysis



Racial inequities in maternal mortality in the U.S. continue to be stark.


The 2015–2018, 4-year total population, county-level, pregnancy-related mortality ratio (PRM; deaths per 100,000 live births; National Center for Health Statistics (NCHS), restricted use mortality file) was linked with the Public Health Exposome (PHE). Using data reduction techniques, 1591 variables were extracted from over 62,000 variables for use in this analysis, providing information on the relationships between PRM and the social, health and health care, natural, and built environments. Graph theoretical algorithms and Bayesian analysis were applied to PHE/PRM linked data to identify latent networks.


PHE variables most strongly correlated with total population PRM were years of potential life lost and overall life expectancy. Population-level indicators of PRM were overall poverty, smoking, lack of exercise, heat, and lack of adequate access to food.


In this high-dimensional analysis, overall life expectancy, poverty indicators, and health behaviors were found to be the strongest predictors of pregnancy-related mortality. This provides strong evidence that maternal death is part of a broader constellation of both similar and unique health behaviors, social determinants and environmental exposures as other causes of death.

Peer Review reports


Study Question: What can analyzing the exposome tell us about maternal mortality?

What’s already known: Several individual-level medical and social factors raise the risk for maternal mortality. Analyzing a large number of predictors simultaneously is computationally difficult

What this study adds: A high-dimensional analysis found county-level overall life expectancy, poverty indicators, and health behaviors to be the strongest predictors of pregnancy-related mortality when an nationwide analysis of the public health exposome was conducted.


Maternal mortality (MM) has become a topic of intense interest in recent years. Statistics suggesting an increase in maternal mortality in the U.S. at a time when most of the world was seeing a decline [1], as well as the stark racial divides [2], has created renewed interest. In 2020, the National Center for Health Statistics (NCHS) reported a national MM ratio (MMR, defined as deaths of women while pregnant or within 42 days of being pregnant, from any cause related to or aggravated by the pregnancy or its management (but not from accidental or incidental causes) at 17.4 per 100,000 live births for 2018, up from to 12.7 in 2007 [3], an increase of 37%. Similarly, the pregnancy-related mortality ratio (PRM, defined as maternal death during pregnancy or within 1 year from any cause related to or aggravated by the pregnancy or its management) has increased 238%, from 7.2 deaths per 100,000 live births in 1987 to 17.2 deaths per 100,000 live births in 2011–2015 [4]. Almost a third of pregnancy-related deaths occur during pregnancy, approximately half of all occur after the day of delivery, and the remainder occur on the day of delivery [5, 6]. Additionally, wide racial inequities in maternal death have been found, with non-Hispanic (NH) Black women being two to four times more likely to die from a pregnancy-related complication than NH White women [3, 4, 6].

While individuals’ biological, behavioral, and demographic characteristics occur within PRM and are inseparable from the social and economic context of the environments in which they live, efforts to understand maternal death are just beginning to extend beyond consideration of individual-level risk factors. The literature on place-based/contextual determinants of maternal mortality remains very sparse, especially in comparison to other reproductive health outcomes, largely due to the rarity of this outcome and unavailability of administrative data with geographic identifiers to achieve adequate sample sizes. Research on this outcome has had a historically narrow focus on clinical factors and medical interventions; a recent systematic review of social determinants of maternal morbidity and mortality concluded that “relatively little attention has been given to examining social determinants of health and maternal health outcomes in the US”, and identified only 5 studies from 1990 through 2018 that investigated associations between “area-level characteristics” and maternal mortality [7]. Ecological systems theory suggests that MM is affected by a range of factors at the societal, policy, community, interpersonal, and individual levels [8]. A recent analysis of county- and zip code-level predictors found that social factors such as income, education, and food access were most strongly associated with pregnancy-related mortality [9]. Most risk factors in that analysis were similar between Black and White PRM but the absolute effect was often stronger for the Black PRM.

Disentangling the web of causation underlying MM and MM disparities is a challenging problem, but may be better understood with an exposome lens. An exposome approach considers how cumulative and lifetime environmental exposures across different domains are related to MM (Fig. 1). The “public health exposome” (PHE) in particular allows for the simultaneous evaluation of individual characteristics and behaviors within the context of exposures from a number of different domains, including the natural, built, social and policy environments, and health and health care. Social determinants such as racial residential segregation in US communities and other features of structural racism in US governance, institutions, cultural norms and values [10], have been identified as fundamental causes of racial health inequities, including maternal health [11, 12]. While the built environment has not been directly studied with respect to maternal mortality, a number of factors have been found to affect pregnancy health [13,14,15,16], including lack of transportation, access/availability of healthy foods, access to quality health care, and safe physical activity spaces, undergirded by social determinants such as high rates of poverty and low educational attainment within communities. Exposure to chemicals in the natural environment, especially air pollutants such as particulate matter and nitrogen oxides, previously have been associated with maternal mortality in studies conducted in Asia and Africa [17, 18]. Thus, a model that allows for the simultaneous consideration of complex, real world exposures from the social, healthcare, built, natural, and policy environments and PRM is warranted.

Fig. 1
figure 1

Conceptual model for disparities in maternal mortality

While incidence of PRM is unacceptably high, maternal deaths are rare in absolute terms. Many state-wide reviews include less than 100 cases [19,20,21,22,23,24]. Analyzing a concept like the exposome – the totality of exposures – even at the ecological level requires methods that can handle thousands of exposures simultaneously. When combined with a small number of cases and the totality of interconnected, multi-level exposures, contexts, and pathways occurring across the life course that shape a person’s risk for death during reproduction, the analytic problem becomes particularly complicated. Complex computational methods provide an alternative methodology for examining large datasets, patterns and places of occurrence, and contributing factors to PRM. A main aim of this study was to apply a high-dimensional exposome-wide analysis to pinpoint novel patterns of associations of risk and protective factors of PRM across counties in the United States in an effort to identify and elucidate associations across the natural built, social and policy environments, health and healthcare. This analysis uses novel methods of data reduction, graph-theoretical algorithms, and Bayesian analysis to examine an extensive dataset with a wide conceptual and geographical range.


Outcome data

The primary outcome in this analysis was the county-level pregnancy-related mortality ratio (PRM; deaths per 100,000 live births). (The analysis uses county equivalents in states that have different administrative units; for brevity, “county” is used in this paper as a generic term.) All cases of pregnancy-related death were identified in the 2015–2018, restricted use, maternal mortality file obtained through the NCHS using International Classification of Disease-10th revision codes for underlying cause of death (A34, and all codes in the Pregnancy, Childbirth, and the Puerperium Classification [O00-O99]). The restricted use maternal mortality file includes death records with geographic identifiers for county of residence for every decedent in the United States and applies the 2018 vital records coding scheme for identifying maternal deaths. The revised coding scheme mitigates misclassification errors resulting from the adoption of a standardized pregnancy-status checkbox on revised death certificates [3]. To maximize rate stability, counts of deaths were summed over four years in each county where a minimum of 1000 live births occurred in the same time period (n = 1709; without this exclusion, a few counties with a low number of births had artifactually high rates; Table S1); no exclusions were made with respect to the number of maternal deaths (Fig. 2). Compared to the included counties, excluded counties with low numbers of births (Table S1) had smaller populations; somewhat lower median income (45,400 vs. 51,300 USD); and a lower percentage NH Black population (median 0.9% vs. 4.2%), consistent with rural counties largely located in the Midwest and West (Fig. 2). Counts of live births by county were aggregated from the NCHS restricted use natality files and summed from 2015–2018 in order to estimate county-level PRM as deaths per 100,000 live births.

Fig. 2
figure 2

Counties included in the analysis (N = 1709)


A large data repository [25] has been collected, collated and curated, providing information on the natural, built, social and policy environments and health and health care, incorporating data from local, national, and global sources (Table 1). Variables were geocoded at scales ranging from 1-km grid cell, to census tract, and county, areas, linked by cross-walks. The PHE data repository contains measures of both health and environmental exposures for 3,141 counties and county equivalents, and spans over 15 years (2003–2018). It currently stores over 62,000 variables that have been geocoded and harmonized at an annual and county level in order to provide the spatial–temporal, contextual environmental data that can be linked to residential addresses and used to analyze environmental context at the county level. The PHE is fully curated with meta-data and a searchable data dictionary. These data were restricted to 2014–2018 in order to coincide with the time frame of the PRM data used in this study.

Table 1 Public health exposome database

Computational analysis

An overview of the analysis flow is provided in Fig. 3. Idiosyncrasies of information capture and variable overlap tend to make data such as this rife with repeated and collinear measures, which can clog the analysis and obfuscate potential results. An exposome autocorrelation guideline [26] was therefore adopted by which variables may be assumed to provide the same information whenever their correlation is at least 0.9. Pearson’s coefficient was chosen because all variables passed the Wilk-Shapiro [27] normality test after correcting for multiple tests using the Benjamini and Hochberg false discovery rate method [28]. Graph theoretical screening algorithms were then applied as follows. A finite, simple, undirected graph was constructed with a vertex for every variable and with edges only between autocorrelates as just defined. Integer Linear Programming via the Python MIP Optimizer [29] was applied using default software settings to compute a minimum dominating set [30]. Vertices not in this set were removed, reducing the focus to 3,903 variables. Confidence limit and margin of error variables were also eliminated, leaving 1591 variables for subsequent analysis.

Fig. 3
figure 3

Overview of analysis structure

County-level total population, NH White, NH Black, and Hispanic PRM ratio variables were then linked to the PHE by Federal Information Processing System (FIPS) codes for state and county. To this dataset, scalable computational methods [31] were applied in both a supervised, targeted fashion and an unsupervised, data science approach.

Supervised analysis was focused exclusively on the four PRM ratio variables. Features for each of these variables were first restricted to county equivalents (counties, parishes, and boroughs) with 1000 or more live births for the respective demographic. This limited the analysis to 1709 county equivalents for total population, 1442 for NH White, 361 for NH Black, and 396 for Hispanic population (Fig. 2). These four feature sets were then used to compute four sorted lists of Pearson correlation coefficients, one list for each PRM ratio variable and all its PHE correlates. Each list was then subdivided into five sublists based on PHE domains of interest (natural, built, social, policy, and health). The result was 24 lists (five domain lists and one combined list for each of the four PRM ratio variables). Correlates with r > 0.15 from each list were then retained for detailed scrutiny. Domains are defined by data source. Thus, for example, the built domain may include variables about income as it relates to housing programs, while the health and health care domain may include variables from morbidity, mortality, and health behavior and health access data.

Unsupervised analysis was more complex and performed in an effort to discover unseen variable networks or relationships in which PRM may play a role. Paraclique analysis is resilient to noise, has high density guarantees, and is robust [32] across multiple parameter settings. It began with the construction of a complete, finite, simple, undirected graph whose vertices denoted PHE and PRM variables and whose edge weights recorded pairwise Pearson correlation coefficients computed over the 1709 county equivalents with 1000 or more total population live births. Completeness and edge weights were eliminated by thresholding. Data inhomogeneity ruled out the use of a single threshold, as is commonly done for example with spectral methods [33], and so absolute values of the correlation coefficient of 0.15, 0.2, 0.3, 0.5, 0.7, and 0.8 were chosen for study. Dense subgraphs in the form of paracliques [34] were then extracted using glom terms [35] ranging from one to three for each of these five graphs A paraclique was retained for further analysis whenever it contained at least one PRM ratio variable.

Both supervised and unsupervised analyses produced variable sets of manageable size, to which traditional but unscalable methods were then applied. Bayesian network analysis was performed using the bnlearn R package [36], with a greedy hill-climbing method employed, again with default settings, to search for a network with minimal Bayesian information criterion score. A sample Bayesian network for the total population PRM ratio variable and its top 20 correlates from each of the PHE domains is illustrated with the finite, simple, directed graph shown in Figure S1.


When total population PRM was considered across the PHE in the supervised analysis, the variable most strongly correlated was years of potential life lost (Pearson correlation coefficient r = 0.254), followed by an inverse correlation with overall life expectancy (r = -0.249). Other correlated indicators were smoking, overall poverty, heat, lack of exercise, and lack of adequate access to food (Table 2). When stratified by race/ethnicity, several similar factors emerged (Table 3). For PRM in the NH Black population, the most highly correlated variables were publicly supplied per capita water use (r = 0.206); an inverse association with births to unmarried, White NH women (r = -0.206); households headed by younger people (r = 0.194); lack of sleep (r = 0.192); poor physical health (r = 0.184); and mental health (0.164); and several indicators of poverty, heat, and sun exposure. The publicly supplied per capita water use correlation seemed to be due to a single parish in Louisiana, which had notably high water use and NH Black PRM; if eliminated from the analysis, the correlation was much lower (r < 0.15). PRM in the Hispanic population was most strongly correlated with lower income in households living in public housing (r = 0.167) birth rates to those aged 15–17 (r = 0.167), and proportion of residents with a high school degree (r = -0.157). The NH White population PRM was associated with overall life expectancy (r = -0.188), household income and poverty (r = 0.181), lack of physical activity (r = 0.180), lack of health insurance (r = 0.176), and population with less than a high school degree (r = 0.174).

Table 2 Factors most strongly associated with national 4-year total population pregnancy-related mortality ratio per 100,000 live births (2015–2018). Correlations are with respect to the relevant PRM ratio variable
Table 3 Factors from all PHE domains most strongly associated with national 4-year pregnancy-related mortality ratio per 100,000 live births (2015–2018), stratified by race/ethnicity. Correlations are with respect to the relevant PRM ratio variable

In domain-specific analyses of total population PRM (Table 4), an array of variables was found to be consistently correlated with overall PRM. For the social domain, various indicators of poverty, lack of access to food, and low educational level were among the strongest correlates. In the health domain, overall life expectancy, smoking (r = 0.246), lack of activity, and diabetes (r = 0.184) were associated with total population PRM. There was also a positive correlation with obesity (r = 0.179), firearm deaths (r = 0.174), rate of AIDS among men who have sex with men (r = 0.168), and an inverse correlation with heavy drinking (r = -0.168). For the built environment, healthy food environment (r = -0.228), several indicators of income related to housing programs, and percentage of mobile-homes (r = 0.155) were correlated with PRM. For the natural and policy environment, most correlations were below our cut-offs. The strongest correlations for the natural environment were all with extreme heat (r = 0.13–0.19). In the policy environment, the only correlations greater than r = 0.1 were with programs to address child poverty (free school lunch eligibility, r = 0.184). Domain-specific analyses of correlations with maternal racial/ethnic specific PRM are presented in supplemental tables S2-S5.

Table 4 Domain-specific factors most strongly associated with national 4-year total population pregnancy-related mortality ratio per 100,000 live births (2015–2018). Correlations are with respect to the relevant PRM ratio variable

In the unsupervised analysis, two paracliques (threshold 15, glom 3, n = 77 and n = 12 variables, Table S5) contained the overall and race/ethnic-specific PRM variables. They also contained several indicators of health care availability, crime arrests (mostly for minor crimes such as bookmaking, curfew, and numbers & lottery, but also embezzlement and offenses against children and adults), population indicators such as birth rates among adolescents, some factors related to water use, and poverty and income inequality. The paraclique containing Hispanic PRM contained variables related to poverty, public housing, and Medicare quality.

In the Bayesian analysis (Table 5), within the factors most correlated with PRM (r > 0.15), the upstream indicators of overall PRM was access to healthy food. Within the built environment, two factors were indicators of overall PRM: access to healthy food; in the health domain, years of potential life lost, proportion of adults reporting no leisure-time physical activity, heart attacks, and chlamydia were most highly correlated; in the natural environment, extreme heat days; and in the social environment, poverty and summary low SES had the highest correlations with overall PRM.

Table 5 Summary of direct predictors based on Bayesian network analysis, pregnancy-related mortality, by domain and race/ethnicity

Study limitations

A limitation of the data used in this study is the relatively small number of PRM cases. Only county equivalents with at least 1000 live births were considered in order to increase rate stability. However, this limited the number of counties that could be examined, particularly for race/ethnicity-specific analyses, and results were still vulnerable to outliers. The outcome could not be analyzed at a lower level, such as neighborhood, due to the rarity of the event and the need to harmonize across data sources. This analysis relied on death records alone to identify cases of PRM. In addition, the method used to identify cases of PRM was limited to information drawn from death records along (ICD-10 code for underlying cause of death). This differs from the CDC’s identification of PRM, which utilizes multiple sources of data and expert reviews in order to determine pregnancy-relatedness in each case. As such, comparisons between these measures cannot be made. Temporal or lagged effects were not considered, because of the approach used to determine PRM (combined 4-year ratio). Finally, while the PHE dataset includes over 62,000 variables to capture area-level features of built, policy, social, and natural environments, the PHE may still fail to capture the totality of potential forces driving PRM in this country. For instance, while a few measures of residential segregation and income inequality are among the PHE variables (but did not turn out to be strong predictors relative to other factors), other measures of structural racism such as racial inequities in education, income, homeownership, or criminal justice [37, 38]; or historical policies like redlining [39] are not included. The policies included in the exposome are relatively limited, and do not, for instance, include those around drugs, sentencing, abortion, or Medicaid eligibility. Causality cannot be explored within an ecological analysis, and inferences cannot be drawn about individual-level associations (it is not known, for example, whether women who died were in poverty or lacked health insurance).

Relying solely on correlation to determine autocorrelates may eliminate potentially highly correlated independent variables from the study. The reduction in autocorrelates, however, produced clearer results. When selecting top correlates to PRM, choosing the 20 most correlated variables may leave out potentially revealing associations. Twenty variables struck a reasonable balance between too many to visualize and too few for exploratory analysis. While dense subgraphs containing PRM via paraclique were found, there is no guarantee that all latent networks that include PRM were detected. This be due to the low overall correlation between PRM and the PHE, requiring a threshold of 0.15. The low cluster thresholds indicate highly subtle variable associations. Finally, all results in this study rely on the details of the methods used and different choices of parameters might affect the results (not unlike alpha = 0.05 or 0.10 for choosing statistical predictors). For this reason, when available, other than for paracliques, we used default software parameter settings for all methods used to aid reproducibility.


In this high-dimensional analysis of over 3900 potential predictors, years of potential life lost was found to be the single, most consistent predictor of total population PRM. The number of maternal deaths is very small relative to other contributors to total years of potential life lost [40], so reverse causality is an unlikely cause of this finding. Since maternal deaths are unusual in occurring in young, often otherwise healthy women, it is not necessarily obvious that they should correlate with overall indicators of mortality, which are driven primarily by older sick people. Reproductive and pregnancy health is often treated as separate, or secondary, to other aspects of health, but this analysis suggests that the major risk factors driving maternal mortality are the same as those driving health more generally. Other related factors by race/ethnicity included overall life expectancy, poverty, smoking, lack of physical activity, and lack of access to healthy food. This paints a clear picture of maternal death as part and parcel of an adverse, social environmental context, rather than caused by a limited, proximal set of risk factors. Maternal deaths are highest in areas where life expectancy is already low. While poverty is widely associated with health [41], it is still striking that the cluster of poverty-related factors rather than potentially more proximal factors such as obstetric services or hypertension rates remain the strongest predictors even after examining thousands of other exposures. The complete list of variables for the built environment domain, for instance, started with about 8700 variables and includes factors like land use; commute time; crowding and density, and rurality; housing age and quality; heating type; amount of housing units, size, and availability; and number of private nonfarm and farm establishments and their characteristics. However, the strongest correlations with PRM within that domain described the social characteristics of the people living in that built environment. Many of the variables in the PHE are related, interconnected, and may be linked to mortality via multiple, multilevel, and complex pathways. In such a case identifying individual risk factors may be less useful than capturing the broader context, represented by multiple indicators, which was accomplished in this analysis.

Race/ethnicity-specific analyses confirmed this overall pattern of interconnected social factors and general health indicators being the strongest predictors of PRM, with identified factors indicating socioeconomic status (education, unemployment) or general population health (adults in fair/poor health, overall life expectancy) as the most important correlates with PRM. While individual variables differed by race/ethnicity, overall, they were part of a larger constellation representing social deprivation (low education, low income, unemployment, lack of health insurance). This finding supports the notion that racism, not race, underlies racial inequities in maternal mortality [11] as with other population health inequities [42, 43]. Historical and contemporary racism restricts access to health-promoting resources such as income, wealth, and education, which results in disproportionate morbidity and mortality [44, 45]. Finally, a few specific health conditions of the population also were identified (diabetes, obesity, heart attacks), which is consistent with the high proportion of pregnancy-related deaths from underlying cardiovascular disease conditions [46].

The graph theoretical methods used in this analysis have several advantages, and we examine our results to determine whether such advantages were realized. First, these methods scale to extreme numbers of variables without issues of multiple comparisons. In this case, several thousand variables and their associations were studied simultaneously. Second, the methods used to elucidate intercorrelated subgroups (paracliques) can provide new insights when compared to traditional methods that can only study associations between variable pairs. In this analysis, various latent associations were thereby identified. Third, these methods allow for the discovery of novel associations or entire groups of factors that are highly intercorrelated. While it is not a surprise to see associations of PRM with lack of physical activity, unhealthy food environment, or low mental health, our first hypotheses might have concerned access to obstetric care or underlying medical causes such as hypertension instead. Thus, the methods employed in the current study fulfilled the goal of identifying novel associations. Finally, use of these methods in the PHE specifically allowed for finding joint effects of the physical and social environments. While most associations identified were from the social and health care domains, some physical environment factors, particularly heat, were also identified as important.

By comparison, our previous work used hypothesis-driven modeling, and found associations between maternal mortality and income inequality [47], levels of violent crime [48], and access to maternity care [49]. These more standard analyses concentrated on specific factors of interest with control for covariates, a model that better supports causal inference. That research, however, was limited to specific social factors, analyzed one at a time. Attempting this for thousands of predictors would be plagued by multiple comparisons, while examining all possible variable combinations would quickly become computationally intractable. The graph theoretical methods we have employed are often best viewed as hypothesis generators. With ecological data especially, they can provide a solid basis for more detailed study incorporating multi-level data, confounders, and mechanistic analysis.

This is an exploratory, ecological study of aggregate-level population data. The results are broadly consistent with a recent hypothesis-based investigation of county-level predictors of PRM in the U.S. which found that county-level unemployment and food insecurity were risk factors and higher income, education, and owner-occupied household were protective factors [9]. Many public health measures are administered at the county level, and for many of the variables included in the Public Health Exposome, existing literature identifies the county level as a relevant geographic resolution for their association with health outcomes [50, 51]. The processes underlying PRM are multi-level, and our focus on the county level does not negate the importance of considering state or more localized levels as well. Indeed, existing evidence identifies state-level contexts, structures, and policies that contribute to the vastly different PRM rates that occur across states [47], but these analyses mask the variation occurring within states. The analysis also pulled out the strongest correlations, not necessarily those that are causal. For instance, food environment could be causally related to maternal death, or it could simply be strongly correlated with income and other health-promoting resources. Other factors, such as rate of HIV among men who have sex with men, are likely a combination of social indicator and indicator of related health factors, such as access to contraception and sexual health services. Most of the data used are national, so should be measured consistently across counties. Ecological studies can be particularly vulnerable to cross-level confounding and interactions [52, 53]; the results of our study are consistent with what is known about individual-level associations, so we think the more important concern is to address multi-level influences and the interactions among effects occurring at different levels.

Strengths of this study include the use of a nationwide, 2015–2018, maternal mortality data set and the wide-ranging environmental and health risk factors represented by the PHE that were considered in the analysis. As can be seen in the tables, the most predictive factors were drawn from a range of databases, and did not necessarily center on a single one, highlighting the value of the full exposome database. If all the predictive results came from a single data source, we should more effectively and efficiently focus our efforts on that source. As the factors we found to be most important are drawn from several of the dozens of databases included in the PHE, it is serving its purpose of describing the range of exposures that may be influencing public health.


Recent years have seen a heightened concern about pregnancy-related mortality in general and disparities in PRM in particular. The examination of 1591 potential associated variables in our analyses did not find unusual or novel risk factors; rather, it confirms that this issue is not separate from the United States’ socioeconomic inequalities and general health problems, but is embedded in them. Study findings contribute to the growing body of evidence documenting the impact of adverse social and environmental factors on maternal death [47, 54, 55]. From a methodological standpoint, we continue to develop “big data” methods that can deal with social and environmental factors as rigorously and completely as the methods for addressing genomics or transcriptomics data. Scalable combinatorial methods have achieved impressive successes in the analysis of high throughput, high fidelity -omics data. In domains such as social science and public health, however, their reach is sometimes limited by data quality and type. Prohibitive levels of noise can occur in many forms, including data inaccuracies, collinearities, conflicts, and uneven sampling. This motivates the sustained development of robust, noise-resilient methods such as the paraclique algorithm and effective feature selection and preprocessing techniques such as MDS-90. From a content standpoint, efforts to reduce maternal mortality must expand beyond individual and clinical risk factors to include the effect of harmful social and environmental contexts, with a particular focus on the experience of women from different racial and ethnic groups, and those who are socio-economically disadvantaged. Further investigation of the intricate links between environmental exposures and maternal death is crucial to elucidate the mechanisms and pathways through which contextual factors contribute to the detrimental health effects of exposed individuals.

Availability of data and materials

Maternal mortality data were acquired from NCHS ( The public health exposome is compiled from publicly-available data sources and maintained by Meharry Medical College Health Disparities Research Center of Excellence. Paul Juarez can be contacted for more information (



Maternal mortality


Maternal mortality ratio


National Center for Health Statistics




Public health exposome


Pregnancy-related mortality


Years of potential life lost


  1. Creanga AA. Maternal mortality in the united states: a review of contemporary data and their limitations. Clin Obstet Gynecol. 2018;61(2):296–306.

    Article  PubMed  Google Scholar 

  2. Jain JA, Temming LA, D’Alton ME, Gyamfi-Bannerman C, Tuuli M, Louis JM, Srinivas SK, Caughey AB, Grobman WA, Hehir M, et al. SMFM Special Report: Putting the “M” back in MFM: Reducing racial and ethnic disparities in maternal morbidity and mortality: A call to action. Am J Obstet Gynecol. 2018;218(2):B9-b17.

    Article  PubMed  Google Scholar 

  3. Hoyert DL, Miniño AM. Maternal mortality in the United States: changes in coding, publication, and data release, 2018. Natl Vital Stat Reports. 2020;69(2):1–18.

    Google Scholar 

  4. Petersen EE, Davis NL, Goodman D, Cox S, Mayes N, Johnston E, Syverson C, Seed K, Shapiro-Mendoza CK, Callaghan WM, et al. Vital signs: pregnancy-RELATED Deaths, United States, 2011–2015, and Strategies for Prevention, 13 States, 2013–2017. MMWR Morb Mortal Wkly Rep. 2019;68(18):423–9.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Creanga AA, Berg CJ, Syverson C, Seed K, Bruce FC, Callaghan WM. Pregnancy-related mortality in the United States, 2006–2010. Obstet Gynecol. 2015;125(1):5–12.

    Article  PubMed  Google Scholar 

  6. Creanga AA, Syverson C, Seed K, Callaghan WM. Pregnancy-Related Mortality in the United States, 2011–2013. Obstet Gynecol. 2017;130(2):366–73.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Wang E, Glazer KB, Howell EA, Janevic TM. Social determinants of pregnancy-related mortality and morbidity in the United States: a systematic review. Obstet Gynecol. 2020;135(4):896–915.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Noursi S, Saluja B, Richey L. Using the ecological systems theory to understand black/white disparities in maternal morbidity and mortality in the United States. J Racial Ethn Health Disparities. 2021;8(3):661–9.

    Article  PubMed  Google Scholar 

  9. Barrera CM, Kramer MR, Merkt PT, Petersen EE, Brantley MD, Eckhaus L, et al. County-level associations between pregnancy-related mortality ratios and contextual sociospatial indicators. Obstet Gynecol. 2022;139(5):855–65.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Krieger N, Waterman PD, Spasojevic J, Li W, Maduro G, Van Wye G. Public health monitoring of privilege and deprivation with the index of concentration at the extremes. Am J Public Health. 2016;106(2):256–63.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Crear-Perry J, Correa-de-Araujo R, Lewis Johnson T, McLemore MR, Neilson E, Wallace M. Social and Structural Determinants of Health Inequities in Maternal Health. J Women’s Health (2002). 2021;30(2):230–5.

  12. Williams DR, Lawrence JA, Davis BA. Racism and health: evidence and needed research. Annu Rev Public Health. 2019;40:105–25.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Woods N, Gilliland J, Seabrook JA. The influence of the built environment on adverse birth outcomes. Neonatal-Perinatal Med. 2017;10(3):233–48.

    Article  CAS  Google Scholar 

  14. Porter AK, Rodríguez DA, Frizzelle BG, Evenson KR. The Association between Neighborhood Environments and Physical Activity from Pregnancy to Postpartum: a Prospective Cohort Study. J Urban Health. 2019;96(5):703–19.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Nieuwenhuijsen MJ, Agier L, Basagaña X, Urquiza J, Tamayo-Uria I, Giorgis-Allemand L, Robinson O, Siroux V, Maitre L, de Castro M, et al. Influence of the Urban Exposome on Birth Weight. Environ Health Perspect. 2019;127(4):47007.

    Article  PubMed  Google Scholar 

  16. Bastain TM, Chavez T, Habre R, Girguis MS, Grubbs B, Toledo-Corral C, Amadeus M, Farzan SF, Al-Marayati L, Lerner D, et al. Study Design, Protocol and Profile of the Maternal And Developmental Risks from Environmental and Social Stressors (MADRES) Pregnancy Cohort: a Prospective Cohort Study in Predominantly Low-Income Hispanic Women in Urban Los Angeles. BMC Pregnancy Childbirth. 2019;19(1):189.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Lien WH, Owili PO, Muga MA, Lin TH. Ambient particulate matter exposure and under-five and maternal deaths in Asia. Int J Environ Res Public Health. 2019;16(20):3855.

    Article  CAS  PubMed Central  Google Scholar 

  18. Owili PO, Lien WH, Muga MA, Lin TH. The associations between types of ambient PM(2.5) and under-five and maternal mortality in Africa. Int J Environ Res Public Health. 2017;14(4):359.

  19. Baeva S, Saxton DL, Ruggiero K, Kormondy ML, Hollier LM, Hellerstedt J, Hall M, Archer NP. Identifying Maternal Deaths in Texas Using an Enhanced Method, 2012. Obstet Gynecol. 2018;131(5):762–9.

    Article  PubMed  Google Scholar 

  20. Burch D, Noell D, Hill WC, Delke I. Pregnancy-associated mortality review: the Florida experience. Semin Perinatol. 2012;36(1):31–6.

    Article  PubMed  Google Scholar 

  21. Geller SE, Koch AR, Martin NJ, Rosenberg D, Bigger HR. Assessing preventability of maternal mortality in Illinois: 2002–2012. Am J Obstet Gynecol. 2014;211(6):698.e691-611.

    Article  Google Scholar 

  22. Hernandez LE, Sappenfield WM, Harris K, Burch D, Hill WC, Clark CL, Delke I. Pregnancy-Related Deaths, Florida, 1999–2012: Opportunities to Improve Maternal Outcomes. Matern Child Health J. 2018;22(2):204–15.

    Article  PubMed  Google Scholar 

  23. Schellpfeffer MA, Gillespie KH, Rohan AM, Blackwell SP. A Review of Pregnancy-Related Maternal Mortality in Wisconsin, 2006–2010. WMJ. 2015;114(5):202–7.

    PubMed  Google Scholar 

  24. Shellhaas C, Conrey E. State-based Review of Maternal Deaths: The Ohio Experience. Clin Obstet Gynecol. 2018;61(2):332–9.

    Article  PubMed  Google Scholar 

  25. Juarez PD, Matthews-Juarez P, Hood DB, Im W, Levine RS, Kilbourne BJ, Langston MA, Al-Hamdan MZ, Crosson WL, Estes MG, et al. The Public Health Exposome: A Population-Based, Exposure Science Approach to Health Disparities Research. Int J Environ Res Public Health. 2014;11:12866–95.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Santos S, Maitre L, Warembourg C, Agier L, Richiardi L, Basagaña X, Vrijheid M. Applying the exposome concept in birth cohort research: a review of statistical approaches. Eur J Epidemiol. 2020;35(3):193–204.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Shapiro SS, Wilk MB. An analysis of variance test for normality (complete samples). Biometrika. 1965;52(3–4):591–611.

    Article  Google Scholar 

  28. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J Roy Stat Soc: Ser B (Methodol). 1995;57(1):289–300.

    Google Scholar 

  29. Túlio A. M. Toffolo HGS: Python MIP.; 2021

  30. Garey MR, Johnson DS. Computers and intractability: a guide to the theory of NP-Completeness. W. H. Freeman and Company; 1979.

  31. Langston MA, Levine RS, Kilbourne BJ, Rogers GL, Kershenbaum AD, Baktash SH, Coughlin SS, Saxton AM, Agboto VA, Hood DB, et al. Scalable Combinatorial Tools for Health Disparities Research. Int J Environ Res Public Health. 2014;11(10):10419–43.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Lu Y, Phillips CA, Langston MA. A robustness metric for biological data clustering algorithms. BMC Bioinformatics. 2019;20(15):503.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Perkins AD, Langston MA: Threshold Selection in Gene Co-Expression Networks Using Spectral Graph Theory Techniques. BMC Bioinformatics 2009, 10.

  34. Chesler EJ, Langston MA. Combinatorial Genetic Regulatory Network Analysis Tools for High Throughput Transcriptomic Data. In: Eskin E, editor. Systems Biology and Regulatory Genomics. New York City: Springer; 2006. p. 150–65 Volume 4023.

    Chapter  Google Scholar 

  35. Hagan RD, Langston MA, Wang K. Lower Bounds on Paraclique Density. Discret Appl Math. 2016;204:208–12.

    Article  Google Scholar 

  36. Scutari M. Learning Bayesian Networks with the bnlearn R Package. J Stat Softw. 2010;35:1–22.

    Article  Google Scholar 

  37. Alson JG, Robinson WR, Pittman L, Doll KM. Incorporating Measures of Structural Racism into Population Studies of Reproductive Health in the United States: A Narrative Review. Health Equity. 2021;5(1):49–58.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Agénor M, Perkins C, Stamoulis C, Hall RD, Samnaliev M, Berland S, et al. Developing a database of structural racism-related state laws for health equity research and practice in the United States. Public Health Rep. 2021;136(4):428–40. Washington, DC: 1974.

  39. Dennis AC, Chung EO, Lodge EK, Martinez RA, Wilbur RE. Looking Back to Leap Forward: A Framework for Operationalizing the Structural Racism Construct in Minority Health Research. Ethn Dis. 2021;31(Suppl 1):301–10.

    Article  PubMed  PubMed Central  Google Scholar 

  40. National Center for Injury Prevention and Control. Years of Potential Life Lost (YPLL) before age 65. Atlanta: Centers for Disease Control and Prevention; 2019.

  41. Chokshi DA. Income, Poverty, and Health Inequality. JAMA. 2018;319(13):1312–3.

    Article  PubMed  Google Scholar 

  42. Bailey ZD, Feldman JM, Bassett MT. How Structural Racism Works - Racist Policies as a Root Cause of U.S. Racial Health Inequities. N Engl J Med. 2021;384(8):768–73.

    Article  PubMed  Google Scholar 

  43. Cooper RS. Race in biological and biomedical research. Cold Spring Harb Perspect Med. 2013;3(11):a008573.

    Article  PubMed  PubMed Central  Google Scholar 

  44. National Academies of Sciences E, and Medicine, Health and Medicine Division, Board on Population Health and Public Health Practice, Committee on Community-Based Solutions to Promote Health Equity in the United States. The state of health disparities in the United States. In: Baciu A, Negussie YAG, editors. Communities in action: pathways to health equity, Vol. 2. Washington, DC: National Academies Press; 2017.

  45. Minehart RD, Bryant AS, Jackson J, Daly JL. Racial/Ethnic Inequities in Pregnancy-Related Morbidity and Mortality. Obstet Gynecol Clin North Am. 2021;48(1):31–51.

    Article  PubMed  Google Scholar 

  46. Varagic J, Desvigne-Nickens P, Gamble-George J, Hollier L, Maric-Bilkan C, Mitchell M, et al. Maternal morbidity and mortality: are we getting to the “Heart” of the matter? J Womens Health (Larchmt). 2021;30(2):178–86.

    Article  Google Scholar 

  47. Vilda D, Wallace M, Dyer L, Harville E, Theall K. Income inequality and racial disparities in pregnancy-related mortality in the US. SSM - Popul Health. 2019;9:100477.

  48. Wallace ME, Friar N, Herwehe J, Theall KP. Violence As a Direct Cause of and Indirect Contributor to Maternal Death. J Women’s Health (2002). 2020;29(8):1032–8.

    Article  Google Scholar 

  49. Wallace M, Dyer L, Felker-Kantor E, Benno J, Vilda D, Harville E, et al. Maternity care deserts and pregnancy-associated mortality in Louisiana. Women’s Health Issues. 2021;31(2):122–9.

    Article  PubMed  Google Scholar 

  50. Juarez PD, Hood DB, Song MA, Ramesh A. Use of an Exposome Approach to Understand the Effects of Exposures From the Natural, Built, and Social Environments on Cardio-Vascular Disease Onset, Progression, and Outcomes. Front Public Health. 2020;8:379.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Cifuentes P, Reichard J, Im W, Smith S, Colen C, Giurgescu C, Williams KP, Gillespie S, Juarez PD, Hood DB. Application of the Public Health Exposome Framework to Estimate Phenotypes of Resilience in a Model Ohio African-American Women’s Cohort. J Urban Health. 2019;96(Suppl 1):57–71.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Blakely TA, Woodward AJ. Ecological effects in multi-level studies. J Epidemiol Community Health. 2000;54(5):367–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Greenland S. Ecologic versus individual-level sources of bias in ecologic estimates of contextual health effects. Int J Epidemiol. 2001;30(6):1343–50.

    Article  CAS  PubMed  Google Scholar 

  54. Eliason EL. Adoption of Medicaid Expansion Is Associated with Lower Maternal Mortality. Womens Health Issues. 2020;30(3):147–52.

    Article  PubMed  Google Scholar 

  55. Nelson DB, Moniz MH, Davis MM. Population-level factors associated with maternal mortality in the United States, 1997–2012. BMC Public Health. 2018;18(1):1007.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


This research has been supported in part by the National Institute of Child Health and Human Development under grant R01HD092653-03S1. It has used resources of the National Energy Research Scientific Computing Center (NERSC), a U.S. Department of Energy Office of Science User Facility located at Lawrence Berkeley National Laboratory, operated under Contract No. DE-AC02-05CH11231.

Author information

Authors and Affiliations



EWH conceived and managed the project, and drafted the manuscript. SKG conducted many of the computational analyses and wrote much of methods section. MAL supervised computational analyses and provided access to computational resources. PJ developed the PHE and advised on its interpretation. DV organized the MM data and advised on background and interpretation of social factors. MEW supervised the parent grant, acquired data, and advised on content. All authors read, edited, and approved the final manuscript.

Corresponding author

Correspondence to E. W. Harville.

Ethics declarations

Ethics approval and consent to participate

This analysis of de-identified aggregate data was ruled exempt by the Tulane Institutional Review Board. All methods were carried out in accordance with relevant guidelines and regulations.

Consent for publication

NA. This analysis of de-identified aggregate data was ruled exempt by the Tulane Institutional Review Board and so informed consent was not required.

Competing interests

None of the authors have a conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Harville, E.W., Grady, S., Langston, M. et al. The public health exposome and pregnancy-related mortality in the United States: a high-dimensional computational analysis. BMC Public Health 22, 2097 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: