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.     

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.     

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.     

SARS-CoV-2 Variants of Concern Infect the Respiratory Tract and Induce Inflammatory Response in Wild-Type Laboratory Mice

The emergence of new severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) variants of concern pose a major threat to public health, due to possible enhanced virulence, transmissibility and immune escape. These variants may also adapt to new hosts, in part through mutations in the spike protein. In this study, we evaluated the infectivity and pathogenicity of SARS-CoV-2 variants of concern in wild-type C57BL/6 mice. Six-week-old mice were inoculated intranasally with a representative virus from the original B.1 lineage, or the emerging B.1.1.7 and B.1.351 lineages. We also infected a group of mice with a mouse-adapted SARS-CoV-2 (MA10). Viral load and mRNA levels of multiple cytokines and chemokines were analyzed in the lung tissues on day 3 after infection. Our data show that unlike the B.1 virus, the B.1.1.7 and B.1.351 viruses are capable of infecting C57BL/6 mice and replicating at high concentrations in the lungs. The B.1.351 virus replicated to higher titers in the lungs compared with the B.1.1.7 and MA10 viruses. The levels of cytokines (IL-6, TNF-α, IL-1β) and chemokine (CCL2) were upregulated in response to the B.1.1.7 and B.1.351 infection in the lungs. In addition, robust expression of viral nucleocapsid protein and histopathological changes were detected in the lungs of B.1.351-infected mice. Overall, these data indicate a greater potential for infectivity and adaptation to new hosts by emerging SARS-CoV-2 variants.     

Clusters of SARS-CoV-2 Lineage B.1.1.7 Infection after Vaccination with Adenovirus-Vectored and Inactivated Vaccines

A SARS-CoV-2 B.1.1.7 variant of concern (VOC) has been associated with increased transmissibility, hospitalization, and mortality. This study aimed to explore the factors associated with B.1.1.7 VOC infection in the context of vaccination. On March 2021, we detected SARS-CoV-2 RNA in nasopharyngeal samples from 14 of 22 individuals vaccinated with a single-dose of ChAdOx1 (outbreak A, n = 26), and 22 of 42 of individuals with two doses of the CoronaVac vaccine (outbreak B, n = 52) for breakthrough infection rates for ChAdOx1 of 63.6% and 52.4% for CoronaVac. The outbreaks were caused by two independent clusters of the B.1.1.7 VOC. The serum of PCR-positive symptomatic SARS-CoV-2-infected individuals had ~1.8–3.4-fold more neutralizing capacity against B.1.1.7 compared to the serum of asymptomatic individuals. These data based on exploratory analysis suggest that the B.1.1.7 variant can infect individuals partially immunized with a single dose of an adenovirus-vectored vaccine or fully immunized with two doses of an inactivated vaccine, although the vaccines were able to reduce the risk of severe disease and death caused by this VOC, even in the elderly.     

Increased transmissibility and global spread of SARS-CoV-2 variants of concern as at June 2021

We present a global analysis of the spread of recently emerged SARS-CoV-2 variants and estimate changes in effective reproduction numbers at country-specific level using sequence data from GISAID. Nearly all investigated countries demonstrated rapid replacement of previously circulating lineages by the World Health Organization-designated variants of concern, with estimated transmissibility increases of 29% (95% CI: 24–33), 25% (95% CI: 20–30), 38% (95% CI: 29–48) and 97% (95% CI: 76–117), respectively, for B.1.1.7, B.1.351, P.1 and B.1.617.2.     

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.     

Analyses of Long-Term Epidemic Trends and Evolution Characteristics of Haplotype Subtypes Reveal the Dynamic Selection on SARS-CoV-2

The scale of SARS-CoV-2 infection and death is so enormous that further study of the molecular and evolutionary characteristics of SARS-CoV-2 will help us better understand and respond to SARS-CoV-2 outbreaks. The present study analyzed the epidemic and evolutionary characteristics of haplotype subtypes or regions based on 1.8 million high-quality SARS-CoV-2 genomic data. The estimated ratio of the rates of non-synonymous to synonymous changes (Ka/Ks) in North America and the United States were always more than 1.0, while the Ka/Ks in other continents and countries showed a sharp decline, then a slow increase to 1.0, and a dramatic increase over time. H1 (B.1) with the highest substitution rate has become the most dominant haplotype subtype since March 2020 and has evolved into multiple haplotype subtypes with smaller substitution rates. Many evolutionary characteristics of early SARS-CoV-2, such as H3 being the only early haplotype subtype that existed for the shortest time, the global prevalence of H1 and H1-5 (B.1.1) within a month after being detected, and many high divergent genome sequences early in February 2020, indicate the missing of early SARS-CoV-2 genomic data. SARS-CoV-2 experienced dynamic selection from December 2019 to August 2021 and has been under strong positive selection since May 2021. Its transmissibility and the ability of immune escape may be greatly enhanced over time. This will bring greater challenges to the control of 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.     

