Genomic and proteomic mutation landscapes of SARS-CoV-2

A systematic mutation analysis of 13 major SARS-CoV-2 variants

SARS-CoV-2 evolves constantly with various novel mutations. Due to their enhanced infectivity, transmissibility and immune evasion, a comprehensive understanding of the association between these mutations and the respective functional changes is crucial. However, previous mutation studies of major SARS-CoV-2 variants remain limited. Here, we performed systematic analyses of full-length amino acids mutation, phylogenetic features, protein physicochemical properties, molecular dynamics and immune escape as well as pseudotype virus infection assays among thirteen major SARS-CoV-2 variants. We found that Omicron exhibited the most abundant and complex mutation sites, higher indices of hydrophobicity and flexibility than other variants. The results of molecular dynamics simulation suggest that Omicron has the highest number of hydrogen bonds and strongest binding free energy between the S protein and ACE2 receptor. Furthermore, we revealed 10 immune escape sites in 13 major variants, some of them were reported previously, but four of which (i.e. 339/373/477/496) are first reported to be specific to Omicron, whereas 462 is specific to Epslion. The infectivity of these variants was confirmed by the pseudotype virus infection assays. Our findings may help us understand the functional consequences of the mutations within various variants and the underlying mechanisms of the immune escapes conferred by the S proteins.

SARS-CoV-2 Accessory Protein ORF8 Targets the Dimeric IgA Receptor pIgR

SARS-CoV-2 is a highly pathogenic respiratory virus that successfully initiates and establishes its infection at the respiratory mucosa. However, little is known about how SARS-CoV-2 antagonizes the host's mucosal immunity. Recent findings have shown a marked reduction in the expression of the polymeric Ig receptor (pIgR) in COVID-19 patients. This receptor maintains mucosal homeostasis by transporting the dimeric IgA (dIgA) and pentameric IgM (pIgM) across mucosal epithelial cells to neutralize the invading respiratory pathogens. By studying the interaction between pIgR and SARS-CoV-2 proteins, we discovered that the viral accessory protein Open Reading Frame 8 (ORF8) potently downregulates pIgR expression and that this downregulation activity of ORF8 correlates with its ability to interact with pIgR. Importantly, the ORF8-mediated downregulation of pIgR diminishes the binding of dIgA or pIgM, and the ORF8 proteins of the variants of concern of SARS-CoV-2 preserve the function of downregulating pIgR, indicating the importance of this conserved activity of ORF8 in SARS-CoV-2 pathogenesis. We further observed that the secreted ORF8 binds to cell surface pIgR, but that this interaction does not trigger the cellular internalization of ORF8, which requires the binding of dIgA to pIgR. These findings suggest the role of ORF8 in SARS-CoV-2 mucosal immune evasion.

Development of Biomarkers and Prognosis Model of Mortality Risk in Patients with COVID-19

Background

As of 30 April 2023, the COVID-19 pandemic has resulted in over 6.9 million deaths worldwide. The virus continues to spread and mutate, leading to continuously evolving pathological and physiological processes. It is imperative to reevaluate predictive factors for identifying the risk of early disease progression.

Methods

A retrospective study was conducted on a cohort of 1379 COVID-19 patients who were discharged from Xin Hua Hospital affiliated with Shanghai Jiao Tong University School of Medicine between 15 December 2022 and 15 February 2023. Patient symptoms, comorbidities, demographics, vital signs, and laboratory test results were systematically documented. The dataset was split into testing and training sets, and 15 different machine learning algorithms were employed to construct prediction models. These models were assessed for accuracy and area under the receiver operating characteristic curve (AUROC), and the best-performing model was selected for further analysis.

Results

AUROC for models generated by 15 machine learning algorithms all exceeded 90%, and the accuracy of 10 of them also surpassed 90%. Light Gradient Boosting model emerged as the optimal choice, with accuracy of 0.928 ± 0.0006 and an AUROC of 0.976 ± 0.0028. Notably, the factors with the greatest impact on in-hospital mortality were growth stimulation expressed gene 2 (ST2,19.3%), interleukin-8 (IL-8,17.2%), interleukin-6 (IL-6,6.4%), age (6.1%), NT-proBNP (5.1%), interleukin-2 receptor (IL-2R, 5%), troponin I (TNI,4.6%), congestive heart failure (3.3%) in Light Gradient Boosting model.

Conclusion

ST-2, IL-8, IL-6, NT-proBNP, IL-2R, TNI, age and congestive heart failure were significant predictors of in-hospital mortality among COVID-19 patients.

Vaccination Shortens the Negative Nucleic Acid Conversion Time of the Older Population: A Retrospective Cohort Study of 73,456 Asymptomatic and Mild Patients with COVID-19 in Shanghai

Purpose

To explore the influencing factors to predict the negative nucleic acid conversion time and ORF1ab gene CT value changes in patients with asymptomatic and mild COVID-19.

Patients and Methods

A total of 73,456 patients with asymptomatic and mild COVID-19 admitted to the Mobile Cabin Hospital in Shanghai from April 3 to April 23, 2022 were selected as the research objects. Epidemiological, clinical, and laboratory data were collected. Correlation analysis was performed.

Results

In patients <18 years old and ≥65 years old, COVID-19 vaccine can shorten the negative nucleic acid conversion time, which is reflected in the lower median or 75% quantile (P<0.001, P<0.05). In patients with underlying diseases, the negative nucleic acid conversion time of booster vaccination and complete vaccination was lower than that of non-vaccinated group (P<0.001, P<0.05). In patients ≤18 years of age or >65 years of age, patients with comorbidity and patients with symptoms, compared with patients 18-65 years of age, patients without comorbidity and patients without symptoms, there was a greater difference in the rate of rise of CT values between vaccinated and unvaccinated patients (P<0.05).

Conclusion

The time of nucleic acid conversion to negative in patients with asymptomatic and mild COVID-19 is affected by age, comorbidity, and first nucleic acid CT value. Vaccination could shorten the negative nucleic acid conversion time of the older population, those with complications or symptoms. The vaccination of older patients does not increase the risk of symptoms.

The RdRp genotyping of SARS-CoV-2 isolated from patients with different clinical spectrum of COVID-19

BACKGROUND

The evolution of SARS-CoV-2 has been observed from the very beginning of the fight against COVID-19, some mutations are indicators of potentially dangerous variants of the virus. However, there is no clear association between the genetic variants of SARS-CoV-2 and the severity of COVID-19. We aimed to analyze the genetic variability of RdRp in correlation with different courses of COVID-19.

RESULTS

The prospective study included 77 samples of SARS-CoV-2 isolated from outpatients (1st degree of severity) and hospitalized patients (2nd, 3rd and 4th degree of severity). The retrospective analyses included 15,898,266 cases of SARS-CoV-2 genome sequences deposited in the GISAID repository. Single-nucleotide variants were identified based on the four sequenced amplified fragments of SARS-CoV-2. The analysis of the results was performed using appropriate statistical methods, with p < 0.05, considered statistically significant. Additionally, logistic regression analysis was performed to predict the strongest determinants of the observed relationships. The number of mutations was positively correlated with the severity of the COVID-19, and older male patients. We detected four mutations that significantly increased the risk of hospitalization of COVID-19 patients (14676C > T, 14697C > T, 15096 T > C, and 15279C > T), while the 15240C > T mutation was common among strains isolated from outpatients. The selected mutations were searched worldwide in the GISAID database, their presence was correlated with the severity of COVID-19.

CONCLUSION

Identified mutations have the potential to be used to assess the increased risk of hospitalization in COVID-19 positive patients. Experimental studies and extensive epidemiological data are needed to investigate the association between individual mutations and the severity of COVID-19.

