CoV-Spectrum: analysis of globally shared SARS-CoV-2 data to identify and characterize new variants

TRIM7 ubiquitinates SARS-CoV-2 membrane protein to limit apoptosis and viral replication

SARS-CoV-2 is a highly transmissible virus that causes COVID-19 disease. Mechanisms of viral pathogenesis include excessive inflammation and viral-induced cell death, resulting in tissue damage. We identified the host E3-ubiquitin ligase TRIM7 as an inhibitor of apoptosis and SARS-CoV-2 replication via ubiquitination of the viral membrane (M) protein. Trim7 -/- mice exhibited increased pathology and virus titers associated with epithelial apoptosis and dysregulated immune responses. Mechanistically, TRIM7 ubiquitinates M on K14, which protects cells from cell death. Longitudinal SARS-CoV-2 sequence analysis from infected patients revealed that mutations on M-K14 appeared in circulating variants during the pandemic. The relevance of these mutations was tested in a mouse model. A recombinant M-K14/K15R virus showed reduced viral replication, consistent with the role of K15 in virus assembly, and increased levels of apoptosis associated with the loss of ubiquitination on K14. TRIM7 antiviral activity requires caspase-6 inhibition, linking apoptosis with viral replication and pathology.

A bispecific antibody exhibits broad neutralization against SARS-CoV-2 Omicron variants XBB.1.16, BQ.1.1 and sarbecoviruses

The Omicron subvariants BQ.1.1, XBB.1.5, and XBB.1.16 of SARS-CoV-2 are known for their adeptness at evading immune responses. Here, we isolate a neutralizing antibody, 7F3, with the capacity to neutralize all tested SARS-CoV-2 variants, including BQ.1.1, XBB.1.5, and XBB.1.16. 7F3 targets the receptor-binding motif (RBM) region and exhibits broad binding to a panel of 37 RBD mutant proteins. We develop the IgG-like bispecific antibody G7-Fc using 7F3 and the cross-neutralizing antibody GW01. G7-Fc demonstrates robust neutralizing activity against all 28 tested SARS-CoV-2 variants and sarbecoviruses, providing potent prophylaxis and therapeutic efficacy against XBB.1 infection in both K18-ACE and BALB/c female mice. Cryo-EM structure analysis of the G7-Fc in complex with the Omicron XBB spike (S) trimer reveals a trimer-dimer conformation, with G7-Fc synergistically targeting two distinct RBD epitopes and blocking ACE2 binding. Comparative analysis of 7F3 and LY-CoV1404 epitopes highlights a distinct and highly conserved epitope in the RBM region bound by 7F3, facilitating neutralization of the immune-evasive Omicron variant XBB.1.16. G7-Fc holds promise as a potential prophylactic countermeasure against SARS-CoV-2, particularly against circulating and emerging variants.

Replication kinetics and infectivity of SARS-CoV-2 Omicron variant sublineages recovered in the Republic of Korea

BACKGROUND

We analyzed the correlation between the infectivity and transmissibility of the severe acute respiratory syndrome coronavirus 2 Omicron sublineages BA.1, BA. 2, BA.4, and BA.5.

METHODS

We assessed viral replication kinetics and infectivity at the cellular level. Nasopharyngeal and oropharyngeal specimens were obtained from patients with coronavirus disease 2019, confirmed using whole-genome sequencing to be caused by the Omicron sublineages BA.1, BA.2, BA.4, or BA.5. These specimens were used to infect Vero E6 cells, derived from monkey kidneys, for the purpose of viral isolation. Viral stocks were then passaged in Vero E6 cells at a multiplicity of infection of 0.01, and culture supernatants were harvested at 12-hour intervals for 72 hours. To evaluate viral replication kinetics, we determined the cycle threshold values of the supernatants using real-time reverse transcription polymerase chain reaction and converted these values to genome copy numbers.

RESULTS

The viral load was comparable between BA.2, BA.4, and BA.5, whereas BA.1 exhibited a lower value. The peak infectious load of BA.4 was approximately 3 times lower than that of BA.2 and BA.5, while the peak load of BA.2 and BA.5 was about 7 times higher than that of BA.1. Notably, BA.1 demonstrated the lowest infectivity over the entire study period.

CONCLUSION

Our results suggest that the global BA.5 wave may have been amplified by the higher viral replication and infectivity of BA.5 compared to other Omicron sublineages. These analyses could support the rapid assessment of the impact of novel variants on case incidence.

Social contacts in Switzerland during the COVID-19 pandemic: Insights from the CoMix study

To mitigate the spread of SARS-CoV-2, the Swiss government enacted restrictions on social contacts from 2020 to 2022. In addition, individuals changed their social contact behavior to limit the risk of COVID-19. In this study, we aimed to investigate the changes in social contact patterns of the Swiss population. As part of the CoMix study, we conducted a survey consisting of 24 survey waves from January 2021 to May 2022. We collected data on social contacts and constructed contact matrices for the age groups 0-4, 5-14, 15-29, 30-64, and 65 years and older. We estimated the change in contact numbers during the COVID-19 pandemic to a synthetic pre-pandemic contact matrix. We also investigated the association of the largest eigenvalue of the social contact and transmission matrices with the stringency of pandemic measures, the effective reproduction number (Re), and vaccination uptake. During the pandemic period, 7084 responders reported an average number of 4.5 contacts (95% confidence interval, CI: 4.5-4.6) per day overall, which varied by age and survey wave. Children aged 5-14 years had the highest number of contacts with 8.5 (95% CI: 8.1-8.9) contacts on average per day and participants that were 65 years and older reported the fewest (3.4, 95% CI: 3.2-3.5) per day. Compared with the pre-pandemic baseline, we found that the 15-29 and 30-64 year olds had the largest reduction in contacts. We did not find statistically significant associations between the largest eigenvalue of the social contact and transmission matrices and the stringency of measures, Re, or vaccination uptake. The number of social contacts in Switzerland fell during the COVID-19 pandemic and remained below pre-pandemic levels after contact restrictions were lifted. The collected social contact data will be critical in informing modeling studies on the transmission of respiratory infections in Switzerland and to guide pandemic preparedness efforts.

Differences and similarities between innate immune evasion strategies of human coronaviruses

So far, seven coronaviruses have emerged in humans. Four recurring endemic coronaviruses cause mild respiratory symptoms. Infections with epidemic Middle East respiratory syndrome-related coronavirus or severe acute respiratory syndrome coronavirus (SARS-CoV)-1 are associated with high mortality rates. SARS-CoV-2 is the causative agent of the coronavirus disease 2019 pandemic. To establish an infection, coronaviruses evade restriction by human innate immune defenses, such as the interferon system, autophagy and the inflammasome. Here, we review similar and distinct innate immune manipulation strategies employed by the seven human coronaviruses. We further discuss the impact on pathogenesis, zoonotic emergence and adaptation. Understanding the nature of the interplay between endemic/epidemic/pandemic coronaviruses and host defenses may help to better assess the pandemic potential of emerging coronaviruses.

Explaining the impact of mutations on quantification of SARS-CoV-2 in wastewater

Wastewater surveillance is an effective tool for monitoring community spread of COVID-19 and other diseases. Quantitative PCR (qPCR) analysis for wastewater surveillance is more susceptible to mutations in target genome regions than binary PCR analysis for clinical surveillance. The SARS-CoV-2 concentrations in wastewater estimated by N1 and N2 qPCR assays started to diverge around July 2022 in data from different sampling sites, analytical methods, and analytical laboratories in Japan. On the basis of clinical genomic surveillance data and experimental data, we demonstrate that the divergence is due to two mutations in the N1 probe region, which can cause underestimation of viral concentrations. We further show that this inaccuracy can be alleviated if the qPCR data are analyzed with the second derivative method or the Cy0 method instead of the crossing point method.

Investigating the cause of a 2021 winter wave of COVID-19 in a border region in eastern Germany: a mixed-methods study, August to November 2021

It is so far unclear how the COVID-19 winter waves started and what should be done to prevent possible future waves. In this study, we deciphered the dynamic course of a winter wave in 2021 in Saxony, a state in Eastern Germany neighbouring the Czech Republic and Poland. The study was carried out through the integration of multiple virus genomic epidemiology approaches to track transmission chains, identify emerging variants and investigate dynamic changes in transmission clusters. For identified local variants of interest, functional evaluations were performed. Multiple long-lasting community transmission clusters have been identified acting as driving force for the winter wave 2021. Analysis of the dynamic courses of two representative clusters indicated a similar transmission pattern. However, the transmission cluster caused by a locally occurring new Delta variant AY.36.1 showed a distinct transmission pattern, and functional analyses revealed a replication advantage of it. This study indicated that long-lasting community transmission clusters starting since early autumn caused by imported or locally occurring variants all contributed to the development of the 2021 winter wave. The information we achieved might help future pandemic prevention.

