Abstract
Host immune responses are a key source of selective pressure driving pathogen evolution. Emergence of many SARS-CoV-2 lineages has been associated with improvements in their ability to evade population immunity resulting from both vaccination and infection. Here we show diverging trends of escape from vaccine-derived and infection-derived immunity for the emerging XBB/XBB.1.5 Omicron lineage. Among 31,739 patients tested in ambulatory settings in Southern California from December, 2022 to February, 2023, adjusted odds of prior receipt of 2, 3, 4, and ≥5 COVID-19 vaccine doses were 10% (95% confidence interval: 1-18%), 11% (3-19%), 13% (3-21%), and 25% (15-34%) lower, respectively, among cases infected with XBB/XBB.1.5 than among cases infected with other co-circulating lineages. Similarly, prior vaccination was associated with greater protection against progression to hospitalization among cases with XBB/XBB.1.5 than among non-XBB/XBB.1.5 cases (70% [30-87%] and 48% [7-71%], respectively, for recipients of ≥4 doses). In contrast, cases infected with XBB/XBB.1.5 had 17% (11-24%) and 40% (19-65%) higher adjusted odds of having experienced 1 and ≥2 prior documented infections, respectively, including with pre-Omicron variants. As immunity acquired from SARS-CoV-2 infection becomes increasingly widespread, fitness costs associated with enhanced vaccine sensitivity in XBB/XBB.1.5 may be offset by increased ability to evade infection-derived host responses.
Host immune responses are a key source of selective pressure influencing the evolutionary dynamics of pathogens.1–3 In the context of expanding population immunity, successive SARS-CoV-2 lineages have shown increasing capacity to evade both vaccine-derived and infection-derived immune responses throughout the course of the COVID-19 pandemic.4,5 Whereas incremental reductions in protection against Epsilon, Gamma, Delta, and other early variants were generally found to be modest,6–9 the Omicron BA.1 lineage was associated with ≥1.8-fold higher risk of breaking through infection-derived immunity in comparison to Delta,10–13 as well as markedly lower effectiveness and duration of COVID-19 vaccine-derived protection.14–16 Although epidemiologic studies did not find strong evidence of differential vaccine-derived or infection-derived protection against the BA.2 lineage in comparison to BA.1, the subsequent BA.4 and BA.5 lineages were associated with notably increased risk of re-infection relative to earlier Omicron lineages,17–20 as well as reduced vaccine effectiveness.21,22 Monitoring the ability of emerging lineages to evade immunity from vaccination or prior infection is central to ongoing efforts aimed at mitigating the burden of SARS-CoV-2, similar to experience with vaccines against influenza,23 pneumococcus,24 and other infectious disease agents.25,26
The XBB/XBB.1.5 Omicron lineages emerged via recombination of BA.2.10.1 and BA.2.75 sublineages and overtook BQ.1/BQ.1.1, along with other BA.5-related lineages, as the leading cause of new infections within the US by late January, 2023.27 While XBB/XBB.1.5 evades neutralization by infection-derived antibodies,28,29 early observational studies have reported that updated (bivalent) COVID-19 booster vaccination confers substantial protection against symptomatic XBB/XBB.1.5 infection.30 It remains unclear whether XBB/XBB.1.5 differs from BA.5-related lineages in its sensitivity to host responses acquired through prior vaccination or infection.
Study population and setting
To address this question, we analyzed data from 31,739 individuals who tested positive for SARS-CoV-2 infection in outpatient settings during the period from 1 December, 2022 to 23 February, 2023 within the Kaiser Permanente Southern California healthcare system, an integrated care organization providing comprehensive medical services to ∼4.7 million residents of southern California across outpatient, emergency department, inpatient, and virtual settings. We limited our sample to individuals whose specimens were processed using the ThermoFisher TaqPath COVID-19 combo kit (TF) assay (39.2% [31,739/80,894] of all eligible cases diagnosed during the study period; see Methods) in order to compare characteristics and clinical outcomes among cases whose positive tests included detection of the spike (S) gene probe or S-gene target failure (SGTF)—a well-described proxy for distinguishing XBB/XBB.1.5 and other BA.2-origin Omicron lineages (associated with S-gene detection) from other co-circulating lineages descending from BA.4 and BA.5 (associated with SGTF).31 Within a random sample of cases for whom sequencing results were available, the positive and negative predictive values of S-gene detection for ascertainment of XBB/XBB.1.5 infection and non-XBB/XBB.1.5 infection, respectively, were 98.2% (269/274) and 99.7% (1,592/1,597). Among sequenced specimens exhibiting SGTF, 99.2% (1,585/1,597) and 0.4% (7) were descendants of BA.5 and BA.4 lineages, respectively.
New detections of SARS-CoV-2 declined over the study period amid reductions in outpatient testing volume (Figure 1). This decline in outpatient detections exceeded reductions in SARS-CoV-2 detection in inpatient settings, where all newly-admitted patients continued to be screened for SARS-CoV-2 infection at the point of admission throughout the study period. The proportion of TF-tested outpatient cases inferred to be infected with the XBB/XBB.1.5 lineage based on S-gene detection increased from 21.1% (45/213) as of 1 December, 2022 to 77.8% (49/63) as of 23 February, 2023. In total, our sample included 9,869 cases infected with XBB/XBB.1.5 and 21,870 cases infected with other lineages. Characteristics of cases with XBB/XBB.1.5 and non-XBB/XBB.1.5 cases were similar, including their age, sex, and racial/ethnic distributions, as well as their prevalence of comorbid conditions and frequency of healthcare utilization over the preceding year (Table 1).
Vaccination history among cases
Comparing vaccination and infection history among individuals infected with distinct lineages provides a strategy to identify differences in the protective effectiveness of these exposures against each lineage.32,33 Importantly, this analytic approach may mitigate confounding that arises when comparing cases to uninfected controls,34 yielding estimates that indicate the relative degree of protection provided by vaccination or prior infection against each lineage. Within our sample, adjusted odds of having received 2, 3, 4, and ≥5 COVID-19 vaccine doses were 10% (95% confidence interval: 1-18%), 11% (3-19%), 13% (3-21%), and 25% (15-34%) lower among cases with XBB/XBB.1.5 in comparison to non-XBB/XBB.1.5 cases (Table 2), suggesting vaccine effectiveness against XBB/XBB.1.5 infection exceeded vaccine effectiveness against infection with co-circulating lineages. Analyses distinguishing vaccine doses received before or after individuals’ first documented SARS-CoV-2 infection yielded similar findings (Table S1).
We next aimed to investigate whether these differences in vaccination status among cases with XBB/XBB.1.5 and non-XBB/XBB.1.5 cases could be explained by differences in timing of vaccination, or in receipt of bivalent doses, among cases infected with each lineage (Table S2). However, these secondary analyses did not identify evidence that time since receipt of the last vaccine dose, or receipt of bivalent doses, differentially affected individuals’ risk of infection with XBB/XBB.1.5 and other lineages. Adjusted odds of having received ≥4 vaccine doses >90 days prior to the date of testing were 17% (8-25%) lower among cases with XBB/XBB.1.5 than among non-XBB/XBB.1.5 cases. Similarly, adjusted odds of having received ≥4 vaccine doses, but with ≥1 dose received <90 days before testing, were 15% (3-25%) lower among cases with XBB/XBB.1.5 than among non-XBB/XBB.1.5 cases. Consistent with these observations, distributions of times from receipt of the most recent vaccine dose were consistent for cases with XBB/XBB.1.5 and non-XBB/XBB.1.5 cases. As compared to non-XBB/XBB.1.5 cases, cases with XBB/XBB.1.5 had 13% (1-23%) lower adjusted odds of having received a vaccine series ≥4 monovalent doses only, and 18% (9-26%) lower adjusted odds of having received mixed monovalent/bivalent series with ≥4 doses. Taken together, these observations suggest that although COVID-19 vaccination was associated with greater degrees of protection against XBB/XBB.1.5 infection than infection with other co-circulating lineages, the added benefits of recent boosting,35–37 including with bivalent vaccine doses,30 may not differ appreciably for protection against XBB/XBB.1.5 as compared to co-circulating lineages.
History of documented SARS-CoV-2 infection among cases
In a pattern opposite to these observations, cases with XBB/XBB.1.5 had 17% (11-24%) and 40% (19-65%) higher adjusted odds of having experienced 1 and ≥2 prior documented infections, respectively, in comparison to non-XBB/XBB.1.5 cases (Table 2). We undertook additional analyses aiming to determine whether these findings were explained by recent BA.4/BA.5 infection, which may have provided enhanced and specific protection against closely related lineages co-circulating with XBB/XBB.1.5 during the study period (Table S3). Consistent with this hypothesis, cases with XBB/XBB.1.5 had 67% (40% to 99%) higher adjusted odds than non-XBB/XBB.1.5 cases of prior documented infection during the BA.4/BA.5 waves (25 June to 30 November, 2022), but did not have substantially higher adjusted odds of prior documented infection during the BA.1/BA.2 waves (20 December, 2021 to 24 June, 2022; 4% [–3% to 12%] higher adjusted odds). However, cases with XBB/XBB.1.5 also had higher adjusted odds of prior documented infection during periods predominated by pre-Omicron variants, including the Delta variant (18% [3-35%]), Alpha/Epsilon variants (19% [4-36%]), and earlier lineages (29% [18-41%]). These findings suggested that protection associated with recent BA.4/BA.5 infection could only partially account for the observed association of XBB/XBB.1.5 infection with prior documented infection; higher odds of prior infection with pre-Omicron variants among cases with XBB/XBB.1.5 suggest that XBB/XBB.1.5 effectively evades infection-derived immunity associated with other SARS-CoV-2 variants as well.
We also undertook several sensitivity analyses aiming to verify that associations of XBB/XBB.1.5 with prior infection were not an artifact of suboptimal vaccine response among certain recipients. First, we distinguished prior infections acquired after individuals had received ≥2 COVID-vaccine doses from prior infections occurring among individuals who had received only a single dose or who were never vaccinated. Whereas post-vaccination infections could be an indicator of poor vaccine response, the same would not be true of pre-vaccination infections (Figure S1). Within these analyses, cases with XBB/XBB.1.5 had 23% (14-31%) higher adjusted odds of having experienced pre-vaccination infection as well as 13% (5-22%) higher adjusted odds of having experienced post-vaccination infections in comparison to non-XBB/XBB.1.5 cases (Table S4). Further, we did not identify evidence of associations between infecting lineage and immunocompromised or immunosuppressed status (adjusted odds ratio 0.95 [0.89-1.05]), which would be expected to mediate differences in risk of post-vaccination infection among vaccine recipients (Table S5). These findings suggest the enhanced capacity of the XBB/XBB.1.5 lineage to overcome immune responses associated with prior infection is independent of the observed association between prior vaccination and XBB/XBB.1.5 infection. In further confirmation of this finding, point estimates indicated stepwise increases in adjusted odds of having experienced 1 or ≥2 prior documented infections among cases with XBB/XBB.1.5 within strata encompassing recipients of 1, 2, 3, and ≥4 COVID-19 vaccine doses (Table S6). These results support the directional relationship between infecting variant and history of documented SARS-CoV-2 infection identified in primary analyses.
Predictors of disease progression
Protection against progression to severe disease endpoints provides an additional dimension for measuring how immunity from vaccination or prior infection impacts SARS-CoV-2 natural history.34,38 Among cases with XBB/XBB.1.5 within our sample, prior receipt of 2, 3, and ≥4 COVID-19 vaccine doses was associated with 41% (–44% to 76%), 54% (0-79%), and 70% (30-87%) reductions, respectively, in adjusted hazards of progression to hospital admission due to any cause in the 30 days following a positive outpatient test (Table 3; single-dose effects were not estimated due to sparse counts). Corresponding estimates of protection against progression to hospital admission associated with 2, 3, and ≥4 COVID-19 vaccine doses were 6% (–65% to 46%), 46% (7-69%), and 48% (7-71%), respectively, among cases infected with other lineages. Our study was underpowered to estimate associations of vaccination with protection against progression to other severe disease endpoints including intensive care unit admission, mechanical ventilation, or mortality; no deaths occurred among cases infected with either lineage over the course of follow-up (Table 4).
Documented prior infections were not strongly associated with risk of hospital admission among cases with XBB/XBB.1.5 or non-XBB/XBB.1.5 cases (adjusted hazard ratio for hospital admission after ≥2 documented infections equal to 0.73 [0.17-3.15] and 0.72 [0.26-2.01], respectively; Table 3). Infection with XBB/XBB.1.5 or non-XBB/XBB.1.5 lineage was not independently associated with individuals’ likelihood of experiencing hospital admission or intensive care unit admission (adjusted hazard ratios equal to 1.03 [0.76 to 1.38] and 1.45 [0.51 to 4.13], respectively, for comparisons of cases with XBB/XBB.1.5 to non-XBB/XBB.1.5 cases; Table 4).
Discussion
Cases in our study infected with the XBB/XBB.1.5 lineage had received fewer COVID-19 vaccine doses and had higher likelihood of prior documented SARS-CoV-2 infection, in comparison to contemporaneous cases infected with other circulating SARS-CoV-2 lineages. These findings suggest that although the XBB/XBB.1.5 lineage has greater capacity than co-circulating lineages (predominantly descending from BA.5) to evade immune responses triggered by prior infection, including with pre-Omicron variants, XBB/XBB.1.5 is more sensitive to immune responses triggered by COVID-19 vaccination. In further support of this hypothesis, we identify that COVID-19 vaccination is associated with greater degrees of protection against hospital admission among cases with XBB/XBB.1.5 than among non-XBB/XBB.1.5 cases. Reassuringly, these findings support the persistent benefit of vaccination for reducing individuals’ risk of infection and severe disease with XBB/XBB.1.5, despite previously-reported immune-evasive characteristics of this lineage.28,29,39 As capture of prior infections in individuals’ electronic health records was likely incomplete for both cases with XBB/XBB.1.5 and non-XBB/XBB.1.5 cases, our study likely underestimates the true magnitude of association of prior infection and individuals’ relative likelihoods of infection with XBB/XBB.1.5 or other SARS-CoV-2 lineages.40 These results demonstrate the considerable ability of XBB/XBB.1.5 to evade immunity resulting from prior infection, in comparison to SARS-CoV-2 variants that have emerged during earlier phases of the pandemic.
Our finding that enhanced vaccine sensitivity co-occurs with enhanced escape of infection-derived immunity in XBB/XBB.1.5 stands in contrast to observations of other SARS-CoV-2 variants emerging after widespread vaccine implementation. The BA.1 and BA.4/BA.5 Omicron lineages were associated with reductions in the protective effectiveness of both prior vaccination and infection in comparison to the Delta and BA.2 Omicron lineages they outcompeted.10–21 Vaccine effectiveness against Delta variant infection was also mildly weaker than effectiveness against earlier variants,6,7 although emergence and dissemination of the Delta variant preceded widespread vaccination in most countries. While the greatest differences in lineage-specific vaccine protection have been apparent in effects on individuals’ risk of acquiring infection, previous studies disaggregating vaccine effectiveness against infection from vaccine effectiveness against hospital admission38 have demonstrated that vaccine escape may be associated with changes in protection against severe disease progression as well, consistent with findings in the present study. In early 2022, protection against progression from a positive outpatient test to hospital admission associated with prior receipt of ≥3 mRNA vaccine doses was 57% for cases within KPSC infected with the BA.1 Omicron lineage, versus 86% for cases infected with the Delta variant.10 In the present study, receipt of 2, 3, and ≥4 COVID-19 vaccine doses was associated with greater stepwise increases in point estimates of protection against hospital admission for cases with XBB/XBB.1.5 (41%, 54%, and 70%, respectively) in comparison to non-XBB/XBB.1.5 cases (6%, 46%, and 48%, respectively). It is important to note that these estimates represent only the reduction in risk of severe disease progression among cases experiencing breakthrough infection following vaccination, and do not account for additional vaccine-derived protection against acquisition of SARS-CoV-2 infection—an important but separate component of vaccine effectiveness against hospitalized illness.38 Taken together, our observations of differential vaccination history among cases with XBB/XBB.1.5 and non-XBB/XBB.1.5 cases, and differential protection against disease progression, are consistent with the hypothesis that XBB/XBB.1.5 shows enhanced sensitivity to vaccine-derived immune responses in comparison to other co-circulating lineages.
Immunological and evolutionary factors driving this apparent bifurcation in evasion of vaccine-derived and infection-derived responses for XBB/XBB.1.5 merit further investigation. Notably, vaccinations available in the US (mRNA-1273, BNT162b2, Ad.26.COV2.S, and NVX-CoV2373) target only the SARS-CoV-2 spike antigen. In contrast, infection with SARS-CoV-2 induces responses against an array of SARS-CoV-2 antigens, some of which may be independently associated with protection.41 Among US blood donors, seroprevalence of infection-derived (anti-nucleocapsid) antibodies reached 58% by February 2022,42 prior to widespread transmission of the BA.2 and BA.4/BA.5 variants. While US seroprevalence estimates following BA.4/BA.5 emergence are unavailable, studies in other countries have reported >80% prevalence of anti-nucleocapsid antibodies by late 2022, consistent with widespread transmission of these lineages.43,44 As prevalence of prior infection approaches or exceeds the proportion of individuals who have received all recommended primary series and booster doses of COVID-19 vaccines,27 the ability to evade responses against non-spike SARS-CoV-2 antigens might be of increasing importance to the fitness of emerging SARS-CoV-2 lineages. Enhanced infectivity, regardless of immune evasion, likewise contributes to the ability of novel SARS-CoV-2 lineages to achieve widespread transmission.8 However, as prior studies have indicated that host receptor binding affinity of XBB/XBB.1.5 is weaker than that observed in co-circulating BQ.1/BQ.1.1 lineages,29,39 evasion of infection-derived immunity from both BA.4/BA.5 and pre-Omicron lineages may of greater relative importance to the successful establishment of XBB/XBB.1.5.
Our study has at least five limitations. First, our comparison of cases infected with XBB/XBB.1.5 and other co-circulating lineages is observational in nature. Although cases with XBB/XBB.1.5 and non-XBB/XBB.1.5 cases did not differ appreciably in most demographic characteristics or measured risk factors, unmeasured differences between the two case populations remain possible sources of confounding, and demographic or epidemiologic characteristics of the study population may not be generalizable to cases in other settings. Relatedly, our study was restricted to cases who received molecular testing in clinical settings, who may be distinct from individuals who received point-of-care antigen testing or who did not seek testing for their infections. However, selecting on receipt of a positive molecular test result enabled us to mitigate potential sources of confounding present in other study designs comparing infected cases to uninfected controls; fewer factors would be expected to confound comparisons of cases with XBB/XBB.1.5 and non-XBB/XBB.1.5 cases.34 As our study conditions on individuals’ acquisition of SARS-CoV-2 infection, observed associations should be interpreted as measures of the relative effect of vaccination and prior infection on individuals’ risk of acquiring XBB/XBB.1.5, compared to their risk of acquiring other co-circulating lineages.45 Alternative designs are needed to quantify absolute effectiveness of prior infection and vaccination in preventing infection. Third, we relied on infections recorded in cases’ electronic health records to determine prior infection history. This is expected to result in exposure misclassification for both cases with XBB/XBB.1.5 and non-XBB/XBB.1.5 cases, as some prior infections may have gone undiagnosed or may have been diagnosed in settings outside the KPSC healthcare system. Our estimates should thus be interpreted as lower-bound measures of the association between prior infection and XBB/XBB.1.5 or non-XBB/XBB.1.5 infecting lineage. Fourth, our analyses of the association of infecting lineage with prior infection during periods predominated by transmission of BA.4/BA.5, BA.1/BA.2, Delta, and earlier lineages are subject to differing degrees of misclassification associated with changes over time in access to both clinical and at-home tests. While these analyses reveal that immune escape is not limited to evasion of responses associated with recent BA.4/BA.5 infection, comparisons of effect size estimates for prior infections with differing lineages should be made with caution. Relatedly, patients received diverse COVID-19 vaccine series (although Ad.26.COV2.S and NVX-CoV2373 were rarely administered within our study population); our focus on numbers of doses received, time since the final dose, and whether bivalent doses were received may not fully capture variation in vaccine-associated immune protection. Last, our study used all-cause hospital admissions to indicate severe disease progression. As some admissions may occur due to causes unrelated to COVID-19 among both cases with XBB/XBB.1.5 and non-XBB/XBB.1.5 cases, associations of vaccination with protection against disease progression may be underestimated. However, excluding cases tested in hospital settings helped to mitigate bias under our study design, as routine SARS-CoV-2 testing on admission may identify more substantial numbers of incidental COVID-19 admissions related to screening46,47 within studies including hospitalized patients.
Our analysis indicates increased vaccine sensitivity of the emerging SARS-CoV-2 XBB/XBB.1.5 lineage as well as enhanced ability of this lineage to evade immunity associated with prior infection, including with pre-Omicron SARS-CoV-2 lineages. Selection for evasion of immune responses associated with prior infection, such as those targeting non-spike SARS-CoV-2 antigens, might be of growing importance to SARS-CoV-2 evolutionary trajectories as immunity from prior SARS-CoV-2 infection becomes increasingly prevalent. While it is reassuring that prior vaccination is associated with enhanced protection against the XBB/XBB.1.5 lineage in comparison to co-circulating lineages, ongoing escape of infection-derived immunity remains a cause for concern. Continuous monitoring of changes in protection associated with vaccination and prior infection is needed to inform responses to emerging SARS-CoV-2 variants.
Methods
Study setting
We undertook this retrospective observational study within the Kaiser Permanente Southern California (KPSC) healthcare system. As a comprehensive, integrated care organization, KPSC delivers healthcare across telehealth, outpatient, emergency department, and inpatient settings for ∼4.7 million members enrolled through employer-provided, government-sponsored, and pre-paid coverage schemes. Electronic health records (EHRs) across all clinical settings, together with laboratory, pharmacy, and immunization data, provide a complete view into care delivered by KPSC. These observations are augmented by insurance claims for out-of-network diagnoses, prescriptions, and procedures, enabling near-complete capture of healthcare interactions for KPSC members. The KPSC Institutional Review Board reviewed and provided ethical approval for the study.
Eligibility criteria
We included cases who: (1) received a positive molecular test result in any outpatient setting between 1 December, 2022 and 23 February, 2023; (2) had specimens processed using the ThermoFisher TaqPath COVID-19 Combo Kit (as described below); (3) had not received any prior SARS-CoV-2 test result in a clinical setting or prior COVID-19 diagnosis within 90 days before their index test; (4) had been enrolled in KPSC health plans for ≥1 year prior to their index test, allowing for up to a 45-day lapse in membership to account for potential delays in re-enrollment; and (5) were not hospitalized at the time of their index test, and had not been hospitalized at any point in the preceding 7 days. Restricting analyses to outpatient cases was expected to provide several design advantages. First, this strategy helped to ensure cases infected with each lineage were similar to each other in terms of healthcare-seeking behavior.34,48 Second, initiating follow-up from the point of outpatient testing helped to ensure cases were ascertained at similar stages of their clinical course, facilitating unbiased comparisons of subsequent progression to severe disease. Outpatient-diagnosed cases at KPSC were automatically enrolled in a home-based symptom monitoring program with standardized criteria for emergency department referral and inpatient admission as a measure to preserve hospital capacity throughout the study period.49 Thus, hospital admission was considered to represent an internally consistent measure of disease severity within the sample followed from the point of outpatient testing,10 whereas cases first intercepted in hospital settings may have had greater variation in clinical status at the point of testing. Relatedly, this approach helped our study to avoid inclusion of incidental SARS-CoV-2 detections among patients who were tested at the point of inpatient admission for causes unrelated to COVID-19. Last, whereas the TaqPath COVID-19 Combo Kit was the primary assay used at regional reference laboratories for outpatient testing, cases diagnosed in hospital settings may have had tests processed in-house using other assays, without readout enabling lineage determination.
Lineage calling
The TaqPath COVID-19 Combo Kit assay included probes for the spike (S), nucleocapsid (N), and orf1a/b genes. Cases with cycle threshold (cT) values below 37 for ≥2 probes were considered positive for SARS-CoV-2. We interpreted S-gene target failure (SGTF), defined as cT≥37 for the S gene but cT<37 for the N and orf1a/b genes, as a proxy for infection with BA.4/BA.5 sublineages, whereas S-gene detection provided a proxy for infection with XBB/XBB.1.5 lineages, consistent with our validation data (Main text) and prior US studies.30,31
Exposures
We characterized the following attributes of included cases using data from their electronic health records: age (defined in 10-year age bands), biological sex; race/ethnicity (white, black, Hispanic of any race, Asian, Pacific Islander, and other/mixed/unknown race); neighborhood socioeconomic status, measured as the median household income within their Census block (<$30,000, $30,000-59,999, $60,000-89,999, $90,000-119,999, $120,000-149,999, and ≥$150,000 per year); cigarette smoking status (current, former, or never smoker); body mass index (BMI; categorized as underweight, normal weight, overweight, obese, or severely obese); Charlson comorbidity index (0, 1-2, 3-5, and ≥6); prior-year emergency department visits and inpatient admissions (each defined as 0, 1, 2, or ≥3 events); prior-year outpatient visits (0-4, 5-9, 10-14, 15-19, 20-29, or ≥30 events); documented prior SARS-CoV-2 infection; and history of COVID-19 vaccination (receipt of 0, 1, 2, 3, 4, or ≥5 doses, according to manufacturer, type, and time from receipt of each dose to the date of the index test). While our analyses include recipients of BNT162b2, mRNA-1273, Ad.26.COV2.S, and NVX-COV2373, we do not distinguish protection associated with receipt of mRNA or non-mRNA vaccine series, as non-mRNA vaccine doses accounted for only 2.0% of all vaccine doses received within the study population (1,741/86,076). For individuals with multiple prior COVID-19 diagnoses or positive SARS-CoV-2 test results, we considered these infections to be distinct if they were not preceded by any other COVID-19 diagnosis or positive test result within 90 days. For data anonymization, index test dates were jittered by random addition of –1, 0, or 1.
Outcomes
We followed outpatient-diagnosed cases from the point of their index test through death, disenrollment, or March 5, 2023, the date of the database cut (providing ≥10 days of follow-up for individuals who did not die or disenroll). The primary endpoint was hospital admission for any cause within 30 days after the index test. Additional endpoints monitored over 60 days after the index test included intensive care unit (ICU) admission, initiation of mechanical ventilation, and death.
Multiple imputation of missing data
To accommodate missing data on cases’ BMI (n=4,799; 15.1%), cigarette smoking (n=4,229; 13.3%), and neighborhood household income (n=6,476; 20.4%), we generated 10 pseudo-datasets completed by sampling from the conditional distribution of these variables, given all other observed characteristics of cases, via multiple imputation. We conducted complete-case statistical analyses across each of the 10 pseudo-datasets and pooled results across these analyses according to Rubin’s rules.50
Logistic regression analysis
Within this analytic sample of cases testing positive for SARS-CoV-2 infection, potential outcomes were binary (infection with XBB/XBB.1.5 or non-XBB/XBB.1.5 lineages). We estimated adjusted odds ratios of prior vaccination and prior documented infection among cases with XBB/XBB.1.5 and non-XBB/XBB.1.5 cases via logistic regression. Models used in primary analyses controlled for variables as categorized in Table 1, as motivated by a directed acyclic graph, and included fixed intercepts for calendar week of testing. For all analyses, we report unadjusted associations of each exposure with XBB/XBB.1.5 or non-XBB/XBB.1.5 infection, associations accounting only for testing week (“time-adjusted” odds ratios), and associations accounting for all confounders (“adjusted” odds ratios). Models using alternative adjustment strategies for vaccination and infection (Figure S1) provided results similar to those of the primary analysis (Table S1; Table S4; Table S6; Table S7).
Survival analysis
For analyses addressing the association of prior vaccination, infection, and infecting lineage with cases’ risk of hospital admission or ICU admission, we fit Cox proportional hazards models to data on cases’ times to each of these events or censoring (at study end date, end of follow-up at 30 or 60 days, or disenrollment, whichever occurred earliest). A survival analysis framework was motivated by the fact that XBB/XBB.1.5 infections accounted for an increasing share of all diagnosed cases over time, and thus had higher likelihood of censoring within <30 days or <60 days in comparison to non-XBB/XBB.1.5 cases. Models defined strata according to cases’ calendar week of testing to control for potential changes in testing and healthcare-seeking practices. We verified the proportional hazards assumption by testing for non-zero slopes of the Schoenfeld residuals.51
Software
We conducted analyses using R (version 4.0.3; R Foundation for Statistical Computing, Vienna, Austria). We used the survival52 package (version 3.5-3) for time-to-event analyses, and the Amelia II package53 (version 1.81.1) for multiple imputation.
Data availability statement
Individual-level testing and clinical outcomes data reported in this study are not publicly shared. Individuals wishing to access disaggregated data, including data reported in this study, should submit requests for access to the corresponding author (sara.y.tartof@kp.org). De-identified data (including, as applicable, participant data and relevant data dictionaries) will be shared upon approval of analysis proposals with signed data-access agreements in place.
Code availability statement
Analysis code is available from github.com/joelewnard/xbb.
Data Availability
Individual-level testing and clinical outcomes data reported in this study are not publicly shared. Individuals wishing to access disaggregated data, including data reported in this study, should submit requests for access to the corresponding author (sara.y.tartof{at}kp.org). De-identified data (including, as applicable, participant data and relevant data dictionaries) will be shared upon approval of analysis proposals with signed data-access agreements in place.
Author contributions statement
JAL, ML, and SYT contributed to the study concept and design. VH, JSK, SFS, BL, HT, and SYT led acquisition of data. VH and JAL led statistical analysis of data. JAL, ML, and SYT led interpretation of data. JAL drafted the manuscript, and all authors critically revised the manuscript for important intellectual content. SYT obtained funding and provided supervision.
Competing interests statement
JAL has received research grants paid directly to his institution and consulting honoraria unrelated to this study from Pfizer. SYT has received research grants paid directly to her institution unrelated to this study from Pfizer. The remaining authors disclose no competing interests.
Figure captions
Acknowledgments
This work was funded by the US Centers for Disease Control and Prevention (CDC; grant 75D3-121C11520 to SYT). The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the CDC. JAL was supported by grant R01-AI14812701A1 from the National Institute for Allergy and Infectious Diseases (US National Institutes of Health), which had no role in design or conduct of the study, or the decision to submit for publication.