Potential Pathogenicity Determinants Identified from Structural Proteomics of SARS-CoV and SARS-CoV-2

"Outlaw" mutations in quasispecies of SARS-CoV-2 inhibit replication

The evolution of SARS-CoV-2, the agent of COVID-19, has been remarkable for its high mutation potential, leading to the appearance of variants. Some mutations have never appeared in the published genomes, which represent consensus, or bona fide genomes. Here we tested the hypothesis that mutations that did not appear in consensus genomes were, in fact, as frequent as the mutations that appeared during the various epidemic episodes, but were not expressed because lethal. To identify these mutations, we analysed the genomes of 90 nasopharyngeal samples and the quasispecies determined by next-generation sequencing. Mutations observed in the quasispecies and not in the consensus genomes were considered to be lethal, what we called "outlaw" mutations. Among these mutations, we analysed the 21 most frequent. Eight of these "outlaws" were in the RNA polymerase and we were able to use a structural biology model and molecular dynamics simulations to demonstrate the functional incapacity of these mutated RNA polymerases. Three other mutations affected the spike, a major protein involved in the pathogenesis of COVID-19. Overall, by analysing the SARS-CoV-2 quasispecies obtained during sequencing, this method made it possible to identify "outlaws," showing areas that could potentially become the target of treatments.

Role of SARS-CoV-2 mutations in the evolution of the COVID-19 pandemic

OBJECTIVES

The SARS-CoV-2 pandemic and large-scale genomic surveillance provided an exceptional opportunity to analyze mutations that appeared over three years in viral genomes. Here we studied mutations and their epidemic consequences for SARS-CoV-2 genomes from our center.

METHODS

We analyzed 61,397 SARS-CoV-2 genomes we sequenced from respiratory samples for genomic surveillance. Mutations frequencies were calculated using Nextclade, Microsoft Excel, and an in-house Python script.

RESULTS

A total of 22,225 nucleotide mutations were identified, 220 (1.0%) being each at the root of ≥836 genomes, classifying mutations as 'hyperfertile'. Two seeded the European pandemic: P323L in RNA polymerase, associated with an increased mutation rate, and D614G in spike that improved fitness. Most 'hyperfertile' mutations occurred in areas not predicted with increased virulence. Their mean number was 8±6 (0-22) per 1000 nucleotides per gene. They were 3.7-times more frequent in accessory than informational genes (13.8 versus 3.7/1000 nucleotides). Particularly, they were 4.1-times more frequent in ORF8 than in the RNA polymerase gene. Interestingly, stop codons were present in 97 positions, almost only in accessory genes, including ORF8 (21/100 codons).

CONCLUSIONS

most 'hyperfertile' mutations did not predict emergence of a new epidemic, and some were stop codons indicating the existence of so-named 'non-virulence' genes.

New-onset and relapsed Graves' disease following COVID-19 vaccination: a comprehensive review of reported cases

Global Coronavir us disease 2019 (COVID-19) vaccination efforts are being intensified to combat the pandemic. As the frequency of immunization against COVID-19 has increased, some adverse effects related to vaccination have emerged. Within this context, this article reviewed 62 Graves' disease (GD) cases following COVID-19 vaccination, to probe the potential association between the vaccination and the onset of GD. A comprehensive search of the PubMed, Web of Science, and Scopus databases was conducted to collect GD cases following COVID-19 vaccination up to June 7, 2023. Among the 62 GD cases included in this review, there were 33 (53.2%) new-onset GD and 10 (16.1%) relapsed GD patients following mRNA vaccination, 14 (22.6%) new-onset GD and 4 (6.5%) relapsed GD patients following viral vector vaccination, and 1 (1.6%) relapsed GD patients following inactivated vaccination. Median durations to symptoms onset for new-onset and relapsed GD were 12 (range: 1-60) and 21 (range: 5-30) days following mRNA vaccination, while 7 (range: 1-28) and 14 (range: 10-14) days following viral vector vaccination, respectively. While the definitive pathogenesis of GD following COVID-19 vaccination remains unclear, it might be associated with cross-immune responses triggered by molecular mimicry, and an adjuvant-induced autoimmune/inflammatory syndrome. However, due to the limited number of observed GD cases following COVID-19 vaccination and the lack of systematic experimental studies, a causal relationship between COVID-19 vaccination and the onset of GD has not been definitively confirmed. It should be highlighted that most of GD patients following COVID-19 vaccination experienced positive outcomes after treatment. In the broader context of ending the COVID-19 pandemic and reducing mortality rates, the benefits of COVID-19 vaccination significantly outweigh mild risks such as treatable GD. Adherence to the COVID-19 vaccination schedule is therefore imperative in effectively managing the pandemic.

