The contribution of asymptomatic SARS-CoV-2 infections to transmission on the Diamond Princess cruise ship

Abstract
Introduction
Results
Discussion
Materials and methods
Appendix 1
Appendix 2
Data availability
References
Article and author information
Metrics

Abstract

A key unknown for SARS-CoV-2 is how asymptomatic infections contribute to transmission. We used a transmission model with asymptomatic and presymptomatic states, calibrated to data on disease onset and test frequency from the Diamond Princess cruise ship outbreak, to quantify the contribution of asymptomatic infections to transmission. The model estimated that 74% (70–78%, 95% posterior interval) of infections proceeded asymptomatically. Despite intense testing, 53% (51–56%) of infections remained undetected, most of them asymptomatic. Asymptomatic individuals were the source for 69% (20–85%) of all infections. The data did not allow identification of the infectiousness of asymptomatic infections, however low ranges (0–25%) required a net reproduction number for individuals progressing through presymptomatic and symptomatic stages of at least 15. Asymptomatic SARS-CoV-2 infections may contribute substantially to transmission. Control measures, and models projecting their potential impact, need to look beyond the symptomatic cases if they are to understand and address ongoing transmission.

Introduction

The ongoing COVID-19 pandemic has spread rapidly across the globe, and the number of individuals infected with SARS-CoV-2 outstrips the number of reported cases (Wang et al., 2020; Golding et al., 2020). One key reason for this may be that a substantial proportion of cases proceed asymptomatically, that is, they either do not experience or are not aware of symptoms throughout their infection but may still transmit to others. In this sense, asymptomatic infections differ from presymptomatic infections, which describes the portion of the incubation period before symptoms develop during which onward transmission is possible.

While pre- and asymptomatic individuals do not directly contribute to morbidity or mortality in an outbreak, they can contribute to ongoing transmission, as has been shown for COVID-19 (Rothe et al., 2020; Chen et al., 2020; Ganyani et al., 2020) and other diseases (Dean et al., 2016; Slater et al., 2019; Esmail et al., 2018). In particular, purely symptom-based interventions (e.g., self-isolation upon onset of disease) will not interrupt transmission from asymptomatic individuals and hence may be insufficient for outbreak control if a substantial proportion of transmission originates from pre- and asymptomatic infections (Chen et al., 2020).

An estimate of the proportion of infections that progress to symptomatic disease, also known as the case-to-infection ratio, provides an indicator of what proportion of infections will remain undetected by symptom-based case detection (Salomon and COVID-19 Statistics, Policy modeling, and Epidemiology Collective, 2020). Evidence so far have suggested that the proportion of SARS-CoV-2 infections that proceed asymptomatically is likely non-trivial (Li et al., 2020; Liu et al., 2020a; He et al., 2020; Mizumoto et al., 2020; Lavezzo et al., 2020; Bendavid et al., 2020), although empirical data are often difficult to interpret due to opportunistic sampling frames (Nishiura et al., 2020) combined with low (Fontanet et al., 2020) and imbalanced participation from individuals with and without symptoms (Gudbjartsson et al., 2020). While it is likely that transmission from asymptomatic individuals can occur, (Bai et al., 2020) quantitative estimates are effectively absent. Improved understanding of the relative infectiousness of asymptomatic SARS-CoV-2 infection, and its contribution to overall transmission will greatly improve the ability to estimate the impact of intervention strategies (Salomon and COVID-19 Statistics, Policy modeling, and Epidemiology Collective, 2020). What is known is that in the presence of active case-finding, presymptomatic infections and symptomatic cases contribute almost equally to overall transmission, as both modelling and empirical studies have shown (Liu et al., 2020a; He et al., 2020).

Documented outbreaks in a closed population with extensive testing of individuals regardless of symptoms provide unique opportunities for improved insights into the dynamics of an infection, as knowledge of the denominator and true proportion infected are crucial, yet often unavailable in other datasets. Here, we use data from the well-documented outbreak on the Diamond Princess cruise ship to capture the mechanics of COVID-19 in a transmission model with explicit asymptomatic and presymptomatic states to infer estimates for the proportion, infectiousness and contribution to transmission of asymptomatic infections. Available data included the date of symptom onset for symptomatic disease for passengers and crew, the number of symptom-agnostic tests administered each day, and the date of positive tests for asymptomatic and presymptomatic individuals (Mizumoto et al., 2020; Nishiura, 2020; NIID, 2020).

Results

Model calibration

The model reflected the data well (Figure 1), including the differently timed peaks for confirmed symptomatic cases for crew (Figure 1A) and passengers (Figure 1B). In addition, the model matched the expected impact of quarantine of passengers on transmission from February 4th as illustrated by the drop in reproductive number (Figure 1E), followed by a later drop in transmission after February 10th, which was driven by a change in contact pattern in crew. See Figure 1—figure supplements 1–2 for full calibration outputs.

Figure 1 with 2 supplements see all

Download asset Open asset

Data from the Diamond Princess and model calibration.

Figure shows data from the Diamond Princess (points (**A-D**) and bars (F)) and results from model calibration. Red lines = median, shading = 95% posterior plus observational interval (**A-C**) and 95% posterior interval only (**D-E**). Two vertical lines show the date of the first confirmed diagnosis (left) and the start of quarantine measures (right). (**A-B**) show confirmed symptomatic cases among crew (A) and passengers (B) with a reported date of onset; (C) shows confirmed pre- or asymptomatic individuals by test date; (D) shows the prevalence of pre/asymptomatic individuals by test date. Points and error bars show point estimates and 95% confidence intervals; (E) shows the basic reproduction number over time for the ship as a whole, reflecting the drop in contact rates (F) shows the number of tests administered irrespective of symptoms, by test date.

Figure 1—source data 1 Marginal posterior parameter values from model calibration. Median and equal-tailed 95% marginal posterior parameter values estimated using 100,000 samples from the joint posterior distribution found during MCMC model calibration. See Table 2 for parameter definitions and descriptions.: https://cdn.elifesciences.org/articles/58699/elife-58699-fig1-data1-v2.xlsx
Download elife-58699-fig1-data1-v2.xlsx

Asymptomatic infections

We estimated that 74% of infections proceeded asymptomatically (70–78%, 95% Posterior Interval (PI)) (see Figure 2A). The strong identifiability of this parameter is driven by the relative proportions of individuals testing positive with and without symptoms, combined with the time-delay between symptom-based and symptom-agnostic testing. As a result, our model estimated that in total 1304 (1,198–1,416) individuals were infected, representing 35% (32–38%) of the initial total population on the Diamond Princess. Over half of these infections had not been detected at disembarkation on February 21 st (53%, 51–56%) consisting of infected individuals who had recovered and became test negative before they were tested (37%, 34–40%), were yet to be tested (15%, 13–16%), or had recently been exposed and were not yet detectable at that point (1%, 1–3%). Nearly two-thirds of pre- and asymptomatic infections (67%, 66–68%) and 8% (6–9%) of symptomatic infections went undetected up until disembarkation (Figure 2C).

Figure 2 with 1 supplement see all

Download asset Open asset

Proportion of infections that are asymptomatic and their contribution to transmission.

(A) Prior (blue) and posterior (red) probability distribution for the proportion progressing to asymptomatic infections. (B) Prior (blue) and posterior (red) probability distribution for the relative infectiousness of asymptomatic infections. (C) Number of pre- and asymptomatic infections and symptomatic cases detected (dark red) and not detected (light red) during the outbreak. Error bars indicate 95% posterior intervals. (D) Posterior probability distribution for proportion of transmission that is from asymptomatic individuals. Dashed and dotted lines show median and interquartile range, respectively.

In contrast to the strong identifiability of the proportion of infections that were asymptomatic, the model was unable to identify the relative infectiousness of asymptomatic infections from the data, that is, a uniform prior was effectively returned (see Figure 2B). This is because the relative infectiousness of asymptomatic infections was degenerate with the overall contact rate, meaning the data were consistent with either relatively frequent contact with less infectious individuals or relatively infrequent contact with more infectious individuals (see Figure 1—figure supplement 1). Despite this, the estimated proportion of transmission due to asymptomatic infections is 69%, with a wide confidence interval (20–85%) and an interquartile range of 56–76% (Figure 2D). The reason this estimate is not effectively 0–100%, as might be expected by the unidentifiable relative infectiousness, is the combination of the strongly identified, relatively high proportion of infections that are asymptomatic and the non-linear relationship between the relative infectiousness of asymptomatics and their contribution to transmission, which quickly saturates to its maximal value (see Figure 2—figure supplement 1). The result is that only a modest relative infectiousness is required to produce a non-trivial contribution to transmission. The relative infectiousness of presymptomatics was also unidentifiable, however, in all scenarios the remaining transmission was equally distributed between the presymptomatic (14%, 1–44%) and symptomatic (17%, 11–42%) individuals. Figure 3 shows the instantaneous proportion of transmission from symptomatic (A), presymptomatic (B) and asymptomatic (C) individuals over the course of the epidemic.