The Canadian VirusSeq Data Portal & Duotang: open resources for SARS-CoV-2 viral sequences and genomic epidemiology

The COVID-19 pandemic led to a large global effort to sequence SARS-CoV-2 genomes from patient samples to track viral evolution and inform public health response. Millions of SARS-CoV-2 genome sequences have been deposited in global public repositories. The Canadian COVID-19 Genomics Network (CanCOGeN - VirusSeq), a consortium tasked with coordinating expanded sequencing of SARS-CoV-2 genomes across Canada early in the pandemic, created the Canadian VirusSeq Data Portal, with associated data pipelines and procedures, to support these efforts. The goal of VirusSeq was to allow open access to Canadian SARS-CoV-2 genomic sequences and enhanced, standardized contextual data that were unavailable in other repositories and that meet FAIR standards (Findable, Accessible, Interoperable and Reusable). In addition, the Portal data submission pipeline contains data quality checking procedures and appropriate acknowledgement of data generators that encourages collaboration. From inception to execution, the portal was developed with a conscientious focus on strong data governance principles and practices. Extensive efforts ensured a commitment to Canadian privacy laws, data security standards, and organizational processes. This Portal has been coupled with other resources like Viral AI and was further leveraged by the Coronavirus Variants Rapid Response Network (CoVaRR-Net) to produce a suite of continually updated analytical tools and notebooks. Here we highlight this Portal, including its contextual data not available elsewhere, and the 'Duotang', a web platform that presents key genomic epidemiology and modeling analyses on circulating and emerging SARS-CoV-2 variants in Canada. Duotang presents dynamic changes in variant composition of SARS-CoV-2 in Canada and by province, estimates variant growth, and displays complementary interactive visualizations, with a text overview of the current situation. The VirusSeq Data Portal and Duotang resources, alongside additional analyses and resources computed from the Portal (COVID-MVP, CoVizu), are all open-source and freely available. Together, they provide an updated picture of SARS-CoV-2 evolution to spur scientific discussions, inform public discourse, and support communication with and within public health authorities. They also serve as a framework for other jurisdictions interested in open, collaborative sequence data sharing and analyses.

Tracking SARS-CoV-2 variants of concern in wastewater: an assessment of nine computational tools using simulated genomic data

Wastewater-based surveillance (WBS) is an important epidemiological and public health tool for tracking pathogens across the scale of a building, neighbourhood, city, or region. WBS gained widespread adoption globally during the SARS-CoV-2 pandemic for estimating community infection levels by qPCR. Sequencing pathogen genes or genomes from wastewater adds information about pathogen genetic diversity, which can be used to identify viral lineages (including variants of concern) that are circulating in a local population. Capturing the genetic diversity by WBS sequencing is not trivial, as wastewater samples often contain a diverse mixture of viral lineages with real mutations and sequencing errors, which must be deconvoluted computationally from short sequencing reads. In this study we assess nine different computational tools that have recently been developed to address this challenge. We simulated 100 wastewater sequence samples consisting of SARS-CoV-2 BA.1, BA.2, and Delta lineages, in various mixtures, as well as a Delta-Omicron recombinant and a synthetic 'novel' lineage. Most tools performed well in identifying the true lineages present and estimating their relative abundances and were generally robust to variation in sequencing depth and read length. While many tools identified lineages present down to 1 % frequency, results were more reliable above a 5 % threshold. The presence of an unknown synthetic lineage, which represents an unclassified SARS-CoV-2 lineage, increases the error in relative abundance estimates of other lineages, but the magnitude of this effect was small for most tools. The tools also varied in how they labelled novel synthetic lineages and recombinants. While our simulated dataset represents just one of many possible use cases for these methods, we hope it helps users understand potential sources of error or bias in wastewater sequencing analysis and to appreciate the commonalities and differences across methods.

Performance of amplicon and capture based next-generation sequencing approaches for the epidemiological surveillance of Omicron SARS-CoV-2 and other variants of concern

To control the SARS-CoV-2 pandemic, healthcare systems have focused on ramping up their capacity for epidemiological surveillance through viral whole genome sequencing. In this paper, we tested the performance of two protocols of SARS-CoV-2 nucleic acid enrichment, an amplicon enrichment using different versions of the ARTIC primer panel and a hybrid-capture method using KAPA RNA Hypercap. We focused on the challenge of the Omicron variant sequencing, the advantages of automated library preparation and the influence of the bioinformatic analysis in the final consensus sequence. All 94 samples were sequenced using Illumina iSeq 100 and analysed with two bioinformatic pipelines: a custom-made pipeline and an Illumina-owned pipeline. We were unsuccessful in sequencing six samples using the capture enrichment due to low reads. On the other hand, amplicon dropout and mispriming caused the loss of mutation G21987A and the erroneous addition of mutation T15521A respectively using amplicon enrichment. Overall, we found high sequence agreement regardless of method of enrichment, bioinformatic pipeline or the use of automation for library preparation in eight different SARS-CoV-2 variants. Automation and the use of a simple app for bioinformatic analysis can simplify the genotyping process, making it available for more diagnostic facilities and increasing global vigilance.

INSaFLU-TELEVIR: an open web-based bioinformatics suite for viral metagenomic detection and routine genomic surveillance

BACKGROUND

Implementation of clinical metagenomics and pathogen genomic surveillance can be particularly challenging due to the lack of bioinformatics tools and/or expertise. In order to face this challenge, we have previously developed INSaFLU, a free web-based bioinformatics platform for virus next-generation sequencing data analysis. Here, we considerably expanded its genomic surveillance component and developed a new module (TELEVIR) for metagenomic virus identification.

RESULTS

The routine genomic surveillance component was strengthened with new workflows and functionalities, including (i) a reference-based genome assembly pipeline for Oxford Nanopore technologies (ONT) data; (ii) automated SARS-CoV-2 lineage classification; (iii) Nextclade analysis; (iv) Nextstrain phylogeographic and temporal analysis (SARS-CoV-2, human and avian influenza, monkeypox, respiratory syncytial virus (RSV A/B), as well as a "generic" build for other viruses); and (v) algn2pheno for screening mutations of interest. Both INSaFLU pipelines for reference-based consensus generation (Illumina and ONT) were benchmarked against commonly used command line bioinformatics workflows for SARS-CoV-2, and an INSaFLU snakemake version was released. In parallel, a new module (TELEVIR) for virus detection was developed, after extensive benchmarking of state-of-the-art metagenomics software and following up-to-date recommendations and practices in the field. TELEVIR allows running complex workflows, covering several combinations of steps (e.g., with/without viral enrichment or host depletion), classification software (e.g., Kaiju, Kraken2, Centrifuge, FastViromeExplorer), and databases (RefSeq viral genome, Virosaurus, etc.), while culminating in user- and diagnosis-oriented reports. Finally, to potentiate real-time virus detection during ONT runs, we developed findONTime, a tool aimed at reducing costs and the time between sample reception and diagnosis.

CONCLUSIONS

The accessibility, versatility, and functionality of INSaFLU-TELEVIR are expected to supply public and animal health laboratories and researchers with a user-oriented and pan-viral bioinformatics framework that promotes a strengthened and timely viral metagenomic detection and routine genomics surveillance. INSaFLU-TELEVIR is compatible with Illumina, Ion Torrent, and ONT data and is freely available at https://insaflu.insa.pt/ (online tool) and https://github.com/INSaFLU (code).

Surfing the Waves of SARS-CoV-2: Analysis of Viral Genome Variants Using an NGS Survey in Verona, Italy

The availability of new technologies for deep sequencing, including next-generation sequencing (NGS), allows for the detection of viral genome variations. The epidemiological determination of SARS-CoV-2 viral genome changes during the pandemic waves displayed the genome evolution and subsequent onset of variants over time. These variants were often associated with a different impact on viral transmission and disease severity. We investigated, in a retrospective study, the trend of SARS-CoV-2-positive samples collected from the start of the Italian pandemic (January 2020) to June 2023. In addition, viral RNAs extracted from 938 nasopharyngeal swab samples were analyzed using NGS between February 2022 and June 2023. Sequences were analyzed with bioinformatic tools to identify lineages and mutations and for phylogenetic studies. Six pandemic waves were detected. In our samples, we predominantly detected BA.2, BQ.1, BA.5.1, BA.5.2, and, more recently, XBB.1 and its subvariants. The data describe the SARS-CoV-2 genome evolution involved in viral interactions with the host and the dynamics of specific genome mutations and deletions.

