Observed serial intervals of SARS-CoV-2 for the Omicron and Delta variants in Belgium based on contact tracing data, 19 November to 31 December 2021 ==================================================================================================================================================== * Cécile Kremer * Toon Braeye * Kristiaan Proesmans * Emmanuel André * Andrea Torneri * Niel Hens ## Abstract The SARS-CoV-2 Omicron BA.1 variant is rapidly spreading worldwide, possibly outcompeting the Delta strain. We investigated the empirical serial interval for both variants using contact tracing data. Overall, we observed a shorter serial interval for Omicron compared to Delta, suggesting faster transmission. Furthermore, results indicate a relation between the empirical serial interval and the vaccination status for both the Omicron and the Delta variant. Consequently, with the progression of the vaccination campaign, the reasons for and extent of dominance of Omicron over Delta may need further assessment. ## Background The WHO designated the SARS-CoV-2 Omicron BA.1 variant (B.1.1.529) as a variant of concern (VoC) on 26 November 2021 [1]. Omicron shows a fast epidemic growth and has taken over as dominant VoC from the previously dominant Delta variant (B.1.617.2) worldwide. In Belgium, the Omicron variant was the dominant circulating strain during the period from 27 December 2021 to 9 January 2022, identified in 88.5% of the sequenced samples [2]. It has been found that the Omicron variant is more efficient at evading immunity, acquired from previous infection or vaccination [3,4], compared to the Delta variant. Another epidemiological characteristic that may contribute to the rapid spread of Omicron is increased transmissibility, possibly due to an increase in the reproduction number or a shortened serial interval, i.e. the time between symptom onset in an infector and infectee [5]. In this study, we estimate the mean and the standard deviation of the serial interval for both the Omicron and Delta variants using contact tracing data collected in Belgium and assess whether these variants are associated with a different observed serial interval. To gain more insights on the possible impact of vaccination, we also compare the observed serial intervals for different combinations of vaccination status in transmission pairs. ### Data Belgium has a contact tracing system in place, where COVID-19 confirmed cases are asked about their contacts during two days prior to symptom onset until 10 days after. Variant detection was done using genotype sequencing. If a variant was found in a transmission chain, all cases belonging to that chain were assumed to be infected by that variant. We collected transmission pairs that could be linked either to Omicron or Delta infections, where the infector reported first symptoms between 19 November and 31 December 2021. During this period Omicron started to spread in Belgium and took over dominance from the Delta variant [6]. The same non-pharmaceutical interventions were in place throughout the considered time window, with a fairly low stringency index (around 48) compared to neighboring countries.We assumed that the first confirmed case in a reconstructed transmission pair (i.e. the ‘index’ case) was the infector, and hence the contact was the infectee. Transmission pairs for which symptom onset was not available for either the infector or the infectee were excluded, as well as pairs for which the observed serial interval was less than -5 days or more than 15 days in order to ensure biologically plausible serial intervals [7,8]. Vaccination status was assigned to both cases in a transmission pair as (i) unvaccinated (including partially vaccinated), (ii)vaccinated (i.e. completed vaccination cycle), and (iii) vaccinated plus booster. Household status was assigned based on information collected during the interview with the index case or provided by the healthcare worker. A more detailed description of the data and analysis is provided in Supplement S1. ### Overall serial interval Of the 2495 included transmission pairs, 86.61% were linked to transmission of the Omicron variant (Figure S1). In what follows we will report the mean and standard deviation (SD) of the observed serial intervals, the median was 3 days for all stratifications. All reported *p*-values are based on a Mann-Whitney U test. Tables S1 and S2 show the proportion of transmission pairs by household and vaccination status. We identified 2161 Omicron and 334 Delta transmission pairs during the period from 19 November to 31 December 2021. The empirical serial interval distribution for Omicron had a mean of 2.75 days (SD 2.53 days), compared to 3.00 days (SD 2.48 days) for Delta (*p* = .019; Figures 1A-B). Table 1 shows the parameter estimates of the normal distribution fit to both empirical serial interval distributions (Figure 1C). View this table: [Table 1.](http://medrxiv.org/content/early/2022/01/30/2022.01.28.22269756/T1) Table 1. Estimated parameters of a normal distribution for the serial interval of Omicron and Delta variants, by different stratifications. There were no Delta transmission pairs in which both cases had received a booster vaccine. Reported estimates are the posterior median and 95% credible interval (CrI). ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/30/2022.01.28.22269756/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2022/01/30/2022.01.28.22269756/F1) Figure 1. Empirical (A-B) and fitted normal (C) distribution of the serial intervals for Omicron and Delta variants. Onset date infector from 19 November until 31 December 2021. ### Serial interval stratified by household status There was a significant difference between the mean empirical serial interval of both variants within households, but not between households (*p* = .034 for within-household, and *p* = .686 for between-household transmission pairs; Figures 2A-D). Figure 2E represents the fit of a normal distribution to the empirical serial intervals for both variants, and in Table 1 we report the parameters of these fitted distributions. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/30/2022.01.28.22269756/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2022/01/30/2022.01.28.22269756/F2) Figure 2. Empirical (A-D) and fitted normal (E; solid = within-household, dashed = between-household) distribution of the serial intervals for Omicron and Delta variant, by household status (excluding 83 pairs for which no information on household status was available). Onset date infector from 19 November until 31 December 2021. ### Serial interval stratified by vaccination status No difference in empirical serial intervals was found for pairs where both cases were unvaccinated or only partially vaccinated (2.69 vs 2.54 days, *p* = .931; Figure 3A-B). For transmission pairs in which both cases were vaccinated (without booster), the empirical serial interval for Omicron was significantly shorter compared to Delta (2.63 vs 3.38 days, *p* = .004; Figure 3C-D). The mean empirical serial interval for Omicron was longer for pairs that received a booster vaccine compared to pairs that were vaccinated with only two doses (3.34 vs 2.63 days, *p* = .065; Figure 3C-E). The mean empirical serial interval was significantly longer for the Delta variant in transmission pairs in which both cases completed the vaccination cycle, compared to those where both cases were unvaccinated or only partially vaccinated (3.38 vs. 2.54 days, *p* = .045; Figure 3B-D). ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/30/2022.01.28.22269756/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2022/01/30/2022.01.28.22269756/F3) Figure 3. Empirical distribution of the serial intervals for Omicron (A, C, E) and Delta (B, D) variant, for transmission pairs where both cases are unvaccinated (A, B), vaccinated (C, D), or vaccinated with a booster (E, no data for Delta variant). Onset date infector from 19 November until 31 December 2021. ## Discussion Our estimates of the empirical serial interval for Omicron are in line with those previously reported. Lee et al. [9] reported a mean serial interval of 2.8 days. Based on 18 transmission pairs from South Korea, Kim et al. [10] estimated the mean serial interval to be 2.22 days (95% CrI 1.48 - 2.97) with a standard deviation of 1.62 days (95% CrI 0.87 - 2.37). Backer et al. [7] reported a mean serial interval of 3.5 and 3 days in two consecutive weeks for Omicron within-household pairs. In line with our findings, they found this to be shorter than for Delta pairs. Shorter serial intervals suggest a possibly shorter generation time for the Omicron variant compared to Delta, pointing to faster transmission. This could hence explain the rapid growth that is observed for the Omicron variant. However, control measures and asymptomatic transmission may lead to different serial and generation interval distributions [11]. Future studies estimating the generation time for both variants are needed to shed more light on this. We found an indication that the empirical serial interval might depend on vaccination status. Precisely, we observed a shorter empirical serial interval for Omicron (2.65 days) compared to Delta (3.38 days) in transmission pairs where both cases completed the vaccination cycle. In addition we found that for the Delta variant the empirical serial interval was longer in transmission pairs where both cases completed the vaccination cycle (3.38 days), compared to those who are unvaccinated or only partially vaccinated (2.54 days). We also found a longer empirical serial interval for Omicron in pairs where both cases already received a booster vaccine (3.34 days) compared to those pairs that were vaccinated but had not yet received a booster (2.63 days). These results indicate that the empirical mean serial interval increases when the infector and infectee have a higher level of vaccine-induced immunity. However, possibly because of limited sample size, this is not observed for all possible combinations of vaccination status (Figure S2) and more data are needed to properly assess the relation between vaccination and serial interval. Interestingly, the empirical mean serial interval for unvaccinated and vaccinated (without booster) transmission pairs was similar for the Omicron variant. If vaccine-induced immunity and serial interval are positively correlated, this result might be explained by lower vaccine efficacy against Omicron for individuals that have not yet received a booster vaccine. While the vaccination campaign progresses and more and more individuals receive a booster vaccine, the reasons for and extent of Omicron’s dominance over Delta might have to be reassessed once more data are available. This work has several limitations. Date of symptom onset is self-reported by the case and may be subject to recall bias. Likewise, there might be bias in the reporting of contacts and the level of reporting may differ for each index case. We have used all reported transmission pairs, although some of them may have been wrongly assigned. We further assume that directionality of transmission was from index to contact. We also do not explicitly account for right truncation [12], but this is assumed not to affect our estimates since it is unlikely that symptomatic infectees were missed as we used the data available on 17 January 2022, but only including infectors that had symptom onset until 31 December 2021 and limiting the serial interval to at most 15 days. However, because contacts are reconstructed until two days prior to symptom onset of the index case, possible left truncation might lead to exclusion of some transmission pairs. While the transmission risk is higher around symptom onset [8], it is possible that not all the infection events caused by an index case are detected. There is also possible selection bias due to the targeted genotype sequencing of suspected Omicron cases (such as previously infected cases or travelers), whereas genotype sequencing resulting in confirmed Delta cases might have been done on samples from a hospital setting as severe disease was an indication for sequencing during the study period. This analysis does not correct for age and previous infection and hence reinfections might be over-represented among the Omicron cases. ## Supporting information Supplementary Material [[supplements/269756_file03.pdf]](pending:yes) ## Data Availability Transmission pair data (without vaccination status for identifiability reasons) and R code for the analysis are available on https://github.com/cecilekremer/serial_interval [https://github.com/cecilekremer/serial\_interval](https://github.com/cecilekremer/serial_interval) ## Funding N.H. and A.T. acknowledge funding from the European Union’s Horizon 2020 research and innovation programme - project EpiPose (Grant agreement number 101003688). N.H. acknowledges funding from the Antwerp Study Centre for Infectious Diseases (ASCID) and the chair in evidence-based vaccinology at the Methusalem Centre of Excellence consortium VAX-IDEA. This project was supported by the VERDI project (101045989), funded by the European Union. Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union or the Health and Digital Executive Agency. Neither the European Union nor the granting authority can be held responsible for them. ## Ethical and privacy statement The data analyses carried out in this work are among the legal tasks of Sciensano (artikel 4, §4 Wet tot oprichting van Sciensano; article 4, §4 Law establishment Sciensano). Article 4, §4 explicitly states that Sciensano is authorized to collect, validate, analyse, report and archive data of a personal nature concerning public health. Sciensano is further authorized to make these processed data available with approval of the qualified sectoral committees. Such approval for these data was obtained from the Information Security Committee Social Security and Health (ISC). ## Contributions Study conceptualisation: NH, AT, CK, TB; Literature research: CK, AT, NH; Data collection: TB, KP; Data analysis code: CK; Data analysis: CK, AT, NH; Results interpretation: AT, CK, NH, EA, TB, KP; Manuscript writing: CK, AT, NH; Manuscript review: TB, KP, EA, NH, AT, CK; Coordination: NH. CK and TB contributed equally to this work. ## Data availability Transmission pair data (without vaccination status for identifiability reasons) and R code for the analysis are available on [https://github.com/cecilekremer/serial_interval](https://github.com/cecilekremer/serial_interval). ## Acknowledgements Data is provided by Sciensano, the Belgian institute for health. We thank Nicola Low for the insightful comments on an earlier draft of this manuscript. * Received January 28, 2022. * Revision received January 28, 2022. * Accepted January 30, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.World Health Organization. Tracking SARS-CoV-2 variants. 2021. [https://www.who.int/en/activities/tracking-SARS-CoV-2-variants](https://www.who.int/en/activities/tracking-SARS-CoV-2-variants). 2. 2.Sciensano. COVID-19 weekly report. 2022. [https://covid-19.sciensano.be/sites/default/files/Covid19/COVID-19\_Weekly\_report\_NL.pdf](https://covid-19.sciensano.be/sites/default/files/Covid19/COVID-19_Weekly_report_NL.pdf) 3. 3.Cele S, Jackson L, Khoury DS, Khan K, Moyo-Gwete T, Tegally H, et al. SARS-CoV-2 Omicron has extensive but incomplete escape of Pfizer BNT162b2 elicited neutralization and requires ACE2 for infection. medRxiv preprint. 2021; doi: [https://doi.org/10.1101/2021.12.08.21267417](https://doi.org/10.1101/2021.12.08.21267417). 4. 4.Pulliam JRC, van Schalkwyk C, Govender N, von Gottberg A, Cohen C, Groome MJ, et al. Increased risk of SARS-CoV-2 reinfection associated with emergence of the Omicron variant in Sourth Africa. medRxiv preprint. 2021; doi: [https://doi.org/10.1101/2021.11.11.21266068](https://doi.org/10.1101/2021.11.11.21266068). 5. 5.Svensson A. A note on generation times in epidemic models. Math Biosci. 2007;208(1):300–311. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.mbs.2006.10.010&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17174352&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F30%2F2022.01.28.22269756.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248196400017&link_type=ISI) 6. 6.Cuypers L, Maes P, Dellicour S, Van Weyenbergh J, Keyaerts E, Baele G. et al. Genomic surveillance of SARS-CoV-2 in Belgium. 2021; [https://assets.uzleuven.be/files/2021-12/genomic\_surveillance\_update\_211130.pdf](https://assets.uzleuven.be/files/2021-12/genomic_surveillance_update_211130.pdf). 7. 7.Backer JA, Eggink D, Andeweg SP, Veldhuijzen IK, van Maarseveen N, Vermaas K. et al. Shorter serial intervals in SARS-CoV-2 cases with Omicron variant compared to Delta variant in the Netherlands, 13 – 19 December 2021. medRxiv preprint; doi: [https://doi.org/10.1101/2022.01.18.22269217](https://doi.org/10.1101/2022.01.18.22269217). 8. 8.He X, Lau EHY, Wu P, Deng X, Wang J, Hao X. et al. Temporal dynamics in viral shedding and transmissibility of COVID-19. Nat Med. 2020;26(5):672–675. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M20-3012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F30%2F2022.01.28.22269756.atom) 9. 9.Lee JJ, Choe YJ, Jeong H, Kim M, Kim S, Yoo H, et al. Importation and transmission of SARS-CoV-2 B.1.1.529 (Omicron) variant of concern in Korea, November 2021. J Korean Med Sci. 2021;36(50):e346. 10. 10.Kim D, Jo J, Lim JS, Ryu S. Serial interval and basic reproduction number of SARS-CoV-2 Omicron variant in South Korea. medRxiv preprint. 2021; doi: : [https://doi.org/10.1101/2021.12.25.21268301](https://doi.org/10.1101/2021.12.25.21268301). 11. 11.Torneri A, Libin P, Scalia Tomba G, Faes C, Wood J, Hens N. On realized serial and generation intervals given control measures: the COVID-19 pandemic case. PLoS Comput Biol. 2021; doi: [https://doi.org/10.1371/journal.pcbi.1008892](https://doi.org/10.1371/journal.pcbi.1008892). 12. 12.Linton NM, Kobayashi T, Yang Y, Hayashi K, Akhmetzhanov AR, Jung SM, et al. Incubation period and other epidemiological characteristics of 2019 novel coronavirus infections with right truncation: a statistical analysis of publicly available case data. J Clin Med. 2020;9(2):538.