Continued selection on cryptic SARS-CoV-2 observed in Missouri wastewater

Results and Discussion

We use next generation sequencing of SARS-CoV-2’s RBD to monitor variants present in Missouri wastewater. Monitoring of MO45 began in March 2021 and continues to the present with roughly weekly sampling (Fig 1). Initially this sewershed was observed to primarily have the Alpha variant with some ancestral sequences. Beta, Gamma, Delta and Mu/Theta sequences were all observed later with Delta becoming the only variant detected by August 2021. Delta was then rapidly replaced by Omicron in December 2021. Since, various Omicron variants have circulated, generally with newer variants displacing older ones, resulting in a mixture of variants co-circulating in late 2022.

• Download figure
• Open in new tab

Fig 1. Presence of Variants Over Time.

Plot of relative defined SARS-CoV-2 variant abundance detected in MO45 sewershed. Wastewater samples were processed for viral RNA extraction and SARS-CoV-2 RBD sequences were amplified and then deep sequenced. Genotype and number of reads were used to estimate the relative abundance of defined variants present in the sewershed.

In addition to the defined variants, a cryptic lineage has also been sporadically detected, first in June 2021 and last in October 2022 (Fig 2). Initial sequences of this variant had K417T T478K E484A Q493K S494P Q498H amino acid changes relative to the ancestral sequence, with E484A and Q493K only appearing in one of the two first detections. E484A and Q493K were both observed in all subsequent sequences of this cryptic lineage, while S494P was not observed again. K417T and T478K had previously been observed in the Gamma and Delta variants respectively, but the other mutations had not yet appeared in any major VOCs.

• Download figure
• Open in new tab

Fig 2. Cryptic Haplotypes Detected in MO45.

Cryptic haplotypes are shown with the dates of their detection and relative abundance in the sample amplification. Select Omicron haplotypes are shown for comparison. Detection dates followed by ‘alt’ indicate the amplification was carried out with primers designed to exclude Omicron lineages for better detection of cryptic sequences.

Several amino acid changes occurred subsequent to the initial observation and appeared to become fixed in the lineage. On February 2, 2022 N460K was first observed in the cryptic lineage and was thereafter fixed. Likewise S477N and F486V were first observed in the cryptic lineage on April 5, 2022, and N440K on April 26, 2022, and in all detections since. N450D was first observed in the cryptic lineage on May 24, 2022. Though the lineage had N450Y on June 16, 2022 instead, the two subsequent detections of the lineage had N450D again. Several other changes were observed in the lineage over time, though none could be concluded to have become fixed. Of note, T547I and T572I both were observed in 3 samples each, but were not observed in the most recent detection. At the last detection of the MO45 cryptic lineage in October 2022, the lineage appeared remarkably similar to an Omicron lineage with 12 amino acid changes in its RBD that were all identical to, or at the same position as, changes found in Omicron lineages.

Most of the residue changes observed in the cryptic lineage predate the changes observed in Omicron. The initial detections of the cryptic lineage, months before the emergence of Omicron, already had two changes that were to be seen ubiquitously in Omicron lineages, T478K and E484A, and three changes at the same residues as changes common in Omicron lineages, K417T, Q493K and Q498H. Likewise, N460K, which appeared in the cryptic lineage in February of 2022, did not become prevalent in an Omicron background until six months later.

The convergence of the cryptic lineage and Omicron variants suggest similar selection pressures. The origin of Omicron and the origin of the MO45 cryptic lineage are unknown. At least in some cases, cryptic lineages appear to be derived from individuals with chronic SARS-CoV-2 infections. However, as the MO45 cryptic lineage hasn’t been traced, a non-human source can not be ruled out. Since the cryptic lineage in some cases acquired changes prior to Omicron, continued monitoring of waste water for such cryptic lineages may provide foreknowledge of changes, or at least the position of changes, likely to be selected for in the circulating Omicron variants.

Sample collection and RNA extraction