Genomic Surveillance for SARS-CoV-2 Variants: Dominance of XBB Replacement - China, January-June 2023

Introduction

In the first half of 2023, a global shift was observed towards the predominance of XBB variants. China faced a significant epidemic between late 2022 and early 2023 due to Omicron subvariants BA.5.2 and BF.7. This study aims to depict the evolving variant distribution among provincial-level administrative divisions (PLADs) in China and explore the factors driving the predominance of XBB replacement.

Methods

Sequences from local and imported coronavirus disease 2019 (COVID-19) cases recorded between January 1 and June 30, 2023, were included. The study analyzed the changing distribution of viral variants and assessed how the prior dominance of specific variants, XBB subvariants, and imported cases influenced the prevalence of the XBB replacement variant.

Results

A total of 56,486 sequences were obtained from local cases, and 8,669 sequences were from imported cases. Starting in April, there was a shift in the prevalence of XBB from imported to local cases, with varying dominance among PLADs. In PLADs previously high in BF.7, the rise of XBB was delayed. A positive correlation was found between XBB proportions in imported cases from January to March and local cases in April. The distribution pattern of XBB subvariants differed between local and imported cases within the same PLAD. No significant differences were noted in the replacement rates of XBB subvariants.

Conclusions

The timing of XBB dominance differed among various PLADs in China in the first half of 2023, correlating closely with the prevalence of XBB variants among imported cases.

Detection of SARS-COV-2 variants and their proportions in wastewater samples using next-generation sequencing in Finland

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants may have different characteristics, e.g., in transmission, mortality, and the effectiveness of vaccines, indicating the importance of variant detection at the population level. Wastewater-based surveillance of SARS-CoV-2 RNA fragments has been shown to be an effective way to monitor the COVID-19 pandemic at the population level. Wastewater is a complex sample matrix affected by environmental factors and PCR inhibitors, causing insufficient coverage in sequencing, for example. Subsequently, results where part of the genome does not have sufficient coverage are not uncommon. To identify variants and their proportions in wastewater over time, we utilized next-generation sequencing with the ARTIC Network's primer set and bioinformatics pipeline to evaluate the presence of variants in partial genome data. Based on the wastewater data from November 2021 to February 2022, the Delta variant was dominant until mid-December in Helsinki, Finland's capital, and thereafter in late December 2022 Omicron became the most common variant. At the same time, the Omicron variant of SARS-CoV-2 outcompeted the previous Delta variant in Finland in new COVID-19 cases. The SARS-CoV-2 variant findings from wastewater are in agreement with the variant information obtained from the patient samples when visually comparing trends in the sewerage network area. This indicates that the sequencing of wastewater is an effective way to monitor temporal and spatial trends of SARS-CoV-2 variants at the population level.

Tracing the origin of SARS-CoV-2 omicron-like spike sequences detected in an urban sewershed: a targeted, longitudinal surveillance study of a cryptic wastewater lineage

BACKGROUND

The origin of novel SARS-CoV-2 spike sequences found in wastewater, without corresponding detection in clinical specimens, remains unclear. We sought to determine the origin of one such cryptic wastewater lineage by tracking and characterising its persistence and genomic evolution over time.

METHODS

We first detected a cryptic lineage, WI-CL-001, in municipal wastewater in Wisconsin, USA, in January, 2022. To determine the source of WI-CL-001, we systematically sampled wastewater from targeted sub-sewershed lines and maintenance holes using compositing autosamplers. Viral concentrations in wastewater samples over time were measured by RT digital PCR. In addition to using metagenomic 12s rRNA sequencing to determine the virus's host species, we also sequenced SARS-CoV-2 spike receptor binding domains, and, where possible, whole viral genomes to identify and characterise the evolution of this lineage.

FINDINGS

We traced WI-CL-001 to its source at a single commercial building. There we detected the cryptic lineage at concentrations as high as 2·7 × 109 genome copies per L. The majority of 12s rRNA sequences detected in wastewater leaving the identified source building were human. Additionally, we generated over 100 viral receptor binding domain and whole-genome sequences from wastewater samples containing the cryptic lineage collected over the 13 consecutive months this virus was detectable (January, 2022, to January, 2023). These sequences contained a combination of fixed nucleotide substitutions characteristic of Pango lineage B.1.234, which circulated in humans in Wisconsin at low levels from October, 2020, to February, 2021. Despite this, mutations in the spike gene and elsewhere resembled those subsequently found in omicron variants.

INTERPRETATION

We propose that prolonged detection of WI-CL-001 in wastewater indicates persistent shedding of SARS-CoV-2 from a single human initially infected by an ancestral B.1.234 virus. The accumulation of convergent omicron-like mutations in WI-CL-001's ancestral B.1.234 genome probably reflects persistent infection and extensive within-host evolution. People who shed cryptic lineages could be an important source of highly divergent viruses that sporadically emerge and spread.

FUNDING

The Rockefeller Foundation, Wisconsin Department of Health Services, Centers for Disease Control and Prevention, National Institute on Drug Abuse, and the Center for Research on Influenza Pathogenesis and Transmission.

Evolving immune evasion and transmissibility of SARS-CoV-2: The emergence of JN.1 variant and its global impact

The continuous evolution of SARS-CoV-2 variants constitutes a significant impediment to the public health. The World Health Organization (WHO) has designated the SARS-CoV-2 variant JN.1, which has evolved from its progenitor BA.2.86, as a Variant of Interest (VOI) in light of its enhanced immune evasion and transmissibility. The proliferating dissemination of JN.1 globally accentuates its competitive superiority and the potential to instigate fresh surges of infection, notably among cohorts previously infected by antecedent variants. Notably, prevailing evidence does not corroborate an increase in pathogenicity associated with JN.1, and antiviral agents retain their antiviral activity against both BA.2.86 and JN.1. The sustained effectiveness of antiviral agents offers a beacon of hope. Nonetheless, the variant's adeptness at eluding the immunoprotective effects conferred by extant vaccines highlights the imperative for the development of more effective vaccines and therapeutic approaches. Overall, the distinct evolutionary trajectories of BA.2.86 and JN.1 underscore the necessity for ongoing surveillance and scholarly inquiry to elucidate their implications for the pandemic's evolution, which requires the international communities to foster collaboration through the sharing of data, exchange of insights, and collective scientific endeavors.

Results from the SARS-CoV-2 wastewater-based surveillance system in Denmark, July 2021 to June 2022

The microbiological analysis of wastewater samples is increasingly used for the surveillance of SARS-CoV-2 globally. We described the setup process of the national SARS-CoV-2 wastewater-based surveillance system in Denmark, presented its main results during the first year of activities, from July 2021 to June 2022, and discussed their operational significance. The Danish SARS-CoV-2 wastewater-based surveillance system was designed to cover 85 % of the population in Denmark and it entailed taking three weekly samples from 230 sites. Samples were RT-qPCR tested for SARS-CoV-2 RNA, targeting the genetic markers N1, N2 and RdRp, and for two faecal indicators, Pepper Mild Mottle Virus and crAssphage. We calculated the weekly SARS-CoV-2 RNA concentration in the wastewater from each sampling site and monitored it in view of the results from individual testing, at the national and regional levels. We attempted to use wastewater results to identify potential local outbreaks, and we sequenced positive wastewater samples using Nanopore sequencing to monitor the circulation of viral variants in Denmark. The system reached its full implementation by October 2021 and covered up to 86.4 % of the Danish population. The system allowed for monitoring of the national and regional trends of SARS-CoV-2 infections in Denmark. However, the system contribution to the identification of potential local outbreaks was limited by the extensive information available from clinical testing. The sequencing of wastewater samples identified relevant variants of concern, in line with results from sequencing of human samples. Amidst the COVID-19 pandemic, Denmark implemented a nationwide SARS-CoV-2 wastewater-based surveillance system that integrated routine surveillance from individual testing. Today, while testing for COVID-19 at the community level has been discontinued, the system is on the frontline to monitor the occurrence and spread of SARS-CoV-2 in Denmark.

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

The SARS-CoV-2 XBB is a group of highly immune-evasive lineages of the Omicron variant of concern 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.

