VarEPS: an evaluation and prewarning system of known and virtual variations of SARS-CoV-2 genomes

Tailored wastewater surveillance framework uncovered the epidemics of key pathogens in a Northwestern city of China

Wastewater surveillance enables rapid pathogen monitoring and community prevalence estimation. However, how to design an integrated and tailored wastewater surveillance framework to monitor major health threats in metropolises remains a major challenge. In this study, we first analyzed the historical clinical data of Xi'an city and designed a wastewater surveillance framework covering five key endemic viruses, namely, SARS-CoV-2, norovirus, influenza A virus (IAV), influenza B virus (IBV), respiratory syncytial virus (RSV), and hantavirus. Amplicon sequencing of SARS-CoV-2, norovirus and hantavirus was conducted biweekly to determine the prevalent community genotypes circulating in this region. The results showed that from April 2023 to August 2023, Xi'an experienced two waves of SARS-CoV-2 infection, which peaked in the middle of May-2023 and late August-2023. The sewage concentrations of IAV and RSV peaked in early March and early May 2023, respectively, while the sewage concentrations of norovirus fluctuated throughout the study period and peaked in late August. The dynamics of the sewage concentrations of SARS-CoV-2, norovirus, IAV, RSV, and hantavirus were in line with the trends in the sentinel hospital percent positivity data, indicating the role of wastewater surveillance in enhancing the understanding of epidemic trends. Amplicon sequencing of SARS-CoV-2 revealed a transition in the predominant genotype, which changed from DY.1 and FR.1.4 to the XBB and EG.5 subvariants. Amplicon sequencing also revealed that there was only one predominant hantavirus genotype in the local population, while highly diverse genotypes of norovirus GI and GII were found in the wastewater. In conclusion, this study provided valuable insights into the dynamics of infection trends and predominant genotypes of key pathogens in a city without sufficient clinical surveillance, highlighting the role of a tailored wastewater surveillance framework in addressing public health priorities. More importantly, our study provides the first evidence demonstrating the applicability of wastewater surveillance for hantavirus, which is a major health threat locally.

CovEpiAb: a comprehensive database and analysis resource for immune epitopes and antibodies of human coronaviruses

Coronaviruses have threatened humans repeatedly, especially COVID-19 caused by SARS-CoV-2, which has posed a substantial threat to global public health. SARS-CoV-2 continuously evolves through random mutation, resulting in a significant decrease in the efficacy of existing vaccines and neutralizing antibody drugs. It is critical to assess immune escape caused by viral mutations and develop broad-spectrum vaccines and neutralizing antibodies targeting conserved epitopes. Thus, we constructed CovEpiAb, a comprehensive database and analysis resource of human coronavirus (HCoVs) immune epitopes and antibodies. CovEpiAb contains information on over 60 000 experimentally validated epitopes and over 12 000 antibodies for HCoVs and SARS-CoV-2 variants. The database is unique in (1) classifying and annotating cross-reactive epitopes from different viruses and variants; (2) providing molecular and experimental interaction profiles of antibodies, including structure-based binding sites and around 70 000 data on binding affinity and neutralizing activity; (3) providing virological characteristics of current and past circulating SARS-CoV-2 variants and in vitro activity of various therapeutics; and (4) offering site-level annotations of key functional features, including antibody binding, immunological epitopes, SARS-CoV-2 mutations and conservation across HCoVs. In addition, we developed an integrated pipeline for epitope prediction named COVEP, which is available from the webpage of CovEpiAb. CovEpiAb is freely accessible at https://pgx.zju.edu.cn/covepiab/.

VarEPS-Influ:an risk evaluation system of occurred and virtual variations of influenza virus genomes

Influenza viruses undergo frequent genomic mutations, leading to potential cross-species transmission, phenotypic changes, and challenges in diagnostic reagents and vaccines. Accurately evaluating and predicting the risk of such variations remain significant challenges. To address this, we developed the VarEPS-Influ database, an influenza virus variations risk evaluation system (VarEPS-Influ). This database employs a 'multi-dimensional evaluation of mutations' strategy, utilizing various tools to assess the physical and chemical properties, primary, secondary, and tertiary structures, receptor affinity, antibody binding capacity, antigen epitopes, and other aspects of the variation's impact. Additionally, we consider space-time distribution, host species distribution, pedigree analysis, and frequency of mutations to provide a comprehensive risk evaluation of mutations and viruses. The VarEPS-Influ database evaluates both observed variations and virtual variations (variations that have not yet occurred), thereby addressing the time-lag issue in risk predictions. Our current one-stop evaluation system for influenza virus genomic variation integrates 1065290 sequences from 224 927 Influenza A, B and C isolates retrieved from public resources. Researchers can freely access the data at https://nmdc.cn/influvar/.

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/.

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.

