Spatiotemporal dynamics and epidemiological impact of SARS-CoV-2 XBB lineages dissemination in Brazil in 2023
=============================================================================================================

* Ighor Arantes
* Marcelo Gomes
* Kimihito Ito
* Sharbilla Sarafim
* Fábio Miyajima
* Ricardo Khouri
* Tiago Gräf
* Marilda Mendonça Siqueira
* Paola Cristina Resende
* Felipe Gomes Naveca
* Gonzalo Bello
* COVID-19 Fiocruz Genomic Surveillance Network

## ABSTRACT

The SARS-CoV-2 XBB is a group of highly immune-evasive lineages of the Omicron VOC that emerged by recombining BA.2-descendent lineages and spread worldwide during 2023. In this study, we combine SARS-CoV-2 genomic data (*n =* 11,065 sequences) with epidemiological data of Severe Acute Respiratory Infection (SARI) cases collected in Brazil between October 2022 and July 2023 to reconstruct the space-time dynamics and epidemiologic impact of XBB dissemination in the country. Our analyses revealed that the introduction and local emergence of lineages carrying convergent mutations within the Spike protein, especially F486P, F456L, and L455F, propelled the spread of XBB* lineages in Brazil. The average relative instantaneous reproduction numbers of XBB*+F486P, XBB*+F486P+F456L, and XBB*+F486P+ F456L+L455F lineages in Brazil were estimated to be 1.24, 1.33, and 1.48 higher than that of other co-circulating lineages (mainly BQ.1\*/BE\*), respectively. Despite such a growth advantage, the dissemination of these XBB* lineages had a reduced impact on Brazil’s epidemiological scenario concerning previous Omicron subvariants. The peak number of SARI cases from SARS-CoV-2 during the XBB wave was approximately 90%, 80%, and 70% lower than that observed during the previous BA.1*, BA.5*, and BQ.1* waves, respectively. These findings revealed the emergence of multiple XBB lineages with progressively increasing growth advantage, yet with relatively limited epidemiological impact in Brazil throughout 2023. The XBB*+F486P+F456L+L455F lineages stand out for their heightened transmissibility, warranting close monitoring in the months ahead.

## INTRODUCTION

Since its emergence in late 2021, the SARS-CoV-2 variant of concern (VOC) Omicron has not only replaced previous VOCs but has also undergone rapid evolutionary changes, giving rise to a multitude of sub-lineages that spread around the world (1). The XBB lineage is one of these Omicron sub-lineages that emerged around July 2022 by the recombination of BA.2-descendent lineages BJ.1 and BM.1.1.1. Its initial dissemination primarily occurred in Asia, with Singapore witnessing a sharp upsurge in cases attributed to this strain (2–5). Although the XBB lineage demonstrated remarkable immune evasiveness and an effective reproduction number (Re) ∼1.2 times higher than those of parental variants (5–20), it displayed a weaker ACE2 binding affinity than other immune evasive Omicron subvariants that circulated during the latter half of 2022, such as the BQ.1* lineages prevalent in Western countries (5, 21, 22).

In late 2022, a new XBB lineage designated XBB.1.5 arose and rapidly spread worldwide. The XBB.1.5 lineage harbors one additional mutation (S:F486P) in the receptor binding domain (RBD) of the Spike that significantly improves the XBB’s ability to attach to the human ACE2 receptor without eroding the variant’s ability to evade humoral immunity (21, 23–27). The Re of XBB.1.5 was estimated to be ∼1.2 times greater than that of the parental lineage XBB.1 and outcompeted BQ.1* lineages (24). Throughout 2023, multiple XBB lineages carrying mutation S:F486P evolved independently in different countries (1). This convergent evolution is the case of the lineage XBB.1.16 that initially spread in India and displayed similar characteristics to lineage XBB.1.5 regarding host cell entry efficiency, neutralization evasion, and growth advantage (17, 28, 29). Thus, increased immune evasion and high receptor-binding affinity resulted in a substantial growth advantage of XBB*+S:F486P relative to previously circulating Omicron lineages.

In 2023, multiple XBB*+F486P subvariants acquired an additional convergent mutation S:F456L in the RBD of the Spike and displayed a growth advantage relative to the XBB variant precursors (30–32). Mutation S:F456L reduces the ACE2 binding affinity but increases the variant’s immune evasion capabilities, even after previous XBB infection (17, 31, 33–36). The most prevalent lineage with this mutation is EG.5, which evolved from the lineage XBB.1.9.2 (S:F486P), showing one of the fastest growth rates among XBB lineages circulating in North America and Europe in early 2023 (30, 37–39). By late March 2023, another group of XBB*+F486P lineages harboring mutations S:L455F and S:F456L, such as lineages XBB.1.5.70, GK*, JD.1* and HK.3, were detected for the first time and are rapidly spreading (40). These lineages exhibit higher immune evasion than XBB*+F486P lineages and heightened ACE2 binding affinity relative to XBB*+F486P+F456L lineages (41, 26), which may thus lead to further growth advantage.

Brazil was one of the most heavily affected countries in the world by the SARS-CoV-2 pandemic). However, little is known about the molecular evolution trajectory of XBB lineages circulating in the country. In this study, we analyzed the temporal transition and transmission advantage of XBB lineages that spread in Brazil during the first half of 2023. We reconstruct the spatiotemporal spread dynamics of major XBB*+F486P, XBB*+F486P+F456L, and XBB*+F486P+F456L+L455F lineages circulating in Brazil to assess the epidemiological impact of their dissemination by analyzing the hospitalizations related to severe acute respiratory illness (SARI) due COVID-19, over the Omicron waves.

### METHODOLOGY

#### SARS-CoV-2 Brazilian genome sequences

