Role of Q675H Mutation in Improving SARS-CoV-2 Spike Interaction with the Furin Binding Pocket

Genotype screening was implemented in Italy and showed a significant prevalence of new SARS-CoV-2 mutants carrying Q675H mutation, near the furin cleavage site of spike protein. Currently, this mutation, which is expressed on different SARS-CoV-2 lineages circulating worldwide, has not been thoughtfully investigated. Therefore, we performed phylogenetic and biocomputational analysis to better understand SARS-CoV-2 Q675H mutants’ evolutionary relationships with other circulating lineages and Q675H function in its molecular context. Our studies reveal that Q675H spike mutation is the result of parallel evolution because it arose independently in separate evolutionary clades. In silico data show that the Q675H mutation gives rise to a hydrogen-bonds network in the spike polar region. This results in an optimized directionality of arginine residues involved in interaction of spike with the furin binding pocket, thus improving proteolytic exposure of the viral protein. Furin was predicted to have a greater affinity for Q675H than Q675 substrate conformations. As a consequence, Q675H mutation could confer a fitness advantage to SARS-CoV-2 by promoting a more efficient viral entry. Interestingly, here we have shown that Q675H spike mutation is documented in all the VOCs. This finding highlights that VOCs are still evolving to enhance viral fitness and to adapt to the human host. At the same time, it may suggest Q675H spike mutation involvement in SARS-CoV-2 evolution.     

SARS-CoV-2 Variant of Concern 202012/01 Has about Twofold Replicative Advantage and Acquires Concerning Mutations

The novel SARS-CoV-2 Variant of Concern (VOC)-202012/01 (also known as B.1.1.7), first collected in United Kingdom on 20 September 2020, is a rapidly growing lineage that in January 2021 constituted 86% of all SARS-CoV-2 genomes sequenced in England. The VOC has been detected in 40 out of 46 countries that reported at least 50 genomes in January 2021. We have estimated that the replicative advantage of the VOC is in the range 1.83–2.18 [95% CI: 1.71–2.40] with respect to the 20A.EU1 variant that dominated in England in November 2020, and in range 1.65–1.72 [95% CI: 1.46–2.04] in Wales, Scotland, Denmark, and USA. As the VOC strain will likely spread globally towards fixation, it is important to monitor its molecular evolution. We have estimated growth rates of expanding mutations acquired by the VOC lineage to find that the L18F substitution in spike has initiated a fast growing VOC substrain. The L18F substitution is of significance because it has been found to compromise binding of neutralizing antibodies. Of concern are immune escape mutations acquired by the VOC: E484K, F490S, S494P (in the receptor binding motif of spike) and Q677H, Q675H (in the proximity of the polybasic cleavage site at the S1/S2 boundary). These mutants may hinder efficiency of existing vaccines and expand in response to the increasing after-infection or vaccine-induced seroprevalence.     

Emergence and Spread of SARS-CoV-2 Lineages B.1.1.7 and P.1 in Italy

Italy’s second wave of SARS-CoV-2 has hit hard, with more than three million cases and over 100,000 deaths, representing an almost ten-fold increase in the numbers reported by August 2020. Herein, we present an analysis of 6515 SARS-CoV-2 sequences sampled in Italy between 29 January 2020 and 1 March 2021 and show how different lineages emerged multiple times independently despite lockdown restrictions. Virus lineage B.1.177 became the dominant variant in November 2020, when cases peaked at 40,000 a day, but since January 2021 this is being replaced by the B.1.1.7 ‘variant of concern’. In addition, we report a sudden increase in another documented variant of concern—lineage P.1—from December 2020 onwards, most likely caused by a single introduction into Italy. We again highlight how international importations drive the emergence of new lineages and that genome sequencing should remain a top priority for ongoing surveillance in Italy.     

SARS-CoV-2 Lineage A.27: New Data from African Countries and Dynamics in the Context of the COVID-19 Pandemic