RNA barcode segments for SARS-CoV-2 identification from HCoVs and SARSr-CoV-2 lineages

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the pathogen responsible for coronavirus disease 2019 (COVID-19), continues to evolve, giving rise to more variants and global reinfections. Previous research has demonstrated that barcode segments can effectively and cost-efficiently identify specific species within closely related populations. In this study, we designed and tested RNA barcode segments based on genetic evolutionary relationships to facilitate the efficient and accurate identification of SARS-CoV-2 from extensive virus samples, including human coronaviruses (HCoVs) and SARSr-CoV-2 lineages. Nucleotide sequences sourced from NCBI and GISAID were meticulously selected and curated to construct training sets, encompassing 1733 complete genome sequences of HCoVs and SARSr-CoV-2 lineages. Through genetic-level species testing, we validated the accuracy and reliability of the barcode segments for identifying SARS-CoV-2. Subsequently, 75 main and subordinate species-specific barcode segments for SARS-CoV-2, located in ORF1ab, S, E, ORF7a, and N coding sequences, were intercepted and screened based on single-nucleotide polymorphism sites and weighted scores. Post-testing, these segments exhibited high recall rates (nearly 100%), specificity (almost 30% at the nucleotide level), and precision (100%) performance on identification. They were eventually visualized using one and two-dimensional combined barcodes and deposited in an online database (http://virusbarcodedatabase.top/). The successful integration of barcoding technology in SARS-CoV-2 identification provides valuable insights for future studies involving complete genome sequence polymorphism analysis. Moreover, this cost-effective and efficient identification approach also provides valuable reference for future research endeavors related to virus surveillance.

SARS-CoV-2 ORF8 as a Modulator of Cytokine Induction: Evidence and Search for Molecular Mechanisms

Severe cases of SARS-CoV-2 infection are characterized by an immune response that leads to the overproduction of pro-inflammatory cytokines, resulting in lung damage, cardiovascular symptoms, hematologic symptoms, acute kidney injury and multiple organ failure that can lead to death. This remarkable increase in cytokines and other inflammatory molecules is primarily caused by viral proteins, and particular interest has been given to ORF8, a unique accessory protein specific to SARS-CoV-2. Despite plenty of research, the precise mechanisms by which ORF8 induces proinflammatory cytokines are not clear. Our investigations demonstrated that ORF8 augments production of IL-6 induced by Poly(I:C) in human embryonic kidney (HEK)-293 and monocyte-derived dendritic cells (mono-DCs). We discuss our findings and the multifaceted roles of ORF8 as a modulator of cytokine response, focusing on type I interferon and IL-6, a key component of the immune response to SARS-CoV-2. In addition, we explore the hypothesis that ORF8 may act through pattern recognition receptors of dsRNA such as TLRs.

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

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

Comparison of infection and human immune responses of two SARS-CoV-2 strains in a humanized hACE2 NIKO mouse model

The COVID-19 pandemic has sickened millions, cost lives and has devastated the global economy. Various animal models for experimental infection with SARS-CoV-2 have played a key role in many aspects of COVID-19 research. Here, we describe a humanized hACE2 (adenovirus expressing hACE2) NOD-SCID IL2Rγ^-/- (NIKO) mouse model and compare infection with ancestral and mutant (SARS-CoV-2-∆382) strains of SARS-CoV-2. Immune cell infiltration, inflammation, lung damage and pro-inflammatory cytokines and chemokines was observed in humanized hACE2 NIKO mice. Humanized hACE2 NIKO mice infected with the ancestral and mutant SARS-CoV-2 strain had lung inflammation and production of pro-inflammatory cytokines and chemokines. This model can aid in examining the pathological basis of SARS-CoV-2 infection in a human immune environment and evaluation of therapeutic interventions.

The End or a New Era of Development of SARS-CoV-2 Virus: Genetic Variants Responsible for Severe COVID-19 and Clinical Efficacy of the Most Commonly Used Vaccines in Clinical Practice

Although the chief of the World Health Organization (WHO) has declared the end of the coronavirus disease 2019 (COVID-19) as a global health emergency, the disease is still a global threat. To be able to manage such pandemics in the future, it is necessary to develop proper strategies and opportunities to protect human life. The data on the SARS-CoV-2 virus must be continuously analyzed, and the possibilities of mutation and the emergence of new, more infectious variants must be anticipated, as well as the options of using different preventive and therapeutic techniques. This is because the fast development of severe acute coronavirus 2 syndrome (SARS-CoV-2) variants of concern have posed a significant problem for COVID-19 pandemic control using the presently available vaccinations. This review summarizes data on the SARS-CoV-2 variants that are responsible for severe COVID-19 and the clinical efficacy of the most commonly used vaccines in clinical practice. The consequences after the disease (long COVID or post-COVID conditions) continue to be the subject of studies and research, and affect social and economic life worldwide.

Immunogenic cell death-led discovery of COVID-19 biomarkers and inflammatory infiltrates

Immunogenic cell death (ICD) serves a critical role in regulating cell death adequate to activate an adaptive immune response, and it is associated with various inflammation-related diseases. However, the specific role of ICD-related genes in COVID-19 remains unclear. We acquired COVID-19-related information from the GEO database and a total of 14 ICD-related differentially expressed genes (DEGs) were identified. These ICD-related DEGs were closely associated with inflammation and immune activity. Afterward, CASP1, CD4, and EIF2AK3 among the 14 DEGs were selected as feature genes based on LASSO, Random Forest, and SVM-RFE algorithms, which had reliable diagnostic abilities. Moreover, functional enrichment analysis indicated that these feature genes may have a potential role in COVID-19 by being involved in the regulation of immune response and metabolism. Further CIBERSORT analysis demonstrated that the variations in the immune microenvironment of COVID-19 patients may be correlated with CASP1, CD4, and EIF2AK3. Additionally, 33 drugs targeting 3 feature genes had been identified, and the ceRNA network demonstrated a complicated regulative association based on these feature genes. Our work identified that CASP1, CD4, and EIF2AK3 were diagnostic genes of COVID-19 and correlated with immune activity. This study presents a reliable diagnostic signature and offers an overview to investigate the mechanism of COVID-19.

SARS-CoV-2 ORF8: A Rapidly Evolving Immune and Viral Modulator in COVID-19

The COVID-19 pandemic has resulted in upwards of 6.8 million deaths over the past three years, and the frequent emergence of variants continues to strain global health. Although vaccines have greatly helped mitigate disease severity, SARS-CoV-2 is likely to remain endemic, making it critical to understand its viral mechanisms contributing to pathogenesis and discover new antiviral therapeutics. To efficiently infect, this virus uses a diverse set of strategies to evade host immunity, accounting for its high pathogenicity and rapid spread throughout the COVID-19 pandemic. Behind some of these critical host evasion strategies is the accessory protein Open Reading Frame 8 (ORF8), which has gained recognition in SARS-CoV-2 pathogenesis due to its hypervariability, secretory property, and unique structure. This review discusses the current knowledge on SARS-CoV-2 ORF8 and proposes actualized functional models describing its pivotal roles in both viral replication and immune evasion. A better understanding of ORF8's interactions with host and viral factors is expected to reveal essential pathogenic strategies utilized by SARS-CoV-2 and inspire the development of novel therapeutics to improve COVID-19 disease outcomes.

Variation in synonymous evolutionary rates in the SARS-CoV-2 genome

Introduction

Coronavirus disease 2019 is an infectious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Influential variants and mutants of this virus continue to emerge, and more effective virus-related information is urgently required for identifying and predicting new mutants. According to earlier reports, synonymous substitutions were considered phenotypically silent; thus, such mutations were frequently ignored in studies of viral mutations because they did not directly cause amino acid changes. However, recent studies have shown that synonymous substitutions are not completely silent, and their patterns and potential functional correlations should thus be delineated for better control of the pandemic.

Methods

In this study, we estimated the synonymous evolutionary rate (SER) across the SARS-CoV-2 genome and used it to infer the relationship between the viral RNA and host protein. We also assessed the patterns of characteristic mutations found in different viral lineages.

Results

We found that the SER varies across the genome and that the variation is primarily influenced by codon-related factors. Moreover, the conserved motifs identified based on the SER were found to be related to host RNA transport and regulation. Importantly, the majority of the existing fixed-characteristic mutations for five important virus lineages (Alpha, Beta, Gamma, Delta, and Omicron) were significantly enriched in partially constrained regions.

Discussion

Taken together, our results provide unique information on the evolutionary and functional dynamics of SARS-CoV-2 based on synonymous mutations and offer potentially useful information for better control of the SARS-CoV-2 pandemic.

Effects of natural polymorphisms in SARS-CoV-2 RNA-dependent RNA polymerase on its activity and sensitivity to inhibitors in vitro

SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) is the key enzyme required for viral replication and mRNA synthesis. RdRp is one of the most conserved viral proteins and a promising target for antiviral drugs and inhibitors. At the same time, analysis of public databases reveals multiple variants of SARS-CoV-2 genomes with substitutions in the catalytic RdRp subunit nsp12. Structural mapping of these mutations suggests that some of them may affect the interactions of nsp12 with its cofactors nsp7/nsp8 as well as with RNA substrates. We have obtained several mutations of these types and demonstrated that some of them decrease specific activity of RdRp in vitro, possibly by changing RdRp assembly and/or its interactions with RNA. Therefore, natural polymorphisms in RdRp may potentially affect viral replication. Furthermore, we have synthesized a series of polyphenol and diketoacid derivatives based on previously studied inhibitors of hepatitis C virus RdRp and found that several of them can inhibit SARS-CoV-2 RdRp. Tested mutations in RdRp do not have strong effects on the efficiency of inhibition. Further development of more efficient non-nucleoside inhibitors of SARS-CoV-2 RdRp should take into account the existence of multiple polymorphic variants of RdRp.

