Structural topological analysis of spike proteins of SARS-CoV-2 variants of concern highlight distinctive amino acid substitution patterns

Since the onset of pandemic in 2019, SARS-CoV-2 has diverged into numerous variants driven by antigenic and infectivity-oriented selection. Some variants have accumulated fitness-enhancing mutations, evaded immunity and spread despite global vaccination campaigns. The spike (S) glycoprotein of SARS-CoV-2 demonstrated the greatest immunogenicity and amino acid substitution diversity owing to its importance in the interaction with human angiotensin receptor 2 (hACE2). The S protein consistently emerges as an amino acid substitution (AAS) hotspot in all six lineages, however, in Omicron this enrichment is significantly higher. This study attempts to design and validate a method of mapping S-protein substitution profile across variants to identify the conserved and AAS regions. A substitution matrix was created based on publicly available databases, and the substitution localization was illustrated on a cryo-electron microscopy generated S-protein model. Our analyses indicated that the diversity of N-terminal (NTD) and receptor-binding (RBD) domains exceeded that of any other regions but still contained extended low substitution density regions particularly considering significantly broader substitution profiles of Omicron BA.2 and BA.4/5. Finally, the substitution matrix was compared to a random sample alignment of variant sequences, revealing discrepancies. Therefore, it was suggested to improve matrix accuracy by processing a large number of S-protein sequences using an automated algorithm. Several critical immunogenic and receptor-interacting residues were identified in the conserved regions within NTD and RBD. In conclusion, the structural and topological analysis of S proteins of SARS-CoV-2 variants highlight distinctive amino acid substitution patterns which may be foundational in predicting future variants.     

Severity of the Omicron SARS-CoV-2 variant compared with the previous lineages: A systematic review

The Omicron variant was first detected in October 2021, which evolved from the original SARS-CoV-2 strain and was found to possess many mutations. Immune evasion was one of the notable consequences of these mutations. Despite Omicron exhibiting increased transmissibility, the rates of hospitalizations and deaths among patients infected with this variant were substantially lower when compared to other strains. However, concluding that the Omicron variant is less severe than other variants of SARS-CoV-2 requires consideration of multiple factors, including the vaccination status of infected patients as well as any previous infections with other variants. This review compiled data about any reported indicators of severity in Omicron-infected patients, including studies comparing Omicron with other variants while adjusting for confounders. A comprehensive search was conducted using different databases to target any studies about Omicron. In total, 62 studies met our inclusion criteria and were included in this study. Many studies reported a significantly reduced risk of hospitalization, ICU admission, need for oxygenation/ventilation, and death in Omicron-infected patients compared to patients infected with other variants, such as Delta. Some studies, however, reported comparable severity in Omicron infected patients as to other variants emphasizing a substantial risk for severe illness. Furthermore, the COVID-19 vaccines were less effective against Omicron relative to previous lineages, except after receiving the booster dose. One study recommended vaccination during pregnancy, which may help prevent future cases of severe SARS-CoV-2 pneumonia in neonates and young infants due to the transfer of humoral response from the mother.     

Clinical characteristics of 310 SARS-CoV-2 Omicron variant patients and comparison with Delta and Beta variant patients in China

Although the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant has spread, data on the clinical characteristics of infected patients are limited. In this study, the demographic, clinical characteristics, and laboratory data of 310 SARS-CoV-2 Omicron variant patients treated at Haihe Hospital of Tianjin were collected and analyzed. Information on these patients was compared to 96 patients with the Delta variant of concern (VOC) and 326 patients with the Beta VOC during the previous coronavirus disease 2019 (COVID-19) outbreak in Harbin. Of the 310 patients infected with the Omicron variant, the median age was 35 years. Most patients were clinically classified as mild (57.74%), and the most common symptoms were cough (48.71%), fever (39.35%), and sore throat (38.26%). The results for different vaccination groups in the Omicron group showed that the median of “SARS-CoV-2 specific IgG” after 2 or 3 doses of vaccination was higher than the unvaccinated group (all Ps < 0.05). Older age was associated with a higher proportion of moderate cases and lower asymptomatic and mild cases based on clinical classifications. Compared to the Delta and Beta groups, the median age of the Omicron group was younger. The total number of asymptomatic patients and mild patients in the Omicron virus group was higher than the Delta and Beta groups (60.97% vs. 54.17% vs. 47.55%). This study presented the clinical characteristics of the first group of patients infected with the Omicron variant in Tianjin, China, and compared their clinical features with patients infected by the Delta and Beta variants, which would increase our understanding of the characteristics of SARS-CoV-2 Omicron variant.     