IMPORTANCE

Brazil was one the most affected countries by the SARS-CoV-2 pandemic, with more than 700,000 deaths by mid-2023. This study reconstructs the dissemination of the virus in the country in the first half of 2023, a period characterized by the dissemination of descendants of XBB.1, a recombinant of Omicron BA.2 lineages evolved in late 2022. The analysis supports that XBB dissemination was marked by the continuous emergence of indigenous lineages bearing similar mutations in key sites of their Spike protein, a process followed by continuous increments in transmissibility, and without repercussions in the incidence of severe cases. Thus, the results suggest that the epidemiological impact of the spread of a SARS-CoV-2 variant is influenced by an intricate interplay of factors that extend beyond the virus's transmissibility alone. The study also underlines the need for SARS-CoV-2 genomic surveillance that allows the monitoring of its ever-shifting composition.

Exploration of low-frequency allelic variants of SARS-CoV-2 genomes reveals coinfections in Mexico occurred during periods of VOCs turnover

A total of 14 973 alleles in 29 661 sequenced samples collected between March 2021 and January 2023 by the Mexican Consortium for Genomic Surveillance (CoViGen-Mex) and collaborators were used to construct a thorough map of mutations of the Mexican SARS-CoV-2 genomic landscape containing Intra-Patient Minor Allelic Variants (IPMAVs), which are low-frequency alleles not ordinarily present in a genomic consensus sequence. This additional information proved critical in identifying putative coinfecting variants included alongside the most common variants, B.1.1.222, B.1.1.519, and variants of concern (VOCs) Alpha, Gamma, Delta, and Omicron. A total of 379 coinfection events were recorded in the dataset (a rate of 1.28 %), resulting in the first such catalogue in Mexico. The most common putative coinfections occurred during the spread of Delta or after the introduction of Omicron BA.2 and its descendants. Coinfections occurred constantly during periods of variant turnover when more than one variant shared the same niche and high infection rate was observed, which was dependent on the local variants and time. Coinfections might occur at a higher frequency than customarily reported, but they are often ignored as only the consensus sequence is reported for lineage identification.

Genome sequence diversity of SARS-CoV-2 in Serbia: insights gained from a 3-year pandemic study

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for the COVID-19 pandemic, has been evolving rapidly causing emergence of new variants and health uncertainties. Monitoring the evolution of the virus was of the utmost importance for public health interventions and the development of national and global mitigation strategies. Here, we report national data on the emergence of new variants, their distribution, and dynamics in a 3-year study conducted from March 2020 to the end of January 2023 in the Republic of Serbia. Nasopharyngeal and oropharyngeal swabs from 2,398 COVID-19-positive patients were collected and sequenced using three different next generation technologies: Oxford Nanopore, Ion Torrent, and DNBSeq. In the subset of 2,107 SARS-CoV-2 sequences which met the quality requirements, detection of mutations, assignment to SARS-CoV-2 lineages, and phylogenetic analysis were performed. During the 3-year period, we detected three variants of concern, namely, Alpha (5.6%), Delta (7.4%), and Omicron (70.3%) and one variant of interest-Omicron recombinant "Kraken" (XBB1.5) (<1%), whereas 16.8% of the samples belonged to other SARS-CoV-2 (sub)lineages. The detected SARS-CoV-2 (sub)lineages resulted in eight COVID-19 pandemic waves in Serbia, which correspond to the pandemic waves reported in Europe and the United States. Wave dynamics in Serbia showed the most resemblance with the profile of pandemic waves in southern Europe, consistent with the southeastern European location of Serbia. The samples were assigned to sixteen SARS-CoV-2 Nextstrain clades: 20A, 20B, 20C, 20D, 20E, 20G, 20I, 21J, 21K, 21L, 22A, 22B, 22C, 22D, 22E, and 22F and six different Omicron recombinants (XZ, XAZ, XAS, XBB, XBF, and XBK). The 10 most common mutations detected in the coding and untranslated regions of the SARS-CoV-2 genomes included four mutations affecting the spike protein (S:D614G, S:T478K, S:P681H, and S:S477N) and one mutation at each of the following positions: 5'-untranslated region (5'UTR:241); N protein (N:RG203KR); NSP3 protein (NSP3:F106F); NSP4 protein (NSP4:T492I); NSP6 protein (NSP6: S106/G107/F108 - triple deletion), and NSP12b protein (NSP12b:P314L). This national-level study is the most comprehensive in terms of sequencing and genomic surveillance of SARS-CoV-2 during the pandemic in Serbia, highlighting the importance of establishing and maintaining good national practice for monitoring SARS-CoV-2 and other viruses circulating worldwide.

Real-time evaluation of signal accuracy in wastewater surveillance of pathogens with high rates of mutation

Wastewater surveillance of coronavirus disease 2019 (COVID-19) commonly applies reverse transcription-quantitative polymerase chain reaction (RT-qPCR) to quantify severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA concentrations in wastewater over time. In most applications worldwide, maximal sensitivity and specificity of RT-qPCR has been achieved, in part, by monitoring two or more genomic loci of SARS-CoV-2. In Ontario, Canada, the provincial Wastewater Surveillance Initiative reports the average copies of the CDC N1 and N2 loci normalized to the fecal biomarker pepper mild mottle virus. In November 2021, the emergence of the Omicron variant of concern, harboring a C28311T mutation within the CDC N1 probe region, challenged the accuracy of the consensus between the RT-qPCR measurements of the N1 and N2 loci of SARS-CoV-2. In this study, we developed and applied a novel real-time dual loci quality assurance and control framework based on the relative difference between the loci measurements to the City of Ottawa dataset to identify a loss of sensitivity of the N1 assay in the period from July 10, 2022 to January 31, 2023. Further analysis via sequencing and allele-specific RT-qPCR revealed a high proportion of mutations C28312T and A28330G during the study period, both in the City of Ottawa and across the province. It is hypothesized that nucleotide mutations in the probe region, especially A28330G, led to inefficient annealing, resulting in reduction in sensitivity and accuracy of the N1 assay. This study highlights the importance of implementing quality assurance and control criteria to continually evaluate, in near real-time, the accuracy of the signal produced in wastewater surveillance applications that rely on detection of pathogens whose genomes undergo high rates of mutation.

S-Gene Target Failure as an Effective Tool for Tracking the Emergence of Dominant SARS-CoV-2 Variants in Switzerland and Liechtenstein, Including Alpha, Delta, and Omicron BA.1, BA.2, and BA.4/BA.5

During the SARS-CoV-2 pandemic, the Dr. Risch medical group employed the multiplex TaqPathTM COVID-19 CE-IVD RT-PCR Kit for large-scale routine diagnostic testing in Switzerland and the principality of Liechtenstein. The TaqPath Kit is a widely used multiplex assay targeting three genes (i.e., ORF1AB, N, S). With emergence of the B.1.1.7 (Alpha) variant, a diagnostic flaw became apparent as the amplification of the S-gene target was absent in these samples due to a deletion (ΔH69/V70) in the Alpha variant genome. This S-gene target failure (SGTF) was the earliest indication of a new variant emerging and was also observed in subsequent variants such as Omicron BA.1 and BA4/BA.5. The Delta variant and Omicron BA.2 did not present with SGTF. From September 2020 to November 2022, we investigated the applicability of the SGTF as a surrogate marker for emerging variants such as B.1.1.7, B.1.617.2 (Delta), and Omicron BA.1, BA.2, and BA.4/BA.5 in samples with cycle threshold (Ct) values < 30. Next to true SGTF-positive and SGTF-negative samples, there were also samples presenting with delayed-type S-gene amplification (higher Ct value for S-gene than ORF1ab gene). Among these, a difference of 3.8 Ct values between the S- and ORF1ab genes was found to best distinguish between "true" SGTF and the cycle threshold variability of the assay. Samples above the cutoff were subsequently termed partial SGTF (pSGTF). Variant confirmation was performed by whole-genome sequencing (Oxford Nanopore Technology, Oxford, UK) or mutation-specific PCR (TIB MOLBIOL). In total, 17,724 (7.4%) samples among 240,896 positives were variant-confirmed, resulting in an overall sensitivity and specificity of 93.2% [92.7%, 93.7%] and 99.3% [99.2%, 99.5%], respectively. Sensitivity was increased to 98.2% [97.9% to 98.4%] and specificity lowered to 98.9% [98.6% to 99.1%] when samples with pSGTF were included. Furthermore, weekly logistic growth rates (α) and sigmoid's midpoint (t0) were calculated based on SGTF data and did not significantly differ from calculations based on comprehensive data from GISAID. The SGTF therefore allowed for a valid real-time estimate for the introduction of all dominant variants in Switzerland and Liechtenstein.