Quasispecies Analysis of SARS-CoV-2 of 15 Different Lineages during the First Year of the Pandemic Prompts Scratching under the Surface of Consensus Genome Sequences

The tremendous majority of SARS-CoV-2 genomic data so far neglected intra-host genetic diversity. Here, we studied SARS-CoV-2 quasispecies based on data generated by next-generation sequencing (NGS) of complete genomes. SARS-CoV-2 raw NGS data had been generated for nasopharyngeal samples collected between March 2020 and February 2021 by the Illumina technology on a MiSeq instrument, without prior PCR amplification. To analyze viral quasispecies, we designed and implemented an in-house Excel file (“QuasiS”) that can characterize intra-sample nucleotide diversity along the genomes using data of the mapping of NGS reads. We compared intra-sample genetic diversity and global genetic diversity available from Nextstrain. Hierarchical clustering of all samples based on the intra-sample genetic diversity was performed and visualized with the Morpheus web application. NGS mapping data from 110 SARS-CoV-2-positive respiratory samples characterized by a mean depth of 169 NGS reads/nucleotide position and for which consensus genomes that had been obtained were classified into 15 viral lineages were analyzed. Mean intra-sample nucleotide diversity was 0.21 ± 0.65%, and 5357 positions (17.9%) exhibited significant (>4%) diversity, in ≥2 genomes for 1730 (5.8%) of them. ORF10, spike, and N genes had the highest number of positions exhibiting diversity (0.56%, 0.34%, and 0.24%, respectively). Nine hot spots of intra-sample diversity were identified in the SARS-CoV-2 NSP6, NSP12, ORF8, and N genes. Hierarchical clustering delineated a set of six genomes of different lineages characterized by 920 positions exhibiting intra-sample diversity. In addition, 118 nucleotide positions (0.4%) exhibited diversity at both intra- and inter-patient levels. Overall, the present study illustrates that the SARS-CoV-2 consensus genome sequences are only an incomplete and imperfect representation of the entire viral population infecting a patient, and that quasispecies analysis may allow deciphering more accurately the viral evolutionary pathways.

The origins of COVID-19 pandemic: A brief overview

The novel coronavirus disease (COVID-19) outbreak that emerged at the end of 2019 has now swept the world for more than 2 years, causing immeasurable damage to the lives and economies of the world. It has drawn so much attention to discovering how the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) originated and entered the human body. The current argument revolves around two contradictory theories: a scenario of laboratory spillover events and human contact with zoonotic diseases. Here, we reviewed the transmission, pathogenesis, possible hosts, as well as the genome and protein structure of SARS-CoV-2, which play key roles in the COVID-19 pandemic. We believe the coronavirus was originally transmitted to human by animals rather than by a laboratory leak. However, there still needs more investigations to determine the source of the pandemic. Understanding how COVID-19 emerged is vital to developing global strategies for mitigating future outbreaks.

The emergence, spread and vanishing of a French SARS-CoV-2 variant exemplifies the fate of RNA virus epidemics and obeys the Mistigri rule

The nature and dynamics of mutations associated with the emergence, spread, and vanishing of SARS-CoV-2 variants causing successive waves are complex. We determined the kinetics of the most common French variant ("Marseille-4") for 10 months since its onset in July 2020. Here, we analyzed and classified into subvariants and lineages 7453 genomes obtained by next-generation sequencing. We identified two subvariants, Marseille-4A, which contains 22 different lineages of at least 50 genomes, and Marseille-4B. Their average lifetime was 4.1 ± 1.4 months, during which 4.1 ± 2.6 mutations accumulated. Growth rate was 0.079 ± 0.045, varying from 0.010 to 0.173. Most of the lineages exhibited a bell-shaped distribution. Several beneficial mutations at unpredicted sites initiated a new outbreak, while the accumulation of other mutations resulted in more viral heterogenicity, increased diversity and vanishing of the lineages. Marseille-4B emerged when the other Marseille-4 lineages vanished. Its ORF8 gene was knocked out by a stop codon, as reported in SARS-CoV-2 of mink and in the Alpha variant. This subvariant was associated with increased hospitalization and death rates, suggesting that ORF8 is a nonvirulence gene. We speculate that the observed heterogenicity of a lineage may predict the end of the outbreak.