A tool to automatically design multiplex PCR primer pairs for specific targets using diverse templates

Multiplex PCR is an increasingly popular method for identifying species, investigating environmental diversity, and conducting phylogenetic analysis. The complexity and increasing availability of diverse templates necessitate a highly automated approach to design degenerate primer pairs for specific targets with multiple sequences. Existing tools for degenerate primer design suffer from poor maintenance, semi-automation, low adaptability, and low tolerance for gaps. We developed PMPrimer, a Python-based tool for automated design and evaluation of multiplex PCR primer pairs for specific targets using diverse templates. PMPrimer automatically designs optimal multiplex PCR primer pairs using a statistical-based template filter; performs multiple sequence alignment, conserved region identification, and primer design; and evaluates the primers based on template coverage, taxon specificity, and target specificity. PMPrimer identifies conserved regions using Shannon's entropy method, tolerates gaps using a haplotype-based method, and evaluates multiplex PCR primer pairs based on template coverage and taxon specificity. We tested PMPrimer using datasets with diverse levels of conservation, sizes, and applications, including tuf genes of Staphylococci, hsp65 genes of Mycobacteriaceae, and 16S ribosomal RNA genes of Archaea. PMPrimer showed outstanding performance compared with existing tools and experimental validated primers. PMPrimer is available as a Python package at https://github.com/AGIScuipeng/PMPrimer .

The role of SARS-CoV-2 N protein in diagnosis and vaccination in the context of emerging variants: present status and prospects

Despite many countries rapidly revising their strategies to prevent contagions, the number of people infected with Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to surge. The emergent variants that can evade the immune response significantly affect the effectiveness of mainstream vaccines and diagnostic products based on the original spike protein. Therefore, it is essential to focus on the highly conserved nature of the nucleocapsid protein as a potential target in the field of vaccines and diagnostics. In this regard, our review initially discusses the structure, function, and mechanism of action of N protein. Based on this discussion, we summarize the relevant research on the in-depth development and application of diagnostic methods and vaccines based on N protein, such as serology and nucleic acid detection. Such valuable information can aid in designing more efficient diagnostic and vaccine tools that could help end the SARS-CoV-2 pandemic.

Early warning of emerging infectious diseases based on multimodal data

The coronavirus disease 2019 (COVID-19) pandemic has dramatically increased the awareness of emerging infectious diseases. The advancement of multiomics analysis technology has resulted in the development of several databases containing virus information. Several scientists have integrated existing data on viruses to construct phylogenetic trees and predict virus mutation and transmission in different ways, providing prospective technical support for epidemic prevention and control. This review summarized the databases of known emerging infectious viruses and techniques focusing on virus variant forecasting and early warning. It focuses on the multi-dimensional information integration and database construction of emerging infectious viruses, virus mutation spectrum construction and variant forecast model, analysis of the affinity between mutation antigen and the receptor, propagation model of virus dynamic evolution, and monitoring and early warning for variants. As people have suffered from COVID-19 and repeated flu outbreaks, we focused on the research results of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and influenza viruses. This review comprehensively viewed the latest virus research and provided a reference for future virus prevention and control research.

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.

Interpretable and Predictive Deep Neural Network Modeling of the SARS-CoV-2 Spike Protein Sequence to Predict COVID-19 Disease Severity

Through the COVID-19 pandemic, SARS-CoV-2 has gained and lost multiple mutations in novel or unexpected combinations. Predicting how complex mutations affect COVID-19 disease severity is critical in planning public health responses as the virus continues to evolve. This paper presents a novel computational framework to complement conventional lineage classification and applies it to predict the severe disease potential of viral genetic variation. The transformer-based neural network model architecture has additional layers that provide sample embeddings and sequence-wide attention for interpretation and visualization. First, training a model to predict SARS-CoV-2 taxonomy validates the architecture's interpretability. Second, an interpretable predictive model of disease severity is trained on spike protein sequence and patient metadata from GISAID. Confounding effects of changing patient demographics, increasing vaccination rates, and improving treatment over time are addressed by including demographics and case date as independent input to the neural network model. The resulting model can be interpreted to identify potentially significant virus mutations and proves to be a robust predctive tool. Although trained on sequence data obtained entirely before the availability of empirical data for Omicron, the model can predict the Omicron's reduced risk of severe disease, in accord with epidemiological and experimental data.

