Serial Interval Distribution of SARS-CoV-2 Infection in Brazil

Text

The severe acute respiratory syndrome novel coronavirus 2 (SARS-CoV-2) first emerged in China in December 2019 [1] and was declared a pandemic by the World Health Organization on 11^th March 2020. The transmission dynamics, such as the serial interval – the time between symptom onset in a primary case and their secondary cases in a chain of transmission [2] – are still being established. The epidemic growth rate and the serial interval jointly determine the basic reproduction number [3], and the existence of pre-symptomatic transmission can be inferred when the serial interval is shorter than the incubation period.

Data for this analysis were provided by the Brazilian Ministry of Health following ethical approval (CONEP protocol number 30127020.0.0000.0068). All confirmed cases of SARS-CoV-2 infection notified on the REDCap system between 25^th February and 19^th March were analyzed. We received the linked pairs data (suspected and primary case) by the Brazilian Ministry of Health. We measured the serial intervals by computing the difference of dates of symptom onset.

In our dataset there were four negative serial intervals (symptom onset in the infectee preceded the infector) and seven zero-valued serial intervals (symptom onset on the same day for the infector-infectee pair). The mean (standard deviation) was 2.97 (3.29) days with a median of 3 days.

It is common to censor observed serial intervals to include only positive values; however, there is no strong theoretical justification for this. As such, we fit a serial interval probability distribution that allows for negative values by modelling the serial interval as where ⌊.⌋ is the floor operator, t_inf,1 is a random variable uniformly distributed in the interval [0,1], (Δt_lat,1, Δt_lat,2) are independent and identically distributed (i.i.d.) chi-squared random variables of mean 2.79 days and (Δt_int,1, Δt_int,2) are i.i.d. chi-squared random variables of mean 1.31 days. See the Supplemental Material for theoretical motivation and method.

Figure 1 shows the comparison between the measured and modelled serial intervals as well as means ± standard deviations for a number of serial interval probability distributions presented in the literature [4–9]. We conclude that our proposed model for the serial interval probability distribution approximates well the serial intervals measured from data.

• Download figure
• Open in new tab

Figure 1.

Upper panel – Serial interval estimation from Brazilian SARS-CoV-2 confirmed cases compared to the modelled distribution. Lower panel – mean ± standard deviation serial interval estimates from the current literature.

In addition, we provide two other analyses of the data in the Supplemental Material in which we i.) fit a normal distribution, and ii.) censor the data to include only positive values, best fit by a lognormal distribution.

The mean serial interval in Brazil of 2.97 days is the shortest reported to our knowledge, but we emphasize that the serial interval is not usually concentrated around the mean due to its high standard deviation (3.29 days) and non-symmetry. Du et al. [4] took an alternative approach fitting a normal distribution to 468 publicly reported serial intervals in which 12% had negative values – see Supplemental Material for comparison with our raw data. 8.38% of the measured serial intervals are negative, which is allowed by our model. Negative serial intervals are described in the literature [2] and can occur with pre-symptomatic transmission, and a wide range for the incubation period – the classical example being HIV.

In two transmissions clusters in Tianjin and Singapore the mean incubation periods were 9.0 and 7.1 days, respectively [5]. Our serial interval estimate is shorter than the reported incubation period, implying that the infectious period for SARS-CoV-2 begins before symptom onset. This is in line with recent observations of a large proportion of undocumented infections, complicating containment of the virus [10].

Our estimates are based on national case notification data, and this is a strength compared to datasets complied solely from media reports. However, there are at least two important sources of bias. Firstly, there is a tendency for secondary cases to recall more recent contacts – recency or recall bias – resulting in a shorter estimate of the serial interval. Similarly, contact tracing is likely to be most effective for recent contacts. Secondly, self-isolation following symptom onset will remove longer serial intervals that would have occurred due to transmission during the symptomatic phase. This second point highlights the role of contextual factors – isolation practices, population density, location of transmission, proportion of the population infected etc. – in determining the serial interval, which is not simply a biological constant of the virus. Given this consideration, it is important generate multiple estimates in different locations and stages of the epidemic. Our report is valuable as the first estimate from Brazil, and Latin America more generally.

Conflict of Interest

The authors declare no conflict of interest.