Geographic and Genomic Distribution of SARS-CoV-2 Mutations

Dynamic of SARS-CoV-2 variants circulation in Tunisian pediatric population, during successive waves, from March 2020 to September 2022

The emergence of SARS-CoV-2 variants has led to several cases among children. However, limited information is available from North African countries. This study describes the SARS-CoV-2 strains circulating in Tunisian pediatric population during successive waves. A total of 447 complete sequences were obtained from individuals aged from 13 days to 18 years, between March 2020 and September 2022: 369 sequences generated during this study and 78 ones, available in GISAID, previously obtained from Tunisian pediatric patients. These sequences were compared with 354 and 274 ones obtained from Tunisian adults and a global dataset, respectively. The variant circulation dynamics of predominant variants were investigated during the study period using maximum-likelihood phylogenetic analysis. Among the studied population, adolescents were the predominant age group, comprising 55.26% of cases. Twenty-three lineages were identified; seven of which were not previously reported in Tunisia. Phylogenetic analysis showed a close relationship between the sequences from Tunisian adults and children. The connections of sequences from other countries were variable according to variants: close relationships were observed for Alpha, B1.160 and Omicron variants, while independent Tunisian clusters were observed for Delta and B.1.177 lineages. These findings highlight the pivotal role of children in virus transmission and underscore the impact of vaccination on virus spread. Vaccination of children, with booster doses, may be considered for better management of future emergences.

Mapping of long stretches of highly conserved sequences in over 6 million SARS-CoV-2 genomes

We identified 11 conserved stretches in over 6.3 million SARS-CoV-2 genomes including all the major variants of concerns. Each conserved stretch is ≥100 nucleotides in length with ≥99.9% conservation at each nucleotide position. Interestingly, six of the eight conserved stretches in ORF1ab overlapped significantly with well-folded experimentally verified RNA secondary structures. Furthermore, two of the conserved stretches were mapped to regions within the S2-subunit that undergo dynamic structural rearrangements during viral fusion. In addition, the conserved stretches were significantly depleted for zinc-finger antiviral protein (ZAP) binding sites, which facilitated the recognition and degradation of viral RNA. These highly conserved stretches in the SARS-CoV-2 genome were poorly conserved at the nucleotide level among closely related β-coronaviruses, thus representing ideal targets for highly specific and discriminatory diagnostic assays. Our findings highlight the role of structural constraints at both RNA and protein levels that contribute to the sequence conservation of specific genomic regions in SARS-CoV-2.

