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

Applying the digital data and the bioinformatics tools in SARS-CoV-2 research

Bioinformatics has been playing a crucial role in the scientific progress to fight against the pandemic of the coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The advances in novel algorithms, mega data technology, artificial intelligence and deep learning assisted the development of novel bioinformatics tools to analyze daily increasing SARS-CoV-2 data in the past years. These tools were applied in genomic analyses, evolutionary tracking, epidemiological analyses, protein structure interpretation, studies in virus-host interaction and clinical performance. To promote the in-silico analysis in the future, we conducted a review which summarized the databases, web services and software applied in SARS-CoV-2 research. Those digital resources applied in SARS-CoV-2 research may also potentially contribute to the research in other coronavirus and non-coronavirus viruses.

Fast end-to-end surface interpretation of SARS-CoV-2 variants by differentiable molecular surface interaction fingerprinting method

Confronting the challenge of persistent mutations of SARS-CoV-2, researchers have turned to deep learning methods to predict the mutated structures of spike proteins and to hypothesize potential changes in their structures and drug efficacies. However, limited works are focused on the surface learning of spike proteins even though their biological functions are usually defined by the geometric and chemical features of 3D molecular surfaces. In addition, the current used geometric deep learning methods are based on mesh representations of proteins to identify potential binding targets for drugs. However, the use of meshes has limitations and is not applicable for many important tasks in molecular biology. To address these limitations, we adopt the differentiable molecular surface interaction fingerprinting (dMaSIF) method which is based on the 3D point clouds and a novel efficient geometric convolutional layer to fast predict the interaction sites on the protein surface. The different binding site patterns for Delta, Omicron and its subvariants are clearly visualized. We find that Delta and Omicron show the similar surface binding patterns while BA.2, BA.2.13, BA.3 and BA.4 present similar ones. BA.4 possesses higher positive interaction site ratio than the others which may account for its higher transmission and infection among humans. In addition, the positive interaction site ratios of BA.2, BA.2.13, BA.3 are higher than Delta and Omicron, which are accordant with their transmission and infection rates. Hopefully our work offers a new effective route to analyze the protein-protein interaction for the SARS-CoV-2 variants.

Structural basis of white-tailed deer, Odocoileus virginianus, ACE2 recognizing all the SARS-CoV-2 variants of concern with high affinity

SARS-CoV-2 has been expanding its host range, among which the white-tailed deer (WTD), Odocoileus virginianus, became the first wildlife species infected on a large scale and might serve as a host reservoir for variants of concern (VOCs) in case no longer circulating in humans. In this study, we comprehensively assessed the binding of the WTD angiotensin-converting enzyme 2 (ACE2) receptor to the spike (S) receptor-binding domains (RBDs) from the SARS-CoV-2 prototype (PT) strain and multiple variants. We found that WTD ACE2 could be broadly recognized by all of the tested RBDs. We further determined the complex structures of WTD ACE2 with PT, Omicron BA.1, and BA.4/5 S trimer. Detailed structural comparison revealed the important roles of RBD residues on 486, 498, and 501 sites for WTD ACE2 binding. This study deepens our understanding of the interspecies transmission mechanisms of SARS-CoV-2 and further addresses the importance of constant monitoring on SARS-CoV-2 infections in wild animals. IMPORTANCE Even if we manage to eliminate the virus among humans, it will still circulate among wildlife and continuously be transmitted back to humans. A recent study indicated that WTD may serve as reservoir for nearly extinct SARS-CoV-2 strains. Therefore, it is critical to evaluate the binding abilities of SARS-CoV-2 variants to the WTD ACE2 receptor and elucidate the molecular mechanisms of binding of the RBDs to assess the risk of spillback events.

Unraveling the Dynamics of Omicron (BA.1, BA.2, and BA.5) Waves and Emergence of the Deltacton Variant: Genomic Epidemiology of the SARS-CoV-2 Epidemic in Cyprus (Oct 2021-Oct 2022)

Commencing in December 2019 with the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), three years of the coronavirus disease 2019 (COVID-19) pandemic have transpired. The virus has consistently demonstrated a tendency for evolutionary adaptation, resulting in mutations that impact both immune evasion and transmissibility. This ongoing process has led to successive waves of infections. This study offers a comprehensive assessment spanning genetic, phylogenetic, phylodynamic, and phylogeographic dimensions, focused on the trajectory of the SARS-CoV-2 epidemic in Cyprus. Based on a dataset comprising 4700 viral genomic sequences obtained from affected individuals between October 2021 and October 2022, our analysis is presented. Over this timeframe, a total of 167 distinct lineages and sublineages emerged, including variants such as Delta and Omicron (1, 2, and 5). Notably, during the fifth wave of infections, Omicron subvariants 1 and 2 gained prominence, followed by the ascendancy of Omicron 5 in the subsequent sixth wave. Additionally, during the fifth wave (December 2021-January 2022), a unique set of Delta sequences with genetic mutations associated with Omicron variant 1, dubbed "Deltacron", was identified. The emergence of this phenomenon initially evoked skepticism, characterized by concerns primarily centered around contamination or coinfection as plausible etiological contributors. These hypotheses were predominantly disseminated through unsubstantiated assertions within the realms of social and mass media, lacking concurrent scientific evidence to validate their claims. Nevertheless, the exhaustive molecular analyses presented in this study have demonstrated that such occurrences would likely lead to a frameshift mutation-a genetic aberration conspicuously absent in our provided sequences. This substantiates the accuracy of our initial assertion while refuting contamination or coinfection as potential etiologies. Comparable observations on a global scale dispelled doubt, eventually leading to the recognition of Delta-Omicron variants by the scientific community and their subsequent monitoring by the World Health Organization (WHO). As our investigation delved deeper into the intricate dynamics of the SARS-CoV-2 epidemic in Cyprus, a discernible pattern emerged, highlighting the major role of international connections in shaping the virus's local trajectory. Notably, the United States and the United Kingdom were the central conduits governing the entry and exit of the virus to and from Cyprus. Moreover, notable migratory routes included nations such as Greece, South Korea, France, Germany, Brazil, Spain, Australia, Denmark, Sweden, and Italy. These empirical findings underscore that the spread of SARS-CoV-2 within Cyprus was markedly influenced by the influx of new, highly transmissible variants, triggering successive waves of infection. This investigation elucidates the emergence of new waves of infection subsequent to the advent of highly contagious and transmissible viral variants, notably characterized by an abundance of mutations localized within the spike protein. Notably, this discovery decisively contradicts the hitherto hypothesis of seasonal fluctuations in the virus's epidemiological dynamics. This study emphasizes the importance of meticulously examining molecular genetics alongside virus migration patterns within a specific region. Past experiences also emphasize the substantial evolutionary potential of viruses such as SARS-CoV-2, underscoring the need for sustained vigilance. However, as the pandemic's dynamics continue to evolve, a balanced approach between caution and resilience becomes paramount. This ethos encourages an approach founded on informed prudence and self-preservation, guided by public health authorities, rather than enduring apprehension. Such an approach empowers societies to adapt and progress, fostering a poised confidence rooted in well-founded adaptation.

Broad protective RBD heterotrimer vaccines neutralize SARS-CoV-2 including Omicron sub-variants XBB/BQ.1.1/BF.7