A total of 11,065 SARS-CoV-2 complete genome sequences recovered across all 26 Brazilian states plus the country’s Federal District between October 1st, 2022, and July 31st, 2023, were newly generated by the COVID-19 Fiocruz Genomic Surveillance Network (42). All samples had real-time RT-PCR cycling threshold (Ct) below 30, indicating elevated viral load. There was no intended bias considering dense sampling of specific outbreaks, and we do not consider the health status (severe, mild, or asymptomatic) or epidemiological data (gender and age) of individuals as a criterion for sampling. SARS-CoV-2 genome sequences were generated using the Illumina COVIDSeq Test kit as previously described (43). Raw data were converted to FASTQ files at Illumina BaseSpace cloud, and consensus sequences were produced with the most up-to-date version of DRAGEN COVID LINEAGE or ViralFlow 1.0 (44). All genomes were evaluated for mutation calling and quality with the Nextclade 2.14.0 algorithm (45) and were uploaded to the EpiCoV database of GISAID ([https://gisaid.org/](https://gisaid.org/)) (**Supplementary Material**). Additionally, we downloaded all publicly available Brazilian sequences (*n =* 14,822) collected since October 1st, 2022, that were submitted to the EpiCov database until July 31st, 2023 with complete collection date and lineage assignment. Whole-genome consensus sequences were classified using the ‘Phylogenetic Assignment of Named Global Outbreak Lineages’ (PANGOLIN) software v4.3 (pangolin-data version 1.21) (46).

#### Non-Brazilian XBB genome sequences and dataset composition

To elucidate the spatiotemporal dynamics of prominent Brazilian XBB lineages that exhibited the S:F486P, S:F456L, and S:L455F mutations, we retrieved all non-Brazilian sequences pertaining to lineages XBB.1.4.2 (*n =* 16), XBB.1.5.86 + HA.1 (*n =* 478), XBB.1.18* + FE.1* (*n =* 1,772), XBB.1.5.70 + GK* (*n =* 350), and XBB.1.5.102 + JD.1* (*n =* 174) that were submitted to the EpiCov database until July 31st 2023 (**Supplementary Material**). To obtain a final dataset for each studied lineage, all genomes, whether Brazilian or foreign, were filtered based on (i) the presence of all described synapomorphic mutations of their respective lineages, (ii) a complete collection date, and (iii) a small fraction of unidentified positions (N < 5%). To enhance the accuracy of the root position in our phylogenetic reconstructions, we obtained all high coverage (N < 1%) genomes with complete collection dates from ancestral lineages for potential Brazilian ones, namely XBB.1 (*n =* 811), XBB.1.5 (*n =* 3,886), and XBB.1.4 (*n =* 12). To streamline our computational analyses without compromising lineage diversity, we clustered them, when needed, with CD-HIT (47), leading to datasets of much-reduced dimensions (XBB.1, *n =* 3 and XBB.1.5, *n =* 118).

#### Maximum-likelihood (ML) phylogenetic analysis

Seven XBB datasets of 29,421 nts. (from position 1 of ORF1ab to position 117 of ORF10) were produced by a combination of Brazilian and non-Brazilian sequences belonging to the following lineages: (i) XBB.1.18* + FE.1 (*n =* 1,843), (ii) FE.1.1 (*n =* 1,251), (iii) FE.1.2 (*n =* 1,051), (iv) XBB.1.5.102 (*n =* 1,522), (v) XBB.1.4.2 (*n =* 132), (vi) XBB.1.5.86 + HA.1 (*n =* 1,010) (vii) XBB.1.5.70 + GK* + JD.1 (*n =* 716). These were aligned using MAFFT v7.467 (48) and manually curated when needed. The aligned datasets were subject to ML phylogenetic analysis using IQ-TREE v2.2.2.7 (49) under the general time-reversible (GTR) model of nucleotide substitution with a gamma-distributed rate variation among sites, four rate categories (G4), a proportion of invariable sites (I) and empirical base frequencies (F) nucleotide substitution model, as selected by the ModelFinder application embedded in IQ-TREE (50). The approximate likelihood-ratio test assessed the branch support based on the Shimodaira–Hasegawa-like procedure (SH-aLRT) with 1,000 replicates.

#### Bayesian time-scaled phylogenetic and phylogeographic analyses

The temporal signal of the XBB datasets was assessed from the ML trees by performing a regression analysis of the root-to-tip divergence against sampling time using TempEst v1.5.36 (51). Sequences that diverged more than 1.5 interquartile ranges from the root-to-tip regression were considered outliers and removed from the analysis. The significance of the association between the two variables was assessed with a Spearman correlation test implemented in the R programming language v.4.1.2. Time-scaled Bayesian phylogenetic trees were estimated in BEAST v1.10.4 (52) under a strict molecular clock model with a uniform distribution between 5.0E-4 and 1.0E-3, and a non-parametric Bayesian skyline (BSKL) model as the coalescent tree prior (53). Due to the large size of some of our datasets (*n* > 500 sequences), time-scaled trees were reconstructed using a modified version of BEAST ([https://beast.community/thorney_beast](https://beast.community/thorney_beast)) that incorporates the possibility of a fixed ML topology. This method alleviates the computational burden by rescaling the input tree branch lengths based on (i) the clock and (ii) coalescent models. The ML phylogenetic trees were inputted in the BEAST XML file as a starting point, and data trees and analyses were performed as specified above. For phylogeographic analyses, a set of 1,000 trees was randomly selected from the posterior distribution of trees resulting from the temporal BEAST analysis. Sampling locations were used as traits in the phylogeographic model, and the ancestral states were reconstructed using a discrete symmetric model (54). Samples were grouped by continental regions, except Brazilian sequences that were grouped in a separate location. Major Brazilian transmission clusters were defined as highly supported (approximate likelihood-ratio test, a*LRT* > 0.75) monophyletic clades that descend from an MRCA node probably located (posterior state probability, *PSP* > 0.75) in Brazil. MCMC simulations were run sufficiently long to ensure convergence (effective sample size, ESS > 200) in all parameter estimates as assessed in TRACER v1.770 (55). The maximum clade credibility (MCC) trees were summarized with TreeAnnotator v1.10 and visualized using treeio v3.1.7 (56) and ggtree v3.2.1 (57).

#### Relative instantaneous reproduction number (RRI) estimations

We measured the transmission advantage of XBB lineages carrying the S:F486P, S:F486P+F456L, and S:F486P+F456L+L455F mutations compared to other lineages using their *RRI*. The details of the estimation method have been described elsewhere (58). Briefly, the *RRI* is an effective reproduction number representing the average number of people an infected individual at time *t* could be expected to infect, given that conditions remain unchanged (59). Suppose that lineage *A*1, …, *A*n, and *a* are circulating in the population. The lineage *a* is used as the baseline of the comparison. Let *R**A*1 (*t*), …, *R**A*n(*t*), and *R**I*(*t*) be the *RRI* of lineages *A*1, …, *A*n, and *a* at time *t*, respectively. Assuming the rate of *R**Ai*(*t*) to *R**I*(*t*) is constant over time for each lineage*A**i*(1 ≤ *i* ≤ n), the effective reproduction number of *A**i* can be written as: ![Formula][1]</img>  where *ki* is a constant. This constant *ki* is called the *RRI* of lineage *Ai* with respect to (w.r.t.) lineage *a*. It was shown that the time evolution of (relative) frequencies of lineages *A*1, …, *A*n, and *a* among a viral population is well-modeled without the number of infections (60, 61). Let *q**Ai*(*t*) and *q**I*(*t*) be the frequencies of lineages *Ai* and *a* respectively. Then *q**Ai*(*t*)and *q**I*(*t*) are represented by the following formula: ![Formula][2]</img>  where *g*(*j*) is the probability mass function of generation time distribution and *l* is the day when the probability of generating new infections becomes negligibly small. Park *et al.* estimated the Omicron strains’ mean generation time to be 3.0 using within-household transmission data (62). The generation time distribution they estimated was a lognormal distribution with a log mean of 0.98 and a log standard deviation of 0.47, of which the mean is 2.97, and the variance is 2.19. Using these values, the probability mass function of the generation time *g*(*j*) was modeled by discretizing the gamma distribution having the same mean and variance as follows: ![Formula][3]</img>  where α and θ are the shape and scale parameters of gamma distribution (Gamma) and the CDF is the cumulative distribution function up to the given point. Parameters of the gamma distribution were set to be α = 4.03 and θ = 0.737, so that the gamma distribution has the same mean and variance to the lognormal distribution.

We consider XBB* lineages carrying S:F486P, S:F486P+F456L, and S:F486P+F456L+L455F mutations as *A*1, *A*2, and *A*3, respectively, and other lineages (mainly BQ.1\*/BE\*) were used as the baseline lineage *a*. We assume that viruses of these lineage groups in a region were sampled according to a Dirichlet multinomial distribution with the shape parameters *q**A*1 (*t*)*M*, *q**A*2 (*t*)*M*, *q**A*3 (*t*)*M*, and *q**I*(*t*)*M* at time *t* where *M* represent the sum of the Dirichlet distribution parameters. We also assume that the relative reproduction numbers of XBB*+S:F486P, XBB*+S:F486P+F456L, and XBB*+S:F486P+F456L+ L455F w.r.t Others are the same in all regions. Lineage groups observed in each region were counted using weekly bins starting from Monday. By maximizing the likelihood function of the multinomial distribution in all regions, we estimated the relative reproduction numbers *k*1, *k*2, *k*3, the sum of parameters Dirichlet distribution parameters *M*, and hypothetical initial frequencies *q**A*1 (*t*), *q**A*2 (*t*), *q**A*3 (*t*) in each region on October 3, 2022. The estimation of the date when each lineage group was introduced in each region is not the purpose of this specific analysis. Instead, we estimated hypothetical initial frequencies of each one, *q**Ai*(*t*), if it existed in each region on October 3, 2022. Thus, a late introduction of a lineage *A**i*to a region results in a small estimate of *q**Ai*(*t*), in the region. The 95% confidence intervals (CI) of estimates were estimated by profile likelihood. The trajectories of variant frequencies were calculated using the maximum likelihood estimates of the parameters in the model. The 95% CI of the trajectories was calculated using combinations of parameters within the 95% confidence region (63).

#### Data on hospitalizations for severe acute respiratory illness (SARI)

We extracted data about hospitalizations resulting from severe acute respiratory illness (SARI), attributed explicitly to SARS-CoV-2 (SARI-COVID), in Brazil during the period spanning from January 2022 to July 2023. This information was sourced from the Influenza Surveillance Information System (SIVEP-Gripe) database ([https://sivepgripe.saude.gov.br/sivepgripe/login.html?0](https://sivepgripe.saude.gov.br/sivepgripe/login.html?0)). To identify SARI-COVID cases, we utilized a set of four criteria, which required individuals to exhibit: (i) fever, including self-reported cases, (ii) cough or sore throat, (iii) dyspnea or oxygen saturation levels below 95% or experiencing respiratory discomfort and (iv) hospitalization. Once an individual meets these criteria and is admitted to a hospital for SARI, their case must be reported and recorded as a distinct entry in the SIVEP-Gripe database. Confirmation of hospitalization for SARI-COVID was contingent on a positive result from the RT-PCR molecular test for SARS-CoV-2. For our analysis, we considered all the records within the SIVEP-Gripe database that adhered to the criteria for defining a hospitalized SARI case. We excluded records related to non-hospitalized deaths from our examination.

## RESULTS

### Circulation of SARS-CoV-2 Omicron subvariants in Brazil in 2023

We analyzed 25,887 Brazilian SARS-CoV-2 genomes collected across all states between October 1st, 2022, and July 31st, 2023 (**Fig. 1a**). Sequences were mostly classified as XBB* (28%) or BA.5* (72%) lineages. While most BA.5* sequences were from the Southeast region (**Fig. 1b**), most XBB* sequences were collected in the North region (**Fig. 1c**). The XBB* sequences were first detected in Brazil in October 2022 and circulated at a low prevalence (<10%) during that year. However, their proportion experienced a rapid surge in 2023 going from ∼15% to ∼90% between January and February (**Fig. 1d**). The XBB* sequences were categorized in four groups based on mutations at the RBD region of the Spike protein: XBB*, XBB*+F486P, XBB*+F486P+F456L, and XBB*+F486P+F456L+L455F. The temporal analysis revealed a clear pattern of substitution of XBB groups over time (**Fig. 1e**). The XBB* group was first detected in October 2022 and displayed a very low prevalence (∼5%) over the studied period. The dominance of XBB lineages in the Brazilian epidemic in 2023 coincides with the introduction and dissemination of the XBB*+F486P group that was first detected by November 2022 and reached its peak in relative prevalence by February 2023. The XBB*+F486P+F456L and XBB*+F486P+F456L+L455F groups emerged later but progressively assumed dominance in the country. The XBB*+F486P+F456L group was first detected by January 2023 and reached its zenith in relative prevalence by May 2023 whereas the XBB*+F486P+F456L+L455F group was first detected by March 2023 and surpassed the XBB*+F486P and XBB*+F486P+F456L groups three and four months later, respectively. The relative instantaneous reproduction number (*RRI*) of XBB* lineages carrying additional Spike mutations *w.r.t* “Other” lineages circulating in Brazil in 2023 was calculated using the renewal-equation-based model developed by Ito *et al* (60) (**Fig. 1f**). The average *RRI* of XBB*+F486P, XBB*+F486P+F456L and XBB*+F486P+F456L+L455F lineages in Brazil were estimated to be 1.24, 1.33 and 1.48 higher than that of other co-circulating lineages (mainly BQ.1\*/BE\*), respectively. Thus, the average *RRI* of XBB*+F486P+F456L lineages was 1.07 times higher than that of XBB*+F486P lineages and the average *RRI* of XBB*+F486P+ F456L+L455F lineages was 1.19 and 1.11 times higher than that of XBB*+F486P and XBB*+F486P+F456L lineages, respectively. The same pattern of XBB subvariants substitutions over time was consistently observed in the whole country, although the current prevalence of the different XBB groups may vary across Brazilian regions (**Fig. 1g**), and across all age groups (**Fig. 1h**).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/30/2023.10.30.23297466/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2023/10/30/2023.10.30.23297466/F1)

Figure 1. Space-time evolution of XBB subvariants in Brazil between October 2022 and July 2023.
**a.** Monthly) distribution of the 25,887 Brazilian SARS-CoV-2 genomes that constituted the dataset for this study. **b-c**. Spatial distribution of SARS-CoV-2 genomes across the five Brazilian regions of the 18,518 non-XBB (**b**) and the 7,369 XBB genomes (**c**) **d**. Temporal fluctuations in the prevalence of both XBB and non-XBB variants in Brazil. **e**. Evolution of distinct XBB subvariants groups relative prevalence in Brazil. The groups are defined by the mutations S:F486P, S:F456L, and S:L455F, all located in the RBD region of the Spike protein. Genomes of XBB subvariants in which all three mutations were absent and those of non-XBB lineages were grouped in the “Others” category. **f.** Relative Instantaneous Reproductive Number (RRI) of XBB* groups w.r.t the “Others” category. **g-h.** Evolution of distinct XBB subvariants groups relative prevalence across the five Brazilian regions (**g**), and distinct age groups (**h**). t. Maps were generated with QGIS v.3.10.2 software ([http://qgis.org](http://qgis.org)) using public access data downloaded from the GADM v.3.6 database ([https://gadm.org](https://gadm.org)) and shapefiles obtained from the Brazilian Institute of Geography and Statistics ([https://portaldemapas.ibge.gov.br/portal.php#homepage](https://portaldemapas.ibge.gov.br/portal.php#homepage)).

### Emergence and dissemination of major XXB lineages in Brazil

The assignment of dynamic Pango-lineages to the Brazilian XBB genomes revealed a heterogeneous composition within each of the groups and further revealed that mutations S:F486P/F456L/L455F arose independently in multiple lineages (**Fig. 2a**). The XBB*+F486P group was mainly composed of lineages XBB.1.5.102 (34%), XBB.1.18.1 (19%), XBB.1.5.86 (12%), XBB.2.3* (3%), XBB.1.4.2 (2%), and others XBB.1.5-descendant lineages (∼25%) (**Fig. 2b**). The XBB*+F486P+F456L group was mostly composed of XBB.1.18.1-descendant lineages such as FE.1.2 (∼50%), FE.1.1 (∼20%), and FE.1 (∼10%) (**Fig. 2b**). The XBB*+F486P+F456L+L455F group was mostly composed of lineage XBB.1.5.70 (∼50%) and their descendent lineages GK.1 (∼30%) and GK.1.1 (∼5%) and by the XBB.1.5.102-descendent lineage JD.1 (∼5%) and JD.1.2 (∼10%) (**Fig. 2b**). The relative frequency of the different XBB lineages within each group varied across Brazilian regions, particularly in the Northern region (**Fig. 2c**). There, the most prevalent lineage in the first group was the XBB.1.5.102 (∼80%), in the second group was an XBB.1.18.1+F456L lineage distinct of FE.1* (∼40%), while in group three prevails the JD.1* (∼40%). All these three variants were scarcely detected in other regions of the country.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/30/2023.10.30.23297466/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2023/10/30/2023.10.30.23297466/F2)

Figure 2. Major XBB* lineages that arose and spread in Brazil.
**a**. Fluxogram representing the major XBB* Pango lineages that emerged and spread in Brazil, focused on the sequential accumulation of key Spike protein mutations (S:F486P, S:F456L, and S:L455F). Lineages that evolved in Brazil are in dark gray, while foreign lineages or those not included in out phylogeographic analysis are in white with a gray halo. **b-c.** Relative frequency of major Pango XBB* lineages bearing the S:F486P, S:F486P+F456L, and S:F486P+F456L+L455F mutations in Brazil (b) and in different country regions (c). **d.** Temporal dynamics of major Pango XBB* lineages that emerged in Brazil. For each XBB* lineage, we represent the TMRCA (circle) and its 95% HPD interval (transparent polygon), the period of cryptic circulation (thinner line), and the sampling range (thicker line). All represented Pango XBB* lineages were inferred as having their MRCA most probably located in Brazil (PSP > 0.75).

Until the end of August 2023, among SARS-CoV-2 genomes submitted to GISAID ([https://gisaid.org/](https://gisaid.org/)), some of the major XBB* lineages detected in Brazil remained mostly restricted to its borders while others were frequently detected outside of the country. Within the XBB*+F486P group, most genomes (55-85%) of lineages XBB.1.18.1, XBB.1.5.86, XBB.1.5.102, and XBB.1.4.2 were sampled in Brazil, while most genomes (>95%) of lineages XBB.1.5 and XBB.2.3* were sampled outside the country. Within the XBB*+F486P+F456L group, lineages FE.1 and FE.1.2 were equivalently sampled in Brazil and abroad (∼50%), while lineages FE.1.1 and HA.1 were mainly detected outside Brazil (75-80%). Within the XBB*+F486P+ F456L+L455F group, JD.1 was mostly sampled in Brazil (∼85%), XBB.1.5.70, and JD.1.2 were equivalently sampled in Brazil and abroad (∼50%), and lineages GK.1, GK.1.1 and were mostly sampled outside the country (60-90%).

In order to ascertain the most probable geographic origin of major XBB* lineages detected in Brazil, we conducted Bayesian phylogeographic inferences by compiling all SARS-CoV-2 complete genomes pertaining to each aforementioned XBB lineage that were available in EpiFlu GISAID by the end of July 2023. Our phylogeographic analyses pinpointed Brazil as the likely origin (*PSP* > 75%) of lineages XBB.1.4.2, XBB.1.5.70, XBB.1.5.86, XBB.1.5.102, XBB.1.18, XBB.1.18.1, FE.1, FE.1.1, FE.1.2, GK.1, GK.1.1, JD.1 and JD.1.2 **(Fig. S1-4 and Table S1)**.

Brazilian autochthonous lineages carrying the S:F486P mutation arose between mid-October (XBB.1.18.1) and late December (XBB.1.4.2) of 2022, followed by Brazilian lineages harboring S:F486P+F456L mutations from mid-December 2022 (FE.1) to early March 2023 (HA.1), and Brazilian lineages harboring S:F486P+F456L+L455F mutations from mid-February (XBB.1.5.70) to early April (GK.1.1) of 2023 (**Fig. 2d and Table S1**). Some of these XBB lineages spread beyond Brazil’s borders and gave rise to new sub-lineages in other locations, including multiple ones in North America, such as GK.1.2 (*PSP* = 0.90), GK.1.3 (*PSP* = 1.0), GK.2 (*PSP* = 0.70), GK.3 (*PSP* = 1.0), and JD.1.1 (*PSP* = 1.0) (**Fig. S1-S4**).

### Epidemiological impact of dissemination of XXB in Brazil

To assess the epidemiological consequences of the proliferation of XBB lineages within the context of Omicron waves in Brazil, we analyzed the temporal trends of hospitalizations attributed to SARI-COVID from January 2022 to July 2023. The available data reveals a recurring pattern throughout 2022 characterized by three discernible waves occurring approximately six months apart, reaching their peaks in January, July, and December, under the dominance of BA.1*, BA.5*, and BQ* subvariants, respectively (**Fig. 3a**). These waves displayed a trend of diminishing magnitude over time, experiencing reductions of approximately 75% and 35%, respectively, compared to their preceding counterparts. As we enter the phase where XBB constitutes the most prevalent subvariant in the country in early 2023, there was a discrete increase in SARI-COVID cases around early March coinciding with the dissemination of the XBB*+F486P group, followed by a period of decline that persists throughout the dissemination of the XBB*+F486P+F456L and XBB*+F486P+F456L+L455L lineages (**Fig. 3b**). The number of SARI-COVID cases at the peak of the XBB* wave was approximately 90%, 80% and 70% lower than that observed during the peak of BA.1*, BA.5* and BQ.1* waves, respectively. Moreover, the number of SARI-COVID notifications was reduced by 60% from January-June to July-December 2022 and by another 60% from July-December 2022 to January-July 2023.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/30/2023.10.30.23297466/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2023/10/30/2023.10.30.23297466/F3)

Figure 3. Severe Acute Respiratory Infection (SARI) detection in Brazil.
**a.** Evolution of SARI cases (yellow line) in Brazil between January 2022 and July 2023, plotted under the backdrop of SARS-CoV-2 Omicron lineages dissemination in the country, according to the color code at the panel’s bottom. **b-c.** Evolution of SARI cases in Brazil between October 2022 and July 2023, plotted under the backdrop of the distinct XBB* groups emergence and dissemination in Brazil (**b**), its five composing regions (**c**), and distinct age groups (**d**). The non-XBB and XBB* lineages lacking the S:F486P, S:F456L, and S:L455F mutations are aggregated and classified under the category “Others”. In all graphs, genomic and epidemiological data are plotted by epidemiological week.

When examining the incidence of SARI-COVID hospitalizations in the context of the spread of XBB* lineages across different Brazilian regions and age groups, a generally analogous trend was discerned, albeit with some variations (**Fig. 3c and 3d**). In the Northeastern, Central-Western, and Southern regions, the surge of SARI-COVID cases was temporally shifted to May. Moreover, the peak of the XBB* wave in the Northeastern region coincided with the dominance of the XBB*+F486P+F456L lineages, unlike the other regions that follow the national pattern. While the BQ.1* lineages dominated the Omicron wave in late 2022 in most Brazilian regions, SARS-CoV-2 cases in the Northern region were driven by a diverse mixture of BA.5-descendent lineages (BE.9, BQ.1*, and DL.1) and the peaks of both the BE.9/BQ.1*/DL.1 and XBB* waves were roughly similar and considerably lower than the peaks observed during the previous BA.1* and BA.5* waves. Finally, individuals under 15 years old exhibit a higher number of SARI cases during the XBB* dissemination period (n = 2,655) compared to what was observed during the BQ.1* dissemination period (n = 2,090), in contrast to all other age groups that displayed a significant decrease in SARI cases (50-60%). Additionally, the number of SARI cases during the XBB* wave (February to July 2023) in individuals under 15 years old (n = 2,655, 23%) exceeded that of the age groups 15-49 (n = 1,276, 11%) and 50-64 (n = 1,420, 12%) and only fell below the count in individuals over 64 years old (n = 6,348, 54%). These data support that the dissemination of XBB* lineages in Brazil did not bring up a significant worsening of the SARS-CoV-2 epidemiological scenario in the country. It’s noteworthy, however, that the incidence of SARI cases during the XBB* wave was higher among children (under 15) compared to young or middle age adults (aged 15-64), in contrast to previous Omicron waves.

## DISCUSSION

Our study reveals that the SARS-CoV-2 XBB epidemic in Brazil in 2023 was driven by multiple introductions followed by dissemination and evolution of different local lineages. The local evolution of XBB lineages in Brazil was characterized by the stepwise accumulation of mutations in the RBD of the Spike protein that increase the immune escape and affinity to ACE2. It showed remarkable convergence with the molecular evolution pathway followed by other Omicron XBB lineages spreading in different countries worldwide (1, 6).

After being first detected in Brazil in October 2022, XBB* represented a small fraction (<5%) of the country’s SARS-CoV-2 circulating lineages until early 2023 and the first XBB lineages to become prevalent in Brazil were those harboring the mutation S:F486P. The major XBB*+F486P lineages that arose and spread within Brazil were the XBB.1.5.86, XBB.1.5.102, XBB.1.18.1, and XBB.1.4.2. These Brazilian XBB*+F486P lineages probably started to be disseminated in Brazil between mid-October and late December 2022. They became the prevalent SARS-CoV-2 variants between January and February 2023, coinciding with the takeover of XBB* variants in Brazil. On a global scale, the XBB*+F486P genomes became the majority (∼50%) among SARS-CoV-2 genomes submitted to GISAID ([https://gisaid.org/](https://gisaid.org/)) by February 2023, reached a peak of prevalence in May (∼85%) and gradually decreased after that. These findings corroborate that XBB lineages were able to substitute the previously dominant BQ.1\*/BE\* lieages circulating in Brazil and other Western countries only after acquiring the S:F486P mutation that significantly increased the ACE2 binding affinity concerning ancestral XBB variants (21, 23–27). Our findings support that the average *RRI* of the XBB*+F486P lineages in Brazil was about 1.24 times higher than that of other cocirculating lineages, such as BQ.1* and ancestral XBB lineages, which is in line with the described average growth advantage of the XBB.1.5 lineage *w.r.t* XBB.1 and BQ.1* lineages circulating in the USA (1.23-1.26 times-higher) (24).

The following step in the local evolution of XBB variants in Brazil was the acquisition of mutation S:F456L. A new Brazilian lineage carrying the mutation S:F456L, designated as FE.1*, arose in December 2022 from the XBB*+F486P lineage XBB.1.18.1. This lineage was first detected in Brazil in late January 2023 and gradually replaced XBB*+F486P lineages, becoming the most prevalent XBB variant by May 2023. The sequential accumulation of Spike mutations F486P and F456L in the Brazilian lineages XBB.1.18.1 and FE.1* reproduced the pattern observed in other XBB lineages that spread worldwide, such as the lineages XBB.1.9.2 and EG.5*. On a global scale, XBB* genomes carrying mutations S:F486P+F456L were first collected in November 2022 but displayed limited frequency (<5%) until May 2023. Subsequently, a fast surge in their prevalence propelled them to near-majority status (40%) by August 2023. The most prevalent XBB*+F486P+F456L lineage globally was the EG.5* (∼55%), followed by FL* (∼10%), FE.1* (∼10%), and multiple (*n =* 90) others XBB* lineages. Mutation S:F456L increases immune escape while reducing the affinity to ACE2 (17, 31, 33–35). Our findings support that the average *RRI* of the XBB*+F486P+F456L lineages circulating in Brazil was 1.07 times higher than that of XBB*+F486P variants. Similar values of growth advantage were described for the FE.1* (1.07-1.12 times higher) and EG.5.1 (1.15-1.23 times higher) lineages concerning the XBB.1.5 lineage in different countries from Asia (China, Japan, Singapore and South Korea) and North America (Canada and USA) (31).

Between mid-February and early March 2023, XBB.1.5-descendant lineages carrying the Spike mutations F456L plus L455F arose in Brazil. This group’s most prevalent Brazilian lineage was the XBB.1.5.70, along with its GK.1* descendants’ lineages, followed by the lineages JD.1*. These lineages were detected in Brazil in March 2023 and surpassed other XBB* variants by June-July 2023. While mutations F456L and L455 individually reduce Spike’s affinity for the ACE2 receptor, their simultaneous presence within the XBB.1.5 backbone collectively confers both increased ACE2 binding and immune escape properties (41, 64). Our analysis found that the RRI of the XBB*+F486P+F456L+L455F lineages in Brazil was 1.19 and 1.11 times higher than that of XBB*+F486P and XBB*+F486P+F456L lineages, respectively, providing the first estimates of growth advantage for these mutants. As of August 2023, XBB* genomes carrying mutations S:F486P+F456L+L455F still displayed a low global prevalence (∼5%) and most genomes were assigned to the XBB.1.5.70+GK* lineages in South America, North America, Europe and Africa (60-80%) and to the HK.3 (EG.5.1.1.3) lineage in Asia (∼50%). However, considering the ability of XBB*+F486P+F456L+L455F lineages to outcompete other XBB* lineages in Brazil and the convergent patterns of XBB evolution in different settings, it is reasonable to anticipate both the global increase in prevalence of these triply mutant lineages in the upcoming weeks and the emergence of mutations F456L+L455F in multiple XBB+F486P backbones.

The pronounced antigenic distance of XBB* lineages to the ancestral B.1 and other Omicron subvariants significantly reduced the efficacy of previous immunity to prevent XBB* infections (65–67). Thus, dissemination of XBB* lineages drove new epidemic waves in several countries (4, 30, 38, 39, 68). In Brazil, the growth advantage of XBB lineages carrying Spike mutations F486L, F486L+F456L, and F486L+F456L+L455F was sufficient to outcompete the BQ.1\*/BE\* lineages. However, their dissemination has had little epidemiological impact in terms of severe cases. SARI-COVID notifications displayed a modest increase during the spread of XBB* lineages in Brazil, much lower than during previous BA.1*, BA.5*, and BQ.1\*/BE\* waves. This scenario was not the result of a less intrinsic virulence because infection with XBB* subvariants results in similar disease severity to infection with other Omicron subvariants (4, 30, 69). The progressive reduction of severe cases across successive Omicron waves in Brazil since early 2022 indicates that population hybrid immunity has been increasingly efficient in preventing new upsurges of severe cases despite the continuous emergence of highly immune-evasive viral variants (65, 70). Repeated Omicron infections may also broaden the breadth of antibody response and increase the neutralizing titers against highly immune evasive variants (34), which is consistent with the observation that reduction is SARI cases from BQ.1\*/BE\* to XBB* waves in Brazil was higher in adults individuals than in youngest (<15 years old) ones. It is also possible that reinfections due to waning immunity were limited in the XBB* wave by the short time interval of ∼3 months from the previous BQ.1\*/BE\* wave, which was half the interval observed between Omicron waves in Brazil in 2022.

Although the diagnostic of SARS-CoV-2 cases in Brazil and most countries is no longer as exhaustive as it was during the first years of the pandemic, our findings reinforce the concept that the epidemiologic impact of SARS-CoV-2 lineage turnovers may greatly vary across different settings according to the previous immune experience of local populations. Countries such as Singapore and Malaysia did not experience a BQ.1\*/BE\* wave, and the magnitude of the XBB wave in late 2022 was similar to the previous BA.4/BA.5 wave occurring around six months earlier (4, 68). By contrast, Western countries that displayed BA.5* followed by BQ.1* waves during 2022 experienced transient surges in cases, hospital admissions, and deaths due to XBB* infections along 2023, particularly following the emergence of the EG.5* lineage, but the overall absolute numbers remained lower than during previous Omicron waves (30, 39, 39, 71, 72). These modest epidemiological consequences stemming from the substitution of BQ.1\*/BE\* lineages by XBB* in Brazil bears a resemblance to the replacement of the Gamma variant by the Delta variant during mid-2021 that also occurred without a significant upsurge of SARS-CoV-2 cases and deaths (73). Of note, the estimated average *RRI* of Delta *w.r.t* Gamma (1.32) (73) was similar to the average *RRI* of XBB* concerning other co-circulating lineages (1.24-1.48) here estimated. Thus, the epidemiological impact of the spread of a new SARS-CoV-2 variant in a given region is influenced by an intricate interplay of factors that extend beyond the virus’s transmissibility alone.

Predicting the antigenic evolution of XBB lineages in Brazil is a difficult task. By the time of writing this article, however, a new GK.1* and JD.1* sub lineages harboring additional RBD mutations were identified in Brazil since early July. The GK.1.4 lineage harbors the RBD mutation K478R; while the GK.1.2.1, GK.1.8 and JD.1.1 harbors the convergent RBD mutation A475V. Interestingly, compared with XBB.1.5, the mutant XBB.1.5+K478R exhibited a 1.7-fold higher infectivity in CaLu-3 cells and similar antibody evasion properties (17, 41). Meanwhile, the mutant XBB.1.5+L455F+F456L+A475V could evade almost all of Class 1 antibodies and is more immune evasive than XBB.1.5+L455F+F456L, but display lower ACE2 binding (74). Some fast-growing XBB lineages displayed the RBD mutation K478R in combination with F486 (such as XBB.1.16), F486P+F456L (such as XBB.1.16.6 and FL.1.5.1) or F486P+F456L+L455F (such as JF.1). Similarly, multiple XBB lineages other than GK.1* and JD.1.1 also harbors the RBD mutation A475V in combination with mutations F486P+F456L+L455F (such as FL.15.1.1, GK.3.1, GK.4, GK.8.1, GW.5.1.1, and GW.5.3.1). Thus, the XBB lineages with the quadruple RBD mutant profiles L455F+F456L+K478R+F486P and L455F+F456L+A475V+F486P might display some growth advantage compared to double or triple RBD mutant XBB lineages currently prevalent and their dissemination in Brazil and elsewhere should be closely monitored.

The major limitation of this study revolves around the availability of XXB* samples over the ten months analyzed in various Brazilian states. Temporally, the majority of XXB* genomes were collected from February to May, constituting approximately 80% of the dataset (*n =* 5,611). As we progressed into June and July, the availability of sequences dropped significantly, accounting for only about 10% (*n =* 797) of the total XBB* dataset, and only 11 out of the 27 Brazilian states were sampled during July. This time frame corresponded with the period of expansion of the XBB*+F486P+F456L+LF455F group, and the prevalence of those triply mutated XBB* variants in Brazil could thus be underestimated. From a spatial perspective, the states of Amazonas (*n =* 1,634), São Paulo (*n =* 991), Bahia (*n =* 591), and Santa Catarina (*n =* 491) contributed with the most viral genomes, surpassing the combined contribution of the remaining 23 states. This imbalanced geographic sampling hindered our ability to perform phylogeographic inferences at the sub-national level. This fact underscores the pressing need for an ongoing and comprehensive SARS-CoV-2 genomic surveillance system in Brazil that provides timely and spatially representative data on the SARS-CoV-2 genomic composition in various regions. It is important to note that as of mid-2023, this crucial data remains unavailable in several states.

In summary, our analysis shows the ongoing emergence of new XBB* lineages in Brazil during 2023, characterized by the gradual accumulation of mutations within the RBD region of the Spike protein and the steadily growing competitive edge. Mutations S:F486P+F456L+LF455F combined resulted in the XBB* variants with the highest transmissibility detected in Brazil in the first half of 2023. Notably, Brazil is the first country to have a high prevalence of XBB*+F486P+F456L+L455F lineages, and real-world transmissibility data obtained here supports the potential of this new Omicron subvariant to spread globally. Notably, despite the significant growth advantage of XBB* lineages compared to other Omicron subvariants currently circulating, their impact on the epidemiological landscape of severe SARS-CoV-2 cases in Brazil has been less pronounced than that of the previously prevalent Omicron lineages. These findings underscore the importance of sustained molecular surveillance of SARS-CoV-2 cases to ensure rapid detection of new emergent viral variants with high transmissibility, pathogenicity, or both in Brazil and also revealed the importance of considering both virologic characteristics and population immune landscapes to predict the local impact of newly emergent SARS-CoV-2 lineages.

## CONTRIBUTION

The study was conceived and designed by GB, FGN and IA. FM, RK, TG, PCR and FGN contributed to diagnostics and sequencing analysis. K.I. contributed to the instantaneous reproduction number (RRI) estimations. MG worked on the retrieval and analysis of Brazilian epidemiological data. FGN and MMS contributed to laboratory management and obtaining financial support. The bioinformatics analysis was performed by IA and SS. I.A. and G.B. wrote the first draft and all authors contributed and approved the final manuscript.

## Supporting information

Supplementary Material [[supplements/297466_file06.pdf]](pending:yes)

## Data Availability

The findings of this study are based on the analysis of 107,881 SARS-CoV-2 Brazilian genomes collected after January 1st 2022, and submitted to the GISAID database until July 31st 2023, accessible at [https://doi.org/10.55876/gis8.231024oa](https://doi.org/10.55876/gis8.231024oa). In our phylogeographic analysis, we additionally downloaded 3,132 global references from the same period available at [https://doi.org/10.55876/gis8.231024qe](https://doi.org/10.55876/gis8.231024qe).

## ACKNOWLEDGMENTS

We gratefully acknowledge all data contributors, i.e., the Authors and their Originating laboratories responsible for obtaining the specimens, and their Submitting laboratories for generating the genetic sequence and metadata and sharing via the GISAID Initiative, on which this research is based. The authors also wish to thank all the healthcare workers and scientists who have worked hard to deal with this pandemic threat. In addition, we appreciate the support of FIOCRUZ COVID-19 Genomics Surveillance Network members (**Supplementary Material**) and the Respiratory Viruses Genomic Surveillance Network of the General Laboratory Coordination (CGLab) of the Brazilian Ministry of Health (MoH) and Brazilian Central Laboratory States (LACENs). Funding support: FAPEAM F.G.N. (INICIATIVA AMAZÔNIA +10 (grant: 01.02.016301.00439/2023-70); Rede Genômica de Vigilância em Saúde - REGESAM); FAPEAM/INOVA FIOCRUZ INOVAÇÃO NA AMAZÔNIA (Chamada Pública n° 04/2022); NPI EXPAND - U.S. Agency for International Development (USAID) implemented by Palladium (7200AA19CA00015). Centers for Disease Control and Prevention (CDC Grant Award 002174); Departamento de Ciência e Tecnologia (DECIT) of the Brazilian MoH. FAPERJ (Grant number E-26/202.896/2018): G.B. FAPERJ - Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, Grant SEI-260003/019669/2022: I.A.

*   Received October 30, 2023.
*   Revision received October 30, 2023.
*   Accepted October 30, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1.Focosi D, Quiroga R, McConnell S, Johnson MC, Casadevall A. 2023. Convergent evolution in sars-cov-2 spike creates a variant soup from which new COVID-19 waves emerge. International Journal of Molecular Sciences 24:2264.
    
    

2.  2.Velavan TP, Ntoumi F, Kremsner PG, Lee SS, Meyer CG. 2023. Emergence and geographic dominance of Omicron subvariants XBB/XBB.1.5 and BF.7 – the public health challenges. International Journal of Infectious Diseases 128:307–309.
    
    

3.  3.Ngiam J, Al-Mubaarak A, Maurer-Stroh S, Tambyah P. 2023. Does the COVID-19 XBB Omicron subvariant signal the beginning of the end of the pandemic? Singapore Medical Journal :.
    
    

4.  4.Goh AXC, Chae S-R, Chiew CJ, Tang N, Pang D, Lin C, Tan KB, Lee VJ, Ho ZJM. 2023. Characteristics of the omicron XBB subvariant wave in Singapore. The Lancet 401:1261–1262.
    
    

5.  5.Tamura T, Ito J, Uriu K, Zahradnik J, Kida I, Anraku Y, Nasser H, Shofa M, Oda Y, Lytras S, Nao N, Itakura Y, Deguchi S, Suzuki R, Wang L, Begum MM, Kita S, Yajima H, Sasaki J, Sasaki-Tabata K, Shimizu R, Tsuda M, Kosugi Y, Fujita S, Pan L, Sauter D, Yoshimatsu K, Suzuki S, Asakura H, Nagashima M, Sadamasu K, Yoshimura K, Yamamoto Y, Nagamoto T, Schreiber G, Maenaka K, Ito H, Misawa N, Kimura I, Suganami M, Chiba M, Yoshimura R, Yasuda K, Iida K, Ohsumi N, Strange AP, Takahashi O, Ichihara K, Shibatani Y, Nishiuchi T, Kato M, Ferdous Z, Mouri H, Shishido K, Sawa H, Hashimoto R, Watanabe Y, Sakamoto A, Yasuhara N, Suzuki T, Kimura K, Nakajima Y, Nakagawa S, Wu J, Shirakawa K, Takaori-Kondo A, Nagata K, Kazuma Y, Nomura R, Horisawa Y, Tashiro Y, Kawai Y, Irie T, Kawabata R, Motozono C, Toyoda M, Ueno T, Hashiguchi T, Ikeda T, Fukuhara T, Saito A, Tanaka S, Matsuno K, Takayama K, Sato K. 2023. Virological characteristics of the SARS-CoV-2 XBB variant derived from recombination of two Omicron subvariants. Nature Communications 14.
    
    

6.  6.Cao Y, Jian F, Wang J, Yu Y, Song W, Yisimayi A, Wang J, An R, Chen X, Zhang N, Wang Y, Wang P, Zhao L, Sun H, Yu L, Yang S, Niu X, Xiao T, Gu Q, Shao F, Hao X, Xu Y, Jin R, Shen Z, Wang Y, Xie XS. 2022. Imprinted SARS-CoV-2 humoral immunity induces convergent Omicron RBD evolution. Nature.
    
    

7.  7.Kurhade C, Zou J, Xia H, Liu M, Chang HC, Ren P, Xie X, Shi P. 2022. Low neutralization of SARS-CoV-2 Omicron BA.2.75.2, BQ.1.1 and XBB.1 by parental mRNA vaccine or a BA.5 bivalent booster. Nature Medicine 29:344–347.
    
    

8.  8.Arora P, Cossmann A, Schulz SR, Ramos GM, Stankov MV, Jäck H-M, Behrens GMN, Pöhlmann S, Hoffmann M. 2023. Neutralisation sensitivity of the SARS-CoV-2 XBB.1 lineage. The Lancet Infectious Diseases 23:147–148.
    
    

9.  9.Dijokaite-Guraliuc A, Das R, Zhou D, Ginn HM, Liu C, Duyvesteyn HME, Huo J, Nutalai R, Supasa P, Selvaraj M, de Silva TI, Plowright M, Newman TAH, Hornsby H, Mentzer AJ, Skelly D, Ritter TG, Temperton N, Klenerman P, Barnes E, Dunachie SJ, Roemer C, Peacock TP, Paterson NG, Williams MA, Hall DR, Fry EE, Mongkolsapaya J, Ren J, Stuart DI, Screaton GR, Conlon C, Deeks A, Frater J, Gardiner S, Jams’en A, Jeffery K, Malone T, Phillips E, Kronsteiner-Dobramysl B, Abraham P, Bibi S, Lambe T, Longet S, Tipton T, Carrol M, Stafford L. 2023. Rapid escape of new SARS-CoV-2 Omicron variants from BA.2-directed antibody responses. Cell Reports 42:112271.
    
    

10. 10.Ju B, Fan Q, Liu C, Shen S, Wang M, Guo H, Zhou B, Ge X, Zhang Z. 2023. Omicron BQ.1.1 and XBB.1 unprecedentedly escape broadly neutralizing antibodies elicited by prototype vaccination. Cell Reports 42:112532.
    
    

11. 11.Miller J, Hachmann NP, Collier AY, Lasrado N, Mazurek CR, Patio RC, Powers O, Surve N, Theiler J, Korber B, Barouch DH. 2023. Substantial neutralization escape by sars-cov-2 omicron variants BQ.1.1 and XBB.1. New England Journal of Medicine 388:662–664.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2214314&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36652339&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F30%2F2023.10.30.23297466.atom) 

12. 12.Mykytyn AZ, Rosu ME, Kok A, Rissmann M, van Amerongen G, Geurtsvankessel C, de Vries RD, Munnink BBO, Smith DJ, Koopmans MPG, Lamers MM, Fouchier RAM, Haagmans BL. 2023. Antigenic mapping of emerging SARS-CoV-2 omicron variants BM.1.1.1, BQ.1.1, and XBB.1. The Lancet Microbe 4:e294–e295.
    
    

13. 13.Tan CY, Chiew CJ, Pang D, Lee VJ, Ong B, Lye DC, Tan KB. 2023. Protective immunity of SARS-CoV-2 infection and vaccines against medically attended symptomatic omicron BA.4, BA.5, and XBB reinfections in Singapore: A national cohort study. The Lancet Infectious Diseases 23:799–805.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(23)00060-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36924786&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F30%2F2023.10.30.23297466.atom) 

14. 14.Zhu K-L, Jiang X-L, Zhan B-D, Wang X-J, Xia X, Cao G-P, Sun W-K, Huang P-X, Zhang J-Z, Gao Y-L, Dai E-H, Gao H-X, Ma M-J. 2023. Durability of neutralization against Omicron subvariants after vaccination and breakthrough infection. Cell Reports 42:112075.
    
    

15. 15.Akerman A, Milogiannakis V, Jean T, Esneau C, Silva MR, Ison T, Fichter C, Lopez JA, Chandra D, Naing Z, Caguicla J, Li D, Walker G, Amatayakul-Chantler S, Roth N, Manni S, Hauser T, Barnes T, Condylios A, Yeang M, Wong M, Foster CSP, Sato K, Lee S, Song Y, Mao L, Sigmund A, Phu A, Vande More AM, Hunt S, Douglas M, Caterson I, Britton W, Sandgren K, Bull R, Lloyd A, Triccas J, Tangye S, Bartlett NW, Darley D, Matthews G, Stark DJ, Petoumenos K, Rawlinson WD, Murrell B, Brilot F, Cunningham AL, Kelleher AD, Aggarwal A, Turville SG. 2023. Emergence and antibody evasion of BQ, BA.2.75 and SARS-CoV-2 recombinant sub-lineages in the face of maturing antibody breadth at the population level. eBioMedicine 90:104545.
    
    

16. 16.He Q, Wu L, Xu Z, Wang X, Xie Y, Chai Y, Zheng A, Zhou J, Qiao S, Huang M, Shang G, Zhao X, Feng Y, Qi J, Gao GF, Wang Q. 2023. An updated atlas of antibody evasion by SARS-CoV-2 Omicron sub-variants including BQ.1.1 and XBB. Cell Reports Medicine 4:100991.
    
    

17. 17.Faraone JN, Qu P, Zheng Y-M, Carlin C, Jones D, Panchal AR, Saif LJ, Oltz EM, Gumina RJ, Liu S-L. 2023. Continued evasion of neutralizing antibody response by Omicron XBB.1.16. Cell Reports 42:113193.
    
    

18. 18.Zou J, Kurhade C, Patel S, Kitchin N, Tompkins K, Cutler M, Cooper D, Yang Q, Cai H, Muik A, Zhang Y, Lee D-Y, Şahin U, Anderson AS, Gruber WC, Xie X, Swanson KA, Shi P-Y. 2023. Neutralization of BA.4–BA.5, BA.4.6, BA.2.75.2, BQ.1.1, and XBB.1 With bivalent vaccine. New England Journal of Medicine 388:854–857.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2214916&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36734885&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F30%2F2023.10.30.23297466.atom) 

19. 19.Davis-Gardner ME, Lai L, Wali B, Samaha H, Solis D, Lee M, Porter-Morrison A, Hentenaar IT, Yamamoto F, Godbole S, Liu Y, Douek DC, Lee FE-H, Rouphael N, Moreno A, Pinsky BA, Suthar MS. 2023. Neutralization against BA.2.75.2, BQ.1.1, and XBB from mRNA Bivalent Booster. New England Journal of Medicine 388:183–185.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2214293&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36546661&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F30%2F2023.10.30.23297466.atom) 

20. 20.Wang Q, Iketani S, Li Z, Liu L, Guo Y, Huang Y, Bowen AD, Liu M, Wang M, Yu J, Valdez R, Lauring AS, Sheng Z, Wang HH, Gordon A, Liu L, Ho DD. 2023. Alarming antibody evasion properties of rising SARS-CoV-2 BQ and XBB subvariants. Cell 186:279–286.e8.
    
    

21. 21.Yue C, Song W, Wang L, Jian F, Chen X, Gao F, Shen Z, Wang Y, Wang X, Cao Y. 2023. Enhanced transmissibility of XBB.1.5 is contributed by both strong ACE2 binding and antibody evasion. Cold Spring Harbor Laboratory.
    
    

22. 22.Ito J, Suzuki R, Uriu K, Itakura Y, Zahradnik J, Kimura KT, Deguchi S, Wang L, Lytras S, Tamura T, Kida I, Nasser H, Shofa M, Begum MM, Tsuda M, Oda Y, Suzuki T, Sasaki J, Sasaki-Tabata K, Fujita S, Yoshimatsu K, Ito H, Nao N, Asakura H, Nagashima M, Sadamasu K, Yoshimura K, Yamamoto Y, Nagamoto T, Kuramochi J, Schreiber G, Suzuki S, Kato M, Ferdous Z, Mouri H, Shishido K, Misawa N, Kimura I, Kosugi Y, Lin P, Suganami M, Chiba M, Yoshimura R, Yasuda K, Iida K, Ohsumi N, Strange AP, Sauter D, Nakagawa S, Wu J, Watanabe Y, Sakamoto A, Yasuhara N, Nakajima Y, Yajima H, Shirakawa K, Takaori-Kondo A, Nagata K, Kazuma Y, Nomura R, Horisawa Y, Tashiro Y, Kawa Y, Irie T, Kawabata R, Shimizu R, Takahashi O, Ichihara K, Motozono C, Toyoda M, Ueno T, Shibatani Y, Nishiuchi T, Saito A, Matsuno K, Takayama K, Hashiguchi T, Tanaka S, Fukuhara T, Ikeda T, Sato K. 2023. Convergent evolution of SARS-CoV-2 Omicron subvariants leading to the emergence of BQ.1.1 variant. Nature Communications 14.
    
    

23. 23.Uraki R, Ito M, Kiso M, Yamayoshi S, Iwatsuki-Horimoto K, Furusawa Y, Sakai-Tagawa Y, Imai M, Koga M, Yamamoto S, Adachi E, Saito M, Tsutsumi T, Otani A, Kikuchi T, Yotsuyanagi H, Halfmann PJ, Pekosz A, Kawaoka Y. 2023. Antiviral and bivalent vaccine efficacy against an omicron XBB.1.5 isolate. The Lancet Infectious Diseases 23:402–403.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(23)00070-1&link_type=DOI) 

24. 24.Uriu K, Ito J, Zahradnik J, Fujita S, Kosugi Y, Schreiber G, Sato K. 2023. Enhanced transmissibility, infectivity, and immune resistance of the SARS-CoV-2 omicron XBB.1.5 variant. The Lancet Infectious Diseases 23:280–281.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(23)00051-8&link_type=DOI) 

25. 25.Hoffmann M, Arora P, Nehlmeier I, Kempf A, Cossmann A, Schulz SR, Morillas Ramos G, Manthey LA, Jäck H-M, Behrens GMN, Pöhlmann S. 2023. Profound neutralization evasion and augmented host cell entry are hallmarks of the fast-spreading SARS-CoV-2 lineage XBB.1.5. Cellular & Molecular Immunology 20:419–422.
    
    

26. 26.Qu P, Faraone JN, Evans JP, Zheng Y-M, Carlin C, Anghelina M, Stevens P, Fernandez S, Jones D, Panchal AR, Saif LJ, Oltz EM, Zhang B, Zhou T, Xu K, Gumina RJ, Liu S-L. 2023. Enhanced evasion of neutralizing antibody response by Omicron XBB.1.5, CH.1.1, and CA.3.1 variants. Cell Reports 42:112443.
    
    

27. 27.Addetia A, Piccoli L, Case JB, Park Y-J, Beltramello M, Guarino B, Dang H, de Melo GD, Pinto D, Sprouse K, Scheaffer SM, Bassi J, Silacci-Fregni C, Muoio F, Dini M, Vincenzetti L, Acosta R, Johnson D, Subramanian S, Saliba C, Giurdanella M, Lombardo G, Leoni G, Culap K, McAlister C, Rajesh A,  Dellota E Jr, Zhou J, Farhat N, Bohan D, Noack J, Chen A, Lempp FA, Quispe J, Kergoat L, Larrous F, Cameroni E, Whitener B, Giannini O, Cippà P, Ceschi A, Ferrari P, Franzetti-Pellanda A, Biggiogero M, Garzoni C, Zappi S, Bernasconi L, Kim MJ, Rosen LE, Schnell G, Czudnochowski N, Benigni F, Franko N, Logue JK, Yoshiyama C, Stewart C, Chu H,  Bourhy He, Schmid MA, Purcell LA, Snell G, Lanzavecchia A, Diamond MS, Corti D, Veesler D. 2023. Neutralization, effector function and immune imprinting of Omicron variants. Nature 621:592–601.
    
    

28. 28.Yamasoba D, Uriu K, Plianchaisuk A, Kosugi Y, Pan L, Zahradnik J, Ito J, Sato K. 2023. Virological characteristics of the SARS-CoV-2 omicron XBB.1.16 variant. The Lancet Infectious Diseases 23:655–656.
    
    

29. 29.Nehlmeier I, Kempf A, Arora P, Cossmann A, Dopfer-Jablonka A, Stankov MV, Schulz SR, Jäck H-M, Behrens GMN, Pöhlmann S, Hoffmann M. 2023. Host cell entry and neutralisation sensitivity of the SARS-CoV-2 XBB.1.16 lineage. Cellular & Molecular Immunology 20:969–971.
    
    

30. 30.ECDC. 2023. ECDC classifies XBB.1.5-like lineages with the amino acid change F456L as variants of interest following an increase in SARS-CoV-2 transmission in EU/EEA countries and abroad.
    
    

31. 31.Kaku Y, Kosugi Y, Uriu K, Ito J,  Hinay AA Jr, Kuramochi J, Sadamasu K, Yoshimura K, Asakura H, Nagashima M, Sato K. 2023. Antiviral efficacy of the SARS-CoV-2 XBB breakthrough infection sera against omicron subvariants including EG.5. The Lancet Infectious Diseases 23:e395–e396.
    
    

32. 32.Tsujino S, Deguchi S, Nomai T, Padilla-Blanco M, Plianchaisuk A, Wang L, Begum MM, Uriu K, Mizuma K, Nao N, Kojima I, Tsubo T, Li J, Matsumura Y, Nagao M, Oda Y, Tsuda M, Anraku Y, Kita S, Yajima H, Tabata K, Guo Z,  Hinay AA Jr, Yoshimatsu K, Yamamoto Y, Nagamoto Y, Asakura H, Nagashima M, Sadamasu K, Yoshimura K, Nasser H, Jonathan M, Putri O, Kim Y, Chen L, Suzuki R, Tamura T, Maenaka K, Irie T, Matsuno K, Tanaka S, Ito J, Ikeda T, Takayama K, Zahradnik J, Hashiguchi T, Fukuhara T, Sato K. 2023. Virological characteristics of the SARS-CoV-2 Omicron EG.5.1 variant. Cold Spring Harbor Laboratory.
    
    

33. 33.Wang Q, Ye S, Zhou Z,  Song Al, Zhu X, Peng J, Liang R, Yang C, Yu X, Huang X, Yu J, Qiu Y, Ge X. 2022. Key mutations in the spike protein of SARSCoV2 affecting neutralization resistance and viral internalization. Journal of Medical Virology 95.
    
    

34. 34.Yisimayi A, Song W, Wang J, Jian F, Yu Y, Chen X, Xu Y, Yang S, Niu X, Xiao T, Wang J, Zhao L, Sun H, An R, Zhang N, Wang Y, Wang P, Yu L, Lv Z, Gu Q, Shao F, Jin R, Shen Z, Xie XS, Wang Y, Cao Y. 2023. Repeated Omicron exposures override ancestral SARS-CoV-2 immune imprinting. Cold Spring Harbor Laboratory.
    
    

35. 35.Zhang L, Kempf A, Nehlmeier I, Cossmann A, Dopfer-Jablonka A, Stankov MV, Schulz SR, Jack H-M, Behrens GMN, Pohlmann S, Hoffmann M. 2023. Neutralisation sensitivity of SARS-CoV-2 lineages EG.5.1 and XBB.2.3. The Lancet Infectious Diseases 23:e391–e392.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(23)00547-9&link_type=DOI) 

36. 36.Wang Q, Guo Y, Zhang RM, Ho J, Mohri H, Valdez R, Manthei DM, Gordon A, Liu L, Ho DD. 2023. Antibody neutralisation of emerging SARS-CoV-2 subvariants: EG.5.1 and XBC.1.6. The Lancet Infectious Diseases 23:e397–e398.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(23)00555-8&link_type=DOI) 

37. 37.WHO. EG.5 Initial Risk Evaluation, 9 August 2023.
    
    

38. 38.Dyer O. 2023. Covid-19: Infections climb globally as EG.5 variant gains ground. BMJ p1900.
    
    

39. 39.Parums DV. 2023. Editorial: A Rapid Global Increase in COVID-19 is Due to the Emergence of the EG.5 (Eris) Subvariant of Omicron SARS-CoV-2. Medical Science Monitor 29.
    
    

40. 40.Chen C, Nadeau S, Yared M, Voinov P, Xie N, Roemer C, Stadler T. 2021. CoV-Spectrum: Analysis of globally shared SARS-CoV-2 data to identify and characterize new variants. Bioinformatics 38:1735–1737.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btab856&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34954792&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F30%2F2023.10.30.23297466.atom) 