Structural and functional characterization of NEMO cleavage by SARS-CoV-2 3CLpro

In addition to its essential role in viral polyprotein processing, the SARS-CoV-2 3C-like protease (3CLpro) can cleave human immune signaling proteins, like NF-κB Essential Modulator (NEMO) and deregulate the host immune response. Here, in vitro assays show that SARS-CoV-2 3CLpro cleaves NEMO with fine-tuned efficiency. Analysis of the 2.50 Å resolution crystal structure of 3CLpro C145S bound to NEMO226-234 reveals subsites that tolerate a range of viral and host substrates through main chain hydrogen bonds while also enforcing specificity using side chain hydrogen bonds and hydrophobic contacts. Machine learning- and physics-based computational methods predict that variation in key binding residues of 3CLpro-NEMO helps explain the high fitness of SARS-CoV-2 in humans. We posit that cleavage of NEMO is an important piece of information to be accounted for, in the pathology of COVID-19.

SARS-CoV-2-positive patients display considerable differences in proteome diversity in urine, nasopharyngeal, gargle solution and bronchoalveolar lavage fluid samples

Proteome profile changes post-severe acute respiratory syndrome coronavirus 2 (post-SARS-CoV-2) infection in different body sites of humans remains an active scientific investigation whose solutions stand a chance of providing more information on what constitutes SARS-CoV-2 pathogenesis. While proteomics has been used to understand SARS-CoV-2 pathogenesis, there are limited data about the status of proteome profile in different human body sites infected by the SARS-CoV-2 virus. To bridge this gap, our study aims to characterize the proteins secreted in urine, bronchoalveolar lavage fluid (BALF), gargle solution, and nasopharyngeal samples and assess the proteome differences in these body samples collected from SARS-CoV-2-positive patients. We downloaded publicly available proteomic data from (https://www.ebi.ac.uk/pride/). The data we downloaded had the following identifiers: (i) PXD019423, n = 3 from Charles Tanford Protein Center in Germany. (ii) IPX0002166000, n = 15 from Beijing Proteome Research Centre, China. (iii) IPX0002429000, n = 5 from Huazhong University of Science and Technology, China, and (iv) PXD022889, n = 18 from Department of Laboratory Medicine and Pathology, Mayo Clinic, Rochester, MN 55905 USA. MaxQuant was used for the human peptide spectral matching using human and SARS-CoV-2 proteome database which we downloaded from the UniProt database (access date 13^th October 2021). The individuals infected with SARS-CoV-2 viruses displayed a different proteome diversity from the different body sites we investigated. Overally, we identified 1809 proteins across the four sample types we compared. Urine and BALF samples had significantly more abundant SARS-CoV-2 proteins than the other body sites we compared. Urine samples had 257(33.7%) unique proteins, followed by nasopharyngeal with 250(32.8%) unique proteins. Gargle solution and BALF had 38(5%) and 73(9.6%) unique proteins respectively. Urine, gargle solution, nasopharyngeal, and bronchoalveolar lavage fluid samples have different protein diversity in individuals infected with SARS-CoV-2. Moreover, our data also demonstrated that a given body site is characterized by a unique set of proteins in SARS-CoV-2 seropositive individuals.

Safety and Efficacy of the Common Vaccines against COVID-19

The worldwide pandemic of coronavirus disease 2019 (COVID-19) has imposed a challenge on human health worldwide, and vaccination represents a vital strategy to control the pandemic. To date, multiple COVID-19 vaccines have been granted emergency use authorization, including inactivated vaccines, adenovirus-vectored vaccines, and nucleic acid vaccines. These vaccines have different technical principles, which will necessarily lead to differences in safety and efficacy. Therefore, we aim to implement a systematic review by synthesizing clinical experimental data combined with mass vaccination data and conducting a synthesis to evaluate the safety and efficacy of COVID-19 vaccines. Compared with other vaccines, adverse reactions after vaccination with inactivated vaccines are relatively low. The efficacy of inactivated vaccines is approximately 60%, adenovirus-vectored vaccines are 65%, and mRNA vaccines are 90%, which are always efficient against asymptomatic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, symptomatic COVID-19, COVID-19 hospitalization, severe or critical hospitalization, and death. RNA-based vaccines have a number of advantages and are one of the most promising vaccines identified to date and are particularly important during a pandemic. However, further improvements are required. In time, all the antibody levels weaken gradually, so a booster dose is needed to maintain immunity. Compared with homologous prime-boost immunization, heterologous prime-boost immunization prompts more robust humoral and cellular immune responses.

First cases of infection with the 21L/BA.2 Omicron variant in Marseille, France

The SARS‐CoV‐2 21K/BA.1, 21L/BA.2, and BA.3 Omicron variants have recently emerged worldwide. To date, the 21L/BA.2 Omicron variant has remained very minority globally but became predominant in Denmark instead of the 21K/BA.1 variant. Here, we describe the first cases diagnosed with this variant in south‐eastern France. We identified 13 cases using variant‐specific qPCR and next‐generation sequencing between 28/11/2021 and 31/01/2022, the first two cases being diagnosed in travelers returning from Tanzania. Overall, viral genomes displayed a mean (±standard deviation) number of 65.9 ± 2.5 (range, 61–69) nucleotide substitutions and 31.0 ± 8.3 (27–50) nucleotide deletions, resulting in 49.6 ± 2.2 (45–52) amino acid substitutions (including 28 in the spike protein) and 12.4 ± 1.1 (12–15) amino acid deletions. Phylogeny showed the distribution in three different clusters of these genomes, which were most closely related to genomes from England and South Africa, from Singapore and Nepal, or from France and Denmark. Structural predictions highlighted a significant enlargement and flattening of the surface of the 21L/BA.2 N‐terminal domain of the spike protein compared to that of the 21K/BA.1 Omicron variant, which may facilitate initial viral interactions with lipid rafts. Close surveillance is needed at global, country, and center scales to monitor the incidence and clinical outcome of the 21L/BA.2 Omicron variant.

Antiviral Strategies Against SARS-CoV-2: A Systems Biology Approach

The unprecedented scientific achievements in combating the COVID-19 pandemic reflect a global response informed by unprecedented access to data. We now have the ability to rapidly generate a diversity of information on an emerging pathogen and, by using high-performance computing and a systems biology approach, we can mine this wealth of information to understand the complexities of viral pathogenesis and contagion like never before. These efforts will aid in the development of vaccines, antiviral medications, and inform policymakers and clinicians. Here we detail computational protocols developed as SARS-CoV-2 began to spread across the globe. They include pathogen detection, comparative structural proteomics, evolutionary adaptation analysis via network and artificial intelligence methodologies, and multiomic integration. These protocols constitute a core framework on which to build a systems-level infrastructure that can be quickly brought to bear on future pathogens before they evolve into pandemic proportions.

Mass Spectrometry-Based Proteomics Research to Fight COVID-19: An Expert Review on Hopes and Challenges

The COVID-19 pandemic caused by the severe acute respiratory syndrome (SARS)-CoV-2 infection is a systemic disease and a major planetary health burden. While SARS-CoV-2 impacts host biology extensively, our knowledge of these alterations from a systems perspective remains incomplete. Moreover, there is currently only a limited description of this systemic disease. For precision diagnosis and treatment of SARS-CoV-2, multiomics technologies and systems science research offer significant prospects. This expert review offers a critical analysis of the prospects and challenges of the emerging mass spectrometry-based proteomics approaches to the study of COVID-19 as seen through a systems medicine lens. We also discuss the ways in which proteomics is poised to offer hope for diagnostics and therapeutics innovation on SARS-CoV-2 infection as the disease transitions from a pandemic to an endemic disease, and thus further challenging the health systems and services worldwide in the coming decade. Proteomics is an important high-throughput technology platform to achieve a functional overview of the ways in which COVID-19 changes host biology, and hence, can help identify possible points of entry for innovation in medicines and vaccines, among others.

A comprehensive review of the analysis and integration of omics data for SARS-CoV-2 and COVID-19

Since the first report of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in December 2019, over 100 million people have been infected by COVID-19, millions of whom have died. In the latest year, a large number of omics data have sprung up and helped researchers broadly study the sequence, chemical structure and function of SARS-CoV-2, as well as molecular abnormal mechanisms of COVID-19 patients. Though some successes have been achieved in these areas, it is necessary to analyze and mine omics data for comprehensively understanding SARS-CoV-2 and COVID-19. Hence, we reviewed the current advantages and limitations of the integration of omics data herein. Firstly, we sorted out the sequence resources and database resources of SARS-CoV-2, including protein chemical structure, potential drug information and research literature resources. Next, we collected omics data of the COVID-19 hosts, including genomics, transcriptomics, microbiology and potential drug information data. And subsequently, based on the integration of omics data, we summarized the existing data analysis methods and the related research results of COVID-19 multi-omics data in recent years. Finally, we put forward SARS-CoV-2 (COVID-19) multi-omics data integration research direction and gave a case study to mine deeper for the disease mechanisms of COVID-19.

A vision of immuno-oncology: the Siena think tank of the Italian network for tumor biotherapy (NIBIT) foundation

BACKGROUND

The yearly Think Tank Meeting of the Italian Network for Tumor Biotherapy (NIBIT) Foundation, brings together in Siena, Tuscany (Italy), experts in immuno-oncology to review the learnings from current immunotherapy treatments, and to propose new pre-clinical and clinical investigations in selected research areas. MAIN: While immunotherapies in non-small cell lung cancer and melanoma led to practice changing therapies, the same therapies had only modest benefit for patients with other malignancies, such as mesothelioma and glioblastoma. One way to improve on current immunotherapies is to alter the sequence of each combination agent. Matching the immunotherapy to the host's immune response may thus improve the activity of the current treatments. A second approach is to combine current immunotherapies with novel agents targeting complementary mechanisms. Identifying the appropriate novel agents may require different approaches than the traditional laboratory-based discovery work. For example, artificial intelligence-based research may help focusing the search for innovative and most promising combination partners.

CONCLUSION

Novel immunotherapies are needed in cancer patients with resistance to or relapse after current immunotherapeutic drugs. Such new treatments may include targeted agents or monoclonal antibodies to overcome the immune-suppressive tumor microenvironment. The mode of combining the novel treatments, including vaccines, needs to be matched to the patient's immune status for achieving the maximum benefit. In this scenario, specific attention should be also paid nowadays to the immune intersection between COVID-19 and cancer.

Proteomics-Based Insights Into the SARS-CoV-2-Mediated COVID-19 Pandemic: A Review of the First Year of Research

In late 2019, a virus subsequently named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in China and led to a worldwide pandemic of the disease termed coronavirus disease 2019. The global health threat posed by this pandemic led to an extremely rapid and robust mobilization of the scientific and medical communities as evidenced by the publication of more than 10,000 peer-reviewed articles and thousands of preprints in the first year of the pandemic alone. With the publication of the initial genome sequence of SARS-CoV-2, the proteomics community immediately joined this effort publishing, to date, more than 100 peer-reviewed proteomics studies and submitting many more preprints to preprint servers. In this review, we focus on peer-reviewed articles published on the proteome, glycoproteome, and glycome of SARS-CoV-2. At a basic level, proteomic studies provide valuable information on quantitative aspects of viral infection course; information on the identities, sites, and microheterogeneity of post-translational modifications; and, information on protein-protein interactions. At a biological systems level, these studies elucidate host cell and tissue responses, characterize antibodies and other immune system factors in infection, suggest biomarkers that may be useful for diagnosis and disease-course monitoring, and help in the development or repurposing of potential therapeutics. Here, we summarize results from selected early studies to provide a perspective on the current rapidly evolving literature.

Comprehensive Intrinsic Disorder Analysis of 6108 Viral Proteomes: From the Extent of Intrinsic Disorder Penetrance to Functional Annotation of Disordered Viral Proteins

Much of our understanding of proteins and proteomes comes from the traditional protein structure-function paradigm. However, in the last 2 decades, both computational and experimental studies have provided evidence that a large fraction of functional proteomes across different domains of life consists of intrinsically disordered proteins, thus triggering a quest to unravel and decipher protein intrinsic disorder. Unlike structured/ordered proteins, intrinsically disordered proteins/regions (IDPs/IDRs) do not possess a well-defined structure under physiological conditions and exist as highly dynamic conformational ensembles. In spite of this peculiarity, these proteins have crucial roles in cell signaling and regulation. To date, studies on the abundance and function of IDPs/IDRs in viruses are rather limited. To fill this gap, we carried out an extensive and thorough bioinformatics analysis of 283 000 proteins from 6108 reference viral proteomes. We analyzed protein intrinsic disorder from multiple perspectives, such as abundance of IDPs/IDRs across diverse virus types, their functional annotations, and subcellular localization in taxonomically divergent hosts. We show that the content of IDPs/IDRs in viral proteomes varies broadly as a function of virus genome types and taxonomically divergent hosts. We have combined the two most commonly used and accurate IDP predictors' results with charge-hydropathy (CH) versus cumulative distribution function (CDF) plots to categorize the viral proteins according to their IDR content and physicochemical properties. Mapping of gene ontology on the disorder content of viral proteins reveals that IDPs are primarily involved in key virus-host interactions and host antiviral immune response downregulation, which are reinforced by the post-translational modifications tied to disorder-enriched viral proteins. The present study offers detailed insights into the prevalence of the intrinsic disorder in viral proteomes and provides appealing targets for the design of novel therapeutics.

SARS-CoV-2 main protease suppresses type I interferon production by preventing nuclear translocation of phosphorylated IRF3

Suppression of type I interferon (IFN) response is one pathological outcome of the infection of highly pathogenic human coronaviruses. To effect this, severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2 encode multiple IFN antagonists. In this study, we reported on the IFN antagonism of SARS-CoV-2 main protease NSP5. NSP5 proteins of both SARS-CoV and SARS-CoV-2 counteracted Sendai virus-induced IFN production. NSP5 variants G15S and K90R commonly seen in circulating strains of SARS-CoV-2 retained the IFN-antagonizing property. The suppressive effect of NSP5 on IFN-β gene transcription induced by RIG-I, MAVS, TBK1 and IKKϵ suggested that NSP5 likely acts at a step downstream of IRF3 phosphorylation in the cytoplasm. NSP5 did not influence steady-state expression or phosphorylation of IRF3, suggesting that IRF3, regardless of its phosphorylation state, might not be the substrate of NSP5 protease. However, nuclear translocation of phosphorylated IRF3 was severely compromised in NSP5-expressing cells. Taken together, our work revealed a new mechanism by which NSP5 proteins encoded by SARS-CoV and SARS-CoV-2 antagonize IFN production by retaining phosphorylated IRF3 in the cytoplasm. Our findings have implications in rational design and development of antiviral agents against SARS-CoV-2.

Potentially adaptive SARS-CoV-2 mutations discovered with novel spatiotemporal and explainable AI models

BACKGROUND

A mechanistic understanding of the spread of SARS-CoV-2 and diligent tracking of ongoing mutagenesis are of key importance to plan robust strategies for confining its transmission. Large numbers of available sequences and their dates of transmission provide an unprecedented opportunity to analyze evolutionary adaptation in novel ways. Addition of high-resolution structural information can reveal the functional basis of these processes at the molecular level. Integrated systems biology-directed analyses of these data layers afford valuable insights to build a global understanding of the COVID-19 pandemic.

RESULTS

Here we identify globally distributed haplotypes from 15,789 SARS-CoV-2 genomes and model their success based on their duration, dispersal, and frequency in the host population. Our models identify mutations that are likely compensatory adaptive changes that allowed for rapid expansion of the virus. Functional predictions from structural analyses indicate that, contrary to previous reports, the Asp614Gly mutation in the spike glycoprotein (S) likely reduced transmission and the subsequent Pro323Leu mutation in the RNA-dependent RNA polymerase led to the precipitous spread of the virus. Our model also suggests that two mutations in the nsp13 helicase allowed for the adaptation of the virus to the Pacific Northwest of the USA. Finally, our explainable artificial intelligence algorithm identified a mutational hotspot in the sequence of S that also displays a signature of positive selection and may have implications for tissue or cell-specific expression of the virus.

CONCLUSIONS

These results provide valuable insights for the development of drugs and surveillance strategies to combat the current and future pandemics.