Characterization of SARS-CoV-2 isolated from a patient in Iran compared to SARS-CoV-2 different variants

Since the onset of the global epidemic of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), whole genome sequencing of virus in all countries has been considered to track and predict virus transmission and variation patterns. In the current study we reported a novel complete genome sequence of SARS-CoV-2 isolated from Iran. Genomics variations and protein sequences were evaluated for the isolated sequence and seven Iranian complete genome sequences of SARS-CoV-2 from NCBI using the reference genome of the SARS-CoV-2 Wuhan-Hu-1. The results showed six nucleotide substitutions. The multiple sequence alignment of the spike protein of the Wuhan-Hu-1 strain and the emerging variants indicated similar its residue pattern in the current sequence to the Wuhan-Hu-1 strain. There were relatively similar binding affinity and residues involved in the interactions of the spike receptor-binding domain (RBD) of the Wuhan-Hu-1 strain, the variants and Hormozgan With angiotensin-converting enzyme 2 (ACE2). Tracing the phylogeny of virus indicated distinct clustering of Iranian variants in branches close to the Asian countries. The mutation effect study on the function of proteins predicted neutral impact of all six nucleotide substitutions. However, the free energy calculations indicated a decreasing the protein stability related to the mutations. This data, consistent with similar studies, showed that despite the high similarity in the nucleotide sequence of the SARS-CoV-2, the mutation pattern varies from country to country. Therefore, any country can benefit from these studies to track and find appropriate strategies for treating and controlling the epidemic.Communicated by Ramaswamy H. Sarma.

SARS-CoV-2 Mutations Lead to a Decrease in the Number of Tissue-Specific MicroRNA-Binding Regions in the Lung

RNA interference in vertebrates acts as an antiviral mechanism only in undifferentiated embryonic stem cells and is mediated by microRNAs. In somatic cells, host microRNAs also bind to the genomes of RNA viruses, regulating their translation and replication. It has been shown that viral (+)RNA can evolve under the influence of host cell miRNAs. In more than two years of the pandemic, the SARS-CoV-2 virus has mutated significantly. It is quite possible that some mutations could be retained in the virus genome under the influence of miRNAs produced by alveolar cells. We demonstrated that microRNAs in human lung tissue exert evolutionary pressure on the SARS-CoV-2 genome. Moreover, a significant number of sites of host microRNA binding with the virus genome are located in the NSP3-NSP5 region responsible for autoproteolysis of viral polypeptides.

The Spread of SARS-CoV-2 Omicron Variant in CALABRIA: A Spatio-Temporal Report of Viral Genome Evolution

We investigated the evolution of SARS-CoV-2 spread in Calabria, Southern Italy, in 2022. A total of 272 RNA isolates from nasopharyngeal swabs of individuals infected with SARS-CoV-2 were sequenced by whole genome sequencing (N = 172) and/or Sanger sequencing (N = 100). Analysis of diffusion of Omicron variants in Calabria revealed the prevalence of 10 different sub-lineages (recombinant BA.1/BA.2, BA.1, BA.1.1, BA.2, BA.2.9, BA.2.10, BA.2.12.1, BA.4, BA.5, BE.1). We observed that Omicron spread in Calabria presented a similar trend as in Italy, with some notable exceptions: BA.1 disappeared in April in Calabria but not in the rest of Italy; recombinant BA.1/BA.2 showed higher frequency in Calabria (13%) than in the rest of Italy (0.02%); BA.2.9, BA.4 and BA.5 emerged in Calabria later than in other Italian regions. In addition, Calabria Omicron presented 16 non-canonical mutations in the S protein and 151 non-canonical mutations in non-structural proteins. Most non-canonical mutations in the S protein occurred mainly in BA.5 whereas non-canonical mutations in non-structural or accessory proteins (ORF1ab, ORF3a, ORF8 and N) were identified in BA.2 and BA.5 sub-lineages. In conclusion, the data reported here underscore the importance of monitoring the entire SARS-CoV-2 genome.

HLA-A*01:01 allele diminishing in COVID-19 patients population associated with non-structural epitope abundance in CD8+ T-cell repertoire

In mid-2021, the SARS-CoV-2 Delta variant caused the third wave of the COVID-19 pandemic in several countries worldwide. The pivotal studies were aimed at studying changes in the efficiency of neutralizing antibodies to the spike protein. However, much less attention was paid to the T-cell response and the presentation of virus peptides by MHC-I molecules. In this study, we compared the features of the HLA-I genotype in symptomatic patients with COVID-19 in the first and third waves of the pandemic. As a result, we could identify the diminishing of carriers of the HLA-A*01:01 allele in the third wave and demonstrate the unique properties of this allele. Thus, HLA-A*01:01-binding immunoprevalent epitopes are mostly derived from ORF1ab. A set of epitopes from ORF1ab was tested, and their high immunogenicity was confirmed. Moreover, analysis of the results of single-cell phenotyping of T-cells in recovered patients showed that the predominant phenotype in HLA-A*01:01 carriers is central memory T-cells. The predominance of T-lymphocytes of this phenotype may contribute to forming long-term T-cell immunity in carriers of this allele. Our results can be the basis for highly effective vaccines based on ORF1ab peptides.

Protein post-translational modification in SARS-CoV-2 and host interaction

SARS-CoV-2 can cause lung diseases, such as pneumonia and acute respiratory distress syndrome, and multi-system dysfunction. Post-translational modifications (PTMs) related to SARS-CoV-2 are conservative and pathogenic, and the common PTMs are glycosylation, phosphorylation, and acylation. The glycosylation of SARS-CoV-2 mainly occurs on spike (S) protein, which mediates the entry of the virus into cells through interaction with angiotensin-converting enzyme 2. SARS-CoV-2 utilizes glycans to cover its epitopes and evade the immune response through glycosylation of S protein. Phosphorylation of SARS-CoV-2 nucleocapsid (N) protein improves its selective binding to viral RNA and promotes viral replication and transcription, thereby increasing the load of the virus in the host. Succinylated N and membrane(M) proteins of SARS-CoV-2 synergistically affect virus particle assembly. N protein regulates its affinity for other proteins and the viral genome through acetylation. The acetylated envelope (E) protein of SARS-CoV-2 interacts with bromodomain-containing protein 2/4 to influence the host immune response. Both palmitoylation and myristoylation sites on S protein can affect the virus infectivity. Papain-like protease is a domain of NSP3 that dysregulates host inflammation by deubiquitination and impinges host IFN-I antiviral immune responses by deISGylation. Ubiquitination of ORF7a inhibits host IFN-α signaling by blocking STAT2 phosphorylation. The methylation of N protein can inhibit the formation of host stress granules and promote the binding of N protein to viral RNA, thereby promoting the production of virus particles. NSP3 macrodomain can reverse the ADP-ribosylation of host proteins, and inhibit the cascade immune response with IFN as the core, thereby promoting the intracellular replication of SARS-CoV-2. On the whole, PTMs have fundamental roles in virus entry, replication, particle assembly, and host immune response. Mutations in various SARS-CoV-2 variants, which lead to changes in PTMs at corresponding sites, cause different biological effects. In this paper, we mainly reviewed the effects of PTMs on SARS-CoV-2 and host cells, whose application is to inform the strategies for inhibiting viral infection and facilitating antiviral treatment and vaccine development for COVID-19.

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.

Mutations in SARS-CoV-2 structural proteins: a global analysis

BACKGROUND

Emergence of new variants mainly variants of concerns (VOC) is caused by mutations in main structural proteins of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Therefore, we aimed to investigate the mutations among structural proteins of SARS-CoV-2 globally.

METHODS