Mobilisation and analyses of publicly available SARS-CoV-2 data for pandemic responses

The COVID-19 pandemic has seen large-scale pathogen genomic sequencing efforts, becoming part of the toolbox for surveillance and epidemic research. This resulted in an unprecedented level of data sharing to open repositories, which has actively supported the identification of SARS-CoV-2 structure, molecular interactions, mutations and variants, and facilitated vaccine development and drug reuse studies and design. The European COVID-19 Data Platform was launched to support this data sharing, and has resulted in the deposition of several million SARS-CoV-2 raw reads. In this paper we describe (1) open data sharing, (2) tools for submission, analysis, visualisation and data claiming (e.g. ORCiD), (3) the systematic analysis of these datasets, at scale via the SARS-CoV-2 Data Hubs as well as (4) lessons learnt. This paper describes a component of the Platform, the SARS-CoV-2 Data Hubs, which enable the extension and set up of infrastructure that we intend to use more widely in the future for pathogen surveillance and pandemic preparedness.

A framework for automated scalable designation of viral pathogen lineages from genomic data

Pathogen lineage nomenclature systems are a key component of effective communication and collaboration for researchers and public health workers. Since February 2021, the Pango dynamic lineage nomenclature for SARS-CoV-2 has been sustained by crowdsourced lineage proposals as new isolates were sequenced. This approach is vulnerable to time-critical delays as well as regional and personal bias. Here we developed a simple heuristic approach for dividing phylogenetic trees into lineages, including the prioritization of key mutations or genes. Our implementation is efficient on extremely large phylogenetic trees consisting of millions of sequences and produces similar results to existing manually curated lineage designations when applied to SARS-CoV-2 and other viruses including chikungunya virus, Venezuelan equine encephalitis virus complex and Zika virus. This method offers a simple, automated and consistent approach to pathogen nomenclature that can assist researchers in developing and maintaining phylogeny-based classifications in the face of ever-increasing genomic datasets.

Molecular Epidemiology of Omicron CH.1.1 Lineage: Genomic and Phenotypic Data Perspective

BACKGROUND

The Omicron variant (B.1.1.529 lineage) of SARS-CoV-2 represents a substantial global health challenge due to its high transmissibility and potential resistance to immunity from vaccines or previous infections. Among the rapidly evolving Omicron lineages, the BA.2.75 and the emerging CH.1.1 have garnered attention. While BA.2.75 is marked by mutations that may enhance immune evasion, CH.1.1 is distinguished by the S: L452R mutation, linked to increased pathogenicity and transmission. Initially identified in India by the end of 2021, these variants have exhibited global dissemination, signaling an urgent need to track and analyze their progression.

METHODS

In this study, the genomic and geographical distribution data of CH.1.1 were collected from the Global Initiative on Sharing Avian Influenza Data (GISAID), PANGOLIN, CoV-Spectrum, and NextStrain databases. Due to the unavailability of epidemiological and genomic data of the CH.1.1 lineage, PubMed and ScienceDirect were used as sources of the phenotypic data of the lineage variations. Amino acid variations utilized in the data mining included S: R346T, S: K444T, S: L452R, and S: F486S.

RESULTS

The current epidemiological data indicate that CH.1.1 is more likely to become one of the dominant spreading lineages in the United Kingdom, New Zealand, Australia, and the United States based on a 32% growth advantage, present CH.1.1 lineage cases number, and the amino acid variation's impact.

CONCLUSION

A significant increase in the newly detected lineage CH.1.1 is highly anticipated. The rise in the detected sequences number from 13,231 on January 21, 2023, to 23,181 on February 6, 2023, supports the prediction and growth advantage of the lineage detected cases. Increases in viral transmissibility caused by higher affinity to ACE2 receptors and immune evasion are deduced from amino acid variations analyzed in the study.

Evolution of enhanced innate immune suppression by SARS-CoV-2 Omicron subvariants

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) human adaptation resulted in distinct lineages with enhanced transmissibility called variants of concern (VOCs). Omicron is the first VOC to evolve distinct globally dominant subvariants. Here we compared their replication in human cell lines and primary airway cultures and measured host responses to infection. We discovered that subvariants BA.4 and BA.5 have improved their suppression of innate immunity when compared with earlier subvariants BA.1 and BA.2. Similarly, more recent subvariants (BA.2.75 and XBB lineages) also triggered reduced innate immune activation. This correlated with increased expression of viral innate antagonists Orf6 and nucleocapsid, reminiscent of VOCs Alpha to Delta. Increased Orf6 levels suppressed host innate responses to infection by decreasing IRF3 and STAT1 signalling measured by transcription factor phosphorylation and nuclear translocation. Our data suggest that convergent evolution of enhanced innate immune antagonist expression is a common pathway of human adaptation and link Omicron subvariant dominance to improved innate immune evasion.

Sotrovimab: A Review of Its Efficacy against SARS-CoV-2 Variants

Among the anti-Spike monoclonal antibodies (mAbs), the S-309 derivative sotrovimab was the most successful in having the longest temporal window of clinical use, showing a high degree of resiliency to SARS-CoV-2 evolution interrupted only by the appearance of the BA.2.86* variant of interest (VOI). This success undoubtedly reflects rational selection to target a highly conserved epitope in coronavirus Spike proteins. We review here the efficacy of sotrovimab against different SARS-CoV-2 variants in outpatients and inpatients, discussing both randomized controlled trials and real-world evidence. Although it could not be anticipated at the time of its development and introduction, sotrovimab's use in immunocompromised individuals who harbor large populations of variant viruses created the conditions for its eventual demise, as antibody selection and viral evolution led to its eventual withdrawal due to inefficacy against later variant lineages. Despite this, based on observational and real-world data, some authorities have continued to promote the use of sotrovimab, but the lack of binding to newer variants strongly argues for the futility of continued use. The story of sotrovimab highlights the power of modern biomedical science to generate novel therapeutics while also providing a cautionary tale for the need to devise strategies to minimize the emergence of resistance to antibody-based therapeutics.

Identification of a highly conserved neutralizing epitope within the RBD region of diverse SARS-CoV-2 variants

The constant emergence of SARS-CoV-2 variants continues to impair the efficacy of existing neutralizing antibodies, especially XBB.1.5 and EG.5, which showed exceptional immune evasion properties. Here, we identify a highly conserved neutralizing epitope targeted by a broad-spectrum neutralizing antibody BA7535, which demonstrates high neutralization potency against not only previous variants, such as Alpha, Beta, Gamma, Delta and Omicron BA.1-BA.5, but also more recently emerged Omicron subvariants, including BF.7, CH.1.1, XBB.1, XBB.1.5, XBB.1.9.1, EG.5. Structural analysis of the Omicron Spike trimer with BA7535-Fab using cryo-EM indicates that BA7535 recognizes a highly conserved cryptic receptor-binding domain (RBD) epitope, avoiding most of the mutational hot spots in RBD. Furthermore, structural simulation based on the interaction of BA7535-Fab/RBD complexes dissects the broadly neutralizing effect of BA7535 against latest variants. Therapeutic and prophylactic treatment with BA7535 alone or in combination with BA7208 protected female mice from the circulating Omicron BA.5 and XBB.1 variant infection, suggesting the highly conserved neutralizing epitope serves as a potential target for developing highly potent therapeutic antibodies and vaccines.

CD4+ and CD8+ T cells are required to prevent SARS-CoV-2 persistence in the nasal compartment

SARS-CoV-2 is the causative agent of COVID-19 and continues to pose a significant public health threat throughout the world. Following SARS-CoV-2 infection, virus-specific CD4+ and CD8+ T cells are rapidly generated to form effector and memory cells and persist in the blood for several months. However, the contribution of T cells in controlling SARS-CoV-2 infection within the respiratory tract are not well understood. Using C57BL/6 mice infected with a naturally occurring SARS-CoV-2 variant (B.1.351), we evaluated the role of T cells in the upper and lower respiratory tract. Following infection, SARS-CoV-2-specific CD4+ and CD8+ T cells are recruited to the respiratory tract and a vast proportion secrete the cytotoxic molecule Granzyme B. Using antibodies to deplete T cells prior to infection, we found that CD4+ and CD8+ T cells play distinct roles in the upper and lower respiratory tract. In the lungs, T cells play a minimal role in viral control with viral clearance occurring in the absence of both CD4+ and CD8+ T cells through 28 days post-infection. In the nasal compartment, depletion of both CD4+ and CD8+ T cells, but not individually, results in persistent and culturable virus replicating in the nasal compartment through 28 days post-infection. Using in situ hybridization, we found that SARS-CoV-2 infection persisted in the nasal epithelial layer of tandem CD4+ and CD8+ T cell-depleted mice. Sequence analysis of virus isolates from persistently infected mice revealed mutations spanning across the genome, including a deletion in ORF6. Overall, our findings highlight the importance of T cells in controlling virus replication within the respiratory tract during SARS-CoV-2 infection.

