Replacement of the Gamma by the Delta variant in Brazil: impact of lineage displacement on the ongoing pandemic
===============================================================================================================

* Marta Giovanetti
* Vagner Fonseca
* Eduan Wilkinson
* Houriiyah Tegally
* Emmanuel James San
* Christian L. Althaus
* Joilson Xavier
* Svetoslav Nanev Slavov
* Vincent Louis Viala
* Alex Ranieri Jerônimo Lima
* Gabriela Ribeiro
* Jayme A. Souza-Neto
* Heidge Fukumasu
* Luiz Lehmann Coutinho
* Rivaldo Venancio da Cunha
* Carla Freitas
* Carlos F Campelo de A e Melo
* Wildo Navegantes
* Rodrigo Fabiano do Carmo Said
* Maria Almiron
* Tulio de Oliveira
* Sandra Coccuzzo Sampaio
* Maria Carolina Elias
* Dimas Tadeu Covas
* Edward C. Holmes
* Jose Lourenço
* Simone Kashima
* Luiz Carlos Junior de Alcantara

## Abstract

The COVID-19 epidemic in Brazil was driven mainly by the spread of Gamma (P.1), a locally emerged Variant of Concern (VOC) that was first detected in early January 2021. This variant was estimated to be responsible for more than 96% of cases reported between January and June 2021, being associated with increased transmissibility and disease severity, a reduction in neutralization antibodies and effectiveness of treatments or vaccines, as well as diagnostic detection failure. Here we show that, following several importations predominantly from the USA, the Delta variant rapidly replaced Gamma after July 2021. However, in contrast to what was seen in other countries, the rapid spread of Delta did not lead to a large increase in the number of cases and deaths reported in Brazil. We suggest that this was likely due to the relatively successful early vaccination campaign coupled with natural immunity acquired following prior infection with Gamma. Our data reinforces reports of the increased transmissibility of the Delta variant and, considering the increasing concern due to the recently identified Omicron variant, argues for the necessity to strengthen genomic monitoring on a national level to quickly detect and curb the emergence and spread of other VOCs that might threaten global health.

## Main text

Since late 2020 the evolution of SARS-CoV-2 has been characterised by the appearance of mutations on the Spike protein leading to the emergence of so-called ‘variants of concern’ (VOC) some of which have spread globally1,2. In particular, some of the mutations in VOCs are suggested to impact viral transmissibility3, resistance to neutralizing antibodies4, and virulence. The identification of such variants has recently challenged public health authorities with respect to tracking transmission and mitigating the impact in the ongoing pandemic. To date, the most important VOCs documented are Alpha, Beta, Gamma, Delta and Omicron (described in late November 2021), first detected in United Kingdom, South Africa, Brazil, India, and South Africa/Botswana respectively5–8.

In the context of continuous surveillance of SARS-CoV-2 clinical samples on behalf of the National Pandemic Alert Network based in the State of São Paulo (Brazil) since early January 2021, our team has been monitoring the proportion of circulating variants in Brazil. Data from this initiative shows that after July 2021 the Delta variant has become highly prevalent. It has been suggested that the Delta variant might be more transmissible9, having led to public health emergencies in other countries due to overwhelming increases in the number of cases, hospitalizations and deaths when compared to previous circulating variants. Here we describe how Delta became the predominant SARS-CoV-2 variant in Brazil, rapidly replacing the previously dominant Gamma, and how this displacement was not associated with an increase in reported case numbers or deaths.

