Quantifying the impact of immune history and variant on SARS-CoV-2 viral kinetics and infection rebound: A retrospective cohort study

  1. James A Hay
  2. Stephen M Kissler
  3. Joseph R Fauver
  4. Christina Mack
  5. Caroline G Tai
  6. Radhika M Samant
  7. Sarah Connolly
  8. Deverick J Anderson
  9. Gaurav Khullar
  10. Matthew MacKay
  11. Miral Patel
  12. Shannan Kelly
  13. April Manhertz
  14. Isaac Eiter
  15. Daisy Salgado
  16. Tim Baker
  17. Ben Howard
  18. Joel T Dudley
  19. Christopher E Mason
  20. Manoj Nair
  21. Yaoxing Huang
  22. John DiFiori
  23. David D Ho
  24. Nathan D Grubaugh  Is a corresponding author
  25. Yonatan H Grad  Is a corresponding author
  1. Harvard TH Chan School of Public Health, United States
  2. Yale School of Public Health, United States
  3. University of Nebraska Medical Center, United States
  4. IQVIA, Real World Solutions, United States
  5. Duke Center for Antimicrobial Stewardship and Infection Prevention, United States
  6. Tempus Labs, United States
  7. Vagelos College of Physicians and Surgeons, Columbia University, United States
  8. Hospital for Special Surgery, United States
  9. National Basketball Association, United States
21 figures, 12 tables and 1 additional file

Figures

PCR Ct value trajectories for confirmed and suspected infections.

(A) PCR Ct value trajectories for each acute Delta (red), Omicron BA.1 (blue), and other (black) infection. Individuals are grouped by the gap between detection and their most recent negative or …

Posterior estimates from a generalized linear model predicting probability of Ct value <30 with spline terms for the interaction between days since detection with age group and the interaction between days since detection with vaccination status and variant.

Shown are the marginal effects of (A) vaccination status and (B) age group on the proportion of individuals with Ct <30 on each day post detection after conditioning on being an Omicron BA.1 …

Effect of booster status, antibody titer against the founder SARS-CoV-2 strain and age group on the probability of Ct value <30 over days post detection.

(A) Proportion of Omicron BA.1 infections with Ct value <30 on each day post detection, stratified by age group and the interaction of booster status at the time of infection and antibody titer …

Estimated viral trajectories by variant and titer.

Points depict measured Ct values, lines depict the estimated population mean viral trajectories, and shaded regions depict the 95% credible intervals for the estimated population viral trajectories. …

Appendix 1—figure 1
Summary of cohort.

Top row describes cohort demographics and data on immune histories. Middle row describes infection data. Bottom row provides additional information on the infection data.

Appendix 1—figure 2
Summary of all infection and vaccination events over time, in addition to serum sample collection dates, included in the dataset.

(A) Histogram showing the distribution of vaccination dates (note that most first doses were administered prior to 2021-06-25), showing when in time the majority of individuals were vaccinated. (B) …

Appendix 1—figure 3
Frequency and proportion of SARS-CoV-2 variant detections.

(A) Frequency of sequenced and unsequenced detected infections over time by week. Vertical dashed lines and shaded backgrounds demarcate periods of variant dominance.

Appendix 1—figure 4
Correlation between authentic virus neutralization assay (ID50) and the Diasorin antibody titer against (A) wildtype and (B) Delta.

Horizontal yellow bar shows an ID50 titer of 50 and 100 respectively, Diagonal lines and shaded regions show mean and 95% confidence intervals (CI) for a linear regression between the Diason …

Appendix 1—figure 5
Distribution of delays from detection to symptom onset among individuals with known symptom status.

Dashed lines mark the median delay between detection and symptom onset. Solid lines mark the day of detection (0).

Appendix 1—figure 6
Distribution of delays from symptom onset to peak Ct values among individuals with known symptom status.

Dashed lines mark the median delay between detection and symptom onset. Solid lines mark the day of symptom onset (0).