We analyzed samples of amino-acid sequences (AASs) for envelope (E), membrane (M), nucleocapsid (N), and spike (S) proteins from the declaration of the coronavirus 2019 (COVID-19) as pandemic to January 2022. The presence and location of mutations were then investigated by aligning the sequences to the reference sequence and categorizing them based on frequency and continent. Finally, the related human genes with the viral structural genes were discovered, and their interactions were reported.

RESULTS

The results indicated that the most relative mutations among the E, M, N, and S AASs occurred in the regions of 7 to 14, 66 to 88, 164 to 205, and 508 to 635 AAs, respectively. The most frequent mutations in E, M, N, and S proteins were T9I, I82T, R203M/R203K, and D614G. D614G was the most frequent mutation in all six geographical areas. Following D614G, L18F, A222V, E484K, and N501Y, respectively, were ranked as the most frequent mutations in S protein globally. Besides, A-kinase Anchoring Protein 8 Like (AKAP8L) was shown as the linkage unit between M, E, and E cluster genes.

CONCLUSION

Screening the structural protein mutations can help scientists introduce better drug and vaccine development strategies.

Progress of cGAS-STING signaling in response to SARS-CoV-2 infection

Coronavirus disease 2019 (COVID-19) is an epidemic respiratory disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that can cause infections in millions of individuals, who can develop lung injury, organ failure, and subsequent death. As the first line of host defense, the innate immune system is involved in initiating the immune response to SARS-CoV-2 infection and the hyperinflammatory phenotype of COVID-19. However, the interplay between SARS-CoV-2 and host innate immunity is not yet well understood. It had become known that the cGAS-STING pathway is involved in the detection of cytosolic DNA, which elicits an innate immune response involving a robust type I interferon response against viral and bacterial infections. Nevertheless, several lines of evidence indicate that SARS-CoV-2, a single-stranded positive-sense RNA virus, triggered the cGAS-STING signaling pathway. Therefore, understanding the molecular and cellular details of cGAS-STING signaling upon SARS-CoV-2 infection is of considerable biomedical importance. In this review, we discuss the role of cGAS-STING signaling in SARS-CoV-2 infection and summarize the potential therapeutics of STING agonists as virus vaccine adjuvants.

Temporal changes in the accessory protein mutations of SARS-CoV-2 variants and their predicted structural and functional effects

Emerging variants enable the continuous spread of SARS-CoV-2 in humans. The factors contributing to behavioral differences in variants remain elusive despite associations with several Spike protein mutations. Exploring accessory proteins may provide a wider understanding of these differences since these proteins may affect viral processes that occur beyond infection. Various bioinformatics tools were utilized to identify significant accessory protein mutations and determine their structural and functional effects over time. The ViruClust web application was used to retrieve accessory protein amino acid sequences and determine mutation frequencies in these sequences across time. The structural and functional effects of the mutations were determined using Missense3D and PROVEAN, respectively. The accessory and Spike protein mutations were compared using mutation densities. Q57H and T151I of ORF3a; T21I and W27L of ORF6; G38V, V82A, and T120I of ORF7a; S31P and T40I of ORF7b; and R52I, C61F, and I121L of ORF8 were highly frequent in most variants of concern and were within known functional domains. Thus, these are good candidates for further experimental evaluation. Among the accessory proteins, ORF6 and ORF8 were highlighted because of their strong and weak correlation with Spike protein mutations, respectively.

Analysis of the mutation dynamics of SARS-CoV-2 genome in the samples from Georgia State of the United States

BACKGROUND

The COVID-19 is caused by a novel coronavirus SARS-CoV-2, which started from China. It spread rapidly throughout the world and was later declared a pandemic by the WHO. Over the course of time, SARS-CoV-2 has mutated for survival advantages, and this led to multiple variants. Multiple studies on mutations identification in SARS-CoV2 have been published covering extensive sample areas. The purpose of this study was to limit the sample area to the Georgia state in the U.S. and to analyze the genome sequences for mutation profiling across the genome and origin of variants.

METHODS

The genome sequences (n = 3,970) were obtained from the NCBI database as of June 12, 2021, with the filter of being complete sequenced genomes, homo-sapiens host, and only from Georgia State of the U.S. NextClade, an online tool was used for the analysis of the sequences using Wuhan-Hu-1/2019 as a reference genome. The algorithm was sequence alignment, translation, mutation calling, phylogenetic placement, clade assignment, and quality control (QC). Thirty-six samples with bad QC were removed from the mutational analysis.

RESULTS

A total 117,743 mutations in the nucleotides were identified (averaging 31.5 mutations per sample). The mutations A23403G, C3037T, C241T, and C14408T were detected in 98% of the samples. Also, a total of 75,517 mutations in the amino acid were identified (averaging 20.2 mutations per sample). The mutations D614G and P314L were identified in >97% samples whereas R203K, G204R, P681H, and N501Y were detected in >50% samples. Analysis also revealed 16 different clades with 20I (49.6%). Clades 20G (24.2%) and 20A (5.5%) being the most abundant, showed that SARS-CoV-2 in the Georgia State originated mainly from Southeast England, other parts of the U.S., and several countries in Western Europe.

CONCLUSION

Looking at the three most common variants in Georgia State of the U.S., we could determine the primary locations of transmission or origin for the virus, and our analyses indicates that majority of the cases originated from Southeast England (Clade 20I), the U.S. itself (Clade 20G), and from Western Europe (Clade 20C).

Immunological study of COVID-19 vaccine candidate based on recombinant spike trimer protein from different SARS-CoV-2 variants of concern

The emergency of new SARS-CoV-2 variants that feature increased immune escape marks an urgent demand for better vaccines that will provide broader immunogenicity. Here, we evaluated the immunogenic capacity of vaccine candidates based on the recombinant trimeric spike protein (S) of different SARS-CoV-2 variants of concern (VOC), including the ancestral Wuhan, Beta and Delta viruses. In particular, we assessed formulations containing either single or combined S protein variants. Our study shows that the formulation containing the single S protein from the ancestral Wuhan virus at a concentration of 2µg (SW2-Vac 2µg) displayed in the mouse model the highest IgG antibody levels against all the three (Wuhan, Beta, and Delta) SARS-CoV-2 S protein variants tested. In addition, this formulation induced significantly higher neutralizing antibody titers against the three viral variants when compared with authorized Gam-COVID-Vac-rAd26/rAd5 (Sputnik V) or ChAdOx1 (AstraZeneca) vaccines. SW2-Vac 2µg was also able to induce IFN-gamma and IL-17, memory CD4 populations and follicular T cells. Used as a booster dose for schedules performed with different authorized vaccines, SW2-Vac 2µg vaccine candidate also induced higher levels of total IgG and IgG isotypes against S protein from different SARS-CoV-2 variants in comparison with those observed with homologous 3-dose schedule of Sputnik V or AstraZeneca. Moreover, SW2-Vac 2µg booster induced broadly strong neutralizing antibody levels against the three tested SARS-CoV-2 variants. SW2-Vac 2µg booster also induced CD4+ central memory, CD4+ effector and CD8+ populations. Overall, the results demonstrate that SW2-Vac 2 µg is a promising formulation for the development of a next generation COVID-19 vaccine.

Incorporating MOOC and COVID-19-Related Scientific Papers into Veterinary Microbiology Teaching to Enhance Students' Learning Performance and Professional Recognition

The COVID-19 pandemic has exerted a huge adverse influence on global teaching activities and students' psychological status. Veterinary microbiology is mainly concerned with bacterial and viral diseases, including coronavirus diseases. An innovative online-to-offline teaching approach for this course was established to stimulate students' learning initiative and mitigate their anxiety about COVID-19. A well-established massive open online course (MOOC) was first adopted as preview material before class, followed by in-person teaching. Additionally, COVID-19-related scientific papers were also used as pre-class reading material in veterinary microbiology and were further explained in class. The effect of this innovative teaching mode was systematically evaluated by final examination scores and questionnaires. The average score (81.75) and excellence score rating (> 85 scores, 37.3%) resulting from this blended teaching mode were not statistically higher than those of the online-only (79.19, p = .115; 28.6%, p = .317) or offline-only (79.47, p = .151; 27.9%, p = .269) teaching modes. This may be due to the sample size investigated; however, the results indicate that the innovative teaching mode did not decrease teaching quality. Additionally, most subjects (72.9%) were satisfied with the blended mode and supported its future use. Intriguingly, the introduction of COVID-19-related scientific papers helped students understand virology, relieve their anxiety, and increase their professional identity. Collectively, the innovative approach to teaching veterinary microbiology in this study provides a beneficial reference for other teachers to maintain and improve teaching quality.