The rapid emergence of multiple sublineages of Omicron (B.1.1.529) variant: Dynamic profiling via molecular phylogenetics and mutational landscape studies

PURPOSE

The recent Omicron (B.1.1.529) variant poses a significant threat to global health. This variant has spread worldwide, and several sublineages have rapidly emerged. Study tried to analyze the microevolution of this variant.

METHODS

We studied the molecular phylogenetics, divergence, geographical distributions, frequencies, risk mutations for antibody affinity, and mutational landscape for Omicron sublineages using in silico analysis and statistical models. The risk mutation of spike for nAb affinity was analyzed and illustrated by statistical plots. Finally, the mutational properties of the spike mutations and their stability were predicted and demonstrated.

RESULTS

First, we studied the microevolutionary Omicron sublineages using molecular phylogenetics. Simultaneously, we revealed divergence events of the Omicron sublineages and observed the lowest minimum divergence of 51 in clade 21K and the highest maximum divergence of 90 in clade 21L. We have demonstrated cluster analyses, geographical distributions, frequencies of Omicron and its sublineages. Finally, we evaluated the mutational landscape of the Omicron sublineages. In this mutational study, we performed a genome-wide analysis of general mutations, mutations in the non-spike genome, and spike mutations of Omicron sublineages. The risk mutation of S-glycoprotein for nAb affinity has been analyzed for Omicron sublineages. Here, we found that some sublineages have all four significant highly destabilizing mutations. Such sublineages are BA.1 (G446S, E484A, T95I, and D614G), BA.2 (H655Y, Q493R, G493S, and D614G), BA.4 (N501Y, Y505H, N969K, and D614G), and BA.2.75 (Q454H, T547K, N764K, D614G and G446S). Finally, from the mutation stability prediction through ΔΔG, we observed that BA.1 and BA.4 had two destabilizing and two stabilizing mutations. Similarly, BA.2, BA.5, and BA.2.12.1 have one destabilizing and three stabilizing mutations. However, all four mutations in BA.2.75 are stabilizing mutations.

CONCLUSIONS

Our molecular phylogenetic studies provided a deeper understanding of the microevolution of sublineages and the creation of Omicron. Similarly, this study might help scientists develop pan-coronavirus vaccines that consider their mutational properties.

A comprehensive analysis of the mutational landscape of the newly emerging Omicron (B.1.1.529) variant and comparison of mutations with VOCs and VOIs

The Omicron variant is spreading rapidly throughout several countries. Thus, we comprehensively analyzed Omicron's mutational landscape and compared mutations with VOC/VOI. We analyzed SNVs throughout the genome, and AA variants (NSP and SP) in VOC/VOI, including Omicron. We generated heat maps to illustrate the AA variants with high mutation prevalence (> 75% frequency) of Omicron, which demonstrated eight mutations with > 90% prevalence in ORF1a and 29 mutations with > 75% prevalence in S-glycoprotein. A scatter plot for Omicron and VOC/VOI's cluster evaluation was computed. We performed a risk analysis of the antibody-binding risk among four mutations (L452, F490, P681, D614) and observed three mutations (L452R, F490S, D614G) destabilized antibody interactions. Our comparative study evaluated the properties of 28 emerging mutations of the S-glycoprotein of Omicron, and the ΔΔG values. Our results showed K417N with minimum and Q954H with maximum ΔΔG value. Furthermore, six important RBD mutations (G339D, S371L, N440K, G446S, T478K, Q498R) were chosen for comprehensive analysis for stabilizing/destabilizing properties and molecular flexibility. The G339D, S371L, N440K, and T478K were noted as stable mutations with 0.019 kcal/mol, 0.127 kcal/mol, 0.064 kcal/mol, and 1.009 kcal/mol. While, G446S and Q498R mutations showed destabilizing results. Simultaneously, among six RBD mutations, G339D, G446S, and Q498R mutations increased the molecular flexibility of S-glycoprotein. This study depicts the comparative mutational pattern of Omicron and other VOC/VOI, which will help researchers to design and deploy novel vaccines and therapeutic antibodies to fight against VOC/VOI, including Omicron.