Appendix 1—figure 7
Proportion of infections with Ct <30 on each day post detection stratified by detection speed.

Solid colored lines and shaded ribbons show posterior means and 95% credible intervals from a generalized linear model predicting probability of Ct <30 as a function of days since detection. Dashed …

Appendix 1—figure 8
All viral trajectories classified as rebound shown in Figure 1B.

Subplots are colored by the most stringent definition for rebound. To be included here, individuals must have 2+consecutive days of Ct ≥30 after an initial Ct <30. The vertical red dotted line marks …

Appendix 1—figure 9
Proportion of infections with Ct <30 on each day post detection by confirmed or suspected variant, vaccination status and detection group, conditional on being in the <30 years age group.

Solid colored lines and shaded ribbons are posterior estimates from a generalized linear model predicting probability of Ct <30 as a function of days since detection and vaccination status, showing …

Appendix 1—figure 10
Identical to Appendix 1—figure 9, but after excluding data from all players.
Appendix 1—figure 11
Proportion of Omicron BA.1 infections, stratified by symptom status, with Ct <30 on each day post detection by booster status and detection group.

Solid colored lines and shaded ribbons are posterior estimates from a generalized linear model predicting probability of Ct <30 as a function of days since detection, variant (only BA.1 results …

Appendix 1—figure 12
Proportion of Omicron BA.1 infections with Ct <30 on each day post detection stratified by vaccination status, age group and detection group.

Solid colored lines and shaded ribbons are posterior estimates from a generalized linear model predicting probability of Ct <30 as a function of days since detection, vaccination status and variant, …

Appendix 1—figure 13
Individual antibody titers to the ancestral SARS-CoV-2 spike over time.

(A) Measured antibody titers by date of sample collection. Lines show longitudinal samples from the same individual, colored by the most recent exposure at the time of sample collection. Lines going …

Appendix 1—figure 14
Histogram of time between (A) second vaccine dose and antibody titer measurement and (B) booster dose and antibody titer measurement.

Dashed line marks the median lag (162 days). 1 individual was infected between receiving their second vaccine dose and having a titer measurement taken (Delta infection). 42 individuals were …

Appendix 1—figure 15
Distribution of antibody titers among Omicron BA.1-infected individuals (colored points) stratified by age group and vaccination status at the time of infection.

Shown are mean titers (large black point) and bootstrapped 95% confidence intervals for the mean (black bars). Note that stratification is by infection and not individual, and that antibody titers …

Appendix 1—figure 16
Proportion of BA.1 Omicron infections with Ct <30 on each day post detection stratified by the interaction of booster status at the time of infection and antibody titer group (note, not stratified by age group here due to small subgroup size).

Shown are posterior estimates from a generalized linear model predicting probability of Ct <30 with a spline term on days since detection, conditional on titer/vaccination status category. Solid …

Appendix 1—figure 17
Schematic diagram of the likelihood function for viral kinetic inference.

The plot depicts the likelihood as a function of ΔCt, the difference between the observed Ct value and the limit of detection, so that ΔCt = 0 (the origin) represents observations at the limit of …

Tables

Table 1
Characteristics of the documented infections.

Counts (N) correspond to numbers of infections.

VariableCategoryNPercent
Total1280100
VariantDelta18014.1
Omicron BA.187868.6
Other22217.3
Vaccination StatusUnvaccinated22817.8
First dose60.5
Second dose42032.8
Boosted62648.9
Antibody Titer13–25047337.0
250–80050439.4
Unknown30323.7
SymptomaticNo25720.1
Yes66451.9
Unknown35928
Detection SpeedDelayed detection87768.5
Frequent testing40331.5
Age Group0–3055643.4
31–5056844.4
50+15512.1
Unknown10.1
Cumulative Infection Number1112888.1
214911.6
320.2
410.1
Days Since Previous ExposureNaïve (no prior exposure)22017.2
<1 month27321.3
1–3 months40331.5
>3 months38430.0
Table 2
Number of rebound infections classified under different definitions for initial clearance and subsequent rebound.
Initial clearance duration (consecutive days with Ct ≥30)Rebound duration (days above Ct value threshold)Ct value threshold of reboundReboundsTotalPercentage
≥4≥4Ct <3007490.00%
≥4≥3Ct <3017490.13%
≥4≥2Ct <3047490.53%
≥3≥4Ct <3029990.20%
≥3≥3Ct <3079990.70%
≥3≥2Ct <30169991.60%
≥2≥4Ct <30713340.53%
≥2≥3Ct <301813341.35%
≥2≥2Ct <304013343.00%
≥4≥4Ct <2507490.00%
≥4≥3Ct <2507490.00%
≥4≥2Ct <2507490.00%
≥3≥4Ct <2519990.10%
≥3≥3Ct <2519990.10%
≥3≥2Ct <2529990.20%
≥2≥4Ct <25113340.08%
≥2≥3Ct <25213340.15%
≥2≥2Ct <25513340.38%
Appendix 1—table 1
Number of identified rebounds stratified by variant, either confirmed through sequencing or assumed based on detection date.

Rebounds are defined here as any trajectory with an initial Ct <30, followed by a sequence of two or more consecutive negative tests or tests with Ct ≥30, and subsequently followed by two or more …

VariantReboundsTotal infectionsPercentage rebounded
Omicron368774.10%
Delta11780.562%
Other32791.08%
Appendix 1—table 2
Number of identified rebounds stratified by vaccination status.

Rebounds are defined here as any trajectory with an initial Ct <30, followed by a sequence of two or more consecutive negative tests or tests with Ct ≥30, and subsequently followed by two or more …

Vaccination statusReboundsTotal infectionsPercentage rebounded
Boosted324946.48%
No record53981.26%
Second dose33230.929%
Unvaccinated21591.26%
Appendix 1—table 3
Comparison of linear logistic regression models predicting probability of Ct <30 on each day since detection among individuals in the frequent testing group.

Models are ranked based on their expected log pointwise predictive density (ELPD), where a lower ELPD implies better predictive accuracy. Model weight refers to the weight of each model in a …

ModelELPD differenceSE differenceWeightAUCClassification accuracyAccuracy (≥30)Accuracy (<30)
Cumulative number of exposures, variant and age0.0000.0000.1570.8970.8280.8900.647
Days since previous exposure, variant and age–1.81310.0110.2990.8960.8230.8780.664
Vaccination status, variant and age–7.4409.8280.0890.8960.8240.8810.656
Cumulative number of exposures and age–27.9358.8230.1040.8940.8230.8870.637
Cumulative number of exposures and variant–31.5369.4280.0000.8940.8240.8910.630
Vaccination status and age–33.93312.3890.0000.8930.8220.8820.646
Vaccination status and variant–35.27413.6930.1560.8930.8240.8860.645
Days since previous exposure and variant–39.75014.2060.0000.8910.8220.8830.644
Days since previous exposure and age–40.92912.5030.0000.8920.8200.8770.654
Variant and age–52.43711.6810.0000.8900.8170.8840.624
Cumulative number of exposures–53.66112.3440.1430.8920.8220.8950.610
Vaccination status–58.50315.1620.0000.8900.8220.8880.632
Age–60.16012.6880.0000.8900.8180.8820.632
Days since previous exposure–69.84215.5470.0450.8880.8170.8840.624
Variant–87.42715.5540.0070.8870.8180.8890.611
Baseline–93.08715.8430.0000.8860.8170.8900.606
Appendix 1—table 4
Comparison of linear logistic regression models predicting probability of Ct <30 on each day since detection among individuals in the delayed testing group.

Models are ranked based on their expected log pointwise predictive density (ELPD), where a lower ELPD implies better predictive accuracy. Model weight refers to the weight of each model in a …

ModelELPD differenceSE differenceWeightAUCClassification accuracyAccuracy (≥30)Accuracy (<30)
Cumulative number of exposures, variant and age0.0000.0000.1540.9140.8500.8930.712
Vaccination status, variant and age–1.8177.8070.0000.9130.8490.8900.720
Vaccination status and age–6.0618.8990.1340.9130.8490.8880.724
Cumulative number of exposures and age–7.9004.7100.0070.9130.8490.8930.712
Days since previous exposure, variant and age–8.0468.9290.2880.9130.8490.8900.721
Days since previous exposure and age–20.1909.3600.0000.9120.8500.8900.725
Variant and age–22.4578.6280.0810.9120.8500.8880.728
Vaccination status and variant–59.93614.9860.1620.9090.8470.8880.717
Cumulative number of exposures and variant–60.13912.0980.0070.9100.8460.8910.703
Age–63.49512.8340.0520.9100.8490.8820.743
Cumulative number of exposures–65.89312.8510.1140.9100.8460.8910.707
Vaccination status–67.52015.4180.0000.9090.8470.8880.717
Days since previous exposure and variant–82.79815.9860.0000.9090.8460.8880.714
Days since previous exposure–93.53716.0580.0000.9090.8460.8870.718
Variant–107.36215.8000.0000.9070.8450.8830.723
Baseline–148.35018.6080.0000.9050.8420.8800.722
Appendix 1—table 5
Posterior estimates of viral trajectory attributes by variant and vaccination status.

Estimates are posterior means with 95% credible intervals.

Trajectory valueVariant/vaccination statusEstimate
Peak viral load (Ct)Other: Unvaccinated25.0 (24.2, 25.9)
Delta: 1–2 doses22.4 (21.4, 23.5)
Omicron: 1–2 doses25.6 (25.0, 26.2)
Omicron: Boosted25.7 (25.4, 26.1)
Peak viral load (log10 copies/ml)Other: Unvaccinated6.8 (6.6, 7.0)
Delta: 1–2 doses7.5 (7.2, 7.8)
Omicron: 1–2 doses6.6 (6.5, 6.8)
Omicron: Boosted6.6 (6.5, 6.7)
Proliferation time (days)Other: Unvaccinated3.5 (3.1, 3.9)
Delta: 1–2 doses3.6 (3.1, 4.1)
Omicron: 1–2 doses3.6 (3.3, 4.0)
Omicron: Boosted4.0 (3.8, 4.3)
Clearance time (days)Other: Unvaccinated9.9 (9.2, 10.6)
Delta: 1–2 doses7.6 (7.0, 8.3)
Omicron: 1–2 doses6.2 (5.8, 6.6)
Omicron: Boosted8.4 (8.0, 8.7)
Appendix 1—table 6
Posterior estimates of viral trajectory attributes by variant and antibody titer.

Estimates are posterior means with 95% credible intervals.

Trajectory attributeVariant/immune statusEstimate
Peak viral load (Ct)Other: Unexposed24.8 (24.0, 25.6)
Delta: Exposed,≤250 AU22.1 (20.6, 23.4)
Delta: Exposed,>250 AU24.4 (22.7, 26.1)
Omicron: Exposed,≤250 AU25.2 (24.7, 25.6)
Omicron: Exposed,>250 AU26.2 (25.7, 26.6)
Peak viral load (log10 copies/ml)Other: Unexposed6.9 (6.6, 7.1)
Delta: Exposed,≤250 AU7.6 (7.3, 8.0)
Delta: Exposed,>250 AU7.0 (6.5, 7.4)
Omicron: Exposed,≤250 AU6.8 (6.7, 6.9)
Omicron: Exposed,>250 AU6.5 (6.4, 6.6)
Proliferation time (days)Other: Unexposed3.5 (3.2, 3.9)
Delta: Exposed,≤250 AU3.6 (2.8, 4.5)
Delta: Exposed,>250 AU4.2 (3.3, 5.3)
Omicron: Exposed,≤250 AU3.9 (3.6, 4.2)
Omicron: Exposed,>250 AU3.7 (3.5, 4.0)
Clearance time (days)Other: Unexposed9.8 (9.1, 10.6)
Delta: Exposed,≤250 AU7.7 (6.9, 8.7)
Delta: Exposed,>250 AU7.6 (6.4, 8.8)
Omicron: Exposed,≤250 AU8.4 (8.0, 8.8)
Omicron: Exposed,>250 AU6.9 (6.6, 7.3)
Appendix 1—table 7
Posterior estimates of Omicron BA.1 viral trajectory attributes by symptom and vaccination status.

Estimates are posterior means with 95% credible intervals.

Trajectory attributeVariant/immune statusEstimate
Peak viral load (Ct)Omicron: 1–2 doses, no symptoms26.7 (25.7, 27.7)
Omicron: 1–2 doses, symptoms25.3 (24.5, 26.1)
Omicron: boosted, no symptoms26.3 (25.7, 27)
Omicron: boosted, symptoms25.4 (25, 25.8)
Peak viral load (log10 copies/ml)Omicron: 1–2 doses, no symptoms6.3 (6.1, 6.6)
Omicron: 1–2 doses, symptoms6.7 (6.5, 6.9)
Omicron: boosted, no symptoms6.4 (6.3, 6.6)
Omicron: boosted, symptoms6.7 (6.6, 6.8)
Proliferation time (days)Omicron: 1–2 doses, no symptoms3.7 (3.2, 4.2)
Omicron: 1–2 doses, symptoms3.7 (3.2, 4.2)
Omicron: boosted, no symptoms4.3 (3.9, 4.8)
Omicron: boosted, symptoms3.8 (3.5, 4.1)
Clearance time (days)Omicron: 1–2 doses, no symptoms5.9 (5.2, 6.6)
Omicron: 1–2 doses, symptoms6.3 (5.8, 6.8)
Omicron: boosted, no symptoms7.4 (6.8, 8.0)
Omicron: boosted, symptoms8.7 (8.3, 9.1)
Appendix 1—table 8
Posterior estimates of Omicron BA.1 viral trajectory attributes by symptom and antibody titer.

Estimates are posterior means with 95% credible intervals.

Trajectory attributeVariant/immune statusEstimate
Peak viral load (Ct)Omicron: Low titer, no symptoms26.1 (25.2, 26.9)
Omicron: Low titer, symptoms24.8 (24.3, 25.4)
Omicron: High titer, no symptoms26.7 (25.9, 27.5)
Omicron: High titer, symptoms25.9 (25.4, 26.5)
Peak viral load (log10 copies/ml)Omicron: Low titer, no symptoms6.5 (6.3, 6.8)
Omicron: Low titer, symptoms6.9 (6.7, 7.0)
Omicron: High titer, no symptoms6.3 (6.1, 6.6)
Omicron: High titer, symptoms6.6 (6.4, 6.7)
Proliferation time (days)Omicron: Low titer, no symptoms4.3 (3.8, 4.9)
Omicron: Low titer, symptoms3.6 (3.3, 4.0)
Omicron: High titer, no symptoms3.8 (3.4, 4.3)
Omicron: High titer, symptoms3.8 (3.4, 4.2)
Clearance time (days)Omicron: Low titer, no symptoms7.4 (6.7, 8.1)
Omicron: Low titer, symptoms8.7 (8.2, 9.3)
Omicron: High titer, no symptoms6.4 (5.8, 7.1)
Omicron: High titer, symptoms7.1 (6.6, 7.5)
Appendix 1—table 9
Posterior estimates of Omicron BA.1 viral RNA clearance times by age and vaccination status.

Estimates are posterior means with 95% credible intervals. Low titer is ≤250 AU, high titer is >250 AU.

Trajectory attributeVariant/immune statusEstimate
Peak viral load (Ct)Omicron:<30, unboosted25.6 (25, 26.3)
Omicron:<30, boosted25.8 (25.1, 26.4)
Omicron: 30–50, unboosted25.3 (24.2, 26.5)
Omicron: 30–50, boosted25.7 (25.2, 26.2)
Omicron: 50+, unboosted25.5 (22.8, 27.9)
Omicron: 50+, boosted25.5 (24.6, 26.3)
Peak viral load (log10 copies/ml)Omicron:<30, unboosted6.6 (6.4, 6.8)
Omicron:<30, boosted6.6 (6.4, 6.8)
Omicron: 30–50, unboosted6.7 (6.4, 7)
Omicron: 30–50, boosted6.6 (6.5, 6.7)
Omicron: 50+, unboosted6.7 (6, 7.4)
Omicron: 50+, boosted6.7 (6.5, 6.9)
Proliferation time (days)Omicron:<30, unboosted3.7 (3.4, 4.1)
Omicron:<30, boosted4.2 (3.8, 4.7)
Omicron: 30–50, unboosted3.6 (3, 4.4)
Omicron: 30–50, boosted4 (3.6, 4.3)
Omicron: 50+, unboosted4.4 (3, 6.3)
Omicron: 50+, boosted3.9 (3.4, 4.4)
Clearance time (days)Omicron:<30, unboosted6 (5.6, 6.5)
Omicron:<30, boosted7.2 (6.7, 7.7)
Omicron: 30–50, unboosted6.4 (5.7, 7.2)
Omicron: 30–50, boosted8.6 (8.1, 9.1)
Omicron: 50+, unboosted8.7 (6.6, 11.2)
Omicron: 50+, boosted9.6 (8.8, 10.6)
Appendix 1—table 10
Posterior estimates of Omicron BA.1 viral RNA clearance times by age and titer.

Estimates are posterior means with 95% credible intervals. Low titer is ≤250 AU, high titer is >250 AU.

Trajectory attributeVariant/immune statusEstimate
Peak viral load (Ct)Omicron:<30, low titer24.8 (24, 25.5)
Omicron:<30, high titer26.2 (25.6, 26.8)
Omicron: 30–50, low titer25.4 (24.8, 26)
Omicron: 30–50, high titer26.1 (25.4, 26.8)
Omicron: 50+, low titer24.9 (23.8, 25.9)
Omicron: 50+, high titer26.3 (24.8, 27.6)
Peak viral load (log10 copies/ml)Omicron:<30, low titer6.9 (6.7, 7.1)
Omicron:<30, high titer6.5 (6.3, 6.7)
Omicron: 30–50, low titer6.7 (6.5, 6.9)
Omicron: 30–50, high titer6.5 (6.3, 6.7)
Omicron: 50+, low titer6.8 (6.6, 7.1)
Omicron: 50+, high titer6.5 (6.1, 6.9)
Proliferation time (days)Omicron:<30, low titer3.9 (3.4, 4.3)
Omicron:<30, high titer4 (3.6, 4.4)
Omicron: 30–50, low titer4 (3.6, 4.4)
Omicron: 30–50, high titer3.6 (3.2, 4)
Omicron: 50+, low titer3.8 (3.2, 4.5)
Omicron: 50+, high titer3.7 (2.9, 4.7)
Clearance time (days)Omicron:<30, low titer7.2 (6.7, 7.8)
Omicron:<30, high titer6.2 (5.8, 6.7)
Omicron: 30–50, low titer8.7 (8.1, 9.3)
Omicron: 30–50, high titer7.4 (6.8, 8)
Omicron: 50+, low titer9.9 (8.7, 11.1)
Omicron: 50+, high titer9.5 (8, 11.1)

Additional files

Download links