SARS-CoV-2 variants with severe immune evasion are a major challenge for COVID-19 prevention, especially the circulating Omicron XBB/BQ.1.1/BF.7 strains. Thus, the next-generation of broad-spectrum vaccines are urgently needed. Previously, we developed a COVID-19 protein subunit vaccine, ZF2001, based on the RBD-homodimer as the immunogen. To adapt SARS-CoV-2 variants, we developed chimeric RBD-heterodimers to induce broad immune responses. In this study, we further explored the concept of tandem RBD homotrimer and heterotrimer. Prototype SARS-CoV-2 RBD-homotrimer, prototype-Delta-BA.1 (PDO) RBD-heterotrimer and Delta-BA.2-BA.5 (DBA2BA5) RBD-heterotrimer were designed. Biochemical and cryo-EM structural characterization demonstrated total epitope exposure of the RBD-trimers. In mouse experiments, PDO and DBA2BA5 elicited broad SARS-CoV-2 neutralization. Potent protection against SARS-CoV-2 variants was observed in challenge assays and was correlated with neutralizing antibody titer. This study validated the design strategy of tandem RBD-heterotrimers as multivalent immunogens and presented a promising vaccine candidate, DBA2BA5, eliciting broad-spectrum immune responses, including against the circulating XBB/BF.7/BQ.1.1.

Dynamical modelling of viral infection and cooperative immune protection in COVID-19 patients

Once challenged by the SARS-CoV-2 virus, the human host immune system triggers a dynamic process against infection. We constructed a mathematical model to describe host innate and adaptive immune response to viral challenge. Based on the dynamic properties of viral load and immune response, we classified the resulting dynamics into four modes, reflecting increasing severity of COVID-19 disease. We found the numerical product of immune system's ability to clear the virus and to kill the infected cells, namely immune efficacy, to be predictive of disease severity. We also investigated vaccine-induced protection against SARS-CoV-2 infection. Results suggested that immune efficacy based on memory T cells and neutralizing antibody titers could be used to predict population vaccine protection rates. Finally, we analyzed infection dynamics of SARS-CoV-2 variants within the construct of our mathematical model. Overall, our results provide a systematic framework for understanding the dynamics of host response upon challenge by SARS-CoV-2 infection, and this framework can be used to predict vaccine protection and perform clinical diagnosis.

Selection for immune evasion in SARS-CoV-2 revealed by high-resolution epitope mapping and sequence analysis

Here, we exploit a deep serological profiling strategy coupled with an integrated, computational framework for the analysis of SARS-CoV-2 humoral immune responses. Applying a high-density peptide array (HDPA) spanning the entire proteomes of SARS-CoV-2 and endemic human coronaviruses allowed identification of B cell epitopes and relate them to their evolutionary and structural properties. We identify hotspots of pre-existing immunity and identify cross-reactive epitopes that contribute to increasing the overall humoral immune response to SARS-CoV-2. Using a public dataset of over 38,000 viral genomes from the early phase of the pandemic, capturing both inter- and within-host genetic viral diversity, we determined the evolutionary profile of epitopes and the differences across proteins, waves, and SARS-CoV-2 variants. Lastly, we show that mutations in spike and nucleocapsid epitopes are under stronger selection between than within patients, suggesting that most of the selective pressure for immune evasion occurs upon transmission between hosts.

Digital PCR Discriminates between SARS-CoV-2 Omicron Variants and Immune Escape Mutations

As severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to evolve, mutations arise that will allow the virus to evade immune defenses and therapeutics. Assays that can identify these mutations can be used to guide personalized patient treatment plans. Digital PCR (dPCR) is a fast and reliable complement to whole-genome sequencing that can be used to discriminate single nucleotide polymorphisms (SNPs) in template molecules. Here, we developed a panel of SARS-CoV-2 dPCR assays and demonstrate its applications for typing variant lineages and therapeutic monoclonal antibody resistance. We first designed multiplexed dPCR assays for SNPs located at residue 3395 in the orf1ab gene that differentiate the Delta, Omicron BA.1, and Omicron BA.2 lineages. We demonstrate their effectiveness on 596 clinical saliva specimens that were sequence verified using Illumina whole-genome sequencing. Next, we developed dPCR assays for spike mutations R346T, K444T, N460K, F486V, and F486S, which are associated with host immune evasion and reduced therapeutic monoclonal antibody efficacy. We demonstrate that these assays can be run individually or multiplexed to detect the presence of up to 4 SNPs in a single assay. We perform these dPCR assays on 81 clinical saliva SARS-CoV-2-positive specimens and properly identify mutations in Omicron subvariants BA.2.75.2, BM.1.1, BN.1, BF.7, BQ.1, BQ.1.1, and XBB. Thus, dPCR could serve as a useful tool to determine if clinical specimens contain therapeutically relevant mutations and inform patient treatment. IMPORTANCE Spike mutations in the SARS-CoV-2 genome confer resistance to therapeutic monoclonal antibodies. Authorization for treatment options is typically guided by general trends of variant prevalence. For example, bebtelovimab is no longer authorized for emergency use in the United States due to the increased prevalence of antibody-resistant BQ.1, BQ.1.1, and XBB Omicron subvariants. However, this blanket approach limits access to life-saving treatment options to patients who are otherwise infected with susceptible variants. Digital PCR assays targeting specific mutations can complement whole-genome sequencing approaches to genotype the virus. In this study, we demonstrate the proof of concept that dPCR can be used to type lineage defining and monoclonal antibody resistance-associated mutations in saliva specimens. These findings show that digital PCR could be used as a personalized diagnostic tool to guide individual patient treatment.

Clinical features of COVID-19 in Italian outpatient children and adolescents during Parental, Delta, and Omicron waves: a prospective, observational, cohort study

Introduction

COVID-19 features changed with the Omicron variant of SARS-CoV-2 in adults. This study aims to describe COVID-19 symptoms in children and adolescents during the Parental, Delta, and Omicron eras.

Methods

A single-centre, prospective observational study was conducted on individuals aged 0-20 years attending the University Hospital of Padua (Italy) from April 2020 to December 2022. COVID-19 cases were defined by positive SARS-CoV-2 molecular detection and/or serology; patient/family symptoms and virological positivity were considered to determine the infection onset. Variables were summarized and compared using appropriate tests of descriptive statistics.

Results

A total of 509 cases [46% female, median age eight years (IQR: 4-12)] were studied. Three-hundred-eighty-seven (76%), 52 (10%), and 70 (14%) subjects experienced COVID-19 during the Parental, Delta, and Omicron waves, respectively. All subjects developed an asymptomatic/mild COVID-19. Overall, the most frequent symptoms were fever (47%) and rhinitis (21%), which showed a significant increasing incidence from the Parental to Omicron waves (p < 0.001). Conversely, diarrhea was most common during the pre-Omicron eras (p = 0.03). Stratifying symptoms according to the age group, fever, rhinitis, and skin rashes were observed more frequently among infants/toddlers; conversely, fatigue was more common in children older than five years. The duration of symptoms was similar across different SARS-CoV-2 variants of concern (VOCs); conversely, the number of symptoms varied according to the age group (p < 0.0001).

Discussion

This study showed differences in COVID-19 clinical presentation among infants, children, and adolescents and confirmed Omicron infection is more likely to be associated with upper respiratory symptoms. However, further population-based studies are needed to support these findings. In addition, active surveillance will play a crucial role in assessing the disease severity of future VOCs.

Cardiac Function Evaluation after SARS-CoV-2 mRNA Vaccination in Children and Adolescents: A Prospective Speckle-Tracking Echocardiography Study

BACKGROUND

Possible cardiac impairment after SARS-CoV-2 mRNA vaccination is a common driver of parental vaccine hesitancy. We performed a comprehensive echocardiographic evaluation of biventricular function in vaccinated children with or without previous COVID-19 compared to healthy controls.

METHODS

We conducted a single-center, prospective, case-control study enrolling children and adolescents aged 5-18 years attending the pediatric clinic of the University Hospital of Padua from April to June 2022. Three months after receiving the primary mRNA vaccination or booster dose, the patients underwent a cardiac assessment, including standard echocardiography and speckle-tracking echocardiography (STE). A pre-pandemic historical cohort of age- and gender-matched healthy children were used as a control.

RESULTS