Second round of an interlaboratory comparison of SARS-CoV2 molecular detection assays used by 45 veterinary diagnostic laboratories in the United States

The COVID-19 pandemic presents a continued public health challenge. Veterinary diagnostic laboratories in the United States use RT-rtPCR for animal testing, and many laboratories are certified for testing human samples; hence, ensuring that laboratories have sensitive and specific SARS-CoV2 testing methods is a critical component of the pandemic response. In 2020, the FDA Veterinary Laboratory Investigation and Response Network (Vet-LIRN) led an interlaboratory comparison (ILC1) to help laboratories evaluate their existing RT-rtPCR methods for detecting SARS-CoV2. All participating laboratories were able to detect the viral RNA spiked in buffer and PrimeStore molecular transport medium (MTM). With ILC2, Vet-LIRN extended ILC1 by evaluating analytical sensitivity and specificity of the methods used by participating laboratories to detect 3 SARS-CoV2 variants (B.1; B.1.1.7 [Alpha]; B.1.351 [Beta]) at various copy levels. We analyzed 57 sets of results from 45 laboratories qualitatively and quantitatively according to the principles of ISO 16140-2:2016. More than 95% of analysts detected the SARS-CoV2 RNA in MTM at ≥500 copies for all 3 variants. In addition, for nucleocapsid markers N1 and N2, 81% and 92% of the analysts detected ≤20 copies in the assays, respectively. The analytical specificity of the evaluated methods was >99%. Participating laboratories were able to assess their current method performance, identify possible limitations, and recognize method strengths as part of a continuous learning environment to support the critical need for the reliable diagnosis of COVID-19 in potentially infected animals and humans.

Emphasis on heat strain to the ocular surface: A functional and clinical study of a modified goggle

Purpose

The limitations of conventional goggles have caused immense inconvenience, and even damage, to the physical and mental health of healthcare workers. Hence, this study aimed to build a modified goggle (MG) with better physical performance. The temperature-humidity index (THI) was used as an indicator to investigate the impact of goggle-related heat strain on the ocular surface.

Methods

The basic functions of antifog, anti-ultraviolet (UV), and anti-blue-light radiation capabilities were evaluated. Furthermore, the clinical impact on noninvasive keratography tear film break-up time (NIKBUT), intraocular pressure, central corneal thickness, Schirmer test I, and the Dry Eye-related Quality of life Score (DEQS) were assessed in 40 healthcare workers by comparing MG with standard goggles (SG). The relationships between THI and the above parameters were explored.

Results

MG had a significantly longer antifog time than SG (212.75 ± 23.95 vs. 138.35 ± 5.54 min, p < 0.05), stronger antiultraviolet ability at 400 nm (99.99 vs. 45.55%), and optimal anti-blue-light performance at 440 nm (33.32 vs. 13.31%). Tear film stability after wearing the goggle was significantly worse than that before wearing them (p < 0.05). Both goggles achieved moderate to strong heat strain, with a THI of >80 at all timepoints. The MG group showed lower THI and DEQS and higher NIKBUT than the SG group (p < 0.05). THI was significantly correlated with DEQS, NIKBUT, and real fogging time (r = 0.876, −0.532, −0.406; all p < 0.05).

Conclusion

Wearing goggles for a long time may cause heat strain to the eyes, thereby leading to eye discomfort and changes in the microenvironment of the ocular surface. Our MG exhibited better antifog, antiultraviolet, and optimal anti-blue-light performance and lower heat strain than SG, thus making it ideally suited for healthcare workers.

Small-Molecule Thioesters as SARS-CoV-2 Main Protease Inhibitors: Enzyme Inhibition, Structure-Activity Relationships, Antiviral Activity, and X-ray Structure Determination

The main protease (Mpro, 3CLpro) of SARS-CoV-2 is an attractive target in coronaviruses because of its crucial involvement in viral replication and transcription. Here, we report on the design, synthesis, and structure-activity relationships of novel small-molecule thioesters as SARS-CoV-2 Mpro inhibitors. Compounds 3w and 3x exhibited excellent SARS-CoV-2 Mpro inhibition with kinac/Ki of 58,700 M-1 s-1 (Ki = 0.0141 μM) and 27,200 M-1 s-1 (Ki = 0.0332 μM), respectively. In Calu-3 and Vero76 cells, compounds 3h, 3i, 3l, 3r, 3v, 3w, and 3x displayed antiviral activity in the nanomolar range without host cell toxicity. Co-crystallization of 3w and 3af with SARS-CoV-2 Mpro was accomplished, and the X-ray structures showed covalent binding with the catalytic Cys145 residue of the protease. The potent SARS-CoV-2 Mpro inhibitors also inhibited the Mpro of other beta-coronaviruses, including SARS-CoV-1 and MERS-CoV, indicating that they might be useful to treat a broader range of coronaviral infections.

VIR-CRISPR: Visual in-one-tube ultrafast RT-PCR and CRISPR method for instant SARS-CoV-2 detection

Until now, corona virus disease 2019 (COVID-19) remained to be an enormous threat for global health. As one viral illness induced by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), versatile, rapid and sensitive method for SARS-CoV-2 detection in early stage is urgently needed. Here, we reported an ultrasensitive and visual in-one-tube detection method which could be accomplished within half an hour from sampling-to-result. By integrating all reactions in one tube, liquid handling steps were omitted and amplicon contamination could be totally avoided. Magnetic beads were employed to achieve the fast extraction of viral nucleic acid and increase the sensitivity. Using portable thermocycler and blue light, the fluorescent results could be directly observed by naked eyes. The proposed method is of higher specificity and sensitivity, nearly at single molecule level. More important, results demonstrated 100% positive detection rate for 40 clinical samples, which was consistent with standard RT-PCR. Thus, our method is considerably simple, rapid, sensitive and accurate, holding great promise for the instant detecting of viruses including SARS-CoV-2 and the next generation of molecular diagnosis.

A bioinformatic approach of targeting SARS-CoV-2 replication by silencing a conserved alternative reserve of the orf8 gene using host miRNAs

The causative agent of the COVID-19 pandemic, the SARS-CoV-2 virus has yielded multiple relevant mutations, many of which have branched into major variants. The Omicron variant has a huge similarity with the original viral strain (first COVID-19 strain from Wuhan). Among different genes, the highly variable orf8 gene is responsible for crucial host interactions and has undergone multiple mutations and indels. The sequence of the orf8 gene of the Omicron variant is, however, identical with the gene sequence of the wild type. orf8 modulates the host immunity making it easier for the virus to conceal itself and remain undetected. Variants seem to be deleting this gene without affecting the viral replication. While analyzing, we came across the conserved orf7a gene in the viral genome which exhibits a partial sequence homology as well as functional similarity with the SARS-CoV-2 orf8. Hence, we have proposed here in our hypothesis that, orf7a might be an alternative reserve of orf8 present in the virus which was compensating for the lost gene. A computational approach was adopted where we screened various miRNAs targeted against the orf8 gene. These miRNAs were then docked onto the orf8 mRNA sequences. The same set of miRNAs was then used to check for their binding affinity with the orf7a reference mRNA. Results showed that miRNAs targeting the orf8 had favorable shape complementarity and successfully docked with the orf7a gene as well. These findings provide a basis for developing new therapeutic approaches where both orf8 and orf7a can be targeted simultaneously.

Exploring COVID-19 pathogenesis on command-line: A bioinformatics pipeline for handling and integrating omics data

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first identified in late 2019 in Wuhan, China, and has proven to be highly pathogenic, making it a global public health threat. The immediate need to understand the mechanisms and impact of the virus made omics techniques stand out, as they can offer a holistic and comprehensive view of thousands of molecules in a single experiment. Mastering bioinformatics tools to process, analyze, integrate, and interpret omics data is a powerful knowledge to enrich results. We present a robust and open access computational pipeline for extracting information from quantitative proteomics and transcriptomics public data. We present the entire pipeline from raw data to differentially expressed genes. We explore processes and pathways related to mapped transcripts and proteins. A pipeline is presented to integrate and compare proteomics and transcriptomics data using also packages available in the Bioconductor and providing the codes used. Cholesterol metabolism, immune system activity, ECM, and proteasomal degradation pathways increased in infected patients. Leukocyte activation profile was overrepresented in both proteomics and transcriptomics data. Finally, we found a panel of proteins and transcripts regulated in the same direction in the lung transcriptome and plasma proteome that distinguish healthy and infected individuals. This panel of markers was confirmed in another cohort of patients, thus validating the robustness and functionality of the tools presented.