The COVID-19 epidemic in Brazil can generally be characterized by two epidemic waves accounting for more than 22 million cases and 616,251 deaths until early December 202110. The first epidemic wave was characterised by the circulation of multiple SARS-CoV-2 lineages (among them the B.1.1.28 and the B.1.1.33), as a direct consequence of multiple independent introduction events between February and March 2020 (Fig.1A)11. By the end of October 2020, even with the implementation of non-pharmaceutical interventions (NPIs), a second wave associated with a dramatic resurgence in cases and death numbers took place. This wave was fuelled by the emergence and circulation of several Variants under Monitoring (VUM), such as P.2 (i.e., Zeta), and some VOCs including Alpha and later Gamma (i.e., P.1), which became widespread by January 2021 and dominated the viral population for nearly 8 consecutive months (Figure 1A). Despite the national vaccination rollout beginning on January 17th 2021, the COVID-19 death toll in the country steadily rose in March 2021, reaching a peak in April 2021 (Fig. 1A). This was followed by a decrease in the number of daily cases and deaths by April (Fig. 1A and Fig. S1), likely as a consequence of a gradual increase in population immunity. Zeta (P.2) mostly dominated the first epidemic wave persisting up to March 2021 when it was replaced by Gamma (P.1) (Fig. 1C). This period was characterized by an upsurge in the number of total cases with a peak registered between February and June 2021. The Alpha variant was also detected from January 2021 onwards, but it remained at a very low frequency nationally (less than 2%)12.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/29/2021.12.27.21268309/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2021/12/29/2021.12.27.21268309/F1)

Figure 1. Replacement of Gamma by Delta variant in Brazil.
(A) Dynamics of the SARS-CoV-2 epidemic in Brazil showing the number of daily COVID-19 deaths and the progression in the proportion of circulating variants in the country over time, with the rapid replacement of the Gamma by the Delta variant. (B-C) Modelled proportion of SARS-CoV-2 variants over time in Brazil in a linear and in a logarithmic scales, showing that Gamma became the dominant variant in Brazil by the beginning of 2021 and was rapidly outcompeted by Delta from June 2021. Model fits are based on a multinomial logistic regression. The size of dots corresponds to the weekly sample size. (D) Prevalence maps following the progression in the monthly average of daily number of cases and proportions of relevant variants per region (North, Northeast, Midwest, Southeast and South) in Brazil from October 2020 to October 2021.

In April 2021, through intensified sampling of likely imported cases associated with returning travellers, the Delta variant was detected. Delta first spread in the southeast, the most populous region of Brazil harbouring the largest urban centres and airports with highest national and international travel flows, that continued to operate throughout 202113. Throughout July 2021, Delta rapidly displaced Gamma, becoming the dominant variant circulating nationally during August 2021 (Fig. 1A and Fig. S1). This variant displacement was not associated with a concurrent increase in COVID-19 incidence (Fig. 1C). We estimate that Delta had a growth advantage of 0.064 (95% confidence interval [CI] 0.058-0.071) per day compared to Gamma (Fig. 1B). Assuming that the variants have the same generation time of 5.2 days14, this corresponds to a growth advantage of 33% (95% CI 30-37%) per generation of viral transmission, which is in good agreement with earlier findings based on sequence data from multiple countries15.

Although the precise cause of these trends is unknown, we hypothesize that a relatively successful early vaccination campaign in the region coupled with a reasonable percentage of the population with natural immunity acquired by prior infection with Gamma contributed to the decreasing case rate around the time of introduction of the Delta variant. In contrast, when Gamma dominated transmission, the observed increasing in the number of cases would have been attributed to a relatively naïve population to infection with any SARS-CoV-2 variant. This epidemiological scenario of rapid switching of new variants without an increase in the number of cases seems not to be uncommon and has also been observed elsewhere16–20.

We further estimated phylogenetic trees to explore the relationship of the sequenced Brazilian genome to those of other isolates across the world. For this purpose, we retrieved 11,147 Delta genomes from Brazil, from which 6,626 were generated by our National Pandemic Alert Network and a globally representative set of other Delta genomes (n=13,261). Our time-stamped phylogeny revealed that the Brazilian Delta isolates are scattered throughout the phylogeny suggesting multiple independent introductions (Fig.2A, Figure S2).

Our analysis revealed at least 124 independent introductions of the Delta variant into Brazil between October 2020 and October 2021 (Fig 2b), with the majority (30%) originating from North America, followed by India (17%), European countries (mostly UK 13%), and other South American countries.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/29/2021.12.27.21268309/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2021/12/29/2021.12.27.21268309/F2)

Figure 2. Phylogenetic analysis of the Delta variant in Brazil.
(A) Time-resolved maximum clade credibility phylogeny of 11,147 SARS-CoV-2 Delta sequences from Brazil, from which 6,626 were generated by our National Network for Pandemic Alert of SARS-CoV-2 along with global Delta sequence (grey). Colours indicate different Brazilian regions (North in green; Midwest in orange; Northeast in yellow; Southeast in red and South in light pink); (B) Inferred locations of importations of the Delta variant into Brazil. (C) Inferred locations of exportation of the Delta variant from Brazil.