Figure 3

Download asset Open asset

Instantaneous proportion of transmission from symptomatic, presymptomatic and asymptomatic individuals.

Instantaneous proportion of transmission from symptomatic (A), presymptomatic (B) and asymptomatic (C) individuals over the course of the epidemic, following the introduction of a single symptomatic individual on 20th Jan. Red lines = median, shading = 95% posterior interval. Two vertical lines show the date of the first confirmed diagnosis (left) and the start of quarantine measures (right).

Because of the non-identifiability of the relative infectiousness of asymptomatic infections we investigated marginal posterior estimates (Table 1). We find that low relative infectiousness of asymptomatic infections (0–25% compared to symptomatic individuals) would need to be compensated by a net reproduction number for individuals during their presymptomatic and symptomatic phase of 15.5–29.1.

Table 1

Model outputs by relative infectiousness of asymptomatic individuals.

Relative infectiousness expressed as proportion compared to symptomatic individuals. All values are 95% posterior ranges from model scenarios. Net reproduction number represents the typical number of infections generated by a single infected individual during their presymptomatic and symptomatic stages.

Range of relative infectiousness of asymptomatic individual	Model output
Range of relative infectiousness of asymptomatic individual	Transmission from asymptomatic individuals (%)	Net reproduction number for presymptomatic passengers	Basic reproduction number
0–1%	0–3	22.7–29.1	6.7–7.6
1–25%	7–58	15.5–25.5	7.0–8.8
25–50%	44–75	11.1–17.6	8.0–9.6
50–75%	60–82	8.7–13.6	8.7–10.2
75–99%	68–86	7.2–11.4	9.3–10.8
99–100%	72–87	6.7–10.2	9.5–10.9

Relative infectiousness expressed as proportion compared to symptomatic individuals. All values are 95% posterior ranges from model scenarios. Net reproduction number represents the typical number of infections generated by a single infected individual during their presymptomatic and symptomatic stages.

Sensitivity analyses

Without an asymptomatic state the model was unable to reconstruct the dynamics of the outbreak (Appendix 2—figures 1–3, Deviance Information Criterion (DIC) = 974 vs 329 for the primary analysis). Moreover, adjusting the relative value for mixing between crew and passengers did not have a qualitative effect on the results (Appendix 2—figures 4–11).

Biased symptom-agnostic testing had a greater impact on the results. Those testing negative having a 50% higher probability of being tested compared to the primary analysis led to a corresponding greater proportion of infections that were asymptomatic (84%, 82–87%), total number infected (2,097, 1,914-2,292) and contribution of asymptomatics to transmission (78%, 36–91%) (see Appendix 2—figures 12–15). Conversely, those testing positive having a 50% higher probability of being tested compared to the primary analysis led to a corresponding smaller proportion of infections that were asymptomatic (66%, 61–71%), total number infected (1,051, 965-1,161) and contribution of asymptomatics to transmission (59%, 9–79%) (see Appendix 2—figures 16–19).

When we assumed a fixed age-specific ratio for the proportion of infections that progress asymptomatically, the model was able to fit the data, although the number of correlated parameters was high. Overall results were similar to the main analysis, with a proportion asymptomatic of 42% (41–44%) and 89% (85–91%) for passengers and crew, respectively. The proportion of all transmission from asymptomatics was 69% (IQR = 59–74%). Relative infectiousness was again unidentifiable. See Appendix 2—figures 20–23 for details.

A longer latent period provided a poorer fit to the data (DIC = 361) (Appendix 2—figures 24–27). Adjusting the duration of the asymptomatic state to half or double the sum of the presymptomatic and symptomatic states made little qualitative difference to the results (Appendix 2—figures 28–35), although the shorter asymptomatic period was a marginally poorer fit to the data (DIC = 338). Finally, recalibrating the model assuming the 35 confirmed pre/asymptomatic cases where a test date was not available were allocated to the last feasible day (13th Feb) made no qualitative difference to our results (see Appendix 2—figures 36–39).

Discussion

Summary

We find that in this well-documented outbreak in a closed population, 74% (70–78%) of infections proceeded asymptomatically, equaling a 1:3.8 (1:3.3-1:4.4) case-to-infection ratio. The majority of asymptomatic infections remained undetected, but may have contributed substantially to ongoing transmission. While the relative infectiousness of asymptomatic infections could not be identified, low infectiousness (e.g. 0–25% compared to symptomatic individuals) would have required a very high net reproduction number for individuals during their presymptomatic and symptomatic stages of (15.5–29.1).

Interpretation

Our results are strongly informed by data, which show that when extensive symptom-agnostic testing was ramped up, substantial numbers of pre- or asymptomatic infections were identified. Given the clear suppression of transmission through quarantining, as indicated by the drop in incident symptomatic disease, this finding is most likely explained by a large proportion of undetected asymptomatic individuals.

The model and data were unable to identify a value for the relative infectiousness, although we showed how different ranges for this key parameter required specific trade-offs, as reflected in the net reproduction number for infected individuals who will develop symptomatic disease. One can argue that a net reproduction number for presymptomatic passengers at the start of the outbreak of over 20 in this population, as required if asymptomatic individuals are effectively unable to transmit (range for relative infectiousness of 0–1%) is unlikely. Such high reproductive numbers are not usually seen, exceeding for example values found for norovirus outbreaks on cruise ships (Gaythorpe et al., 2018). While SARS-CoV-2 has been shown to survive on surfaces (van Doremalen et al., 2020), this does not seem to be the primary mode of transmission. In combination with growing evidence around viral load in asymptomatic infections and their involvement in transmission chains (Lavezzo et al., 2020), as well as anecdotal evidence about transmission from asymptomatic individuals (Bai et al., 2020; Chau et al., 2020) including in closed populations (Arons et al., 2020), it is reasonable to assume that asymptomatic infections play some role in ongoing SARS-CoV transmission. In our model, asymptomatic infections were responsible for more than half of all transmission in 83% of the scenarios compatible with the data.

It is important to note that our conclusion that asymptomatic infections may have contributed substantially to ongoing transmission is critically dependent on the setting. In this case symptomatic infections were quickly identified and removed from the ship before symptom-agnostic testing began, thereby leaving asymptomatic infections to dominate transmission. Such dominance should therefore not be interpreted as a constant of nature, but instead an important consideration in settings where prompt isolation of symptomatic infections is already in place but with little to no consideration for asymptomatic infections.

Comparison to other studies

Our estimated proportion of asymptomatic infections in this outbreak is higher than previous studies, which relied on diagnosed cases only (Mizumoto et al., 2020). As we have shown, a substantial number of infections were not detected, which would explain some of the difference. Other empirical studies have found usually lower values, while some found similar ranges. While underestimation in other estimates due to low (Fontanet et al., 2020) and imbalanced participation from individuals with and without symptoms (Gudbjartsson et al., 2020) will be part of the explanation, there remains scope for unexplained variation from more complete samples (Lavezzo et al., 2020). In addition, it is possible that PCR-based testing has a lower sensitivity for asymptomatic individuals, which would further increase the proportion of infections that were asymptomatic (Chau et al., 2020).

A sensitivity analysis showed that our results were robust to age-specific probabilities of progressing to asymptomatic infections, as well as other assumptions made in the model, and driven by trends in the data.

Our estimated substantial contribution to transmission from asymptomatic infections confirms a hypothesis from Nishiura after analysing symptomatic cases occurring post-disembarkation (Nishiura, 2020). Our initial reproduction number of 9.3 (7.4–10.6) reflects the high transmission environment expected on cruise ships via increased contacts in a confined space, although lower than the value found in an earlier analysis by Rocklöv et al., 2020.

Our finding of similar contribution to transmission from presymptomatic and symptomatic individuals also matches findings by others (Ganyani et al., 2020; Liu et al., 2020a; He et al., 2020). In line with this, it is clear that having symptoms, or at least being aware of them, is not required in the transmission of SARS-CoV-2 (Ganyani et al., 2020; Liu et al., 2020a; He et al., 2020; Kimball et al., 2020; Wei et al., 2020). Although cough is often considered essential for transmission of respiratory infections (Patterson and Wood, 2019), work in tuberculosis, influenza and other coronaviruses has shown that while a cough may increase spread, it is not a requirement. Transmission from breathing, talking and sneezing is also possible, as well as transmission from contaminated surfaces (van Doremalen et al., 2020; Leung et al., 2020; Asadi et al., 2019; Williams et al., 2020).

Limitations

Additional data, in particular on the distribution of asymptomatic infections across crew and passengers, by age and shared quarantine environments would have benefited the model and potentially enabled us to estimate a range for the relative infectiousness of asymptomatic infections. A serological survey of the population, and the date and testing history of individuals who developed symptoms post-disembarkation, would also have likely informed more precise model estimates. In addition, better evidence on performance of the test used, and the associated likelihood of false-negative or false-positive results would help refine estimates. As more data become available, future model analyses of SARS-CoV-2 dynamics in closed populations should further inform the key questions we have looked to address here.