41. 41.Jian F, Yang S, Yu Y, Song W, Yisimayi A, Chen X, Xu Y, Wang P, Yu L, Wang J, Liu L, Niu X, Wang J, Xiao T, An R, Wang Y, Gu Q, Shao F, Jin R, Shen Z, Wang Y, Cao Y. 2023. Convergent evolution of SARS-CoV-2 XBB lineages on receptor-binding domain 455-456 enhances antibody evasion and ACE2 binding. Cold Spring Harbor Laboratory.
    
    

42. 42.FIOCRUZ. SARS-CoV-2 Genomic Surveillance in Brazil. FIOCRUZ.
    
    

43. 43.Resende PC, Delatorre E, Gräf T, Mir D, Motta FC, Appolinario LR, da Paixão ACD, da Mendonça AC F, Ogrzewalska M, Caetano B, Wallau GL, Docena C assia, dos Santos MC,  de Almeida Ferreira J,  Sousa Junior EC, da Silva SP, Fernandes SB, Vianna LA, da Souza L C, Ferro JFG, Nardy VB, Santos CA, Riediger I, do Carmo Debur M,  Croda Júlio, Oliveira WK,  Abreu Ae, Bello G, Siqueira MM. 2021. Evolutionary dynamics and dissemination pattern of the sars-cov-2 lineage B.1.1.33 during the early pandemic phase in Brazil. Frontiers in Microbiology 11.
    
    