Following introductions into Brazil, in line with recent findings21, the Delta variant appears to have further spread within Brazilian regions (Fig.2A), highlighting complex local transmission dynamics maintained by travel. Importantly, this variant appears to have been introduced into each region of Brazil through multiple ports of entry making it challenging to accurately reconstruct transmission pathways across the country. Our analysis further shows that the Delta variant subsequently spread from Brazil into countries in South America, and also to North America, Europe and Africa (Fig 2C)22.

The first months of the Brazilian epidemic were fuelled by the circulation of multiple lineages, as a direct consequence of multiple viral independent introductions from overseas. As the epidemic progressed, the observed large-scale community transmission led to local emergence of VOCs and VUMs23. These variants dominated the end of the first and second epidemic wave in the country, resulting in an exponential increase in the number of daily cases and deaths, making Brazil one of the countries hardest hit by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic worldwide.

In this study we analysed the displacement of Gamma by the Delta variant in Brazil, showing also how this replacement was not associated with an increase in reported cases numbers. Given that during the displacement NPIs remained relatively relaxed locally, it seems likely that a significant proportion of the Brazilian population has developed immunity (either through natural infection or vaccination) which helped to prevent a rise in case numbers. We estimated that in Brazil, Delta had a transmission advantage of 33% (95% CI 30-37%) compared to Gamma. The spread of future variants will be possible if they exhibit increased transmissibility or immunity evasiveness, both of which seem likely in the case of the Omicron variant8. However, as shown here, displacement does not necessarily equate to new epidemic waves of reported cases or deaths, since the outcome of the spread of a new variant will depend on a complex interplay between local interventions, infrastructure and the immunity landscape already present in the population due to the circulation of previous variants. Considering that past replacement events have taken place so rapidly in Brazil, future epidemiological assessment of new VOC must be conducted rapidly and regularly. The capacity to detect and respond to new variants requires continued support and funding for molecular surveillance and sequencing capacity more generally. Moreover, enhanced sampling efforts in countries like Brazil are needed to ensure better geographical representativeness of available SARS-CoV-2 sequences and the rapid detection of emerging variants when their frequencies are still low. These factors will be key to detect, understand and respond to the likely upcoming Omicron wave which has already been detected in more than 87 countries, including Brazil, and which will possibly alter the landscape of variants currently circulating in the country.

## Supporting information

Supplementary\_Figure\_1 [[supplements/268309_file04.pdf]](pending:yes)

Supplementary\_Figure\_2 [[supplements/268309_file05.pdf]](pending:yes)

Supplementary\_Table\_1 [[supplements/268309_file06.xlsx]](pending:yes)

## Data Availability