Emergence and Spread of a B.1.1.28-Derived P.6 Lineage with Q675H and Q677H Spike Mutations in Uruguay

Uruguay controlled the viral dissemination during the first nine months of the SARS-CoV-2 pandemic. Unfortunately, towards the end of 2020, the number of daily new cases exponentially increased. Herein, we analyzed the country-wide genetic diversity of SARS-CoV-2 between November 2020 and April 2021. We identified that the most prevalent viral variant during the first epidemic wave in Uruguay (December 2020–February 2021) was a B.1.1.28 sublineage carrying Spike mutations Q675H + Q677H, now designated as P.6, followed by lineages P.2 and P.7. P.6 probably arose around November 2020, in Montevideo, Uruguay’s capital department, and rapidly spread to other departments, with evidence of further local transmission clusters; it also spread sporadically to the USA and Spain. The more efficient dissemination of lineage P.6 with respect to P.2 and P.7 and the presence of mutations (Q675H and Q677H) in the proximity of the key cleavage site at the S1/S2 boundary suggest that P.6 may be more transmissible than other lineages co-circulating in Uruguay. Although P.6 was replaced by the variant of concern (VOC) P.1 as the predominant lineage in Uruguay since April 2021, the monitoring of the concurrent emergence of Q675H + Q677H in VOCs should be of worldwide interest.     

Multiple Early Introductions of SARS-CoV-2 to Cape Town,South Africa

Cape Town was the first city in South Africa to experience the full impact of the coronavirus disease 2019 (COVID-19) pandemic. We acquired samples from all suspected cases and their contacts during the first month of the pandemic from Tygerberg Hospital. Nanopore sequencing generated SARS-CoV-2 whole genomes. Phylogenetic inference with maximum likelihood and Bayesian methods were used to determine lineages that seeded the local epidemic. Three patients were known to have travelled internationally and an outbreak was detected in a nearby supermarket. Sequencing of 50 samples produced 46 high-quality genomes. The sequences were classified as lineages: B, B.1, B.1.1.1, B.1.1.161, B.1.1.29, B.1.8, B.39, and B.40. All the sequences from persons under investigation (PUIs) in the supermarket outbreak (lineage B.1.8) fall within a clade from the Netherlands with good support (p > 0.9). In addition, a new mutation, 5209A>G, emerged within the Cape Town cluster. The molecular clock analysis suggests that this occurred around 13 March 2020 (95% confidence interval: 9–17 March). The phylogenetic reconstruction suggests at least nine early introductions of SARS-CoV-2 into Cape Town and an early localized transmission in a shopping environment. Genomic surveillance was successfully used to investigate and track the spread of early introductions of SARS-CoV-2 in Cape Town.     

Genomic monitoring of SARS-CoV-2 variants using sentinel SARI hospital surveillance

Background

To support the COVID-19 pandemic response, many countries, including Belgium, implemented baseline genomic surveillance (BGS) programs aiming to early detect and characterize new SARS-CoV-2 variants. In parallel, Belgium maintained a sentinel network of six hospitals that samples patients with severe acute respiratory infections (SARI) and integrated SARS-CoV-2 detection within a broader range of respiratory pathogens. We evaluate the ability of the SARI surveillance to monitor general trends and early signals of viral genetic evolution of SARS-CoV-2 and compare it with the BGS as a reference model.

Methods

Nine-hundred twenty-five SARS-CoV-2 positive samples from patients fulfilling the Belgian SARI definition between January 2020 and December 2022 were sequenced using the ARTIC Network amplicon tiling approach on a MinION platform. Weekly variant of concern (VOC) proportions and types were compared to those that were circulating between 2021 and 2022, using 96,251 sequences of the BGS.

Results

SARI surveillance allowed timely detection of the Omicron (BA.1, BA.2, BA.4, and BA.5) and Delta (B.1.617.2) VOCs, with no to 2 weeks delay according to the start of their epidemic growth in the Belgian population. First detection of VOCs B.1.351 and P.1 took longer, but these remained minor in Belgium. Omicron BA.3 was never detected in SARI surveillance. Timeliness could not be evaluated for B.1.1.7, being already major at the start of the study period.

Conclusions

Genomic surveillance of SARS-CoV-2 using SARI sentinel surveillance has proven to accurately reflect VOCs detected in the population and provides a cost-effective solution for long-term genomic monitoring of circulating respiratory viruses.     

