Skip to main content

ORIGINAL RESEARCH article

Front. Med., 19 June 2020
Sec. Infectious Diseases – Surveillance, Prevention and Treatment
This article is part of the Research Topic Coronavirus Disease (COVID-19): Pathophysiology, Epidemiology, Clinical Management and Public Health Response View all 400 articles

Strongly Heterogeneous Transmission of COVID-19 in Mainland China: Local and Regional Variation

\nYuke Wang
Yuke Wang*Peter TeunisPeter Teunis
  • Center for Global Safe Water, Sanitation, and Hygiene, Hubert Department of Global Health, Rollins School of Public Health, Emory University, Atlanta, GA, United States

Background: The outbreak of novel coronavirus disease 2019 (COVID-19) started in the city of Wuhan, China, with a period of rapid initial spread. Transmission on a regional and then national scale was promoted by intense travel during the holiday period of the Chinese New Year. We studied the variation in transmission of COVID-19, locally in Wuhan, as well as on a larger spatial scale, among different cities and even among provinces in mainland China.

Methods: In addition to reported numbers of new cases, we have been able to assemble detailed contact data for some of the initial clusters of COVID-19. This enabled estimation of the serial interval for clinical cases, as well as reproduction numbers for small and large regions.

Findings: We estimated the average serial interval was 4.8 days. For early transmission in Wuhan, any infectious case produced as many as four new cases, transmission outside Wuhan was less intense, with reproduction numbers below two. During the rapid growth phase of the outbreak the region of Wuhan city acted as a hot spot, generating new cases upon contact, while locally, in other provinces, transmission was low.

Interpretation: COVID-19 is capable of spreading very rapidly. The sizes of outbreak in provinces of mainland China mainly depended on the numbers of cases imported from Wuhan as the local reproduction numbers were low. The COVID-19 epidemic should be controllable with appropriate interventions (suspension of public transportation, cancellation of mass gatherings, implementation of surveillance and testing, and promotion of personal hygiene and face mask use).

1. Introduction

In December 2019, several cases of severe pneumonia appeared in Wuhan, the capital city of Hubei province in China. The outbreak was caused by a novel coronavirus: Severe Acute Respiratory Syndrome Coronavirus 2, or SARS-CoV-2 (1) and the disease (COVID-19) started to spread rapidly (2) in Wuhan. As of 28 May 2020, globally 5,593,631 cases have been confirmed with 353,334 deaths in 227 countries (3). In the early outbreak in China, by 19 February 2020, 83.2% (62,031/74,576) of the confirmed cases were located in Hubei province and 60.4% (45,027/74,576) were located in the city of Wuhan, where the outbreak originated (4, 5). The “pneumonia of unknown etiology” appeared in Wuhan from 8 December 2019. Many early cases have been reported to be linked to the Huanan (Southern China) Seafood Wholesale Market (hereafter referred to as “the Market”) (6). By 2 January 2020, 41 initial cases were confirmed as COVID-19 (7).

As the novel coronavirus epidemic was spreading within Wuhan, the Chinese New Year/Spring Festival (25 January 2020), the most important holiday in China, was approaching. In 2019, 2.99 billion of people traveled by bus, train, and plane during 40 days around Chinese New Year (8). Wuhan, with a population size of 11 million, is one of the four most important railway hubs in China. With billions of people traveling and lots of family and friends gathering, there was greatly increased risk of rapidly spreading this newly emerging infectious disease, nationally, and even globally. On 19 January 2020, the first confirmed COVID-19 case outside of Wuhan appeared in Shenzhen, Guangdong (9). As of 23 January 2020, confirmed COVID-19 cases have been reported in 29 provinces of mainland China and nine countries and areas outside of mainland China (10, 11). On the same day, the central government of the People's Republic of China initiated a lockdown in Wuhan and two nearby prefectural cities, Huanggang and Ezhou, in Hubei province, to prevent spreading of the COVID-19 outbreak (12, 13). However, as many cases had “escaped” from Wuhan before the lockdown, COVID-19 has spread to most provinces. The numbers of cases imported from Wuhan for different provinces heavily depend on their connectivities with Wuhan. At the same time, municipalities responded differently, regarding timeliness and adequacy of measures (including declarations of public health emergency, holiday extension, event cancelation, and surveillance using infrared thermometers in public spaces.