Whilst symptom-agnostic testing provided valuable insights into the pre- and asymptomatic states, such testing was not necessarily random, as was assumed in our primary analysis. Indeed, it is known that individuals were generally screened in reverse-age order (NIID, 2020). Sensitivity analyses considering biased testing still produced non-trivial results for the proportion of infections that were asymptomatic however.

Our model also assumed that the infectiousness of all transmissible states was constant over time. If instead symptomatic individuals are most infectious immediately after symptom onset (He et al., 2020), an estimate of their contribution to transmission would in principle increase. Given the likely heightened awareness of symptoms on board however, such an increase is likely to be marginal. Indeed, since our assumption of an average one-day, exponentially distributed delay between symptom onset and removal from the ship is likely to be an overestimate, a prompter removal distribution would at least in part offset such an increase.

A similar simplification was made by assuming that the probability of an individual progressing to either a presymptomatic or asymptomatic infection was independent of who infected whom. It is possible, however, that transmission from a symptomatic infection may be more likely to ultimately result in another symptomatic infection, owing to a higher infecting dose for example.

Conclusion

Asymptomatic SARS-CoV-2 infections may contribute substantially to transmission. This is essential to consider for countries when assessing the potential effectiveness of ongoing control measures to contain COVID-19.

Materials and methods

Data

Data from the Diamond Princess outbreak have been widely reported. (Mizumoto et al., 2020; Nishiura, 2020; NIID, 2020) On January 20th, the Diamond Princess cruise ship departed from Yokohama on a tour of Southeast Asia. A passenger that disembarked on January 25th in Hong Kong subsequently tested positive for SARS-CoV-2 on February 1st, reporting the date of symptom onset as January 23rd.

After arriving back in Yokohama on February 3rd, all passengers and crew were screened for symptoms, and those screening positive were then tested. The ship began quarantine on February 5th with all passengers confined to their cabins and crew undertaking essential activities only. At the start of quarantine there were 3711 individuals on board (2666 passengers and 1045 crew) with a median age of 65 (45–75 interquartile range).

Testing capacity was limited until February 11th and before then the majority of individuals tested had reported symptoms, referred to here as ‘symptom-based testing’. All individuals with a positive test at any stage were promptly removed from the ship and isolated. After February 11th, testing capacity increased and the testing of individuals irrespective of symptoms, referred to here as ‘symptom-agnostic testing’, was scaled up. In total, 314 symptomatic and 320 pre- or asymptomatic infections were reported before disembarkation was principally completed on February 21st.

We extracted the following data from Mizumoto et al., 2020; Nishiura, 2020; NIID, 2020 (see Figure 1). Firstly, the number of symptomatic cases per day (i.e. those testing positive having reported symptoms) by date of symptom onset, separately for passengers and crew. The date of symptom onset was not available for 115 cases, which we accounted for in our model structure by assuming they were distributed over time proportional to those cases with a reported date of symptom onset (see Appendix 1—table 1). Secondly, we extracted the number of pre- or asymptomatic infections identified per day (i.e. individuals testing positive having not reported symptoms) by date of test. The test date was not available for 35 pre- or asymptomatic individuals between the February 6-14th, which we assumed were distributed over time proportional to the daily number of tests performed amongst individuals not reporting symptoms. No data were available on how many individuals that tested positive in the absence of symptoms became symptomatic after disembarkation. Finally, we extracted the number of tests performed per day amongst individuals not reporting symptoms (see Appendix 1—table 2).

Model

Request a detailed protocol

We built a deterministic, compartmental model to capture transmission, disease development and the effect of interventions on board the Diamond Princess. Following exposure, after which an individual is assumed to test negative for SARS-CoV-2 for the duration of the latent phase (see Table 2), a proportion of individuals proceed asymptomatically with the remainder becoming presymptomatic. This proportion equates to a universal probability of becoming either presymptomatic or asymptomatic, independent of who infected whom. Individuals in the presymptomatic, asymptomatic or symptomatic state are assumed to test positive and have independent infectiousness, expressed relative to those with symptomatic disease.

Table 2

Model parameters and priors/values.

Parameter	Description	Prior/value	Source/Notes
$\bar{β}$	Overall contact rate (1/days)	Estimated: Uniform(0,100)
$c^{(c c)}$	Relative initial contact rate between crew/crew	Fixed: 1
$c^{(p p)}$	Relative initial contact rate between passengers/passengers	Estimated: Uniform(0,100)
$c^{(p c)}$	Relative initial contact rate between passengers/crew	Fixed relative to $c^{(p p)}$
X	Ratio: $\frac{c^{(p c)}}{c^{(p p)}}$	Fixed: 0.1	Assumed. Varied in sensitivity analyses
b₁	Percentage reduction in all initial contact rates (%)	Estimated: Uniform(0,100)
b₂	Rate of change of all contact rates (1/days)	Fixed: 10	Assumed. Transitions completed over approximately one day
$τ^{(p p)}, τ^{(p c)}$	Time of transition for contacts between passengers/passengers and passengers/crew (days)	Estimated: Uniform(0,32)	Assumed to be equal to each other
$τ^{(c c)}$	Time of transition for contacts between crew/crew (days)	Estimated: Uniform(0,32)
$θ_{p}$	Relative infectiousness of presymptomatic state	Estimated: Uniform(0,1)	Relative to symptomatic state
$θ_{a}$	Relative infectiousness of asymptomatic state	Estimated: Uniform(0,1)	Relative to symptomatic state
$χ$	Proportion of infections that proceed to asymptomatic state	Estimated: Uniform(0,1)
$\frac{1}{v}$	Latent period (days)	Fixed: 4.3	Derived from Backer et al., 2020
$\frac{1}{γ_{a}}$	Mean duration in asymptomatic state (days)	Fixed: 5.0	Assumed. Sum of mean durations in presymptomatic and symptomatic states. Varied in sensitivity analyses.
$\frac{1}{γ_{p}}$	Mean duration in presymptomatic state (days)	Fixed: 2.1	Derived from Backer et al., 2020
$\frac{1}{γ_{s}}$	Mean duration in infectious symptomatic state (days).	Fixed: 2.9	From Liu et al., 2020b Applicable only until quarantine starts on 5th Feb
$\frac{1}{μ}$	Mean delay between onset of symptomatic disease and symptom-based testing and removal (days).	Fixed: 1	Assumed. Applicable only after quarantine starts on 5th Feb.
$\frac{1}{η}$	Mean duration of test positivity following recovery (days)	Fixed: 7	From Woelfel et al., 2020
$ϕ$	Proportion of symptomatic cases with a reported onset date	Fixed: 0.661 (199/314)	From Mizumoto et al., 2020; Nishiura, 2020
$f (t)$	Rate of symptom-agnostic testing and removal (1/days)	Fixed: Calculated	From Mizumoto et al., 2020 Calculated using the number of tests administered per day amongst individuals not reporting symptoms (see Appendix 1)
$N^{(p)}$	Total number of passengers on the ship as at start of quarantine on 5th Feb	Fixed: 2666	From NIID, 2020
$N^{(c)}$	Total number of crew on the ship as at start of quarantine on 5th Feb	Fixed: 1045	From NIID, 2020

Individuals with presymptomatic infection are either detected through symptom-agnostic testing before being removed from the ship, or develop symptomatic disease. Once symptomatic disease starts, individuals can either recover undetected on the ship or, following the start of quarantine on February 5th, be detected through symptom-based testing and removed from the ship with an average delay of one day following symptom onset. We allowed for individuals to test positive after their infectious period for an average of seven days (Woelfel et al., 2020). After this, we assume they would test negative.

Individuals with asymptomatic infections either recover undetected on the ship, or are detected by symptom-agnostic testing before being removed from the ship. See Appendix 1—figure 1 for a diagram of the model.

Symptom-agnostic testing was assumed to have been random amongst those not reporting symptoms and no delay was introduced between testing and removal of those that tested positive from the ship. As such, the number of people that tested positive through symptom-agnostic testing before being removed from the ship per day was calculated using the number of tests performed per day (Figure 1F) and the proportion of individuals that were either presymptomatic, asymptomatic or recovered but continued to test positive for up to seven days, amongst all individuals on the ship not reporting symptoms. All testing was assumed to have 100% sensitivity and specificity.

Crew and passengers were modelled separately, using stratified data on the number of confirmed symptomatic cases (Figure 1A–B). We estimated the within-crew and within-passenger contact rates through calibration to the data, but assumed that the between-group contact rate was a fixed factor of 1/10th of the within-passenger rate, and explored the impact of this assumption in sensitivity analyses. We enabled the model to capture potential changes in contact behaviour between individuals by representing contact rates as sigmoid functions over time, reflecting any reductions in contact. The dates and extent of the changes were determined solely through model calibration to the data.

