Tracking SARS-CoV-2 genetic diversity is strongly indicated because diversifying selection may lead to the emergence of novel variants resistant to naturally acquired or vaccine-induced immunity. To monitor New York City (NYC) for the presence of novel variants, we deep sequence most of the receptor binding domain coding sequence of the S protein of SARS-CoV-2 isolated from the New York City wastewater. Here we report detecting increasing frequencies of novel cryptic SARS-CoV-2 lineages not recognized in GISAID’s EpiCoV database. These lineages contain mutations that had been rarely observed in clinical samples, including Q493K, Q498Y, E484A, and T572N and share many mutations with the Omicron variant of concern. Some of these mutations expand the tropism of SARS-CoV-2 pseudoviruses by allowing infection of cells expressing the human, mouse, or rat ACE2 receptor. Finally, pseudoviruses containing the spike amino acid sequence of these lineages were resistant to different classes of receptor binding domain neutralizing monoclonal antibodies. We offer several hypotheses for the anomalous presence of these lineages, including the possibility that these lineages are derived from unsampled human COVID-19 infections or that they indicate the presence of a non-human animal reservoir.

The first 2 years of the COVID-19 pandemic were mainly characterized by recurrent mutations of SARS-CoV-2 Spike protein at residues K417, L452, E484, N501 and P681 emerging independently across different variants of concern (Alpha, Beta, Gamma, and Delta). Such homoplasy is a marker of convergent evolution. Since Spring 2022 and the third year of the pandemic, with the advent of Omicron and its sublineages, convergent evolution has led to the observation of different lineages acquiring an additional group of mutations at different amino acid residues, namely R346, K444, N450, N460, F486, F490, Q493, and S494. Mutations at these residues have become increasingly prevalent during Summer and Autumn 2022, with combinations showing increased fitness. The most likely reason for this convergence is the selective pressure exerted by previous infection- or vaccine-elicited immunity. Such accelerated evolution has caused failure of all anti-Spike monoclonal antibodies, including bebtelovimab and cilgavimab. While we are learning how fast coronaviruses can mutate and recombine, we should reconsider opportunities for economically sustainable escape-proof combination therapies, and refocus antibody-mediated therapeutic efforts on polyclonal preparations that are less likely to allow for viral immune escape.

The severe acute respiratory syndrome coronavirus 2 (SARS-COV-2) Spike glycoprotein is solely responsible for binding to the host cell receptor and facilitating fusion between the viral and host membranes. The ability to generate viral particles pseudotyped with SARS-COV-2 Spike is useful for many types of studies, such as characterization of neutralizing antibodies or development of fusion-inhibiting small molecules. Here, we characterized the use of a codon-optimized SARS-COV-2 Spike glycoprotein for the generation of pseudotyped HIV-1, murine leukemia virus (MLV), and vesicular stomatitis virus (VSV) particles. The full-length Spike protein functioned inefficiently with all three systems but was enhanced over 10-fold by deleting the last 19 amino acids of the cytoplasmic tail. Infection of 293FT target cells was possible only if the cells were engineered to stably express the human angiotensin-converting enzyme 2 (ACE2) receptor, but stably introducing an additional copy of this receptor did not further enhance susceptibility. Stable introduction of the Spike-activating protease TMPRSS2 further enhanced susceptibility to infection by 5- to 10-fold. Replacement of the signal peptide of the Spike protein with an optimal signal peptide did not enhance or reduce infectious particle production. However, modifications D614G and R682Q further enhanced infectious particle production. With all enhancing elements combined, the titer of pseudotyped HIV-1 particles reached almost 10⁶ infectious particles/ml. Finally, HIV-1 particles pseudotyped with SARS-COV-2 Spike were successfully used to detect neutralizing antibodies in plasma from coronavirus disease 2019 (COVID-19) patients, but not in plasma from uninfected individuals. IMPORTANCE In work with pathogenic viruses, it is useful to have rapid quantitative tests for viral infectivity that can be performed without strict biocontainment restrictions. A common way of accomplishing this is to generate viral pseudoparticles that contain the surface glycoprotein from the pathogenic virus incorporated into a replication-defective viral particle that contains a sensitive reporter system. These pseudoparticles enter cells using the glycoprotein from the pathogenic virus, leading to a readout for infection. Conditions that block entry of the pathogenic virus, such as neutralizing antibodies, will also block entry of the viral pseudoparticles. However, viral glycoproteins often are not readily suited for generating pseudoparticles. Here, we describe a series of modifications that result in the production of relatively high-titer SARS-COV-2 pseudoparticles that are suitable for the detection of neutralizing antibodies from COVID-19 patients.