A measles-vectored vaccine candidate expressing prefusion-stabilized SARS-CoV-2 spike protein brought to phase I/II clinical trials: protection of African green monkeys from COVID-19 disease

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged at the end of 2019 and is responsible for the largest human pandemic in 100 years. Thirty-four vaccines are currently approved for use worldwide, and approximately 67% of the world population has received a complete primary series of one, yet countries are dealing with new waves of infections, variant viruses continue to emerge, and breakthrough infections are frequent secondary to waning immunity. Here, we evaluate a measles virus (MV)-vectored vaccine expressing a stabilized prefusion SARS-CoV-2 spike (S) protein (MV-ATU3-S2PΔF2A; V591) with demonstrated immunogenicity in mouse models (see companion article [J. Brunet, Z. Choucha, M. Gransagne, H. Tabbal, M.-W. Ku et al., J Virol 98:e01693-23, 2024, https://doi.org/10.1128/jvi.01693-23]) in an established African green monkey model of disease. Animals were vaccinated with V591 or the control vaccine (an equivalent MV-vectored vaccine with an irrelevant antigen) intramuscularly using a prime/boost schedule, followed by challenge with an early pandemic isolate of SARS-CoV-2 at 56 days post-vaccination. Pre-challenge, only V591-vaccinated animals developed S-specific antibodies that had virus-neutralizing activity as well as S-specific T cells. Following the challenge, V591-vaccinated animals had lower infectious virus and viral (v) RNA loads in mucosal secretions and stopped shedding virus in these secretions earlier. vRNA loads were lower in these animals in respiratory and gastrointestinal tract tissues at necropsy. This correlated with a lower disease burden in the lungs as quantified by PET/CT at early and late time points post-challenge and by pathological analysis at necropsy.IMPORTANCESevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the largest human pandemic in 100 years. Even though vaccines are currently available, countries are dealing with new waves of infections, variant viruses continue to emerge, breakthrough infections are frequent, and vaccine hesitancy persists. This study uses a safe and effective measles vaccine as a platform for vaccination against SARS-CoV-2. The candidate vaccine was used to vaccinate African green monkeys (AGMs). All vaccinated AGMs developed robust antigen-specific immune responses. After challenge, these AGMs produced less virus in mucosal secretions, for a shorter period, and had a reduced disease burden in the lungs compared to control animals. At necropsy, lower levels of viral RNA were detected in tissue samples from vaccinated animals, and the lungs of these animals lacked the histologic hallmarks of SARS-CoV-2 disease observed exclusively in the control AGMs.

Intrinsic RNA Targeting Triggers Indiscriminate DNase Activity of CRISPR-Cas12a

The CRISPR-Cas12a system has emerged as a powerful tool for next-generation nucleic acid-based molecular diagnostics. However, it has long been believed to be effective only on DNA targets. Here, we investigate the intrinsic RNA-enabled trans-cleavage activity of AsCas12a and LbCas12a and discover that they can be directly activated by full-size RNA targets, although LbCas12a exhibits weaker trans-cleavage activity than AsCas12a on both single-stranded DNA and RNA substrates. Remarkably, we find that the RNA-activated Cas12a possesses higher specificity in recognizing mutated target sequences compared to DNA activation. Based on these findings, we develop the "Universal Nuclease for Identification of Virus Empowered by RNA-Sensing" (UNIVERSE) assay for nucleic acid testing. We incorporate a T7 transcription step into this assay, thereby eliminating the requirement for a protospacer adjacent motif (PAM) sequence in the target. Additionally, we successfully detect multiple PAM-less targets in HIV clinical samples that are undetectable by the conventional Cas12a assay based on double-stranded DNA activation, demonstrating unrestricted target selection with the UNIVERSE assay. We further validate the clinical utility of the UNIVERSE assay by testing both HIV RNA and HPV 16 DNA in clinical samples. We envision that the intrinsic RNA targeting capability may bring a paradigm shift in Cas12a-based nucleic acid detection and further enhance the understanding of CRISPR-Cas biochemistry.

Wastewater sequencing as a powerful tool to reveal SARS-CoV-2 variant introduction and spread in French Guiana, South America

The origin of introduction of a new pathogen in a country, the evolutionary dynamics of an epidemic within a country, and the role of cross-border areas on pathogen dynamics remain complex to disentangle and are often poorly understood. For instance, cross-border areas represent the ideal location for the sharing of viral variants between countries, with international air travel, land travel and waterways playing an important role in the cross-border spread of infectious diseases. Unfortunately, monitoring the point of entry and the evolutionary dynamics of viruses in space and time within local populations remain challenging. Here we tested the efficiency of wastewater-based epidemiology and genotyping in monitoring Covid-19 epidemiology and SARS-CoV-2 variant dynamics in French Guiana, a tropical country located in South America. Our results suggest that wastewater-based epidemiology and genotyping are powerful tools to monitor variant introduction and disease evolution within a tropical country but the inclusion of both clinical and wastewater samples could still improve our understanding of genetic diversity co-circulating. Wastewater sequencing also revealed the cryptic transmission of SARS-CoV-2 variants within the country. Interestingly, we found some amino acid changes specific to the variants co-circulating in French Guiana, suggesting a local evolution of the SARS-CoV-2 variants after their introduction. More importantly, our results showed that the proximity to bordering countries was not the origin of the emergence of the French Guianese B.1.160.25 variant, but rather that this variant emerged from an ancestor B.1.160 variant introduced by European air plane travelers, suggesting thus that air travel remains a significant risk for cross-border spread of infectious diseases. Overall, we suggest that wastewater-based epidemiology and genotyping provides a cost effective and non-invasive approach for pathogen monitoring and an early-warning tool for disease emergence and spread within a tropical country.

A benchmark for evaluation of structure-based online tools for antibody-antigen binding affinity

The prediction of binding affinity changes caused by missense mutations can elucidate antigen-antibody interactions. A few accessible structure-based online computational tools have been proposed. However, selecting suitable software for particular research is challenging, especially research on the SARS-CoV-2 spike protein with antibodies. Therefore, benchmarking of the mutation-diverse SARS-CoV-2 datasets is critical. Here, we collected the datasets including 1216 variants about the changes in binding affinity of antigens from 22 complexes for SARS-CoV-2 S proteins and 22 monoclonal antibodies as well as applied them to evaluate the performance of seven binding affinity prediction tools. The tested tools' Pearson correlations between predicted and measured changes in binding affinity were between -0.158 and 0.657, while accuracy in classification tasks on predicting increasing or decreasing affinity ranged from 0.444 to 0.834. These tools performed relatively better on predicting single mutations, especially at epitope sites, whereas poor performance on extremely decreasing affinity. The tested tools were relatively insensitive to the experimental techniques used to obtain structures of complexes. In summary, we constructed a list of datasets and evaluated a range of structure-based online prediction tools that will explicate relevant processes of antigen-antibody interactions and enhance the computational design of therapeutic monoclonal antibodies.

Prediction of Prospective Mutational Landscape of SARS-CoV-2 Spike ssRNA and Evolutionary Basis of Its Host Interaction

Booster doses are crucial against severe COVID-19, as rapid virus mutations and variant emergence prolong the pandemic crisis. The virus's quick evolution, short generation-time, and adaptive changes impact virulence and evolvability, helping predictions about variant of concerns' (VOCs') landscapes. Here, in this study, we used a new computational algorithm, to predict the mutational pattern in SARS-CoV-2 ssRNA, proteomics, structural identification, mutation stability, and functional correlation, as well as immune escape mechanisms. Interestingly, the sequence diversity of SARS Coronavirus-2 has demonstrated a predominance of G- > A and C- > U substitutions. The best validation statistics are explored here in seven homologous models of the expected mutant SARS-CoV-2 spike ssRNA and employed for hACE2 and IgG interactions. The interactome profile of SARS-CoV-2 spike with hACE2 and IgG revealed a strong correlation between phylogeny and divergence time. The systematic adaptation of SARS-CoV-2 spike ssRNA influences infectivity and immune escape. Data suggest higher propensity of Adenine rich sequence promotes MHC system avoidance, preferred by A-rich codons. Phylogenetic data revealed the evolution of SARS-CoV-2 lineages' epidemiology. Our findings may unveil processes governing the genesis of immune-resistant variants, prompting a critical reassessment of the coronavirus mutation rate and exploration of hypotheses beyond mechanical aspects.

SARS-CoV-2 outbreak: role of viral proteins and genomic diversity in virus infection and COVID-19 progression

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection is the cause of coronavirus disease 2019 (COVID-19); a severe respiratory distress that has emerged from the city of Wuhan, Hubei province, China during December 2019. COVID-19 is currently the major global health problem and the disease has now spread to most countries in the world. COVID-19 has profoundly impacted human health and activities worldwide. Genetic mutation is one of the essential characteristics of viruses. They do so to adapt to their host or to move to another one. Viral genetic mutations have a high potentiality to impact human health as these mutations grant viruses unique unpredicted characteristics. The difficulty in predicting viral genetic mutations is a significant obstacle in the field. Evidence indicates that SARS-CoV-2 has a variety of genetic mutations and genomic diversity with obvious clinical consequences and implications. In this review, we comprehensively summarized and discussed the currently available knowledge regarding SARS-CoV-2 outbreaks with a fundamental focus on the role of the viral proteins and their mutations in viral infection and COVID-19 progression. We also summarized the clinical implications of SARS-CoV-2 variants and how they affect the disease severity and hinder vaccine development. Finally, we provided a massive phylogenetic analysis of the spike gene of 214 SARS-CoV-2 isolates from different geographical regions all over the world and their associated clinical implications.

Early predictors of intensive care unit admission among COVID-19 patients in Qatar

Background

COVID-19 is associated with significant morbidity and mortality. This study aimed to explore the early predictors of intensive care unit (ICU) admission among patients with COVID-19.

Methods

This was a case-control study of adult patients with confirmed COVID-19. Cases were defined as patients admitted to ICU during the period February 29-May 29, 2020. For each case enrolled, one control was matched by age and gender.

Results

A total of 1,560 patients with confirmed COVID-19 were included. Each group included 780 patients with a predominant male gender (89.7%) and a median age of 49 years (interquartile range = 18). Predictors independently associated with ICU admission were cardiovascular disease (adjusted odds ratio (aOR) = 1.64, 95% confidence interval (CI): 1.16-2.32, p = 0.005), diabetes (aOR = 1.52, 95% CI: 1.08-2.13, p = 0.016), obesity (aOR = 1.46, 95% CI: 1.03-2.08, p = 0.034), lymphopenia (aOR = 2.69, 95% CI: 1.80-4.02, p < 0.001), high AST (aOR = 2.59, 95% CI: 1.53-4.36, p < 0.001), high ferritin (aOR = 1.96, 95% CI: 1.40-2.74, p < 0.001), high CRP (aOR = 4.09, 95% CI: 2.81-5.96, p < 0.001), and dyspnea (aOR = 2.50, 95% CI: 1.77-3.54, p < 0.001).

Conclusion

Having cardiovascular disease, diabetes, obesity, lymphopenia, dyspnea, and increased AST, ferritin, and CRP were independent predictors for ICU admission in patients with COVID-19.

Global and genetic diversity of SARS-CoV-2 in wastewater

The analysis of SARS-CoV-2 in wastewater has enabled us to better understand the spread and evolution of the virus worldwide. To deepen our understanding of its epidemiological and genomic characteristics, we analyzed 10,147 SARS-CoV-2 sequences from 5 continents and 21 countries that were deposited in the GISAID database up until January 31, 2023. Our results revealed over 100 independent lineages of the virus circulating in water samples from March 2020 to January 2023, including variants of interest and concern. We observed four clearly defined periods of global distribution of these variants over time, with one variant being replaced by another. Interestingly, we found that SARS-CoV-2 water-borne sequences from different countries had a close phylogenetic relationship. Additionally, 40 SARS-CoV-2 water-borne sequences from Europe and the USA did not show any phylogenetic relationship with SARS-CoV-2 human sequences. We also identified a significant number of non-synonymous mutations, some of which were detected in previously reported cryptic lineages. Among the countries analyzed, France and the USA showed the highest degree of sequence diversity, while Austria reported the highest number of genomes (6,296). Our study provides valuable information about the epidemiological and genomic diversity of SARS-CoV-2 in wastewater, which can be employed to support public health initiatives and preparedness.

Strategy to develop broadly effective multivalent COVID-19 vaccines against emerging variants based on Ad5/35 platform

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron strain has evolved into highly divergent variants with several sub-lineages. These newly emerging variants threaten the efficacy of available COVID-19 vaccines. To mitigate the occurrence of breakthrough infections and re-infections, and more importantly, to reduce the disease burden, it is essential to develop a strategy for producing updated multivalent vaccines that can provide broad neutralization against both currently circulating and emerging variants. We developed bivalent vaccine AdCLD-CoV19-1 BA.5/BA.2.75 and trivalent vaccines AdCLD-CoV19-1 XBB/BN.1/BQ.1.1 and AdCLD-CoV19-1 XBB.1.5/BN.1/BQ.1.1 using an Ad5/35 platform-based non-replicating recombinant adenoviral vector. We compared immune responses elicited by the monovalent and multivalent vaccines in mice and macaques. We found that the BA.5/BA.2.75 bivalent and the XBB/BN.1/BQ.1.1 and XBB.1.5/BN.1/BQ.1.1 trivalent vaccines exhibited improved cross-neutralization ability compared to their respective monovalent vaccines. These data suggest that the developed multivalent vaccines enhance immunity against circulating Omicron subvariants and effectively elicit neutralizing antibodies across a broad spectrum of SARS-CoV-2 variants.

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

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

Novel Omicron Variants Enhance Anchored Recognition of TMEM106B: A New Pathway for SARS-CoV-2 Cellular Invasion

The recent discovery that TMEM106B serves as a receptor mediating ACE2-independent SARS-CoV-2 entry into cells deserves attention, especially in the background of the frequent emergence of mutant strains. Here, the structure-dynamic features of this novel pathway are dissected deeply. Our investigation revealed that the large loop (RBD@471-491) could anchor TMEM106B, which was then firmly locked by another loop (RBD@444-451). The novel and widely disseminated Omicron variants (BA.2.86/EG.5.1) affect the anchoring recognition of proteins, with BA.2.86 being more likely to impact cells with limited or undetectable ACE2 expression. The large loop of the EG.5.1 variant captures TMEM106B poorly due to impaired electrostatic complementarity. Furthermore, we emphasize that antibody design against these two loops could enhance the protection of ACE2 low-expressing cells according to the alanine scanning mutagenesis of multiple antibodies. We hope this study will provide a novel perspective for the prevention and treatment against this new viral invasion pathway.

Identifying novel inhibitors targeting Exportin-1 for the potential treatment of COVID-19

The nuclear export protein 1 (XPO1) mediates the nucleocytoplasmic transport of proteins and ribonucleic acids (RNAs) and plays a prominent role in maintaining cellular homeostasis. XPO1 has emerged as a promising therapeutic approach to interfere with the lifecycle of many viruses. In our earlier study, we proved the inhibition of XPO1 as a therapeutic strategy for managing SARS-COV-2 and its variants. In this study, we have utilized pharmacophore-assisted computational methods to identify prominent XPO1 inhibitors. After several layers of screening, a few molecules were shortlisted for further experimental validation on the in vitro SARS-CoV-2 cell infection model. It was observed that these compounds reduced spike positivity, suggesting inhibition of SARS-COV-2 infection. The outcome of this study could be considered further for developing novel antiviral therapeutic strategies against SARS-CoV-2.

Motifs in SARS-CoV-2 evolution

We present a novel framework enhancing the prediction of whether novel lineage poses the threat of eventually dominating the viral population. The framework is based purely on genomic sequence data, without requiring prior established biological analysis. Its building blocks are sets of coevolving sites in the alignment (motifs), identified via coevolutionary signals. The collection of such motifs forms a relational structure over the polymorphic sites. Motifs are constructed using distances quantifying the coevolutionary coupling of pairs and manifest as coevolving clusters of sites. We present an approach to genomic surveillance based on this notion of relational structure. Our system will issue an alert regarding a lineage, based on its contribution to drastic changes in the relational structure. We then conduct a comprehensive retrospective analysis of the COVID-19 pandemic based on SARS-CoV-2 genomic sequence data in GISAID from October 2020 to September 2022, across 21 lineages and 27 countries with weekly resolution. We investigate the performance of this surveillance system in terms of its accuracy, timeliness, and robustness. Lastly, we study how well each lineage is classified by such a system.

Structure-Based High-Throughput Virtual Screening and Molecular Dynamics Simulation for the Discovery of Novel SARS-CoV-2 NSP3 Mac1 Domain Inhibitors

Since the emergence of SARS-CoV-2, many genetic variations within its genome have been identified, but only a few mutations have been found in nonstructural proteins (NSPs). Among this class of viral proteins, NSP3 is a multidomain protein with 16 different domains, and its largest domain is known as the macrodomain or Mac1 domain. In this study, we present a virtual screening campaign in which we computationally evaluated the NCI anticancer library against the NSP3 Mac1 domain, using Molegro Virtual Docker. The top hits with the best MolDock and Re-Rank scores were selected. The physicochemical analysis and drug-like potential of the top hits were analyzed using the SwissADME data server. The binding stability and affinity of the top NSC compounds against the NSP3 Mac1 domain were analyzed using molecular dynamics (MD) simulation, using Desmond software, and their interaction energies were analyzed using the MM/GBSA method. In particular, by applying subsequent computational filters, we identified 10 compounds as possible NSP3 Mac1 domain inhibitors. Among them, after the assessment of binding energies (ΔGbind) on the whole MD trajectories, we identified the four most interesting compounds that acted as strong binders of the NSP3 Mac1 domain (NSC-358078, NSC-287067, NSC-123472, and NSC-142843), and, remarkably, it could be further characterized for developing innovative antivirals against SARS-CoV-2.

Novel Mutations in the Non-Structure Protein 2 of SARS-CoV-2

Introduction

Mutation in the genome of SARS-CoV-2 may play a role in immune evasion, pathogenicity, and speed of its transmission. Our investigation aimed to evaluate the mutations that exist in the NSP2.

Materials and Method

RNA was extracted from nasopharyngeal swabs from 100 COVID-19 patients. RT-PCR was performed on all samples using NSP2-specific primers. Following gel electrophoresis, the bands were cut, purified, and sequenced using the Sanger method. After sequencing, 90 sequences could be used for further analysis. Bioinformatics analysis was conducted to investigate the effect of mutations on protein structure, stability, prediction of homology models, and phylogeny tree.

Results

The patients' mean age was 51.08. The results revealed that 8 of the 17 NSP2 mutations (R207C, T224I, G262V, T265I, K337D, N348S, G392D, and I431M) were missense. One deletion was also found in NSP2. Among NSP2 missense mutations studied, K337D and G392D increased structural stability while the others decreased it. The homology-designed models demonstrated that the homologies were comparable to the sequences of the Wuhan-HU-1 virus.

Conclusion

Our study suggested that the mutations K337D and G392D modulate the stability of NSP2, and tracking viral evolution should be implemented and vaccine development updated.

COVID-19 in 16 West African Countries: An Assessment of the Epidemiology and Genetic Diversity of SARS-CoV-2 after Four Epidemic Waves

West Africa faced the COVID-19 pandemic in early March 2020 and, as of March 31, 2022, had more than 900,000 confirmed cases and more than 12,000 deaths. During this period, SARS-CoV-2 genomes evolved genetically, resulting in the emergence of distinct lineages. This review was conducted to provide the epidemiological profile of COVID-19, the mutational profile of SARS-CoV-2, and the dynamics of its lineages in the 16 west African countries by analyzing data from 33 studies and seven situation reports. For a more complete representation of the epidemiology and genetic diversity of SARS-CoV-2, we used reliable public data in addition to eligible studies. As of March 31, 2022, the 16 west African countries experienced four epidemic waves with variable intensities. Higher mortality was noted during the third wave with a case fatality rate (CFR) of 1.9%. After these four epidemic waves, Liberia recorded the highest CFR (4.0%), whereas Benin had the lowest CFR (0.6%). Through mutational analysis, a high genetic heterogeneity of the genomes was observed, with a predominance of mutations in the spike protein. From this high mutational rate, different lineages emerged. Our analysis of the evolutionary diversity allowed us to count 205 lineages circulating in west Africa. This study has provided a good representation of the mutational profile and the prevalence of SARS CoV-2 lineages beyond the knowledge of the global epidemiology of the 16 African countries.

Multiple thermocycles followed by LAMP with only two primers for ultrasensitive colorimetric viral RNA testing and tracking at single-base resolution

Rapid, accurate and high throughput measurement of infectious viruses is an urgent need to prevent viral transmission. Loop-mediated isothermal amplification (LAMP) is an attractive isothermal amplification method for nucleic acid detection, especially for point-of-care (POC) testing, but it needs at least four primers and its sensitivity is also limited when integrating with visual detection methods. Herein, by designing only two primers to precisely recognize the four regions of the target, we developed a multiple thermocycles-based LAMP method (MTC-LAMP) for sensitive and specific testing and tracking of viral RNA. We also introduced a novel SYBR Green I (SG)-assisted stable colorimetric assay induced by the amplification products through the charge neutralization effect of positively charged SG toward gold nanoparticles (AuNPs). The ultralow nonspecific background of the double exponential amplification improved the detection sensitivity to near single-molecule level (1 aM, 3 copies in 5 μL solution), which was higher than RT-PCR and RT-LAMP. After adding AuNPs, a significant color difference between target and blank was immediately observed by naked eye. By introducing a peptide nucleic acid (PNA) clamp into our colorimetric MTC-LAMP assay, the specific distinguish of virus variants at single-base resolution was observed without the requirement of any equipment. This assay shows great potential for large-scale screening and tracking of the threatening viruses with ultrahigh sensitivity and pronounced colorimetric output, which is of great importance for pandemic control.

CovidShiny: An Integrated Web Tool for SARS-CoV-2 Mutation Profiling and Molecular Diagnosis Assay Evaluation In Silico

The coronavirus disease 2019 (COVID-19) pandemic is still ongoing, with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continuing to evolve and accumulate mutations. While various bioinformatics tools have been developed for SARS-CoV-2, a well-curated mutation-tracking database integrated with in silico evaluation for molecular diagnostic assays is currently unavailable. To address this, we introduce CovidShiny, a web tool that integrates mutation profiling, in silico evaluation, and data download capabilities for genomic sequence-based SARS-CoV-2 assays and data download. It offers a feasible framework for surveilling the mutation of SARS-CoV-2 and evaluating the coverage of the molecular diagnostic assay for SARS-CoV-2. With CovidShiny, we examined the dynamic mutation pattern of SARS-CoV-2 and evaluated the coverage of commonly used assays on a large scale. Based on our in silico analysis, we stress the importance of using multiple target molecular diagnostic assays for SARS-CoV-2 to avoid potential false-negative results caused by viral mutations. Overall, CovidShiny is a valuable tool for SARS-CoV-2 mutation surveillance and in silico assay design and evaluation.

Correlative analysis of wastewater trends with clinical cases and hospitalizations through five dominant variant waves of COVID-19

Wastewater-based epidemiology (WBE) has been utilized to track community infections of Coronavirus Disease 2019 (COVID-19) by detecting RNA of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), within samples collected from wastewater. The correlations between community infections and wastewater measurements of the RNA can potentially change as SARS-CoV-2 evolves into new variations by mutating. This study analyzed SARS-CoV-2 RNA, and indicators of human waste in wastewater from two sewersheds of different scales (University of Miami (UM) campus and Miami-Dade County Central District wastewater treatment plant (CDWWTP)) during five internally defined COVID-19 variant dominant periods (Initial, Pre-Delta, Delta, Omicron and Post-Omicron wave). SARS-CoV-2 RNA quantities were compared against COVID-19 clinical cases and hospitalizations to evaluate correlations with wastewater SARS-CoV-2 RNA. Although correlations between documented clinical cases and hospitalizations were high, prevalence for a given wastewater SARS-CoV-2 level varied depending upon the variant analyzed. The correlative relationship was significantly steeper (more cases per level found in wastewater) for the Omicron-dominated period. For hospitalization, the relationships were steepest for the Initial wave, followed by the Delta wave with flatter slopes during all other waves. Overall results were interpreted in the context of SARS-CoV-2 virulence and vaccination rates among the community.

Genome analysis of SARS-CoV-2 isolates from a population reveals the rapid selective sweep of a haplotype carrying many pre-existing and new mutations

To understand the mechanism underlying the evolution of SARS-CoV-2 in a population, we sequenced 92 viral genomes from Assam, India. Analysis of these and database sequences revealed a complete selective sweep of a haplotype in Assam carrying 13 pre-existing variants, including a high leap in frequency of a variant on ORF8, which is involved in immune evasion. A comparative study between sequences of same lineage and similar time frames in and outside Assam showed that 10 of the 13 pre-existing variants had a frequency ranging from 96 to 99%, and the remaining 3 had a low frequency outside Assam. Using a phylogenetic approach to infer sequential occurrences of variants we found that the variant Phe120del on ORF8, which had a low frequency (1.75%) outside Assam, is at the base of the phylogenetic tree of variants and became totally fixed (100%) in Assam population. Based on this observation, we inferred that the variant on ORF8 had a selective advantage, so it carried the haplotype to reach the100% frequency. The haplotype also carried 32 pre-existing variants at a frequency from 1.00 to 80.00% outside Assam. Those of these variants that are more closely linked to the S-protein locus, which often carries advantageous mutations and is tightly linked to the ORF8 locus, retained higher frequencies, while the less tightly linked variants showed lower frequencies, likely due to recombination among co- circulating variants in Assam. The ratios of non-synonymous substitutions to synonymous substitutions suggested that some genes such as those coding for the S-protein and non-structural proteins underwent positive selection while others were subject to purifying selection during their evolution in Assam. Furthermore, we observed negative correlation of the Ct value of qRT-PCR of the patients with abundant ORF6 transcripts, suggesting that ORF6 can be used as a marker for estimating viral titer. In conclusion, our in-depth analysis of SARS-CoV-2 genomes in a regional population reveals the mechanism and dynamics of viral evolution.

Unraveling the impact of ORF3a Q57H mutation on SARS-CoV-2: insights from molecular dynamics

ORF3a is a conserved accessory protein of SARS-CoV-2, linked to viral infection and pathogenesis, with acquired mutations at various locations. Previous studies have shown that the occurrence of the Q57H mutation is higher in comparison to other positions in ORF3a. This mutation is known to induce conformational changes, yet the extent of structural alteration and its role in the viral adaptation process remain unknown. Here we performed molecular dynamics (MD) simulations of wt-ORF3a, Q57H, and Q57A mutants to analyze structural changes caused by mutations compared to the native protein. The MD analysis revealed that Q57H and Q57A mutants show significant structural changes in the dimer conformation than the wt-ORF3a. This dimer conformer narrows down the ion channel cavity, which reduces Na + or K + permeability leading to decrease the antigenic response that can help the virus to escape the host immune system. Non-bonding interaction analysis shows the Q57H mutant has more interacting residues, resulting in more stability within dimer conformation than the wt-ORF3a and Q57A. Moreover, both mutant dimers (Q57H and Q57A) form a novel salt-bridge interaction at the same position between A:Asp142 and B:Lys61, whereas such an interaction is absent in the wt-ORF3a dimer. We have also noticed that the TM3 domain's flexibility in Q57H is increased because of strong inter-domain interactions of TM1 and TM2 within the dimer conformation. These unusual interactions and flexibility of Q57H mutant can have significant impacts on the SARS-CoV-2 adaptations, virulence, transmission, and immune system evasion. Our findings are consistent with the previous experimental data and provided details information on the structural perturbation in ORF3a caused by mutations, which can help better understand the structural change at the molecular level as well as the reason for the high virulence properties of this variant.Communicated by Ramaswamy H. Sarma.

Timeline of SARS-CoV-2 Transmission in Sabah, Malaysia: Tracking the Molecular Evolution

BACKGROUND

The COVID-19 pandemic poses an unprecedented public health challenge in Malaysia. The impact of COVID-19 varies between countries, including geographically divided states within a country. The deadly transmission of COVID-19 has taken a heavy toll in Sabah, Malaysia's third most populous state, contributing nearly 10% to the recorded national death toll as of 31 December 2022. Although several SARS-CoV-2 genome sequences have been analysed in Malaysia, molecular epidemiology data from Sabah focusing on the diversity and evolution of SARS-CoV-2 variants are still lacking. This study examines the major SARS-CoV-2 variants and emerging mutations from Sabah, the Malaysian Borneo, which is geographically divided from West Malaysia by the South China Sea.

METHODS

A total of 583 COVID-19 samples were subjected to whole genome sequencing and analysed with an additional 1123 Sabah COVID-19 sequences retrieved from the GISAID EpiCoV consortium. Nextclade and Pangolin were used to classify these sequences according to the clades and lineages. To determine the molecular evolutionary characteristics, Bayesian time-scaled phylogenetic analysis employing the maximum likelihood algorithm was performed on selected SARS-CoV-2 genome sequences, using the Wuhan-Hu-1 sequence as a reference.

RESULTS

Sabah was affected starting from the second COVID-19 wave in Malaysia, and the early sequences were classified under the O clade. The clade was gradually replaced during subsequent waves by G, GH, GK and GRA, with the latter being dominant as of December 2022. Phylogenetically, the Delta isolates in this study belong to the three main subclades 21A, 21J and 21I, while Omicron isolates belong to 21M, 21L and 22B. The time-scaled phylogeny suggested that SARS-CoV-2 introduced into Sabah originated from Peninsular Malaysia in early March 2020, and phylodynamic analysis indicated that increased viral spread was observed in early March and declined in late April, followed by an evolutionary stationary phase in June 2020.

CONCLUSION

Continuous molecular epidemiology of SARS-CoV-2 in Sabah will provide a deeper understanding of the emergence and dominance of each variant in the locality, thus facilitating public health intervention measures.

Comparison of genome replication fidelity between SARS-CoV-2 and influenza A virus in cell culture

Since the emergence of COVID-19, several SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) variants have emerged and spread widely. These variants are produced through replication errors of the viral genome by viral RNA-dependent RNA polymerase (RdRp). Seasonal epidemics of influenza are also known to occur because of new variants of influenza A virus (IAV), which are generated by the introduction of mutations by viral RdRp with low fidelity. Variants with different antigenicities appear because of mutations in envelope glycoproteins. In this study, we calculated and compared the mutation rates in genome replication of IAV and SARS-CoV-2. Average mutation rates per passage were 9.01 × 10-5 and 3.76 × 10-6 substitutions/site for IAV and SARS-CoV-2, respectively. The mutation rate of SARS-CoV-2 was 23.9-fold lower than that of IAV because of the proofreading activity of the SARS-CoV-2 RdRp complex. Our data could be useful in establishing effective countermeasures against COVID-19.

ACE2 receptor polymorphism in humans and animals increases the risk of the emergence of SARS-CoV-2 variants during repeated intra- and inter-species host-switching of the virus

Like other coronaviruses, SARS-CoV-2 has ability to spread through human-to-human transmission and to circulate from humans to animals and from animals to humans. A high frequency of SARS-CoV-2 mutations has been observed in the viruses isolated from both humans and animals, suggesting a genetic fitness under positive selection in both ecological niches. The most documented positive selection force driving SARS-CoV-2 mutations is the host-specific immune response. However, after electrostatic interactions with lipid rafts, the first contact between the virus and host proteins is the viral spike-cellular receptor binding. Therefore, it is likely that the first level of selection pressure impacting viral fitness relates to the virus's affinity for its receptor, the angiotensin I converting enzyme 2 (ACE2). Although sufficiently conserved in a huge number of species to support binding of the viral spike with enough affinity to initiate fusion, ACE2 is highly polymorphic both among species and within a species. Here, we provide evidence suggesting that when the viral spike-ACE2 receptor interaction is not optimal, due to host-switching, mutations can be selected to improve the affinity of the spike for the ACE2 expressed by the new host. Notably, SARS-CoV-2 is mutation-prone in the spike receptor binding domain (RBD), allowing a better fit for ACE2 orthologs in animals. It is possibly that this may also be true for rare human alleles of ACE2 when the virus is spreading to billions of people. In this study, we present evidence that human subjects expressing the rare E₃₂₉G allele of ACE2 with higher allele frequencies in European populations exhibit a improved affinity for the SARS-CoV-2 spike N₅₀₁Y variant of the virus. This may suggest that this viral N₅₀₁Y variant emerged in the human population after SARS-CoV-2 had infected a human carrying the rare E₃₂₉G allele of ACE2. In addition, this viral evolution could impact viral replication as well as the ability of the adaptive humoral response to control infection with RBD-specific neutralizing antibodies. In a shifting landscape, this ACE2-driven genetic drift of SARS-CoV-2 which we have named the 'boomerang effect', could complicate the challenge of preventing COVID with a SARS-CoV-2 spike-derived vaccine.

The NSP4 T492I mutation increases SARS-CoV-2 infectivity by altering non-structural protein cleavage

The historically dominant SARS-CoV-2 Delta variant and the currently dominant Omicron variants carry a T492I substitution within the non-structural protein 4 (NSP4). Based on in silico analyses, we hypothesized that the T492I mutation increases viral transmissibility and adaptability, which we confirmed with competition experiments in hamster and human airway tissue culture models. Furthermore, we showed that the T492I mutation increases the replication capacity and infectiveness of the virus and improves its ability to evade host immune responses. Mechanistically, the T492I mutation increases the cleavage efficiency of the viral main protease NSP5 by enhancing enzyme-substrate binding, which increases production of nearly all non-structural proteins processed by NSP5. Importantly, the T492I mutation suppresses viral-RNA-associated chemokine production in monocytic macrophages, which may contribute to the attenuated pathogenicity of Omicron variants. Our results highlight the importance of NSP4 adaptation in the evolutionary dynamics of SARS-CoV-2.

Circulating Variants of SARS-CoV-2 Among Macedonian COVID - 19 Patients in the First Year of Pandemic

Genomic epidemiology has proven to be a useful tool for investigating pandemic outbreaks and tracking pathogen spread and evolution. This study describes the circulation of SARS-CoV-2 strains in N. Macedonia during a period of one year, encompassing three waves of the COVID-19 pandemic. A certain percentage (2-3%) of positive cases were continuously selected and analyzed by whole genome sequencing (WGS) technology. Using this approach, a total of 337 SARS-CoV-2 genomes were sequenced and 26 different lineages belonging to 7 clades were detected. During the first wave of the pandemic, the most dominant lineage was B.1.1, followed by B.1.1.70, which became the most dominant in the second wave. The B.1.1.7 lineage completely overpassed all other variants in the third wave. Our study strengthens the notion that the progression of COVID-19 pandemic is associated with emergence of new SARS-CoV-2 variants with increased virulence. The measure of the impact of this viral dynamic on the spread of the pandemic should be evaluated in association with other factors that might influence the transmission.

Sequence Analysis of Hot Spot Regions of Spike and RNA-dependent-RNA polymerase (RdRp) Genes of SARS-CoV-2 in Kerman, Iran

Background

Mutations in the SARS-CoV-2 genome might influence pathogenicity, transmission rate, and evasion of the host immune system. Therefore, the purpose of the present study was to investigate the genetic alteration as well as assess their effects on the receptor binding domain (RBD) of the spike and the putative RNA binding site of the RdRp genes of SARS-CoV-2 using bioinformatics tools.

Materials and Method

In this cross-sectional study, 45 confirmed COVID-19 patients using qRT-PCR were included and divided into mild, severe, and critical groups based on the severity of the disease. RNA was extracted from nasopharyngeal swab samples using a commercial kit. RT-PCR was performed to amplify the target sequences of the spike and RdRp genes and sequence them by the Sanger method. Clustal OMEGA, MEGA 11 software, I-mutant tools, SWISS-MODEL, and HDOCK web servers were used for bioinformatics analyses.

Results

The mean age of the patients was 50.68±2.73. The results showed that four of six mutations (L452R, T478K, N501Y, and D614G) in RBD and three of eight in the putative RNA binding site (P314L, E1084D, V1883T) were missense. In the putative RNA binding site, another deletion was discovered. Among missense mutations, N501Y and V1883T were responsible for increasing structural stability, while others were responsible for decreasing it. The various homology models designed showed that these homologies were like the Wuhan model. The molecular docking analysis revealed that the T478K mutation in RBD had the highest binding affinity. In addition, 35 RBD samples (89.7%) and 33 putative RNA binding site samples (84.6%) were similar to the Delta variant.

Conclusion

Our results indicated that double mutations (T478K and N501Y) in the S protein might increase the binding affinity of SARS-CoV-2 to human ACE2 compared to the wild-type (WT) strain. Moreover, variations in the spike and RdRp genes might influence the stability of encoded proteins.

Functional Foods: A Promising Strategy for Restoring Gut Microbiota Diversity Impacted by SARS-CoV-2 Variants

Natural herbs and functional foods contain bioactive molecules capable of augmenting the immune system and mediating anti-viral functions. Functional foods, such as prebiotics, probiotics, and dietary fibers, have been shown to have positive effects on gut microbiota diversity and immune function. The use of functional foods has been linked to enhanced immunity, regeneration, improved cognitive function, maintenance of gut microbiota, and significant improvement in overall health. The gut microbiota plays a critical role in maintaining overall health and immune function, and disruptions to its balance have been linked to various health problems. SARS-CoV-2 infection has been shown to affect gut microbiota diversity, and the emergence of variants poses new challenges to combat the virus. SARS-CoV-2 recognizes and infects human cells through ACE2 receptors prevalent in lung and gut epithelial cells. Humans are prone to SARS-CoV-2 infection because their respiratory and gastrointestinal tracts are rich in microbial diversity and contain high levels of ACE2 and TMPRSS2. This review article explores the potential use of functional foods in mitigating the impact of SARS-CoV-2 variants on gut microbiota diversity, and the potential use of functional foods as a strategy to combat these effects.

The Mutational Landscape of SARS-CoV-2

Mutation research is crucial for detecting and treating SARS-CoV-2 and developing vaccines. Using over 5,300,000 sequences from SARS-CoV-2 genomes and custom Python programs, we analyzed the mutational landscape of SARS-CoV-2. Although almost every nucleotide in the SARS-CoV-2 genome has mutated at some time, the substantial differences in the frequency and regularity of mutations warrant further examination. C>U mutations are the most common. They are found in the largest number of variants, pangolin lineages, and countries, which indicates that they are a driving force behind the evolution of SARS-CoV-2. Not all SARS-CoV-2 genes have mutated in the same way. Fewer non-synonymous single nucleotide variations are found in genes that encode proteins with a critical role in virus replication than in genes with ancillary roles. Some genes, such as spike (S) and nucleocapsid (N), show more non-synonymous mutations than others. Although the prevalence of mutations in the target regions of COVID-19 diagnostic RT-qPCR tests is generally low, in some cases, such as for some primers that bind to the N gene, it is significant. Therefore, ongoing monitoring of SARS-CoV-2 mutations is crucial. The SARS-CoV-2 Mutation Portal provides access to a database of SARS-CoV-2 mutations.

4’-fluorouridine as a potential COVID-19 oral drug?: a review

The available antiviral drugs against coronavirus disease 2019 (COVID-19) are limited. Oral drugs that can be prescribed to non-hospitalized patients are required. The 4′-fluoruridine, a nucleoside analog similar to remdesivir, is one of the promising candidates for COVID-19 oral therapy due to its ability to stall viral RdRp. Available data suggested that 4'-fluorouridine has antiviral activity against the respiratory syncytial virus, hepatitis C virus, lymphocytic choriomeningitis virus, and other RNA viruses, including SARS-CoV-2. In vivo study revealed that SARS-CoV-2 is highly susceptible to 4'-fluorouridine and was effective with a single daily dose versus molnupiravir administered twice daily. Although  4'-fluorouridine is considered as strong candidates, further studies are required to determine its efficacy in the patients and  it’s genetic effects on humans. In this review, we the antiviral activity of 4′-fluorouridine is reviewed and compared it to other drugs currently in development. The current literature on 4′-fluorouridine's antiviral activity against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is compiled and discussed.

Density estimation of SARS-CoV2 spike proteins using super pixels segmentation technique

The worldwide outbreak of COVID-19 disease was caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV 2). The existence of spike proteins, which allow these viruses to infect host cells, is one of the distinctive biological traits of various prior viruses. As a result, the process by which these viruses infect people is largely dependent on spike proteins. The density of SARS-CoV-2 spike proteins must be estimated to better understand and develop diagnostics and vaccines against the COVID-19 pandemic. CT scans and X-rays have three issues: frosted glass, consolidation, and strange roadway layouts. Each of these issues can be graded separately or together. Although CT scan is sensitive to COVID-19, it is not very specific. Therefore, patients who obtain these results should have more comprehensive clinical and laboratory tests to rule out other probable reasons. This work collected 586 SARS-CoV 2 transmission electron microscopy (TEM) images from open source for density estimation of virus spike proteins through a segmentation approach based on the superpixel technique. As a result, the spike density means of SARS-CoV2 and SARS-CoV were 21,97 nm and 22,45 nm, respectively. Furthermore, in the future, we aim to include this model in an intelligent system to enhance the accuracy of viral detection and classification. Moreover, we can remotely connect hospitals and public sites to conduct environmental hazard assessments and data collection.