Systematic Genomic Surveillance of SARS-CoV-2 at Xiamen International Airport and the Port of Xiamen Reveals the Importance of Incoming Travelers in Lineage Diversity

Sever Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) is still a threat to human health globally despite the World Health Organization (WHO) announcing the end of the COVID-19 pandemic. Continued surveillance of SARS-CoV-2 at national borders would be helpful in understanding the epidemics of novel imported variants and updating local strategies for disease prevention and treatment. This study focuses on the surveillance of imported SARS-CoV-2 variants among travelers entering Xiamen International Airport and the Port of Xiamen from February to August 2023. A total of 97 imported SARS-CoV-2 sequences among travelers from 223 cases collected from 12 different countries and regions were identified by real-time RT-PCR. Next-generation sequencing was used to generate high-quality complete sequences for phylogenetic and population dynamic analysis. The study revealed a dominant shift in variant distribution, in which the XBB subvariant (XBB.1.5, XBB.1.16, XBB.1.9, XBB.2.3, and EG.5.1) accounted for approximately 88.8% of the sequenced samples. In detail, clades 23D and 23E accounted for 26.2% and 21.4% of the sequenced samples, respectively, while clades 23B (13.6%) and 23F (10.7%) took the third and fourth spots in the order of imported sequences, respectively. Additionally, the XBB.2.3 variants were first identified in imported cases from the mainland of Xiamen, China on 27 February 2023. The spatiotemporal analyses of recent viral genome sequences from a limited number of travelers into Xiamen provide valuable insights into the situation surrounding SARS-CoV-2 and highlight the importance of sentinel surveillance of SARS-CoV-2 variants in the national border screening of incoming travelers, which serves as an early warning system for the presence of highly transmissible circulating SARS-CoV-2 lineages.

Stronger and durable SARS-CoV-2 immune response to mRNA vaccines in 5-11 years old children with prior COVID-19

BACKGROUND AND OBJECTIVES

mRNA vaccines elicit a durable humoral response to SARS-CoV-2 in adults, whereas evidence in children is scarce. This study aimed to assess the early and long-term immune response to the mRNA vaccine in children with or without previous SARS-CoV-2 infection.

METHODS

In a multicentre prospective observational study, we profiled the immune response to the Pfizer BioNTech (BNT162b2) vaccine in 5-11-year-old children attending the University Pediatric Hospital of Padua and Bambino-Gesù Hospital in Rome (Italy) from December-2021 to February-2023. Blood samples were collected pre-, 1-, and 6-months after vaccination. Neutralizing antibodies (NAbs) and anti-spike-receptor-binding-domain (anti-S-RBD) IgG titers were analyzed through Plaque Reduction Neutralization Test (PRNT) and chemiluminescent immune-enzymatic assay (CLIA), respectively. Immune cell phenotypes were analyzed by flow cytometry.

RESULTS

Sixty children (26 [43 %] female, median age = 8 years [IQR = 7-10.7]) were enrolled in the study, including 46 children with a laboratory-confirmed previous COVID-19 (SARS-CoV-2-recovered) and 14 SARS-CoV-2-naïve participants defined as the absence of antigen-specific antibodies before vaccination. SARS-CoV-2-recovered participants recorded higher anti-S-RBD IgG and Wild-type and Omicron BA.2 NAbs titers than SARS-CoV-2-naïve participants at both 1- and 6-months after vaccination. Antibody titers correlated with T (Tregs) and B (Bregs) regulatory cell frequencies in SARS-CoV-2-recovered children. Both SARS-CoV-2-recovered and SARS-CoV-2-naïve participants decreased antibody titers by approximately 100 to 250 % from 1 to 6 months. While children with immunocompromising underlying conditions developed immune responses comparable to those of healthy children, solid organ transplant recipients exhibited lower levels of NAbs and anti-S-RBD IgG titers, as well as reduced frequencies of Tregs and Bregs.

CONCLUSIONS

mRNA vaccination triggered a higher production of specific anti-SARS-CoV-2 antibodies along with increased levels of regulatory cells in children with previous SARS-CoV-2 infection up to the following 6 months. These findings provide insights into boosting pre-existing immunity.

V-pipe 3.0: a sustainable pipeline for within-sample viral genetic diversity estimation

The large amount and diversity of viral genomic datasets generated by next-generation sequencing technologies poses a set of challenges for computational data analysis workflows, including rigorous quality control, scaling to large sample sizes, and tailored steps for specific applications. Here, we present V-pipe 3.0, a computational pipeline designed for analyzing next-generation sequencing data of short viral genomes. It is developed to enable reproducible, scalable, adaptable, and transparent inference of genetic diversity of viral samples. By presenting 2 large-scale data analysis projects, we demonstrate the effectiveness of V-pipe 3.0 in supporting sustainable viral genomic data science.

Divergent trajectory of replication and intrinsic pathogenicity of SARS-CoV-2 Omicron post-BA.2/5 subvariants in the upper and lower respiratory tract

BACKGROUND

Earlier Omicron subvariants including BA.1, BA.2, and BA.5 emerged in waves, with a subvariant replacing the previous one every few months. More recently, the post-BA.2/5 subvariants have acquired convergent substitutions in spike that facilitated their escape from humoral immunity and gained ACE2 binding capacity. However, the intrinsic pathogenicity and replication fitness of the evaluated post-BA.2/5 subvariants are not fully understood.

METHODS

We systemically investigated the replication fitness and intrinsic pathogenicity of representative post-BA.2/5 subvariants (BL.1, BQ.1, BQ.1.1, XBB.1, CH.1.1, and XBB.1.5) in weanling (3-4 weeks), adult (8-10 weeks), and aged (10-12 months) mice. In addition, to better model Omicron replication in the human nasal epithelium, we further investigated the replication capacity of the post-BA.2/5 subvariants in human primary nasal epithelial cells.

FINDINGS

We found that the evaluated post-BA.2/5 subvariants are consistently attenuated in mouse lungs but not in nasal turbinates when compared with their ancestral subvariants BA.2/5. Further investigations in primary human nasal epithelial cells revealed a gained replication fitness of XBB.1 and XBB.1.5 when compared to BA.2 and BA.5.2.

INTERPRETATION

Our study revealed that the post-BA.2/5 subvariants are attenuated in lungs while increased in replication fitness in the nasal epithelium, indicating rapid adaptation of the circulating Omicron subvariants in the human populations.

FUNDING

The full list of funding can be found at the Acknowledgements section.

The Era of the FLips: How Spike Mutations L455F and F456L (and A475V) Are Shaping SARS-CoV-2 Evolution

Convergent evolution of the SARS-CoV-2 Spike protein has been mostly driven by immune escape, in particular by escape to the viral infection-neutralizing antibodies (nAbs) elicited by previous infections and/or vaccinations [...].

Molecular biology of SARS-CoV-2 and techniques of diagnosis and surveillance

The World Health Organization (WHO) declared coronavirus disease 2019 (COVID-19), a disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a global pandemic in March 2020. Reverse transcription-polymerase chain reaction (RT-PCR) is the reference technique for molecular diagnosis of SARS-CoV-2 infection. The SARS-CoV-2 virus is constantly mutating, and more transmissible variants have emerged, making genomic surveillance a crucial tool for investigating virus transmission dynamics, detecting novel genetic variants, and assessing mutation impact. The S gene, which encodes the spike protein, is frequently mutated, and it plays an important role in transmissibility. Spike protein mutations affect infectivity and vaccine effectiveness. SARS-CoV-2 variants are tracked using whole genome sequencing (WGS) and S-gene analysis. WGS, Sanger sequencing, and many S-gene-targeted RT-PCR methods have been developed. WGS and Sanger sequencing are standard methods for detecting mutations and can be used to identify known and unknown mutations. Melting curve analysis, endpoint genotyping assay, and S-gene target failure are used in the RT-PCR-based method for the rapid detection of specific mutations in SARS-CoV-2 variants. Therefore, these assays are suitable for high-throughput screening. The combinatorial use of RT-PCR-based assays, Sanger sequencing, and WGS enables rapid and accurate tracking of SARS-CoV-2 variants. In this review, we described RT-PCR-based detection and surveillance techniques for SARS-CoV-2.