SARS-CoV-2 is constantly evolving with lineages emerging and others eclipsing. Some lineages have an important epidemiological impact and are known as variants of interest (VOIs), variants under monitoring (VUMs) or variants of concern (VOCs). Lineage A.27 was first defined as a VUM since it holds mutations of concern. Here, we report additional lineage A.27 data and sequences from five African countries and describe the molecular characteristics, and the genetic history of this lineage worldwide. Based on the new sequences investigated, the most recent ancestor (tMRCA) of lineage A.27 was estimated to be from April 2020 from Niger. It then spread to Europe and other parts of the world with a peak observed between February and April 2021. The detection rate of A.27 then decreased with only a few cases reported during summer 2021. The phylogenetic analysis revealed many sub-lineages. Among them, one was defined by the substitution Q677H in the spike (S) gene, one was defined by the substitution D358N in the nucleoprotein (N) gene and one was defined by the substitution A2143V in the ORF1b gene. This work highlights the importance of molecular characterization and the timely submission of sequences to correctly describe the circulation of particular strains in order to be proactive in monitoring the pandemic.     

A Potential SARS-CoV-2 Variant of Interest (VOI) Harboring Mutation E484K in the Spike Protein Was Identified within Lineage B.1.1.33 Circulating in Brazil

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epidemic in Brazil was dominated by two lineages designated as B.1.1.28 and B.1.1.33. The two SARS-CoV-2 variants harboring mutations at the receptor-binding domain of the Spike (S) protein, designated as lineages P.1 and P.2, evolved from lineage B.1.1.28 and are rapidly spreading in Brazil. Lineage P.1 is considered a Variant of Concern (VOC) because of the presence of multiple mutations in the S protein (including K417T, E484K, N501Y), while lineage P.2 only harbors mutation S:E484K and is considered a Variant of Interest (VOI). On the other hand, epidemiologically relevant B.1.1.33 deriving lineages have not been described so far. Here we report the identification of a new SARS-CoV-2 VOI within lineage B.1.1.33 that also harbors mutation S:E484K and was detected in Brazil between November 2020 and February 2021. This VOI displayed four non-synonymous lineage-defining mutations (NSP3:A1711V, NSP6:F36L, S:E484K, and NS7b:E33A) and was designated as lineage N.9. The VOI N.9 probably emerged in August 2020 and has spread across different Brazilian states from the Southeast, South, North, and Northeast regions.     

A Comprehensive Molecular Epidemiological Analysis of SARS-CoV-2 Infection in Cyprus from April 2020 to January 2021: Evidence of a Highly Polyphyletic and Evolving Epidemic

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) resulted in an extraordinary global public health crisis. In early 2020, Cyprus, among other European countries, was affected by the SARS-CoV-2 epidemic and adopted lockdown measures in March 2020 to limit the initial outbreak on the island. In this study, we performed a comprehensive retrospective molecular epidemiological analysis (genetic, phylogenetic, phylodynamic and phylogeographic analyses) of SARS-CoV-2 isolates in Cyprus from April 2020 to January 2021, covering the first ten months of the SARS-CoV-2 infection epidemic on the island. The primary aim of this study was to assess the transmissibility of SARS-CoV-2 lineages in Cyprus. Whole SARS-CoV-2 genomic sequences were generated from 596 clinical samples (nasopharyngeal swabs) obtained from community-based diagnostic testing centers and hospitalized patients. The phylogenetic analyses revealed a total of 34 different lineages in Cyprus, with B.1.258, B.1.1.29, B.1.177, B.1.2, B.1 and B.1.1.7 (designated a Variant of Concern 202012/01, VOC) being the most prevalent lineages on the island during the study period. Phylodynamic analysis showed a highly dynamic epidemic of SARS-CoV-2 infection, with three consecutive surges characterized by specific lineages (B.1.1.29 from April to June 2020; B.1.258 from September 2020 to January 2021; and B.1.1.7 from December 2020 to January 2021). Genetic analysis of whole SARS-CoV-2 genomic sequences of the aforementioned lineages revealed the presence of mutations within the S protein (L18F, ΔH69/V70, S898F, ΔY144, S162G, A222V, N439K, N501Y, A570D, D614G, P681H, S982A and D1118H) that confer higher transmissibility and/or antibody escape (immune evasion) upon the virus. Phylogeographic analysis indicated that the majority of imports and exports were to and from the United Kingdom (UK), although many other regions/countries were identified (southeastern Asia, southern Europe, eastern Europe, Germany, Italy, Brazil, Chile, the USA, Denmark, the Czech Republic, Slovenia, Finland, Switzerland and Pakistan). Taken together, these findings demonstrate that the SARS-CoV-2 infection epidemic in Cyprus is being maintained by a continuous influx of lineages from many countries, resulting in the establishment of an ever-evolving and polyphyletic virus on the island.     