HaploCoV: unsupervised classification and rapid detection of novel emerging variants of SARS-CoV-2

Accurate and timely monitoring of the evolution of SARS-CoV-2 is crucial for identifying and tracking potentially more transmissible/virulent viral variants, and implement mitigation strategies to limit their spread. Here we introduce HaploCoV, a novel software framework that enables the exploration of SARS-CoV-2 genomic diversity through space and time, to identify novel emerging viral variants and prioritize variants of potential epidemiological interest in a rapid and unsupervised manner. HaploCoV can integrate with any classification/nomenclature and incorporates an effective scoring system for the prioritization of SARS-CoV-2 variants. By performing retrospective analyses of more than 11.5 M genome sequences we show that HaploCoV demonstrates high levels of accuracy and reproducibility and identifies the large majority of epidemiologically relevant viral variants - as flagged by international health authorities - automatically and with rapid turn-around times.Our results highlight the importance of the application of strategies based on the systematic analysis and integration of regional data for rapid identification of novel, emerging variants of SARS-CoV-2. We believe that the approach outlined in this study will contribute to relevant advances to current and future genomic surveillance methods.

Analysis of SARS-CoV-2 variants from patient specimens in Nevada from October 2020 to August 2021

In early 2020, the emergence and spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the human population quickly developed into a global pandemic. SARS-CoV-2 is the etiological agent of coronavirus disease 2019 (COVID-19) which has a broad range of respiratory illnesses. As the virus circulates, it acquires nucleotide changes. These mutations are potentially due to the inherent differences in the selection pressures within the human population compared to the original zoonotic reservoir of SARS-CoV-2 and formerly naïve humans. The acquired mutations will most likely be neutral, but some may have implications for viral transmission, disease severity and resistance to therapies or vaccines. This is a follow-up study from our early report (Hartley et al. J Genet Genomics. 01202021;48(1):40-51) which detected a rare variant (nsp12, RdRp P323F) circulating within Nevada in mid 2020 at high frequency. The primary goals of the current study were to determine the phylogenetic relationship of the SARS-CoV-2 genomes within Nevada and to determine if there are any unusual variants within Nevada compared to the current database of SARS-CoV-2 sequences. Whole genome sequencing and analysis of SARS-CoV-2 from 425 positively identified nasopharyngeal/nasal swab specimens were performed from October 2020 to August 2021 to determine any variants that could result in potential escape from current therapeutics. Our analysis focused on nucleotide mutations that generated amino acid variations in the viral Spike (S) protein, Receptor binding domain (RBD) and the RNA-dependent RNA-polymerase (RdRp) complex. The data indicate that SARS-CoV-2 sequences from Nevada did not contain any unusual variants that had not been previously reported. Additionally, we did not detect the previously identified the RdRp P323F variant in any of the samples. This suggests that the rare variant we detected before was only able to circulate because of the stay-at-home orders and semi-isolation experience during the early months of the pandemic. IMPORTANCE: SARS-COV-2 continues to circulate in the human population. In this study, SARS-CoV-2 positive nasopharyngeal/nasal swab samples were used for whole genome sequencing to determine the phylogenetic relationship of SARS-CoV-2 sequences within Nevada from October 2020 to August 2021. The resulting data is being added to a continually growing database of SARS-CoV-2 sequences that will be important for understanding the transmission and evolution of the virus as it spreads around the globe.