Rise and Fall of SARS-CoV-2 Lineage A.27 in Germany

Here, we report on the increasing frequency of the SARS-CoV-2 lineage A.27 in Germany during the first months of 2021. Genomic surveillance identified 710 A.27 genomes in Germany as of 2 May 2021, with a vast majority identified in laboratories from a single German state (Baden-Wuerttemberg, n = 572; 80.5%). Baden-Wuerttemberg is located near the border with France, from where most A.27 sequences were entered into public databases until May 2021. The first appearance of this lineage based on sequencing in a laboratory in Baden-Wuerttemberg can be dated to early January ’21. From then on, the relative abundance of A.27 increased until the end of February but has since declined—meanwhile, the abundance of B.1.1.7 increased in the region. The A.27 lineage shows a mutational pattern typical of VOIs/VOCs, including an accumulation of amino acid substitutions in the Spike glycoprotein. Among those, L18F, L452R and N501Y are located in the epitope regions of the N-terminal- (NTD) or receptor binding domain (RBD) and have been suggested to result in immune escape and higher transmissibility. In addition, A.27 does not show the D614G mutation typical for all VOIs/VOCs from the B lineage. Overall, A.27 should continue to be monitored nationally and internationally, even though the observed trend in Germany was initially displaced by B.1.1.7 (Alpha), while now B.1.617.2 (Delta) is on the rise.     

Serum Neutralizing Activity against B.1.1.7, B.1.351, and P.1 SARS-CoV-2 Variants of Concern in Hospitalized COVID-19 Patients

The recent spreading of new SARS-CoV-2 variants, carrying several mutations in the spike protein, could impact immune protection elicited by natural infection or conferred by vaccination. In this study, we evaluated the neutralizing activity against the viral variants that emerged in the United Kingdom (B.1.1.7), Brazil (P.1), and South Africa (B.1.351) in human serum samples from hospitalized patients infected by SARS-CoV-2 during the first pandemic wave in Italy in 2020. Of the patients studied, 59.5% showed a decrease (≥2 fold) in neutralizing antibody titer against B.1.1.7, 83.3% against P.1, and 90.5% against B.1.351 with respect to the original strain. The reduction in antibody titers against all analyzed variants, and in particular P.1 and B.1.351, suggests that previous symptomatic infection might be not fully protective against exposure to SARS-CoV-2 variants carrying a set of relevant spike mutations.     

First Detection of SARS-CoV-2 B.1.1.7 Variant of Concern in an Asymptomatic Dog in Spain

Natural SARS-CoV-2 infection in pets has been widely documented during the last year. Although the majority of reports suggested that dogs’ susceptibility to the infection is low, little is known about viral pathogenicity and transmissibility in the case of variants of concern, such as B.1.1.7 in this species. Here, as part of a large-scale study on SARS-CoV-2 prevalence in pets in Spain, we have detected the B.1.1.7 variant of concern (VOC) in a dog whose owners were infected with SARS-CoV-2. The animal did not present any symptoms, but viral loads were high in the nasal and rectal swabs. In addition, viral isolation was possible from both swabs, demonstrating that the dog was shedding infectious virus. Seroconversion occurred 23 days after the first sampling. This study documents the first detection of B.1.1.7 VOC in a dog in Spain and emphasizes the importance of performing active surveillance and genomic investigation on infected animals.     

Replication Kinetics of B.1.351 and B.1.1.7 SARS-CoV-2 Variants of Concern Including Assessment of a B.1.1.7 Mutant Carrying a Defective ORF7a Gene

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the etiological agent of COVID-19, is a readily transmissible and potentially deadly pathogen which is currently re-defining human susceptibility to pandemic viruses in the modern world. The recent emergence of several genetically distinct descendants known as variants of concern (VOCs) is further challenging public health disease management, due to increased rates of virus transmission and potential constraints on vaccine effectiveness. We report the isolation of SARS-CoV-2 VOCs imported into Australia belonging to the B.1.351 lineage, first described in the Republic of South Africa (RSA), and the B.1.1.7 lineage originally reported in the United Kingdom, and directly compare the replication kinetics of these two VOCs in Vero E6 cells. In this analysis, we also investigated a B.1.1.7 VOC (QLD1516/2021) carrying a 7-nucleotide deletion in the open reading frame 7a (ORF7a) gene, likely truncating and rendering the ORF7a protein of this virus defective. We demonstrate that the replication of the B.1.351 VOC (QLD1520/2020) in Vero E6 cells can be detected earlier than the B.1.1.7 VOCs (QLD1516/2021 and QLD1517/2021), before peaking at 48 h post infection (p.i.), with significantly higher levels of virus progeny. Whilst replication of the ORF7a defective isolate QLD1516/2021 was delayed longer than the other viruses, slightly more viral progeny was produced by the mutant compared to the unmutated isolate QLD1517/2021 at 72 h p.i. Collectively, these findings contribute to our understanding of SARS-CoV-2 replication and evolutionary dynamics, which have important implications in the development of future vaccination, antiviral therapies, and epidemiological control strategies for COVID-19.     