Evidence for a mouse origin of the SARS-CoV-2 Omicron variant

The rapid accumulation of mutations in the SARS-CoV-2 Omicron variant that enabled its outbreak raises questions as to whether its proximal origin occurred in humans or another mammalian host. Here, we identified 45 point mutations that Omicron acquired since divergence from the B.1.1 lineage. We found that the Omicron spike protein sequence was subjected to stronger positive selection than that of any reported SARS-CoV-2 variants known to evolve persistently in human hosts, suggesting a possibility of host-jumping. The molecular spectrum of mutations (i.e., the relative frequency of the 12 types of base substitutions) acquired by the progenitor of Omicron was significantly different from the spectrum for viruses that evolved in human patients but resembled the spectra associated with virus evolution in a mouse cellular environment. Furthermore, mutations in the Omicron spike protein significantly overlapped with SARS-CoV-2 mutations known to promote adaptation to mouse hosts, particularly through enhanced spike protein binding affinity for the mouse cell entry receptor. Collectively, our results suggest that the progenitor of Omicron jumped from humans to mice, rapidly accumulated mutations conducive to infecting that host, then jumped back into humans, indicating an inter-species evolutionary trajectory for the Omicron outbreak.     

In-silico genomic landscape characterization and evolution of SARS-CoV-2 variants isolated in India shows significant drift with high frequency of mutations

In-silico studies on SARS-CoV-2 genome are considered important to identify the significant pattern of variations and its possible effects on the structural and functional characteristics of the virus. The current study determined such genetic variations and their possible impact among SARS-CoV-2 variants isolated in India. A total of 546 SARS-CoV-2 genomic sequences (India) were retrieved from the gene bank (NCBI) and subjected to alignment against the Wuhan variant (NC_045512.2), the corresponding amino acid changes were analyzed using NCBI Protein-BLAST. These 546 variants revealed 841 mutations; most of these were non-synonymous 464/841 (55.1%), there was no identical variant compared to the original strain. All genes; coding and non-coding showed nucleotide changes, most of the structural genes showed frequent nonsynonymous mutations. The most affected genes were ORF1a/b followed by the S gene which showed 515/841 (61.2%) and 120/841 (14.3%) mutations, respectively. The most frequent non-synonymous mutation 486/546 (89.01%) occurred in the S gene (structural gene) at position 23,403 where A changed to G leading to the replacement of aspartic acid by glycine in position (D614G). Interestingly, four variants also showed deletion. The variants MT800923 and MT800925 showed 12 consecutive nucleotide deletion in position 21982–21993 resulting in 4 consecutive amino acid deletions that were leucine, glycine, valine, and tyrosine in positions 141, 142, 143, and 144 respectively. The present study exhibited a higher mutations rate per variant compared to other studies carried out in India.     

The Global Landscape of SARS-CoV-2 Genomes,Variants, and Haplotypes in 2019nCoVR

On January 22, 2020, China National Center for Bioinformation (CNCB) released the 2019 Novel Coronavirus Resource (2019nCoVR), an open-access information resource for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). 2019nCoVR features a comprehensive integration of sequence and clinical information for all publicly available SARS-CoV-2 isolates, which are manually curated with value-added annotations and quality evaluated by an automated in-house pipeline. Of particular note, 2019nCoVR offers systematic analyses to generate a dynamic landscape of SARS-CoV-2 genomic variations at a global scale. It provides all identified variants and their detailed statistics for each virus isolate, and congregates the quality score, functional annotation, and population frequency for each variant. Spatiotemporal change for each variant can be visualized and historical viral haplotype network maps for the course of the outbreak are also generated based on all complete and high-quality genomes available. Moreover, 2019nCoVR provides a full collection of SARS-CoV-2 relevant literature on the coronavirus disease 2019 (COVID-19), including published papers from PubMed as well as preprints from services such as bioRxiv and medRxiv through Europe PMC. Furthermore, by linking with relevant databases in CNCB, 2019nCoVR offers data submission services for raw sequence reads and assembled genomes, and data sharing with NCBI. Collectively, SARS-CoV-2 is updated daily to collect the latest information on genome sequences, variants, haplotypes, and literature for a timely reflection, making 2019nCoVR a valuable resource for the global research community. 2019nCoVR is accessible at https://bigd.big.ac.cn/ncov/.     

Mutations in spike protein T cell epitopes of SARS-COV-2 variants: Plausible influence on vaccine efficacy