All SARS-CoV-2 whole genome sequences produced by the National Network for Pandemic Alert of SARS-CoV-2 have been deposited in the GISAID sequence database and are publicly available subject to the terms and conditions of the GISAID database. The GISAID accession numbers of sequences used in the phylogenetic analysis, are provided in the Supplementary Table S1. [https://www.gisaid.org/](https://www.gisaid.org/)

[https://www.gisaid.org/](https://www.gisaid.org/) 

## Data availability

All SARS-CoV-2 whole genome sequences produced by the National Network for Pandemic Alert of SARS-CoV-2 have been deposited in the GISAID sequence database and are publicly available sub- ject to the terms and conditions of the GISAID database. The GISAID accession numbers of sequenc- es used in the phylogenetic analysis, are provided in the **Supplementary Table S1**.

## Code availability

All input files along with all resulting output files and scripts used in the present study will be made available upon request and publicly shared on GitHub at final publication.

## Funding

This work was supported in part through National Institutes of Health USA grant U01 AI151698 for the United World Arbovirus Research Network (UWARN), the CRP- ICGEB RESEARCH GRANT 2020 Project CRP/BRA20-03, Contract CRP/20/03, the Oswaldo Cruz Foundation VPGDI-027-FIO- 20-2-2-30, the Brazilian Ministry of Health (SCON2021-00180), the CNPq (426559/2018-5) and the Faperj E-26/202.930/2016. MG and LCJA are supported by Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ) (E-26/202.248/2018(238504); E26/202.665/2019(247400)). CA received funding from the European Union’s Horizon 2020 research and innovation programme – project EpiPose (No 101003688). Research reported in this publication/article was also supported by the South African Department of Science and Innovation (DSI) and the South African Medical Research Council (SAMRC) under the BRICS JAF #2020/049. The content and findings reported/illustrated herein are the sole deduction, view and responsibility of the researcher/s and do not reflect the official position and sentiments of the funders.

## Competing interests statement

The authors declare no competing interests.

## Author Contributions

Molecular screening and produced SARS-CoV-2 genomic data: MG, JX, SNS, VLV, ARJ, GR, JASN, HF, LLC, SCS, MCE, and SK; Collected samples and curated metadata: SCS, MCE, VF, DTC, SK, LCJA, CF, CFCAM, WN, RFCS, MA; Analysed the data: MG, VF, EW, HT, EJS, CLA; Helped with study design and data interpretation: MG, VF, EW, HT, EJS, CLA;, ECH, JL, and LCJA; Wrote the initial manuscript, which was reviewed by all authors: MG, CLA, EH, JL and LCJA.

## Materials and Methods

### Ethics statement

This research was approved by the Ethics Review Committee of the Pan American Health Organization (PAHOERC.0344.01) and by the Federal University of Minas Gerais (CEP/CAAE: 32912820.6.1001.5149). The availability of these samples for research purposes during outbreaks of national concern is allowed under the terms of the 510/2016 Resolution of the National Ethical Committee for Research – Brazilian Ministry of Health (CONEP - Comissão Nacional de Ética em Pesquisa, Ministério da Saúde) that authorize, without the necessity of an informed consent, the use of clinical samples collected in the Brazilian Central Public Health Laboratories to accelerate knowledge building and contribute to surveillance and outbreak response. The samples processed in this study were obtained anonymously from material exceeding the routine diagnosis in Brazilian public health laboratories that belong to the public network within BrMoH.

### Epidemiological data

We analysed daily cases of SARS-CoV-2 in Brazil up to 2nd November 2021 from the COVIDA network available at [https://github.com/wcota/covid19br](https://github.com/wcota/covid19br). For convenience, the geographical locations were aggregated by Brazilian macro regions: North, Northeast, Southeast, South, and Midwest. The National COVIDA network releases daily updates on the number of confirmed new cases, deaths, and recoveries, with a breakdown by states and regions.

### Sample collection and molecular diagnostic assays

As part of the National Pandemic Alert Network, since early January 2021, our team has been monitoring the proportion of circulating variants in Brazil. For this purpose, convenience clinical samples, from public laboratories in Brazil were received and randomly selected for sequencing every week. Depending on the partner institution, library preparation and sequencing was done either on the Illumina and/or Oxford Nanopore Platform. Viral RNA was extracted from nasopharyngeal swabs using an automated protocol and tested for SARS-CoV-2 by multiplex real-time PCR assays: (i) the Allplex 2019-nCoV Assay (Seegene) targeting the envelope (E), the RNA dependent RNA polymerase (RdRp) and the nucleocapsid (N) genes; (ii) the Charité: SARS-CoV2 (E/RP) assay (Bio- Manguinhos/Fiocruz) targeting the E gene, and (iii) the GeneFinder COVID-19 Plus RealAmp Kit (Osang Healthcare, South Korea) supplied by the Brazilian Ministry of Health (BrMoH), Butantan Institute and the Pan-American Health Organization (OPAS).

### cDNA synthesis and whole genome sequencing

Samples were selected for sequencing based on the Ct value (≤30) and availability of epidemiological metadata, such as date of sample collection, sex, age and municipality of residence. The preparation of SARS-CoV-2 genomic libraries was performed using two different strategies: (i) the Illumina COVIDSeq test following the manufacturer’s instructions11; and (ii) the Oxford Nanopore sequencing using the ARTIC Network primal scheme24. The normalized libraries were loaded for the Illumina sequencing onto a 300-cycle MiSeq Reagent Kit v2 and run on the Illumina MiSeq instrument (Illumina, San Diego, CA, USA) and for the Nanopore strategy into a R9.4 flow cell (Oxford Nanopore Technologies) as previously described25. All experiments were performed in a biosafety level-2 cabinet. In each sequencing run, we used negative controls to prevent and check for possible contamination with less than 2% mean coverage.

### Genome assembly

Sequences generated on the Illumina and nanopore platforms were assembled using Genome Detective 1.132/326.

### Estimating relative transmission advantage

We analyzed 11,147 SARS-CoV-2 Brazilian sequences, from which n=6,626 were generated by our National Pandemic Alert Network, that have been uploaded to GISAID (Global Initiative On Sharing All Influenza Data) from 26 April 2021 to 23 October 2021. We used a multinomial logistic regression model to estimate the growth advantage of Delta compared to Gamma in Brazil15,27 We added splines to account for time-varying growth rates in the model fit and estimated the overall growth advantage without splines. We fitted the model using the *multinom* function of the *nnet* package (Venables & Ripley)28 in R.

### Phylogenetic analysis

We analyzed 11,147 Delta variants from Brazil, publicly available on GISAID21 as at the 4 November 2021, of which n=6,626 were generated by our National Pandemic Alert Network and used in this study. These were put in their phylogenetic context through comparison with a globally representative (n=13,261) set of other SARS-CoV-2 Delta variants from around the world, sampled from October 23th 2020 to October 23th 2021. The full set of sequences were aligned with NextAlign29 to obtain a good codon quality alignment against the Wuhan-Hu-1 universal reference sequence. The subsequent alignment was then used to infer a Maximum Likelihood tree topology in IQTREE230 (-m GTR, -b 100). Transfer bootstrap support for splits in the topology was inferred using Booster31. The resulting consensus ML tree topology was assessed for molecular clock signal in TempEst32. Potential outlier sequences or sequences lacking required metadata (e.g. date and location of sampling) were pruned off the topology with the ape package33 in R prior to dating. The branches in the ML-tree topology were then converted into units of calendar time in TreeTime34 using a constant rate of 0.0008 substitutions/site/year35 with a clock standard deviation of 0.0004 substitutions/site/year. Following the dating of the phylogeny we annotated the tips and internal nodes using the “mugration” package extension of TreeTime and then counted the state changes from one country to another and their inferred time points. This gave us the number and time of SARS-CoV-2 Delta isolates entering Brazil over time. To obtain a measure of confidence in the time and source of viral transitions, we performed the discrete ancestral state reconstruction on 10 bootstrap replicate trees.

## Additional information

**Supplementary Figure 1**. Dynamics of the SARS-CoV-2 epidemic in each Brazilian showing the number of daily COVID-19 deaths and the progression in the proportion of circulating variants in the country over time, showing the rapid replacement of Gamma by Delta.

**Supplementary Figure 2**. Progressive distribution of the top 10 Delta PANGO sub-lineages in Brazil over time.

**Supplementary Table 1**. GISAID accession numbers of sequences used in the phylogenetic analysis.

## Acknowledgments

The authors acknowledge the National Network for Pandemic Alert of SARS-CoV-2 and the contribution of all employees of General Coordination of Public Health Laboratories and professionals of Public Health Laboratories of Brazil for their contribution towards the sequencing effort and for their commitment and work during the fight of the COVID-19 pandemic. We also thank all the authors who have kindly deposited and shared genome data on GISAID.

*   Received December 27, 2021.
*   Revision received December 27, 2021.
*   Accepted December 29, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Safari, I. & Elahi, E. Evolution of the SARS-CoV-2 genome and emergence of variants of concern. Arch Virol (2021) doi:10.1007/s00705-021-05295-5.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00705-021-05295-5&link_type=DOI) 

2.  2.The biological and clinical significance of emerging SARS-CoV-2 variants | Nature Reviews Genetics. [https://www.nature.com/articles/s41576-021-00408-x](https://www.nature.com/articles/s41576-021-00408-x).
    
    

3.  3.Nonaka, C. K. V. et al. Genomic Evidence of SARS-CoV-2 Reinfection Involving E484K Spike Mutation, Brazil - Volume 27, Number 5—May 2021 - Emerging Infectious Diseases journal - CDC. doi:10.3201/eid2705.210191.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2705.210191&link_type=DOI) 

4.  4.Greaney, A. J. et al. Comprehensive mapping of mutations in the SARS-CoV-2 receptor-binding domain that affect recognition by polyclonal human plasma antibodies. Cell Host & Microbe 29, 463-476.e6 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.CHOM.2021.02.003&link_type=DOI) 