covSampler: A Subsampling Method with Balanced Genetic Diversity for Large-Scale SARS-CoV-2 Genome Data Sets

Phylogenetic analysis has been widely used to describe, display and infer the evolutionary patterns of viruses. The unprecedented accumulation of SARS-CoV-2 genomes has provided valuable materials for the real-time study of SARS-CoV-2 evolution. However, the large number of SARS-CoV-2 genome sequences also poses great challenges for data analysis. Several methods for subsampling these large data sets have been introduced. However, current methods mainly focus on the spatiotemporal distribution of genomes without considering their genetic diversity, which might lead to postsubsampling bias. In this study, a subsampling method named covSampler was developed for the subsampling of SARS-CoV-2 genomes with consideration of both their spatiotemporal distribution and their genetic diversity. First, covSampler clusters all genomes according to their spatiotemporal distribution and genetic variation into groups that we call divergent pathways. Then, based on these divergent pathways, two kinds of subsampling strategies, representative subsampling and comprehensive subsampling, were provided with adjustable parameters to meet different users’ requirements. Our performance and validation tests indicate that covSampler is efficient and stable, with an abundance of options for user customization. Overall, our work has developed an easy-to-use tool and a webserver (https://www.covsampler.net) for the subsampling of SARS-CoV-2 genome sequences.

Omicron variant (B.1.1.529) of SARS-CoV-2: understanding mutations in the genome, S-glycoprotein, and antibody-binding regions

The Omicron variant has been detected in nearly 150 countries. We analyzed the mutational landscape of Omicron throughout the genome, focusing the S-glycoprotein. We also evaluated mutations in the antibody-binding regions and observed some important mutations overlapping those of previous variants including N501Y, D614G, H655Y, N679K, and P681H. Various new receptor-binding domain mutations were detected, including Q493K, G496S, Q498R, S477N, G466S, N440K, and Y505H. New mutations were found in the NTD (Δ143-145, A67V, T95I, L212I, and Δ211) including one new mutation in fusion peptide (D796Y). There are several mutations in the antibody-binding region including K417N, E484A, Q493K, Q498R, N501Y, and Y505H and several near the antibody-binding region (S477N, T478K, G496S, G446S, and N440K). The impact of mutations in regions important for the affinity between spike proteins and neutralizing antibodies was evaluated. Furthermore, we examined the effect of significant antibody-binding mutations (K417N, T478K, E484A, and N501Y) on antibody affinity, stability to ACE2 interaction, and possibility of amino acid substitution. All the four mutations destabilize the antibody-binding affinity. This study reveals future directions for developing neutralizing antibodies against the Omicron variant.

The 2022 Nucleic Acids Research database issue and the online molecular biology database collection

The 2022 Nucleic Acids Research Database Issue contains 185 papers, including 87 papers reporting on new databases and 85 updates from resources previously published in the Issue. Thirteen additional manuscripts provide updates on databases most recently published elsewhere. Seven new databases focus specifically on COVID-19 and SARS-CoV-2, including SCoV2-MD, the first of the Issue's Breakthrough Articles. Major nucleic acid databases reporting updates include MODOMICS, JASPAR and miRTarBase. The AlphaFold Protein Structure Database, described in the second Breakthrough Article, is the stand-out in the protein section, where the Human Proteoform Atlas and GproteinDb are other notable new arrivals. Updates from DisProt, FuzDB and ELM comprehensively cover disordered proteins. Under the metabolism and signalling section Reactome, ConsensusPathDB, HMDB and CAZy are major returning resources. In microbial and viral genomes taxonomy and systematics are well covered by LPSN, TYGS and GTDB. Genomics resources include Ensembl, Ensembl Genomes and UCSC Genome Browser. Major returning pharmacology resource names include the IUPHAR/BPS guide and the Therapeutic Target Database. New plant databases include PlantGSAD for gene lists and qPTMplants for post-translational modifications. The entire Database Issue is freely available online on the Nucleic Acids Research website (https://academic.oup.com/nar). Our latest update to the NAR online Molecular Biology Database Collection brings the total number of entries to 1645. Following last year's major cleanup, we have updated 317 entries, listing 89 new resources and trimming 80 discontinued URLs. The current release is available at http://www.oxfordjournals.org/nar/database/c/.