With emerging SARS-CoV-2 variants, vaccines approved so far are under scrutiny for long term effectiveness against the circulating strains. There is a prevalent obsession with humoral immunity as in vitro studies have indicated diminished effects of vaccine-induced neutralizing antibodies. However, this need not clinically translate to vaccine resistance as immune response against all forms of present vaccine preparations is T dependent unlike that against native viral particles which can induce T independent immune responses. Thus, we focused on this major correlate of protection against infections, T cell response. Using bioinformatics tools, we analyzed SARS-CoV-2 Spike protein T cell epitopes and their diversity across Delta plus/B.1.617.2.1, Gamma/P.1 (variant of concern), B.1.1.429, Zeta/P.2 and Mink cluster 5/B.1.1.298 variants as well as Omicron/B.1.1.529 (variant of concern). We also compared HLA restriction profiles of the mutant epitopes with that of the native epitopes (from Wuhan_hu_1 strain, used in vaccine formulations). Our observations show ~90% conservation of CD4+ and CD8+ epitopes across Delta plus/B.1.617.2.1, Gamma/P.1 (variant of concern), B.1.1.429, Zeta/P.2 and Mink cluster 5/B.1.1.298. For the Omicron/B.1.1.529 variant, ~75% of CD4+ and ~ 87% CD8+ epitopes were conserved. Majority of the mutated CD4+ and CD8+ epitopes of this variant were predicted to retain the HLA restriction pattern as their native epitopes. The results of our bioinformatics analysis suggest largely conserved T cell responses across the studied variants, ability of T cells to tackle new SARS-CoV-2 variants and aid in protection from COVID-19 post vaccination. In conclusion, the results suggest that current vaccines may not be rendered completely ineffective against new variants.     

Mutational sensitivity of D614G in spike protein of SARS-CoV-2 in Jordan

BackgroundSpike protein is the surface glycoprotein of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) necessary for the entry of the virus via the transmembrane receptors of the human respiratory cells causing COVID-19 disease.AimHere, we aimed to predict the three-dimensional monomer structure of spike protein of SARS-CoV-2 from 20 Jordanian nasopharyngeal samples and to determine the percentage of single amino acid variants (SAV) in the spike protein of SARS-CoV-2.MethodsThe output of the Protein Homology/analogY Recognition Engine V 2.0 (Phyre2) found four single amino acid variants in the spike gene.ResultsThe first variant represented by 5% of samples that showed tyrosine deletion at Y144 located in the N terminal domain. The second and the dominant variant, represented by 62%, showed aspartate a coil amino acid substitution to glycine an extracellular amino acid at D614G located in the spike recognition binding site. The third variant, represented by 5%, showed aspartate substitution to tyrosine at D1139Y, and the fourth variant, represented by 5% glycine substitution to serine at G1167S.ConclusionOur results have shown low mutational sensitivity in all variants except to D614G the one with the most likely neutral mutational sensitivity that all variants might not explicitly affect the function of spike glycoprotein. However, D614G might change the viral conformational plasticity and hence a potential viral fitness gain but one must be cautious about drawing any concrete conclusions about the severity of symptoms and viral transmission from genomic data only.General significanceStudying mutations such as D614G in deep is essential to control the pandemic in terms of immune systems, antibodies, or even vaccines.     

The humoral and cellular immune evasion of SARS-CoV-2 Omicron and sub-lineages

The recently discovered SARS-CoV-2 variant Omicron (B.1.1.529) has rapidly become a global public health issue. The substantial mutations in the spike protein in this new variant have raised concerns about its ability to escape from pre-existing immunity established by natural infection or vaccination. In this review, we give a summary of current knowledge concerning the antibody evasion properties of Omicron and its subvariants (BA.2, BA.2.12.1, BA.4/5, and BA.2.75) from therapeutic monoclonal antibodies and the sera of SARS-CoV-2 vaccine recipients or convalescent patients. We also summarize whether vaccine-induced cellular immunity (memory B cell and T cell response) can recognize Omicron specifically. In brief, the Omicron variants demonstrated remarkable antibody evasion, with even more striking antibody escape seen in the Omicron BA.4 and BA.5 sub-lineages. Luckily, the third booster vaccine dose significantly increased the neutralizing antibodies titers, and the vaccine-induced cellular response remains conserved and provides second-line defense against the Omicron.     

High-resolution melting analysis to discriminate between the SARS-CoV-2 Omicron variants BA.1 and BA.2

