Hostname: page-component-586b7cd67f-tf8b9 Total loading time: 0 Render date: 2024-11-22T17:49:30.888Z Has data issue: false hasContentIssue false

Intercontinental transmission and local demographic expansion of SARS-CoV-2

Published online by Cambridge University Press:  13 April 2021

Hong-yin Hu
Affiliation:
State Key Laboratory of Grassland Agro-Ecosystem, Institute of Innovation Ecology & School of Life Sciences, Lanzhou University, Lanzhou, China
Fang Yan
Affiliation:
College of Life Science, Yunnan University, Kunming650500, China
Jia-ming Zhu
Affiliation:
State Key Laboratory of Grassland Agro-Ecosystem, Institute of Innovation Ecology & School of Life Sciences, Lanzhou University, Lanzhou, China
Alex Plimo Karuno
Affiliation:
State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, 650223, China
Wei-wei Zhou*
Affiliation:
State Key Laboratory of Grassland Agro-Ecosystem, Institute of Innovation Ecology & School of Life Sciences, Lanzhou University, Lanzhou, China State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, 650223, China
*
Author for correspondence: Wei-wei Zhou, E-mail: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

The global outbreak of coronavirus disease 2019 (COVID-19) is greatly threatening the public health in the world. We reconstructed global transmissions and potential demographic expansions of severe acute respiratory syndrome coronavirus 2 based on genomic information. We found that intercontinental transmissions were rare in January and early February but drastically increased since late February. After world-wide implements of travel restrictions, the transmission frequencies decreased to a low level in April. We identified a total of 88 potential demographic expansions over the world based on the star-radiative networks and 75 of them were found in Europe and North America. The expansion numbers peaked in March and quickly dropped since April. These findings are highly concordant with epidemic reports and modelling results and highlight the significance of quarantine validity on the global spread of COVID-19. Our analyses indicate that the travel restrictions and social distancing measures are effective in containing the spread of COVID-19.

Type
Original Paper
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Author(s), 2021. Published by Cambridge University Press

Introduction