The increase in SARS-CoV-2 lineages during 2020-2022 in a state in the Brazilian Northeast is associated with a number of cases

SARS-CoV-2 has caused a high number of deaths in several countries. In Brazil, there were 37 million confirmed cases of COVID-19 and 700,000 deaths caused by the disease. The population size and heterogeneity of the Brazilian population should be considered in epidemiological surveillance due to the varied tropism of the virus. As such, municipalities and states must be factored in for their unique specificities, such as socioeconomic conditions and population distribution. Here, we investigate the spatiotemporal dispersion of emerging SARS-CoV-2 lineages and their dynamics in each microregion from Sergipe state, northeastern Brazil, in the first 3 years of the pandemic. We analyzed 586 genomes sequenced between March 2020 and November 2022 extracted from the GISAID database. Phylogenetic analyses were carried out for each data set to reconstruct evolutionary history. Finally, the existence of a correlation between the number of lineages and infection cases by SARS-CoV-2 was evaluated. Aracaju, the largest city in northeastern Brazil, had the highest number of samples sequenced. This represented 54.6% (320) of the genomes, and consequently, the largest number of lineages identified. Studies also analyzed the relationship between mean lineage distributions and mean monthly infections, daily cases, daily deaths, and hospitalizations of vaccinated and unvaccinated patients. For this, a correlation matrix was created. Results revealed that the increase in the average number of SARS-CoV-2 variants was related to the average number of SARS-CoV-2 cases in both unvaccinated and vaccinated individuals. Thus, our data indicate that it is necessary to maintain epidemiological surveillance, especially in capital cities, since they have a high rate of circulation of resident and non-resident inhabitants, which contributes to the dynamics of the virus.

Low-sugar universal mRNA vaccine against coronavirus variants with deletion of glycosites in the S2 or stem of SARS-CoV-2 spike messenger RNA (mRNA)

Since the outbreak of Severe Acute Respiratory Syndrome Virus-2 (SARS-CoV-2) in 2019, more than 15 million spike protein sequences have been identified, raising a new challenge for the development of a broadly protective vaccine against the various emerging variants. We found that the virus, like most other human viruses, depends on host-made glycans to shield the conserved epitopes on spike protein from immune response and demonstrated that deletion of the glycan shields exposed highly conserved epitopes and elicited broadly protective immune responses. In this study, we identified 17 conserved epitopes from 14 million spike protein sequences and 11 of the conserved epitopes are in the S2 domain, including the six most conserved epitopes in the stem region. We also demonstrated that deletion of the glycosites in the spike messenger RNA (mRNA) S2 domain or the stem region exposed the highly conserved epitopes and elicited broadly protective immune responses, particularly CD-8+ T cell response against various SARS-CoV-2 variants, and other human coronaviruses including MERS, SARS viruses, and those causing common cold.

Immune evasion and membrane fusion of SARS-CoV-2 XBB subvariants EG.5.1 and XBB.2.3

Immune evasion by SARS-CoV-2 paired with immune imprinting from monovalent mRNA vaccines has resulted in attenuated neutralizing antibody responses against Omicron subvariants. In this study, we characterized two new XBB variants rising in circulation - EG.5.1 and XBB.2.3, for their neutralization and syncytia formation. We determined the neutralizing antibody titers in sera of individuals that received a bivalent mRNA vaccine booster, BA.4/5-wave infection, or XBB.1.5-wave infection. Bivalent vaccination-induced antibodies neutralized ancestral D614G efficiently, but to a much less extent, two new EG.5.1 and XBB.2.3 variants. In fact, the enhanced neutralization escape of EG.5.1 appeared to be driven by its key defining mutation XBB.1.5-F456L. Notably, infection by BA.4/5 or XBB.1.5 afforded little, if any, neutralization against EG.5.1, XBB.2.3 and previous XBB variants - especially in unvaccinated individuals, with average neutralizing antibody titers near the limit of detection. Additionally, we investigated the infectivity, fusion activity, and processing of variant spikes for EG.5.1 and XBB.2.3 in HEK293T-ACE2 and CaLu-3 cells but found no significant differences compared to earlier XBB variants. Overall, our findings highlight the continued immune evasion of new Omicron subvariants and, more importantly, the need to reformulate mRNA vaccines to include XBB spikes for better protection.

SARS-CoV-2 Omicron variants BA.4 and BA.5 dominated the fifth COVID-19 epidemiological wave in Mexico

In Mexico, the BA.4 and BA.5 Omicron variants dominated the fifth epidemic wave (summer 2022), superseding BA.2, which had circulated during the inter-wave period. The present study uses genome sequencing and statistical and phylogenetic analyses to examine these variants' abundance, distribution, and genetic diversity in Mexico from April to August 2022. Over 35 % of the sequenced genomes in this period corresponded to the BA.2 variant, 8 % to the BA.4 and 56 % to the BA.5 variant. Multiple subvariants were identified, but the most abundant, BA.2.9, BA.2.12.1, BA.5.1, BA.5.2, BA.5.2.1 and BA.4.1, circulated across the entire country, not forming geographical clusters. Contrastingly, other subvariants exhibited a geographically restricted distribution, most notably in the Southeast region, which showed a distinct subvariant dynamic. This study supports previous results showing that this region may be a significant entry point and contributed to introducing and evolving novel variants in Mexico. Furthermore, a differential distribution was observed for certain subvariants among specific States through time, which may have contributed to the overall increased diversity observed during this wave compared to the previous ones. This study highlights the importance of sustaining genomic surveillance to identify novel variants that may impact public health.

Convergent evolution of SARS-CoV-2 XBB lineages on receptor-binding domain 455-456 synergistically enhances antibody evasion and ACE2 binding

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) XBB lineages have achieved dominance worldwide and keep on evolving. Convergent evolution of XBB lineages on the receptor-binding domain (RBD) L455F and F456L is observed, resulting in variants with substantial growth advantages, such as EG.5, FL.1.5.1, XBB.1.5.70, and HK.3. Here, we show that neutralizing antibody (NAb) evasion drives the convergent evolution of F456L, while the epistatic shift caused by F456L enables the subsequent convergence of L455F through ACE2 binding enhancement and further immune evasion. L455F and F456L evade RBD-targeting Class 1 public NAbs, reducing the neutralization efficacy of XBB breakthrough infection (BTI) and reinfection convalescent plasma. Importantly, L455F single substitution significantly dampens receptor binding; however, the combination of L455F and F456L forms an adjacent residue flipping, which leads to enhanced NAbs resistance and ACE2 binding affinity. The perturbed receptor-binding mode leads to the exceptional ACE2 binding and NAb evasion, as revealed by structural analyses. Our results indicate the evolution flexibility contributed by epistasis cannot be underestimated, and the evolution potential of SARS-CoV-2 RBD remains high.

Comparative study showed that children faced a 78% higher risk of new-onset conditions after they had COVID-19

AIM

Children have largely been unaffected by severe COVID-19 compared to adults, but data suggest that they may have experienced new conditions after developing the disease. We compared outcomes in children who had experienced COVID-19 and healthy controls.

METHODS

A retrospective nested cohort study assessed the incidence rate of new-onset conditions after COVID-19 in children aged 0-14 years. Data were retrieved from an Italian paediatric primary care database linked to Veneto Region registries. Exposed children with a positive nasopharyngeal swab were matched 1:1 with unexposed children who had tested negative. Conditional Cox regression was fitted to estimate the adjusted hazard ratios (aHR) and 95% confidence intervals (CI) for the exposure and outcome associations after adjusting for covariates.

RESULTS

We compared 1656 exposed and 1656 unexposed children from 1 February 2020 to 30 November 2021. The overall excess risk for new-onset conditions after COVID-19 was 78% higher in the exposed than unexposed children. We found significantly higher risks for some new conditions in exposed children, including mental health issues (aHR 1.8, 95% CI 1.1-3.0) and neurological problems (aHR 2.4, 95% CI 1.4-4.1).

CONCLUSION

Exposed children had a 78% higher risk of developing new conditions of interest after COVID-19 than unexposed children.