A total of 39 post-VACCINE cases (24, 61% female), mean age 12.6 ± 2.6 years (range 8-17), were enrolled in the study. Ninety percent (N = 35) of patients were previously healthy. No differences in left ventricular diameters, left ventricular ejection fraction (LVEF), and tricuspid annular plane systolic excursion (TAPSE) were observed between cases and controls. Global longitudinal strain (GLS) was in the normal range in all individuals, with no differences between post-VACCINE cases and controls (-21.7 ± 2.3% vs. 21.2 ± 1.8%; p = 0.338). However, GLS was found to be slightly but significantly reduced in post-VACCINE children with a previous COVID-19 compared to naïve-vaccinated individuals (post-VACCINE+COVID-19: -19.9 ± 1.1% vs. post-VACCINE-only: -22.0 ± 2.3%; p = 0.002).

CONCLUSIONS

We did not observe an impairment in GLS or in other indices of LV structure or function after mRNA COVID-19 vaccination.

Severity outcomes associated with SARS-CoV-2 XBB variants, an observational analysis

The rapidity with which SARS-CoV-2 XBB variants rose to predominance has been alarming. We used a large cohort of patients diagnosed with Omicron infections between September 2022 and mid-February 2023 to evaluate the likelihood of admission or need for supplemental oxygen in patients infected with XBB variants. Our data showed no significant association between XBB or XBB.1.5 infections and admissions. Older age groups, lack of vaccination, immunosuppression and underlying heart, kidney, and lung disease showed significant associations with hospitalization.

Towards a post-pandemic future for global pathogen genome sequencing

Pathogen genome sequencing has become a routine part of our response to active outbreaks of infectious disease and should be an important part of our preparations for future epidemics. In this Essay, we discuss the innovations that have enabled routine pathogen genome sequencing, as well as how genome sequences can be used to understand and control the spread of infectious disease. We also explore the impact of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) pandemic on the field of pathogen genomics and outline the challenges we must address to further improve the utility of pathogen genome sequencing in the future.

Epidemic history and evolution of an emerging threat of international concern, the severe acute respiratory syndrome coronavirus 2

This comprehensive review focuses on the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its impact as the cause of the COVID-19 pandemic. Its objective is to provide a cohesive overview of the epidemic history and evolutionary aspects of the virus, with a particular emphasis on its emergence, global spread, and implications for public health. The review delves into the timelines and key milestones of SARS-CoV-2's epidemiological progression, shedding light on the challenges encountered during early containment efforts and subsequent waves of transmission. Understanding the evolutionary dynamics of the virus is crucial in monitoring its potential for adaptation and future outbreaks. Genetic characterization of SARS-CoV-2 is discussed, with a focus on the emergence of new variants and their implications for transmissibility, severity, and immune evasion. The review highlights the important role of genomic surveillance in tracking viral mutations linked to establishing public health interventions. By analyzing the origins, global spread, and genetic evolution of SARS-CoV-2, valuable insights can be gained for the development of effective control measures, improvement of pandemic preparedness, and addressing future emerging infectious diseases of international concern.

Rapid and Technically Simple Detection of SARS-CoV-2 Variants Using CRISPR Cas12 and Cas13

The worldwide proliferation of the SARS-CoV-2 virus in the past 3 years has allowed the virus to accumulate numerous mutations. Dangerous lineages have emerged one after another, each leading to a new wave of the pandemic. In this study, we have developed the THRASOS pipeline to rapidly discover lineage-specific mutation signatures and thus advise the development of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-based diagnostic tests. We also optimized a strategy to modify loop-mediated isothermal amplification amplicons for downstream use with Cas12 and Cas13 for future multiplexing. The close ancestry of the BA.1 and BA.2 variants of SARS-CoV-2 (Omicron) made these excellent candidates for the development of a first test using this workflow. With a quick turnaround time and low requirements for laboratory equipment, the test we have created is ideally suited for settings such as mobile clinics lacking equipment such as Next-Generation Sequencers or Sanger Sequencers and the personnel to run these devices.

Importation of Alpha and Delta variants during the SARS-CoV-2 epidemic in Switzerland: Phylogenetic analysis and intervention scenarios

The SARS-CoV-2 pandemic has led to the emergence of various variants of concern (VoCs) that are associated with increased transmissibility, immune evasion, or differences in disease severity. The emergence of VoCs fueled interest in understanding the potential impact of travel restrictions and surveillance strategies to prevent or delay the early spread of VoCs. We performed phylogenetic analyses and mathematical modeling to study the importation and spread of the VoCs Alpha and Delta in Switzerland in 2020 and 2021. Using a phylogenetic approach, we estimated between 383-1,038 imports of Alpha and 455-1,347 imports of Delta into Switzerland. We then used the results from the phylogenetic analysis to parameterize a dynamic transmission model that accurately described the subsequent spread of Alpha and Delta. We modeled different counterfactual intervention scenarios to quantify the potential impact of border closures and surveillance of travelers on the spread of Alpha and Delta. We found that implementing border closures after the announcement of VoCs would have been of limited impact to mitigate the spread of VoCs. In contrast, increased surveillance of travelers could prove to be an effective measure for delaying the spread of VoCs in situations where their severity remains unclear. Our study shows how phylogenetic analysis in combination with dynamic transmission models can be used to estimate the number of imported SARS-CoV-2 variants and the potential impact of different intervention scenarios to inform the public health response during the pandemic.

VariantHunter: a method and tool for fast detection of emerging SARS-CoV-2 variants

With the progression of the COVID-19 pandemic, large datasets of SARS-CoV-2 genome sequences were collected to closely monitor the evolution of the virus and identify the novel variants/strains. By analyzing genome sequencing data, health authorities can 'hunt' novel emerging variants of SARS-CoV-2 as early as possible, and then monitor their evolution and spread. We designed VariantHunter, a highly flexible and user-friendly tool for systematically monitoring the evolution of SARS-CoV-2 at global and regional levels. In VariantHunter, amino acid changes are analyzed over an interval of 4 weeks in an arbitrary geographical area (continent, country, or region); for every week in the interval, the prevalence is computed and changes are ranked based on their increase or decrease in prevalence. VariantHunter supports two main types of analysis: lineage-independent and lineage-specific. The former considers all the available data and aims to discover new viral variants. The latter evaluates specific lineages/viral variants to identify novel candidate designations (sub-lineages and sub-variants). Both analyses use simple statistics and visual representations (diffusion charts and heatmaps) to track viral evolution. A dataset explorer allows users to visualize available data and refine their selection. VariantHunter is a web application free to all users. The two types of supported analysis (lineage-independent and lineage-specific) allow user-friendly monitoring of the viral evolution, empowering genomic surveillance without requiring any computational background. Database URL http://gmql.eu/variant_hunter/.

CoVigator-A Knowledge Base for Navigating SARS-CoV-2 Genomic Variants

BACKGROUND

The outbreak of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) resulted in the global COVID-19 pandemic. The urgency for an effective SARS-CoV-2 vaccine has led to the development of the first series of vaccines at unprecedented speed. The discovery of SARS-CoV-2 spike-glycoprotein mutants, however, and consequentially the potential to escape vaccine-induced protection and increased infectivity, demonstrates the persisting importance of monitoring SARS-CoV-2 mutations to enable early detection and tracking of genomic variants of concern.

RESULTS

We developed the CoVigator tool with three components: (1) a knowledge base that collects new SARS-CoV-2 genomic data, processes it and stores its results; (2) a comprehensive variant calling pipeline; (3) an interactive dashboard highlighting the most relevant findings. The knowledge base routinely downloads and processes virus genome assemblies or raw sequencing data from the COVID-19 Data Portal (C19DP) and the European Nucleotide Archive (ENA), respectively. The results of variant calling are visualized through the dashboard in the form of tables and customizable graphs, making it a versatile tool for tracking SARS-CoV-2 variants. We put a special emphasis on the identification of intrahost mutations and make available to the community what is, to the best of our knowledge, the largest dataset on SARS-CoV-2 intrahost mutations. In the spirit of open data, all CoVigator results are available for download. The CoVigator dashboard is accessible via covigator.tron-mainz.de.