Model parameterisation

Request a detailed protocol

We used data from the literature to inform the natural history of COVID-19, in particular for the duration of presymptomatic and symptomatic phases (see Table 2).

Model calibration

Request a detailed protocol

The model was calibrated in a Bayesian framework. We fitted to the daily incidence of confirmed symptomatic cases with a known onset date, separately for passengers and crew, assuming a Poisson distribution in the likelihood. We simultaneously fitted to the daily number of confirmed pre- and asymptomatic infections for passengers and crew combined by using the number of tests administered per day and the prevalence of presymptomatic, asymptomatic and post-infection test-positive individuals, assuming a binomial distribution in the likelihood. We used uniform priors for the parameters to be estimated (see Table 2) and sampled the posterior of the model parameters using sequential Markov Chain Monte-Carlo (MCMC). A burn in phase during which the proposal distributions were adapted in both scale and shape to provide optimal sampling efficiency was discarded, leaving chains with 1.5 million iterations. The resultant MCMC chains were visually inspected for convergence.

Model outputs

Request a detailed protocol

Model outputs were calculated by randomly sampling 100,000 parameter values from the posterior distribution. Model trajectories were generated and compared to the data in Figure 1A–C to inspect model fit. The basic reproduction number was also calculated over time using the next-generation matrix (Diekmann et al., 2010), as a measure of ongoing transmission. We estimated the proportion of infections that become asymptomatic and the relative infectiousness of asymptomatic infections using their respective marginal posterior parameter values. Finally, the contribution of asymptomatic infections to overall transmission, as well as the net reproduction number for presymptomatic passengers at the beginning of the outbreak (i.e. the typical number of infections generated by a single presymptomatic individual) were estimated, both overall and by specific ranges of relative infectiousness. We report the median and 95% equal-tailed posterior intervals throughout.

Sensitivity analyses

Request a detailed protocol

We recalibrated the model for a number of alternative scenarios to assess model sensitivity. Firstly, we assessed the impact of removing the asymptomatic phase (i.e. 100% of infections progressed to symptomatic disease). Secondly, we explored the impact of assuming different values for the relative mixing between crew and passengers as well as shorter and longer durations of asymptomatic infection. Thirdly, we considered the impact of biased symptom-agnostic testing. Specifically, we first assumed that those that would test positive were 50% more likely to be tested, before then assuming that those that would test negative were 50% more likely to be tested, both compared to purely random testing as per the primary analysis. We also explored the impact of assuming a different proportion of asymptomatic infections for crew and passengers based on their distinct median ages (36 years for crew, 69 years for passengers), using a fixed ratio for the two proportions taken from the results of a model fitted to epidemic data in six countries by Davies et al., 2020. In addition, we explored a longer latent period given the relatively high age in our population (Jiang et al., 2020). Finally, we recalibrated the model assuming the 35 confirmed pre/asymptomatic cases where a test date was not available were allocated to the last feasible day (13th Feb) instead of proportionate to the overall number of tests over the period February 6-14th(see Appendix 2 for further details).

All analyses were conducted using R version 3.5.0 (R Development Core Team, 2014). Bayesian calibration was performed in LibBi (Murray, 2013) using RBi (Funk, 2019) as an interface. Replication data and analyses scripts are available on GitHub at https://github.com/thimotei/covid19_asymptomatic_trans (Emery et al., 2020; copy archived at https://github.com/elifesciences-publications/covid19_asymptomatic_trans).

Role of funding source

Request a detailed protocol

The funder of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report. The corresponding author had full access to all the data in the study and had final responsibility for the decision to submit for publication.

Appendix 1

Extended methods

Data

Data for confirmed symptomatic cases were extracted from Nishiura, 2020. Appendix 1—table 1 shows n = 199 confirmed symptomatic cases by date of symptom onset for passengers and crew separately. Symptom onset dates were unavailable for a further n = 115 confirmed symptomatic cases. These were accounted for in the model structure (see Appendix 1—figure 1) by assuming they were distributed over time proportional to those cases with a reported date of symptom onset. The data itself were not augmented.

Appendix 1—table 1

Confirmed symptomatic cases (n = 199) by date of symptom onset for passengers and crew separately, extracted from Nishiura, 2020.

A further n = 115 confirmed symptomatic cases without symptom onset dates are not included in the table.

Date of symptom onset	Confirmed symptomatic cases
Date of symptom onset	Passengers	Crew	Total
20-Jan	2	0	2
21-Jan	0	0	0
22-Jan	0	0	0
23-Jan	1	0	1
24-Jan	0	0	0
25-Jan	0	0	0
26-Jan	0	0	0
27-Jan	0	0	0
28-Jan	0	0	0
29-Jan	1	0	1
30-Jan	1	0	1
31-Jan	0	0	0
01-Feb	4	0	4
02-Feb	4	0	4
03-Feb	4	0	4
04-Feb	6	0	6
05-Feb	12	0	12
06-Feb	15	2	17
07-Feb	29	2	31
08-Feb	17	2	19
09-Feb	19	5	24
10-Feb	7	3	10
11-Feb	11	8	19
12-Feb	5	7	12
13-Feb	9	8	17
14-Feb	2	5	7
15-Feb	1	3	4
16-Feb	0	3	3
17-Feb	0	1	1
18-Feb	0	0	0
19-Feb	0	0	0
20-Feb	0	0	0
Total	150	49	199

Data for confirmed pre/asymptomatic cases and symptom-agnostic testing were extracted from Mizumoto et al., 2020. Appendix 1—table 2 shows n = 2749 symptom-agnostic tests and n = 320 confirmed pre/asymptomatic cases by date of test for passengers and crew combined, since stratification by passenger/crew was unavailable. The number of symptom-agnostic tests was inferred from the total number of tests each day, minus the number of positive results in individuals reporting symptoms in Mizumoto et al., 2020. Test dates were not available for n = 35 confirmed pre/asymptomatic cases between 5th-14th Feb. These were distributed proportional to the total number of tests (symptom-based and symptom-agnostic) on those days. An alternative scenario where all n = 35 confirmed pre/asymptomatic cases are allocated to the last possible day (13th Feb) is explored in sensitivity analyses.

Appendix 1—table 2

Confirmed pre/asymptomatic cases (n = 320) and symptom-agnostic tests (n = 2749) by date of test for passengers and crew combined, extracted from Mizumoto et al., 2020.

⁺Test dates were not available for n = 35 confirmed pre/asymptomatic cases between 5th-14th Feb. These were distributed proportional to the total number of tests (symptom-based and symptom-agnostic) on those days.

Date of test	Number of symptom agnostic tests	Number of confirmed pre/asymptomatic cases
20-Jan	0	0
21-Jan	0	0
22-Jan	0	0
23-Jan	0	0
24-Jan	0	0
25-Jan	0	0
26-Jan	0	0
27-Jan	0	0
28-Jan	0	0
29-Jan	0	0
30-Jan	0	0
31-Jan	0	0
01-Feb	0	0
02-Feb	0	0
03-Feb	0	0
04-Feb	0	0
05-Feb⁺	23	2
06-Feb⁺	64	3
07-Feb⁺	138	8
08-Feb⁺	3	0
09-Feb⁺	54	3
10-Feb⁺	43	5
11-Feb⁺	0	0
12-Feb⁺	17	3
13-Feb⁺	188	11
14-Feb⁺	0	0
15-Feb	188	38
16-Feb	257	38
17-Feb	475	70
18-Feb	658	65
19-Feb	596	68
20-Feb	45	6
Total	2749	320

Model

The model described in the main article is shown in detail in Appendix 1—figure 1, where passengers ( $i = p$ ) and crew ( $i = c$ ) are modelled separately and the annotated parameters are described in Table 2 in the main article.

Appendix 1—figure 1

Download asset Open asset

Model diagram for the outbreak onboard the Diamond Princess cruise ship described in the main paper.

The annotated transition parameters are defined in Table 2 in the main article and detailed further, below. The model is stratified by i = passengers or crew. The asymptomatic, presymptomatic and symptomatic states are all assumed to be infectious and individuals would test positive during symptom-based or symptom-agnostic testing. Individuals that recover are also assumed to test positive for an average of 1 week after they are no longer infectious.

The force of infection is given by

λ^{(i)} (t) = \sum_{j = p, c} β^{(i j)} (t) \frac{(θ_{a} I_{a}^{(j)} + θ_{p} I_{p}^{(j)} + I_{s k}^{(j)} + I_{s u}^{(j)})}{N^{(j)}}

where the time dependent contact parameters are described by sigmoid functions

β^{(i j)} (t) = \bar{β} c^{(i j)} (1 - \frac{b_{1}}{1 + e^{- b_{2} (t - τ^{(i j)})}})