SARS-CoV-2 nsp3 and nsp4 are minimal constituents of a pore spanning replication organelle

Coronavirus replication is associated with the remodeling of cellular membranes, resulting in the formation of double-membrane vesicles (DMVs). A DMV-spanning pore was identified as a putative portal for viral RNA. However, the exact components and the structure of the SARS-CoV-2 DMV pore remain to be determined. Here, we investigate the structure of the DMV pore by in situ cryo-electron tomography combined with subtomogram averaging. We identify non-structural protein (nsp) 3 and 4 as minimal components required for the formation of a DMV-spanning pore, which is dependent on nsp3-4 proteolytic cleavage. In addition, we show that Mac2-Mac3-DPUP-Ubl2 domains are critical for nsp3 oligomerization and crown integrity which influences membrane curvature required for biogenesis of DMVs. Altogether, SARS-CoV-2 nsp3-4 have a dual role by driving the biogenesis of replication organelles and assembly of DMV-spanning pores which we propose here to term replicopores.

Analysis of SARS-CoV-2 mutations associated with resistance to therapeutic monoclonal antibodies that emerge after treatment

The mutation rate of the Omicron sublineage has led to baseline resistance against all previously authorized anti-Spike monoclonal antibodies (mAbs). Nevertheless, in case more antiviral mAbs will be authorized in the future, it is relevant to understand how frequently treatment-emergent resistance has emerged so far, under different combinations and in different patient subgroups. We report the results of a systematic review of the medical literature for case reports and case series for treatment-emergent immune escape, which is defined as emergence of a resistance-driving mutation in at least 20% of sequences in a given host at a given timepoint. We identified 32 publications detailing 216 cases that included different variants of concern (VOC) and found that the incidence of treatment emergent-resistance ranged from 10% to 50%. Most of the treatment-emergent resistance events occurred in immunocompromised patients. Interestingly, resistance also emerged against cocktails of two mAbs, albeit at lower frequencies. The heterogenous therapeutic management of those cases doesn't allow inferences about the clinical outcome in patients with treatment-emergent resistance. Furthermore, we noted a temporal correlation between the introduction of mAb therapies and a subsequent increase in SARS-CoV-2 sequences across the globe carrying mutations conferring resistance to that mAb, raising concern as to whether these had originated in mAb-treated individuals. Our findings confirm that treatment-emergent immune escape to anti-Spike mAbs represents a frequent and concerning phenomenon and suggests that these are associated with mAb use in immunosuppressed hosts.

CoVEx: SARS-CoV-2 Mutation Explorer for genomic surveillance

Effective management of emerging diseases relies on timely pathogen identification and monitoring. The emergence of COVID-19 in December 2019, rapidly evolved into a global pandemic, with millions of cases and deaths reported worldwide. The accumulation of SARS-CoV-2 genomes provided unprecedented opportunities for studying the virus's evolutionary dynamics, understanding the impact of mutations, and identifying emerging Variants of Interest (VOIs) and Variants of Concern (VOCs). During the COVID-19 pandemic, health systems faced challenges in promptly detecting such variants and timely notifying. To facilitate the continuous monitoring of mutations, various initiatives and open-source pipelines have been established. However, these platforms often lack integration for conducting user sequence analysis and comparing it with publicly reported data on platforms like GISAID. Here, we present CoVEx, an easy-to-use tool for analyzing and visualizing SARS-CoV-2 variant sequences obtained using Illumina sequencing technology. CoVEx integrates quality control, alignment, genome annotation, lineage designation, and mutation analysis tools. Implemented in Python, CoVEx also has a mutation explorer feature that generates interactive graphs summarising identified mutations in an intuitive manner. Similarly, it leverages the Outbreak.info package to create heatmaps highlighting the mutations associated with designated Pangolin lineages. Furthermore, by comparing mutation profiles against GISAID data, CoVEx offers valuable insights into the prevalence and distribution of mutations worldwide. We validated CoVEx using raw sequence data (n = 108) and demonstrated its accuracy in assembling sequences and predicting Pangolin and Nextclade Pango lineages. Notably, the tool revealed the emergence of a previously unreported mutation, ORF1a:I2501T, within the Costa Rica GN.1 lineage. This finding highlights CoVEx's capability to identify novel mutations in the different lineages, providing valuable information to researchers and public health decision makers. CoVEx and documentation are freely available on GitLab: https://gitlab.com/CNCA_CeNAT/covex.

Spike Protein Mutation-Induced Changes in the Kinetic and Thermodynamic Behavior of Its Receptor Binding Domains Explain Their Higher Propensity to Attain Open States in SARS-CoV-2 Variants of Concern

Spike (S) protein opening in SARS-CoV-2 controls the accessibility of its receptor binding domains (RBDs) to host receptors and immune recognition. Along the evolution of SARS-CoV-2 to its variants of concern (VOC)-alpha, beta, gamma, delta, and omicron-their S proteins showed a higher propensity to attain open states. Deciphering how mutations in S protein can shape its conformational dynamics will contribute to the understanding of viral host tropism. Here using microsecond-scale multiple molecular dynamics simulations (MDS), we provide insights into the kinetic and thermodynamic contributions of these mutations to RBD opening pathways in S proteins of SARS-CoV-2 VOCs. Mutational effects were analyzed using atomistic (i) equilibrium MDS of closed and open states of S proteins and (ii) nonequilibrium MDS for closed-to-open transitions. In MDS of closed or open states, RBDs in S proteins of VOCs showed lower thermodynamic stability with higher kinetic fluctuations, compared to S proteins of ancestral SARS-CoV-2. For closed-to-open transitions in S proteins of VOCs, we observed apparently faster RBD opening with a 1.5-2-fold decrease in the thermodynamic free-energy barrier (ΔGclosed→open). Saturation mutagenesis studies highlighted S protein mutations that may control its conformational dynamics and presentation to host receptors.

The return of the "Mistigri" (virus adaptative gain by gene loss) through the SARS-CoV-2 XBB.1.5 chimera that predominated in 2023

Severe acute respiratory syndrome coronavirus 2 XBB.1.5 is the first recombinant lineage to predominate at the country and global scales. Very interestingly, like the Marseille-4B subvariant (or B.1.160) and the pandemic variant B.1.1.7 (or Alpha) previously, it has its ORF8 gene inactivated by a stop codon. We aimed here to study the distribution of stop codons in ORF8 of XBB.1.5 and non-XBB.1.5 genomes. We identified that a stop codon was present at 89 (74%) ORF8 codons in ≥1 of 15 222 404 genomes available in GISAID. The mean proportion of genomes with a stop codon per codon was 0.11% (range, 0%-7.8%). In addition, a stop codon was detected at 15 (12%) codons in at least 1000 genomes. These 15 codons are notably located on seven stem-loop hairpin regions and in the signal peptide region for the case of the XBB.1.5 lineage (codon 8). Thus, it is very likely that stop codons in ORF8 gene contributed on at least three occasions and independently during the pandemic to the evolutionary success of a lineage that became transiently predominant. Such association of gene loss with evolutionary success, which suits the recently described Mistigri rule, is an important biological phenomenon very unknown in virology while largely described in cellular organisms.

RCoV19: A One-stop Hub for SARS-CoV-2 Genome Data Integration, Variant Monitoring, and Risk Pre-warning

The Resource for Coronavirus 2019 (RCoV19) is an open-access information resource dedicated to providing valuable data on the genomes, mutations, and variants of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In this updated implementation of RCoV19, we have made significant improvements and advancements over the previous version. Firstly, we have implemented a highly refined genome data curation model. This model now features an automated integration pipeline and optimized curation rules, enabling efficient daily updates of data in RCoV19. Secondly, we have developed a global and regional lineage evolution monitoring platform, alongside an outbreak risk pre-warning system. These additions provide a comprehensive understanding of SARS-CoV-2 evolution and transmission patterns, enabling better preparedness and response strategies. Thirdly, we have developed a powerful interactive mutation spectrum comparison module. This module allows users to compare and analyze mutation patterns, assisting in the detection of potential new lineages. Furthermore, we have incorporated a comprehensive knowledgebase on mutation effects. This knowledgebase serves as a valuable resource for retrieving information on the functional implications of specific mutations. In summary, RCoV19 serves as a vital scientific resource, providing access to valuable data, relevant information, and technical support in the global fight against COVID-19. The complete contents of RCoV19 are available to the public at https://ngdc.cncb.ac.cn/ncov/.