CONCLUSIONS

With increasing demand worldwide in genome surveillance for tracking the spread of SARS-CoV-2, CoVigator will be a valuable resource of an up-to-date list of mutations, which can be incorporated into global efforts.

Genomic Insights Into the Evolution and Demographic History of the SARS-CoV-2 Omicron Variant: Population Genomics Approach

Background

A thorough understanding of the patterns of genetic subdivision in a pathogen can provide crucial information that is necessary to prevent disease spread. For SARS-CoV-2, the availability of millions of genomes makes this task analytically challenging, and traditional methods for understanding genetic subdivision often fail.

Objective

The aim of our study was to use population genomics methods to identify the subtle subdivisions and demographic history of the Omicron variant, in addition to those captured by the Pango lineage.

Methods

We used a combination of an evolutionary network approach and multivariate statistical protocols to understand the subdivision and spread of the Omicron variant. We identified subdivisions within the BA.1 and BA.2 lineages and further identified the mutations associated with each cluster. We further characterized the overall genomic diversity of the Omicron variant and assessed the selection pressure for each of the genetic clusters identified.

Results

We observed concordant results, using two different methods to understand genetic subdivision. The overall pattern of subdivision in the Omicron variant was in broad agreement with the Pango lineage definition. Further, 1 cluster of the BA.1 lineage and 3 clusters of the BA.2 lineage revealed statistically significant signatures of selection or demographic expansion (Tajima's D<-2), suggesting the role of microevolutionary processes in the spread of the virus.

Conclusions

We provide an easy framework for assessing the genetic structure and demographic history of SARS-CoV-2, which can be particularly useful for understanding the local history of the virus. We identified important mutations that are advantageous to some lineages of Omicron and aid in the transmission of the virus. This is crucial information for policy makers, as preventive measures can be designed to mitigate further spread based on a holistic understanding of the variability of the virus and the evolutionary processes aiding its spread.

Clinical and Epidemiological Profiles of Primary Healthcare Professionals with COVID-19 Infection and Long COVID: An Observational Study

Health professionals have been one of the groups most affected by the SARS-CoV-2 virus. Currently, there is little scientific evidence on the similarities and differences between COVID-19 infection and the development of long COVID in primary care (PC) workers. Therefore, it is necessary to analyse their clinical and epidemiological profiles in depth. This study was observational and descriptive, including PC professionals who were divided into three comparison groups based on the diagnostic test for acute SARS-CoV-2 infection. The responses were analysed using descriptive and bivariate analysis to examinate the relationship between independent variables and the presence or not of long COVID. Binary logistic regression analysis was also conducted, with each symptom as the dependent variable and each group as the independent variable. The results describe the sociodemographic characteristics of these population groups, revealing that women in the health sector are the most affected by long COVID and that being in this group is associated with its development. Furthermore, individuals with long COVID exhibited the highest number of symptoms and pathologies. Certain symptoms were found to be associated with long COVID development in this population, including an altered sense of smell, pneumonia, fever, and sore throat, among others. Similarly, altered senses of smell and taste, chest tightness, and joint pain, among others, were found to be associated with acute COVID-19 infection. Additionally, patients with pre-existing overweight or obesity were more likely to experience acute COVID-19 and develop long COVID. The data obtained can be crucial for improving the detection, diagnosis, and treatment of long COVID patients, ultimately leading to an enhancement in their quality of life.

LAPIS is a fast web API for massive open virus sequencing data

BACKGROUND

Recent epidemic outbreaks such as the SARS-CoV-2 pandemic and the mpox outbreak in 2022 have demonstrated the value of genomic sequencing data for tracking the origin and spread of pathogens. Laboratories around the globe generated new sequences at unprecedented speed and volume and bioinformaticians developed new tools and dashboards to analyze this wealth of data. However, a major challenge that remains is the lack of simple and efficient approaches for accessing and processing sequencing data.

RESULTS

The Lightweight API for Sequences (LAPIS) facilitates rapid retrieval and analysis of genomic sequencing data through a REST API. It supports complex mutation- and metadata-based queries and can perform aggregation operations on massive datasets. LAPIS is optimized for typical questions relevant to genomic epidemiology. Using a newly-developed in-memory database engine, it has a high speed and throughput: between 25 January and 4 February 2023, the SARS-CoV-2 instance of LAPIS, which contains 14.5 million sequences, processed over 20 million requests with a mean response time of 411 ms and a median response time of 1 ms. LAPIS is the core engine behind our dashboards on genspectrum.org and we currently maintain public LAPIS instances for SARS-CoV-2 and mpox.

CONCLUSIONS

Powered by an optimized database engine and available through a web API, LAPIS enhances the accessibility of genomic sequencing data. It is designed to serve as a common backend for dashboards and analyses with the potential to be integrated into common database platforms such as GenBank.

Sequential intrahost evolution and onward transmission of SARS-CoV-2 variants

Persistent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections have been reported in immune-compromised individuals and people undergoing immune-modulatory treatments. Although intrahost evolution has been documented, direct evidence of subsequent transmission and continued stepwise adaptation is lacking. Here we describe sequential persistent SARS-CoV-2 infections in three individuals that led to the emergence, forward transmission, and continued evolution of a new Omicron sublineage, BA.1.23, over an eight-month period. The initially transmitted BA.1.23 variant encoded seven additional amino acid substitutions within the spike protein (E96D, R346T, L455W, K458M, A484V, H681R, A688V), and displayed substantial resistance to neutralization by sera from boosted and/or Omicron BA.1-infected study participants. Subsequent continued BA.1.23 replication resulted in additional substitutions in the spike protein (S254F, N448S, F456L, M458K, F981L, S982L) as well as in five other virus proteins. Our findings demonstrate not only that the Omicron BA.1 lineage can diverge further from its already exceptionally mutated genome but also that patients with persistent infections can transmit these viral variants. Thus, there is, an urgent need to implement strategies to prevent prolonged SARS-CoV-2 replication and to limit the spread of newly emerging, neutralization-resistant variants in vulnerable patients.

Tracking the first SARS-CoV-2 Omicron BA.5.1.3 outbreak in China

The SARS-CoV-2 is still undergoing rapid evolution, resulting in the emergence of several variants of concern, especially the Omicron variants (B.1.1.529), which are surging worldwide. In this study, we tracked Omicron subvariant BA.5.1.3 as the causative agent in the Hainan Province wave in China, which started on 1 August 2022. This was China's first case of Omicron subvariant BA.5.1.3 and led to an indefinite total lockdown in Hainan with more than 8,500 confirmed cases. We obtained 391 whole genomes from positive nasopharyngeal swab samples in the city of Sanya in Hainan Province, which was the center of this outbreak. More than half of the infected cases were female (58%, 227/391) with a median age of 37.0 years (IQR 23.0-53.0). Median Ct values were 24.9 (IQR 22.6-27.3) and 25.2 (IQR 22.9-27.6) for ORF1ab and N genes, respectively. The total single-nucleotide polymorphism (SNP) numbers of Omicron BA.5.1.3 sampled in Sanya (median 69.0, IQR = 69.0-70.0) compared to those worldwide (median 63.0, IQR = 61.0-64.0) showed a significant difference (p < 0.05). Unique core mutations, including three non-synonymous mutations in ORF1ab (Y1064N, S2844G, and R3574K) and one synonymous mutation in ORF3a (S74S), were found. Phylogenetic analysis showed that virus from Sanya formed an independent sub-clade within the BA.5.1.3 subvariant, and could be divided into 15 haplotypes based on the S gene. The most recent common ancestor for the virus from Sanya was estimated as appearing on 5 July 2022, with 95% HPD ranging from 15 May to 20 September 2022. Thanks to our results, we were also able to delineate the mutational profile of this outbreak and highlight the importance of global genomic surveillance and data sharing.

Web Resources for SARS-CoV-2 Genomic Database, Annotation, Analysis and Variant Tracking

The SARS-CoV-2 genomic data continue to grow, providing valuable information for researchers and public health officials. Genomic analysis of these data sheds light on the transmission and evolution of the virus. To aid in SARS-CoV-2 genomic analysis, many web resources have been developed to store, collate, analyze, and visualize the genomic data. This review summarizes web resources used for the SARS-CoV-2 genomic epidemiology, covering data management and sharing, genomic annotation, analysis, and variant tracking. The challenges and further expectations for these web resources are also discussed. Finally, we highlight the importance and need for continued development and improvement of related web resources to effectively track the spread and understand the evolution of the virus.

Varying Cellular Immune Response against SARS-CoV-2 after the Booster Vaccination: A Cohort Study from Fukushima Vaccination Community Survey, Japan

Booster vaccination reduces the incidence of severe cases and mortality related to COVID-19, with cellular immunity playing an important role. However, little is known about the proportion of the population that has achieved cellular immunity after booster vaccination. Thus, we conducted a Fukushima cohort database and assessed humoral and cellular immunity in 2526 residents and healthcare workers in Fukushima Prefecture in Japan through continuous blood collection every 3 months from September 2021. We identified the proportion of people with induced cellular immunity after booster vaccination using the T-SPOT.COVID test, and analyzed their background characteristics. Among 1089 participants, 64.3% (700/1089) had reactive cellular immunity after booster vaccination. Multivariable analysis revealed the following independent predictors of reactive cellular immunity: age < 40 years (adjusted odds ratio: 1.81; 95% confidence interval: 1.19-2.75; p-value: 0.005) and adverse reactions after vaccination (1.92, 1.19-3.09, 0.007). Notably, despite IgG(S) and neutralizing antibody titers of ≥500 AU/mL, 33.9% (349/1031) and 33.5% (341/1017) of participants, respectively, did not have reactive cellular immunity. In summary, this is the first study to evaluate cellular immunity at the population level after booster vaccination using the T-SPOT.COVID test, albeit with several limitations. Future studies will need to evaluate previously infected subjects and their T-cell subsets.

A Web Tool to Estimate Baseline Anti-Spike Monoclonal Antibody Efficacy Based on Regional Genomic Surveillance

Drug appropriateness is a pillar of modern evidence-based medicine, but the turnaround times of genomic sequencing are not compatible with the urgent need to deliver treatments against microorganisms. Massive worldwide genomic surveillance has created an unprecedented landscape for exploiting viral sequencing for therapeutic purposes. When it comes to therapeutic antiviral antibodies, using IC50 against specific polymorphisms of the target antigen can be calculated in vitro, and a list of mutations leading to drug resistance (immune escape) can be compiled. The author encountered this type of knowledge (available from the Stanford University Coronavirus Antiviral Resistance Database,) in a publicly accessible repository of SARS-CoV-2 sequences. The author used a custom function of the CoV-Spectrum.org web portal to deliver up-to-date, regional prevalence estimates of baseline efficacy for each authorized anti-spike mAb across all co-circulating SARS-CoV-2 sublineages at a given time point. This publicly accessible tool can inform therapeutic choices that would otherwise be blind.

Neutralization escape of Omicron XBB, BR.2, and BA.2.3.20 subvariants

New Omicron subvariants continue to emerge throughout the world. In particular, the XBB subvariant, which is a recombinant virus between BA.2.10.1.1 and BA.2.75.3.1.1.1, as well as the BA.2.3.20 and BR.2 subvariants that contain mutations distinct from BA.2 and BA.2.75, are currently increasing in proportion of variants sequenced. Here we show that antibodies induced by 3-dose mRNA booster vaccination as well as BA.1- and BA.4/5-wave infection effectively neutralize BA.2, BR.2, and BA.2.3.20 but have significantly reduced efficiency against XBB. In addition, the BA.2.3.20 subvariant exhibits enhanced infectivity in the lung-derived CaLu-3 cells and in 293T-ACE2 cells. Overall, our results demonstrate that the XBB subvariant is highly neutralization resistant, which highlights the need for continued monitoring of the immune escape and tissue tropism of emerging Omicron subvariants.

Enhanced evasion of neutralizing antibody response by Omicron XBB.1.5, CH.1.1, and CA.3.1 variants

Omicron subvariants continuingly challenge current vaccination strategies. Here, we demonstrate nearly complete escape of the XBB.1.5, CH.1.1, and CA.3.1 variants from neutralizing antibodies stimulated by three doses of mRNA vaccine or by BA.4/5 wave infection, but neutralization is rescued by a BA.5-containing bivalent booster. CH.1.1 and CA.3.1 show strong immune escape from monoclonal antibody S309. Additionally, XBB.1.5, CH.1.1, and CA.3.1 spike proteins exhibit increased fusogenicity and enhanced processing compared with BA.2. Homology modeling reveals the key roles of G252V and F486P in the neutralization resistance of XBB.1.5, with F486P also enhancing receptor binding. Further, K444T/M and L452R in CH.1.1 and CA.3.1 likely drive escape from class II neutralizing antibodies, whereas R346T and G339H mutations could confer the strong neutralization resistance of these two subvariants to S309-like antibodies. Overall, our results support the need for administration of the bivalent mRNA vaccine and continued surveillance of Omicron subvariants.

Future trajectory of SARS-CoV-2: Constant spillover back and forth between humans and animals

SARS-CoV-2, known as severe acute respiratory syndrome coronavirus 2, is causing a massive global public health dilemma. In particular, the outbreak of the Omicron variants of SARS-CoV-2 in several countries has aroused the great attention of the World Health Organization (WHO). As of February 1st, 2023, the WHO had counted 671,016,135 confirmed cases and 6,835,595 deaths worldwide. Despite effective vaccines and drug treatments, there is currently no way to completely and directly eliminate SARS-CoV-2. Moreover, frequent cases of SARS-CoV-2 infection in animals have also been reported. In this review, we suggest that SARS-CoV-2, as a zoonotic virus, may be frequently transmitted between animals and humans in the future, which provides a reference and warning for rational prevention and control of COVID-19.

Zooanthroponotic transmission of SARS-CoV-2 and host-specific viral mutations revealed by genome-wide phylogenetic analysis

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is a generalist virus, infecting and evolving in numerous mammals, including captive and companion animals, free-ranging wildlife, and humans. Transmission among non-human species poses a risk for the establishment of SARS-CoV-2 reservoirs, makes eradication difficult, and provides the virus with opportunities for new evolutionary trajectories, including the selection of adaptive mutations and the emergence of new variant lineages. Here, we use publicly available viral genome sequences and phylogenetic analysis to systematically investigate the transmission of SARS-CoV-2 between human and non-human species and to identify mutations associated with each species. We found the highest frequency of animal-to-human transmission from mink, compared with lower transmission from other sampled species (cat, dog, and deer). Although inferred transmission events could be limited by sampling biases, our results provide a useful baseline for further studies. Using genome-wide association studies, no single nucleotide variants (SNVs) were significantly associated with cats and dogs, potentially due to small sample sizes. However, we identified three SNVs statistically associated with mink and 26 with deer. Of these SNVs, ~⅔ were plausibly introduced into these animal species from local human populations, while the remaining ~⅓ were more likely derived in animal populations and are thus top candidates for experimental studies of species-specific adaptation. Together, our results highlight the importance of studying animal-associated SARS-CoV-2 mutations to assess their potential impact on human and animal health.

Mechanism of a rabbit monoclonal antibody broadly neutralizing SARS-CoV-2 variants

Due to the continuous evolution of SARS-CoV-2, the Omicron variant has emerged and exhibits severe immune evasion. The high number of mutations at key antigenic sites on the spike protein has made a large number of existing antibodies and vaccines ineffective against this variant. Therefore, it is urgent to develop efficient broad-spectrum neutralizing therapeutic drugs. Here we characterize a rabbit monoclonal antibody (RmAb) 1H1 with broad-spectrum neutralizing potency against Omicron sublineages including BA.1, BA.1.1, BA.2, BA.2.12.1, BA.2.75, BA.3 and BA.4/5. Cryo-electron microscopy (cryo-EM) structure determination of the BA.1 spike-1H1 Fab complexes shows that 1H1 targets a highly conserved region of RBD and avoids most of the circulating Omicron mutations, explaining its broad-spectrum neutralization potency. Our findings indicate 1H1 as a promising RmAb model for designing broad-spectrum neutralizing antibodies and shed light on the development of therapeutic agents as well as effective vaccines against newly emerging variants in the future.

Analysis of the ARTIC V4 and V4.1 SARS-CoV-2 primers and their impact on the detection of Omicron BA.1 and BA.2 lineage-defining mutations

The ARTIC protocol uses a multiplexed PCR approach with two primer pools tiling the entire SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) genome. Primer pool updates are necessary for accurate amplicon sequencing of evolving SARS-CoV-2 variants with novel mutations. The suitability of the ARTIC V4 and updated V4.1 primer scheme was assessed using whole genome sequencing of Omicron from clinical samples using Oxford Nanopore Technology. Analysis of Omicron BA.1 genomes revealed that 93.22 % of clinical samples generated improved genome coverage at 50× read depth with V4.1 primers when compared to V4 primers. Additionally, the V4.1 primers improved coverage of BA.1 across amplicons 76 and 88, which resulted in the detection of the variant-defining mutations G22898A, A26530G and C26577G. The Omicron BA.2 sub-variant (VUI-22JAN-01) replaced BA.1 as the dominant variant by March 2022, and analysis of 168 clinical samples showed reduced coverage across amplicons 15 and 75. Upon further interrogation of primer binding sites, a mutation at C4321T [present in 163/168 (97 %) of samples] was identified as a possible cause of complete dropout of amplicon 15. Furthermore, two mutations were identified within the primer binding regions for amplicon 75: A22786C (present in 90 % of samples) and C22792T (present in 12.5 % of samples). Together, these mutations may result in reduced coverage of amplicon 75, and further primer updates would allow the identification of the two BA.2-defining mutations present in amplicon 75: A22688G and T22679C. This work highlights the need for ongoing surveillance of primer matches as circulating variants evolve and change.

