Occurrence of a substitution or deletion of SARS-CoV-2 spike amino acid 677 in various lineages in Marseille, France

Within-host evolutionary dynamics and tissue compartmentalization during acute SARS-CoV-2 infection

The global evolution of SARS-CoV-2 depends in part upon the evolutionary dynamics within individual hosts with varying immune histories. To characterize the within-host evolution of acute SARS-CoV-2 infection, we sequenced saliva and nasal samples collected daily from vaccinated and unvaccinated individuals early during infection. We show that longitudinal sampling facilitates high-confidence genetic variant detection and reveals evolutionary dynamics missed by less-frequent sampling strategies. Within-host dynamics in both unvaccinated and vaccinated individuals appeared largely stochastic; however, in rare cases, minor genetic variants emerged to frequencies sufficient for forward transmission. Finally, we detected significant genetic compartmentalization of viral variants between saliva and nasal swab sample sites in many individuals. Altogether, these data provide a high-resolution profile of within-host SARS-CoV-2 evolutionary dynamics.IMPORTANCEWe detail the within-host evolutionary dynamics of SARS-CoV-2 during acute infection in 31 individuals using daily longitudinal sampling. We characterized patterns of mutational accumulation for unvaccinated and vaccinated individuals, and observed that temporal variant dynamics in both groups were largely stochastic. Comparison of paired nasal and saliva samples also revealed significant genetic compartmentalization between tissue environments in multiple individuals. Our results demonstrate how selection, genetic drift, and spatial compartmentalization all play important roles in shaping the within-host evolution of SARS-CoV-2 populations during acute infection.

SPLASH: A statistical, reference-free genomic algorithm unifies biological discovery

Today's genomics workflows typically require alignment to a reference sequence, which limits discovery. We introduce a unifying paradigm, SPLASH (Statistically Primary aLignment Agnostic Sequence Homing), which directly analyzes raw sequencing data, using a statistical test to detect a signature of regulation: sample-specific sequence variation. SPLASH detects many types of variation and can be efficiently run at scale. We show that SPLASH identifies complex mutation patterns in SARS-CoV-2, discovers regulated RNA isoforms at the single-cell level, detects the vast sequence diversity of adaptive immune receptors, and uncovers biology in non-model organisms undocumented in their reference genomes: geographic and seasonal variation and diatom association in eelgrass, an oceanic plant impacted by climate change, and tissue-specific transcripts in octopus. SPLASH is a unifying approach to genomic analysis that enables expansive discovery without metadata or references.

[WITHDRAWN] SPLASH: a statistical, reference-free genomic algorithm unifies biological discovery

The authors have withdrawn this manuscript due to a duplicate posting of manuscript number BIORXIV/2022/497555. Therefore, the authors do not wish this work to be cited as reference for the project. If you have any questions, please contact the corresponding author. The correct preprint can be found at doi: https://doi.org/10.1101/2022.06.24.497555.

O-Linked Sialoglycans Modulate the Proteolysis of SARS-CoV-2 Spike and Likely Contribute to the Mutational Trajectory in Variants of Concern

The emergence of a polybasic cleavage motif for the protease furin in SARS-CoV-2 spike has been established as a major factor for human viral transmission. The region N-terminal to that motif is extensively mutated in variants of concern (VOCs). Besides furin, spikes from these variants appear to rely on other proteases for maturation, including TMPRSS2. Glycans near the cleavage site have raised questions about proteolytic processing and the consequences of variant-borne mutations. Here, we identify that sialic acid-containing O-linked glycans on Thr678 of SARS-CoV-2 spike influence furin and TMPRSS2 cleavage and posit O-linked glycosylation as a likely driving force for the emergence of VOC mutations. We provide direct evidence that the glycosyltransferase GalNAc-T1 primes glycosylation at Thr678 in the living cell, an event that is suppressed by mutations in the VOCs Alpha, Delta, and Omicron. We found that the sole incorporation of N-acetylgalactosamine did not impact furin activity in synthetic O-glycopeptides, but the presence of sialic acid reduced the furin rate by up to 65%. Similarly, O-glycosylation with a sialylated trisaccharide had a negative impact on TMPRSS2 cleavage. With a chemistry-centered approach, we substantiate O-glycosylation as a major determinant of spike maturation and propose disruption of O-glycosylation as a substantial driving force for VOC evolution.

Computational Analysis Predicts Correlations among Amino Acids in SARS-CoV-2 Proteomes

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a serious global challenge requiring urgent and permanent therapeutic solutions. These solutions can only be engineered if the patterns and rate of mutations of the virus can be elucidated. Predicting mutations and the structure of proteins based on these mutations have become necessary for early drug and vaccine design purposes in anticipation of future viral mutations. The amino acid composition (AAC) of proteomes and individual viral proteins provide avenues for exploitation since AACs have been previously used to predict structure, shape and evolutionary rates. Herein, the frequency of amino acid residues found in 1637 complete proteomes belonging to 11 SARS-CoV-2 variants/lineages were analyzed. Leucine is the most abundant amino acid residue in the SARS-CoV-2 with an average AAC of 9.658% while tryptophan had the least abundance of 1.11%. The AAC and ranking of lysine and glycine varied in the proteome. For some variants, glycine had higher frequency and AAC than lysine and vice versa in other variants. Tryptophan was also observed to be the most intolerant to mutation in the various proteomes for the variants used. A correlogram revealed a very strong correlation of 0.999992 between B.1.525 (Eta) and B.1.526 (Iota) variants. Furthermore, isoleucine and threonine were observed to have a very strong negative correlation of -0.912, while cysteine and isoleucine had a very strong positive correlation of 0.835 at p < 0.001. Shapiro-Wilk normality test revealed that AAC values for all the amino acid residues except methionine showed no evidence of non-normality at p < 0.05. Thus, AACs of SARS-CoV-2 variants can be predicted using probability and z-scores. AACs may be beneficial in classifying viral strains, predicting viral disease types, members of protein families, protein interactions and for diagnostic purposes. They may also be used as a feature along with other crucial factors in machine-learning based algorithms to predict viral mutations. These mutation-predicting algorithms may help in developing effective therapeutics and vaccines for SARS-CoV-2.

QTQTN motif upstream of the furin-cleavage site plays a key role in SARS-CoV-2 infection and pathogenesis

The furin cleavage site (FCS), an unusual feature in the SARS-CoV-2 spike protein, has been spotlighted as a factor key to facilitating infection and pathogenesis by increasing spike processing. Similarly, the QTQTN motif directly upstream of the FCS is also an unusual feature for group 2B coronaviruses (CoVs). The QTQTN deletion has consistently been observed in in vitro cultured virus stocks and some clinical isolates. To determine whether the QTQTN motif is critical to SARS-CoV-2 replication and pathogenesis, we generated a mutant deleting the QTQTN motif (ΔQTQTN). Here, we report that the QTQTN deletion attenuates viral replication in respiratory cells in vitro and attenuates disease in vivo. The deletion results in a shortened, more rigid peptide loop that contains the FCS and is less accessible to host proteases, such as TMPRSS2. Thus, the deletion reduced the efficiency of spike processing and attenuates SARS-CoV-2 infection. Importantly, the QTQTN motif also contains residues that are glycosylated, and disruption of its glycosylation also attenuates virus replication in a TMPRSS2-dependent manner. Together, our results reveal that three aspects of the S1/S2 cleavage site-the FCS, loop length, and glycosylation-are required for efficient SARS-CoV-2 replication and pathogenesis.