Molecular Evolution and Epidemiological Characteristics of SARS COV-2 in (Northwestern) Poland

The emergence of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) evolved into a worldwide outbreak, with the first Polish cases in February/March 2020. This study aimed to investigate the molecular epidemiology of the circulating virus lineages between March 2020 and February 2021. We performed variant identification, spike mutation pattern analysis, and phylogenetic and evolutionary analyses for 1106 high-coverage whole-genome sequences, implementing maximum likelihood, multiple continuous-time Markov chain, and Bayesian birth–death skyline models. For time trends, logistic regression was used. In the dataset, virus B.1.221 lineage was predominant (15.37%), followed by B.1.258 (15.01%) and B.1.1.29 (11.48%) strains. Three clades were identified, being responsible for 74.41% of infections over the analyzed period. Expansion in variant diversity was observed since September 2020 with increasing frequency of the number in spike substitutions, mainly H69V70 deletion, P681H, N439K, and S98F. In population dynamics inferences, three periods with exponential increase in infection were observed, beginning in March, July, and September 2020, respectively, and were driven by different virus clades. Additionally, a notable increase in infections caused by the B.1.1.7 lineage since February 2021 was noted. Over time, the virus accumulated mutations related to optimized transmissibility; therefore, faster dissemination is reflected by the second wave of epidemics in Poland.     

In-Flight Transmission of a SARS-CoV-2 Lineage B.1.617.2 Harbouring the Rare S:E484Q Immune Escape Mutation

We describe a flight-associated infection scenario of seven individuals with a B.1.617.2 (Delta) lineage, harbouring an S:E484Q point mutation. In Sweden, at least 10% of all positive SARS-CoV-2 samples were sequenced in each county; the B.1.717.2 + S:E484Q combination was not detected in Sweden before and was imported within the scenario described in this report. The high transmission rate of the delta lineage combined with the S:E484Q mutation, associated with immune escape in other lineages, makes this specific genetic combination a possible threat to the global fight against the COVID-19 pandemic. Even within the Omicron wave, the B.1.617.2 + S:E484Q variant appeared in community samples in Sweden, as it seems that this combination has an evolutionary gain compared to other B.1.617.2 lineages. The here described genomic combination was not detectable with the common fasta file-based Pango-lineage analysis, hence increasing the probability of the true global prevalence to be higher.     

Phylodynamic Dispersal of SARS-CoV-2 Lineages Circulating across Polish–German Border Provinces

Introduction: The emergence of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has evolved into a worldwide outbreak, with significant molecular evolution over time. Large-scale phylodynamic studies allow to map the virus spread and inform preventive strategies. Aim: This study investigates the extent of binational dispersal and dynamics of SARS-CoV-2 lineages between seven border provinces of the adjacent countries of Poland and Germany to reconstruct SARS-CoV-2 transmission networks. Methods: Following three pandemic waves from March 2020 to the end of May 2021, we analysed a dataset of 19,994 sequences divided into B.1.1.7|Alpha and non-Alpha lineage groups. We performed phylogeographic analyses using the discrete diffusion models to identify the pathways of virus spread. Results: Based on population dynamics inferences, in total, 673 lineage introductions (95% HPD interval 641–712) for non-Alpha and 618 (95% HPD interval 599–639) for B.1.1.7|Alpha were identified in the area. For non-Alpha lineages, 5.05% binational, 86.63% exclusively German, and 8.32% Polish clusters were found, with a higher frequency of international clustering observed for B.1.1.7|Alpha (13.11% for binational, 68.44% German and 18.45% Polish, p < 0.001). We identified key transmission hubs for the analysed lineages, namely Saxony, West Pomerania and Lower Silesia. Conclusions: Clustering patterns between Poland and Germany reflect the viral variant transmission dynamics at the international level in the borderline area. Tracing the spread of the virus between two adjacent large European countries may provide a basis for future intervention policies in cross-border cooperation efforts against the spread of the pandemics.     