Outbreak.info genomic reports: scalable and dynamic surveillance of SARS-CoV-2 variants and mutations

In response to the emergence of SARS-CoV-2 variants of concern, the global scientific community, through unprecedented effort, has sequenced and shared over 11 million genomes through GISAID, as of May 2022. This extraordinarily high sampling rate provides a unique opportunity to track the evolution of the virus in near real-time. Here, we present outbreak.info , a platform that currently tracks over 40 million combinations of Pango lineages and individual mutations, across over 7,000 locations, to provide insights for researchers, public health officials and the general public. We describe the interpretable visualizations available in our web application, the pipelines that enable the scalable ingestion of heterogeneous sources of SARS-CoV-2 variant data and the server infrastructure that enables widespread data dissemination via a high-performance API that can be accessed using an R package. We show how outbreak.info can be used for genomic surveillance and as a hypothesis-generation tool to understand the ongoing pandemic at varying geographic and temporal scales.

A pseudovirus system enables deep mutational scanning of the full SARS-CoV-2 spike

A major challenge in understanding SARS-CoV-2 evolution is interpreting the antigenic and functional effects of emerging mutations in the viral spike protein. Here, we describe a deep mutational scanning platform based on non-replicative pseudotyped lentiviruses that directly quantifies how large numbers of spike mutations impact antibody neutralization and pseudovirus infection. We apply this platform to produce libraries of the Omicron BA.1 and Delta spikes. These libraries each contain ∼7,000 distinct amino acid mutations in the context of up to ∼135,000 unique mutation combinations. We use these libraries to map escape mutations from neutralizing antibodies targeting the receptor-binding domain, N-terminal domain, and S2 subunit of spike. Overall, this work establishes a high-throughput and safe approach to measure how ∼105 combinations of mutations affect antibody neutralization and spike-mediated infection. Notably, the platform described here can be extended to the entry proteins of many other viruses.