And $τ^{(p p)} = τ^{(p c)} = τ^{(c p)}$ (i.e. contact between passengers/passengers and passengers/crew is reduced at the same time, which can differ from contact between crew/crew).

The transition from exposed to presymptomatic or asymptomatic is modelled as an erlang distribution using two compartments (i.e. a shape parameter $k = 2$ ), each with a mean duration of $1 / 2 v$ .

The rate of symptom agnostic testing and removal of individuals not reporting symptoms is given by the total number of symptom agnostic tests administered per day divided by the total number of individuals not presenting symptoms being tested on that day

f (t) = \frac{N^{t e s t s}}{(S + E + I_{a} + I_{p} + T + C)}

Where $N^{t e s t s}$ is taken from the data in Appendix 1—table 2 and variables without indices represent the totals among passengers and crew (e.g. $S = S^{(p)} + S^{(c)}$ )

To reflect heightened symptom awareness following quarantine, the transition rate from symptomatic infection to recovered on the ship is constant before quarantine and zero afterwards, whilst the rate of removal of individuals reporting symptoms is zero before quarantine and a constant afterwards

\begin{array}{ll} γ_{s} (t) = γ_{s}, μ (t) = 0 & f o r t < 5 t h F e b \\ γ_{s} (t) = 0, μ (t) = μ & f o r t \geq 5 t h F e b \end{array}

All other model transitions are exponentially distributed.

The model is initialised with a single symptomatic passenger with a known onset date on 20th Jan, with all other individuals susceptible

I_{s k}^{(p)} (0) = 1, S^{(p)} = N^{(p)} - 1, S^{(c)} = N^{(c)}

Model calibration

The model was calibrated in a Bayesian framework to fit to the two sets of empirical observations from the ship (Appendix 1—tables 1–2). We used a Poisson likelihood for the incident symptomatic cases with a known onset date for crew and passengers separately. We used a Binomial likelihood for the number of confirmed pre- and asymptomatic infections for passengers and crew combined, using the number of tests administered per day and the prevalence of presymptomatic, asymptomatic and post-infection test-positive individuals. The complete likelihood is given by

L = (\prod_{k = 1}^{K} P o i s s o n (Z_{k}^{(c)} | m e a n = z_{k}^{(c)})) (\prod_{k = 1}^{K} P o i s s o n (Z_{k}^{(p)} | m e a n = z_{k}^{(p)})) (\prod_{k = 1}^{K} B i n o m (Y_{k} | N_{k}^{t e s t s}, m e a n = y_{k}))

where $Z_{k}^{(i)}$ is the observed incidence of symptomatic cases with a known date of onset on day k for passengers p or crew $c, z_{k}^{(i)}$ is the model predicted incidence, $Y_{k}$ is the observed prevalence of presymptomatic, asymptomatic and post-infection test-positive individuals (passengers and crew combined) amongst $N_{k}^{t e s t s}$ symptom-agnostic tests, and $y_{k}$ is the model predicted prevalence

y (t) = \frac{I_{a} + I_{p} + T}{S + E + I_{a} + I_{p} + T + C}

We used uniform priors for the parameters to be estimated (see Table 2 in the main article).

Model outputs

The basic reproduction number as a function of time $R_{0} (t)$ was calculated by first constructing the next generation matrix (NGM) at each time point using the relevant Jacobian matrices (Diekmann et al., 2010). The basic reproduction number is then given by the absolute value of the dominant eigenvalue of the NGM.

The net reproduction number for a presymptomatic infection (i.e. the typical number of secondary infections caused by a single presymptomatic individual throughout both their presymptomatic and symptomatic periods) at the beginning of the outbreak is given by the respective entry in the NGM evaluated at $t = 0$ .

The instantaneous proportion of transmission from either symptomatic, presymptomatic or asymptomatic individuals was calculated by dividing the number of infections generated by the respective infected state in the previous timestep by the total number of new infections in the previous timestep. The overall proportion of transmission from asymptomatic individuals was given by the cumulative number of infections caused by asymptomatics divided by the cumulative number of total infections, evaluated at the end of the outbreak.

Appendix 2

Sensitivity analyses results

1. Presymptomatic infection only

Assumes the proportion of infections that are asymptomatic and their relative infectiousness are zero ( $χ = 0$ and $θ_{a} = 0$ ). The latent period $1 / v$ is estimated with a uniform prior between 1 and 21 days.

Appendix 2—figure 1

Download asset Open asset

Appendix 2—figure 2

Download asset Open asset

Parameter correlation plot containing parameter values from 10,000 samples of the joint posterior distribution found during MCMC model calibration.

See Table 2 in the main article for parameter definitions and descriptions.

Appendix 2—figure 3

Download asset Open asset

Parameter trace plot showing all 1.5 million samples from the MCMC model calibration sequentially.

See Table 2 in the main article for parameter definitions and descriptions.

2. Relative passenger-crew contact rate: X = 0.02

Assumes the contact rate between passengers and crew is 1/50th of contacts between passengers and passengers $(c^{(p c)} / c^{(p p)} = 0.02)$ .

Appendix 2—figure 4

Download asset Open asset

Appendix 2—figure 5

Download asset Open asset

Appendix 2—figure 6

Download asset Open asset

Appendix 2—figure 7

Download asset Open asset

3. Relative passenger-crew contact rate: X = 0.5

Assumes the contact rates between passengers and crew is half that of between passengers and passengers $(c^{(p c)} / c^{(p p)} = 0.5)$ .

Appendix 2—figure 8

Download asset Open asset

Appendix 2—figure 9

Download asset Open asset

Appendix 2—figure 10

Download asset Open asset

Appendix 2—figure 11

Download asset Open asset

4. Biased symptom-agnostic testing: test-negative individuals more likely to be tested

Assumes individuals that would test negative during symptom-agnostic testing are 50% more likely to be tested compared to the primary analysis, where testing is random.

Appendix 2—figure 12

Download asset Open asset

Appendix 2—figure 13

Download asset Open asset

Appendix 2—figure 14

Download asset Open asset

Appendix 2—figure 15

Download asset Open asset

5. Biased symptom-agnostic testing: test-positive individuals more likely to be tested

Assumes individuals that would test positive during symptom-agnostic testing are 50% more likely to be tested compared to the primary analysis, where testing is random.

Appendix 2—figure 16

Download asset Open asset

Appendix 2—figure 17

Download asset Open asset

Appendix 2—figure 18

Download asset Open asset

Appendix 2—figure 19

Download asset Open asset

6. Age dependent proportion asymptomatic

Assumes separate asymptomatic proportions for crew $(χ^{(c)})$ and passengers $(χ^{(p)})$ to reflect their different age demographics (median ages of 36 and 69 respectively), compared to a single asymptomatic proportion in the primary analysis. The ratio $χ^{(p)} / χ^{(c)}$ was fixed at 0.48 using the results for asymptomatic proportion by age from a model fitted to epidemic data in six countries by Davies et al., 2020.

Appendix 2—figure 20

Download asset Open asset

Appendix 2—figure 21

Download asset Open asset

Appendix 2—figure 22

Download asset Open asset

Appendix 2—figure 23

Download asset Open asset

7. Duration of latent period: 8.8 days

Assumes the average duration for the latent period is $1 / v = 8.8$ days (Jiang et al., 2020), compared to 4.3 in the primary analysis.

Appendix 2—figure 24

Download asset Open asset

Appendix 2—figure 25

Download asset Open asset

Appendix 2—figure 26

Download asset Open asset

Appendix 2—figure 27

Download asset Open asset

8. Duration of asymptomatic infection: 2.5 days

Assumes the average duration of asymptomatic infection is $1 / γ_{a} = 2.5$ days, compared to 5 days in the primary analysis.

Appendix 2—figure 28

Download asset Open asset

Appendix 2—figure 29

Download asset Open asset

Appendix 2—figure 30

Download asset Open asset

Appendix 2—figure 31

Download asset Open asset

9. Duration of asymptomatic infection: 10 days

Assumes the average duration of asymptomatic infection is $1 / γ_{a} = 10$ days, compared to 5 days in the primary analysis.

Appendix 2—figure 32

Download asset Open asset

Appendix 2—figure 33

Download asset Open asset

Appendix 2—figure 34

Download asset Open asset

Appendix 2—figure 35

Download asset Open asset

10. Alternative distribution of n = 35 confirmed pre/asymptomatic cases

Assumes that n = 35 confirmed pre/asymptomatic cases without a test are apportioned to the last possible day (13th Feb), compared to proportional to the total number of tests administered over 6th-13th Feb in the primary analysis.

Appendix 2—figure 36

Download asset Open asset

Appendix 2—figure 37

Download asset Open asset

Appendix 2—figure 38

Download asset Open asset

Appendix 2—figure 39

Download asset Open asset

Data availability

All data analysed during this study are included in the manuscript and supporting files. Model code and data are available through github at https://github.com/thimotei/covid19_asymptomatic_trans (copy archived athttps://github.com/elifesciences-publications/covid19_asymptomatic_trans).