CoronaVac and ChAdOx1 Vaccination and Gamma Infection Elicited Neutralizing Antibodies against the SARS-CoV-2 Delta Variant

The emergence of new SARS-CoV-2 variants represents a constant threat to world public health. The SARS-CoV-2 Delta variant was identified in late 2020 in India; since then, it has spread to many other countries, replacing other predominant lineages and raising concerns about vaccination efficiency. We evaluated the sensitivity of the Delta variant to antibodies elicited by COVID-19 vaccinated (CoronaVac and ChAdOx1) and convalescent individuals previously infected by earlier lineages and by the Gamma variant. No reduction in the neutralizing efficacy of the Delta variant was observed when compared to B lineage and a reduced neutralization was observed for the Gamma variant. Our results indicate that neutralization of the Delta variant is not compromised in individuals vaccinated by CoronaVac or ChAdOx1; however, a reduction in neutralization efficacy is expected for individuals infected by the Gamma variant, highlighting the importance of continuous vaccination even for previously infected individuals.     

Genetic Analysis of SARS-CoV-2 Variants in Mexico during the First Year of the COVID-19 Pandemic

During the first year of the SARS-CoV-2 pandemic in Mexico, more than two million people were infected. In this study, we analyzed full genome sequences from 27 February 2020 to 28 February 2021 to characterize the geographical and temporal distribution of SARS-CoV-2 lineages and identify the most common circulating lineages during this period. We defined six different geographical regions with particular dynamics of lineage circulation. The Northeast and Northwest regions were the ones that exhibited the highest lineage diversity, while the Central south and South/Southeast regions presented less diversity with predominance of a certain lineage. Additionally, by late February 2021, lineage B.1.1.519 represented more than 89% of all circulating lineages in the country.     

Genotype Uniformity,Wild Bird-to-Poultry Transmissions,and Farm-to-Farm Carryover during the Spread of the Highly Pathogenic Avian Influenza H5N8 in the Czech Republic in 2021

In 2020–2021, the second massive dissemination of a highly pathogenic avian influenza of the H5Nx subtype occurred in Europe. During this period, the virus caused numerous outbreaks in poultry, including in the Czech Republic. In the present study, we provide an insight into the genetic variability of the Czech/2021 (CZE/2021) H5N8 viruses to determine the relationships between strains from wild and domestic poultry and to infer transmission routes between the affected flocks of commercial poultry. For this purpose, whole genome sequencing and phylogenetic analysis of 70 H5N8 genomes representing 79.7% of the cases were performed. All CZE/2021 H5N8 viruses belonged to the 2.3.4.4b H5 lineage and circulated without reassortment, retaining the A/chicken/Iraq/1/2020 H5N8-like genotype constellation. Phylogenetic analysis suggested the frequent local transmission of H5N8 from wild birds to backyard poultry and extensive spread among commercial poultry farms. In addition, the analysis suggested one cross-border transmission event. Indirect transmission via contaminated materials was considered the most likely source of infection. Improved biosecurity and increased collaboration between field veterinarians and the laboratory are essential to limit the local spread of the virus and to reveal and interrupt critical routes of infection.     

The SARS-CoV-2 Lambda variant and its neutralisation efficiency following vaccination with Comirnaty,Israel, April to June 2021

The SARS-CoV-2 Lambda (Pango lineage designation C.37) variant of interest, initially identified in Peru, has spread to additional countries. First detected in Israel in April 2021 following importations from Argentina and several European countries, the Lambda variant infected 18 individuals belonging to two main transmission chains without further spread. Micro-neutralisation assays following Comirnaty (BNT162b2 mRNA, BioNTech-Pfizer) vaccination demonstrated a significant 1.6-fold reduction in neutralising titres compared with the wild type virus, suggesting increased susceptibility of vaccinated individuals to infection.     