44. 44.Dezordi FZ, da Neto AM S,  de Campos Tulio L, Jeronimo PMC, Aksenen CF, Almeida SP, Wallau GL. 2022. ViralFlow: A versatile automated workflow for sars-cov-2 genome assembly, lineage assignment, mutations and intrahost variant detection. Viruses 14:217.
    
    

45. 45.Aksamentov I, Roemer C, Hodcroft E, Neher R. 2021. Nextclade: Clade assignment, mutation calling and quality control for viral genomes. Journal of Open Source Software 6:3773.
    
    

46. 46.O’Toole Áine, Scher E, Underwood A, Jackson B, Hill V, McCrone JT, Colquhoun R, Ruis C, Abu-Dahab K, Taylor B, Yeats C, du Plessis L, Maloney D, Medd N, Attwood SW, Aanensen DM, Holmes EC, Pybus OG, Rambaut A. 2021. Assignment of epidemiological lineages in an emerging pandemic using the pangolin tool. Virus Evolution 7.
    
    

47. 47.Li W, Jaroszewski L, Godzik A. 2001. Clustering of highly homologous sequences to reduce the size of large protein databases. Bioinformatics 17:282–283.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/17.3.282&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11294794&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F30%2F2023.10.30.23297466.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000168053800009&link_type=ISI) 

