Estimating the establishment of local transmission and the cryptic phase of the COVID-19 pandemic in the USA

Onset of local transmission

Stochastic simulations of the worldwide epidemic spread yield international/domestic infection importations, incidence of infections, and deaths per subpopulation at a daily resolution in the continental US. In Fig.1B we show the model estimates for the median daily incidence of new infections up to March 1, 2020 in the continental US. There is a stark contrast between the model output and the number of officially reported cases by the same date, highlighting the significant number of transmission events that may have already occurred before many states had implemented testing strategies independent of travel history. For model validation we compare our projections of the number of infections during the week of March 8, 2020, to the number of cases reported during that week within 48 of the continental states that reported at least 1 case (shown in Fig.1B inset). We see a strong correlation between the reported cases and our model’s projected number of infections, (Pearson’s correlation coefficient on log-values 0.79, p < 0.001), although many fewer cases had actually been reported by that time. If we assume that the number of reported cases and simulated infections are related through a simple binomial stochastic sampling process, we find that on average 5 in 1,000 infections (90%CI [1 – 25 per 1,000]) were detected by March 8, 2020. The ascertainment rate grows quickly as testing capacity increases, and our estimate doubled to detecting 11 in 1,000 infections (90%CI [2 – 36 per 1, 000]) by March 15, 2020. By April 1, 2020 our model suggests a detection rate of 96 cases per 1,000 infections (90%CI: [20- 289 per 1,000]). The Centers for Disease Control and Prevention (CDC) estimates that 42 to 166 in 1,000 infections were detected during March and April, 2020, at different locations in the US, with estimated numbers of infections at least 10 times greater than the number of reported cases in locations across seven different states (37, 38). These SARS-CoV-2 infections are also distributed heterogeneously across the US. In Fig.1C we show the model-based median daily number of new infections on March 1, 2020.

In Fig. 2 we plot the posterior distribution for the earliest date within each state when at least 10 new infections per day occurred in the community. California, New York and New Jersey are the first states with a probability larger than 50% to have experienced local transmission by mid January (California) or beginning of February (New York and New Jersey), 2020. However, the wide distribution of dates suggests that SARS-CoV-2 cryptic transmission may have started as early as December, 2019. The posterior distribution for the timing of the onset of local transmission peaks by mid February in Florida, Illinois, Maryland, Massachusetts, Texas, Washington, Arizona, Colorado, Connecticut, Georgia, Indiana, Michigan, Minnesota, Nevada, New Hampshire, North Carolina, Ohio, Oregon, Pennsylvania, Tennessee, and Utah Virginia. From the posterior distribution of Fig. 2, all states have a median date of the onset of local transmission by early March, with the large majority of them in February, 2020, a critical month for the cryptic spread of SARS-CoV-2 in the continental US. However, during that time, testing in the US was still focused on returning travelers from China.

• Download figure
• Open in new tab

Figure 2: Timing of the onset of local transmission.

Posterior distributions of the week when each state first reached 10 locally generated SARS-CoV-2 transmission events per day.

International and domestic seeding

As the model allows the recording of the origin and destination of SARS-CoV-2 carriers at the global scale, we can study the possible sources of infection importation for each state. In particular, we are able to record the flow of latent and infectious individuals through international and domestic flight connections. However, the model also considers the effect of possible infections through commuters that may spend a few hours in a neighboring subpopulation (7). It is important to stress that the model’s realizations explore the many possible paths of the epidemic. Thus, the analysis provided here must be considered as a statistical description of the potential sources of SARS-CoV-2 importations, rather than providing a specific, single causal chain of events.

To identify seeding importation events relevant to the onset of local transmission in each stochastic realization of the model, we record the number of importations (of latent and infectious individuals) before the local transmission chains were established (defined as 10 daily local transmission events). We visualize the origin of the seeding importations relevant for establishing local transmission by aggregating the importation sources, considering some key geographical regions (e.g. Europe and Asia) while keeping the US and China separate and aggregating all the other countries (i.e. Others) in Fig. 3. It is worth clarifying that seeding importations are different from the actual number of times the virus has been introduced to each state with subsequent onward transmission. Even after local transmission has started, future importation events may give rise to additional onward transmission forming independently-introduced transmission lineages of the virus as seen in the United Kingdom (39). Statistics for importations through March 1, 2020, are reported in the SM file.