Genomic Surveillance of Circulating SARS-CoV-2 in South East Italy: A One-Year Retrospective Genetic Study

In order to provide insights into the evolutionary and epidemiological viral dynamics during the current COVID-19 pandemic in South Eastern Italy, a total of 298 genomes of SARS-CoV-2 strains collected in the Apulia and Basilicata regions, between March 2020 and January 2021, were sequenced. The genomic analysis performed on the draft genomes allowed us to assign the genetic clades and lineages of belonging to each sample and provide an overview of the main circulating viral variants. Our data showed the spread in Apulia and Basilicata of SARS-CoV-2 variants which have emerged during the second wave of infections and are being currently monitored worldwide for their increased transmission rate and their possible impact on vaccines and therapies. These results emphasize the importance of genome sequencing for the epidemiological surveillance of the new SARS-CoV-2 variants’ spread.     

Genomic Epidemiology of SARS-CoV-2 in Tocantins State and the Diffusion of P.1.7 and AY.99.2 Lineages in Brazil

Tocantins is a state in the cross-section between the Central-West, North and Northeast regions of Brazilian territory; it is a gathering point for travelers and transportation from the whole country. In this study, 9493 genome sequences, including 241 local SARS-CoV-2 samples (collected from 21 December 2020, to 16 December 2021, and sequenced in the MinION platform) were analyzed with the following aims: (i) identify the relative prevalence of SARS-CoV-2 lineages in the state of Tocantins; (ii) analyze them phylogenetically against global SARS-CoV-2 sequences; and (iii) hypothesize the viral dispersal routes of the two most abundant lineages found in our study using phylogenetic and phylogeographic approaches. The performed analysis demonstrated that the majority of the strains sequenced during the period belong to the Gamma P.1.7 (32.4%) lineage, followed by Delta AY.99.2 (27.8%), with the first detection of VOC Omicron. As expected, there was mainly a dispersion of P.1.7 from the state of São Paulo to Tocantins, with evidence of secondary spreads from Tocantins to Goiás, Mato Grosso, Amapá, and Pará. Rio de Janeiro was found to be the source of AY.99.2 and from then, multiple cluster transmission was observed across Brazilian states, especially São Paulo, Paraiba, Federal District, and Tocantins. These data show the importance of trade routes as pathways for the transportation of the virus from Southeast to Northern Brazil.     

From the Cover: Multiple spillovers from humans and onward transmission of SARS-CoV-2 in white-tailed deer