48. 48.Katoh K, Rozewicki J, Yamada KD. 2017. MAFFT online service: Multiple sequence alignment, interactive sequence choice and visualization. Briefings in Bioinformatics 20:1160–1166.
    
    

49. 49.Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, Lanfear R. 2020. IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. Molecular Biology and Evolution 37:1530– 1534.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msaa015&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32556291&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F30%2F2023.10.30.23297466.atom) 

50. 50.Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. 2017. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nature Methods 14:587–589.
    
    

51. 51.Rambaut A, Lam TT, Max Carvalho L, Pybus OG. 2016. Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen). Virus Evolution 2:vew007.
    
    [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%2F2023%2F10%2F30%2F2023.10.30.23297466.atom) 

52. 52.Suchard MA, Lemey P, Baele G, Ayres DL, Drummond AJ, Rambaut A. 2018. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evolution 4.
    
    

53. 53.Drummond AJ. 2005. Bayesian coalescent inference of past population dynamics from molecular sequences. Molecular Biology and Evolution 22:1185–1192.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msi103&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15703244&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F30%2F2023.10.30.23297466.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000228700200004&link_type=ISI) 

54. 54.Lemey P, Rambaut A, Drummond AJ, Suchard MA. 2009. Bayesian phylogeography finds its roots. PLoS Computational Biology 5:e1000520.
    
    