• Download figure
• Open in new tab

Figure 3:

Importation sources. Each state is displayed in a clockwise order with respect to the start of the local outbreak (as seen in Fig. 2). Importation flows are directed and weighted. We normalize links considering the total in-flow for each state so that the sum of importations flows, for each state, is one. In the SM we report the complete list of countries contributing, as importation sources, in each group (i.e., geographical region).

Cryptic spreading phase

From late January to early March, SARS-CoV-2 had been spreading across the US mostly undetected. We estimate that cities such as Los Angeles, New York, Chicago, Seattle, and San Francisco have possibly experienced local transmission beginning in the first half of February (Fig. 4A). In the time span of one additional month, through March 15, 2020, most large cities in the US had sizeable ongoing outbreaks. California most likely generated 10 local infections by late January (see Fig. 2), around the same time as the adoption of draconian containment measures and international travel restrictions in mainland China. The model also allows us to estimate possible COVID-19 related deaths. By March 1, 2020, we estimate a median of 63 [90% CI10 – 629] cumulative deaths in the US while only 1 death was reported. Although some states launched investigations in search of evidence that COVID-19 was the cause of death as far back as December 2019, it is likely that most deaths were missed because COVID-19 testing guidelines were based on travel history (40).

• Download figure
• Open in new tab

Figure 4: Comparing model projections to surveillance and air traffic data.

(A) Map of the US showing the date where regions observed at least an average of 100 infections. (B) The correlation between the ordering of each state to reach 100 infections in the model projections and to reach 100 reported cases in the surveillance data. The correlation is computed using the Kendall rank correlation coefficient, τ. (C) The correlation between the ordering of each state considering the time needed to reach 100 reported cases in the surveillance data and the ranking of the combined international and domestic air traffic. (D) The order correlation of the onset of excess deaths due to pneumonia and influenza and the order of the model projection for the date of 10 transmissions per day. Circle size in A, B, and C corresponds to the population size of each state. (E) The cumulative distribution of the probability that the US reached 100 locally generated transmissions per day by a given day for both a scenario using the real timelines of travel restrictions and a counterfactual scenario where the travel restrictions to/from mainland China are shifted one week before. The median dates for the real and counterfactual timelines are February 2, 2020, and February 4, 2020, respectively. We show a horizontal line at a density of 0.5

To provide a model consistency check with respect to surveillance data concerning the epidemic, we compared the model-based estimates with observed surveillance data. Particularly, early on in the epidemic, surveillance data were known to be highly unreliable because of under-detection. For each state we compared the order in which they surpassed 100 infections in the model and in the surveillance data gathered from the John Hopkins University Coronavirus Resource Centre (41). In Fig. 4B we plot the ordering for states and compute the Kendall rank correlation coefficient τ (see SM for details). The correlation is positive (τ = 0.69, p < 0.001) indicating that, despite the detection and testing issues, the expected patterns of epidemic diffusion of the states is largely described by the model. As mentioned before, one possible major driver of the diffusion pattern is air traffic. We compare the ordering of states according to their air travel volume to their epidemic order as previously defined (Fig. 4C). We consider both national and international traffic and find a positive correlation (τ = 0.67 with p < 0.001) between the epidemic ordering derived from surveillance data and air traffic, suggesting the passenger volume of both international and national traffic are key factors driving the early phases of outbreak across the country. Similar observations have been reported in China, where the initial spreading of the virus outside Hubei was strongly correlated with the traffic to/from the province (42). The correlation between the air travel ordering and the simulations is high (τ = 0.87, p < 0.001) as actual airline data are used in the model to simulate the mobility of individuals (see SM). It is worth remarking that population size is also correlated with both the traveling flows (τ = 0.7, p < 0.001) and the epidemic order of each state (τ = 0.68, p < 0.001) as discussed in the SM. In our model it is not possible to exclude increased contacts in highly populated places before social distancing interventions and disentangle this effect from increased seeding due to the correlation between travel volume and population size. As yet another independent test for the model, in Fig. 4D we show the positive correlation between excess deaths, as estimated in (43), and the order in which the states reached the threshold of 10 daily transmissions (τ = 0.51, p < 0.001).