References

1. Arons MM
2. Hatfield KM
3. Reddy SC
4. Kimball A
5. James A
6. Jacobs JR
7. Taylor J
8. Spicer K
9. Bardossy AC
10. Oakley LP
11. Tanwar S
12. Dyal JW
13. Harney J
14. Chisty Z
15. Bell JM
16. Methner M
17. Paul P
18. Carlson CM
19. McLaughlin HP
20. Thornburg N
21. Tong S
22. Tamin A
23. Tao Y
24. Uehara A
25. Harcourt J
26. Clark S
27. Brostrom-Smith C
28. Page LC
29. Kay M
30. Lewis J
31. Montgomery P
32. Stone ND
33. Clark TA
34. Honein MA
35. Duchin JS
36. Jernigan JA
37. Public Health–Seattle and King County and CDC COVID-19 Investigation Team
(2020) Presymptomatic SARS-CoV-2 infections and transmission in a skilled nursing facility
New England Journal of Medicine 382:2081–2090.

https://doi.org/10.1056/NEJMoa2008457
- PubMed
- Google Scholar
1. Asadi S
2. Wexler AS
3. Cappa CD
4. Barreda S
5. Bouvier NM
6. Ristenpart WD
(2019) Aerosol emission and superemission during human speech increase with voice loudness
Scientific Reports 9:2348.

https://doi.org/10.1038/s41598-019-38808-z
- PubMed
- Google Scholar
(2020) Incubation period of 2019 novel coronavirus (2019-nCoV) infections among travellers from Wuhan, China, 20–28 January 2020
Eurosurveillance 25:2000062.

https://doi.org/10.2807/1560-7917.ES.2020.25.5.2000062
- Google Scholar
1. Bai Y
2. Yao L
3. Wei T
4. Tian F
5. Jin D-Y
6. Chen L
7. Wang M
(2020) Presumed asymptomatic carrier transmission of COVID-19
Jama 323:1406–1407.

https://doi.org/10.1001/jama.2020.2565
- Google Scholar
Website
1. Bendavid E
2. Mulaney B
3. Sood N
4. Shah S
5. Ling E
6. Bromley-Dulfano R
(2020) COVID-19 antibody seroprevalence in santa clara county, California
Epidemiology. Accessed April 18, 2020.

http://medrxiv.org/lookup/doi/10.1101/2020.04.14.20062463
Preprint
1. Chau NVV
2. Lam VT
3. Dung NT
4. Yen LM
5. Minh NNQ
6. Hung LM
(2020) The natural history and transmission potential of asymptomatic SARS-CoV-2 infection
medRxiv.

https://doi.org/10.1101/2020.04.27.20082347
- Google Scholar
1. Chen Y
2. Wang A
3. Yi B
4. Keqin D
5. Wang H
6. Wang J
(2020) The epidemiological characteristics of infection in close contacts of COVID-19 in Ningbo City
Chin J Epidemiol 41:668–672.

https://doi.org/10.3760/cma.j.cn112338-20200304-00251
- Google Scholar
Preprint
(2020) Age-dependent effects in the transmission and control of COVID-19 epidemics
medRxiv.

https://doi.org/10.1101/2020.03.24.20043018
- Google Scholar
1. Dean NE
2. Halloran ME
3. Yang Y
4. Longini IM
(2016) Transmissibility and pathogenicity of ebola virus: a systematic review and Meta-analysis of household secondary attack rate and asymptomatic infection
Clinical Infectious Diseases 62:1277–1286.

https://doi.org/10.1093/cid/ciw114
- PubMed
- Google Scholar
(2010) The construction of next-generation matrices for compartmental epidemic models
Journal of the Royal Society Interface 7:873–885.

https://doi.org/10.1098/rsif.2009.0386
- Google Scholar
Software
(2020) The contribution of asymptomatic SARS-CoV-2 infections to transmission - a model-based analysis of the Diamond Princess outbreak, version 116b196
GitHub.

https://github.com/thimotei/covid19_asymptomatic_trans
(2018) Tuberculosis transmission during the subclinical period: could unrelated cough play a part?
The Lancet Respiratory Medicine 6:244–246.

https://doi.org/10.1016/S2213-2600(18)30105-X
- PubMed
- Google Scholar
Preprint
1. Fontanet A
2. Tondeur L
3. Madec Y
4. Grant R
5. Besombes C
6. Jolly N
(2020) Cluster of COVID-19 in northern France: a retrospective closed cohort study
medRxiv.

https://doi.org/10.1101/2020.04.18.20071134
- Google Scholar
Software
1. Funk S
(2019) Sbfnk/RBi, version 43e77da
Sbfnk/RBi.

https://github.com/sbfnk/RBi
Preprint
1. Ganyani T
2. Kremer C
3. Chen D
4. Torneri A
5. Faes C
6. Wallinga J
(2020) Estimating the generation interval for COVID-19 based on symptom onset data
medRxiv.

https://doi.org/10.1101/2020.03.05.20031815
- Google Scholar
(2018) Norovirus transmission dynamics: a modelling review
Epidemiology and Infection 146:147–158.

https://doi.org/10.1017/S0950268817002692
- PubMed
- Google Scholar
Preprint
1. Golding N
2. Russell TW
3. Abbott S
4. Hellewell J
5. Pearson CAB
6. van Zandvoort K
7. Jarvis CI
8. Gibbs H
9. Liu Y
10. Eggo RM
11. Edmunds JW
12. Kucharski A
(2020) Reconstructing the global dynamics of under-ascertained COVID-19 cases and infections
medRxiv.

https://doi.org/10.1101/2020.07.07.20148460
- Google Scholar
(2020) Spread of SARS-CoV-2 in the icelandic population
New England Journal of Medicine 382:2302–2315.

https://doi.org/10.1056/NEJMoa2006100
- PubMed
- Google Scholar
1. He X
2. Lau EHY
3. Wu P
4. Deng X
5. Wang J
6. Hao X
7. Lau YC
8. Wong JY
9. Guan Y
10. Tan X
11. Mo X
12. Chen Y
13. Liao B
14. Chen W
15. Hu F
16. Zhang Q
17. Zhong M
18. Wu Y
19. Zhao L
20. Zhang F
21. Cowling BJ
22. Li F
23. Leung GM
(2020) Temporal dynamics in viral shedding and transmissibility of COVID-19
Nature Medicine 26:672–675.

https://doi.org/10.1038/s41591-020-0869-5
- PubMed
- Google Scholar
Preprint
1. Jiang AB
2. Lieu R
3. Quenby S
(2020) Significantly longer Covid-19 incubation times for the elderly from a case study of 136 patients throughout China
medRxiv.

https://doi.org/10.1101/2020.04.14.20065896
- Google Scholar
1. Kimball A
2. Hatfield KM
3. Arons M
4. James A
5. Taylor J
6. Spicer K
7. Bardossy AC
8. Oakley LP
9. Tanwar S
10. Chisty Z
11. Bell JM
12. Methner M
13. Harney J
14. Jacobs JR
15. Carlson CM
16. McLaughlin HP
17. Stone N
18. Clark S
19. Brostrom-Smith C
20. Page LC
21. Kay M
22. Lewis J
23. Russell D
24. Hiatt B
25. Gant J
26. Duchin JS
27. Clark TA
28. Honein MA
29. Reddy SC
30. Jernigan JA
31. Public Health – Seattle & King County
32. CDC COVID-19 Investigation Team
(2020) Asymptomatic and presymptomatic SARS-CoV-2 infections in residents of a Long-Term care skilled nursing facility - King county, Washington, march 2020
MMWR. Morbidity and Mortality Weekly Report 69:377–381.

https://doi.org/10.15585/mmwr.mm6913e1
- PubMed
- Google Scholar
Website
(2020) Suppression of COVID-19 outbreak in the municipality of vo, Italy
Epidemiology. Accessed April 19, 2020.

https://www.medrxiv.org/content/10.1101/2020.04.17.20053157v1
1. Leung NHL
2. Chu DKW
3. Shiu EYC
4. Chan K-H
5. McDevitt JJ
6. Hau BJP
7. Yen H-L
8. Li Y
9. Ip DKM
10. Peiris JSM
11. Seto W-H
12. Leung GM
13. Milton DK
14. Cowling BJ
(2020) Respiratory virus shedding in exhaled breath and efficacy of face masks
Nature Medicine 26:676–680.

https://doi.org/10.1038/s41591-020-0843-2
- Google Scholar
1. Li R
2. Pei S
3. Chen B
4. Song Y
5. Zhang T
6. Yang W
7. Shaman J
(2020) Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2)
Science 368:489–493.

https://doi.org/10.1126/science.abb3221
- Google Scholar
(2020a) The contribution of pre-symptomatic infection to the transmission dynamics of COVID-2019
Wellcome Open Research 5:58.