The Future of Epidemic and Pandemic Vaccines to Serve Global Public Health Needs

This Review initiates a wide-ranging discussion over 2023 by selecting and exploring core themes to be investigated more deeply in papers submitted to the Vaccines Special Issue on the "Future of Epidemic and Pandemic Vaccines to Serve Global Public Health Needs". To tackle the SARS-CoV-2 pandemic, an acceleration of vaccine development across different technology platforms resulted in the emergency use authorization of multiple vaccines in less than a year. Despite this record speed, many limitations surfaced including unequal access to products and technologies, regulatory hurdles, restrictions on the flow of intellectual property needed to develop and manufacture vaccines, clinical trials challenges, development of vaccines that did not curtail or prevent transmission, unsustainable strategies for dealing with variants, and the distorted allocation of funding to favour dominant companies in affluent countries. Key to future epidemic and pandemic responses will be sustainable, global-public-health-driven vaccine development and manufacturing based on equitable access to platform technologies, decentralised and localised innovation, and multiple developers and manufacturers, especially in low- and middle-income countries (LMICs). There is talk of flexible, modular pandemic preparedness, of technology access pools based on non-exclusive global licensing agreements in exchange for fair compensation, of WHO-supported vaccine technology transfer hubs and spokes, and of the creation of vaccine prototypes ready for phase I/II trials, etc. However, all these concepts face extraordinary challenges shaped by current commercial incentives, the unwillingness of pharmaceutical companies and governments to share intellectual property and know-how, the precariousness of building capacity based solely on COVID-19 vaccines, the focus on large-scale manufacturing capacity rather than small-scale rapid-response innovation to stop outbreaks when and where they occur, and the inability of many resource-limited countries to afford next-generation vaccines for their national vaccine programmes. Once the current high subsidies are gone and interest has waned, sustaining vaccine innovation and manufacturing capability in interpandemic periods will require equitable access to vaccine innovation and manufacturing capabilities in all regions of the world based on many vaccines, not just "pandemic vaccines". Public and philanthropic investments will need to leverage enforceable commitments to share vaccines and critical technology so that countries everywhere can establish and scale up vaccine development and manufacturing capability. This will only happen if we question all prior assumptions and learn the lessons offered by the current pandemic. We invite submissions to the special issue, which we hope will help guide the world towards a global vaccine research, development, and manufacturing ecosystem that better balances and integrates scientific, clinical trial, regulatory, and commercial interests and puts global public health needs first.

