Abstract
Importance Seasonal influenza hospitalizations pose a considerable burden in the United States, with BIPOC (Black, Indigenous, and other People of Color) communities being disproportionately affected.
Objective To determine and quantify the effects of different types of mitigation strategies on inequities in influenza outcomes (symptomatic infections and hospitalizations).
Design In this simulation study, we fit a race-stratified agent-based model of influenza transmission to demographic and hospitalization data of the United States.
Participants We consider five racial-ethnic groups: non-Hispanic White persons, non-Hispanic Black persons, non-Hispanic Asian persons, non-Hispanic American Indian or Alaska Native persons, and Hispanic or Latino persons.
Setting We tested five idealized equity-promoting interventions to determine their effectiveness in reducing inequity in influenza outcomes. The interventions assumed (i) equalized vaccination rates, (ii) equalized comorbidities, (iii) work-risk distribution proportional to the distribution of the population, (iv) reduced work contacts for all, or (v) a combination of equalizing vaccination rates and comorbidities and reducing work contacts.
Main Outcomes and Measures Reduction in symptomatic or hospitalization risk ratios, defined as the ratio of the number of symptomatic infections (hospitalizations respectively) in each age- and racial-ethnic group and their corresponding white counterpart. We also evaluated the reduction in the absolute mean number of symptomatic infections or hospitalizations in each age- and racial-ethnic group compared to the fitted scenario (baseline).
Results Our analysis suggests that symptomatic infections were equalized and reduced (by up to 17% in BIPOC adults aged 18-49) by strategies reducing work contacts or equalizing vaccination rates. Reducing comorbidities resulted in significant decreases in hospitalizations, with a reduction of over 40% in BIPOC groups. All tested interventions reduced the inequity in influenza hospitalizations in all racial-ethnic groups, but interventions reducing comorbidities in marginalized populations were the most effective. Notably, these interventions resulted in better outcomes across all racial-ethnic groups, not only those prioritized by the interventions.
Conclusions and Relevance In this simulation modeling study, equalizing vaccination rates and reducing number of work contacts (which are relatively simple strategies to implement) reduced the both the inequity in hospitalizations and the absolute number of symptomatic infections and hospitalizations in all age and racial-ethnic groups.
Question What interventions might be effective in reducing the observed race-based disparities in seasonal influenza hospitalizations in the United States (US)?
Findings Simulated interventions equalizing comorbidity rates by in marginalized racial-ethnic groups would prevent up to 40% hospitalizations in these groups. Increasing vaccination rates among marginalized communities or reducing work contacts would reduce hospitalizations by over 10% and 20%, respectively, in most racial-ethnic groups and would similarly reduce symptomatic infections.
Meaning Increasing vaccination rates among marginalized communities, decreasing comorbidities in these groups, or reducing work contacts would help reduce the race-based disparities in seasonal influenza hospitalizations in the US.
1. Introduction
According to the Centers for Disease Control and Prevention (CDC), approximately 8% of Americans are infected with seasonal influenza each year [1]. The burden of seasonal influenza in the United States (US) varies widely yearly, with 360,000 hospitalizations and 21,000 deaths estimated for the 2022-23 season, which was typical [2]. Not all Americans, however, face the same risk of infection or severe outcomes. Young children (aged 0-4) and older adults (aged 65 and older) are more likely to experience severe outcomes than any other age group [3]. Significant racial and ethnic inequities also exist in seasonal influenza outcomes in the US, where the BIPOC (Black, Indigenous, and other People of Color) population faces a greater burden of exposure to infection and less access to prevention and treatment, resulting in higher age-adjusted rates of hospitalization and ICU admission [4, 5, 6, 3].
Disparities in disease outcomes between racial and ethnic groups are due to a variety of complex factors. Social determinants of health (SDOH) including socioeconomic status, house-hold structure, access to healthcare, education, and occupation [7] are known to affect influenza outcomes [8, 9]. Systemic inequities lead to higher rates of underlying medical conditions [10], which may worsen influenza outcomes. Furthermore, existing racial and ethnic disparities in vaccine uptake exacerbate differences in the rates of severe outcomes. These disparities might be due to a complex mix of factors, including mistrust in the healthcare system, lack of transportation or childcare, inability to take time off or to pay for the vaccine or the co-pay [11, 12, 13, 14]. BIPOC populations are, therefore, more frequently exposed to influenza, more susceptible to infection, more likely to develop symptomatic disease, and more likely to progress to hospitalization or death [15, 16].
In this paper, we simulated influenza transmission using a race-structured approach to investigate the racial and ethnic inequities in influenza outcomes, allowing us to understand and quantify the contribution of different social and structural factors in observed inequities in influenza hospitalizations. We tested five equity-promoting interventions aimed at reducing either structural inequities (represented by different levels of exposure based on race and ethnicity at the workplace) or health-related inequities (represented by different levels of comorbidities and/or susceptibility to infection). Projecting the effects of equity-promoting interventions into realistic synthetic populations may help inform policy decisions aimed at increasing equity while also increasing the social and economic benefits for influenza and other respiratory infectious diseases.
2. Methods
We constructed an agent-based model of influenza transmission adapted from COVASIM (COVID-19 Agent-based Simulator) [17, 18]. Agents interact with each other through a contact network in four locations: households, schools, work and community. Based on CDC’s definitions [19], we included five racial-ethnic groups in our model (referred below as race groups for succinctness): non-Hispanic white persons, non-Hispanic Black persons, non-Hispanic Asian persons, non-Hispanic American Indian or Alaska Native (AI/AN) persons, and Hispanic or Latino persons. We calibrated our network to US-based data to incorporate differences by race and ethnicity in each location. Households were stratified to match both size and race distributions according to the US Census data [20, 21]. We assumed three types of workplaces, low-, medium-, or high-risk workplace, according to size, and calibrated them to US-based workplace data [22, 23, 24]. Using data from [25, 26, 27] we modeled two types of schools: majority-white and racially diverse. Figure 1 shows the structure of race-based contacts in the model. Supplemental Figure S1 shows the fits to the data.
Disease progression probabilities were stratified by age and race in the model, according to [3]. We included vaccination in the simulation by assuming that vaccinated people will have a lower probability of symptomatic infection than those who are unvaccinated [28, 29]. The proportion of vaccinated people in each age and racial-ethnic group was based on historic vaccination rates from [15]. We assumed a vaccine efficacy against symptomatic infection of 50%, in agreement with estimates from typical influenza seasons [30]. We simulate a single influenza season with five million people and 100 runs. We report median disease outcomes with corresponding interquartile ranges. A complete description of the model is given in the Supplement.
2.1 Equity-promoting interventions under study
We tested five equity promoting interventions: 1.) “Equal Vaccination” scenario: the age-stratified vaccination rates are equalized for all race groups, choosing the highest vaccination rate as the reference (Fig. 2 A), 2.) “Equal Comorbidities” scenario: the age-stratified risks of disease progression given infection (hospitalization, ICU, and death) are equalized for all race groups (Fig. 2 B), 3.) “Equal Work-Risk” scenario: the race distribution of workers is proportional to that of the population regardless of workplace size (Fig. 2 C), 4.) “Low-Risk Workplaces” scenario: all workplaces are assumed to have less contacts (Fig. 2 D), and 5.) “Equal Disease Progression and Low-Risk Workplaces” a scenario where the probabilities of vaccinations and severe outcomes are equalized and workplaces are low-risk. Each of these scenarios targets different challenges faced by marginalized populations. The first two scenarios target underlying health inequities that affect the outcomes of an influenza infection, resulting in a higher probability of developing symptoms due for example, to reduced vaccination rates in the BIPOC populations, (targeted in scenario 1) or a higher rate of severe outcomes due, for example, to more prevalent comorbidities in marginalized populations (targeted in scenario 2). Scenarios 3 and 4 aim at addressing structural inequities that stem from the overall societal structure in which a substantially higher proportion of BIPOC populations are employed in lower-income frontline occupations. Scenario 5 is a combination of scenarios 1, 2, and 4.
Obviously, we do not pretend to solve the complex problem of systemic inequity in health in the US with these simple scenarios. Rather, the goal of this project is to characterize the potential impact of five equity-promoting interventions to the inequities in hospitalization rates for seasonal influenza.
3. Results
3.1 Results from Baseline Model
The model was adjusted to match the number of hospitalizations and the age- and race-stratified hospitalization rate ratios (HRRs) from a typical influenza season [3]. The hospitalization rate ratio was computed as the ratio between the number of hospitalizations in a given race- and age-group and the number of hospitalizations in the corresponding white age-group.
Our model closely matches the reported average hospitalizations and the HRRs for all age and race subgroups [3], Figures 3A and 3B. Because of the relatively small size of the AI/AN population, there was more variability observed for all ages in this race group.
Mild symptomatic influenza infections are not usually reported [31, 32], and marginalized race groups, who are likely to have less access to health care [10, 33, 34] and less ability to take time off when sick, are less likely to report their mild infections [31]. Hence, rates of symptomatic influenza infection are difficult to determine from data and are likely to be biased. Using our simulations, we estimate the number of symptomatic infections in each race- and age-group and corresponding symptomatic-infection rate ratios (SRRs), (calculated similarly to the HRRs, Figure 3C).
In our model, the Black, AI/AN, and Hispanic populations experience higher rates of symptomatic infection than those of the white population. This is especially true for the adult age-groups, for which the rates of symptomatic infection were between about 1.2 to 1.4 times higher (Fig. 3C). The Asian population, however, experiences rates of symptomatic infection closer to those of the white population for all age-groups (Fig. 3C).
3.2 Results from Equity-Promoting Interventions
Using the results in Figure 3B-C as a baseline, we can compare changes in the inequities in both symptomatic infections and hospitalizations that result from the equity-promoting interventions described in Section 2.1.
3.2.1 Effects on Symptomatic Infection
For the populations facing the most inequities, scenarios impacting work contacts improved in-equities the most relative to the baseline (Fig. 4A-B). Scenario 4 (low-risk workplaces) reduced the SRRs the most for all groups and races, with SSRs closest to equity (SRR = 1) for all race groups in most of the adult age groups (Fig. 4A, magenta boxplots). This intervention would also reduce symptomatic infections the most with a maximum reduction of 17%, 12%, and 18% in the 50-64 year-old Black, Hispanic, and AI/AN groups, respectively (Fig. 4B, fourth column). An intervention aiming at achieving a workplace distribution similar to the overall population distribution (denoted Equal Work-Risk) also reduced the inequities observed in the SRRs (Fig. 4A, orange boxplots), but it had little impact in the overall number of symptomatic infections averted. Achieving equity in vaccination rates (scenario 1, equal vaccination) resulted in significant gains in the number of symptomatic infections averted, with a maximum reduction of 17% infections averted in the Black young-adult age group (18-49 years old, Fig. 4B, first column), but had a more modest reduction of the SRRs (Fig. 4A, green boxplots). As expected, removing inequities in severe outcomes, did not improve inequities in symptomatic infection.
Finally, the combination scenario with both equal disease progression and low-risk work-places (scenario 5, equal disease progression and low-risk workplaces) reduced the SRRs the most across all age-groups for the three most marginalized populations (the Black, Hispanic and AI/AN adults, Fig. 4A, brown bars), and it averted the most symptomatic infections (maximum of 21% in the 50-64 year-old Hispanic population, Fig. 4B, fifth column). This is expected, as this scenario is a combination of the best single intervention scenarios.
Importantly, scenarios 1 (equal vaccination), 4 (low-risk workplaces), and 5 (equal disease progression and low-risk workplaces) resulted in significantly less symptomatic infections over-all, with less symptomatic infections in all ages and race groups. With these scenarios, there was a maximum of approximately 40% less symptomatic infections in marginalized populations (for scenario 5 in the Black population aged 18-49), and up to 26% and 23% less infections in the white (aged 18-49) and Asian (aged 50-64) populations, respectively (for scenario 5, equal disease progression and low-risk workplaces, in Fig. 4B and Table 2 in the Supplementary Materials).
3.2.2 Effects on Hospitalizations
For most age and race groups, these interventions improved the observed inequities when compared to the baseline (Fig. 5A). When considering single interventions alone, the intervention achieving the most equity varied by age and race. The equal comorbidities scenario performed the best for reducing both inequity and hospitalizations in children. Additionally, this scenario also performed the best for Black adults under 75 years old. For adults in the Hispanic population, however, the low-risk workplaces scenario outperformed the equal comorbidities scenario for adults 50 and older. These results suggest that the inequities in severe outcomes in the His-panic population are primarily due to inequities in occupational risk, while those in the Black population are primarily due to differences in comorbidity rates. Therefore, our results emphasize the need for different interventions to decrease inequity in different groups. Consistent with the results for the SRRs, the combined scenario 5 (equal disease progression and low-risk work-places scenario) achieved the most equitable results for most groups. This scenario was closest to equity for most age groups for the Black, AI/AN and Hispanic populations. Furthermore, this scenario resulted in the biggest decrease in hospitalizations for all races and all age groups, including the white population, with over 40% fewer influenza hospitalizations compared to baseline for all white adult groups (Table 3 in the SI and Fig. 5B, last column).
We then compared the effect of each scenario on the total number of hospitalizations to the baseline scenario (Fig. 5B and Table 3 in the SI). The largest gain for a single intervention was observed with the equal comorbidities scenario, with most marginalized groups seeing an over 40% reduction in hospitalizations. For example, if the Black adults aged 18-64 had similar disease progression as their white counterparts, they would experience over 70% fewer hospitalizations when compared to baseline. Under this scenario, the AI/AN and Hispanic groups saw over 50% and 37% reduction in influenza hospitalizations, respectively, in age groups younger than 65 years old. Interestingly, under this scenario, young children (0-4 years old) in the Asian group experienced 37.5% fewer influenza hospitalizations. Reducing contacts in the workplace (scenario 4) led to a modest decrease in hospitalizations, with the maximum reduction observed in the AI/AN group (50% reduction in the 50-65 age group). Equalizing the vaccination rates across ages and race groups led to a small decrease in total hospitalizations for all groups but those aged 5-17 and those over 75 years old in the AI/AN group (33% and 50% less hospitalizations, respectively). As before, combining two interventions resulted in important reductions in the total number of hospitalizations in all groups (maximum reduction: 75% in young children and adults aged 50-64 in the AI/AN group). As with the symptomatic infections, it is key to note that equity promoting interventions did not result in less hospitalizations for those in the BIPOC communities exclusively. All the interventions but one (equal work risk) resulted in fewer influenza hospitalizations in all groups. For example, equalizing vaccination rates resulted in at least 10% less infections in all white age-groups.
4. Discussion
As the recent COVID-19 pandemic has highlighted, the burden of infectious diseases disproportionately falls on marginalized populations [35, 5, 36, 37]. There is currently new momentum to close these disparity gaps with interventions aimed to increase prevention, vaccination, access to healthcare, and treatment [38, 39, 40, 41, 42]. To better support those in marginalized communities facing current outbreaks and to prepare for future pandemics, it is important to quantify the potential impact of interventions aimed to reduce inequities [43]. In this project, we constructed an agent-based mathematical model of influenza transmission that captures both structural differences between racial and ethnic groups in the US (household composition, work-related risk, and school contacts) and differences in health outcomes (rates of disease progression) deriving from long-held inequities and systemic racism. Our goal is to provide a framework for studying racial inequities in disease outcomes and to evaluate the effects of potential interventions.
We explored the effects of five equity-promoting interventions aimed to reduce different social determinants of health that affect influenza outcomes. These determinants can be structural inequities leading to differences in exposure (represented by work contacts and type of occupation), health inequities (represented by differences in disease progression rates), or both. Strategies for reducing inequities in disease outcomes could include vaccination campaigns that are tailored to reach marginalized groups, interventions intended to reduce comorbidities in these groups, or more systemic approaches aimed at reducing inequities in health outcomes (e.g. improve health care coverage and access in marginalized populations). Strategies for reducing inequities in work contacts could include social-distancing measures in workplaces that reduce the probability of transmission and exposure, incorporating better air filtration systems or more systemic cultural and workplace changes. Evidently, systemic changes are more difficult to attain and will require major commitment from all sectors of society [44]. Nevertheless, we considered it important to consider them and to quantify the potential impact such strategies could have on influenza hospitalizations.
Our results suggest that tailored campaigns to increase vaccination rates among marginalized communities may have an equalizing effect in symptomatic infection rates, with younger adults (aged 18-49 years) in Black, AI/AN and Hispanic communities benefiting the most from this intervention, with up to 17% less symptomatic infections compared to the baseline scenario. This is important because these populations have generally less access to medical care [8, 10] and work benefits (e.g. working from home [22] or paid sick leave [45]) and are more likely to be employed in frontline occupations [22, 23]. Interestingly, their corresponding white counterparts would also experience fewer symptomatic infections, with up to 15% less compared to baseline. In addition, this strategy would lessen the significant economic burden [46] posed by symptomatic influenza every year. Moreover, our results suggest that policies that reduce race-based comorbidities would have a large impact on the inequities in influenza hospitalizations. In addition, the scenario implementing fewer contacts in the workplace (called low-risk workplaces) had a much larger effect, and benefited more groups, than the one where contacts were maintained at the same rates, but the distribution of contacts was proportional to that of the population. This suggests that implementing (relatively) small changes in workplaces (such as better air filtration systems, social distancing measures) can have a large impact in influenza outcomes. Finally, it is important to note that our results suggest that improving equity and fairness in health outcomes and reducing contacts in workplaces resulted in significant reductions in the number of symptomatic influenza infections and hospitalizations for all racial-groups, not only those in marginalized groups targeted by these interventions.
This study, like any mathematical modeling analysis, has limitations. Our model is a simplification of a very complex problem: health inequities are the result of a multi-factorial problem rooted in a long history of systemic racism, socioeconomic deprivation, and a fragmented healthcare system that makes it difficult to align data, resources, and interventions, optimally [47]. As such, we included in this model SDOH for which data was readily available while maintaining the tractability of the model. We considered race-stratified household size, school composition, workplace and disease progression, but additional SDOH could be included (e.g. specific influenza comorbidities such as heart disease, diabetes, and asthma [48, 49]). We did not include interracial households, which account for 10% of households in the United States, which could result in epidemics with more assortative mixing. However, agents in our model contact people from other racial-ethnic groups in other locations, minimizing this effect.
We use racial identity as a proxy for exposure to racism and to study inequities in influenza cases and hospitalizations, but racial identity is known to be a poor proxy for studying health inequities [50]. However, other SDOH (e.g. occupation, access to health care, socioeconomic status) are usually not collected along with case information. This highlights the pressing need to collect more informative social determinants of health data, so that inequities can be better understood and analyzed.
Researchers and decision makers are increasingly recognizing the need to incorporate fairness considerations into analyses aiming to evaluate public health interventions [51, 52, 53, 54, 55], including analysis using mathematical models of infectious diseases. Previous work has focused on reducing racial and ethnic inequities in COVID-19 outcomes [56, 57, 58] assessing the impact of COVID-19 mitigation strategies on different racial and ethnic groups [59, 60], optimizing equitable vaccine distribution for COVID-19 and influenza [61, 56, 62] studying inequity between subgroups for influenza and other diseases [63, 64], and promoting equity in resource allocation between geographic regions for Ebola [65]. We hope that the present work will be a helpful addition to the existing literature and will promote further discussion and use of quantitative methods to evaluate equity-promoting interventions in public health.
Funding Statement
This work was partially supported by grants from the National Institutes of Health (UM1AI068635). L.M. and D.D. were also supported by a grant from Centers for Disease Control and Prevention (NU38OT000297-02) through their cooperative agreement with the Council of State and Territorial Epidemiologists. L.M. was partially funded through the IDCRC Early Career Investigator Pilot Award subaward #A925844 (NIH # 5UM1AI148684-05). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Data Availability
All data produced will be available online.
https://github.com/laura-matrajt/inequities_flu_hospitalizations
Supplementary Material
Supplementary Methods
In our simulation, individuals in the population are represented in our model by agents who interact with each other through a contact network, similarly to [1, 2]. Influenza transmission occurs when agents in the network contact each other. At the beginning of a simulation (representing the beginning of an influenza season), most agents in the network are susceptible. Upon exposure, susceptible individuals can become infected, but are not immediately infectious. Once infectious, individuals are considered either asymptomatic (they can infect others but they have no symptoms) or presymptomatic. Presymptomatic individuals develop mild symptoms and can either recover or progress to severe or critical disease. In our model, children are assumed to be more likely to develop symptoms [3, 4]. Severe cases are those individuals who require hospitalization, and critical cases are those who are admitted to the ICU. Critical cases can either recover or die. In our model, children are assumed to be more likely to develop symptoms [3, 4]. The probability that an individual becomes symptomatic, and has mild, severe, or critical symptoms depends on both age and race/ethnicity. Transition times between these states are sampled from a log-normal distribution with mean times given in Table 1. A full list of disease parameters used in the model is also given in Table 1.
Our contact network was adapted from a previously developed network model, SynthPops [5], to include both age-specific and race-and-ethnicity specific contact patterns in household, school, workplace, and community settings. SynthPops constructs synthetic networks that are calibrated to data such as age distributions, household sizes, employment rates, workplace sizes, school enrollment rates, school sizes, etc. In addition, we calibrated our network to US-based data to incorporate differences by race and ethnicity as described in the following sections.
A more detailed description of the base models, COVASIM and SynthPops, is given by Kerr et al. [1] and Mistry et al. [5].
Incorporating racial and ethnic differences in the model
Using the CDC race and ethnicity definitions [6], we included five categories defining racial and ethnic groups (referred below as race groups) in our model: non-Hispanic white persons, non-Hispanic Black persons, non-Hispanic Asian persons, non-Hispanic American Indian or Alaska Native (AI/AN) persons, and Hispanic or Latino persons.
Household stratification by race group
Individuals were assigned to a given race group during household creation as follows: After generating households based on given household-size and age distributions mimicking the US household distribution [7, 8], we assigned a race group to each household using a probability distribution based on US Census data on household size by race and ethnicity [7]. It is important to note that this method does not capture interracial households. The age and race distributions of households calibrated in the model are in agreement with Census data (Figure S1(A-B)).
Work stratification by race group
We used US-based workplace data [9, 10] to stratify workplaces and work distribution by race group. To do this, we first create each workplace with a size determined according to a distribution of workplaces sizes from US Census data [11]. We used the workplace size as a proxy to determine the level of infection risk at the workplace, divided into three categories: low-, medium-, or high-risk environment. We considered low-risk occupations to be non-frontline work and assumed that those in low-risk occupations have the fewest work contacts. Medium-risk and high-risk occupations were considered to be higher-income and lower-income frontline work, respectively. Those in medium-risk occupations have more work contacts than those in low-risk occupations, but fewer than those in high-risk occupations. Workplaces with at most 50 workers are deemed to be low risk, those with between 50 and 100 workers to be medium risk, and those with more than 100 workers to be high risk. We assigned working individuals (determined from age-stratified employment rates) to workplaces based on racial and ethnic distributions for both the total workforce [9] and low-income, high-income, and non-frontline work [10]. Figure S1(E) shows the fit of our model workplace distribution to the data.
School stratification by race group
We assign students to schools, which include pre-schools through university, based on school racial-composition data [12, 13, 14]. Because schools composition varies greatly by region, we consider two types of schools: majority-white schools and schools that are racially diverse. We adapt examples of these two types of schools from [12] to reflect the general US population using US Census data [13, 14] and use them as probability distributions for assigning students to schools. To create a new school, we first draw an age and race from probability distributions gathered from Census data. The age determines the school type (pre-school, elementary school, middle school, high school, or university), and the race determines the racial distribution of the school (majority-white if white, and racially diverse, otherwise). We use school size distributions to determine the size of the school, and we assign students based on the race-stratified probability distributions. Figure S1(C-D) shows the average school compositions for our simulated population compared to the example data.
Influenza outcomes stratification by race group
To capture differences between racial groups in influenza outcomes, we informed transition probabilities from developing symptomatic infection through hospitalization, ICU admission to recovery or death from available clinical data [15]. Since we do not have age and race stratified data on the rates of symptomatic infection, we run the model with the parameters given in Table 1 to estimate these rates. Influenza vaccination has been shown to reduce the risk of symptomatic infection [16, 17], but influenza vaccination rates vary by race and age, with decreased vaccination coverage in marginalized populations [18, 19]. While our model does not incorporate vaccination explicitly, we included vaccination by assuming that people who got vaccinated in the current influenza season will have a lower probability of symptomatic infection than those who are unvaccinated. We determined the proportion of vaccinated people in each age and racial-ethnic group based on race-stratified historic vaccination rates from [19]. We assumed a vaccine efficacy against symptomatic infection of 50%, in agreement with estimates from typical influenza seasons [20].
Supplementary Figures and Tables
Footnotes
Extended the abstract, moved methods section to the supplemental material