The model is consistent with the picture emerging from the genomic epidemiology (44–46) of an early start of the COVID-19 epidemic in coastal states, followed by the propagation, dominated by the domestic traveling patterns, to the less globally connected regions in the US. The model suggests that COVID-19 spread across the US in about 7 weeks, and that by the end of February many states were experiencing sustained SARS-CoV-2 local transmission. The median time for the onset of local transmission in coastal states is as early as February 2020, and it raises the question of what would the unfolding of the epidemic have looked like in the case of an earlier warning and issuing of travel restrictions to/from China. For this reason we have performed simulations of a counterfactual scenario in which the timeline of all the travel restrictions to/from mainland China and interventions is shifted one week earlier compared to the actual one. In Fig. 4E, we show the distribution of the probability that the US experiences sustained local transmission by a given date in both the counterfactual scenario and the real timeline of interventions. We define sustained national transmission as 100 new infections per day. Interventions implemented one week earlier amount to a proportional delay of the onset of local transmission in the continental US. In particular we find that in the real timeline intervention analysis a 50% cumulative probability of reaching the 100 locally generated transmissions mark is reached by February 2, 2020, compared to February 4, 2020, in the counterfactual scenario. This is in agreement with the evidence provided by several studies (12, 47, 48) that a considerable number of infections had already traveled from mainland China to international destinations before mid January, thus potentially seeding multiple epidemic outbreaks across the world, and leading to the international expansion of the COVID-19 epidemic, despite the mainland China travel ban. Our analysis however does not consider spontaneous behavioral changes that people might have adopted before the official national and local guidelines were announced. While certainly some individuals might have taken precautions as a result of the news from China, evidence from surveys of public concern from several countries in Europe suggest that in late February only a very limited fraction of people considered COVID-19 as a concrete threat (49). Our model also does not contain any calibration or constraint on the trajectory of the outbreak in the months of March and April. We provide this analysis in the SM showing the consistency of the results.

Discussion

Our study characterizes the cryptic transmission phase during which SARS-CoV-2 spread largely undetected in the US. The results suggest that the first sustained local transmission chains took place as early as mid January, and by the end of February the infection was spreading to many other domestic locations. This timeline is shifted several weeks ahead with respect to the detection of cases in surveillance data. This is consistent with the fact that in January and February no country had the capacity to do mass testing. Countries adopted a policy of testing symptomatic individuals with a travel history linked to China, thus, generally missing the cryptic transmission occurring domestically. We find that the order in which the virus initially progressed across states according to our model is highly correlated with the official record. The model highlights that the geographical heterogenities in the observed spreading patterns are explained by the features of the air transportation network and population distributions. The results also indicate that the source of introduction of SARS-CoV-2 infections into the US changed substantially and rapidly through time. While early importations were from international sources, most introductions occurred during February and March 2020. Our results indicate that many states were seeded from domestic sources rather than international. The presented results could be of potential interest in combination with sequencing data of SARS-CoV-2 genomes to reconstruct in greater detail the early epidemic history of the US COVID-19 epidemic. The estimated SARS-CoV-2 importation pattern and the cryptic transmission phase dynamic are of potential use when planning and modelling public health policies in the context of international travel.

Author Contributions

Author contributions: J.T.D., M.C., N.P. and A.V. designed research; M.C., J.T.D., N.P., K.M., A.P.P, M.A., N.E.D., C.G., M.L., S.M.,L.R., K.S., X.X., M.E.H., I.M.L., C.V., and A.V. per-formed research; M.C., J.T.D., N.P., A.P.P., K.M. and A.V. analyzed data; and M.C., J.T.D., N.P., K.M., A.P.P., M.A., N.E.D., C.G., M.L., S.M., L.R., K.S., X.X., M.E.H., I.M.L., C.V., and A.V. wrote and edited the paper.

Competing Interests

M.E.H. reports grants from National Institute of General Medical Sciences, during the conduct of the study; A.V. reports grants and personal fees from Metabiota inc., outside the submitted work; M.C. and A.P.P. report grants from Metabiota inc., outside the submitted work. No other relationships or activities that could appear to have influenced the submitted work.

Data and materials availability

Proprietary airline data are commercially available from Official Aviation Guide (OAG) and IATA databases. The GLEAM model is publicly available at http://www.gleamviz.org/.