Exploration of potent antiviral phytomedicines from Lauraceae family plants against SARS-CoV-2 RNA-dependent RNA polymerase

RNA-dependent RNA polymerase, also known as RdRp, is a possible therapeutic target that could be used to suppress the proliferation of RNA viruses such as SARS-CoV-2. This protein has two major functional sites (a) catalytic and (b) substrate entry, which regulate the natural substrate entry and its corresponding interaction with the protein. In this study, a computational drug design pipeline was applied to investigate potential inhibitors against SARS-CoV-2 RdRp from Lauraceae plants, and five top hits were selected based on the docked score (< -7 kcal/mol). The docking study suggested that the Glochidioboside had a minimum binding score of -7.8 kcal/mol. This compound showed total five hydrogen bonds while two of them were with catalytic residues Asp618 and Asp760. However, another compound, Sitogluside showed a binding score of -7.3 kcal/mol with four hydrogen bonds targeting three functional residues (Arg555, Ser759, and Asp760). Later, 100 ns explicit solvent molecular dynamics (MD) simulation was performed to evaluate the stability of the protein-ligand docked system. These compounds translocated their positions from the catalytic site to the substrate entry site, as observed in the MD simulation trajectory. However, translocation did not affect the binding strength of these compounds, and they retained the strong binding affinity (ΔG < -11.5 kcal/mol), estimated using the MM/GBSA method. In general, the findings of this study indicated the potential therapeutic compounds that may be used targeting SARS-CoV-2 RdRp. However, these compounds still need to be validated by experimentation in order to determine their inhibitory function.Communicated by Ramaswamy H. Sarma.