Association of COVID-19 with Comorbidities: An Update

Coronavirus disease (COVID-19) is caused by severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) which was identified in Wuhan, China in December 2019 and jeopardized human lives. It spreads at an unprecedented rate worldwide, with serious and still-unfolding health conditions and economic ramifications. Based on the clinical investigations, the severity of COVID-19 appears to be highly variable, ranging from mild to severe infections including the death of an infected individual. To add to this, patients with comorbid conditions such as age or concomitant illnesses are significant predictors of the disease's severity and progression. SARS-CoV-2 enters inside the host cells through ACE2 (angiotensin converting enzyme2) receptor expression; therefore, comorbidities associated with higher ACE2 expression may enhance the virus entry and the severity of COVID-19 infection. It has already been recognized that age-related comorbidities such as Parkinson's disease, cancer, diabetes, and cardiovascular diseases may lead to life-threatening illnesses in COVID-19-infected patients. COVID-19 infection results in the excessive release of cytokines, called "cytokine storm", which causes the worsening of comorbid disease conditions. Different mechanisms of COVID-19 infections leading to intensive care unit (ICU) admissions or deaths have been hypothesized. This review provides insights into the relationship between various comorbidities and COVID-19 infection. We further discuss the potential pathophysiological correlation between COVID-19 disease and comorbidities with the medical interventions for comorbid patients. Toward the end, different therapeutic options have been discussed for COVID-19-infected comorbid patients.

Rise of the BQ.1.1.37 SARS-CoV-2 Sublineage, Italy

BQ.1.1 has dominated the Europe and Americas COVID-19 wave across the 2022-2023 winter, and further viral evolution is expected to escape the consolidating immune responses. We report here the emergence of the BQ.1.1.37 variant in Italy, peaking in January 2022 before suffering competition by XBB.1.*. We attempted to correlate the potential fitness of BQ.1.1.37 with a unique two-amino acid insertion within the Spike protein.

Exhaled Breath Aerosol Shedding of Highly Transmissible Versus Prior Severe Acute Respiratory Syndrome Coronavirus 2 Variants

BACKGROUND

Aerosol inhalation is recognized as the dominant mode of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission. Three highly transmissible lineages evolved during the pandemic. One hypothesis to explain increased transmissibility is that natural selection favors variants with higher rates of viral aerosol shedding. However, the extent of aerosol shedding of successive SARS-CoV-2 variants is unknown. We aimed to measure the infectivity and rate of SARS-CoV-2 shedding into exhaled breath aerosol (EBA) by individuals during the Delta and Omicron waves and compared those rates with those of prior SARS-CoV-2 variants from our previously published work.

METHODS

Individuals with coronavirus disease 2019 (COVID-19) (n = 93; 32 vaccinated and 20 boosted) were recruited to give samples, including 30-minute breath samples into a Gesundheit-II EBA sampler. Samples were quantified for viral RNA using reverse-transcription polymerase chain reaction and cultured for virus.

RESULTS

Alpha (n = 4), Delta (n = 3), and Omicron (n = 29) cases shed significantly more viral RNA copies into EBAs than cases infected with ancestral strains and variants not associated with increased transmissibility (n = 57). All Delta and Omicron cases were fully vaccinated and most Omicron cases were boosted. We cultured virus from the EBA of 1 boosted and 3 fully vaccinated cases.

CONCLUSIONS

Alpha, Delta, and Omicron independently evolved high viral aerosol shedding phenotypes, demonstrating convergent evolution. Vaccinated and boosted cases can shed infectious SARS-CoV-2 via EBA. These findings support a dominant role of infectious aerosols in transmission of SARS-CoV-2. Monitoring aerosol shedding from new variants and emerging pathogens can be an important component of future threat assessments and guide interventions to prevent transmission.