Ex vivo and in vivo suppression of SARS-CoV-2 with combinatorial AAV/RNAi expression vectors

Despite rapid development and deployment of vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), clinically relevant modalities to curb the pandemic by directly attacking the virus on a genetic level remain highly desirable and are urgently needed. Here we comprehensively illustrate the capacity of adeno-associated virus (AAV) vectors co-expressing a cocktail of three short hairpin RNAs (shRNAs; RNAi triggers) directed against the SARS-CoV-2 RdRp and N genes as versatile and effective antiviral agents. In cultured monkey cells and human gut organoids, our most potent vector, SAVIOR (SARS virus repressor), suppressed SARS-CoV-2 infection to background levels. Strikingly, in control experiments using single shRNAs, multiple SARS-CoV-2 escape mutants quickly emerged from infected cells within 24-48 h. Importantly, such adverse viral adaptation was fully prevented with the triple-shRNA AAV vector even during long-term cultivation. In addition, AAV-SAVIOR efficiently purged SARS-CoV-2 in a new model of chronically infected human intestinal cells. Finally, intranasal AAV-SAVIOR delivery using an AAV9 capsid moderately diminished viral loads and/or alleviated disease symptoms in hACE2-transgenic or wild-type mice infected with human or mouse SARS-CoV-2 strains, respectively. Our combinatorial and customizable AAV/RNAi vector complements ongoing global efforts to control the coronavirus disease 2019 (COVID-19) pandemic and holds great potential for clinical translation as an original and flexible preventive or therapeutic antiviral measure.

SARS-CoV-2 Point Mutation and Deletion Spectra and Their Association with Different Disease Outcomes

Mutant spectra of RNA viruses are important to understand viral pathogenesis and response to selective pressures. There is a need to characterize the complexity of mutant spectra in coronaviruses sampled from infected patients. In particular, the possible relationship between SARS-CoV-2 mutant spectrum complexity and disease associations has not been established. In the present study, we report an ultradeep sequencing (UDS) analysis of the mutant spectrum of amplicons from the nsp12 (polymerase)- and spike (S)-coding regions of 30 nasopharyngeal isolates (diagnostic samples) of SARS-CoV-2 of the first COVID-19 pandemic wave (Madrid, Spain, April 2020) classified according to the severity of ensuing COVID-19. Low-frequency mutations and deletions, counted relative to the consensus sequence of the corresponding isolate, were overwhelmingly abundant. We show that the average number of different point mutations, mutations per haplotype, and several diversity indices was significantly higher in SARS-CoV-2 isolated from patients who developed mild disease than in those associated with moderate or severe disease (exitus). No such bias was observed with RNA deletions. Location of amino acid substitutions in the three-dimensional structures of nsp12 (polymerase) and S suggest significant structural or functional effects. Thus, patients who develop mild symptoms may be a richer source of genetic variants of SARS-CoV-2 than patients with moderate or severe COVID-19. IMPORTANCE The study shows that mutant spectra of SARS-CoV-2 from diagnostic samples differ in point mutation abundance and complexity and that significantly larger values were observed in virus from patients who developed mild COVID-19 symptoms. Mutant spectrum complexity is not a uniform trait among isolates. The nature and location of low-frequency amino acid substitutions present in mutant spectra anticipate great potential for phenotypic diversification of SARS-CoV-2.

Both simulation and sequencing data reveal coinfections with multiple SARS-CoV-2 variants in the COVID-19 pandemic

SARS-CoV-2 is a single-stranded RNA betacoronavirus with a high mutation rate. The rapidly emerging SARS-CoV-2 variants could increase transmissibility and diminish vaccine protection. However, whether coinfection with multiple SARS-CoV-2 variants exists remains controversial. This study collected 12,986 and 4,113 SARS-CoV-2 genomes from the GISAID database on May 11, 2020 (GISAID20May11), and Apr 1, 2021 (GISAID21Apr1), respectively. With single-nucleotide variant (SNV) and network clique analyses, we constructed single-nucleotide polymorphism (SNP) coexistence networks and discovered maximal SNP cliques of sizes 16 and 34 in the GISAID20May11 and GISAID21Apr1 datasets, respectively. Simulating the transmission routes and SNV accumulations, we discovered a linear relationship between the size of the maximal clique and the number of coinfected variants. We deduced that the COVID-19 cases in GISAID20May11 and GISAID21Apr1 were coinfections with 3.20 and 3.42 variants on average, respectively. Additionally, we performed Nanopore sequencing on 42 COVID-19 patients and discovered recurrent heterozygous SNPs in twenty of the patients, including loci 8,782 and 28,144, which were crucial for SARS-CoV-2 lineage divergence. In conclusion, our findings reported SARS-CoV-2 variants coinfection in COVID-19 patients and demonstrated the increasing number of coinfected variants.

Dynamic Allostery Highlights the Evolutionary Differences between the CoV-1 and CoV-2 Main Proteases

The SARS-CoV-2 coronavirus has become one of the most immediate and widely-studied systems since its identification and subsequent global outbreak from 2019-2021. In an effort to understand the biophysical changes as a result of mutations, the mechanistic details of multiple different proteins within the SARS-CoV-2 virus have been studied and compared with SARS-CoV-1. Focusing on the main protease (mPro), we first explored the long-range dynamics using the Dynamic Coupling Index (DCI) to investigate the dynamic coupling between the catalytic site residues and the rest of the protein, both inter and intra chain, for the CoV-1 and CoV-2 mPro. We found that there is significant cross-chain coupling between these active sites and specific distal residues in the CoV-2 mPro not present in CoV-1. The enhanced long distance interactions, particularly between the two chains, suggest subsequently enhanced cooperativity for CoV-2. A further comparative analysis of the dynamic flexibility using the Dynamic Flexibility Index (DFI) between the CoV-1 and CoV-2 mPros shows that the inhibitor binding near active sites induces change in flexibility to a distal region of the protein, opposite in behavior between the two systems; this region becomes more flexible upon inhibitor binding in CoV-1 while it becomes less flexible in the CoV-2 mPro. Upon inspection, we show that, on average, the dynamic flexibility of the sites substituted from CoV-1 to CoV-2 changes significantly less than the average calculated across all residues within the structure, indicating that the differences in behaviors between the two systems is likely the result of allosteric influence, where the new substitutions in CoV-2 induce flexibility and dynamical changes elsewhere in the structure.

A screening strategy for identifying the dominant variant of SARS-COV-2 in the fifth peak of Kurdistan- Iran population using HRM and Probe-based RT-PCR assay

By the emergence of SARS CoV-2 variants, many studies were developed to deal with it. The high transmissibility and mortality rate of some variants, in particular developing countries have caused the operation of simple diagnostic tests for genomic surveillance. In this study, we developed two assays of High Resolution Melting (HRM) and Probe-based RT-PCR as simple and inexpensive methods to identify the variants. We screened the mutations of del69–70, E484K, E484Q, D614G, L452R, and T478K in 100 cases from SARS-COV-2 positive patients in Kurdistan- Iran population. In general, the result of the two methods overlapped each other, nevertheless, we suggested HRM results be confirmed with a standard assay (Whole-Genome Sequencing). This work indicated that HRM as the rapid and inexpensive method could identify and categorize the variants of SARS CoV-2 and reduce the costs for carrying out sequencing.

Discovery of SARS-CoV-2 main protease covalent inhibitors from a DNA-encoded library selection

Covalent inhibitors targeting the main protease (M^pro, or 3CLpro) of SARS-CoV-2 have shown promise in preclinical investigations. Herein, we report the discovery of two new series of molecules that irreversibly bind to SARS-CoV-2 M^pro. These acrylamide containing molecules were discovered using our covalent DNA-encoded library (DEL) screening platform. Following selection against SARS-CoV-2 M^pro, off-DNA compounds were synthesized and investigated to determine their inhibitory effects, the nature of their binding, and to generate preliminary structure-activity relationships. LC-MS analysis indicates a 1:1 (covalent) binding stoichiometry between our hit molecules and SARS-CoV-2 M^pro. Fluorescent staining assay for covalent binding in the presence of cell lysate suggests reasonable selectivity for SARS-CoV-2 M^pro. And lastly, inhibition of enzymatic activity was also observed against a panel of 3CLpro enzymes from different coronavirus strains, with IC₅₀ values ranging from inactive to single digit micromolar. Our results indicate that DEL selection is a useful approach for identifying covalent inhibitors of cysteine proteases.