Fidelity of Ribonucleotide Incorporation by the SARS-CoV-2 Replication Complex

The SARS-CoV-2 coronavirus has caused a global pandemic. Despite the initial success of vaccines at preventing infection, genomic variation has led to the proliferation of variants capable of higher infectivity. Mutations in the SARS-CoV-2 genome are the consequence of replication errors, highlighting the importance of understanding the determinants of SARS-CoV-2 replication fidelity. The RNA-dependent RNA polymerase (RdRp) is the central catalytic subunit for SARS-CoV-2 RNA replication and genome transcription. Here, we report the fidelity of ribonucleotide incorporation by SARS-CoV-2 RdRp (nsp12), along with its co-factors nsp7/nsp8, using steady-state kinetic analysis. Our analysis suggests that in the absence of the proofreading subunit (nsp14), the nsp12/7/8 complex has a surprisingly low base substitution fidelity (10^-1-10^-3). This is orders of magnitude lower than the fidelity reported for other coronaviruses (10^-6-10^-7), highlighting the importance of proofreading for faithful SARS-CoV-2 replication. We performed a mutational analysis of all reported SARS-CoV-2 genomes and identified mutations in both nsp12 and nsp14 that appear likely to lower viral replication fidelity through mechanisms that include impairing the nsp14 exonuclease activity or its association with the RdRp. Our observations provide novel insight into the mechanistic basis of replication fidelity in SARS-CoV-2 and the potential effect of nsp12 and nsp14 mutations on replication fidelity, informing the development of future antiviral agents and SARS-CoV-2 vaccines.