Collection and processing of samples were as previously described. Twenty-four-hour composite samples were collected at the MO45 wastewater treatment facility and maintained at 4 °C until they were delivered to the analysis lab, generally within 24 h of collection. Samples were then centrifuged at 3000× g for 10 min and followed by filtration through a 0.22 µM polyethersolfone membrane (Millipore, Burlington, MA, USA). Approximately 37.5 mL of wastewater was mixed with 12.5 mL solution containing 50% (w/vol) polyethylene glycol 8000 and 1.2 M NaCl, mixed, and incubated at 4 °C for at least 1 h. Samples were then centrifuged at 12,000× g for 2 h at 4 °C. Supernatant was decanted and RNA was extracted from the remaining pellet (usually not visible) with the QIAamp Viral RNA Mini Kit (Qiagen, Germantown, MD, USA) using the manufacturer’s instructions. RNA was extracted in a final volume of 60 µL. MiSeq Similar to our previous protocol, the primary RBD RT-PCR was performed using the Superscript IV One-Step RT-PCR System (Thermo Fisher Scientific,12594100). Primary RT-PCR amplification was performed as follows: 25°C (2:00) + 50°C (20:00) + 95°C (2:00) + [95°C (0:15) + 55°C (0:30) + 72°C (1:00)] × 25 cycles using the MiSeq primary PCR primers CTGCTTTACTAATGTCTATGCAGATTC and NCCTGATAAAGAACAGCAACCT. Secondary PCR (25 μL) was performed on RBD amplifications using 5 μL of the primary PCR as template with MiSeq nested gene specific primers containing 5′ adapter sequences (0.5 μM each) acactctttccctacacgacgctcttccgatctGTRATGAAGTCAGMCAAATYGC and gtgactggagttcagacgtgtgctcttccgatctATGTCAAGAATCTCAAGTGTCTG, dNTPs (100 μM each) (New England Biolabs, N0447L) and Q5 DNA polymerase (New England Biolabs, M0541S).

Secondary PCR amplification was performed as follows: 95°C (2:00) + [95°C (0:15) + 55°C (0:30) + 72°C (1:00)] × 20 cycles.

For Omicron exclusion amplification, the primary RBD RT-PCR was performed using the MiSeq primary PCR primers ATTCTGTCCTATATAATTCCGCAT and CCCTGATAAAGAACAGCAACCT (the first primer was changed to TATATAATTCCGCATCATTTTCCAC starting in May, 2022 to adapt to changing Omicron lineages) and secondary PCR used MiSeq nested gene specific primers containing 5′ adapter sequences (0.5 μM each) acactctttccctacacgacgctcttccgatctGTGATGAAGTCAGACAAATCGC and gtgactggagttcagacgtgtgctcttccgatctATGTCAAGAATCTCAAGTGTCTG.

A tertiary PCR (50 μL) was performed to add adapter sequences required for Illumina cluster generation with forward and reverse primers (0.2 μM each), dNTPs (200 μM each) (New England Biolabs, N0447L) and Phusion High-Fidelity or (KAPA HiFi for CA samples) DNA Polymerase (1U) (New England Biolabs, M0530L). PCR amplification was performed as follows: 98°C (3:00) + [98°C (0:15) + 50°C (0:30) + 72°C (0:30)] × 7 cycles +72°C (7:00). Amplified product (10 μl) from each PCR reaction is combined and thoroughly mixed to make a single pool. Pooled amplicons were purified by addition of Axygen AxyPrep MagPCR Clean-up beads (Axygen, MAG-PCR-CL-50) or in a 1.0 ratio to purify final amplicons. The final amplicon library pool was evaluated using the Agilent Fragment Analyzer automated electrophoresis system, quantified using the Qubit HS dsDNA assay (Invitrogen), and diluted according to Illumina’s standard protocol. The Illumina MiSeq instrument was used to generate paired-end 300 base pair reads. Adapter sequences were trimmed from output sequences using Cutadapt.