key: cord-0941598-32v44sw9 authors: Wu, Joseph T; Leung, Kathy; Leung, Gabriel M title: Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study date: 2020-03-06 journal: The Lancet DOI: 10.1016/s0140-6736(20)30260-9 sha: a99b24730ab1768a968a582caa998ab84b96f007 doc_id: 941598 cord_uid: 32v44sw9 Summary Background Since Dec 31, 2019, the Chinese city of Wuhan has reported an outbreak of atypical pneumonia caused by the 2019 novel coronavirus (2019-nCoV). Cases have been exported to other Chinese cities, as well as internationally, threatening to trigger a global outbreak. Here, we provide an estimate of the size of the epidemic in Wuhan on the basis of the number of cases exported from Wuhan to cities outside mainland China and forecast the extent of the domestic and global public health risks of epidemics, accounting for social and non-pharmaceutical prevention interventions. Methods We used data from Dec 31, 2019, to Jan 28, 2020, on the number of cases exported from Wuhan internationally (known days of symptom onset from Dec 25, 2019, to Jan 19, 2020) to infer the number of infections in Wuhan from Dec 1, 2019, to Jan 25, 2020. Cases exported domestically were then estimated. We forecasted the national and global spread of 2019-nCoV, accounting for the effect of the metropolitan-wide quarantine of Wuhan and surrounding cities, which began Jan 23–24, 2020. We used data on monthly flight bookings from the Official Aviation Guide and data on human mobility across more than 300 prefecture-level cities in mainland China from the Tencent database. Data on confirmed cases were obtained from the reports published by the Chinese Center for Disease Control and Prevention. Serial interval estimates were based on previous studies of severe acute respiratory syndrome coronavirus (SARS-CoV). A susceptible-exposed-infectious-recovered metapopulation model was used to simulate the epidemics across all major cities in China. The basic reproductive number was estimated using Markov Chain Monte Carlo methods and presented using the resulting posterior mean and 95% credibile interval (CrI). Findings In our baseline scenario, we estimated that the basic reproductive number for 2019-nCoV was 2·68 (95% CrI 2·47–2·86) and that 75 815 individuals (95% CrI 37 304–130 330) have been infected in Wuhan as of Jan 25, 2020. The epidemic doubling time was 6·4 days (95% CrI 5·8–7·1). We estimated that in the baseline scenario, Chongqing, Beijing, Shanghai, Guangzhou, and Shenzhen had imported 461 (95% CrI 227–805), 113 (57–193), 98 (49–168), 111 (56–191), and 80 (40–139) infections from Wuhan, respectively. If the transmissibility of 2019-nCoV were similar everywhere domestically and over time, we inferred that epidemics are already growing exponentially in multiple major cities of China with a lag time behind the Wuhan outbreak of about 1–2 weeks. Interpretation Given that 2019-nCoV is no longer contained within Wuhan, other major Chinese cities are probably sustaining localised outbreaks. Large cities overseas with close transport links to China could also become outbreak epicentres, unless substantial public health interventions at both the population and personal levels are implemented immediately. Independent self-sustaining outbreaks in major cities globally could become inevitable because of substantial exportation of presymptomatic cases and in the absence of large-scale public health interventions. Preparedness plans and mitigation interventions should be readied for quick deployment globally. Funding Health and Medical Research Fund (Hong Kong, China). The National Health Commission of China has developed a case-definition system to facilitate the classification of patients (panel). To mitigate the spread of the virus, the Chinese Government has pro gressively implemented metro politan-wide quarantine of Wuhan and several nearby cities since Jan 23-24, 2020. Numerous domestic See Online for appendix airports and train stations, as well as international airports, have adopted temperature screening measures to detect individuals with fever. Two other novel coronaviruses (CoVs) have emerged as major global health threats since 2002, namely severe acute respiratory syndrome coronavirus (SARS-CoV; in 2002) that spread to 37 countries, and Middle East respiratory syndrome coronavirus (MERS-CoV; in 2012) that spread to 27 countries. SARS-CoV caused more than 8000 infections and 800 deaths, and MERS-CoV infected 2494 individuals and caused 858 deaths worldwide to date. Both are zoonotic viruses and epidemiologically similar, except that SARS-CoV has virtually no subclinical manifestation, whereas MERS-CoV behaves more similarly to the other four commonly circulating human CoVs, with a substantial proportion of asymptomatic infections (table 1). Symptomatic cases of both viruses usually present with moderate-to-severe respiratory symptoms that often progress to severe pneumonia. A notable common characteristic of both SARS-CoV and MERS-CoV is that they have low potential for sustained community transmission (ie, low basic reproductive number). 15, 28, 29 However, the most worrisome aspect is the ability of the viruses to cause unusually large case clusters via superspreading, which can exceed 100 individuals and are apparently seeded by a single index case. 5, [28] [29] [30] [31] In this study, we provide a nowcast of the probable size of the epidemic, recognising the challenge of complete ascertainment of a previously unknown pathogen with an unclear clinical spectrum and testing capacity, even after identification of the aetiological cause. More importantly, from a public health control viewpoint, we then forecast the probable course of spread domes tically and inter nationally, first by assuming similar transmissibility as the initial phase in Wuhan (ie, little or no mitigation interventions), then accounting for the potential impact of the various social and personal nonpharmaceutical interventions that have been progressively and quickly implemented in January, 2020. In this modelling study, we first inferred the basic reproductive number of 2019-nCoV and the outbreak size in Wuhan from Dec 1, 2019, to Jan 25, 2020, on the basis of the number of cases exported from Wuhan to cities outside mainland China. We then estimated the number of cases that had been exported from Wuhan to other cities in mainland China. Finally, we forecasted the spread of 2019-nCoV within and outside mainland China, accounting for the Greater Wuhan region quarantine implemented since Jan 23-24, 2020, and other public health interventions. Wuhan is the major air and train transportation hub of central China (figure 1). We estimated the daily number of outbound travellers from Wuhan by air, train, and road with data from three sources (see appendix p 1 for details): (1) the monthly number of global flight bookings to Wuhan for January and February, 2019, obtained from the Official Aviation Guide (OAG); (2) the daily number of domestic passengers by means of transportation recorded by the location-based services of the Tencent (Shenzhen, China) database from Wuhan to more than 300 prefecture-level Evidence before this study In central China, Wuhan is investigating an outbreak of atypical pneumonia caused by the zoonotic 2019 novel coronavirus (2019-nCoV). Few data and details are currently available. Two other novel coronaviruses (CoVs) have emerged as major global epidemics in recent decades; in 2002, severe acute respiratory syndrome coronavirus (SARS-CoV) spread to 37 countries and caused more than 8000 cases and almost 800 deaths, and in 2012, Middle East respiratory syndrome coronavirus (MERS-CoV) spread to 27 countries, causing 2494 cases and 858 deaths worldwide to date. These CoVs are both zoonotic and have low potential for sustained community transmission. However, larger clusters involving superspreading events were observed for SARS in 2003 and MERS in 2014-17. As of Jan 25, 2020, the scale and geographical extent of the 2019-nCoV outbreak both within and outside mainland China are highly uncertain. National and global spread of this disease is particularly concerning given that chunyun, a 40-day period with extremely high air and train traffic across China because of the lunar new year Spring Festival, began on Jan 10, 2020. We searched PubMed and preprint archives for articles published up to Jan 25, 2020, that contained information about the Wuhan outbreak using the terms "coronavirus", "CoV", "2019-nCoV", "Wuhan", "transmission", "China", "superspreading", and "Chinese New Year". We found six studies that reported the relative risks of case exportation from Wuhan to areas outside mainland China. In the absence of a robust and complete line list for characterising the epidemiology of this novel pathogen, we inferred the outbreak size of 2019-nCoV in Wuhan from the number of confirmed cases that have been exported to cities outside mainland China. We used this outbreak size estimate to project the number of cases that have been exported to other Chinese cities. We forecasted the spread of 2019-nCoV both within and outside of mainland China. 32 On Jan 19, 2020, the Chinese Center for Disease Control and Prevention (CDC) reported that only 43 (22%) of the 198 confirmed cases in its outbreak investigation had been exposed to the Huanan seafood wholesale market, 33 the most probable index source of zoonotic 2019-nCoV infections, which was closed and disinfected on Jan 1, 2020. As such, and because of the difficulties of tracing all infections, we assumed that during Dec 1-31, 2019, the epidemic in Wuhan was seeded by a constant zoonotic force of infection that caused 86 cases (ie, twice the 43 confirmed cases with zoonotic exposure) in our baseline scenario. For the sensitivity analysis, we assumed 129 and 172 cases (50% and 100% higher than the baseline scenario value). Given that both 2019-nCoV and SARS-CoV could cause self-sustaining human-tohuman transmission in the community, we assumed that the serial interval of 2019-nCoV was the same as that of SARS-CoV in Hong Kong (mean 8·4 days; table 1). 4 We assumed that the incubation period of 2019-nCoV was similar to that of SARS-CoV and MERS-CoV (mean 6 days; table 1). These assumptions are consistent with preliminary estimates of the serial interval (mean 7·5 days) and incubation period (mean 6·1 days) using line-list data from China CDC. 34 We assumed that the catchment population size of the Wuhan Tianhe International Airport at Wuhan was 19 million (ie, Greater Wuhan region with 11 million people from Wuhan city plus 8 million from parts of several neighbouring cities). We estimated the number of cases in Greater Wuhan on the basis of the number of confirmed cases exported to cities outside mainland China whose symptom onset date had been reported to fall from Dec 25, 2019, to Jan 19, 2020 (appendix pp 2-4). The start date of this period (day D s ) corresponded to one mean incubation period before the Wuhan outbreak was announced on Dec 31, 2019, whereas the end date (day D e ) was 7 days before the time of writing (Jan 26, 2020). This end date was chosen to minimise the effect of lead time bias on case confirmation (the time between onset and case confirmation was ≤7 days in 46 [80%] of 55 cases exported to cities outside mainland China; see appendix pp 2-4). We let χ d to be the number of such case exportation on day d. Our 2019 OAG data indicated that for cities outside mainland China excluding Hong Kong, the daily average number of international outbound air passengers was L W,I =3633 and that of international inbound air passengers was L I,W =3546 in Greater Wuhan during January-Feb ru ary, 2019 (table 2) . We excluded Hong Kong in this estimate because travel from mainland China to Hong Kong had dropped sharply since August, 2019, because of social unrest in Hong Kong. Our calibrated Tencent mobility data indicated that for cities in mainland China, the daily number of all domestic outbound travellers was L W,C (t)=502 013 and that of all domestic inbound travellers was L C,W (t)=487 310 in Wuhan at time t before chunyun (Jan 10). During chunyun, these estimates were L W,C (t)=717 226 and L C,W (t)=810 500. We used the following susceptible-exposed-infectiousrecovered (SEIR) model to simulate the Wuhan epidemic since it was established in December, 2019: where S(t), E(t), I(t), and R(t) were the number of susceptible, latent, infectious, and removed individuals at time t; D E and D I were the mean latent (assumed to be the same as incubation) and infectious period (equal to the serial interval minus the mean latent period 4 ); R 0 was the basic reproductive number; z(t) was the zoonotic The case definition of 2019-nCoV differs depending on the context in which it is used. A suspected or probable case is defined as a case that meets: (1) three clinical criteria or (2) two clinical criteria and one epidemiological criterion. Clinical criteria are: fever; radiographic evidence of pneumonia or acute respiratory distress syndrome; and low or normal white blood cell count or low lymphocyte count. Epidemiological criteria are: living in Wuhan or travel history to Wuhan within 14 days before symptom onset; contact with patients with fever and symptoms of respiratory infection within 14 days before symptom onset; and a link to any confirmed cases or clusters of suspected cases. The definition of a confirmed case, for the first case in a province, is a suspected or probable case with detection of viral nucleic acid at the city CDC and provincial CDC. For the second case and all subsequent cases, the definition is a suspected or probable case with detection of virus nucleic acid at the city CDC. We defined cases as individuals who were symptomatic, who could be detected by temperature screening at international borders, or who had a disease severity requiring hospital admission, or both, plus travel history to Wuhan. force of infection equal to 86 cases per day in the baseline scenario before market closure on Jan 1, 2020, and equal to 0 thereafter. The cumulative number of infections and cases that had occurred in Greater Wuhan up to time t was obtained from the SEIR model. We assumed that travel behaviour was not affected by disease and hence international case exportation occurred according to a non-homogeneous process with rate As such, the expected number of international case exportation on day d was Taken together, the likelihood function was We estimated R 0 using Markov Chain Monte Carlo methods with Gibbs sampling and non-informative flat prior. For a given R 0 from the resulting posterior distribution, we used the same SEIR model to estimate the corresponding outbreak size in Wuhan and the probability distribution of the number of cases that had been exported domestically to other cities in mainland China (on the basis of the destination distribution in our calibrated Tencent mobility data). Point estimates were presented using posterior means, and statistical uncertainty was presented using 95% credible intervals (CrIs). We extended the above SEIR model into a SEIR-metapopulation model to simulate the spread of 2019-nCoV across mainland China, assuming the trans missibility of 2019-nCoV was similar across all cities. 35 The movements of individuals between more than 300 prefecture-level cities were modelled using the daily average traffic volumes in our calibrated Tencent mobility data. Given that 2019-nCoV has caused widespread outbreak awareness not only among public health professionals (ie, WHO and government health authorities), but also among the general public in China and other countries, the transmissibility of the epidemic might be reduced compared with its nascent stage at Wuhan because of community-wide social distancing measures and other non-pharmaceutical interventions (eg, use of face masks and improved personal hygiene). Previous studies sug gested that non-pharmaceutical interventions might be able to reduce influenza transmission by up to 50%. 36 As such, we simulated local epidemics across mainland China assuming that the transmissibility of 2019-nCoV was reduced by 0%, 25%, and 50% after Wuhan was quarantined on Jan 23, 2020. The epidemics would fade out if transmissibility could be reduced by 1-1/R 0 . Further more, we considered 50% reduction in intercity mobility. Finally, we hypothesised that citywide population quarantine at Wuhan had negligible effect on the epidemic trajectories of the rest of the country and tested this hypothesis by making the extreme assumption that all inbound and outbound mobility at Wuhan were eliminated indefinitely on Jan 23, 2020. 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. Figure 2 summarises our estimates of the basic reproductive number R 0 and the outbreak size of 2019-nCoV in Wuhan as of Jan 25, 2020. In our baseline scenario, we estimated that R 0 was 2·68 (95% CrI 2·47-2·86) with an epidemic doubling time of 6·4 days (95% CrI 5·8-7·1; figure 2 ). We estimated that 75 815 individuals (95% CrI 37 304-130 330) individuals had been infected in figure 3) . Beijing, Shanghai, Guangzhou, and Shenzhen were the mainland Chinese cities that together accounted for 53% of all outbound international air travel from China and 69% of international air travel outside Asia, whereas Chongqing is a large metropolis that has a population of 32 million and very high ground traffic volumes with Wuhan. Substantial epidemic take-off in these cities would thus contribute to the spread of 2019-nCoV within and outside mainland China. If the zoonotic force of infection that initiated the Wuhan epidemic was 50% and 100% higher than the baseline scenario value, then R 0 would be 2·53 (95% CrI 2·32-2·71) and 2·42 (2·22-2·60), respectively. The corresponding estimate of the number of infections in Wuhan would be 38% and 56% lower than baseline. The number of exported cases and infections in Chongqing, Beijing, Shanghai, Guangzhou, and Shenzhen would be similarly reduced in magnitude ( figure 3) . Figure 4 shows the epidemic curves for Wuhan, Chongqing, Beijing, Shanghai, Guangzhou, and Shenzhen with a R 0 of 2·68, assuming 0%, 25%, or 50% decrease in transmissibility across all cities, together with 0% or 50% reduction in inter-city mobility after Wuhan was quarantined on Jan 23, 2020. The epidemics would fade out if transmissibility was reduced by more than 1-1/R 0 =63%. Our estimates suggested that a 50% reduction in inter-city mobility would have a negligible effect on epidemic dynamics. We estimated that if there was no reduction in trans missibility, the Wuhan epidemic would peak around April, 2020, and local epidemics across cities in mainland China would lag by 1-2 weeks. If transmissibility was reduced by 25% in all cities domestically, then both the growth rate and magnitude of local epidemics would be substantially reduced; the epidemic peak would be delayed by about 1 month and its magnitude reduced by about 50%. A 50% reduction in transmissibility would push the viral reproductive number to about 1·3, in which case the epidemic would grow slowly without peaking during the first half of 2020. However, our simulation suggested that wholesale quarantine of population movement in Greater Wuhan would have had a negligible effect on the forward City C h o n g q i n g B e i j i n g S h a n g h a i G u a n g z h o u S h e n z h e n C h e n g d u Z h e n g z h o u C h a n g s h a N a n c h a n g H e f e i N a n j i n g H a n g z h o u W e n z h o u K u n m i n g X i a m e n H a i k o u T i a n j i n N a n n i n g X i 'a n Q i n g d a o Base case 50% higher zoonotic FOI 100% higher zoonotic FOI Estimates are as of Jan 25, 2020. Cases corresponded to infections that were symptomatic or infectious. The number of cases was smaller than the number of infections because some individuals with the infection were still in the incubation period. We assumed that infected individuals were not infectious during the incubation period (ie, similar to severe acute respiratory syndrome-related coronavirus 37 ). PDF=probability density function. FOI=force of infection. Base case 50% higher zoonotic FOI 100% higher zoonotic FOI trajectories of the epidemic because multiple major Chinese cities had already been seeded with more than dozens of infections each (results not shown because they are visually indistinguishable from figure 4) . The probability that the chain of transmission initiated by an infected case would fade out without causing exponential epidemic growth decreases sharply as R 0 increases (eg, <0·2 when R 0 >2). 1, 38 As such, given the substantial volume of case importation from Wuhan (figure 3), local epidemics are probably already growing exponentially in multiple major Chinese cities. Given that Beijing, Shanghai, Guangzhou, and Shenzhen together accounted for more than 50% of all outbound interna tional air travel in mainland China, other countries would likely be at risk of experiencing 2019-nCoV epidemics during the first half of 2020. During the period of an epidemic when human-to-human transmission is established and reported case numbers are rising exponentially, nowcasting and forecasting are of crucial importance for public health planning and control domestically and internationally. 39, 40 In this study, we have estimated the outbreak size of 2019-nCoV thus far in Wuhan and the probable extent of disease spread to other cities domestically. Our findings suggest that independent self-sustaining human-to-human spread is already present in multiple major Chinese cities, many of which are global transport hubs with huge numbers of both inbound and outbound passengers (eg, Beijing, Shanghai, Guangzhou, and Shenzhen). Therefore, in the absence of substantial public health interventions that are immediately applied, further international seeding and subsequent local establishment of epidemics might become inevitable. On the present trajectory, 2019-nCoV could be about to become a global epidemic in the absence of mitigation. Nevertheless, it might still be possible to secure containment of the spread of infection such that initial imported seeding cases or even early local transmission does not lead to a large epidemic in locations outside Wuhan. To possibly succeed, substantial, even draconian measures that limit population mobility should be seriously and immediately considered in affected areas, as should strategies to drastically reduce within-population contact rates through cancellation of mass gatherings, school closures, and instituting work-from-home arrange ments, for example. Precisely what and how much should be done is highly contextually specific and there is no one-size-fits-all set of prescriptive interventions that would be appropriate across all settings. Should containment fail and local transmission is established, mitigation measures according to plans that had been drawn up and executed during previous major outbreaks, such as those of SARS, MERS, or pandemic influenza, could serve as useful reference templates. The overriding epidemiological priority to inform public health control would be to compile and release a line list of suspected, possible, probable, and confirmed cases and close contacts that is updated daily and linked to clinical outcomes and laboratory test results. A robust line list is essential for the generation of accurate and precise epidemiological parameters as inputs into transmission models to inform situational awareness and optimising the responses to the epidemic. 41 Additionally, given the extent of spread and level of public concern it has already generated, the clinical spectrum and severity profile of 2019-nCoV infections needs rapid ascertainment by unbiased and reliable methods in unselected samples of cases, especially those with mild or subclinical presentations. The modelling techniques that we used in this study are very similar to those used by other researchers who are working towards the same goal of characterising the epidemic dynamics of 2019-nCoV (Zhanwei Du, University of Texas at Austin, personal communication). [42] [43] [44] [45] The consensus on our methodology provides some support for the validity of our nowcasts and forecasts. An additional strength of our study is that our model is parameterised with the latest mobility data from OAG and Tencent. Nonetheless, our study has several major limitations. First, we assumed that travel behaviour was not affected by disease status and that all infections eventually have symptoms (albeit possibly very mild). We would have underestimated the outbreak size in Greater Wuhan if individuals with increased risk of infection (eg, confounded by socioeconomic status) were less likely to travel internationally or if the proportion of asymptomatic infections were substantial. Second, our estimate of transmis sibility and outbreak size was somewhat sensitive to our assumption regarding the zoonotic mechanism that initiated the epidemic at Wuhan. However, our overall conclusion regarding the extent of case exportation in major Chinese cities would remain the same even for our lowest estimate of transmissibility ( figure 3) . Third, our epidemic forecast was based on inter-city mobility data from 2019 that might not necessarily reflect mobility patterns in 2020, especially in the presence of current public vigilance and response regarding the health threat posed by 2019-nCoV (appendix p 5). Fourth, little is known regarding the seasonality of coro navirus transmission. If 2019-nCoV, similar to influenza, has strong seasonality in its transmission, our epidemic forecast might not be reliable. Identifying and eliminating the zoonotic source remains an important task to prevent new animal-tohuman seeding events. The renewal of a complete ban on market trading and sale of wild game meat in China on Jan 26 can provide only temporary suspension of demand, even if completely adhered to. 46 Vaccine platforms should be accelerated for real-time deployment in the event of a second wave of infections. Above all, for health protection within China and internationally, especially those locations with the closest travel links with major Chinese ports, preparedness plans should be readied for deployment at short notice, including securing supply chains of pharmaceuticals, personal pro tective equipment, hospital supplies, and the necessary human resources to deal with the consequences of a global outbreak of this magnitude. JTW, GML, and KL designed the experiments. KL collected data. JTW and KL analysed data. KL, JTW, and GML interpreted the results and wrote the manuscript. We declare no competing interests. We collated epidemiological data from publicly available data sources (news articles, press releases, and published reports from public health agencies). All the epidemiological information that we used is documented in the Article. We purchased the 2019 global flight bookings data from the Official Aviation Guide (OAG). Our data purchase agreement with OAG prohibits us from sharing these data with third parties, but interested parties can contact OAG to make the same data purchase. We purchased the January-February, 2019, intercity mobility data from Tencent. Our data purchase agreement with Tencent prohibits us from sharing these data with third parties, but interested parties can contact Tencent to make the same data purchase. Superspreading and the effect of individual variation on disease emergence Model parameters and outbreak control for SARS Transmission dynamics of the etiological agent of SARS in Hong Kong: impact of public health interventions Transmission dynamics and control of severe acute respiratory syndrome Transmission characteristics of MERS and SARS in the healthcare setting: a comparative study The role of superspreading in Middle East respiratory syndrome coronavirus (MERS-CoV) transmission Unraveling the drivers of MERS-CoV transmission MERS transmission and risk factors: a systematic review Clinical impact of human coronaviruses 229E and OC43 infection in diverse adult populations Epidemiology characteristics of human coronaviruses in patients with respiratory infection symptoms and phylogenetic analysis of HCoV-OC43 during 2010-2015 in Guangzhou Epidemiology and clinical presentations of the four human coronaviruses 229E, HKU1, NL63, and OC43 detected over 3 years using a novel multiplex real-time PCR method The epidemiology of severe acute respiratory syndrome in the 2003 Hong Kong epidemic: an analysis of all 1755 patients A comparative epidemiologic analysis of SARS in Hong Kong, Beijing and Taiwan Comparison of incubation period distribution of human infections with MERS-CoV in South Korea and Saudi Arabia Complete genomic sequence of human coronavirus OC43: molecular clock analysis suggests a relatively recent zoonotic coronavirus transmission event Incubation periods of acute respiratory viral infections: a systematic review Murray and Nadel's textbook of respiratory medicine, 5th edn Preliminary epidemiological assessment of MERS-CoV outbreak in South Korea Seroprevalence of IgG antibodies to SARS-coronavirus in asymptomatic or subclinical population groups Comparative serological study for the prevalence of anti-MERS coronavirus antibodies in high-and low-risk groups in Qatar Middle East respiratory syndrome coronavirus (MERS-CoV) neutralising antibodies in a high-risk human population Development and application of an enzyme immunoassay for coronavirus OC43 antibody in acute respiratory illness Examination of seroprevalence of coronavirus HKU1 infection with S protein-based ELISA and neutralization assay against viral spike pseudotyped virus Middle East respiratory syndrome: what we learned from the 2015 outbreak in the Republic of Korea An outbreak of human coronavirus OC43 infection and serological cross-reactivity with SARS coronavirus Human coronavirus alone or in co-infection with rhinovirus C is a risk factor for severe respiratory disease and admission to the pediatric intensive care unit: a one-year study in Southeast Brazil Case fatality of SARS in mainland China and associated risk factors Interhuman transmissibility of Middle East respiratory syndrome coronavirus: estimation of pandemic risk Middle East respiratory syndrome coronavirus: quantification of the extent of the epidemic, surveillance biases, and transmissibility The SARS epidemic in Hong Kong: what lessons have we learned? SARS in healthcare facilities People's Daily Online. Wuhan is expected to send 15 million passengers during Spring Festival in 2020 The 2019-nCoV Outbreak Joint Field Epidemiology Investigation Team, Li Q. An outbreak of NCIP (2019-nCoV) infection in China-Wuhan Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia Characterizing the dynamics underlying global spread of epidemics The effect of public health measures on the 1918 influenza pandemic in U.S. cities Frequently asked questions about SARS Mathematical epidemiology of infectious diseases: model building, analysis and interpretation Real-time epidemic forecasting: challenges and opportunities Outbreak analytics: a developing data science for informing the response to emerging pathogens Improving the evidence base for decision making during a pandemic: the example of 2009 influenza A/H1N1 MRC Centre for Global Infectious Disease Analysis: Wuhan coronavirus reports 1-3 Novel coronavirus 2019-nCoV: early estimation of epidemiological parameters and epidemic predictions Transmission dynamics of 2019 novel coronavirus (2019-nCoV) Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: a data-driven analysis in the early phase of the outbreak Announcement of the State Administration of Market Supervision, Ministry of Agriculture and Rural Affairs, and the National Forestry and Grass Bureau on banning wildlife trading (in Chinese) We thank Chi-Kin Lam and Miky Wong from School of Public Health, The University of Hong Kong (Hong Kong, China) for technical support. This research was supported by a commissioned grant from the Health and Medical Research Fund from the Government of the Hong Kong Special Administrative Region. Editorial note: the Lancet Group takes a neutral position with respect to territorial claims in published maps and institutional affiliations.