The outbreak of coronavirus disease 2019 (COVID-19) [Reference Wang1] caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) [Reference Gonzalez-Reiche2] has been and is still threatening global public health. The origin of SARS-CoV-2 is still under debate. It is difficult to find the root cause of the disease and different approaches yield contrasting results [Reference Pipes3]. Meanwhile, two studies based on network analyses identified ancestral haplogroup by using sequences of coronavirus from bats as an outgroup [Reference Forster4, Reference Yu5]. However, the reliability of this result is under debate [Reference Sánchez-Pacheco6Reference Mavian8]. After an outbreak in east Asia, COVID-19 dramatically spread worldwide. The first cases were confirmed in Europe, North America and Oceania in January 2020. In February, the cases were officially reported in all continents. Until May 2020, the officially confirmed cases were reported in more than 200 countries. Up to 1 July 2020, the world-wide outbreak of COVID-19 had caused more than 510 thousand deaths (https://coronavirus.jhu.edu/map.html). Studies based on epidemic models evaluated the transmissions and expansions of SARS-CoV-2 in the world [Reference Chinazzi9Reference Siedner14], suggesting the likely complex demographic dynamics. However, the direct evidence from population genomic data is lacking. This is critical to trace transmission events and identify community-based outbreaks based on the demographic expansions [Reference Croucher and Didelot15, Reference Volz, Koelle and Bedford16]. After the outbreak of COVID-19, a large number of viral genomes were sequenced and uploaded to the public database (https://www.gisaid.org/). These open resources provide us a chance to examine evolutionary history of the SARS-CoV-2. We can trace the intercontinental transmission and potential demographic expansions against timescales and accompanying travel restrictions in the world. In this study, we evaluated the impact of travel restrictions on the spread of COVID-19 based on genomic information. Meanwhile, we checked if the strong social distancing measures are effective in reducing local demographic expansions of virus.

Materials and methods

Sequences selection

To reconstruct the evolutionary history of SARS-CoV-2, we downloaded 31 820 full-length genome sequences on 26 May 2020. The sequences were downloaded from the Global Initiative on Sharing Avian Influenza Data (GISAID) (Table S1). We also collected the sampling time of the sequences. As the random sequencing errors could generate a large amount of single site mutation, which would affect the accurate definition of the haplotype and network, we removed all low-quality sequences following the reference [Reference Zhao17]. We filtered the data based on mutation density and numbers of unknown bases (N), degenerate bases and gaps. Mutation density was defined by mutations/slide window size (20nt). The gaps were identified by comparing with reference sequence (MN908947). The unknown bases should be less than 15, the degenerate bases should be less than 50 and the gaps should be less than three. If a sequence contained high mutation density regions (mutation density larger than or equal to 0.25), it would also be treated as low quality.

Construction of phylogenetic tree and network

The genome sequences were aligned using MAFFT [Reference Kuraku18]. Haplotypes were determined using the R package ‘pegas’ [Reference Paradis19]. We used NETWORK 4.10 to build a median-joining network [Reference Bandelt, Forster and Röhl20]. The phylogenetic tree was constructed by using IQtree [Reference Nguyen21]. Support values were estimated by using ultrafast bootstrap approximation [Reference Hoang22, Reference Minh, Nguyen and Von Haeseler23] with 1000 bootstrap replicates. Haplotype diversity and nucleotide diversity were calculated using the R package ‘pegas’ [Reference Paradis19].

Transmission pattern reconstruction

Based on the sampling date of the genome sequences, we obtained the potential origin time for each haplotype. We set the earliest sampling date as origin time of the haplotype and the sampling dates of sequences in other continents as transmission time. We corrected origin time and place of haplotypes based on the epidemic information, especially the travelling history of the first few confirmed cases of each country. For some haplotypes, the first reported sequences were isolated from patients who had an international travelling history. Therefore, we corrected the origin place as where they may be infected. In total, the origin place of eight haplotypes were corrected (H10, H121, H3196, H2714, H2354, H179, H1565 and H1035). We assumed that the same haplotype or two directly connected haplotypes in the network found in the different continents represented one intercontinental transmission. The links with more than 10 mutations were not involved in the analyses, as it may be caused by sequencing errors.

Local demographic expansion

The rapid population expansion of the virus through transmission will lead to the production of the numerous haplotypes with short step (one or two mutations) from one central haplotype [Reference Bandelt24]. Except seriously sampling biases, the frequency of these haplotypes should be lower than the central haplotype. This local demographic expansion from the potential virus outbreak could be identified by the star-like phylogeny of the haplotypes in the network [Reference Bandelt24]. We set haplotypes that linked with 10 or more haplotypes in the network as central haplotype. Then we extracted the haplotypes directly connected with the central haplotypes. The samples of these haplotypes together with the central haplotype were used to define local demographic expansion because of the outbreaks of SARS-CoV-2. We examined the distributions of each expansion in the continent and the total occurrence frequency of expansions. Analyses were repeated under different thresholds to see if we can get similar patterns. We also assigned sequences to clades following the approach of Rambaut et al. [Reference Rambaut25]. Then we checked the frequency of outbreaks for each clade.

Results

Genome sequence variation and network construction

In total, we downloaded 31 820 full-length genome sequences and 14 206 sequences passed the quality controls. After trimming the ends, the alignment we used was 29 182 bps. The alignment included 6275 variable sites, in which 2757 are potential phylogenetic informative sites. In addition, we identified a total of 7177 haplotypes from all the reported genomes (Table S1). We constructed the phylogenetic tree (Fig. 1a) and the network based on sequences of haplotypes. The overall nucleotide diversity was 0.000098 and the haplotype diversity was 0.996.

Fig. 1. (a) A phylogeny tree based on SARS-CoV-2 genomes by IQtree. Haplotypes associated with the local demographic expansions identified by network analyses are indicated in red. (b) Intercontinental transmissions against timescale. x-axis stands for time and y-axis stands for number of transmissions. (c) The illustrated intercontinental transmissions between different continents and the line thickness correspond to transmission frequencies.

Global transmission pattern of COVID-19

When a haplotype or two directly connected haplotypes on the network were found in different continents, one intercontinental transmission was defined. In total, we identified 412 intercontinental transmissions (Table S2). The transmissions occurred rarely in January and February, but drastically increased since the end of February (Fig. 1b). The transmission frequencies peaked in mid-March, which was about 10 times higher than those in January. Since April, the transmission frequencies dropped down to a low level as in January when travel restrictions were implemented in the world.

Local demographic expansion

Rapid population expansion because of frequent transmissions of virus will lead to the fixations of rare and random mutation due to founder effect [Reference Croucher and Didelot15, Reference Campbell26]. This will result in the star-phylogeny network, in which rare haplotypes with short mutation are connected to a central haplotype with high frequency [Reference Bandelt24]. We defined a subnetwork with 10 or more haplotypes with few mutations connected to one high-frequency central haplotype as one potential expansion. In total, we identified 88 potential demographic expansions of COVID-19 (Table S3 and Fig. 2a) by this approach. Temporal pattern of the expansion, which yielded from the earliest sampling time of central haplotypes, indicated that most expansions emerged since the end of February. Six expansions occurred before February, 18 in February, 61 in March and 2 since April. For the 18 ones in February, 16 of them took place after 20 February. Following the clade assignment, we identified 23 expansions in clade 19A, 13 ones in clade 19B, 25 ones in clade 20A, 10 ones in clade 20B and 17 ones in clade 20C (Tables S3 and S4). The expansion frequencies may change when different thresholds were used to define each expansion, but the total trend remains similar (Table S5).

Fig. 2. (a) The demographic expansions against timescale. x-axis indicates time and y-axis shows the number of expansions. (b) Examples of the star-radiative structures in the median-joining subnetwork of the SARS-CoV-2 genomes. Colours stand for haplotype frequency.

Discussion and conclusion

Based on the 14 206 sequences that passed the filtering, we identified 412 intercontinental transmissions. Our analyses suggested that the intercontinental transmissions occurred at a lower frequency in January and February (Fig. 1b). However, since the end of February, such intercontinental transmissions increased drastically at a high frequency (Fig. 1b). Most identified transmissions are related with Europe and North America. The intercontinental transmissions started to decrease since mid-March, which is roughly concordant with the time that large-scale travel bans started to be implemented in the world. The frequency returned to the low level since April. Undoubtedly these transmissions only comprised a part of the actual transmissions. Considering the incubation period, the transmission times we estimated may be later than the real time, but the general pattern of transmissions was not impacted by such bias. It should be noted that although travel bans reduced the transmission frequencies, these policies failed to prevent the continuing epidemic of COVID-19. A recent study presented similar patterns [Reference Rito27]. The European clade originated from east Asian clade. However, the subsequent worldwide spread since March 2020 were mainly related to the European clade [Reference Rito27]. As suggested by the epidemic model fitting [Reference Chinazzi9], very few transmissions may have led to severe epidemics and secondary outbreaks when social distancing measures or strict community quarantine are lacking.

This study further identified 88 potential demographic expansions of SARS-CoV-2. Most expansions were reported in Europe and North America. The temporal pattern indicated that multiple expansions occurred when quarantine policies to reduce social contacts had not been implemented [Reference Thunström28]. The local expansions occurred at a lower frequency as the intercontinental transmissions in January and February (Fig. 2a). However, since the end of February, the frequency of expansions drastically increased and peaked in March. After strong social distancing measures were strictly imposed [Reference Thunström28], the number of expansions quickly dropped (Fig. 2a and Table S3). Only two expansions were identified since April. Similar patterns were observed when applying different thresholds (Table S5). This is concordant with other studies [Reference Courtemanche13] and the trends of daily new cases (https://coronavirus.jhu.edu/map.html). For example, in Europe and North America, daily new cases peaked in mid-March or early April, then the curve reached a plateau or started to drop in most countries (https://coronavirus.jhu.edu/map.html).

Overall, the present study examined intercontinental transmission and local demographic expansion of the SARS-CoV-2 in the world based on a big dataset of the genomes. These findings inferred directly from the genome sequences are largely concordant with the epidemic reports and epidemic modelling [Reference Zhao and Chen11, Reference Hsiang12]. Our network analyses illustrated history and relationships of populations without any prior assumption on demographic histories. These results can be used as guidance in subsequent modelling and other studies. Both sample numbers and sampling bias may affect our estimations. For example, transmissions and expansions could not be estimated based on our methods when very few sequences are available. Dense samplings will identify more transmissions and expansion events. But our conclusions are not impacted by such bias, as the genome sequences were continuously increasing in April and May (Table S1). In this way, the decrease of frequency of international transmissions and number of expansions since April could not be an artefact because of fewer samples. Meanwhile, our results are concordant with other studies [Reference Rito27], which indicate the results are reliable. Despite these caveats, we provided the genomic evidence for quarantine validity on the spread of COVID-19 between continents across the world.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0950268821000777.

Acknowledgements

We owe thanks to the website of GISAID for gathering data in our study.

Author contributions

Hong-yin Hu, Fang Yan and Wei-wei Zhou designed the study. Hong-yin Hu, Jia-ming Zhu and Wei-wei Zhou contributed to data collection. Hong-yin Hu and Alex Plimo Karuno carried out the data curation. Hong-yin Hu and Wei-wei Zhou wrote the paper. All authors read and approved the final version of the manuscript.

Conflict of interest

The authors declare no competing financial interests.

Data

The data that support the findings of this study are available from GISAID (https://www.gisaid.org/).

Footnotes

*

These authors contributed equally to this work.

References

Wang, C et al. (2020) A novel coronavirus outbreak of global health concern. The Lancet 395, 470473.CrossRefGoogle ScholarPubMed
Gonzalez-Reiche, AS et al. (2020) Introductions and early spread of SARS-CoV-2 in the New York city area. Science (New York, N.Y.) 369, 297301.CrossRefGoogle ScholarPubMed
Pipes, L et al. (2021) Assessing uncertainty in the rooting of the SARS-CoV-2 phylogeny. Molecular Biology and Evolution 38, 15371543.CrossRefGoogle ScholarPubMed
Forster, P et al. (2020) Phylogenetic network analysis of SARS-CoV-2 genomes. Proceedings of the National Academy of Sciences of the United States of America 117, 92419243.CrossRefGoogle ScholarPubMed
Yu, WB et al. (2020) Decoding the evolution and transmissions of the novel pneumonia coronavirus (SARS-CoV-2 / HCoV-19) using whole genomic data. Zoological Research 41, 247257.CrossRefGoogle Scholar
Sánchez-Pacheco, SJ et al. (2020) Median-joining network analysis of SARS-CoV-2 genomes is neither phylogenetic nor evolutionary. Proceedings of the National Academy of Sciences of the United States of America 117, 1251812519.CrossRefGoogle ScholarPubMed
Forster, P et al. (2020) Reply to Sánchez-Pacheco, et al., Chookajorn, and Mavian et al.: explaining phylogenetic network analysis of SARS-CoV-2 genomes. Proceedings of the National Academy of Sciences of the United States of America 117, 1252412525.CrossRefGoogle ScholarPubMed
Mavian, C et al. (2020) Sampling bias and incorrect rooting make phylogenetic network tracing of SARS-COV-2 infections unreliable. Proceedings of the National Academy of Sciences of the United States of America 117, 1252212523.CrossRefGoogle ScholarPubMed
Chinazzi, M et al. (2020) The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak. Science 368, 395400.CrossRefGoogle ScholarPubMed
Wells, CR et al. (2020) Impact of international travel and border control measures on the global spread of the novel 2019 coronavirus outbreak. Proceedings of the National Academy of Sciences of the United States of America 117, 75047509.CrossRefGoogle ScholarPubMed
Zhao, S and Chen, H (2020) Modeling the epidemic dynamics and control of COVID-19 outbreak in China. Quantitative Biology 8, 1119.CrossRefGoogle Scholar
Hsiang, S et al. (2020) The effect of large-scale anti-contagion policies on the COVID-19 pandemic. Nature 584, 262267.CrossRefGoogle ScholarPubMed
Courtemanche, C et al. (2020) Strong social distancing measures in the United States reduced the COVID-19 growth rate. Health Affairs 39, 12371246.CrossRefGoogle ScholarPubMed
Siedner, MJ et al. (2020) Social distancing to slow the US COVID-19 epidemic: longitudinal pretest-posttest comparison group study. PLoS Medicine 17, e1003244.CrossRefGoogle ScholarPubMed
Croucher, NJ and Didelot, X (2015) The application of genomics to tracing bacterial pathogen transmission. Current Opinion in Microbiology 23, 6267.CrossRefGoogle ScholarPubMed
Volz, EM, Koelle, K and Bedford, T (2013) Viral phylodynamics. PLoS Computational Biology 9, e1002947.CrossRefGoogle ScholarPubMed
Zhao, WM et al. (2020) The 2019 novel coronavirus resource. Hereditas 42, 212221.Google ScholarPubMed
Kuraku, S et al. (2013) aLeaves facilitates on-demand exploration of metazoan gene family trees on MAFFT sequence alignment server with enhanced interactivity. Nucleic Acids Research 41, W22W28.CrossRefGoogle ScholarPubMed
Paradis, E (2010) Pegas: an R package for population genetics with an integrated-modular approach. Bioinformatics 26, 419420.CrossRefGoogle Scholar
Bandelt, HJ, Forster, P and Röhl, A (1999) Median-joining networks for inferring intraspecific phylogenies. Molecular Biology and Evolution 16, 3748.CrossRefGoogle ScholarPubMed
Nguyen, LT et al. (2015) IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Molecular Biology and Evolution 32, 268274.CrossRefGoogle ScholarPubMed
Hoang, DT et al. (2018) UFBoot2: improving the ultrafast bootstrap approximation. Molecular Biology and Evolution 35, 518522.CrossRefGoogle ScholarPubMed
Minh, BQ, Nguyen, MAT and Von Haeseler, A (2013) Ultrafast approximation for phylogenetic bootstrap. Molecular Biology and Evolution 30, 11881195.CrossRefGoogle ScholarPubMed
Bandelt, HJ et al. (1995) Mitochondrial portraits of human populations using median networks. Genetics 16, 3748.Google Scholar
Rambaut, A et al. (2020) A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology. Nature Microbiology 5, 14031407.CrossRefGoogle ScholarPubMed
Campbell, F et al. (2018) When are pathogen genome sequences informative of transmission events? PLoS Pathogens 14, e1006885.CrossRefGoogle ScholarPubMed
Rito, T et al. (2020) Phylogeography of 27,000 SARS-CoV-2 genomes: Europe as the major source of the COVID-19 pandemic. Microorganisms 8, 1678.CrossRefGoogle ScholarPubMed
Thunström, L, et al. (2020) The benefits and costs of using social distancing to flatten the curve for COVID-19. Journal of Benefit-Cost Analysis 11, 117.CrossRefGoogle Scholar
Figure 0

Fig. 1. (a) A phylogeny tree based on SARS-CoV-2 genomes by IQtree. Haplotypes associated with the local demographic expansions identified by network analyses are indicated in red. (b) Intercontinental transmissions against timescale. x-axis stands for time and y-axis stands for number of transmissions. (c) The illustrated intercontinental transmissions between different continents and the line thickness correspond to transmission frequencies.

Figure 1

Fig. 2. (a) The demographic expansions against timescale. x-axis indicates time and y-axis shows the number of expansions. (b) Examples of the star-radiative structures in the median-joining subnetwork of the SARS-CoV-2 genomes. Colours stand for haplotype frequency.

Supplementary material: File

Hu et al. supplementary material

Hu et al. supplementary material

Download Hu et al. supplementary material(File)
File 922.7 KB