Many animal species are susceptible to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and could act as reservoirs; however, transmission in free-living animals has not been documented. White-tailed deer, the predominant cervid in North America, are susceptible to SARS-CoV-2 infection, and experimentally infected fawns can transmit the virus. To test the hypothesis that SARS-CoV-2 is circulating in deer, 283 retropharyngeal lymph node (RPLN) samples collected from 151 free-living and 132 captive deer in Iowa from April 2020 through January of 2021 were assayed for the presence of SARS-CoV-2 RNA. Ninety-four of the 283 (33.2%) deer samples were positive for SARS-CoV-2 RNA as assessed by RT-PCR. Notably, following the November 2020 peak of human cases in Iowa, and coinciding with the onset of winter and the peak deer hunting season, SARS-CoV-2 RNA was detected in 80 of 97 (82.5%) RPLN samples collected over a 7-wk period. Whole genome sequencing of all 94 positive RPLN samples identified 12 SARS-CoV-2 lineages, with B.1.2 (n = 51; 54.5%) and B.1.311 (n = 19; 20%) accounting for ∼75% of all samples. The geographic distribution and nesting of clusters of deer and human lineages strongly suggest multiple human-to-deer transmission events followed by subsequent deer-to-deer spread. These discoveries have important implications for the long-term persistence of the SARS-CoV-2 pandemic. Our findings highlight an urgent need for a robust and proactive “One Health” approach to obtain enhanced understanding of the ecology, molecular evolution, and dissemination of SARS-CoV-2.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the cause of COVID-19 in humans, is a novel coronavirus in the genus Betacoronavirus (subgenus Sarbecovirus) (1). SARS-CoV-2 was first identified in Wuhan, China, toward the end of 2019 (2), and has caused a pandemic with more than 250 million COVID-19 cases and over 5 million deaths globally as of November 2021 (3). The virus continues to evolve and spread, with a growing concern for the emergence of new variants. SARS-CoV-2 uses the host angiotensin-converting enzyme 2 (ACE-2) receptor to enter cells (4). ACE-2 receptors are well conserved across vertebrate species, including humans (5), and computational analyses predict high binding affinities of SARS-CoV-2 to the ACE-2 receptor in multiple animal species, indicating potential susceptibility to infection (5). Included among these are three species of cervids: the Père David’s deer (Elaphurus davidianus), reindeer (Rangifer tarandus), and white-tailed deer (Odocoileus virginianus) (6).Global dissemination of SARS-CoV-2 among humans provides opportunities for spillover transmission into nonhuman hosts (5, 7). Consistent with this idea, SARS‐CoV‐2 infections have been documented in dogs, cats, zoo animals, and farmed mink (8, 9). In principle, SARS-CoV-2 infection of a nonhuman animal host might result in establishing a reservoir that can further drive the emergence of novel variants with the potential for spillback to humans. This type of transmission cycle has been described among workers on mink farms (10). However, widespread SARS-CoV-2 transmission in a free-living animal species has not yet been documented.Our study was prompted by a recent report that 40% of free-living white-tailed deer in the United States had antibodies against SARS-CoV-2 (11). Recent studies have also provided evidence of SARS-CoV-2 transmission among experimentally infected deer (6, 12). To test the hypothesis that infection and subsequent transmission of SARS-CoV-2 in deer occurs in nature, we assayed 283 retropharyngeal lymph node (RPLN) samples collected from free-living and captive deer in Iowa from April 2020 through January 2021. Sequence analysis of SARS-CoV-2 genomes in all 94 positive samples discovered multiple lineages corresponding to viral genotypes circulating contemporaneously in humans. In the aggregate, our results are consistent with a model of multiple independent human-to-deer transmission events followed by subsequent deer-to-deer transmission. The findings raise the concerning possibility of reverse zoonoses, especially in exurban areas with high deer density. Our results also highlight the potential risks and considerable knowledge gaps associated with the continued molecular evolution of SARS-CoV-2 in animal hosts. Implementation of enhanced surveillance programs to define the magnitude of the deer reservoir problem and determine whether other reservoir species exist at the animal–human interface is warranted.     

Emergence of Two Distinct SARS-CoV-2 Gamma Variants and the Rapid Spread of P.1-like-II SARS-CoV-2 during the Second Wave of COVID-19 in Santa Catarina,Southern Brazil

The western mesoregion of the state of Santa Catarina (SC), Southern Brazil, was heavily affected as a whole by the COVID-19 pandemic in early 2021. This study aimed to evaluate the dynamics of the SARS-CoV-2 virus spreading patterns in the SC state from March 2020 to April 2021 using genomic surveillance. During this period, there were 23 distinct variants, including Beta and Gamma, among which the Gamma and related lineages were predominant in the second pandemic wave within SC. A regionalization of P.1-like-II in the Western SC region was observed, concomitant to the increase in cases, mortality, and the case fatality rate (CFR) index. This is the first evidence of the regionalization of the SARS-CoV-2 transmission in SC and it highlights the importance of tracking the variants, dispersion, and impact of SARS-CoV-2 on the public health systems.     

Contrasting Epidemiology and Population Genetics of COVID-19 Infections Defined by Multilocus Genotypes in SARS-CoV-2 Genomes Sampled Globally