Computational Analysis Predicts Correlations among Amino Acids in SARS-CoV-2 Proteomes

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a serious global challenge requiring urgent and permanent therapeutic solutions. These solutions can only be engineered if the patterns and rate of mutations of the virus can be elucidated. Predicting mutations and the structure of proteins based on these mutations have become necessary for early drug and vaccine design purposes in anticipation of future viral mutations. The amino acid composition (AAC) of proteomes and individual viral proteins provide avenues for exploitation since AACs have been previously used to predict structure, shape and evolutionary rates. Herein, the frequency of amino acid residues found in 1637 complete proteomes belonging to 11 SARS-CoV-2 variants/lineages were analyzed. Leucine is the most abundant amino acid residue in the SARS-CoV-2 with an average AAC of 9.658% while tryptophan had the least abundance of 1.11%. The AAC and ranking of lysine and glycine varied in the proteome. For some variants, glycine had higher frequency and AAC than lysine and vice versa in other variants. Tryptophan was also observed to be the most intolerant to mutation in the various proteomes for the variants used. A correlogram revealed a very strong correlation of 0.999992 between B.1.525 (Eta) and B.1.526 (Iota) variants. Furthermore, isoleucine and threonine were observed to have a very strong negative correlation of -0.912, while cysteine and isoleucine had a very strong positive correlation of 0.835 at p < 0.001. Shapiro-Wilk normality test revealed that AAC values for all the amino acid residues except methionine showed no evidence of non-normality at p < 0.05. Thus, AACs of SARS-CoV-2 variants can be predicted using probability and z-scores. AACs may be beneficial in classifying viral strains, predicting viral disease types, members of protein families, protein interactions and for diagnostic purposes. They may also be used as a feature along with other crucial factors in machine-learning based algorithms to predict viral mutations. These mutation-predicting algorithms may help in developing effective therapeutics and vaccines for SARS-CoV-2.

Structural patterns of SARS-CoV-2 variants of concern (alpha, beta, gamma, delta) spike protein are influenced by variant-specific amino acid mutations: A computational study with implications on viral evolution

SARS-CoV-2 (SARS2) regularly mutates resulting to variants of concern (VOC) which have higher virulence and transmissibility rates while concurrently evading available therapeutic strategies. This highlights the importance of amino acid mutations occurring in the SARS2 spike protein structure since it may affect virus biology. However, this was never fully elucidated. Here, network analysis was performed based on the COVID-19 genomic epidemiology network between December 2019-July 2021. Representative SARS2 VOC spike protein models were generated and quality checked, protein model superimposition was done, and common contact based on contact mapping was established. Throughout this study, we found that: (1) certain individual variant-specific amino acid mutations can affect the spike protein structural pattern; (2) certain individual variant-specific amino acid mutations had no affect on the spike protein structural pattern; and (3) certain combination of variant-specific amino acids are putatively epistatic mutations that can potentially influence the VOC spike protein structural pattern. This manuscript was submitted as part of a theme issue on "Modelling COVID-19 and Preparedness for Future Pandemics".

Sequence analysis of SARS-CoV-2 Delta variant isolated from Makassar, South Sulawesi, Indonesia

Introduction

This study aimed to perform mutation and phylogenetic analyses of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Delta variants and analyze the characteristic signs and symptoms of patients infected with SARS-CoV-2 Delta variant originated from Makassar during the Delta outbreak.Methods: We collected samples from patients who were infected with coronavirus disease 2019 (COVID-19) between June and October 2021. We selected the Quantitative Reverse Transcription-Polymerase Chain Reaction (PCR)-positive samples with a cycle threshold value of <30 for whole genome sequencing. Total viral ribonucleic acid (RNA) was isolated from 34 PCR-positive nasopharyngeal swab samples, and whole genome sequencing was performed using the Oxford Nanopore GridlON sequencer. Phylogenetic and maximum clade credibility analyses were performed using the Bayesian Markov chain Monte Carlo method.

Results

It was found that 33 patients were infected with the SARS-CoV-2 Delta variant in this cohort study, among whom 63.6% (21) patients were female. According to the clinical data, 24 (72.7%), 7 (21.2%), and 2 (6.1%) patients had mild, moderate, and severe COVID-19 infections. Phylogenetic analysis based on the spike and RNA-dependent RNA polymerase (RdRp) genes showed that the collected samples were clustered in the main lineage of B.1.617.2 (Delta variant). The Delta variants had a high frequency of distinct mutations in the spike protein region, including T19R (94.12%), L452R (88.23%), T478K (91.17%), D614G (97%), P681R (97%), and D950 N (97%). Other unique mutations found in a smaller frequency in our samples were present in the N-terminal domain, including A27T (2.94%) and A222V (14.70%), and in the receptor-binding domain, including Q414K (5.88%), G446V (2.94%), and T470 N (2.94%).

Conclusion

This study revealed the unique mutations in the S protein region of Delta variants. T19R, L452R, T478K/T478R, D614G, P681R, and D950 N were the most common substitutions in Makassar's Delta variant.

Energetic and frustration analysis of SARS-CoV-2 nucleocapsid protein mutations

The ongoing COVID-19 spreads worldwide with the ability to evolve in diverse human populations. The nucleocapsid (N) protein is one of the mutational hotspots in the SARS-CoV-2 genome. The N protein is an abundant RNA-binding protein critical for viral genome packaging. It comprises two large domains including the N-terminal domain (NTD) and the C-terminal domain (CTD) linked by the centrally located linker region. Mutations in N protein have been reported to increase the severity of disease by modulating viral transmissibility, replication efficiency as well as virulence properties of the virus in different parts of the world. To study the effect of N protein missense mutations on protein stability, function, and pathogenicity, we analyzed 228 mutations from each domain of N protein. Further, we have studied the effect of mutations on local residual frustration changes in N protein. Out of 228 mutations, 11 mutations were predicted to be deleterious and destabilized. Among these mutations, R32C, R191C, and R203 M mutations fall into disordered regions and show significant change in frustration state. Overall, this work reveals that by altering the energetics and residual frustration, N protein mutations might affect the stability, function, and pathogenicity of the SARS-CoV-2.

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.

Expanding repertoire of SARS-CoV-2 deletion mutations contributes to evolution of highly transmissible variants

The emergence of highly transmissible SARS-CoV-2 variants and vaccine breakthrough infections globally mandated the characterization of the immuno-evasive features of SARS-CoV-2. Here, we systematically analyzed 2.13 million SARS-CoV-2 genomes from 188 countries/territories (up to June 2021) and performed whole-genome viral sequencing from 102 COVID-19 patients, including 43 vaccine breakthrough infections. We identified 92 Spike protein mutations that increased in prevalence during at least one surge in SARS-CoV-2 test positivity in any country over a 3-month window. Deletions in the Spike protein N-terminal domain were highly enriched for these 'surge-associated mutations' (Odds Ratio = 14.19, 95% CI 6.15-32.75, p value = 3.41 × 10-10). Based on a longitudinal analysis of mutational prevalence globally, we found an expanding repertoire of Spike protein deletions proximal to an antigenic supersite in the N-terminal domain that may be one of the key contributors to the evolution of highly transmissible variants. Finally, we generated clinically annotated SARS-CoV-2 whole genome sequences from 102 patients and identified 107 unique mutations, including 78 substitutions and 29 deletions. In five patients, we identified distinct deletions between residues 85-90, which reside within a linear B cell epitope. Deletions in this region arose contemporaneously on a diverse background of variants across the globe since December 2020. Overall, our findings based on genomic-epidemiology and clinical surveillance suggest that the genomic deletion of dispensable antigenic regions in SARS-CoV-2 may contribute to the evasion of immune responses and the evolution of highly transmissible variants.

