Taking stock of the mutations in human SARS-CoV-2 spike proteins: From early days to nearly the end of COVID-19 pandemic

In depth sequencing of a serially sampled household cohort reveals the within-host dynamics of Omicron SARS-CoV-2 and rare selection of novel spike variants

SARS-CoV-2 has undergone repeated and rapid evolution to circumvent host immunity. However, outside of prolonged infections in immunocompromised hosts, within-host positive selection has rarely been detected. Here we combine daily longitudinal sampling of individuals with replicate sequencing to increase the accuracy of and lower the threshold for variant calling. We sequenced 577 specimens from 105 individuals in a household cohort during the BA.1/BA.2 variant period. Individuals exhibited extremely low viral diversity, and we estimated a low within-host evolutionary rate. Within-host dynamics were dominated by genetic drift and purifying selection. Positive selection was rare but highly concentrated in spike. A Wright Fisher Approximate Bayesian Computational model identified positive selection at 14 loci with 7 in spike, including S:448 and S:339. This detectable immune-mediated selection is unusual in acute respiratory infections and may be caused by the relatively narrow antibody repertoire in individuals during the early Omicron phase of the SARS-CoV-2 pandemic.

Conformational and Stability Analysis of SARS-CoV-2 Spike Protein Variants by Molecular Simulation

We performed a comprehensive structural analysis of the conformational space of several spike (S) protein variants using molecular dynamics (MD) simulations. Specifically, we examined four well-known variants (Delta, BA.1, XBB.1.5, and JN.1) alongside the wild-type (WT) form of SARS-CoV-2. The conformational states of each variant were characterized by analyzing their distributions within a selected space of collective variables (CVs), such as inter-domain distances between the receptor-binding domain (RBD) and the N-terminal domain (NTD). Our primary focus was to identify conformational states relevant to potential structural transitions and to determine the set of native contacts (NCs) that stabilize these conformations. The results reveal that genetically more distant variants, such as XBB.1.5, BA.1, and JN.1, tend to adopt more compact conformational states compared to the WT. Additionally, these variants exhibit novel NC profiles, characterized by an increased number of specific contacts distributed among ionic, polar, and nonpolar residues. We further analyzed the impact of specific mutations, including T478K, N500Y, and Y504H. These mutations not only enhance interactions with the human host receptor but also alter inter-chain stability by introducing additional NCs compared to the WT. Consequently, these mutations may influence the accessibility of certain protein regions to neutralizing antibodies. Overall, these findings contribute to a deeper understanding of the structural and functional variations among S protein variants.

Phylogeogenomic analysis of the earliest reported sequences of SARS-CoV-2 from 161 countries

The SARS-CoV-2 is the causative agent of COVID-19 whose evolutionary path with geographical context forms the focus of present study. The first reported sequence from each of the 161 countries was downloaded from the GISAID database. Multiple sequence alignment was performed using MAFFT v.7, and a TCS-based network was constructed using PopART v.1.7. A total of 27 proteins were analyzed including structural and non-structural proteins. NSP3 and NSP12, responsible for viral replication and RNA synthesis, respectively, had the highest mutation incidence and frequency among non-structural proteins. The spike (S) protein, critical for viral attachment and entry, had the highest prevalence and frequency of mutations. ORF3a had the highest mutation incidence and frequency among accessory proteins. The phylogeogenomic network identified six haplogroups containing 35 sequences, while the remaining sequences belonged to different haplotypes. The virus's genetic distinctiveness was higher in European genomes, with four haplogroups dominated by Europe-linked sequences. The triangular-shaped pattern observed in the virus's evolutionary path suggests that it spread to different continents from Asia. Multiple transmission pathways connecting different countries affirm the virus's ability to emerge in multiple countries by early 2020. The possibility of new species emergence through "saltation" due to the pandemic is also discussed.

In depth sequencing of a serially sampled household cohort reveals the within-host dynamics of Omicron SARS-CoV-2 and rare selection of novel spike variants