High-resolution melting (HRM) analysis was conducted to discriminate between SARS-CoV-2 Omicron variant BA.1 (B.1.1.529.1) and subvariant BA.2 (B.1.1.529.2). We performed two-step PCR consisting of the first PCR and the second nested PCR to prepare the amplicon for HRM analysis, which detected G339D, N440K, G446S and D796Y variations in the SARS-CoV-2 spike protein. The melting temperatures (Tms) of the amplicons from the cDNA of the Omicron variant BA.1 and subvariant BA.2 receptor binding domain (RBD) in spike protein were the same: 75.2 °C (G339D variation) and 73.4 °C (D796Y variation). These Tms were distinct from those of SARS-CoV-2 isolate Wuhan-Hu-1, and were specific to the Omicron variant. In HRM analyses that detected the N440K and G446S variations, the Tms of amplicons from the cDNA of the Omicron variant BA.1 and subvariant BA.2 RBDs were 73.0 °C (N440K and G446S variations) and 73.5 °C (G446S variation). This difference indicates that the SARS-CoV-2 Omicron variants BA.1 and BA.2 can be clearly discriminated. Our study demonstrates the usefulness of HRM analysis after two-step PCR for the discrimination of SARS-CoV-2 variants.     

SARS-CoV-2 VOCs,Mutational diversity and clinical outcome: Are they modulating drug efficacy by altered binding strength?

The global COVID-19 pandemic continues due to emerging Severe Acute Respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOC). Here, we performed comprehensive analysis of in-house sequenced SARS-CoV-2 genome mutations dynamics in the patients infected with the VOCs - Delta and Omicron, within Recovered and Mortality patients. Statistical analysis highlighted significant mutations - T4685A, N4992N, and G5063S in RdRp; T19R in NTD spike; K444N and N532H in RBD spike, associated with Delta mortality. Mutations, T19I in NTD spike, Q493R and N440K in the RBD spike were significantly associated with Omicron mortality. We performed molecular docking for possible effect of significant mutations on the binding of Remdesivir. We found that Remdesivir showed less binding efficacy with the mutant Spike protein of both Delta and Omicron mortality compared to recovered patients. This indicates that mortality associated mutations could have a modulatory effect on drug binding which could be associated with disease outcome.     

Genomic Epidemiology of SARS-CoV-2 in Pakistan

COVID-19 has swept globally and Pakistan is no exception. To investigate the initial introductions and transmissions of the SARS-CoV-2 in Pakistan, we performed the largest genomic epidemiology study of COVID-19 in Pakistan and generated 150 complete SARS-CoV-2 genome sequences from samples collected from March 16 to June 1, 2020. We identified a total of 347 mutated positions, 31 of which were over-represented in Pakistan. Meanwhile, we found over 1000 intra-host single-nucleotide variants (iSNVs). Several of them occurred concurrently, indicating possible interactions among them or coevolution. Some of the high-frequency iSNVs in Pakistan were not observed in the global population, suggesting strong purifying selections. The genomic epidemiology revealed five distinctive spreading clusters. The largest cluster consisted of 74 viruses which were derived from different geographic locations of Pakistan and formed a deep hierarchical structure, indicating an extensive and persistent nation-wide transmission of the virus that was probably attributed to a signature mutation (G8371T in ORF1ab) of this cluster. Furthermore, 28 putative international introductions were identified, several of which are consistent with the epidemiological investigations. In all, this study has inferred the possible pathways of introductions and transmissions of SARS-CoV-2 in Pakistan, which could aid ongoing and future viral surveillance and COVID-19 control.     

Characteristics of replication and pathogenicity of SARS-CoV-2 Alpha and Delta isolates

The continuously arising of SARS-CoV-2 variants has been posting a great threat to public health safety globally, from B.1.17 (Alpha), B.1.351 (Beta), P.1 (Gamma), B.1.617.2 (Delta) to B.1.1.529 (Omicron). The emerging or reemerging of the SARS-CoV-2 variants of concern is calling for the constant monitoring of their epidemics, pathogenicity and immune escape. In this study, we aimed to characterize replication and pathogenicity of the Alpha and Delta variant strains isolated from patients infected in Laos. The amino acid mutations within the spike fragment of the isolates were determined via sequencing. The more efficient replication of the Alpha and Delta isolates was documented than the prototyped SARS-CoV-2 in Calu-3 and Caco-2 cells, while such features were not observed in Huh-7, Vero E6 and HPA-3 cells. We utilized both animal models of human ACE2 (hACE2) transgenic mice and hamsters to evaluate the pathogenesis of the isolates. The Alpha and Delta can replicate well in multiple organs and cause moderate to severe lung pathology in these animals. In conclusion, the spike protein of the isolated Alpha and Delta variant strains was characterized, and the replication and pathogenicity of the strains in the cells and animal models were also evaluated.     

