Abstract
Background SARS-CoV-2 vaccines have been shown to be safe and effective against infection and severe COVID-19 disease worldwide. Certain co-morbid conditions cause immune dysfunction and may reduce immune response to vaccination. In contrast, those with co-morbidities may practice infection prevention strategies. Thus, the real-world clinical impact of co-morbidities on SARS-CoV-2 infection in the recent post-vaccination period is not well established. We performed this study to understand the epidemiology of Omicron breakthrough infection and evaluate associations with number of comorbidities in a vaccinated and boosted population.
Methods and Findings We performed a retrospective clinical cohort study utilizing the Northwestern Medicine Enterprise Data Warehouse. Our study population was identified as fully vaccinated adults with at least one booster. The primary risk factor of interest was the number of co-morbidities. Our primary outcome was incidence and time to first positive SARS-CoV-2 molecular test in the Omicron predominant era. We performed multivariable analyses stratified by calendar time using Cox modeling to determine hazard of SARS-CoV-2. In total, 133,191 patients were analyzed. Having 3+ comorbidities was associated with increased hazard for breakthrough (HR=1.2 CI 1.2-1.6). During the second half of the study, having 2 comorbidities (HR= 1.1 95% CI 1.02-1.2) and having 3+ comorbidities (HR 1.7, 95% CI 1.5-1.9) were associated with increased hazard for Omicron breakthrough. Older age was associated with decreased hazard in the first 6 months of follow-up. Interaction terms for calendar time indicated significant changes in hazard for many factors between the first and second halves of the follow-up period.
Conclusions Omicron breakthrough is common with significantly higher risk for our most vulnerable patients with multiple co-morbidities. Age related behavioral factors play an important role in breakthrough infection with the highest incidence among young adults. Our findings reflect real-world differences in immunity and exposure risk behaviors for populations vulnerable to COVID-19.
Introduction
Vaccines against SARS-COV-2, have been developed and shown in numerous studies to be safe and highly effective at reducing SARS-CoV-2 infection and COVID-19 disease [1–3]. However, clinical trials and population-based observational studies excluded or did not compare certain groups at highest risk for severe outcomes of COVID-19. The impact of age and burden of immune disorders or certain chronic conditions on vaccine effectiveness (including boosting) in terms of acquisition of infection, has not been adequately studied in our current era dominated by omicron SARS-CoV-2 subvariants.
Although there have been impressive advances in our understanding of protective immunity against COVID-19 after vaccination and natural infection, risk of infection is not completely understood for our most vulnerable patients including individuals who are immune compromised (age, HIV, malignancies, solid organ transplant, stem cell transplant) or have chronic illnesses (diabetes, obesity, chronic liver disease, and chronic kidney disease). SARS-CoV-2 vaccines stimulate both B and T cell responses to virus spike protein to elicit an effective immune response [4]. Those with dysfunctional immunity have been observed to have lower responses to vaccination with antibody titers as indicators of immunogenicity [5]. Questions remain about how lower immunogenicity translates into diminished clinical effectiveness of COVID-19 vaccines in real world populations with different co-morbidities. Individuals who have chronic disease, advanced age, and/or immunodeficiencies may be at higher risk for breakthrough infection due to poor vaccine response (defined as COVID-19 infection after completion of all require doses with a typical 2-week lag period) [6]; however, they may also practice better infection prevention such as mask wearing, avoiding travel or large gatherings, and social distancing [7].
The incidence of SARS-CoV-2 breakthrough infection is an increasingly important issue worldwide, with vulnerable populations at high risk of infection at a time when vaccine-induced immunity may not be fully optimized. To gain insight into this issue and to inform public health decision making, our study aimed and was able to determine the incidence and risk factors associated with SARS-CoV-2 breakthrough infection in the Omicron-variant era among vaccinated and boosted individuals.
Methods
Study Design
We performed a clinical cohort study of breakthrough infection in the first year of SARS-CoV-2 Omicron-variant era (January 1, 2022, until December 31, 2022) among fully vaccinated and boosted adults (18 years and older) as per CDC/FDA vaccine guidelines for COVID-19 [8], All demographic, lab, vaccine, and comorbidity data were collected from the Northwestern Medicine (NM) Enterprise Data Warehouse (EDW)[9]. Each participant had a unique study identifier and protected health information was stored separately with access limited to the principal investigator.
Population & Definitions
We included adults that were boosted prior to the Omicron era, which we defined as those who received at least 3 (first dose mRNA) or 2 (first dose J&J) SARS-CoV-2 vaccine doses prior to December 15, 2021. To include a representative population that was likely to have SARS-CoV-2 testing performed within the NM system, we only included patients if they had at least two medical system visits (including inpatient, outpatient, telemedicine, and lab testing) at least 180 days apart between January 1, 2020, and December 15, 2021.
We observed individuals in this cohort from January 1, 2022, to December 31, 2022, for incident breakthrough infection with SARS-CoV-2 as our primary outcome defined as the first positive SARS-CoV-2 PCR test performed at an NM facility after January 1, 2022. Patients who tested positive for SARS-CoV-2 after their most recent booster dose but prior to January 1, 2022, were excluded. Data was accessed throughout this study period. Comorbidities of interest (diabetes, obesity, solid organ transplant, stem cell transplant, HIV, hematologic malignancy, chronic liver disease, and chronic kidney disease) were identified using ICD9/ICD10 coding.
Statistical Analysis
Descriptive statistics, including median (IQR) and counts (%), were calculated for patient characteristics and compared between those with and without breakthrough SARS-CoV-2 infection during the study period. We calculated cumulative incidence as the number of breakthrough infections divided by the total number of those at risk (no prior breakthrough infection during study and not right censored) and plotted these curves with 95% confidence intervals assuming normality over the study period. Patients were right censored in cases of death, additional SARS-CoV-2 vaccines, or loss to follow-up. Loss to follow-up was defined as 90 days without a visit at NM (outpatient visit, hospital admission, or laboratory testing). We performed Cox regression modeling to determine proportional hazards of breakthrough infection with the following covariates: age, sex, race, ethnicity, time from booster dose to study period start, and number of comorbidities. We included the overall number of comorbidities rather than the specific comorbidities themselves as this study aimed to measure associations with overall health, not individual disease epidemiology. The proportional hazards assumption was assessed graphically with cumulative incidence curves for each covariate. For covariates that violated this assumption, interaction terms for calendar time were considered. Multicollinearity was assessed with generalized variance inflation factors (GVIFs). If (GVIF1/(2⋅DF))2 > 3, removal of those variables from the model was considered. Linearity of log hazards for continuous covariates was examined with martingale residuals from the fitted model. The car, survival and tidycmprsk packages in R 4.2.3 software were used for statistical analysis and plot production[10–13].
Results
Clinical and demographic characteristics of the cohort included and analyzed are presented in Table 1. In total, there were 133,191 patients in the cohort with a median (IQR) age of 61 years (47, 72) 63% female sex, 84% white, 95 % Non-Hispanic or Latino ethnicity, and 77% with any comorbid condition.
Overall, 99% of individuals received mRNA SARS-CoV-2 booster. One booster dose vaccine was administered to 99.2% of individuals and 0.8% had received more than one booster dose.
Of the total cohort, 7,904 (5.9%) tested positive for SARS-Cov-2 infection by PCR during the study period in a median (IQR) of 135 (34, 196) days. The infections occurred in a median (IQR) of 196 (115, 260) days after a booster dose of vaccine. Approximately 31% of the total cohort was right-censored due to loss to follow-up, 21% due to additional boosters, and 0.4% due to death.
The proportional hazards assumption in Cox regression was violated for the following covariates: age, race, days since vaccination, and number of comorbidities. Graphical assessment indicated that significant changes in incidence rate trends began near the 6-month (halfway) mark of the study period. To account for this, we introduced an interaction term to calculate separate hazard ratios for periods 1 (Jan 1, 2022, to June 30, 2022) and 2 (July 1, 2022, to Dec 31, 2022) for those 4 covariates.
In Cox regression multivariable analysis (Table 2), increasing age was associated with reduced hazard during this first period when controlling for other covariates. The greatest reduction in hazard compared to the reference age group (18-39) was seen in the 80+ age group (HR 0.41 95% CI 0.36-0.46). Similar but lower reductions were observed in the 60-79 (HR 0.44 95% CI 0.4-0.48) and 40-59 (HR 0.71 95% CI 0.66-0.76) age groups. A significant interaction between calendar time and the hazards for age was found in each age group compared to the reference group (p < .001 for each). In the second half of the study period, only the 60-79 age group demonstrated a significantly different hazard than the reference group (HR 0.86, 95% CI 0.76-0.99). Trends in cumulative incidence stratified by number of comorbidities in are observed in Fig. 1.
Cox regression modeling (Table 2) also revealed an increased hazard for those of Asian race compared to the reference group (white race) in both periods (period 1: HR 1.32 95% CI 1.19-1.46; period 2 HR 1.23 95% CI 1.03-1.48). Significant associations with hazard were seen in the Black or African American race group compared to the reference group, though the directions of these associations differed between study periods (period 1: HR 0.83 95% CI 0.73-0.94; period 2: HR 1.16 95% CI 1.03-1.48). People of Hispanic, Latino, or Spanish ethnicity also had a higher hazard for COVID-19 breakthrough infection compared to those not in that ethnic group (HR 1.48 95% CI 1.35-1.62). Days from booster to start of study showed an increased hazard of 1.03 per week (95% CI 1.03-1.04) and later 1.01 per week (95% CI 1.00-1.02). The interaction term for calendar time and this variable was significant (p < .001).
Compared to having none, having one or two comorbidities was not significantly associated with increased hazard for Omicron infection during period 1 when controlling for other variables. In contrast, having three or more comorbidities was associated with increased hazard for breakthrough (HR=1.16 CI 1.08-1.26). During the second half of the study, having one comorbidity (HR = 1.27 95% CI 1.11-1.46), having two comorbidities (HR= 1.45 95% CI 1.26-1.67) and having three or more comorbidities (HR 1.73, 95% CI 1.51-1.97) were all associated with increased hazard for Omicron breakthrough compared to those with no comorbidities. Cumulative incidence over time between co-morbidity groups and stratified by age category are presented in Fig. 2. Cumulative incidence over time between age categories and stratified by co-morbidities groups are presented in Fig. 1.
Discussion
In this large cohort study across an urban healthcare system in Chicago, we found that omicron breakthrough infections were common with a one-year cumulative incidence of 5.1% for all vaccinated and boosted individuals. We observed a higher incidence of breakthrough infections among individuals with multiple comorbidities and those who had a longer time since vaccine booster. However, age was also an important factor with increasing age independently associated with reduced risk of Omicron breakthrough likely due to less risk taking and lower exposure among the highest age groups. We found racial and Ethnic differences in infection incidence as demonstrated by increased risk of Omicron breakthrough for individuals identifying as Asian race or Hispanic/Latino ethnicity, again likely due to differences in exposure from social behaviors rather than inherent immunogenetic differences. Perhaps our most intriguing finding was that across all age groups, those with three or more co-morbidities had the highest incidence of SARS-CoV-2 infection, and this factor did not appear to alter individual behavior to lower exposure as much as age. The immune effects of co-morbidity might have overwhelmed any decrease in social behavior risk related exposure. Thus, acquisition of Omicron SARS-CoV-2 was complex with behavioral and individual risk or immunologic factors driving SARS-CoV-2 infection despite vaccination.
Our findings are consistent with previous studies that have reported a significant burden of breakthrough infections despite vaccination efforts. For example, a study by Stouten et al. [14] found a similar cumulative incidence of 11.2% among vaccinated individuals in a different urban setting. Similarly, Sun et al.[15] reported an incidence of 7.1% in those without immune dysfunction and slightly higher rates for those with specific diseases affecting the immune system. Our cumulative incidence is lower than the ones reported in these studies. This likely relates to several factors including study population, methodologic differences, and SAR-CoV-2 variant circulating during the study period. The predominant circulating variant of concern in prior studies was the Delta variant, which has been the most studied in terms of breakthrough and consistently demonstrated high incidence rates [16]. Additionally, studies conducted in rural populations have demonstrated lower incidence of breakthrough [17, 18]. One study in a small cohort in New York City reported an incidence <1% [19], but this discrepancy is likely due to study time period and significantly smaller sample size compared to our study and others with higher incidence. As expected, and generally observed, evidence suggests that breakthrough infections are common across many different populations and healthcare settings and may vary with overtime with shifts in vaccine coverage, population immunity, and circulating SARS-CoV-2 variant strain.
Our study revealed an intriguing finding in that having three or more comorbid diseases heightened the risk of breakthrough infections across all age groups. Having more comorbidities increases risk of infection greater than one would expect of the individual comorbidities. Surprisingly, this heightened risk did not seem to prompt significant alterations in individual behavior. It is plausible that the impact of these comorbidities on immunity overshadowed any potential behavioral changes. Alternatively, it could be that comorbidities are not as influential in driving behavioral and lifestyle changes as age. However, the latter would contrast with prior studies suggesting that perceived vulnerability and severity, along with self-efficacy and intention, are drivers of COVID-19 protective behaviors[20–23]. Our findings align with those of Smits et al. [24], who examined four specific comorbidities and observed that patients with two or more of these conditions also faced a greater risk than would be expected based on the individual effects of each comorbidity. This contrasts with the findings of Walmsley et al. [25], who did not identify significant variations in the rate of breakthrough infections among individuals with underlying comorbidities. Theses investigators acknowledged that this discrepancy may be due to the low number of participants with comorbid diseases in their study, which could have hindered the identification of a clear association.
Immune compromising co-morbidities have been shown to increase risk of SARS-CoV-2 infection [26];[27–29], however, we did not find any significant increased risk among those in our immunocompromised groups of stem cell or solid organ transplant recipients as has been previously reported[28]. Whether due to the impact on the immune system or lack of an effect on behaviors, our findings reveal that independent of age or type of comorbidity, an increasing amount of comorbidity burden increased the risk for Omicron SARS-CoV-2 infection after vaccination within a population most vulnerable to worse outcomes and severe COVID-19 disease. Thus, public health messing should continue to emphasize the importance of infection prevention measures within this key group.
We noted that older individuals became at higher risk for infection at later time points, suggesting reduction in vaccine response or increase in risk taking behaviors over time. This finding is consistent with a previous study that demonstrated an inverse association between age and antibody titers only three months after mRNA vaccination [30]. Similarly, a study focusing on immunocompromised individuals found that although vaccination was associated with a modest risk reduction, older individuals continued to have higher rates of breakthrough SARS-CoV-2 infections [15]. Research has shown that as the immune system ages, the number of naïve T and B cells decreases, which can lead to reduced vaccine efficiency and a predisposition to breakthrough infections. While the quality of antibodies remains unaffected, aging is associated with a decrease in the quantity of antibodies produced after vaccination, rendering individuals more susceptible to infection [31, 32]. Our results support prior evidence establishing that, independent of comorbidities or immunocompromised status, increased age is a strong risk factor for breakthrough infection likely due to declining immunity over time. Overall, this evidence suggest that the diminished quantity of antibodies produced in older individuals puts them at greater risk compared to the longer-lasting protection seen in younger people. Our findings also support recent recommendations for SARS-CoV-2 vaccine boosting approximately twice per year among individuals over 65 years of age [33].
In contrast to our findings regarding older individuals, multiple studies have identified younger age as a risk factor for COVID-19 breakthrough infection[14, 25, 34] . It’s important to note that these studies varied in design, cohort size, and follow-up duration, which may have influenced their ability to capture current evolving trends. These findings align with our observations shortly after vaccination, where younger age emerged as an important risk factor. We attribute this phenomenon to social behaviors playing a larger role than age during earlier time periods, when immunity has not yet significantly declined and vaccine response more protective. Young adults engage in more social interactions and risk-taking behaviors for respiratory virus acquisition compared to the elderly, potentially increasing their exposure to SARS-CoV-2through daily activities [35, 36]. Certain occupations, such as healthcare work, transit operation, and retail roles, have been identified as particularly high-risk for breakthrough infection despite full vaccination [37, 38]. Additionally, younger individuals may perceive themselves as less vulnerable to severe COVID-19 disease and therefore be more likely to disregard public health precautions such as social distancing and mask-wearing, or they may be less adept at recognizing the signs and symptoms of the virus, leading to further spread of the disease [39, 40]. Thus, it appears that prior to a significant decline in immunity and vaccine efficacy, behaviors among younger individuals likely drive and risk for breakthrough SARS-Co-V2 infections.
Time since last booster dose and certain racial/ethnic group also emerged as independent risk factors. A study described a similar finding in relation to the time since second vaccine dose [41]. This finding is consistent with waning vaccine immunity that has been previously described [42–46]. Our study also highlights the increased risk of Omicron SARS-CoV-2 infection among vaccinated individuals from Hispanic, Asian, and Black ethnic/racial groups. We believe it is unlikely that there is a biological basis to these associations– rather due to well-established socioeconomic factors known to play significant roles in influencing the occurrence and outcomes of SARS-CoV-2 and COVID disease [47–49]. Prior studies examining breakthrough SARS-CoV-2 infection do not comment on demographic information, likely limited geographical regions being studied. [14, 16]. However, reports in immunocompromised populations did not find any significant racial or ethnic differences [15, 50]. Disparities in healthcare access and quality may contribute to certain demographic groups facing a higher risk of SARS-CoV-2 infection[47, 51] .
Engagement in higher-risk occupations, larger household size, lower income level, distrust in healthcare, lack of health insurance, and unequal access to healthcare services are among the key contributors to these disparities [52]. Efforts to address these disparities and improve healthcare outcomes for high-risk demographic groups are needed and should include interventions aimed at expanding access to healthcare and insurance, establishing more equitable care models, and addressing the underlying social determinants of health.
Our study has several limitations that should be considered when interpreting the results. Firstly, our research was conducted using electronic medical records within a single urban healthcare system, which likely did not capture all breakthrough SARS-CoV-2 infections. Patients who sought care outside of our NM system or opted for at-home testing were not always identified in our analysis, potentially leading to an underestimation of the true cumulative incidence of breakthrough SARS-CoV-2 infections. Additionally, mild, or asymptomatic breakthrough infections may be underreported in our study, as individuals with less severe symptoms might be less likely to seek medical care or testing. Secondly, given our study population selection criteria (see above), we may have included individuals who were generally more proactive about seeking healthcare services and thus may be overrepresented. This may affect the generalizability of our findings to a broader less engaged population. Additionally, the findings of this study may not be generalizable to other populations or healthcare settings. Factors such as population demographics, vaccination rates, and healthcare infrastructure can vary widely between different regions and may impact the incidence and risk factors for breakthrough SARS-CoV-2 infections. Lastly, our study was conducted during a specific timeframe and primarily focused on the Omicron variant. The future incidence and risk factors for breakthrough SARS-Co-2 infections will likely vary with evolving population immunity (natural and/or updated vaccines) and changes in circulating immune evasive Omicron sub-variants. Despite these limitations, we studied a large population with nearly 8,000 breakthrough infection events during the Omicron era and were able to assess associations with key demographic and clinical characteristics driving SARS-CoV-2 acquisition using strong statistical methods.
Conclusion
Our study sheds light on the complex interplay of factors influencing Omicron breakthrough infections in a large urban healthcare system in Chicago among a population of SARS-CoV-2 vaccinated and boosted individuals. We found a substantial one-year cumulative incidence of 5.1%, highlighting the ongoing challenges posed by SARS-CoV-2 despite vaccination efforts. These findings reflect real-world differences in immunity and exposure risk behaviors for populations vulnerable to Omicron variants of SARS-CoV-2 worldwide. By identifying key risk factors and disparities, our findings can inform targeted public health interventions to mitigate the impact of breakthrough infections in vaccinated populations. Public health messages should continue to emphasize the importance of considering both co-morbidities and age as critical factors in understanding and mitigating the risk of SARS-CoV-2 acquisition, ensuring that interventions are tailored to address the specific needs of vulnerable populations. As people make decisions about booster vaccinations, our findings provide more information for patients to consider personal risk in their decision making. Ongoing public health surveillance and research are crucial to understand the long-term effectiveness of vaccines against SARS-CoV-2. Future research should focus on understanding mechanisms of declining immunity, immune evasion by SARS-CoV-2 viruses, drivers of acquisition behavior, and optimizing protective vaccine.
Financial Disclosures
Research reported in this publication was supported, in part, by the National Institutes of Health’s National Center for Advancing Translational Sciences (grant number UL1TR001422) which supports the Northwestern Medicine Enterprise Data Warehouse. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. Research reported in this publication was also supported by the Division of Infectious Diseases emerging and re-emerging pathogens program (EREPP) at the Northwestern University Feinberg School of Medicine
Supporting Information
S1 Dataset. Vaccine Cohort Dataset.
S2 Supporting Information. Vaccine Cohort Dataset Codebook
Data Availability
All relevant data are within the manuscript and its Supporting Information files.