Wastewater-based epidemiology (WBE) is an effective way of tracking the appearance and spread of SARS-COV-2 lineages through communities. Beginning in early 2021, we implemented a targeted approach to amplify and sequence the receptor binding domain (RBD) of SARS-COV-2 to characterize viral lineages present in sewersheds. Over the course of 2021, we reproducibly detected multiple SARS-COV-2 RBD lineages that have never been observed in patient samples in 9 sewersheds located in 3 states in the USA. These cryptic lineages contained between 4 to 24 amino acid substitutions in the RBD and were observed intermittently in the sewersheds in which they were found for as long as 14 months. Many of the amino acid substitutions in these lineages occurred at residues also mutated in the Omicron variant of concern (VOC), often with the same substitutions. One of the sewersheds contained a lineage that appeared to be derived from the Alpha VOC, but the majority of the lineages appeared to be derived from pre-VOC SARS-COV-2 lineages. Specifically, several of the cryptic lineages from New York City appeared to be derived from a common ancestor that most likely diverged in early 2020. While the source of these cryptic lineages has not been resolved, it seems increasingly likely that they were derived from long-term patient infections or animal reservoirs. Our findings demonstrate that SARS-COV-2 genetic diversity is greater than what is commonly observed through routine SARS-CoV-2 surveillance. Wastewater sampling may more fully capture SARS-CoV-2 genetic diversity than patient sampling and could reveal new VOCs before they emerge in the wider human population.

Sequencing Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) from wastewater has become a useful tool in monitoring the spread of viral variants. Approaches to this task have been varied, relying on differing sequencing methods and computational analyses. We used a novel computation workflow based on amplicon sequencing of SARS-CoV-2 spike domains in order to track viral populations in wastewater. As part of this workflow, we developed a program, SAM Refiner, that has a variety of outputs, including novel variant reporting as well as functions designed to remove polymerase chain reaction (PCR) generated chimeric sequences. With these methods, we were able to track viral population dynamics over time. We report here on the emergence of two variants of concern, B.1.1.7 (Alpha) and P.1 (Gamma), and their displacement of the D614G B.1 variant in a Missouri sewershed.

Retrovirus assembly is driven by the multidomain structural protein Gag. Interactions between the capsid domains (CA) of Gag result in Gag multimerization, leading to an immature virus particle that is formed by a protein lattice based on dimeric, trimeric, and hexameric protein contacts. Among retroviruses the inter- and intra-hexamer contacts differ, especially in the N-terminal sub-domain of CA (CA_NTD). For HIV-1 the cellular molecule inositol hexakisphosphate (IP6) interacts with and stabilizes the immature hexamer, and is required for production of infectious virus particles. We have used in vitro assembly, cryo-electron tomography and subtomogram averaging, atomistic molecular dynamics simulations and mutational analyses to study the HIV-related lentivirus equine infectious anemia virus (EIAV). In particular, we sought to understand the structural conservation of the immature lentivirus lattice and the role of IP6 in EIAV assembly. Similar to HIV-1, IP6 strongly promoted in vitro assembly of EIAV Gag proteins into virus-like particles (VLPs), which took three morphologically highly distinct forms: narrow tubes, wide tubes, and spheres. Structural characterization of these VLPs to sub-4Å resolution unexpectedly showed that all three morphologies are based on an immature lattice with preserved key structural components, highlighting the structural versatility of CA to form immature assemblies. A direct comparison between EIAV and HIV revealed that both lentiviruses maintain similar immature interfaces, which are established by both conserved and non-conserved residues. In both EIAV and HIV-1, IP6 regulates immature assembly via conserved lysine residues within the CA_CTD and SP. Lastly, we demonstrate that IP6 stimulates in vitro assembly of immature particles of several other retroviruses in the lentivirus genus, suggesting a conserved role for IP6 in lentiviral assembly. The structural polyprotein Gag is conserved among all retroviruses and mediates virus assembly via oligomerization into incomplete lattices that are stabilized by dimeric, trimeric and hexameric contacts. Despite a high degree of conservation at the secondary and tertiary structure level, the quaternary interactions between the CA domains of retroviral Gag vary. Recently, the small cellular molecule IP6 was identified as an assembly co-factor of the lentivirus HIV-1. To better understand the structural variability of retroviruses and to determine if IP6 is an assembly cofactor of other lentiviruses we determined the structure of the HIV-1 related retrovirus EIAV. Using cryo-electron tomography and subtomogram averaging, in vitro assembly, mutation analysis, and molecular dynamic simulations, we determined and characterized the structure of the EIAV immature lattice. Furthermore, we found that IP6 is an assembly cofactor of EIAV, and other lentiviruses.