严重急性呼吸综合征冠状病毒2 Omicron变异株研究进展

世界卫生组织(World Health Organization, WHO)于2021年11月26日将首次在南非报告的新型冠状病毒 B．1.1.529 变异株列为受关注变种(variant of concern, VOC),并将其命名为奥密克戎(Omicron)。该变异株存在约50个突变,仅在刺突蛋白区域就有至少30个突变,远远超过其他流行株的突变位点数量。根据对突变位点的分析以及初步实验证实,该毒株可能具有极强的传染性以及免疫逃逸能力。Omicron变异株会怎样影响新冠疫情的走向引起了各国的广泛关注,本文将从Omicron变异株的基本特征、检测、致病性、传染性、免疫逃逸等方面进行综述。     

Host microRNAs exhibit differential propensity to interact with SARS-CoV-2 and variants of concern

A significant number of SARS-CoV-2-infected individuals naturally overcome viral infection, suggesting the existence of a potent endogenous antiviral mechanism. As an innate defense mechanism, microRNA (miRNA) pathways in mammals have evolved to restrict viruses, besides regulating endogenous mRNAs. In this study, we systematically examined the complete repertoire of human miRNAs for potential binding sites on SARS-CoV-2 Wuhan-Hu-1, Beta, Delta, and Omicron. Human miRNA and viral genome interaction were analyzed using RNAhybrid 2.2 with stringent parameters to identify highly bonafide miRNA targets. Using publicly available data, we filtered for miRNAs expressed in lung epithelial cells/tissue and oral keratinocytes, concentrating on the miRNAs that target SARS-CoV-2 S protein mRNAs. Our results show a significant loss of human miRNA and SARS-CoV-2 interactions in Omicron (130 miRNAs) compared to Wuhan-Hu-1 (271 miRNAs), Beta (279 miRNAs), and Delta (275 miRNAs). In particular, hsa-miR-3150b-3p and hsa-miR-4784 show binding affinity for S protein of Wuhan strain but not Beta, Delta, and Omicron. Loss of miRNA binding sites on N protein was also observed for Omicron. Through Ingenuity Pathway Analysis (IPA), we examined the experimentally validated and highly predicted functional role of these miRNAs. We found that hsa-miR-3150b-3p and hsa-miR-4784 have several experimentally validated or highly predicted target genes in the Toll-like receptor, IL-17, Th1, Th2, interferon, and coronavirus pathogenesis pathways. Focusing on the coronavirus pathogenesis pathway, we found that hsa-miR-3150b-3p and hsa-miR-4784 are highly predicted to target MAPK13. Exploring miRNAs to manipulate viral genome/gene expression can provide a promising strategy with successful outcomes by targeting specific VOCs.     

COVID-19: Omicron – the latest,the least virulent,but probably not the last variant of concern of SARS-CoV-2

The Omicron variant rapidly became the dominant SARS-CoV-2 strain in South Africa and elsewhere. This review explores whether this rise was due to an increased transmission of the variant or its escape from population immunity by an extensively mutated spike protein. The mutations affected the structure of the spike protein leading to the loss of neutralization by most, but not all, therapeutic monoclonal antibodies. Omicron also shows substantial immune escape from serum antibodies in convalescent patients and vaccinees. A booster immunization increased, however, the titre and breadth of antiviral antibody response. The cellular immune response against Omicron was largely preserved explaining a satisfying protection of boosted vaccinees against severe infections. Clinicians observed less severe infection with Omicron, but other scientists warned that this must not necessarily reflect less intrinsic virulence. However, in animal experiments with mice and hamsters, Omicron infections also displayed a lesser virulence than previous VOCs and lung functions were less compromised. Cell biologists demonstrated that Omicron differs from Delta by preferring the endocytic pathway for cell entry over fusion with the plasma membrane which might explain Omicron’s distinct replication along the respiratory tract compared with Delta. Omicron represents a distinct evolutionary lineage that deviated from the mainstream of evolving SARS-CoV-2 already in mid-2020 raising questions about where it circulated before getting widespread in December 2021. The role of Omicron for the future trajectory of the COVID-19 pandemic is discussed.