5.  5.Meng, B. et al. Recurrent emergence of SARS-CoV-2 spike deletion H69/V70 and its role in the Alpha variant B.1.1.7. Cell Reports 35, 109292 (2021).
    
    

6.  6.Tegally, H. et al. Detection of a SARS-CoV-2 variant of concern in South Africa. Nature 592, 438–443 (2021).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F29%2F2021.12.27.21268309.atom) 

7.  7.Faria, N. R. et al. Genomics and epidemiology of the P.1 SARS-CoV-2 lineage in Manaus, Brazil. Science 372, 815–821 (2021).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNzIvNjU0NC84MTUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMi8yOS8yMDIxLjEyLjI3LjIxMjY4MzA5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

8.  8.Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa | medRxiv. [https://www.medrxiv.org/content/10.1101/2021.12.19.21268028v1](https://www.medrxiv.org/content/10.1101/2021.12.19.21268028v1).
    
    

9.  9.Planas, D. et al. Reduced sensitivity of SARS-CoV-2 variant Delta to antibody neutralization. Nature 596, 276–280 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/2021.05.26.445838&link_type=DOI) 

10. 10.WHO | World Health Organization. [https://www.who.int/](https://www.who.int/).
    
    

11. 11.Genomic epidemiology reveals how restriction measures shaped the SARS-CoV-2 epidemic in Brazil | medRxiv. [https://www.medrxiv.org/content/10.1101/2021.10.07.21264644v1](https://www.medrxiv.org/content/10.1101/2021.10.07.21264644v1).
    
    

12. 12.COVID-19 Situation Reports - PAHO/WHO | Pan American Health Organization. [https://www.paho.org/en/covid-19-situation-reports](https://www.paho.org/en/covid-19-situation-reports).
    
    

13. 13.Rodriguez-Morales, A. J. et al. COVID-19 in Latin America: The implications of the first confirmed case in Brazil. Travel Medicine and Infectious Disease 35, 101613 (2020).
    
    

14. 14.Ganyani, T. et al. Estimating the generation interval for coronavirus disease (COVID-19) based on symptom onset data, March 2020. Eurosurveillance 25, (2020).
    
    

15. 15.Eurosurveillance | Increased transmissibility and global spread of SARS-CoV-2 variants of concern as at June 2021. [https://www.eurosurveillance.org/content/10.2807/1560-7917.ES.2021.26.24.2100509](https://www.eurosurveillance.org/content/10.2807/1560-7917.ES.2021.26.24.2100509).
    
    

16. 16.Rapid replacement of the Beta variant by the Delta variant in South Africa | medRxiv. [https://www.medrxiv.org/content/10.1101/2021.09.23.21264018v1](https://www.medrxiv.org/content/10.1101/2021.09.23.21264018v1).
    
    

17. 17.Replacement of the Alpha variant of SARS-CoV-2 by the Delta variant in Lebanon between April and June 2021 | medRxiv. [https://www.medrxiv.org/content/10.1101/2021.08.10.21261847v1](https://www.medrxiv.org/content/10.1101/2021.08.10.21261847v1).
    
    

18. 18.The SARS-CoV-2 variant Delta displaced the variants Gamma and Gamma plus in Amazonas, Brazil - SARS-CoV-2 coronavirus / nCoV-2019 Genomic Epidemiology - Virological. [https://virological.org/t/the-sars-cov-2-variant-delta-displaced-the-variants-gamma-and-gamma-plus-in-amazonas-brazil/765](https://virological.org/t/the-sars-cov-2-variant-delta-displaced-the-variants-gamma-and-gamma-plus-in-amazonas-brazil/765).
    
    

19. 19.Marques, A. D. et al. SARS-CoV-2 variants associated with vaccine breakthrough in the Delaware Valley through summer 2021. [http://medrxiv.org/lookup/doi/10.1101/2021.10.18.21264623](http://medrxiv.org/lookup/doi/10.1101/2021.10.18.21264623) (2021) doi:10.1101/2021.10.18.21264623.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4xMC4xOC4yMTI2NDYyM3YyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMjkvMjAyMS4xMi4yNy4yMTI2ODMwOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

20. 20.Olsen, R. J. et al. Trajectory of Growth of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Variants in Houston, Texas, January through May 2021, Based on 12,476 Genome Sequences. The American Journal of Pathology 191, 1754–1773 (2021).
    
    

21. 21.Emergence and Spread of the SARS-CoV-2 Variant of Concern Delta Across Different Brazilian Regions | medRxiv. [https://www.medrxiv.org/content/10.1101/2021.11.25.21266251v1](https://www.medrxiv.org/content/10.1101/2021.11.25.21266251v1).
    
    

22. 22.Candido, D. D. S. et al. Routes for COVID-19 importation in Brazil. Journal of Travel Medicine 27, taaa042 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jtm/taaa042&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32211799&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F29%2F2021.12.27.21268309.atom) 

23. 23.Giovanetti, M. et al. Genomic epidemiology reveals how restriction measures shaped the SARS-CoV-2 epidemic in Brazil. 2021.10.07.21264644 [https://www.medrxiv.org/content/10.1101/2021.10.07.21264644v1](https://www.medrxiv.org/content/10.1101/2021.10.07.21264644v1) (2021) doi:10.1101/2021.10.07.21264644.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4xMC4wNy4yMTI2NDY0NHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMjkvMjAyMS4xMi4yNy4yMTI2ODMwOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

24. 24.Artic Network. [https://artic.network/ncov-2019](https://artic.network/ncov-2019).
    
    

25. 25.Xavier, J. et al. The ongoing COVID-19 epidemic in Minas Gerais, Brazil: insights from epidemiological data and SARS-CoV-2 whole genome sequencing. Emerging Microbes & Infections 9, 1824–1834 (2020).
    
    

26. 26.Vilsker, M. et al. Genome Detective: an automated system for virus identification from high-throughput sequencing data. Bioinformatics 35, 871–873 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bty695&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30124794&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F29%2F2021.12.27.21268309.atom) 

27. 27.Davies, N. G. et al. Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England. Science 372, eabg3055 (2021).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiIzNzIvNjUzOC9lYWJnMzA1NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzI5LzIwMjEuMTIuMjcuMjEyNjgzMDkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

28. 28.Modern Applied Statistics with S, 4th ed. [https://www.stats.ox.ac.uk/pub/MASS4/](https://www.stats.ox.ac.uk/pub/MASS4/).
    
    

29. 29.GitHub - nextstrain/nextclade: Viral genome alignment, mutation calling, clade assignment, quality checks and phylogenetic placement. [https://github.com/nextstrain/nextclade](https://github.com/nextstrain/nextclade).
    
    

30. 30.Nguyen, L.-T., Schmidt, H. A., von Haeseler, A. & Minh, B. Q. IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Molecular Biology and Evolution 32, 268–274 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msu300&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25371430&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F29%2F2021.12.27.21268309.atom) 

31. 31.Lemoine, F. et al. Renewing Felsenstein’s phylogenetic bootstrap in the era of big data. Nature 556, 452–456 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-018-0043-0&link_type=DOI) 

32. 32.Rambaut, A., Lam, T. T., Max Carvalho, L. & Pybus, O. G. Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen). Virus Evol 2, vew007 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ve/vew007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27774300&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F29%2F2021.12.27.21268309.atom) 

33. 33.ape 3.0: New tools for distance-based phylogenetics and evolutionary analysis in R | Bioinformatics | Oxford Academic. [https://academic.oup.com/bioinformatics/article/28/11/1536/267314](https://academic.oup.com/bioinformatics/article/28/11/1536/267314).
    
    

34. 34.Sagulenko, P., Puller, V. & Neher, R. A. TreeTime: Maximum-likelihood phylodynamic analysis. Virus Evol 4, vex042 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ve/vex042&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29340210&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F29%2F2021.12.27.21268309.atom) 

35. 35.Wilkinson, E. et al. A year of genomic surveillance reveals how the SARS-CoV-2 pandemic unfolded in Africa. Science , eabj4336.