https://doi.org/10.12688/wellcomeopenres.15788.1
- PubMed
- Google Scholar
Preprint
1. Liu T
2. Hu J
3. Kang M
4. Lin L
5. Zhong H
6. Xiao J
(2020b) Transmission dynamics of 2019 novel coronavirus (2019-nCoV)
bioRxiv.

https://doi.org/10.1101/2020.01.25.919787
- Google Scholar
(2020) Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the diamond princess cruise ship, Yokohama, Japan, 2020
Eurosurveillance 25:180.

https://doi.org/10.2807/1560-7917.ES.2020.25.10.2000180
- Google Scholar
Preprint
1. Murray LM
(2013) Bayesian State-Space modelling on High-Performance hardware using LibBi
arXiv.

https://arxiv.org/abs/1306.3277
- Google Scholar
Website
1. NIID
(2020) Field briefing: diamond princess COVID-19 cases
Accessed March 12, 2020.

https://www.niid.go.jp/niid/en/2019-ncov-e/9417-covid-dp-fe-02.html
Website
1. Nishiura H
2. Kobayashi T
3. Miyama T
4. Suzuki A
5. Jung S
6. Hayashi K
(2020) Estimation of the asymptomatic ratio of novel coronavirus infections (COVID-19) [Internet]
Epidemiology. Accessed March 17, 2020.

http://medrxiv.org/lookup/doi/10.1101/2020.02.03.20020248
1. Nishiura H
(2020) Backcalculating the incidence of infection with COVID-19 on the diamond princess
Journal of Clinical Medicine 9:657.

https://doi.org/10.3390/jcm9030657
- Google Scholar
1. Patterson B
2. Wood R
(2019) Is cough really necessary for TB transmission?
Tuberculosis 117:31–35.

https://doi.org/10.1016/j.tube.2019.05.003
- PubMed
- Google Scholar
Software
1. R Development Core Team
(2014) R: A Language and Environment for Statistical Computing
R Foundation for Statistical Computing, Vienna, Austria.

http://www.R-project.org/
(2020) COVID-19 outbreak on the diamond princess cruise ship: estimating the epidemic potential and effectiveness of public health countermeasures
Journal of Travel Medicine 27:taaa030.

https://doi.org/10.1093/jtm/taaa030
- PubMed
- Google Scholar
1. Rothe C
2. Schunk M
3. Sothmann P
4. Bretzel G
5. Froeschl G
6. Wallrauch C
7. Zimmer T
8. Thiel V
9. Janke C
10. Guggemos W
11. Seilmaier M
12. Drosten C
13. Vollmar P
14. Zwirglmaier K
15. Zange S
16. Wölfel R
17. Hoelscher M
(2020) Transmission of 2019-nCoV infection from an asymptomatic contact in Germany
New England Journal of Medicine 382:970–971.

https://doi.org/10.1056/NEJMc2001468
- PubMed
- Google Scholar
Website
1. Salomon JA
2. COVID-19 Statistics, Policy modeling, and Epidemiology Collective
(2020) Defining high-value information for COVID-19 decision-making [Internet]
Health Policy. Accessed April 20, 2020.

http://medrxiv.org/lookup/doi/10.1101/2020.04.06.20052506
1. Slater HC
2. Ross A
3. Felger I
4. Hofmann NE
5. Robinson L
6. Cook J
7. Gonçalves BP
8. Björkman A
9. Ouedraogo AL
10. Morris U
11. Msellem M
12. Koepfli C
13. Mueller I
14. Tadesse F
15. Gadisa E
16. Das S
17. Domingo G
18. Kapulu M
19. Midega J
20. Owusu-Agyei S
21. Nabet C
22. Piarroux R
23. Doumbo O
24. Doumbo SN
25. Koram K
26. Lucchi N
27. Udhayakumar V
28. Mosha J
29. Tiono A
30. Chandramohan D
31. Gosling R
32. Mwingira F
33. Sauerwein R
34. Paul R
35. Riley EM
36. White NJ
37. Nosten F
38. Imwong M
39. Bousema T
40. Drakeley C
41. Okell LC
(2019) The temporal dynamics and infectiousness of subpatent Plasmodium falciparum infections in relation to parasite density
Nature Communications 10:1433.

https://doi.org/10.1038/s41467-019-09441-1
- PubMed
- Google Scholar
(2020) Aerosol and surface stability of SARS-CoV-2 as compared with SARS-CoV-1
New England Journal of Medicine 382:1564–1567.

https://doi.org/10.1056/NEJMc2004973
- PubMed
- Google Scholar
Website
1. Wang C
2. Liu L
3. Hao X
4. Guo H
5. Wang Q
6. Huang J
(2020) Evolving epidemiology and impact of Non-pharmaceutical interventions on the outbreak of coronavirus disease 2019 in Wuhan, China [Internet]
Epidemiology. Accessed April 10, 2020.

http://medrxiv.org/lookup/doi/10.1101/2020.03.03.20030593
1. Wei WE
2. Li Z
3. Chiew CJ
4. Yong SE
5. Toh MP
6. Lee VJ
(2020) Presymptomatic transmission of SARS-CoV-2 - Singapore, January 23-March 16, 2020
MMWR. Morbidity and Mortality Weekly Report 69:411–415.

https://doi.org/10.15585/mmwr.mm6914e1
- PubMed
- Google Scholar
1. Williams CM
2. Abdulwhhab M
3. Birring SS
4. De Kock E
5. Garton NJ
6. Townsend E
7. Pareek M
8. Al-Taie A
9. Pan J
10. Ganatra R
11. Stoltz AC
12. Haldar P
13. Barer MR
(2020) Exhaled Mycobacterium tuberculosis output and detection of subclinical disease by face-mask sampling: prospective observational studies
The Lancet Infectious Diseases 20:607–617.

https://doi.org/10.1016/S1473-3099(19)30707-8
- PubMed
- Google Scholar
1. Woelfel R
2. Corman VM
3. Guggemos W
4. Seilmaier M
5. Zange S
6. Mueller MA
(2020) Clinical presentation and virological assessment of hospitalized cases of coronavirus disease 2019 in a travel-associated transmission cluster
Infectious Diseases 1:502.

https://doi.org/10.1101/2020.03.05.20030502
- Google Scholar

Article and author information

Author details

Jon C Emery

Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, United Kingdom

Contribution
Conceptualization, Data curation, Software, Formal analysis, Methodology, Writing - original draft, Writing - review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-6644-7604
Timothy W Russell

Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, United Kingdom

Contribution
Software, Formal analysis, Methodology, Writing - review and editing

Competing interests
No competing interests declared
Yang Liu

Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, United Kingdom

Contribution
Data curation, Software, Methodology, Writing - review and editing

Competing interests
No competing interests declared
Joel Hellewell

Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, United Kingdom

Contribution
Software, Formal analysis, Methodology, Writing - review and editing

Competing interests
No competing interests declared
Carl AB Pearson

Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, United Kingdom

Contribution
Software, Methodology, Writing - review and editing

Competing interests
No competing interests declared
CMMID COVID-19 Working Group

Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, United Kingdom

Contribution
Data curation, Writing - review and editing

Competing interests
No competing interests declared
1. Katherine E Atkins
2. Petra Klepac
3. Akira Endo
4. Christopher I Jarvis
5. Nicholas G Davies
6. Eleanor M Rees
7. Sophie R Meakin
8. Alicia Rosello
9. Kevin van Zandvoort
10. James D Munday
11. W John Edmunds
12. Thibaut Jombart
13. Megan Auzenbergs
14. Emily S Nightingale
15. Mark Jit
16. Sam Abbott
17. David Simons
18. Nikos I Bosse
19. Quentin J Leclerc
20. Simon R Procter
21. C Julian Villabona-Arenas
22. Damien C Tully
23. Arminder K Deol
24. Fiona Yueqian Sun
25. Stéphane Hué
26. Anna M Foss
27. Kiesha Prem
28. Graham Medley
29. Amy Gimma
30. Rachel Lowe
31. Samuel Clifford
32. Matthew Quaife
33. Charlie Diamond
34. Hamish P Gibbs
35. Billy J Quilty
36. Kathleen OReilly
Gwenan M Knight

Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, United Kingdom

Contribution
Methodology, Writing - review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-7263-9896
Rosalind M Eggo

Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, United Kingdom

Contribution
Methodology, Writing - review and editing

Competing interests
No competing interests declared
Adam J Kucharski

Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, United Kingdom

Contribution
Methodology, Writing - review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-8814-9421
Sebastian Funk

Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, United Kingdom

Contribution
Software, Methodology, Writing - review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-2842-3406
Stefan Flasche

Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, United Kingdom

Contribution
Conceptualization, Methodology, Writing - review and editing

Competing interests
No competing interests declared
Rein MGJ Houben

Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, United Kingdom

Contribution
Conceptualization, Methodology, Writing - original draft, Writing - review and editing