SARS-CoV-2 Omicron Subvariants Balance Host Cell Membrane, Receptor, and Antibody Docking via an Overlapping Target Site

Variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are emerging rapidly and offer surfaces that are optimized for recognition of host cell membranes while also evading antibodies arising from vaccinations and previous infections. Host cell infection is a multi-step process in which spike heads engage lipid bilayers and one or more angiotensin-converting enzyme 2 (ACE-2) receptors. Here, the membrane binding surfaces of Omicron subvariants are compared using cryo-electron microscopy (cEM) structures of spike trimers from BA.2, BA.2.12.1, BA.2.13, BA.2.75, BA.3, BA.4, and BA.5 viruses. Despite significant differences around mutated sites, they all maintain strong membrane binding propensities that first appeared in BA.1. Both their closed and open states retain elevated membrane docking capacities, although the presence of more closed than open states diminishes opportunities to bind receptors while enhancing membrane engagement. The electrostatic dipoles are generally conserved. However, the BA.2.75 spike dipole is compromised, and its ACE-2 affinity is increased, and BA.3 exhibits the opposite pattern. We propose that balancing the functional imperatives of a stable, readily cleavable spike that engages both lipid bilayers and receptors while avoiding host defenses underlies betacoronavirus evolution. This provides predictive criteria for rationalizing future pandemic waves and COVID-19 transmissibility while illuminating critical sites and strategies for simultaneously combating multiple variants.

Early Introduction and Community Transmission of SARS-CoV-2 Omicron Variant, New York, New York, USA

The Omicron variant of SARS-CoV-2 has become dominant in most countries and has raised significant global health concerns. As a global commerce center, New York, New York, USA, constantly faces the risk for multiple variant introductions of SARS-CoV-2. To elucidate the introduction and transmission of the Omicron variant in the city of New York, we created a comprehensive genomic and epidemiologic analysis of 392 Omicron virus specimens collected during November 25-December 11, 2021. We found evidence of 4 independent introductions of Omicron subclades, including the Omicron subclade BA.1.1 with defining substitution of R346K in the spike protein. The continuous genetic divergence within each Omicron subclade revealed their local community transmission and co-circulation in New York, including both household and workplace transmissions supported by epidemiologic evidence. Our study highlights the urgent need for enhanced genomic surveillance and effective response planning for better prevention and management of emerging SARS-CoV-2 variants.

Imprinted SARS-CoV-2 humoral immunity induces convergent Omicron RBD evolution

Continuous evolution of Omicron has led to a rapid and simultaneous emergence of numerous variants that display growth advantages over BA.5 (ref. 1). Despite their divergent evolutionary courses, mutations on their receptor-binding domain (RBD) converge on several hotspots. The driving force and destination of such sudden convergent evolution and its effect on humoral immunity remain unclear. Here we demonstrate that these convergent mutations can cause evasion of neutralizing antibody drugs and convalescent plasma, including those from BA.5 breakthrough infection, while maintaining sufficient ACE2-binding capability. BQ.1.1.10 (BQ.1.1 + Y144del), BA.4.6.3, XBB and CH.1.1 are the most antibody-evasive strains tested. To delineate the origin of the convergent evolution, we determined the escape mutation profiles and neutralization activity of monoclonal antibodies isolated from individuals who had BA.2 and BA.5 breakthrough infections2,3. Owing to humoral immune imprinting, BA.2 and especially BA.5 breakthrough infection reduced the diversity of the neutralizing antibody binding sites and increased proportions of non-neutralizing antibody clones, which, in turn, focused humoral immune pressure and promoted convergent evolution in the RBD. Moreover, we show that the convergent RBD mutations could be accurately inferred by deep mutational scanning profiles4,5, and the evolution trends of BA.2.75 and BA.5 subvariants could be well foreseen through constructed convergent pseudovirus mutants. These results suggest that current herd immunity and BA.5 vaccine boosters may not efficiently prevent the infection of Omicron convergent variants.

Convergent Evolution in SARS-CoV-2 Spike Creates a Variant Soup from Which New COVID-19 Waves Emerge

The first 2 years of the COVID-19 pandemic were mainly characterized by recurrent mutations of SARS-CoV-2 Spike protein at residues K417, L452, E484, N501 and P681 emerging independently across different variants of concern (Alpha, Beta, Gamma, and Delta). Such homoplasy is a marker of convergent evolution. Since Spring 2022 and the third year of the pandemic, with the advent of Omicron and its sublineages, convergent evolution has led to the observation of different lineages acquiring an additional group of mutations at different amino acid residues, namely R346, K444, N450, N460, F486, F490, Q493, and S494. Mutations at these residues have become increasingly prevalent during Summer and Autumn 2022, with combinations showing increased fitness. The most likely reason for this convergence is the selective pressure exerted by previous infection- or vaccine-elicited immunity. Such accelerated evolution has caused failure of all anti-Spike monoclonal antibodies, including bebtelovimab and cilgavimab. While we are learning how fast coronaviruses can mutate and recombine, we should reconsider opportunities for economically sustainable escape-proof combination therapies, and refocus antibody-mediated therapeutic efforts on polyclonal preparations that are less likely to allow for viral immune escape.

Detection of Circulating SARS-CoV-2 Variants of Concern (VOCs) Using a Multiallelic Spectral Genotyping Assay

Throughout the coronavirus disease 2019 (COVID-19) pandemic, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continuously evolved, resulting in new variants, some of which possess increased infectivity, immune evasion, and virulence. Such variants have been denoted by the World Health Organization as variants of concern (VOC) because they have resulted in an increased number of cases, posing a strong risk to public health. Thus far, five VOCs have been designated, Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2), and Omicron (B.1.1.529), including their sublineages. Next-generation sequencing (NGS) can produce a significant amount of information facilitating the study of variants; however, NGS is time-consuming and costly and not efficient during outbreaks, when rapid identification of VOCs is urgently needed. In such periods, there is a need for fast and accurate methods, such as real-time reverse transcription PCR in combination with probes, which can be used for monitoring and screening of the population for these variants. Thus, we developed a molecular beacon-based real-time RT-PCR assay according to the principles of spectral genotyping. This assay employs five molecular beacons that target ORF1a:ΔS3675/G3676/F3677, S:ΔH69/V70, S:ΔE156/F157, S:ΔΝ211, S:ins214EPE, and S:ΔL242/A243/L244, deletions and an insertion found in SARS-CoV-2 VOCs. This assay targets deletions/insertions because they inherently provide higher discrimination capacity. Here, the design process of the molecular beacon-based real-time RT-PCR assay for detection and discrimination of SARS-CoV-2 is presented, and experimental testing using SARS-CoV-2 VOC samples from reference strains (cultured virus) and clinical patient samples (nasopharyngeal samples), which have been previously classified using NGS, were evaluated. Based on the results, it was shown that all molecular beacons can be used under the same real-time RT-PCR conditions, consequently improving the time and cost efficiency of the assay. Furthermore, this assay was able to confirm the genotype of each of the tested samples from various VOCs, thereby constituting an accurate and reliable method for VOC detection and discrimination. Overall, this assay is a valuable tool that can be used for screening and monitoring the population for VOCs or other emerging variants, contributing to limiting their spread and protecting public health.

Extraordinary Evasion of Neutralizing Antibody Response by Omicron XBB.1.5, CH.1.1 and CA.3.1 Variants