Visualized Genotyping from "Sample to Results" Within 25 Minutes by Coupling Recombinase Polymerase Amplification (RPA) With Allele-Specific Invasive Reaction Assisted Gold Nanoparticle Probes Assembling

A simple and rapid genotyping method with less-instrumentation is essential for realizing point-of-care detection of personalized medicine-related gene biomarkers. Herein, we developed a rapid and visualized genotyping method by coupling recombinase polymerase amplification (RPA) with allele-specific invader reaction assisted gold nanoparticle probes assembling. In the method, the DNA targets were firstly amplified by using RPA, which is a rapid isothermal amplification technology. Then an allele-specific invasion reaction was performed to recognize the single nucleotide polymorphisms (SNPs) site in the amplicons, to produce signal molecules that caused discoloration of gold nanoparticle probes. As a result, genotyping was achieved by observing the color change of the reaction by using naked eye without the requirement for any expensive instrument. In order to achieve rapid genotyping detection, the genomic DNA from oral swab lysate samples were used for the RPA templates amplification. In this way, a visualized genotyping from "samples to results" within 25 min was realized. Two clopidogrel related SNPs CYP2C19*2 and CYP2C19*3 of 56 clinical samples were correctly genotyped by using this rapid visualized genotyping assay. In addition, the feasibility for this pathogen genotyping method was also verified by detecting plasmid DNA containing three SARS-COV-2 gene mutation sites, indicating that this method has the potential for clinical sample detection.

Discovery of 9,10-dihydrophenanthrene derivatives as SARS-CoV-2 3CLpro inhibitors for treating COVID-19

The epidemic coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has now spread worldwide and efficacious therapeutics are urgently needed. 3-Chymotrypsin-like cysteine protease (3CLpro) is an indispensable protein in viral replication and represents an attractive drug target for fighting COVID-19. Herein, we report the discovery of 9,10-dihydrophenanthrene derivatives as non-peptidomimetic and non-covalent inhibitors of the SARS-CoV-2 3CLpro. The structure-activity relationships of 9,10-dihydrophenanthrenes as SARS-CoV-2 3CLpro inhibitors have carefully been investigated and discussed in this study. Among all tested 9,10-dihydrophenanthrene derivatives, C1 and C2 display the most potent SARS-CoV-2 3CLpro inhibition activity, with IC50 values of 1.55 ± 0.21 μM and 1.81 ± 0.17 μM, respectively. Further enzyme kinetics assays show that these two compounds dose-dependently inhibit SARS-CoV-2 3CLprovia a mixed-inhibition manner. Molecular docking simulations reveal the binding modes of C1 in the dimer interface and substrate-binding pocket of the target. In addition, C1 shows outstanding metabolic stability in the gastrointestinal tract, human plasma, and human liver microsome, suggesting that this agent has the potential to be developed as an orally administrated SARS-CoV-2 3CLpro inhibitor.

Distinct Symptoms and Underlying Comorbidities with Latitude and Longitude in COVID-19: A Systematic Review and Meta-Analysis

The coronavirus disease 2019 (COVID-19) pandemic is straining global health resources, and the prevalence of severe disease appears to vary across countries. In accordance with PRISMA guidelines, we performed a systematic review and meta-analysis of clinical features and underlying medical conditions of COVID-19. Eighty-seven studies, involving 1,434,931 COVID-19 patients from the Americas, Asia, Europe, and Oceania, were included. Geographically, the rate of severity was highest in Asia (95% confidence interval (CI) 0.23‒0.30). The rates of comorbidities of COVID-19 patients in the Americas were significantly higher than those in Asia. Most Asian patients had fever (95%CI 0.70‒0.81), and most Oceanian patients had cough (95%CI 0.68‒0.70) as their prevalent symptom. Dyspnea was common in the Americas (95%CI 0.33‒0.64), Europe (95%CI 0.29‒0.64), and high latitude regions (95%CI 0.53‒0.82). European patients exhibited significantly high rates of loss of smell and taste (95%CI 0.60-0.97). In low-latitude regions, cancer (95%CI 14.50‒4.89) had the strongest correlation with illness severity. Comorbid diseases and clinical manifestations of severe COVID-19 patients vary substantially between latitudes and longitudes. Region-specific care should be considered to treat and improve the prognosis of COVID-19 patients.

The mutational dynamics of the SARS-CoV-2 virus in serial passages in vitro

Since its outbreak in 2019, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) keeps surprising the medical community by evolving diverse immune escape mutations in a rapid and effective manner. To gain deeper insight into mutation frequency and dynamics, we isolated ten ancestral strains of SARS-CoV-2 and performed consecutive serial incubation in ten replications in a suitable and common cell line and subsequently analysed them using RT-qPCR and whole genome sequencing. Along those lines we hoped to gain fundamental insights into the evolutionary capacity of SARS-CoV-2 in vitro. Our results identified a series of adaptive genetic changes, ranging from unique convergent substitutional mutations and hitherto undescribed insertions. The region coding for spike proved to be a mutational hotspot, evolving a number of mutational changes including the already known substitutions at positions S:484 and S:501. We discussed the evolution of all specific adaptations as well as possible reasons for the seemingly inhomogeneous potential of SARS-CoV-2 in the adaptation to cell culture. The combination of serial passage in vitro with whole genome sequencing uncovers the immense mutational potential of some SARS-CoV-2 strains. The observed genetic changes of SARS-CoV-2 in vitro could not be explained solely by selectively neutral mutations but possibly resulted from the action of directional selection accumulating favourable genetic changes in the evolving variants, along the path of increasing potency of the strain. Competition among a high number of quasi-species in the SARS-CoV-2 in vitro population gene pool may reinforce directional selection and boost the speed of evolutionary change.

•

Ten genetically similar strains evolved very differently in serial passage in vitro.

•

Observed mutations included substitutions at important spike positions.

•

The three strains with the highest replication rates developed two convergent mutations.

•

Via directional selection favourable genetic changes are accumulated.

•

Competition among many quasi-species boosts the speed of evolutionary change.

Genomic, proteomic and metabolomic profiling of severe acute respiratory syndrome-Coronavirus-2

Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome-Coronavirus-2 (SARS-CoV-2) is one of the worst human health problems faced by humanity in recent centuries. An end to this health crisis relies on our ability to monitor viral transmission dynamics to check spread, develop therapeutics and preventatives for treatment of SARS-CoV-2 infection and understand the pathophysiology of the disease for better management of the patients. Omics technologies have played a crucial part in understanding the different aspects of COVID-19 disease. While whole-genome sequencing of SARS-CoV-2 isolates from across the globe has aided in the development of molecular diagnostic assays and informed about the viral evolution, knowledge of structure and function of viral proteome fueled the development of small molecule and biologicals therapeutics as well as vaccines. Concurrently, metabolomic profiling of samples from COVID-19 patients experiencing a varying level of disease severity has provided a snapshot of the pathophysiology of the disease helping device effective treatment regimen. This chapter deals with genomic, proteomic, and metabolomic profiling of SRAS-CoV-2.

A global analysis of conservative and non-conservative mutations in SARS-CoV-2 detected in the first year of the COVID-19 world-wide diffusion

The ability of SARS-CoV-2 to rapidly mutate represents a remarkable complicancy. Quantitative evaluations of the effects that these mutations have on the virus structure/function is of great relevance and the availability of a large number of SARS-CoV-2 sequences since the early phases of the pandemic represents a unique opportunity to follow the adaptation of the virus to humans. Here, we evaluated the SARS-CoV-2 amino acid mutations and their progression by analyzing publicly available viral genomes at three stages of the pandemic (2020 March 15th and October 7th, 2021 February 7th). Mutations were classified in conservative and non-conservative based on the probability to be accepted during the evolution according to the Point Accepted Mutation substitution matrices and on the similarity of the encoding codons. We found that the most frequent substitutions are T > I, L > F, and A > V and we observe accumulation of hydrophobic residues. These findings are consistent among the three stages analyzed. We also found that non-conservative mutations are less frequent than conservative ones. This finding may be ascribed to a progressive adaptation of the virus to the host. In conclusion, the present study provides indications of the early evolution of the virus and tools for the global and genome-specific evaluation of the possible impact of mutations on the structure/function of SARS-CoV-2 variants.