55. 55.Rambaut A, Drummond AJ, Xie D, Baele G, Suchard MA. 2018. Posterior summarization in bayesian phylogenetics using tracer 1.7. Systematic Biology 67:901–904.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/sysbio/syy032&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29718447&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F30%2F2023.10.30.23297466.atom) 

56. 56.Wang L-G, Lam TT-Y, Xu S, Dai Z, Zhou L, Feng T, Guo P, Dunn CW, Jones BR, Bradley T, Zhu H, Guan Y, Jiang Y, Yu G. 2019. Treeio: An R package for phylogenetic tree input and output with richly annotated and associated data. Molecular Biology and Evolution 37:599–603.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msz240&link_type=DOI) 

57. 57.Yu G. 2020. Using ggtree to Visualize Data on TreeLike Structures. Current Protocols in Bioinformatics 69.
    
    

58. 58.Ito K, Piantham C, Nishiura H. 2021. Predicted dominance of variant Delta of SARS-CoV-2 before Tokyo Olympic Games, Japan, July 2021. Eurosurveillance 26.
    
    

59. 59.Fraser C. 2007. Estimating individual and household reproduction numbers in an emerging epidemic. PLoS ONE 2:e758.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0000758&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17712406&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F30%2F2023.10.30.23297466.atom) 

