skip to main content


Title: Host-Genome Similarity Characterizes the Adaption of SARS-CoV-2 to Humans
The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a high mutation rate and many variants have emerged in the last 2 years, including Alpha, Beta, Delta, Gamma and Omicron. Studies showed that the host-genome similarity (HGS) of SARS-CoV-2 is higher than SARS-CoV and the HGS of open reading frame (ORF) in coronavirus genome is closely related to suppression of innate immunity. Many works have shown that ORF 6 and ORF 8 of SARS-CoV-2 play an important role in suppressing IFN-β signaling pathway in vivo. However, the relation between HGS and the adaption of SARS-CoV-2 variants is still not clear. This work investigates HGS of SARS-CoV-2 variants based on a dataset containing more than 40,000 viral genomes. The relation between HGS of viral ORFs and the suppression of antivirus response is studied. The results show that ORF 7b, ORF 6 and ORF 8 are the top 3 genes with the highest HGS. In the past 2 years, the HGS values of ORF 8 and ORF 7B of SARS-CoV-2 have increased greatly. A remarkable correlation is discovered between HGS and inhibition of antivirus response of immune system, which suggests that the similarity between coronavirus and host gnome may be an indicator of the suppression of innate immunity. Among the five variants (Alpha, Beta, Delta, Gamma and Omicron), Delta has the highest HGS and Omicron has the lowest HGS. This finding implies that the high HGS in Delta variant may indicate further suppression of host innate immunity. However, the relatively low HGS of Omicron is still a puzzle. By comparing the mutations in genomes of Alpha, Delta and Omicron variants, a commonly shared mutation ACT > ATT is identified in high-HGS strain populations. The high HGS mutations among the three variants are quite different. This finding strongly suggests that mutations in high HGS strains are different in different variants. Only a few common mutations survive, which may play important role in improving the adaptability of SARS-CoV-2. However, the mechanism for how the mutations help SARS-CoV-2 escape immunity is still unclear. HGS analysis is a new method to study virus–host interaction and may provide a way to understand the rapid mutation and adaption of SARS-CoV-2.  more » « less
Award ID(s):
2136095
NSF-PAR ID:
10393684
Author(s) / Creator(s):
Date Published:
Journal Name:
Biomolecules
Volume:
12
Issue:
7
ISSN:
2218-273X
Page Range / eLocation ID:
972
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The glycosylation on the spike (S) protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, modulates the viral infection by altering conformational dynamics, receptor interaction and host immune responses. Several variants of concern (VOCs) of SARS-CoV-2 have evolved during the pandemic, and crucial mutations on the S protein of the virus have led to increased transmissibility and immune escape. In this study, we compare the site-specific glycosylation and overall glycomic profiles of the wild type Wuhan-Hu-1 strain (WT) S protein and five VOCs of SARS-CoV-2: Alpha, Beta, Gamma, Delta and Omicron. Interestingly, both N- and O-glycosylation sites on the S protein are highly conserved among the spike mutant variants, particularly at the sites on the receptor-binding domain (RBD). The conservation of glycosylation sites is noteworthy, as over 2 million SARS-CoV-2 S protein sequences have been reported with various amino acid mutations. Our detailed profiling of the glycosylation at each of the individual sites of the S protein across the variants revealed intriguing possible association of glycosylation pattern on the variants and their previously reported infectivity. While the sites are conserved, we observed changes in the N- and O-glycosylation profile across the variants. The newly emerged variants, which showed higher resistance to neutralizing antibodies and vaccines, displayed a decrease in the overall abundance of complex-type glycans with both fucosylation and sialylation and an increase in the oligomannose-type glycans across the sites. Among the variants, the glycosylation sites with significant changes in glycan profile were observed at both theN-terminal domain and RBD of S protein, with Omicron showing the highest deviation. The increase in oligomannose-type happens sequentially from Alpha through Delta. Interestingly, Omicron does not contain more oligomannose-type glycans compared to Delta but does contain more compared to the WT and other VOCs. O-glycosylation at the RBD showed lower occupancy in the VOCs in comparison to the WT. Our study on the sites and pattern of glycosylation on the SARS-CoV-2 S proteins across the VOCs may help to understand how the virus evolved to trick the host immune system. Our study also highlights how the SARS-CoV-2 virus has conserved bothN- andO- glycosylation sites on the S protein of the most successful variants even after undergoing extensive mutations, suggesting a correlation between infectivity/ transmissibility and glycosylation.

     
    more » « less
  2. Due to the emergence of new variants of the SARS-CoV-2 coronavirus, the question of how the viral genomes evolved, leading to the formation of highly infectious strains, becomes particularly important. Three major emergent strains, Alpha, Beta and Delta, characterized by a significant number of missense mutations, provide a natural test field. We accumulated and aligned 4.7 million SARS-CoV-2 genomes from the GISAID database and carried out a comprehensive set of analyses. This collection covers the period until the end of October 2021, i.e., the beginnings of the Omicron variant. First, we explored combinatorial complexity of the genomic variants emerging and their timing, indicating very strong, albeit hidden, selection forces. Our analyses show that the mutations that define variants of concern did not arise gradually but rather co-evolved rapidly, leading to the emergence of the full variant strain. To explore in more detail the evolutionary forces at work, we developed time trajectories of mutations at all 29,903 sites of the SARS-CoV-2 genome, week by week, and stratified them into trends related to (i) point substitutions, (ii) deletions and (iii) non-sequenceable regions. We focused on classifying the genetic forces active at different ranges of the mutational spectrum. We observed the agreement of the lowest-frequency mutation spectrum with the Griffiths–Tavaré theory, under the Infinite Sites Model and neutrality. If we widen the frequency range, we observe the site frequency spectra much more consistently with the Tung–Durrett model assuming clone competition and selection. The coefficients of the fitting model indicate the possibility of selection acting to promote gradual growth slowdown, as observed in the history of the variants of concern. These results add up to a model of genomic evolution, which partly fits into the classical drift barrier ideas. Certain observations, such as mutation “bands” persistent over the epidemic history, suggest contribution of genetic forces different from mutation, drift and selection, including recombination or other genome transformations. In addition, we show that a “toy” mathematical model can qualitatively reproduce how new variants (clones) stem from rare advantageous driver mutations, and then acquire neutral or disadvantageous passenger mutations which gradually reduce their fitness so they can be then outcompeted by new variants due to other driver mutations. 
    more » « less
  3. Abstract Background

    Four severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants predominated in the United States since 2021. Understanding disease severity related to different SARS-CoV-2 variants remains limited.

    Method

    Viral genome analysis was performed on SARS-CoV-2 clinical isolates circulating March 2021 through March 2022 in Cleveland, Ohio. Major variants were correlated with disease severity and patient outcomes.

    Results

    In total 2779 patients identified with either Alpha (n = 1153), Gamma (n = 122), Delta (n = 808), or Omicron variants (n = 696) were selected for analysis. No difference in frequency of hospitalization, intensive care unit (ICU) admission, and death were found among Alpha, Gamma, and Delta variants. However, patients with Omicron infection were significantly less likely to be admitted to the hospital, require oxygen, or admission to the ICU (χ2 = 12.8, P < .001; χ2 = 21.6, P < .002; χ2 = 9.6, P = .01, respectively). In patients whose vaccination status was known, a substantial number had breakthrough infections with Delta or Omicron variants (218/808 [26.9%] and 513/696 [73.7%], respectively). In breakthrough infections, hospitalization rate was similar regardless of variant by multivariate analysis. No difference in disease severity was identified between Omicron subvariants BA.1 and BA.2.

    Conclusions

    Disease severity associated with Alpha, Gamma, and Delta variants is comparable while Omicron infections are significantly less severe. Breakthrough disease is significantly more common in patients with Omicron infection.

     
    more » « less
  4. Kazarinoff, P. (Ed.)
    Disparities in undergraduate STEM degree completions across the United States are a national concern. Undergraduate-level research opportunities are vital for developing future researchers and building their scientific identity. These experiences can help students in community colleges acquire 21st-century skills and build confidence in their ability to do science [1-3]. The development and implementation of guided research experiences provide users with a topic they are familiar with but not necessarily experts in, like SARS-CoV2 infections. In this particular study, the Immune Epitope Database (IEDB) was used to identify amino acid residues located on the immunogenic regions of the spike glycoprotein of SARS-CoV-2 variants: Alpha, Beta, Gamma, Delta, and Omicron. IEDB is a web-based bioinformatics tool that contains published epitope information and prediction aids that can be used as a research platform for studying infectious diseases. The objective of this study aimed to map the immunogenic regions on the spike glycoproteins of the SARS-CoV-2 variants and predict the immune evasion of these variants [4-6]. Identifying the antigenic determinations that bind to the antibodies is essential for designing future candidates for peptide-based vaccines. This study aims to map the immunogenic regions on the spike glycoproteins of the SARS-CoV-2 variants and predict the immune evasion of these variants [4-6]. Identifying the antigenic determinations that bind to the antibodies is essential for designing future candidates for peptide-based vaccines. This research identifies regions where mutations have occurred in the virus, which are important to study as they can affect the virus’s immune evasion and impact available vaccines. Targeting multiple immunogenic regions unaffected by mutations can serve as potential targets for new vaccines, providing better protection against different variants. 
    more » « less
  5. Abstract

    Severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2), the virus causing COVID‐19, has continued to mutate and spread worldwide despite global vaccination efforts. In particular, the Omicron variant, first identified in South Africa in late November 2021, has become the dominant strain worldwide. Compared to the original strain identified in Wuhan, Omicron features 50 genetic mutations, with 15 mutations in the receptor‐binding domain (RBD) of the spike protein, which binds to the human angiotensin‐converting enzyme 2 (ACE2) receptor for viral entry. However, it is not completely understood how these mutations alter the interaction and binding strength between the Omicron RBD and ACE2. In this study, we used a combined steered molecular dynamics (SMD) simulation and experimental microscale thermophoresis (MST) approach to quantify the interaction between Omicron RBD and ACE2. We report that the Omicron brings an enhanced RBD‐ACE2 interface through N501Y, Q498R, and T478K mutations; the changes further lead to unique interaction patterns, reminiscing the features of previously dominated variants, Alpha (N501Y) and Delta (L452R and T478K). Among the Q493K and Q493R, we report that Q493R shows stronger binding to ACE2 than Q493K due to increased interactions. Our MST data confirmed that the Omicron mutations in RBD are associated with a five‐fold higher binding affinity to ACE2 compared to the RBD of the original strain. In conclusion, our results could help explain the Omicron variant's prevalence in human populations, as higher interaction forces or affinity for ACE2 likely promote greater viral binding and internalization, leading to increased infectivity.

     
    more » « less