Evolutionary Dynamics of Indels in SARS-CoV-2 Spike Glycoprotein

SARS-CoV-2, responsible for the current COVID-19 pandemic that claimed over 5.0 million lives, belongs to a class of enveloped viruses that undergo quick evolutionary adjustments under selection pressure. Numerous variants have emerged in SARS-CoV-2, posing a serious challenge to the global vaccination effort and COVID-19 management. The evolutionary dynamics of this virus are only beginning to be explored. In this work, we have analysed 1.79 million spike glycoprotein sequences of SARS-CoV-2 and found that the virus is fine-tuning the spike with numerous amino acid insertions and deletions (indels). Indels seem to have a selective advantage as the proportions of sequences with indels steadily increased over time, currently at over 89%, with similar trends across countries/variants. There were as many as 420 unique indel positions and 447 unique combinations of indels. Despite their high frequency, indels resulted in only minimal alteration of N-glycosylation sites, including both gain and loss. As indels and point mutations are positively correlated and sequences with indels have significantly more point mutations, they have implications in the evolutionary dynamics of the SARS-CoV-2 spike glycoprotein.

Rapid, inexpensive methods for exploring SARS CoV-2 D614G mutation

A common mutation has occurred in the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS CoV-2), known as D614G (A23403G). There are discrepancies in the impact of this mutation on the virus's infectivity. Also, the whole genome sequencings are expensive and time-consuming. This study aims to develop three fast economical assays for prompt identifications of the D614G mutation including Taqman probe-based real-time reverse transcriptase polymerase chain reaction (rRT PCR), an amplification refractory mutation system (ARMS) RT and restriction fragment length polymorphism (RFLP), in nasopharyngeal swab samples. Both rRT and ARMS data showed G614 mutants indicated by the presence of HEX probe and 176 bp, respectively. Additionally, the results of the RFLP data and DNA sequencings confirmed the prevalence of the G614 mutants. These methods will be important, in epidemiological, reinfections and zoonotic aspects, through detecting the G614 mutant in retro-perspective samples to track its origins and future re-emergence of D614 wild type.

Impact of the Nucleic Acid Extraction Method and the RT-qPCR Assay on SARS-CoV-2 Detection in Low-Viral Samples

COVID-19 was initially reported in China at the end of 2019 and soon thereafter, in March 2020, the WHO declared it a pandemic. Until October 2021, over 240 million COVID-19 cases were recorded, with 4.9 mln deaths. In order to stop the spread of this disease, it is crucial to monitor and detect any infected person. The etiologic agent of COVID-19 is a novel coronavirus called SARS-CoV-2. The gold standard for the detection of the virus is the RT-qPCR method. This study evaluated two RNA extraction methods and four commercial RT-qPCR assays routinely used in diagnostic laboratories for detecting SARS-CoV-2 in human specimens from the upper respiratory tract. We analyzed a panel of 70 clinical samples with varying RNA loads. Our study demonstrated the significant impact of the diagnostic methods selected by the laboratory on the SARS-CoV-2 detection in clinical specimens with low viral loads.

A Global Mutational Profile of SARS-CoV-2: A Systematic Review and Meta-Analysis of 368,316 COVID-19 Patients

Since its first detection in December 2019, more than 232 million cases of COVID-19, including 4.7 million deaths, have been reported by the WHO. The SARS-CoV-2 viral genomes have evolved rapidly worldwide, causing the emergence of new variants. This systematic review and meta-analysis was conducted to provide a global mutational profile of SARS-CoV-2 from December 2019 to October 2020. The review was conducted according to the Preferred Reporting Items for Systematic Reviews and Meta-analysis (PRISMA), and a study protocol was lodged with PROSPERO. Data from 62 eligible studies involving 368,316 SARS-CoV-2 genomes were analyzed. The mutational data analyzed showed most studies detected mutations in the Spike protein (n = 50), Nucleocapsid phosphoprotein (n = 34), ORF1ab gene (n = 29), 5′-UTR (n = 28) and ORF3a (n = 25). Under the random-effects model, pooled prevalence of SARS-CoV-2 variants was estimated at 95.1% (95% CI; 93.3–96.4%; I² = 98.952%; p = 0.000) while subgroup meta-analysis by country showed majority of the studies were conducted ‘Worldwide’ (n = 10), followed by ‘Multiple countries’ (n = 6) and the USA (n = 5). The estimated prevalence indicated a need to continuously monitor the prevalence of new mutations due to their potential influence on disease severity, transmissibility and vaccine effectiveness.

RNA Viromics of Southern California Wastewater and Detection of SARS-CoV-2 Single-Nucleotide Variants

Municipal wastewater provides an integrated sample of a diversity of human-associated microbes across a sewershed, including viruses. Wastewater-based epidemiology (WBE) is a promising strategy to detect pathogens and may serve as an early warning system for disease outbreaks. Notably, WBE has garnered substantial interest during the coronavirus disease 2019 (COVID-19) pandemic to track disease burden through analyses of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA. Throughout the COVID-19 outbreak, tracking SARS-CoV-2 in wastewater has been an important tool for understanding the spread of the virus. Unlike traditional sequencing of SARS-CoV-2 isolated from clinical samples, which adds testing burden to the health care system, in this study, metatranscriptomics was used to sequence virus directly from wastewater. Here, we present a study in which we explored RNA viral diversity through sequencing 94 wastewater influent samples across seven wastewater treatment plants (WTPs), collected from August 2020 to January 2021, representing approximately 16 million people in Southern California. Enriched viral libraries identified a wide diversity of RNA viruses that differed between WTPs and over time, with detected viruses including coronaviruses, influenza A, and noroviruses. Furthermore, single-nucleotide variants (SNVs) of SARS-CoV-2 were identified in wastewater, and we measured proportions of overall virus and SNVs across several months. We detected several SNVs that are markers for clinically important SARS-CoV-2 variants along with SNVs of unknown function, prevalence, or epidemiological consequence. Our study shows the potential of WBE to detect viruses in wastewater and to track the diversity and spread of viral variants in urban and suburban locations, which may aid public health efforts to monitor disease outbreaks. IMPORTANCE Wastewater-based epidemiology (WBE) can detect pathogens across sewersheds, which represents the collective waste of human populations. As there is a wide diversity of RNA viruses in wastewater, monitoring the presence of these viruses is useful for public health, industry, and ecological studies. Specific to public health, WBE has proven valuable during the coronavirus disease 2019 (COVID-19) pandemic to track the spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) without adding burden to health care systems. In this study, we used metatranscriptomics and reverse transcription-droplet digital PCR (RT-ddPCR) to assay RNA viruses across Southern California wastewater from August 2020 to January 2021, representing approximately 16 million people from Los Angeles, Orange, and San Diego counties. We found that SARS-CoV-2 quantification in wastewater correlates well with county-wide COVID-19 case data, and that we can detect SARS-CoV-2 single-nucleotide variants through sequencing. Likewise, wastewater treatment plants (WTPs) harbored different viromes, and we detected other human pathogens, such as noroviruses and adenoviruses, furthering our understanding of wastewater viral ecology.

Human bronchial-pulmonary proteomics in coronavirus disease 2019 (COVID-19) pandemic: applications and implications

INTRODUCTION

The outbreak of the newly discovered human coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has disrupted the normal life of almost every civilization worldwide. Studies have shown that the coronavirus disease 2019 (COVID-19) caused by the SARS-CoV-2 can affect multiple human organs and physiological systems, but the respiratory system remains the primary location for viral infection.

AREAS COVERED

We summarize how omics technologies are used in SARS-CoV-2 research and specifically review the current knowledge of COVID-19 from the aspect of human bronchial-pulmonary proteomics. Also, knowledge gaps in COVID-19 that can be fulfilled by proteomics are discussed.

EXPERT OPINION

Overall, human bronchial-pulmonary proteomics plays an important role in revealing the dynamics, functions, tropism, and pathogenicity of SARS-CoV-2, which is crucial for COVID-19 biomarker and therapeutic target discoveries. To more fully understand the impact of COVID-19, research from various angles using multi-omics approaches should also be conducted on the lungs as well as other organs.