60. 60.Piantham C, Ito K. 2022. Predicting the trajectory of replacements of sars-cov-2 variants using relative reproduction numbers. Viruses 14:2556.
    
    

61. 61.Ito K, Piantham C, Nishiura H. 2022. Estimating relative generation times and relative reproduction numbers of Omicron BA.1 and BA.2 with respect to Delta in Denmark. Cold Spring Harbor Laboratory.
    
    

62. 62.Park SW, Sun K, Abbott S, Sender R, Bar-on YM, Weitz JS, Funk S, Grenfell BT, Backer JA, Wallinga J, Viboud C, Dushoff J. 2023. Inferring the differences in incubation-period and generation-interval distributions of the Delta and Omicron variants of SARS-CoV-2. Proceedings of the National Academy of Sciences 120.
    
    

63. 63.Held L, Bove DSes. 2020. Likelihood and bayesian inference: With applications in biology and medicine. Springer Nature.
    
    

64. 64.Qu P, Xu K, Faraone JN, Goodarzi N, Zheng Y-M, Carlin C, Bednash JS, Horowitz JC, Mallampalli RK, Saif LJ, Oltz EM, Jones D, Gumina RJ, Liu S-L. 2023. Immune Evasion, Infectivity, and Fusogenicity of SARS-CoV-2 Omicron BA.2.86 and FLip Variants. Cold Spring Harbor Laboratory.
    
    