Since its emergence in 2019, SARS-CoV-2 has spread and evolved globally, with newly emerged variants of concern (VOCs) accounting for more than 500 million COVID-19 cases and 6 million deaths. Continuous surveillance utilizing simple genetic tools is needed to measure the viral epidemiological diversity, risk of infection, and distribution among different demographics in different geographical regions. To help address this need, we developed a proof-of-concept multilocus genotyping tool and demonstrated its utility to monitor viral populations sampled in 2020 and 2021 across six continents. We sampled globally 22,164 SARS-CoV-2 genomes from GISAID (inclusion criteria: available clinical and demographic data). They comprised two study populations, “2020 genomes” (N = 5959) sampled from December 2019 to September 2020 and “2021 genomes” (N = 16,205) sampled from 15 January to 15 March 2021. All genomes were aligned to the SARS-CoV-2 reference genome and amino acid polymorphisms were called with quality filtering. Thereafter, 74 codons (loci) in 14 genes including orf1ab polygene (N = 9), orf3a, orf8, nucleocapsid (N), matrix (M), and spike (S) met the 0.01 minimum allele frequency criteria and were selected to construct multilocus genotypes (MLGs) for the genomes. At these loci, 137 mutant/variant amino acids (alleles) were detected with eight VOC-defining variant alleles, including N KR_203&204, orf1ab (I₂₆₅, F₃₆₀₆, and L₄₇₁₅), orf3a H₅₇, orf8 S₈₄, and S G₆₁₄, being predominant globally with > 35% prevalence. Their persistence and selection were associated with peaks in the viral transmission and COVID-19 incidence between 2020 and 2021. Epidemiologically, older patients (≥20 years) compared to younger patients (<20 years) had a higher risk of being infected with these variants, but this association was dependent on the continent of origin. In the global population, the discriminant analysis of principal components (DAPC) showed contrasting patterns of genetic clustering with three (Africa, Asia, and North America) and two (North and South America) continental clusters being observed for the 2020 and 2021 global populations, respectively. Within each continent, the MLG repertoires (range 40–199) sampled in 2020 and 2021 were genetically differentiated, with ≤4 MLGs per repertoire accounting for the majority of genomes sampled. These data suggested that the majority of SARS-CoV-2 infections in 2020 and 2021 were caused by genetically distinct variants that likely adapted to local populations. Indeed, four GISAID clade-defined VOCs - GRY (Alpha), GH (Beta), GR (Gamma), and G/GK (Delta variant) were differentiated by their MLG signatures, demonstrating the versatility of the MLG tool for variant identification. Results from this proof-of-concept multilocus genotyping demonstrates its utility for SARS-CoV-2 genomic surveillance and for monitoring its spatiotemporal epidemiology and evolution, particularly in response to control interventions including COVID-19 vaccines and chemotherapies.     

The Impact on Infectivity and Neutralization Efficiency of SARS-CoV-2 Lineage B.1.351 Pseudovirus

A new variant of SARS-CoV-2 B.1.351 lineage (first found in South Africa) has been raising global concern due to its harboring of multiple mutations in the spike that potentially increase transmissibility and yield resistance to neutralizing antibodies. We here tested infectivity and neutralization efficiency of SARS-CoV-2 spike pseudoviruses bearing particular mutations of the receptor-binding domain (RBD) derived either from the Wuhan strains (referred to as D614G or with other sites) or the B.1.351 lineage (referred to as N501Y, K417N, and E484K). The three different pseudoviruses B.1.351 lineage related significantly increased infectivity compared with other mutants that indicated Wuhan strains. Interestingly, K417N and E484K mutations dramatically enhanced cell–cell fusion than N501Y even though their infectivity were similar, suggesting that K417N and E484K mutations harboring SARS-CoV-2 variant might be more transmissible than N501Y mutation containing SARS-CoV-2 variant. We also investigated the efficacy of two different monoclonal antibodies, Casirivimab and Imdevimab that neutralized SARS-CoV-2, against several kinds of pseudoviruses which indicated Wuhan or B.1.351 lineage. Remarkably, Imdevimab effectively neutralized B.1.351 lineage pseudoviruses containing N501Y, K417N, and E484K mutations, while Casirivimab partially affected them. Overall, our results underscore the importance of B.1.351 lineage SARS-CoV-2 in the viral spread and its implication for antibody efficacy.