For correspondence
rein.houben@lshtm.ac.uk

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-4132-7467

Funding

European Research Council (Action Number 757699)

Jon C Emery
Rein MGJ Houben

Bill and Melinda Gates Foundation (INV-003174)

Yang Liu

Bill and Melinda Gates Foundation (NTD Modelling Consortium OPP1184344)

Carl AB Pearson

Department for International Development, UK Government (Epidemic Preparedness Coronavirus research programme 221303/Z/20/Z)

Carl AB Pearson

Horizon 2020 (project EpiPose (101003688))

Yang Liu

HDR UK (MR/S003975/1)

Rosalind M Eggo

National Institute for Health Research (16/137/109)

Yang Liu

Medical Research Council (MC_PC 19065)

Rosalind M Eggo

Medical Research Council (MR/P014658/1)

Gwenan M Knight

Wellcome (206250/Z/17/Z)

Timothy W Russell
Adam J Kucharski

Wellcome (208812/Z/17/Z)

Stefan Flasche

Wellcome (210758/Z/18/Z)

Joel Hellewell
Sebastian Funk

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

The authors would like to thank Dr Taichi Kidani for help in translating primary data sources.

The following funding sources are acknowledged as providing funding for the named authors. This research was partly funded by the Bill and Melinda Gates Foundation (INV-003174: YL; NTD Modelling Consortium OPP1184344: CABP). DFID/Wellcome Trust (Epidemic Preparedness Coronavirus research programme 221303/Z/20/Z: CABP). ERC Starting Grant (#757699: JCE, RMGJH). This project has received funding from the European Union’s Horizon 2020 research and innovation programme - project EpiPose (101003688: YL). HDR UK (MR/S003975/1: RME). This research was partly funded by the National Institute for Health Research (NIHR) using UK aid from the UK Government to support global health research. The views expressed in this publication are those of the author(s) and not necessarily those of the NIHR or the UK Department of Health and Social Care (16/137/109: YL). UK MRC (MC_PC 19065: RME; MR/P014658/1: GMK). Wellcome Trust (206250/Z/17/Z: AJK, TWR; 208812/Z/17/Z: SFlasche; 210758/Z/18/Z: JH, SFunk).

The following funding sources are acknowledged as providing funding for the CMMID COVID-19 working group. Alan Turing Institute (AE). BBSRC LIDP (BB/M009513/1: DS). This research was partly funded by the Bill and Melinda Gates Foundation (INV-003174: KP, MJ; NTD Modelling Consortium OPP1184344: GM; OPP1180644: SRP; OPP1183986: ESN; OPP1191821: KO’R, MA). DFID/Wellcome Trust (Epidemic Preparedness Coronavirus research programme 221303/Z/20/Z: KvZ). Elrha R2HC/UK DFID/Wellcome Trust/This research was partly funded by the National Institute for Health Research (NIHR) using UK aid from the UK Government to support global health research. The views expressed in this publication are those of the author(s) and not necessarily those of the NIHR or the UK Department of Health and Social Care (KvZ). ERC Starting Grant (#757699: MQ). This project has received funding from the European Union’s Horizon 2020 research and innovation programme - project EpiPose (101003688: KP, MJ, PK, WJE). This research was partly funded by the Global Challenges Research Fund (GCRF) project ‘RECAP’ managed through RCUK and ESRC (ES/P010873/1: AG, CIJ, TJ). Nakajima Foundation (AE). NIHR (16/137/109: BJQ, CD, FYS, MJ; Health Protection Research Unit for Modelling Methodology HPRU-2012–10096: NGD, TJ; PR-OD-1017–20002: AR). Royal Society (Dorothy Hodgkin Fellowship: RL; RP\EA\180004: PK). UK DHSC/UK Aid/NIHR (ITCRZ 03010: HPG). UK MRC (LID DTP MR/N013638/1: EMR, QJL). Authors of this research receive funding from UK Public Health Rapid Support Team funded by the United Kingdom Department of Health and Social Care (TJ). Wellcome Trust (208812/Z/17/Z: SC; 210758/Z/18/Z: JDM, NIB, SA, SRM). No funding (AKD, AMF, DCT, SH).

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

6,514

views
639

downloads
58

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Jon C Emery
Timothy W Russell
Yang Liu
Joel Hellewell
Carl AB Pearson
CMMID COVID-19 Working Group
Gwenan M Knight
Rosalind M Eggo
Adam J Kucharski
Sebastian Funk
Stefan Flasche
Rein MGJ Houben

(2020)

The contribution of asymptomatic SARS-CoV-2 infections to transmission on the Diamond Princess cruise ship

eLife 9:e58699.

https://doi.org/10.7554/eLife.58699

Categories and tags

Research organism

Human

Share this article

Cite this article

Data from the Diamond Princess and model calibration.

Figure 1—source data 1

Proportion of infections that are asymptomatic and their contribution to transmission.

Instantaneous proportion of transmission from symptomatic, presymptomatic and asymptomatic individuals.

Model outputs by relative infectiousness of asymptomatic individuals.

Model parameters and priors/values.

Confirmed symptomatic cases (n = 199) by date of symptom onset for passengers and crew separately, extracted from Nishiura, 2020.

Confirmed pre/asymptomatic cases (n = 320) and symptom-agnostic tests (n = 2749) by date of test for passengers and crew combined, extracted from Mizumoto et al., 2020.

Model diagram for the outbreak onboard the Diamond Princess cruise ship described in the main paper.

Data from the Diamond Princess and model calibration.

Parameter correlation plot containing parameter values from 10,000 samples of the joint posterior distribution found during MCMC model calibration.

Parameter trace plot showing all 1.5 million samples from the MCMC model calibration sequentially.

Data from the Diamond Princess and model calibration.

Proportion of infections that are asymptomatic and their contribution to transmission.

Parameter correlation plot containing parameter values from 10,000 samples of the joint posterior distribution found during MCMC model calibration.

Parameter trace plot showing all 1.5 million samples from the MCMC model calibration sequentially.

Data from the Diamond Princess and model calibration.

Proportion of infections that are asymptomatic and their contribution to transmission.

Parameter correlation plot containing parameter values from 10,000 samples of the joint posterior distribution found during MCMC model calibration.

Parameter trace plot showing all 1.5 million samples from the MCMC model calibration sequentially.

Data from the Diamond Princess and model calibration.

Proportion of infections that are asymptomatic and their contribution to transmission.

Parameter correlation plot containing parameter values from 10,000 samples of the joint posterior distribution found during MCMC model calibration.

Parameter trace plot showing all 1.5 million samples from the MCMC model calibration sequentially.

Data from the Diamond Princess and model calibration.

Proportion of infections that are asymptomatic and their contribution to transmission.

Parameter correlation plot containing parameter values from 10,000 samples of the joint posterior distribution found during MCMC model calibration.

Parameter trace plot showing all 1.5 million samples from the MCMC model calibration sequentially.

Data from the Diamond Princess and model calibration.

Proportion of infections that are asymptomatic and their contribution to transmission.

Parameter correlation plot containing parameter values from 10,000 samples of the joint posterior distribution found during MCMC model calibration.

Parameter trace plot showing all 1.5 million samples from the MCMC model calibration sequentially.

Data from the Diamond Princess and model calibration.

Proportion of infections that are asymptomatic and their contribution to transmission.

Parameter correlation plot containing parameter values from 10,000 samples of the joint posterior distribution found during MCMC model calibration.

Parameter trace plot showing all 1.5 million samples from the MCMC model calibration sequentially.

Data from the Diamond Princess and model calibration.

Proportion of infections that are asymptomatic and their contribution to transmission.

Parameter correlation plot containing parameter values from 10,000 samples of the joint posterior distribution found during MCMC model calibration.

Parameter trace plot showing all 1.5 million samples from the MCMC model calibration sequentially.

Data from the Diamond Princess and model calibration.

Proportion of infections that are asymptomatic and their contribution to transmission.

Parameter correlation plot containing parameter values from 10,000 samples of the joint posterior distribution found during MCMC model calibration.

Parameter trace plot showing all 1.5 million samples from the MCMC model calibration sequentially.

Data from the Diamond Princess and model calibration.

Proportion of infections that are asymptomatic and their contribution to transmission.

Parameter correlation plot containing parameter values from 10,000 samples of the joint posterior distribution found during MCMC model calibration.

Parameter trace plot showing all 1.5 million samples from the MCMC model calibration sequentially.

Author details

Jon C Emery

Contribution

Competing interests

Timothy W Russell

Contribution

Competing interests

Yang Liu

Contribution

Competing interests

Joel Hellewell

Contribution

Competing interests

Carl AB Pearson

Contribution

Competing interests

CMMID COVID-19 Working Group

Contribution

Competing interests

Gwenan M Knight

Contribution

Competing interests

Rosalind M Eggo

Contribution

Competing interests

Adam J Kucharski

Contribution

Competing interests

Sebastian Funk

Contribution