High-Throughput Adaptable SARS-CoV-2 Screening for Rapid Identification of Dominant and Emerging Regional Variants ================================================================================================================== * Zita Hubler * Xiao Song * Cameron Norris * Mehul Jani * David Alouani * Maureen Atchley * Lisa Stempak * Sarah Cherian * Christine Schmotzer * Navid Sadri ## ABSTRACT **Objectives** Emerging SARS-CoV-2 variant strains can be associated with increased transmissibility, more severe disease, and reduced effectiveness of treatments. To improve the availability of regional variant surveillance, we describe a variant genotyping system that is rapid, accurate, adaptable, and able to detect new low-level variants built with existing hospital infrastructure. **Methods** We use a tiered high-throughput SARS-CoV-2 screening program to characterizes variants in a supra-regional health system over 76 days. Combining targeted qPCR and selective sequencing, we screen positive SARS-CoV-2 samples from all hospitals within our health care system for genotyping dominant and emerging variants. **Results** The median turnaround for genotyping was two days using the high-throughput qPCR-based screen, allowing us to rapidly characterize the emerging Delta variant. In our population, the Delta variant is associated with a lower CT value, lower age at infection, and increased vaccine breakthrough cases. Detection of low-level and potentially emerging variants highlights the utility of a tiered approach. **Conclusions** These findings underscore the need for fast, low-cost, high-throughput monitoring of regional viral sequences as the pandemic unfolds and the emergence of SARS-CoV-2 variants increases. Combing qPCR-based screening with selective sequencing allows for rapid genotyping of variants and dynamic system improvement. **Key messages** * A tiered approach that uses qPCR-based screening to identify dominant variants and sequencing for unique variants maximizes throughput, turnaround time, and information gleaned from each sample. * In our population, the Delta variant became dominant in less than a month and is associated with lower CT, lower age at infection, and increased breakthrough cases. * We identified low-level variants, including the variant of interest B.1.621 and a Delta variant with an E484K mutation in our population using existing hospital infrastructure. Keywords * SARS-CoV-2 * variant identification * Delta variant * targeted testing * variant of concern * variant screening * E484K ## INTRODUCTION Over the last year, new variants of SARS-CoV-2 have emerged from all corners of the globe and some are associated with increased virulence, differing response to treatment, and the ability to evade vaccines. [1] New variants, combined with increased travel and variable vaccination status, can abruptly strain hospital resources.[2, 3] There is a need for quick identification and characterization of both dominant and low-level SARS-CoV-2 variants on a region-per-region basis to allow for appropriate allocation of resources in preparation for waxing and waning infection rates. Further, this identification system needs to be dynamic to monitor newly emerging variants that may become the dominant strain in a matter of weeks. While it would be ideal, use of sequencing for variant identification can be cost-prohibitive, low-throughput, and not broadly available in many clinical laboratories. By contrast, qPCR is available to most clinical laboratories and represents a speedy and cost-effective alternative to genotype common SARS-CoV-2 variants.[4–7] Several versions of qPCR-based variants screens detect common genetic alterations in the SARS-CoV-2 sequence to differentiate between variants. [4–7] Yet, qPCR-based screens are limited in their availability to detect newly emerging variants that are not identifiable with the existing genomic signatures. This work describes a tiered approach to SARS-CoV-2 variant screening that combines qPCR-based screening with the sequencing of select samples to monitor known variants and detect newly emerging variants quickly, accurately, and with flexible infrastructure investment. ## METHODS ### Samples SARS-Cov-2 positive nasal or nasopharyngeal samples identified via routine clinical testing at one of 12 regional hospitals between May 21, 2021 and August 4, 2021 were included in the variant screening program. The program was developed to support regional and state variant surveillance efforts. The routine clinical testing platforms varied between regional hospitals and included: ID Now COVID-19 (Abbott), Simplexa COVID-19 Direct Kit (DiaSorin), Xpert Xpress SARS-CoV-2 (Cepheid), TaqPath COVID-19 Combo Kit (ThermoFisher Scientific), and Aptima SARS-CoV-2 (Hologic) Assays. ### RNA extraction and qPCR All regional hospitals were instructed to send positive samples to a central lab performing the TaqPath COVID-19 Combo Kit assay to enable processing in a high-throughput 96 well plate format. RNA extraction and real-time qPCRs were performed per manufacturer’s instructions in a semi-automated fashion using a KingerFisher Flex System (ThermoFisher Scientific) and TaqPath COVID-19 Combo Kit (ThermoFisher Scientific) to assess for sufficient levels of ORF1ab, N-gene, and S-gene. The RNA plates were then transported to the sequencing laboratory to perform additional multiplexed qPCRs for regions of interest: ORF1del, N501Y, E484K, L452R. (For primers and more detailed methods, see Supplemental Table 1.) Data from all PCR plates were combined and integrated with patient demographic data. A deep learning model was used to identify any instances of errant cycle-threshold (Ct) calls. [4] Samples with N gene Ct greater than 33 were deemed to have an insufficient viral load for variant calling. Vaccine breakthrough cases are defined as having occurred in patients who completed their vaccination series within the hospital system more than two weeks before the sample collection date. ### Variant preliminary calling A preliminary probability matrix was calculated using publically available data that correlated the detection of specific viral genome markers with major WHO variants and CDC Variants of Concern (VoC)/Variants of Interest (VoI) using the primers described. An example of a probability matrix for publically available data is shown in Supplemental Table 2. The preliminary probability matrix was the basis for a ruleset that guides an automated preliminary calling algorithm. (Supplemental Table 3) The qPCR gene profile of each sample and the initial automated variant call are quickly reviewed and finalized. (Figure 1) Samples that have unclear expression profiles can be sequenced for further variant identification. (Figure 1) ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/08/2021.08.31.21262625/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/09/08/2021.08.31.21262625/F1) Figure 1: Schematic illustrating Dynamic Variant Calling Pipeline. Gray arrows indicate optional pathways, with examples of samples chosen for sequencing. Arrows and boxes coming from variant identification boxes show areas for iterative improvement and quality control of the pipeline. Confirmed positives are samples that have an N-gene Ct of less than 33. ### Sequencing Samples were sequenced for one of three reasons: unclear expression profiles, variant sampling (pipeline validation), or clinically interesting cases (such as vaccine breakthrough). (Figure 1) The number of samples sequenced and sequencing timing is based on availability of clinical laboratory personnel and resources, and case numbers needed to comprise a sequencing batch. We performed whole-genome sequencing using the SARS-CoV-2 Research Panel (Ion Gene Studio S5 system, ThermoFisher Scientific) as previously described. [8] Pangolin ([pangolin.cog-uk.io](http://pangolin.cog-uk.io), v.3.1.3, 8/6/2021) and NextClad ([clades.nextstrain.org](http://clades.nextstrain.org), v.1.5.4) were used for lineage or clade assignment, respectively. Sequencing data was submitted to GeneBank ([ncbi.nlm.nih.gov/genbank](http://ncbi.nlm.nih.gov/genbank)) for public use. ### Data Analysis For this study, variant call from the qPCR-based Variant Screen and the sequencing results were combined and analyzed using Excel, R, and GraphPad Prism. Statistics were performed using GraphPad Prism. Data for the Delta+E484K samples were downloaded from GISAID ([www.gisaid.org/](http://www.gisaid.org/)) on 8/18/21. Data for the Cuyahoga County vaccination rates were downloaded from the Ohio Department of Health COVID-19 Dashboard ([coronavirus.ohio.gov/wps/portal/gov/covid-19/dashboards](http://coronavirus.ohio.gov/wps/portal/gov/covid-19/dashboards)) on 8/6/21. ## RESULTS ### Pipeline Performance Between May 21, 2021 and August 4, 2021, 1150 positive COVID samples were processed in the qPCR-Based Screening Program, representing 2.0% of all COVID tests performed. (Figure 2) 99.5% of the screening results had an identifiable variant expression profile. Of all samples screened, 167 were subsequently sequenced for verification or clarification of variant calling. (Figure 2) Of note, based on our ruleset, the algorithm is not able to differentiate Pangolin lineages of the same variant, such as B.1.617.2 versus AY.3. However, it correctly classified the overall WHO variant as Delta. More granular identification of Pangolin lineages can be accomplished by incorporating novel primers and adjusting the ruleset (Figure 1). As part of our iterative improvement, we monitor the validity of the ruleset given the current variant distribution in our population. In the sequencing confirmed samples, for each variant called in the screen, the identification was correct 100% of the time (Figure 2). These results show that our qPCR screening is accurate and high-throughput, a necessary feature given the increasing rate of SARS-CoV-2 infections. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/08/2021.08.31.21262625/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2021/09/08/2021.08.31.21262625/F2) Figure 2: Flowchart of pipeline throughput over 76 days. Data collected from 5-21-21 to 08-04-21. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/08/2021.08.31.21262625/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2021/09/08/2021.08.31.21262625/F3) Figure 3: Tracking of Variants and Characterization of Delta A) Percentage of each variant per total samples tested per day. Each bar corresponds to one day. Variants are color-coded and indicated in the legends. The totals indicated in the legend are the total count detected via the pipeline over 76 days. “Low” indicates samples with insufficient viral load for screening as measured by a high N gene Ct>33. “Not Screened” samples were not received at the main facility. B) Each dot indicates one sample, the y-axis is patient’s age, x-axis is the variant identified. Median and 95% confidence interval indicated with the bars. **, P<0.005, Mann Whitney two-tailed test. C) Variant versus Cycle threshold (Ct) for core markers: ORF1ab, N gene, and S gene. Each dot represents one sample. Median and 95% confidence interval indicated. \**\*|\*, P<0.0001, One-way Anova with Tukey multiple comparison correction. D) Percentage of each variant that is in patients known to be vaccinated relative to all patients identified with that variant. \**\*|\*, P < 0.0005, two-tailed Fisher exact test. **, P< 0.005, two-tailed Fisher exact test. Another essential feature of SARS-CoV-2 variant surveillance is quick turnaround time. We found the median lag time before screening is two days after sample collection, with a degree of variability which may be attributable to differing transportation times from various facilities. (Supplemental Fig. 1A) By contrast, time to sequencing fluctuates, and may be delayed due to practical limitations such as personnel and equipment availability, as well as the number of available samples for a batch. Consequently, the median time between sample collection and sequencing was 14 days, substantially longer than the initial screen. (Supplemental Fig. 1A) This discrepancy in median screening versus sequencing time underscores the advantages of a rapid preliminary screen paired with subsequent sequencing. ### Variability in Variant Distribution over time We analyzed the data for the variants called within the last 76 days, as this timeframe represents a relative plateau of vaccination initiation and completion in our county, reducing the temporal effects of changing vaccination rates in our community on the analyses performed. (Supplemental Fig. 1 B). We noted that in less than a month, the dominant strains shifted from being Alpha and Gamma to Delta. (Figure 2A) The transition rate between variants highlights the utility of a rapid variant calling algorithm. Overall, this data parallels the international trend; after the wave of Alpha swept our county, Delta arrived and was associated with an increase in cases that dominated all other strains in our region.[9, 10] ### Characterization of Delta Variant Given the ongoing concern for Delta’s increasing viral load, infection of younger patients, and vaccine breakthrough cases, we wanted to characterize the impact of the Delta variant in our population.[9–13] We found a statistically significant decrease in the age of the patients testing positive for Delta compared to Alpha. (Figure 2B). Further, we found that Delta samples had cycle thresholds (Ct) significantly lower than Alpha and Gamma for the three SARS-CoV-2 markers assessed. (Figure 2C)[14] This may suggest an association between the Delta variant and increased viral load, as Ct has been used as an imperfect surrogate measure for viral load. Lastly, there was a statistically significant difference in the proportion of patients with each variant who were known to be vaccinated. [10] However, we noted a temporal association between the Alpha variant being dominant and the occurrence of low viral load vaccine-breakthrough events that were not genotyped. This suggested that difference in Ct values between Alpha and Delta could be confounding the analysis of vaccine-breakthrough. (Supplemental Fig. 1 C) Therefore, we assumed that all high Ct breakthrough cases were due to Alpha and Gamma, and despite this, Delta is associated with more vaccinated breakthrough cases. Together these results suggest that in our population, the Delta strain is infecting younger patients, with higher viral load, and more vaccine breakthrough cases. ### Identification of Novel Delta Variant One advantage of our process is the ability to selectively use sequencing for samples that are most likely to be informative. While variants aside from Alpha, Gamma, and Delta, represented less than 3.6% of identified cases, these low-level variants can provide clues about the next more virulent or endemic strain. [2] (Figure 2A) Unbiased sample sequencing, especially when performed in limited numbers, may miss low-level variants. Our high-throughput qPCR screen identified some emerging low-level variants, such as 5 cases of B.1.621 (Figure 2A). However, novel or unique low-level variants require biased sequencing for further identification [2]. Our identification of a Delta variant with an E484K mutation exemplifies biased sequencing of samples with ambiguous qPCR screening results as a mechanism for identifying low-level variants. The sample had an unclear mutational signature on our qPCR screen that included an E484K and an L452R mutation. Within two days, we used sequencing and clustering to elucidate that the sample is a Delta variant based on its Pangolin lineage B.1.617.2 (ambiguity score: 0.9995) and Nextclade clade 21A. (Figure 4A) Sequencing confirmed the presence of an E484K mutation which is surprising because, this mutation is detected in less than 0.03% of Delta variants in GISAID ([gisaid.org](http://gisaid.org)). Delta containing an E484K mutation is not prevalent in the USA but was detected in 56 Turkish samples, suggesting it may spread regionally and warrants additional surveillance. (Figure 4B) Interestingly, when phylogenetically arranged by Nextclade, our sample did not cluster with the samples from Turkey, making direct spread from Turkey less likely. (Supplemental Fig. 2). Instead was most similar to samples from the USA and Russia with, or without, an E484K mutation. (Supplemental Fig. 2) Overall, our detection of an uncommon Delta variant with an E484K mutation highlights the strength for monitoring low-level variants using a tiered qPCR and sequencing approach. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/08/2021.08.31.21262625/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2021/09/08/2021.08.31.21262625/F4) Figure 4: Identification of Delta + E484K A) Phylogenetic tree showing the clade of the unknown sample (shown with the star). Radial axis is the number of mutations with Delta, Gamma, and Alpha highlighted. B) Number of reported cases of Delta +E484K in GISAID per country. ## DISCUSSION As the COVID pandemic unfolds, the SARS-COV2 variants are evolving to be increasingly transmittable, have varying responses to treatment, and higher risk of evading existing vaccines. [1, 11] Information about newly emerging variants is vital for public health and hospital decision-making as they can overwhelm regional facilities and resources. [15] Therefore, healthcare systems need rapid and cost-effective variant detection to differentiate between subtly distinct and newly emerging variants. While no one molecular technique is ideal for all facilities, combining technologies in a tiered approach that combines qPCR screening and whole-genome sequencing allows for quick and accurate variant detection that is amenable to optimization. qPCR-based testing is well-suited for identifying known variants with established mutational signatures but is limited in detecting novel variants. Clinical laboratories of all resource levels use qPCR-based variant screening due to ease of implementation and cost-efficiency. The widespread uptake of this technology is directly attributable to early work in the area that demonstrated the feasibility of variant identification using a combination of mutational markers, including all or some of the following: L452R, E484K, N501Y, H69_V70del. [5–7, 16, 17] Nevertheless, new variants with unique mutational signatures may not be identifiable with these existing qPCR-based screenings, demonstrating the need for further monitoring of low-level or novel variants. Sequencing is the gold standard for identifying variants and the cornerstone for global monitoring of SARS-CoV-2 viral evolution. In fact, qPCR-based monitoring methods rely on global sequencing efforts to stay abreast with the newest viral strains. Resource-rich facilities can use sequencing for variant surveillance on a representative subset of unbiased samples. However, random sampling for sequencing may not be statistically geared to identify low-level variants on a regional basis. For example, identifying the Delta+E484K variant using unbiased sample selection for sequencing would necessitate sequencing thousands of un-mutated Delta variant samples. Turnaround time can be slower as sequencing runs may be delayed due to demands on laboratory resources when caseloads are high or until batches are filled when caseloads are low, limiting the utility of variant monitoring for real-time outbreak management. To enrich the sequencing data for information, some use targeted sequencing of mutational hotspots or restrict the population to hospitalized patients. [15, 18] Nevertheless, monitoring variants by exclusively sequencing is not feasible for most clinical laboratories due to resource limitations, preventing more widespread variant monitoring. Overall, both qPCR screening and unbiased population sequencing have some capacity to detect novel low-level variants on a regional basis. Using a tiered approach, we amplify the strengths of each, which may accelerate the identification of emerging clinically important strains. We implemented a pipeline with a qPCR-based screen for the most common variants and subsequent sequencing of select samples. Examples of sequenced samples include: ambiguous qPCR expression profiles, clinically interesting cases, or quality control samples. Depending on available resources, sequencing can be performed in-house or commercially. In short, this tiered approach matches the strengths of each technology to the most informative samples. It is not intended to replace sequencing nor qPCR-based screening at all hospitals; it is an alternative approach for facilities to optimize their existing resources. One advantage of the tiered approach is its capacity for dynamic feedback and iterative improvement. As the virus evolves, so too must the tests we use to detect it. In our workflow, the variant calling results can be used to adjust the primers and the screening ruleset to reflect variability in the regional level of variants. For example, if a variant of concern is noted in another country, such as Delta, we incorporate primers to detect it before the variant arrives in our population. Or, if a novel mutation is seen in the quality control samples from a dominant strain, the ruleset can be adjusted to account for the new subtype, such as for the Delta with an E484K. The current setup will not catch all mutations or new variants; instead, it is designed to be adaptable and broadly implementable. Given that the SARS-CoV2 genome is accumulating mutations and new variants are evolving, stagnant screening systems can become obsolete.[19–21] A tiered approach allows for continuous quality control and improvement to match the pace of variant evolution. The merits of this system are exemplified by our timely detection of rare variants while maintaining high-throughput variant identification of the dominant strains. Due to the high-throughput nature of the qPCR screen, we identified the emerging Delta variant and, within the first few weeks, collected sufficient samples to perform preliminary characterization of it in our population. We found the Delta variant is associated with younger patient age, lower CT values, and increased vaccine breakthrough. As a testament to the utility of sequencing, we were able to detect a rare and emerging variant of Delta with an E484K mutation, which had not yet been described in our population and only a handful of times previously in the USA. Therefore, our approach demonstrates the strengths of both qPCR-based screening and sequencing and outlines the practical implementation of these tests using existing and strained clinical laboratory infrastructure during the middle of an evolving pandemic. ## Supporting information Supplemental [[supplements/262625_file06.pdf]](pending:yes) ## Data Availability Data sharing: All data relevant to the study are included in the article or uploaded as supplementary information. Sequencing data have been submitted to Genebank. ## Acknowledgements We are grateful for the technical, medical, and administrative staff at University Hospitals Health System for their unwavering support and resilience of our SAR2- CoV-2 testing and variant identification programs. ## Footnotes * **Competing interests:** No competing interests to disclose. * **Funding:** None * **Data sharing:** All data relevant to the study are included in the article or uploaded as supplementary information. Sequencing data have been submitted to Genebank. * **Patient consent and ethical approval:** The study was reviewed by the University Hospitals Institutional Review Board (Cleveland, Ohio) and ethical approval was waived. * **Reporting guidelines:** The Equator network did not recommend a specific set of guidelines for pathology methods papers. * Discussion section updated. * Received August 31, 2021. * Revision received September 7, 2021. * Accepted September 8, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.CDC. SARS-CoV-2 Variant Classifications and Definitions. [website] 2021 Aug. 17, 2021 [cited 2021 08/18/2021]; Available from: [https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-info.html](https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-info.html). 2. 2.Snell, L.B., et al., Rapid genome sequencing in hospitals to identify potential vaccine-escape SARS-CoV-2 variants. Lancet Infect Dis, 2021. 3. 3.Xia, F., et al., Quantifying competitive advantages of mutant strains in a population involving importation and mass vaccination rollout. Infect Dis Model, 2021. 6: p. 988–996. 4. 4.Alouani, D.J., et al., Specificity of SARS-CoV-2 Real-Time PCR Improved by Deep Learning Analysis. J Clin Microbiol, 2021. 59(6). 5. 5.Korukluoglu, G., et al., 40 minutes RT-qPCR Assay for Screening Spike N501Y and HV69-70del Mutations. bioRxiv, 2021. 6. 6.Vogels, C.B.F., et al., Multiplex qPCR discriminates variants of concern to enhance global surveillance of SARS-CoV-2. PLoS Biol, 2021. 19(5): p. e3001236. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pbio.3001236&link_type=DOI) 7. 7.Wang, H., et al., Multiplex SARS-CoV-2 Genotyping Reverse Transcriptase PCR for Population-Level Variant Screening and Epidemiologic Surveillance. J Clin Microbiol, 2021. 59(8): p. e0085921. 8. 8.Sekulic, M., et al., Molecular Detection of SARS-CoV-2 Infection in FFPE Samples and Histopathologic Findings in Fatal SARS-CoV-2 Cases. Am J Clin Pathol, 2020. 154(2): p. 190–200. 9. 9.Mishra, S., et al., Changing composition of SARS-CoV-2 lineages and rise of Delta variant in England. EClinicalMedicine, 2021. 39: p. 101064. 10. 10.Musser, J.M., et al., Delta variants of SARS-CoV-2 cause significantly increased vaccine breakthrough COVID-19 cases in Houston, Texas. medRxiv, 2021. 11. 11.Kirola, L., Genetic emergence of B.1.617.2 in COVID-19. New Microbes New Infect, 2021. 43: p. 100929. 12. 12.Lopez Bernal, J., et al., Effectiveness of Covid-19 Vaccines against the B.1.617.2 (Delta) Variant. N Engl J Med, 2021. 385(7): p. 585–594. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F08%2F2021.08.31.21262625.atom) 13. 13.Li, B., et al., Viral infection and transmission in a large, well-traced outbreak caused by the SARS-CoV-2 Delta variant. medRxiv, 2021. 14. 14.Alaimo, J.T., et al., Loss-of-function mutations in ISCA2 disrupt 4Fe-4S cluster machinery and cause a fatal leukodystrophy with hyperglycinemia and mtDNA depletion. Hum Mutat, 2018. 39(4): p. 537–549. 15. 15.Page, A.J., et al., Large-scale sequencing of SARS-CoV-2 genomes from one region allows detailed epidemiology and enables local outbreak management. Microb Genom, 2021. 7(6). 16. 16.Zhao, Y., et al., A novel diagnostic test to screen SARS-CoV-2 variants containing E484K and N501Y mutations. Emerg Microbes Infect, 2021. 10(1): p. 994–997. 17. 17.Wang, H., et al., Mutation-Specific SARS-CoV-2 PCR Screen: Rapid and Accurate Detection of Variants of Concern and the Identification of a Newly Emerging Variant with Spike L452R Mutation. J Clin Microbiol, 2021. 59(8): p. e0092621. 18. 18.Jørgensen, T.S., et al., A rapid, cost efficient and simple method to identify current SARS-CoV-2 variants of concern by Sanger sequencing part of the spike protein gene. medRxiv, 2021. 19. 19.Bandoy, D. and B.C. Weimer, Analysis of SARS-CoV-2 genomic epidemiology reveals disease transmission coupled to variant emergence and allelic variation. Sci Rep, 2021. 11(1): p. 7380. 20. 20.Jonsson, H., et al., Molecular benchmarks of a SARS-CoV-2 epidemic. Nat Commun, 2021. 12(1): p. 3633. 21. 21.Sjaarda, C.P., et al., Temporal Dynamics and Evolution of SARS-CoV-2 Demonstrate the Necessity of Ongoing Viral Genome Sequencing in Ontario, Canada. mSphere, 2021. 6(3).