In this study, we used a method for analysing transmission patterns based on the serial interval between clinical cases of COVID-19 (14). Based on dates of symptom onset and contact information for confirmed COVID-19 cases in Tianjin province, it was possible to estimate the serial interval distribution. With this information, it was possible to estimate numbers of new cases caused by any subject infectious with COVID-19, or their reproduction numbers during early stages of the outbreak in Wuhan. More importantly, we examined the heterogeneity in transmission among 30 provinces of mainland China and also among 20 cities in Guangdong province.

2. Methods

2.1. Data Sources and Assumptions

Various data with different levels of detail were collected from different sources.

First, detailed data for 112 confirmed cases between 21 January 2020 and 12 February 2020 in Tianjin province was obtained from the website of Tianjin Health Commission, the website of Tianjin Centers for Disease Control and Prevention, and the official Weibo (China's Twitter equivalent) of a local newspaper (1517). Collected information included demographic information, travel history to Wuhan/Hubei, date of symptom onset, and any contact information that could be found. See Table 1 for a summary.

TABLE 1
www.frontiersin.org

Table 1. Characteristics of COVID-19 outbreak in Tianjin by 12 February 2020.

Second, to evaluate the role of “the Market” in Wuhan and estimate the numbers of cases caused by contact with this source (its reproduction number) in early transmission, data for the first confirmed 425 COVID-19 cases with date of symptom onset and exposure information to “the Market” was extracted from a recently published report (6).

And in third place, we collected epidemic curve data: numbers of cases (total of 6941) by province and date of confirmation from the Wikipedia page of COVID-19 cases in mainland China (18) and combined these data with detailed information, wherever available, about travel history and location (city and district) extracted from announcements of the Health Commissions of the provinces. These data were used to compare transmission of COVID-19 in different provinces and cities. Since the travel history to Wuhan is currently not available for all cases outside of Wuhan, the probability that any confirmed cases were linked with Wuhan was estimated using the 228 confirmed cases with known travel history between 20 January 2020 to 3 February 2020 in Beijing.

2.2. Transmission Analysis

Adopting terminology of Teunis et al. (14) a transmission probability matrix V can be defined where element vi,j is the probability that subject i was infected by another subject j; vi is a vector of transmission probabilities linking case i to any other case. The total number of observed subjects is n. All observed cases could appear twice: as a descendant (infected subject) i and an ancestor (infecting subject) j. Since any case can only have one ancestor the corresponding network must be sparse. Elements of V can be estimated by utilizing a distance kernel κi,j(Xi,j|ij), that defines a pairwise likelihood that subject i was infected by subject j. The distribution of the serial interval (Figure 1A): the distance in time between pairs of cases defines a practical distance kernel, translating the time intervals between symptom onsets in any two cases into a likelihood that these cases were linked as a transmission pair (19).

FIGURE 1
www.frontiersin.org

Figure 1. The serial interval between linked cases. (A) when a subject is infected, there is a latency until they start shedding virus and become infectious to others, then they (may) become symptomatic and can be detected clinically. When subject j infects subject i the time between symptom onset in case j and symptom onset in their descendant case i is the serial interval for onset of clinical symptoms. (B) Distribution of the serial interval for onset of clinical symptoms, estimated from COVID-19 clusters with partially known links (Figure 2). Shown are the prior and the posterior mode distribution (the most probable), with 90% predictive interval.

As outlined in Teunis et al. (14), the elements of the transmission probability matrix V may be estimated in a Markov chain Monte Carlo procedure. The elements of the transmission matrix are subject to constraints. Diagonal elements must be zero (subjects cannot infect themselves) and rows must add to 1 (unless the parent of the corresponding subject was not observed). Additional constraints may be imposed, for instance by preventing links between subjects known to have not been in contact. A mask matrix M (n × n, like V) may be defined, with elements 1 where links are admissible, and 0 where they are not. This mask M may be applied to the matrix of kernels κi,j() by elementwise multiplication. Elements of V representing pairs of subjects with inadmissible links are thus excluded: the corresponding vi,j are set to zero, and they are not updated in MCMC estimation. The mask M can be used to define contacts: whenever subject i is known to have been in contact with another subject j element mi,j = 1. And perhaps more importantly: when subject i is known to not have had contact with subject k then mi,k = 0.

When a sufficient number of infectious contacts is known, the serial interval distribution may be estimated from outbreak data. First, the elements of the transmission probability matrix V are estimated, using (plausible) starting values for the parameters of the serial interval distribution. Then, V is fixed and the serial interval distribution parameters are estimated. Then the serial interval distribution is fixed and V is estimated. This alternating procedure can be repeated until no more improvement (in posterior probability) is found (14).

In the present analyses a special node (0) was defined, that has a uniform kernel κi,0: any subject can be infected by node 0 at any time, within a given time range. Outside that time range κi,0 = 0. Such a node could represent environmental transmission, i.e., from “the Market” to anyone in contact with that environment, or contact with a pool of infectious subjects, i.e., from any infectious subjects within Wuhan/Hubei to subjects outside of Wuhan/Hubei.

As the outbreak progressed, travel to Wuhan became increasingly less likely, especially after the lockdown of the city of Wuhan on 23 January 2020. However, contacts did not cease abruptly: travel records from the confirmed cases in Beijing showed a gradually decreasing probability of contact with Wuhan. Therefore the probability of linking any confirmed cases outside Hubei province to Wuhan was modeled as a logistic function of the date of confirmation

P(linked with Wuhan|t)=eβ0+β1t1+eβ0+β1t    (1)

where t is the number of days after 20 January 2020. The parameters were estimated using the travel history records from the 228 confirmed cases in Beijing (Figure A1). The fitted logistic relation was then used to impute contacts with Wuhan sources for all cases outside Hubei province with unknown travel history.

Once infection has been imported from Wuhan, further transmission depends on local contacts: Since these comprised relatively few cases it is unlikely that all local cases were equally connected. To simulate inhomogeneous transmission local contacts were imputed by distributing cases within a province or city into clusters. Cases within a cluster were assumed fully connected while cases in different clusters were not connected, thus allowing only transmission within a cluster and not between clusters. Cluster size was assumed random, with average size five people, and a high probability (88.8%) of small clusters (average 1–5 people) and a low probability (11.2%) of medium to large clusters (average 10–25 people).

The transmission probability matrix may be used to estimate reproduction numbers by calculating the row sums

Rj=i=1nvi,j    (2)

that represent expected numbers of cases infected by infectious case j (its outdegree), or its reproduction number. Reproduction numbers thus are calculated for individual subjects. Therefore, they can be grouped in any way useful: by time, to illustrate the progression of transmission during the outbreak; by location, to illustrate spatial variation in transmission; or by other characteristics, like age or gender.

In order to translate the transmission probability matrix into a transmission network a multinomial sample of size 1 and probability vector row vi reduces the vector vi to a binary vector with exactly one element 1 and the rest zero. This reduced matrix may be interpreted as an adjacency matrix, equivalent to a directed graph representing a transmission tree of the outbreak (14).

3. Results

For 112 confirmed cases in Tianjin between 21 January 2020 and 12 February 2020, the transmission network could be estimated based on the dates of symptom onset, augmented by contact information between cases. Table 1 shows some summary statistics of demographic and contact information for those cases in Tianjin. As contacts between many of the cases could be identified, so that many elements of the transmission probability matrix were known, joint estimation of the remaining unknown contact probabilities and the serial interval distribution was feasible. With a prior for the serial interval distribution set as Gamma(4,2), 1,000 updates were performed, alternating between probability transmission matrix and serial interval distribution parameters, with 100 iterations for each update and the best fit parameters were chosen as those with the highest posterior probability (14). Figure 2 shows an estimated (posterior mode) transmission network. Using the available contact information, there appeared to be considerable variation in the sizes of clusters of cases. There were two major clusters: one among crew members of Tianjin railway and another one among sales representatives and customers in Baodi department store. Figure 1B shows the posterior mode serial interval distribution as a gamma distribution (shape parameter 3.16, scale parameter 1.52), with a mean of 4.8 days. This best fit serial interval distribution has been used in all following analyses.

FIGURE 2
www.frontiersin.org

Figure 2. (Posterior mode) Transmission tree for Tianjin between 21 January 2020 and 12 February 2020. The nodes are positioned horizontally corresponding to the x-axis representing their date of symptom onset. Vertical positioning is arbitrary (adjusted to prevent crossing of links). Most initial (founding) cases of clusters in Tianjin were imported from Wuhan/Hubei. There were two clusters of cases that could be identified: one among crew members of Tianjin railway and another one among sales representatives and customers in Baodi department store. Both clusters are marked with a box. Many small clusters were among relatives and friends (5 of size 4; 3 of size 3; 5 of size 2) and a proportion of imported cases did not infect anybody (24 of size 1).

Figure 3 shows estimated reproduction numbers by day for the initial 425 confirmed cases in Wuhan. “The market” node as an infectious source was linked to as many as 13.6 cases, on average. For those 425 initial cases, the mean reproduction number was 2.5 (until 31 December 2019).

FIGURE 3
www.frontiersin.org

Figure 3. Reproduction number estimates for the period 8 December 2019 to 21 January 2020 in Wuhan, Hubei province. The x-axis represents the date of symptom onset. The y-axis in the top plot shows the number of cases while the y-axis in the bottom plot shows the reproduction number. The dot and whiskers represent the mean and 2.5th–97.5th percentile range of reproduction number on a specific day, respectively.

The estimated numbers of imported cases from Wuhan varied among provinces in mainland China. Provinces bordering Hubei province include Henan, Hunan, Anhui, and Jiangxi (213, 175, 163, 146), and provinces with close economic ties with Hubei province include Zhejiang and Guangdong (271 and 228, respectively). These provinces all had high estimated numbers of cases imported from Wuhan by 3 February 2020 (Figure 4B).

FIGURE 4
www.frontiersin.org

Figure 4. Transmission in provinces in China: (A) average reproduction numbers and (B) estimated numbers of cases imported (linked to) Wuhan.

Outside Hubei there was less transmission than within Hubei province, with average reproduction numbers between 1 and 2 during the first 5 days after the first confirmed local case. The provinces Guizhou, Hebei, Shandong, and Guangdong had higher average reproduction numbers: 2.25, 2.14, 1.95, and 1.91, respectively (Figure 4A). Though the average reproduction numbers were not very high, there existed large variations in reproduction numbers of individual cases, which indicated some individual cases could have spread the disease to many people (Figure 2).

Within a province, there was also spatial variation in the amount of transmission. Figure 5 shows the disease transmission in 21 cities within Guangdong province. Shenzhen and Guangdong, the two largest and most affluent urban centers in Guangdong, had the highest estimated numbers of cases imported from Wuhan (127 and 124, respectively) by 3 February 2020. Guangzhou, Zhaoqing, Shaoguan, Shenzhen, and Zhuhai had reproduction numbers slightly larger than 1 (1.35, 1.31, 1.30, 1.10, and 1.05, respectively) for first 5 days after the first confirmation of a local case in the city.

FIGURE 5
www.frontiersin.org

Figure 5. Transmission in cities in Guangdong province: (A) average reproduction numbers and (B) estimated numbers of cases imported (linked to) Wuhan.

Figures A2, A3 show boxplots of estimated reproduction numbers for provinces in China and cities in Guandong province, respectively. Figure A4 shows reproduction number estimates by gender and age for Beijing and Tianjin.

4. Discussion

Spreading patterns of COVID-19 at different spatial scales, from local (railway department and department store) to regional (cities within Guangdong Province) to nationwide (Provinces within China), appear similar: new cases arise from contacts with a reservoir, spreading to remote locations through travel from the origin outbreak region. There likely is some secondary spread locally, from cases infected from the reservoir, but local reproduction numbers are low, often insufficient to support sustained transmission.

Knowledge of the serial interval, for symptomatic cases in the outbreak, is essential for analysis of the transmission probabilities (14). Previous studies have estimated the distribution of serial intervals from a set of confirmed transmission pairs, where both the ancestor (who caused infection) and the descendant (who became infected) were known (19, 20). As such information is mostly lacking (6) (and may be hard to obtain) in this COVID-19 outbreak, we have estimated the serial interval distribution from a curated network, using whatever contact information was available. Existing knowledge of contacts (and timing of those contacts) between cases was used to construct a matrix of prior probabilities, to restrict the transmission network to only those links that are possible.

When all transmission links are known, the serial interval distribution can be easily estimated. However, at early stages of an outbreak, information on transmission links is usually incomplete. When a sufficiently large proportion of the transmission links is known, it is possible to estimate both the serial interval distribution, and the transmission probability matrix (14). Starting from a plausible serial interval distribution, the transmission probability matrix V is updated until convergence; then V is frozen and the serial interval parameters are updated, again until obtaining a new optimum; then the serial interval distribution is frozen and the transmission probability matrix is updated again, and so on, until no further improvement can be found. The estimated serial interval distribution, Gamma(3.16, 1.52), has a mean of 4.8 days and good spread out (i.e., large scale parameter) to long serial interval. The variation of serial interval could caused by highly varied incubation period (21). Such a short serial interval, compared to SARS (mean: 8.4 days) (22) and MERS (mean: 6.8 days) (23), gives COVID-19 ability to spread more rapidly. The rapid spread of COVID-19 in South Korea (from 31 cases on 18 February to 2,022 cases on 28 February) and Italy (from 20 cases on 21 February to 650 cases on 28 February), shows how missed infectious subjects may cause rapid transmission within a very short period, due to the combination of a short serial interval and an occasionally high reproduction number (24).

As the incubation period seems to be highly variable, it may be possible that appearance of symptoms in any case precedes symptom onset in its ancestor. When that happens, the serial interval for symptom onset is negative. Details of the negative serial interval can be found in the Supplementary Material. To check whether negative serial intervals would adversely affect analysis we used an alternative distribution, where the serial intervals were shifted leftward by a small amount. A shift of 1 or 2 days had no destructive effect on estimation of V, and the resulting estimates of the effective reproduction numbers did not change substantially.

Early transmission of a newly emerging infectious disease in a population lacking immunity could reveal the basic reproduction number as long as there is no intervention. Currently, the source of COVID-19 is still unclear but many of the early cases were reported to have had contact with the Huanan Seafood Wholesale Market. A recent announcement released by the Chinese Academy of Sciences indicated that the source may not have come from this Seafood Market (25). Any patient zero would have occurred before the initial 425 confirmed cases, but lacking any specific information about early spread of COVID-19, it can be assumed that a reservoir of infectious subjects could function as a super–spreading node causing new infections upon contact at any time within the outbreak. Assuming such a source for the initial cluster in Wuhan (“the Market”) led to an estimated basic reproduction number 13.6. It must be noted that this environmental “super–spreading node” is hypothetical, and represents the production of new cases by the joint presence of infectious subjects connected to the initial source.

There was strong heterogeneity in disease transmission in different provinces of mainland China. The present analysis produces two relevant characteristics: the numbers of cases imported from Wuhan, and the average reproduction number for the first 5 days after any confirmed imported case in each province. When there are many imported cases from Wuhan, these constitute a large base number of cases to start local spreading of the disease. Provinces geographically close to Hubei province are expected to experience mass migration from Wuhan, especially near the holiday season. Provinces with big tier–1 cities like Zhejiang and Guangdong both have very high estimated numbers of cases imported from Wuhan due to their economic connection with Hubei. When early transmission has a high reproduction number, local growth rate in a province will be high. Provinces geographically distant from Hubei, like Hebei, Guizhou, and Shandong, appear to have relatively high reproduction numbers during early local transmission. In provinces distant from Wuhan, the central hub of the outbreak, inhabitants and local governments may have been less cautious (with few public and personal prevention measures) considering high costs of interventions and perceived low risk.

Since Guangdong province has both large numbers of imported cases and high reproduction numbers for early transmission, we further examined the transmission inside of Guangdong by city. Shenzhen and Guangzhou, two of the biggest cities in China, had many imported cases. This is expected considering intense economic connectivity: transportation links, for goods and people including enormous numbers of migrant workers working in those two cities. The numbers of cases imported from Wuhan into other cities in Guangdong province was smaller. This pattern matches the flow of people in public transport (by train and airplane) converging onto the main transportation hubs (Shenzen and Guangzhou) and diverging from there toward other destinations in Guangdong province. For early transmission following import from Wuhan, the reproduction number in most cities in Guangdong had an average reproduction number around 1, during the first 5 days. With appropriate and sustained disease prevention and control measures (e.g., suspension of public transportation, cancellation of mass gatherings, implementation of surveillance and testing, and promotion of personal hygiene and face mask use), the outbreak is unlikely to spread out in those cities.

Among reported cases, gender differences seemed unimportant; most serious illnesses occurred among the elderly, in particular those with health problems prior to infection, while only a small proportion of cases were young of age (21). In the present study, when reproduction number estimates were grouped by gender, no differences were found. Similarly, estimates of transmission from different age groups showed that small children and elderly people were equally likely to transmit infection as any other age group (Figure A4). We also examined the effect of the size of clusters imputed on the results and the difference was trivial between average cluster size three and five.

An important issue in analyzing transmission of COVID-19 is the amount of silent transmission. As mentioned earlier, some infectious subjects may transmit their infection before they become symptomatic. As long as their descendant cases also go on to develop symptoms their transmission link may still be established. However, when infected subjects who remain completely asymptomatic could be infectious to susceptible contacts, the appearance of such contacts, when symptomatic, could not be linked to their immediate ancestors. Given the fact that symptoms seem to be milder in those who are young, such unobservable transmission cannot be excluded (26). Silent transmission has been seen in other infectious diseases (27) where, notably, asymptomatically infected subjects appeared to cause fewer transmission. As asymptomatic infections would lead to an antibody response, serology could be a valuable tool to assess the importance of asymptomatic transmission.

Data Availability Statement

The datasets analyzed for this study can be found at Github (https://github.com/ywan446/covid-19-transmission).

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We were grateful to the alumni of cohort B5 '08 at Shijiazhuang No. 2 Middle School, including Di An, Jiayi Liu, Xun Li, Menghan Hou, Yunpeng Zhang, Fan Huang, and Hongxu Li for gathering and sharing information form various sources about COVID-19.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2020.00329/full#supplementary-material

References

1. Gorbalenya AE, Baker SC, Baric RS, de Groot RJ, Drosten C, Gulyaeva AA, et al. Severe acute respiratory syndrome-related coronavirus: the species and its viruses a statement of the coronavirus study group. bioRxiv. (2020). doi: 10.1101/2020.02.07.937862

CrossRef Full Text | Google Scholar

2. World Health Organization. Novel Coronavirus (2019-nCoV) Situation Report 22: 11 February 2020. Geneva: World Health Organization (2020).

Google Scholar

3. World Health Organization. Coronavirus Disease (COVID-19) Situation Report 129: 28 May 2020. Geneva: World Health Organization (2020).

Google Scholar

4. National Health Commission of the People's Republic of China (2020). Available online at: http://www.nhc.gov.cn/xcs/yqtb/202002/4dcfcb9b74ea4a408fc1d56d4db61f93.shtml (accessed February 12, 2020).

5. Hubei Government (2020). Available online at: http://wjw.hubei.gov.cn/fbjd/dtyw/202002/t20200220_2141274.shtml (accessed February 12, 2020).

6. Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y, et al. Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia. N Engl J Med. (2020) 382:1199–207. doi: 10.1056/NEJMoa2001316

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. (2020) 385:497–506. doi: 10.1016/S0140-6736(20)30183-5

PubMed Abstract | CrossRef Full Text | Google Scholar

8. CNN (2020). Available online at: https://www.cnn.com/travel/article/chunyun-2020-lunar-new-year-travel-rush-china/index.html (accessed February 12, 2020).

9. Health Commission of Guangdong Province (2020). Available online at: http://wsjkw.gd.gov.cn/zwyw_yqxx/content/post_2876057.html (accessed February 12, 2020).

10. National Health Commission of the People's Republic of China (2020). Available online at: http://www.nhc.gov.cn/xcs/yqtb/202001/c5da49c4c5bf4bcfb320ec2036480627.shtml (accessed February 12, 2020).

11. World Health Organization. Novel Coronavirus (2019-nCoV) Situation Report 4: 24 January 2020. Geneva: World Health Organization (2020).

Google Scholar

12. Hubei Government (2020). Available online at: http://www.hubei.gov.cn/zhuanti/2020/gzxxgzbd/zxtb/202001/t20200123_2014402.shtml (accessed February 12, 2020).

13. BBC (2020). Available online at: https://www.bbc.com/news/world-asia-china-51217455 (accessed February 12, 2020).

14. Teunis P, Heijne JCM, Sukhrie F, van Eijkeren J, Koopmans M, Kretzschmar M. Infectious disease transmission as a forensic problem: who infected whom? J R Soc Interface. (2013) 10:20120955. doi: 10.1098/rsif.2012.0955

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Tianjin Health Commission (2020). Available online at: http://www.tj.gov.cn/ (accessed February 12, 2020).

Google Scholar

16. Tianjin Centers for Disease Control and Prevention (2020). Available online at: http://www.cdctj.com.cn/ (accessed February 12, 2020).

17. Weibo (2020). Available online at: https://www.weibo.com/p/1002061960785875/home?from=page_100206&mod=TAB&is_all=1#place (accessed February 12, 2020).

18. Wikipedia Mainland China (2020). Available online at: https://zh.wikipedia.org/wiki/2019%E5%86%A0%E7%8B%80%E7%97%85%E6%AF%92%E7%97%85%E4%B8%AD%E5%9C%8B%E5%A4%A7%E9%99%B8%E7%97%85%E4%BE%8B (accessed February 12, 2020).

19. Wallinga J, Teunis P. Different epidemic curves for Severe Acute Respiratory Syndrome reveal similar impacts of control measures. Am J Epidemiol. (2004) 160:509–16. doi: 10.1093/aje/kwh255

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Heijne JCM, Teunis PFM, Morroy G, Wijkmans C, Oostveen S, Duizer E, et al. Enhanced hygiene measures and norovirus transmission during an outbreak. Emerg Infect Dis. (2009) 15:24–30. doi: 10.3201/eid1501.080299

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Guan W, Ni Z, Hu Y, Liang W, Ou C, He J, et al. Clinical characteristics of 2019 novel coronavirus infection in China. MedRxiv. (2020). doi: 10.1101/2020.02.06.20020974

CrossRef Full Text | Google Scholar

22. Lipsitch M, Cohen T, Cooper B, Robins JM, Ma S, James L, et al. Transmission dynamics and control of severe acute respiratory syndrome. Science. (2003) 300:1966–70. doi: 10.1126/science.1086616

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Cauchemez S, Nouvellet P, Cori A, Jombart T, Garske T, Clapham H, et al. Unraveling the drivers of MERS-CoV transmission. Proc Natl Acad Sci USA. (2016) 113:9081–6. doi: 10.1073/pnas.1519235113

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Wallinga J, Lipsitch M. How generation intervals shape the relationship between growth rates and reproductive numbers. Proc R Soc B Biol Sci. (2007) 274:599–604. doi: 10.1098/rspb.2006.3754

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Yu WB, Tang GD, Zhang L, Corlett RT. Decoding the evolution and transmissions of the novel pneumonia coronavirus (SARS-CoV-2) using whole genomic data. chinaXiv. (2020) chinaXiv:202002.00033v2.

PubMed Abstract | Google Scholar

26. Zhang W, Du RH, Li B, Zheng XS, Yang XL, Hu B, et al. Molecular and serological investigation of 2019-nCoV infected patients: implication of multiple shedding routes. Emerg Microbes Infect. (2020) 9:386–9. doi: 10.1080/22221751.2020.1729071

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Sukhrie FHA, Teunis PFM, Beersma MFC, Copra C, Vennema H, Bogerman J, et al. Nosocomial transmission of norovirus is mainly caused by symptomatic cases. Clin Infect Dis. (2012) 54:931–7. doi: 10.1093/cid/cir971

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: novel coronavirus disease, COVID-19, transmission, serial interval, reproduction number

Citation: Wang Y and Teunis P (2020) Strongly Heterogeneous Transmission of COVID-19 in Mainland China: Local and Regional Variation. Front. Med. 7:329. doi: 10.3389/fmed.2020.00329

Received: 26 March 2020; Accepted: 04 June 2020;
Published: 19 June 2020.

Edited by:

Zhongheng Zhang, Sir Run Run Shaw Hospital, China

Reviewed by:

Thiago DeSouza-Vieira, National Institute of Allergy and Infectious Diseases, National Institutes of Health (NIH), United States
Hong-Qiu Gu, National Clinical Research Center for Neurological Diseases, China
Lu Ke, Medical School of Nanjing University, China
Adriano Barra, University of Salento, Italy

Copyright © 2020 Wang and Teunis. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yuke Wang, yuke.wang@emory.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.