In-depth genetic characterization of the SARS-CoV-2 pandemic in a two-year frame in North Macedonia using second and third generation sequencing technologies

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a persistent negative impact on both the public health and the global economy. To comprehend the origin, transmission routes and discover the mutations that alter the virus’s transmissibility and pathogenicity, full-length SARS-CoV-2 genomes have to be molecularly characterized. Focusing on a two-year time frame (2020-2021), we provide an in-depth virologic and epidemiological overview of the SARS-CoV-2 pandemic in the Republic of North Macedonia by assessing the frequency and distribution of the circulating SARS-CoV-2 variants. Using genetic characterization and phylogenetic analysis we shed light on the molecular evolution of the virus as well as test for a possible connection between specific SARS-CoV-2 haplotypes and the severity of the clinical symptoms. Our results show that one fifth (21.51%) of the tested respiratory samples for SARS-CoV-2 were positive. A noticeable trend in the incidence and severity of the COVID-19 infections was observed in the 60+ age group between males and females. Of the total number of positive cases, the highest incidence of SARS-CoV-2 was noticed in 60+ males (4,170.4/100,000), with a statistically significant (0,0001) difference between the two sexes. Additionally, a 1.8x increase in male mortality and consequentially significantly higher number of death cases was observed compared to females of the same age group (0.001). A total of 327 samples were sequenced in the period March 2020 - August 2021, showing the temporal distribution of SARS-CoV-2 variants circulating in North Macedonia. The phylogenetic analysis showed that most of the viral genomes were closely related and clustered in four distinctive lineages, B.1, B.1.1.7, B.1.351 and B.1.617.2. A statistically significant difference was observed in the 2C_1 haplotype (p=0.0013), where 10.5% of the patients were hospitalized due to severe clinical condition. By employing genetic sequencing, coupled with epidemiological investigations, we investigated viral distribution patterns, identified emerging variants and detected vaccine breakthrough infections. The present work is the first molecular study giving a comprehensive overview of the genetic landscape of circulating SARS-CoV-2 viruses in North Macedonia in a period of two years.

Whole-Genome Comparison of Representatives of All Variants of SARS-CoV-2, Including Subvariant BA.2 and the GKA Clade

Since its discovery at the end of 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has rapidly evolved into many variants, including the subvariant BA.2 and the GKA clade. Genomic clarification is needed for better management of the current pandemic as well as the possible reemergence of novel variants. The sequence of the reference genome Wuhan-Hu-1 and approximately 20 representatives of each variant were downloaded from GenBank and GISAID. Two representatives with no track of in-definitive nucleotides were selected. The sequences were aligned using muscle. The location of insertion/deletion (indel) in the genome was mapped following the open reading frame (ORF) of Wuhan-Hu-1. The phylogeny of the spike protein coding region was constructed using the maximum likelihood method. Amino acid substitutions in all ORFs were analyzed separately. There are two indel sites in ORF1AB, eight in spike, and one each in ORF3A, matrix (MA), nucleoprotein (NP), and the 3'-untranslated regions (3'UTR). Some indel sites and residues/substitutions are not unique, and some are variant-specific. The phylogeny shows that Omicron, Deltacron, and BA2 are clustered together and separated from other variants with 100% bootstrap support. In conclusion, whole-genome comparison of representatives of all variants revealed indel patterns that are specific to SARS-CoV-2 variants or subvariants. Polymorphic amino acid comparison across all coding regions also showed amino acid residues shared by specific groups of variants. Finally, the higher transmissibility of BA.2 might be due at least in part to the 48 nucleotide deletions in the 3'UTR, while the seem-to-be extinction of GKA clade is due to the lack of genetic advantages as a consequence of amino acid substitutions in various genes.

Genome Structure, Life Cycle, and Taxonomy of Coronaviruses and the Evolution of SARS-CoV-2

Coronaviruses have a broad host range and exhibit high zoonotic potential. In this chapter, we describe their genomic organization in terms of encoded proteins and provide an introduction to the peculiar discontinuous transcription mechanism. Further, we present evolutionary conserved genomic RNA secondary structure features, which are involved in the complex replication mechanism. With a focus on computational methods, we review the emergence of SARS-CoV-2 starting with the 2019 strains. In that context, we also discuss the debated hypothesis of whether SARS-CoV-2 was created in a laboratory. We focus on the molecular evolution and the epidemiological dynamics of this recently emerged pathogen and we explain how variants of concern are detected and characterised. COVID-19, the disease caused by SARS-CoV-2, can spread through different transmission routes and also depends on a number of risk factors. We describe how current computational models of viral epidemiology, or more specifically, phylodynamics, have facilitated and will continue to enable a better understanding of the epidemic dynamics of SARS-CoV-2.

Small-Sample Estimation of the Mutational Support and Distribution of SARS-CoV-2

We consider the problem of determining the mutational support and distribution of the SARS-CoV-2 viral genome in the small-sample regime. The mutational support refers to the unknown number of sites that may eventually mutate in the SARS-CoV-2 genome while mutational distribution refers to the distribution of point mutations in the viral genome across a population. The mutational support may be used to assess the virulence of the virus and guide primer selection for real-time RT-PCR testing. Estimating the distribution of mutations in the genome of different subpopulations while accounting for the unseen may also aid in discovering new variants. To estimate the mutational support in the small-sample regime, we use GISAID sequencing data and our state-of-the-art polynomial estimation techniques based on new weighted and regularized Chebyshev approximation methods. For distribution estimation, we adapt the well-known Good-Turing estimator. Our analysis reveals several findings: First, the mutational supports exhibit significant differences in the ORF6 and ORF7a regions (older versus younger patients), ORF1b and ORF10 regions (females versus males) and in almost all ORFs (Asia/Europe/North America). Second, even though the N region of SARS-CoV-2 has a predicted 10% mutational support, mutations fall outside of the primer regions recommended by the CDC.

Mutations of SARS-CoV-2 and their impact on disease diagnosis and severity

Numerous variations of the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2), including D614G, B.1.1.7 (United Kingdom), B.1.1.28 (Brazil P1, P2), CAL.20C (Southern California), B.1.351 (South Africa), B.1.617 (B.1.617.1 Kappa & Delta B.1.617.2) and B.1.1.529, have been reported worldwide. The receptor-binding domain (RBD) of the spike (S) protein is involved in virus-cell binding, where virus-neutralizing antibodies (NAbs) react. Novel variants in the S-protein could maximize viral affinity for the human angiotensin-converting enzyme 2 (ACE2) receptor and increase virus transmission. Molecular detection with false-negative results may refer to mutations in the part of the virus's genome used for virus diagnosis. Furthermore, these changes in S-protein structure alter the neutralizing ability of NAbs, resulting in a reduction in vaccine efficiency. Further information is needed to evaluate how new mutations may affect vaccine efficacy.

Mutations in SARS-CoV-2: Insights on structure, variants, vaccines, and biomedical interventions

COVID-19 is a worldwide pandemic caused by SARS-coronavirus-2 (SARS-CoV-2). Less than a year after the emergence of the Covid-19 pandemic, many vaccines have arrived on the market with innovative technologies in the field of vaccinology. Based on the use of messenger RNA (mRNA) encoding the Spike SARS-Cov-2 protein or on the use of recombinant adenovirus vectors enabling the gene encoding the Spike protein to be introduced into our cells, these strategies make it possible to envisage the vaccination in a new light with tools that are more scalable than the vaccine strategies used so far. Faced with the appearance of new variants, which will gradually take precedence over the strain at the origin of the pandemic, these new strategies will allow a much faster update of vaccines to fight against these new variants, some of which may escape neutralization by vaccine antibodies. However, only a vaccination policy based on rapid and massive vaccination of the population but requiring a supply of sufficient doses could make it possible to combat the emergence of these variants. Indeed, the greater the number of infected individuals, the faster the virus multiplies, with an increased risk of the emergence of variants in these RNA viruses. This review will discuss SARS-CoV-2 pathophysiology and evolution approaches in altered transmission platforms and emphasize the different mutations and how they influence the virus characteristics. Also, this article summarizes the common vaccines and the implication of the mutations and genetic variety of SARS-CoV-2 on the COVID-19 biomedical arbitrations.