Newly emerging Omicron subvariants continue to emerge around the world, presenting potential challenges to current vaccination strategies. This study investigates the extent of neutralizing antibody escape by new subvariants XBB.1.5, CH.1.1, and CA.3.1, as well as their impacts on spike protein biology. Our results demonstrated a nearly complete escape of these variants from neutralizing antibodies stimulated by three doses of mRNA vaccine, but neutralization was rescued by a bivalent booster. However, CH.1.1 and CA.3.1 variants were highly resistant to both monovalent and bivalent mRNA vaccinations. We also assessed neutralization by sera from individuals infected during the BA.4/5 wave of infection and observed similar trends of immune escape. In these cohorts, XBB.1.5 did not exhibit enhanced neutralization resistance over the recently dominant BQ.1.1 variant. Notably, the spike proteins of XBB.1.5, CH.1.1, and CA.3.1 all exhibited increased fusogenicity compared to BA.2, correlating with enhanced S processing. Overall, our results support the administration of new bivalent mRNA vaccines, especially in fighting against newly emerged Omicron subvariants, as well as the need for continued surveillance of Omicron subvariants.

Detection and Molecular Characterization of the SARS-CoV-2 Delta Variant and the Specific Immune Response in Companion Animals in Switzerland

In human beings, there are five reported variants of concern of severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2). However, in contrast to human beings, descriptions of infections of animals with specific variants are still rare. The aim of this study is to systematically investigate SARS-CoV-2 infections in companion animals in close contact with SARS-CoV-2-positive owners ("COVID-19 households") with a focus on the Delta variant. Samples, obtained from companion animals and their owners were analyzed using a real-time reverse transcriptase-polymerase chain reaction (RT-qPCR) and next-generation sequencing (NGS). Animals were also tested for antibodies and neutralizing activity against SARS-CoV-2. Eleven cats and three dogs in nine COVID-19-positive households were RT-qPCR and/or serologically positive for the SARS-CoV-2 Delta variant. For seven animals, the genetic sequence could be determined. The animals were infected by one of the pangolin lineages B.1.617.2, AY.4, AY.43 and AY.129 and between zero and three single-nucleotide polymorphisms (SNPs) were detected between the viral genomes of animals and their owners, indicating within-household transmission between animal and owner and in multi-pet households also between the animals. NGS data identified SNPs that occur at a higher frequency in the viral sequences of companion animals than in viral sequences of humans, as well as SNPs, which were exclusively found in the animals investigated in the current study and not in their owners. In conclusion, our study is the first to describe the SARS-CoV-2 Delta variant transmission to animals in Switzerland and provides the first-ever description of Delta-variant pangolin lineages AY.129 and AY.4 in animals. Our results reinforce the need of a One Health approach in the monitoring of SARS-CoV-2 in animals.

Longitudinal neutralization activities on authentic Omicron variant provided by three doses of BBIBP-CorV vaccination during one year

The majority of people in China have been immunized with the inactivated viral vaccine BBIBP-CorV. The emergence of the Omicron variant raised the concerns about protection efficacy of the inactivated viral vaccine in China. However, longitudinal neutralization data describing protection efficacy against Omicron variant is still lacking. Here we present one-year longitudinal neutralization data of BBIBP-CorV on authentic Omicron, Delta, and wild-type strains using 224 sera collected from 14 volunteers who have finished three doses BBIBP-CorV. The sera were also subjected for monitoring the SARS-CoV-2 specific IgG, IgA, and IgM responses on protein and peptide microarrays. The neutralization titers showed different protection efficacies against the three strains. By incorporating IgG and IgA signals of proteins and Spike protein derived peptide on microarray, panels as potential surrogate biomarkers for rapid estimation of neutralization titers were established. These data support the necessity of the 3^rd dose of BBIBP-CorV vaccination. After further validation and assay development, the panels could be used for reliable, convenient and fast evaluation of the efficacy of vaccination.

Genomic Epidemiology of the SARS-CoV-2 Epidemic in Cyprus from November 2020 to October 2021: The Passage of Waves of Alpha and Delta Variants of Concern

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in December 2019 resulted in the coronavirus disease 2019 (COVID-19) pandemic, which has had devastating repercussions for public health. Over the course of this pandemic, the virus has continuously been evolving, resulting in new, more infectious variants that have frequently led to surges of new SARS-CoV-2 infections. In the present study, we performed detailed genetic, phylogenetic, phylodynamic and phylogeographic analyses to examine the SARS-CoV-2 epidemic in Cyprus using 2352 SARS-CoV-2 sequences from infected individuals in Cyprus during November 2020 to October 2021. During this period, a total of 61 different lineages and sublineages were identified, with most falling into three groups: B.1.258 & sublineages, Alpha (B.1.1.7 & Q. sublineages), and Delta (B.1.617.2 & AY. sublineages), each encompassing a set of S gene mutations that primarily confer increased transmissibility as well as immune evasion. Specifically, these lineages were coupled with surges of new infections in Cyprus, resulting in the following: the second wave of SARS-CoV-2 infections in Cyprus, comprising B.1.258 & sublineages, during late autumn 2020/beginning of winter 2021; the third wave, comprising Alpha (B.1.1.7 & Q. sublineages), during spring 2021; and the fourth wave, comprising Delta (B.1.617.2 & AY. sublineages) during summer 2021. Additionally, it was identified that these lineages were primarily imported from and exported to the UK, Greece, and Sweden; many other migration links were also identified, including Switzerland, Denmark, Russia, and Germany. Taken together, the results of this study indicate that the SARS-CoV-2 epidemic in Cyprus was characterized by successive introduction of new lineages from a plethora of countries, resulting in the generation of waves of infection. Overall, this study highlights the importance of investigating the spatiotemporal evolution of the SARS-CoV-2 epidemic in the context of Cyprus, as well as the impact of protective measures placed to mitigate transmission of the virus, providing necessary information to safeguard public health.

Neutralization sensitivity, fusogenicity, and infectivity of Omicron subvariants

BACKGROUND

The emergence of SARS-CoV-2 Omicron subvariants has raised questions regarding resistance to immunity by natural infection or immunization. We examined the sensitivity of Delta and Omicron subvariants (BA.1, BA.1.1, BA.2, BA.2.12.1, BA.4/5, and BA.3) to neutralizing antibodies from BBIBP-CorV-vaccinated and BBIBP-CorV- or ZF2001-boosted individuals, as well as individuals with Delta and BA.1 breakthrough infections, and determined their fusogenicity and infectivity.

METHODS

In this cross-sectional study, serum samples from two doses of BBIBP-CorV-vaccinated individuals 1 (n = 36), 3 (n = 36), and 7 (n = 37) months after the second dose; BBIBP-CorV- (n = 25) or ZF2001-boosted (n = 30) individuals; and fully vaccinated individuals with Delta (n = 30) or BA.1 (n = 26) infection were collected. The serum-neutralizing reactivity and potency of bebtelovimab were assessed against D614G, Delta, and Omicron subvariants (BA.1, BA.1.1, BA.2, BA.2.12.1, BA.4/5, and BA.3) through a pseudovirus neutralization assay. The fusogenicity and infectivity of D614G, Delta, and Omicron subvariants were determined by cell-cell fusion assay and pseudovirus infection assay, respectively.

RESULTS

Omicron subvariants markedly escaped vaccine-elicited neutralizing antibodies after two doses of BBIBP-CorV with comparable efficiency. A third dose vaccination of BBIBP-CorV or ZF2001 increased neutralizing antibody titers and breadth against Delta and three Omicron subvariants. Delta and BA.1 breakthrough infections induced comparable neutralizing antibody titers against D614G and Delta variants, whereas BA.1 breakthrough infections elicited a stronger and broader antibody response against three Omicron subvariants than Delta breakthrough infections. BA.2.12.1 and BA.4/5 are more resistant to immunity induced by breakthrough infections. Bebtelovimab had no significant loss of potency against the Delta and Omicron subvariants. Cell culture experiments showed Omicron subvariants to be less fusogenic and have higher infectivity than D614G and Delta with comparable efficiency.

CONCLUSIONS

These findings have important public health implications and highlight the importance of repeated exposure to SARS-CoV-2 antigens to broaden the neutralizing antibody response against Omicron subvariants.