SARS-CoV-2 has undergone repeated and rapid evolution to circumvent host immunity. However, outside of prolonged infections in immunocompromised hosts, within-host positive selection has rarely been detected. The low diversity within-hosts and strong genetic linkage among genomic sites make accurately detecting positive selection difficult. Longitudinal sampling is a powerful method for detecting selection that has seldom been used for SARS-CoV-2. Here we combine longitudinal sampling with replicate sequencing to increase the accuracy of and lower the threshold for variant calling. We sequenced 577 specimens from 105 individuals from a household cohort primarily during the BA.1/BA.2 variant period. There was extremely low diversity and a low rate of divergence. Specimens had 0-12 intrahost single nucleotide variants (iSNV) at >0.5% frequency, and the majority of the iSNV were at frequencies <2%. Within-host dynamics were dominated by genetic drift and purifying selection. Positive selection was rare but highly concentrated in spike. Two individuals with BA.1 infections had S:371F, a lineage defining substitution for BA.2. A Wright Fisher Approximate Bayesian Computational model identified positive selection at 14 loci with 7 in spike, including S:448 and S:339. We also detected significant genetic hitchhiking between synonymous changes and nonsynonymous iSNV under selection. The detectable immune-mediated selection may be caused by the relatively narrow antibody repertoire in individuals during the early Omicron phase of the SARS-CoV-2 pandemic. As both the virus and population immunity evolve, understanding the corresponding shifts in SARS-CoV-2 within-host dynamics will be important.

Differential epitope prediction across diverse circulating variants of SARS-COV-2 in Brazil

COVID-19, caused by the SARS-COV-2 virus, induces numerous immunological reactions linked to the severity of the clinical condition of those infected. The surface Spike protein (S protein) present in Sars-CoV-2 is responsible for the infection of host cells. This protein presents a high rate of mutations, which can increase virus transmissibility, infectivity, and immune evasion. Therefore, we propose to evaluate, using immunoinformatic techniques, the predicted epitopes for the S protein of seven variants of Sars-CoV-2. MHC class I and II epitopes were predicted and further assessed for their immunogenicity, interferon-gamma (IFN-γ) inducing capacity, and antigenicity. For B cells, linear and structural epitopes were predicted. For class I MHC epitopes, 40 epitopes were found for the clades of Wuhan, Clade 2, Clade 3, and 20AEU.1, Gamma, and Delta, in addition to 38 epitopes for Alpha and 44 for Omicron. For MHC II, there were differentially predicted epitopes for all variants and eight equally predicted epitopes. These were evaluated for differences in the MHC II alleles to which they would bind. Regarding B cell epitopes, 16 were found in the Wuhan variant, 14 in 22AEU.1 and in Clade 3, 15 in Clade 2, 11 in Alpha and Delta, 13 in Gamma, and 9 in Omicron. When compared, there was a reduction in the number of predicted epitopes concerning the Spike protein, mainly in the Delta and Omicron variants. These findings corroborate the need for updates seen today in bivalent mRNA vaccines against COVID-19 to promote a targeted immune response to the main circulating variant, Omicron, leading to more robust protection against this virus and avoiding cases of reinfection. When analyzing the specific epitopes for the RBD region of the spike protein, the Omicron variant did not present a B lymphocyte epitope from position 390, whereas the epitope at position 493 for MHC was predicted only for the Alpha, Gamma, and Omicron variants.

Development of a highly stable, active small interfering RNA with broad activity against SARS-CoV viruses

Treatment options for COVID-19 remain limited. Here, we report the optimization of an siRNA targeting the highly conserved leader region of SARS-CoV-2. The siRNA was rendered nuclease resistant by the introduction of modified nucleotides without loss of activity. Importantly, the siRNA also retained its inhibitory activity against the emerged omicron sublineage variant BA.2, which occurred after the siRNA was designed and is resistant to other antiviral agents such as antibodies. In addition, we show that a second highly active siRNA designed against the viral 5'-UTR can be applied as a rescue molecule, to minimize the spread of escape mutations. We therefore consider our siRNA-based molecules to be promising broadly active candidates for the treatment of current and future SARS-CoV-2 variants.