65. 65.Tan CY, Chiew CJ, Pang D, Lee VJ, Ong B, Lye DC, Tan KB. 2023. Protective immunity of SARS-CoV-2 infection and vaccines against medically attended symptomatic omicron BA.4, BA.5, and XBB reinfections in Singapore: A national cohort study. The Lancet Infectious Diseases 23:799–805.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(23)00060-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36924786&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F30%2F2023.10.30.23297466.atom) 

66. 66.Shrestha NK, Burke PC, Nowacki AS, Simon JF, Hagen A, Gordon SM. 2023. Effectiveness of the coronavirus disease 2019 bivalent vaccine. Open Forum Infectious Diseases 10.
    
    

67. 67.Lewnard JA, Hong V, Kim JS, Shaw SF, Lewin B, Takhar H, Lipsitch M, Tartof SY. 2023. Increased vaccine sensitivity of an emerging SARS-CoV-2 variant. Nature Communications 14.
    
    

68. 68.Lee TY, Lim WF, Ang GY, Yu CY. 2023. Genomic surveillance of sars-cov-2 in Malaysia during the era of endemic COVID-19. Life 13:1644.
    
    

69. 69.Klein EY, Fall A, Norton JM, Eldesouki RE, Abdullah O, Han L, Yunker M, Mostafa HH. 2023. Severity outcomes associated with SARS-CoV-2 XBB variants, an observational analysis. Journal of Clinical Virology 165:105500.
    
    

70. 70.Lasrado N, Barouch DH. 2023. SARS-CoV-2 hybrid immunity: The best of both worlds. The Journal of Infectious Diseases.
    
    

71. 71.CDC. 2020. COVID Data Tracker. [https://covid.cdc.gov/covid-data-tracker/#trends\_weeklyhospitaladmissions\_select\_00](https://covid.cdc.gov/covid-data-tracker/#trends_weeklyhospitaladmissions_select_00). Retrieved 12 October 2023.
    
    

72. 72.GOV.UK. England Summary. [https://coronavirus.data.gov.uk/](https://coronavirus.data.gov.uk/). Retrieved 12 October 2023.
    
    

73. 73.Arantes I, Bello G, Nascimento V, Souza V, da Silva A, Silva D, Nascimento F, Mejía M,  Brandão M Júlia, Gonçalves L, Silva G, da Costa CF, Abdalla L, Santos João H, Ramos TCA, Piantham C, Ito K, Siqueira MM, Resende PC, Wallau GL, Delatorre E, Gräf T, Naveca FG. 2023. Comparative epidemic expansion of SARS-CoV-2 variants Delta and Omicron in the Brazilian State of Amazonas. Nature Communications 14.
    
    

74. 74.Cao YR. X. [https://twitter.com/yunlong\_cao/status/1714694257175654415](https://twitter.com/yunlong_cao/status/1714694257175654415). Retrieved 23 October 2023.

 [1]: /embed/graphic-1.gif
 [2]: /embed/graphic-2.gif
 [3]: /embed/graphic-3.gif