National Scale Real-Time Surveillance of SARS-CoV-2 Variants Dynamics by Wastewater Monitoring in Israel

In this report, we describe a national-scale monitoring of the SARS-CoV-2 (SC-2) variant dynamics in Israel, using multiple-time sampling of 13 wastewater treatment plants. We used a combination of inclusive and selective quantitative PCR assays that specifically identify variants A19/A20 or B.1.1.7 and tested each sample for the presence and relative viral RNA load of each variant. We show that between December 2020 and March 2021, a complete shift in the SC-2 variant circulation was observed, where the B.1.1.7 replaced the A19 in all examined test points. We further show that the normalized viral load (NVL) values and the average new cases per week reached a peak in January 2021 and then decreased gradually in almost all test points, in parallel with the progression of the national vaccination campaign, during February–March 2021. This study demonstrates the importance of monitoring SC-2 variant by using a combination of inclusive and selective PCR tests on a national scale through wastewater sampling, which is far more amendable for high-throughput monitoring compared with sequencing. This approach may be useful for real-time dynamics surveillance of current and future variants, such as the Omicron (BA.1, BA.2) and other variants.     

Impact of January 2021 curfew measures on SARS-CoV-2 B.1.1.7 circulation in France

Following the spread of the SARS-CoV-2 B.1.1.7 variant, social distancing was strengthened in France in January 2021. Using a two-strain mathematical model calibrated on genomic surveillance, we estimated that curfew measures allowed hospitalisations to plateau by decreasing transmission of the historical strains while B.1.1.7 continued to grow. School holidays appear to have further slowed down progression in February. Without progressively strengthened social distancing, a rapid surge of hospitalisations is expected, despite the foreseen increase in vaccination rhythm.     

The Alpha Variant (B.1.1.7) of SARS-CoV-2 in Children: First Experience from 3544 Nucleic Acid Amplification Tests in a Cohort of Children in Germany

In May 2021, the Alpha variant (B.1.1.7) of SARS-CoV-2 was found in 91% of the SARS-CoV-2 cases in Germany. Not much is known about the symptoms, courses of disease, and infectiousness in pediatric patients with the Alpha variant. Objective: The aim of this retrospective analysis was to gain information on the infection with the Alpha variant in children and adolescents. Methods: Between 12 January 2021 and 3 June 2021, all nucleic acid amplification tests (NAATs) of children who received a swab for SARS-CoV-2 were included. Data were collected on standardized questionnaires. The analysis of data was anonymized and retrospective. Results: We investigated 3544 NAATs; 95 children were tested positive (2.7%) for SARS-CoV-2. For the sub-analysis, 65 children were analyzed. In 59 children, the Alpha variant was found (90.8%), and 54.2% (n = 32/59) were symptomatic. The most common symptoms were fever, cough, and rhinitis. The median Ct value was 24.0 (min 17.0; max 32.7). Conclusions: We can underline early findings that children are still less effected by SARS-CoV-2 infection with the spread of the Alpha variant. We found no evidence that children infected with the Alpha variant showed more severe symptoms or suffered from a more severe clinical course than those infected with the wild type.     

First confirmation of importation and transmission in Spain of the newly identified SARS-CoV-2 B.1.1.7 variant

IntroductionA newly identified SARS-CoV-2 variant, VOC202012/01 originating lineage B.1.1.7, recently emerged in the United Kingdom. The rapid spread in the UK of this new variant has caused other countries to be vigilant.Material and methodsWe based our initial screening of B.1.1.7 on the dropout of the S gene signal in the TaqPath assay, caused by the 69/70 deletion. Subsequently, we confirmed the B.1.1.7 candidates by whole genome sequencing.ResultsWe describe the first three imported cases of this variant from London to Madrid, subsequent post-arrival household transmission to three relatives, and the two first cases without epidemiological links to UK. One case required hospitalization. In all cases, drop-out of gene S was correctly associated to the B.1.1.7 variant, as all the corresponding sequences carried the 17 lineage-marker mutations.ConclusionThe first identifications of the SARS-CoV-2 B.1.1.7 variant in Spain indicate the role of independent introductions from the UK coexisting with post-arrival transmission in the community, since the early steps of this new variant in our country.