A Comparison of Whole Genome Sequencing of SARS-CoV-2 Using Amplicon-Based Sequencing, Random Hexamers, and Bait Capture

Pathogenesis and transmission of SARS-CoV-2 D614G, Alpha, Gamma, Delta, and Omicron variants in golden hamsters

Since the emergence of SARS-CoV-2 in humans, novel variants have evolved to become dominant circulating lineages. These include D614G (B.1 lineage), Alpha (B.1.1.7), Gamma (P.1), Delta (B.1.617.2), and Omicron BA.1 (B.1.1.529) and BA.2 (B.1.1.529.2) viruses. Here, we compared the viral replication, pathogenesis, and transmissibility of these variants. Replication kinetics and innate immune response against the viruses were tested in ex vivo human nasal epithelial cells (HNEC) and induced pluripotent stem cell-derived lung organoids (IPSC-LOs), and the golden hamster model was employed to test pathogenicity and potential for transmission by the respiratory route. Delta, BA.1, and BA.2 viruses replicated more efficiently, and outcompeted D614G, Alpha, and Gamma viruses in an HNEC competition assay. BA.1 and BA.2 viruses, however, replicated poorly in IPSC-LOs compared to other variants. Moreover, BA.2 virus infection significantly increased secretion of IFN-λ1, IFN-λ2, IFN-λ3, IL-6, and IL-1RA in HNECs relative to D614G infection, but not in IPSC-LOs. The BA.1 and BA.2 viruses replicated less effectively in hamster lungs compared to the other variants; and while the Gamma virus reached titers comparable to D614G and Delta viruses, it caused greater lung pathology. Lastly, the Gamma and Delta variants transmitted more efficiently by the respiratory route compared to the other viruses, while BA.1 and BA.2 viruses transmitted less efficiently. These findings demonstrate the ongoing utility of experimental risk assessment as SARS-CoV-2 variants continue to evolve.

Identification of patient demographic, clinical, and SARS-CoV-2 genomic factors associated with severe COVID-19 using supervised machine learning: a retrospective multicenter study

BACKGROUND

Drivers of COVID-19 severity are multifactorial and include multidimensional and potentially interacting factors encompassing viral determinants and host-related factors (i.e., demographics, pre-existing conditions and/or genetics), thus complicating the prediction of clinical outcomes for different severe acute respiratory syndrome coronavirus (SARS-CoV-2) variants. Although millions of SARS-CoV-2 genomes have been publicly shared in global databases, linkages with detailed clinical data are scarce. Therefore, we aimed to establish a COVID-19 patient dataset with linked clinical and viral genomic data to then examine associations between SARS-CoV-2 genomic signatures and clinical disease phenotypes.

METHODS

A cohort of adult patients with laboratory confirmed SARS-CoV-2 from 11 participating healthcare institutions in the Greater Toronto Area (GTA) were recruited from March 2020 to April 2022. Supervised machine learning (ML) models were developed to predict hospitalization using SARS-CoV-2 lineage-specific genomic signatures, patient demographics, symptoms, and pre-existing comorbidities. The relative importance of these features was then evaluated.

RESULTS

Complete clinical data and viral whole genome level information were obtained from 617 patients, 50.4% of whom were hospitalized. Notably, inpatients were older with a mean age of 66.67 years (SD ± 17.64 years), whereas outpatients had a mean age of 44.89 years (SD ± 16.00 years). SHapley Additive exPlanations (SHAP) analyses revealed that underlying vascular disease, underlying pulmonary disease, and fever were the most significant clinical features associated with hospitalization. In models built on the amino acid sequences of functional regions including spike, nucleocapsid, ORF3a, and ORF8 proteins, variants preceding the emergence of variants of concern (VOCs) or pre-VOC variants, were associated with hospitalization.

CONCLUSIONS

Viral genomic features have limited utility in predicting hospitalization across SARS-CoV-2 diversity. Combining clinical and viral genomic datasets provides perspective on patient specific and virus-related factors that impact COVID-19 disease severity. Overall, clinical features had greater discriminatory power than viral genomic features in predicting hospitalization.

Enriching the future of public health microbiology with hybridization bait capture

SUMMARYPublic health microbiology focuses on microorganisms and infectious agents that impact human health. For years, this field has relied on culture or molecular methods to investigate complex samples of public health importance. However, with the increase in accuracy and decrease in sequencing cost over the last decade, there has been a transition to the use of next-generation sequencing in public health microbiology. Nevertheless, many available sequencing methods (e.g., shotgun metagenomics and amplicon sequencing) do not work well in complex sample types, require deep sequencing, or have inherent biases associated with them. Hybridization bait capture, also known as target enrichment, brings in solutions for such limitations. It is an increasingly popular technique to simultaneously characterize many thousands of genetic elements while reducing the amount of sequencing needed (thereby reducing the sequencing costs). Here, we summarize the concept of hybridization bait capture for public health, reviewing a total of 35 bait sets designed in six key topic areas for public health microbiology [i.e., antimicrobial resistance (AMR), bacteria, fungi, parasites, vectors, and viruses], and compare hybridization bait capture to previously relied upon methods. Furthermore, we provide an in-depth comparison of the three most popular bait sets designed for AMR by evaluating each of them against three major AMR databases: Comprehensive Antibiotic Resistance Database, Microbial Ecology Group Antimicrobial Resistance Database, and Pathogenicity Island Database. Thus, this article provides a review of hybridization bait capture for public health microbiologists.

Development of targeted whole genome sequencing approaches for Crimean-Congo haemorrhagic fever virus (CCHFV)

Crimean-Congo haemorrhagic fever (CCHF) is the most prevalent human tick-borne viral disease, with a reported case fatality rate of 30 % or higher. The virus contains a tri-segmented, negative-sense RNA genome consisting of the small (S), medium (M) and large (L) segments encoding respectively the nucleoprotein (NP), the glycoproteins precursor (GPC) and the viral RNA-dependent RNA polymerase (RDRP). CCHFV is one of the most genetically diverse arboviruses, with seven distinct lineages named after the region they were first reported in and based on S segment phylogenetic analysis. Due to the high genetic divergence of the virus, a single targeted tiling PCR strategy to enrich for viral nucleic acids prior to sequencing is difficult to develop, and previously we have developed and validated a tiling PCR enrichment method for the Europe 1 genetic lineage. We have developed a targeted, probe hybridisation capture method and validated its performance on clinical as well as cell-cultured material of CCHFV from different genetic lineages, including Europe 1, Europe 2, Africa 2 and Africa 3. The method produced over 95 % reference coverages with at least 10x sequencing depth. While we were only able to recover a single complete genome sequence from the tested Europe 1 clinical samples with the capture hybridisation protocol, the data provides evidence of its applicability to different CCHFV genetic lineages. CCHFV is an important tick-borne human pathogen with wide geographical distribution. Environmental as well as anthropogenic factors are causing increased CCHFV transmission. Development of strategies to recover CCHFV sequences from genetically diverse lineages of the virus is of paramount importance to monitor the presence of the virus in new areas, and in public health responses for CCHFV molecular surveillance to rapidly detect, diagnose and characterise currently circulating strains.

UnCoVar: a reproducible and scalable workflow for transparent and robust virus variant calling and lineage assignment using SARS-CoV-2 as an example

BACKGROUND

At a global scale, the SARS-CoV-2 virus did not remain in its initial genotype for a long period of time, with the first global reports of variants of concern (VOCs) in late 2020. Subsequently, genome sequencing has become an indispensable tool for characterizing the ongoing pandemic, particularly for typing SARS-CoV-2 samples obtained from patients or environmental surveillance. For such SARS-CoV-2 typing, various in vitro and in silico workflows exist, yet to date, no systematic cross-platform validation has been reported.

RESULTS

In this work, we present the first comprehensive cross-platform evaluation and validation of in silico SARS-CoV-2 typing workflows. The evaluation relies on a dataset of 54 patient-derived samples sequenced with several different in vitro approaches on all relevant state-of-the-art sequencing platforms. Moreover, we present UnCoVar, a robust, production-grade reproducible SARS-CoV-2 typing workflow that outperforms all other tested approaches in terms of precision and recall.

CONCLUSIONS

In many ways, the SARS-CoV-2 pandemic has accelerated the development of techniques and analytical approaches. We believe that this can serve as a blueprint for dealing with future pandemics. Accordingly, UnCoVar is easily generalizable towards other viral pathogens and future pandemics. The fully automated workflow assembles virus genomes from patient samples, identifies existing lineages, and provides high-resolution insights into individual mutations. UnCoVar includes extensive quality control and automatically generates interactive visual reports. UnCoVar is implemented as a Snakemake workflow. The open-source code is available under a BSD 2-clause license at github.com/IKIM-Essen/uncovar.

A novel fast hybrid capture sequencing method for high-efficiency common human coronavirus whole-genome acquisition

Rapid and accurate sequencing of the entire viral genome, coupled with continuous monitoring of genetic changes, is crucial for understanding the epidemiology of coronaviruses. We designed a novel method called micro target hybrid capture system (MT-Capture) to enable whole-genome sequencing in a timely manner. The novel design of probes used in target binding exhibits a unique and synergistic "hand-in-hand" conjugation effect. The entire hybrid capture process is within 2.5 hours, overcoming the time-consuming and complex operation characteristics of the traditional liquid-phase hybrid capture (T-Capture) system. By designing specific probes for these coronaviruses, MT-Capture effectively enriched isolated strains and 112 clinical samples of coronaviruses with cycle threshold values below 37. Compared to multiplex PCR sequencing, it does not require frequent primer updates and has higher compatibility. MT-Capture is highly sensitive and capable of tracking variants.IMPORTANCEMT-Capture is meticulously designed to enable the efficient acquisition of the target genome of the common human coronavirus. Coronavirus is a kind of virus that people are generally susceptible to and is epidemic and infectious, and it is the virus with the longest genome among known RNA viruses. Therefore, common human coronavirus samples are selected to evaluate the accuracy and sensitivity of MT-Capture. This method utilizes innovative probe designs optimized through probe conjugation techniques, greatly shortening the time and simplifying the handwork compared with traditional hybridization capture processes. Our results demonstrate that MT-Capture surpasses multiplex PCR in terms of sensitivity, exhibiting a thousandfold increase. Moreover, MT-Capture excels in the identification of mutation sites. This method not only is used to target the coronaviruses but also may be used to diagnose other diseases, including various infectious diseases, genetic diseases, or tumors.

Performance of amplicon and capture based next-generation sequencing approaches for the epidemiological surveillance of Omicron SARS-CoV-2 and other variants of concern

To control the SARS-CoV-2 pandemic, healthcare systems have focused on ramping up their capacity for epidemiological surveillance through viral whole genome sequencing. In this paper, we tested the performance of two protocols of SARS-CoV-2 nucleic acid enrichment, an amplicon enrichment using different versions of the ARTIC primer panel and a hybrid-capture method using KAPA RNA Hypercap. We focused on the challenge of the Omicron variant sequencing, the advantages of automated library preparation and the influence of the bioinformatic analysis in the final consensus sequence. All 94 samples were sequenced using Illumina iSeq 100 and analysed with two bioinformatic pipelines: a custom-made pipeline and an Illumina-owned pipeline. We were unsuccessful in sequencing six samples using the capture enrichment due to low reads. On the other hand, amplicon dropout and mispriming caused the loss of mutation G21987A and the erroneous addition of mutation T15521A respectively using amplicon enrichment. Overall, we found high sequence agreement regardless of method of enrichment, bioinformatic pipeline or the use of automation for library preparation in eight different SARS-CoV-2 variants. Automation and the use of a simple app for bioinformatic analysis can simplify the genotyping process, making it available for more diagnostic facilities and increasing global vigilance.

Lessons learned: overcoming common challenges in reconstructing the SARS-CoV-2 genome from short-read sequencing data via CoVpipe2

Background

Accurate genome sequences form the basis for genomic surveillance programs, the added value of which was impressively demonstrated during the COVID-19 pandemic by tracing transmission chains, discovering new viral lineages and mutations, and assessing them for infectiousness and resistance to available treatments. Amplicon strategies employing Illumina sequencing have become widely established for variant detection and reference-based reconstruction of SARS-CoV-2 genomes, and are routine bioinformatics tasks. Yet, specific challenges arise when analyzing amplicon data, for example, when crucial and even lineage-determining mutations occur near primer sites.

Methods

We present CoVpipe2, a bioinformatics workflow developed at the Public Health Institute of Germany to reconstruct SARS-CoV-2 genomes based on short-read sequencing data accurately. The decisive factor here is the reliable, accurate, and rapid reconstruction of genomes, considering the specifics of the used sequencing protocol. Besides fundamental tasks like quality control, mapping, variant calling, and consensus generation, we also implemented additional features to ease the detection of mixed samples and recombinants.

Results

We highlight common pitfalls in primer clipping, detecting heterozygote variants, and dealing with low-coverage regions and deletions. We introduce CoVpipe2 to address the above challenges and have compared and successfully validated the pipeline against selected publicly available benchmark datasets. CoVpipe2 features high usability, reproducibility, and a modular design that specifically addresses the characteristics of short-read amplicon protocols but can also be used for whole-genome short-read sequencing data.

Conclusions

CoVpipe2 has seen multiple improvement cycles and is continuously maintained alongside frequently updated primer schemes and new developments in the scientific community. Our pipeline is easy to set up and use and can serve as a blueprint for other pathogens in the future due to its flexibility and modularity, providing a long-term perspective for continuous support. CoVpipe2 is written in Nextflow and is freely accessible from \href{https://github.com/rki-mf1/CoVpipe2}{github.com/rki-mf1/CoVpipe2} under the GPL3 license.

Enrichment of SARS-CoV-2 sequence from nasopharyngeal swabs whilst identifying the nasal microbiome

Simultaneously characterising the genomic information of coronaviruses and the underlying nasal microbiome from a single clinical sample would help characterise infection and disease. Metatranscriptomic approaches can be used to sequence SARS-CoV-2 (and other coronaviruses) and identify mRNAs associated with active transcription in the nasal microbiome. However, given the large sequence background, unenriched metatranscriptomic approaches often do not sequence SARS-CoV-2 to sufficient read and coverage depth to obtain a consensus genome, especially with moderate and low viral loads from clinical samples. In this study, various enrichment methods were assessed to detect SARS-CoV-2, identify lineages and define the nasal microbiome. The methods were underpinned by Oxford Nanopore long-read sequencing and variations of sequence independent single primer amplification (SISPA). The utility of the method(s) was also validated on samples from patients infected seasonal coronaviruses. The feasibility of profiling the nasal microbiome using these enrichment methods was explored. The findings shed light on the performance of different enrichment strategies and their applicability in characterising the composition of the nasal microbiome.

fosA11, a novel chromosomal-encoded fosfomycin resistance gene identified in Providencia rettgeri

This study investigated resistance genes corresponding to the fosfomycin resistance phenotype in clinical isolate Providencia rettgeri W986, as well as characterizing the enzymatic activity of FosA11 and the genetic environment. Antimicrobial susceptibility testing was performed using the agar microdilution method based on the Clinical and Laboratory Standards Institute guidelines. The whole genomic sequence of Providencia rettgeri W986 was obtained using Illumina sequencing and the PacBio platform. The fosA-11 gene was amplified by PCR and cloned into the pUCP20 vector. The recombinant strain pCold1-fosA11-BL21 was expressed to extract the target protein, and absorbance photometry was applied for enzymatic parameter determination. Minimal inhibitory concentration (MIC) tests showed that W986 conferred fosfomycin resistance and was inhibited by phosphonoformate, thereby indicating the presence of a FosA protein. A novel resistance gene designated as fosA11 was identified by whole-genome sequencing and bioinformatics analysis, and it shared 54.41%-64.23% amino acid identity with known FosA proteins. Cloning fosA11 into Escherichia coli obtained a significant increase (32-fold) in the MIC with fosfomycin. Determination of the enzyme kinetics showed that FosA11 had a high catalytic effect on fosfomycin, with Km = 18 ± 4 and Kcat = 56.1 ± 3.2. We also found that fosA11 was located on the chromosome, but the difference in the GC content between the chromosome and fosA11 was dubious, and thus further investigation is required. In this study, we identified and characterized a novel fosfomycin inactivation enzyme called FosA11. The origin and prevalence of the fosA11 gene in other bacteria require further investigation.IMPORTANCEFosfomycin is an effective antimicrobial agent against Enterobacterales strains. However, the resistance rate of fosfomycin is increasing year by year. Therefore, it is necessary to study the deep molecular mechanism of bacterial resistance to fosfomycin. We identified a novel chromosomal fosfomycin glutathione S-transferase, FosA11 from Providencia rettgeri, which shares a very low identity (54.41%-64.23%) with the previously known FosA and exhibits highly efficient catalytic ability against fosfomycin. Analysis of the genetic context and origin of fosA11 displays that the gene and its surrounding environments are widely conserved in Providencia and no mobile elements are discovered, implying that FosA11 may be broadly important in the natural resistance to fosfomycin of Providencia species.

Nosocomial surveillance of multidrug-resistant Acinetobacter baumannii: a genomic epidemiological study

Acinetobacter baumannii is a major opportunistic pathogen causing hospital-acquired infections, and it is imperative to comprehend its evolutionary and epidemiological dynamics in hospitals to prevent and control nosocomial transmission. Here, we present a comprehensive genomic epidemiological study involving the genomic sequencing and antibiotic resistance profiling of 634 A. baumannii strains isolated from seven intensive care units (ICUs) of a Chinese general hospital over 2 consecutive years. Our study reveals that ST2 is highly dominant (90.54%) in the ICUs, with 98.90% of the ST2 exhibiting multidrug resistant or extensively drug resistant. Phylogenetic analyses of newly sequenced genomes and public data suggest that nosocomial isolates originated outside the hospital but evolved inside. The major lineages appear to be stable, with 9 of the 28 identified nosocomial epidemic clones infecting over 60% of the affected patients. However, outbreaks of two highly evolved clones have been observed in different hospitals, suggesting significant inter-hospital transmission chains. By coupling patient medical records and genomic divergence of the ST2, we found that cross-ward patient transfer played a crucial role in pathogen's nosocomial transmission. Additionally, we identified 831 potential adaptive evolutionary loci and 44 associated genes by grouping and comparing the genomes of clones with different prevalence. Overall, our study provides a comprehensive and contemporary survey on the epidemiology and genomic evolution of A. baumannii in a large Chinese general hospital. These findings shed light on the nosocomial evolution and transmission of A. baumannii and offers valuable information for transmission prevention and antibiotic therapy.IMPORTANCEThis study delved into the genomic evolution and transmission of nosocomial Acinetobacter baumannii on a large scale, spanning both an extended time period and the largest sample size to date. Through molecular epidemiological investigations based on genomics, we can directly trace the origin of the pathogen, detecting and monitoring outbreaks of infectious diseases in a timely manner, and ensuring public health safety. In addition, this study also collects a large amount of genomic and antibiotic resistance detection data, which is helpful for phenotype prediction based on genomic sequencing. It enables patients to receive personalized antibiotic treatment quickly, helps doctors select antibiotics more accurately, and contributes to reducing the use of antibiotics and lowering the risk of antibiotic resistance development.

Applying targeted gene hybridization capture to viruses with a focus to SARS-CoV-2

Although next-generation sequencing technologies are advancing rapidly, many research topics often require selective sequencing of genomic regions of interest. In addition, sequencing low-titre viruses is challenging, especially for coronaviruses, which are the largest RNA viruses. Prior to sequencing, enrichment of viral particles can help to significantly increase target sequence information as well as avoid large sequencing efforts and, consequently, can increase sensitivity and reduce sequencing costs. Targeting nucleic acids using capture by hybridization is another efficient method that can be performed by applying complementary probes (DNA or RNA baits) to directly enrich genetic information of interest while removing background non-target material. In studies where sequence capture by hybridization has been applied to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus, most authors agree that this technique is useful to easily access sequence targets in complex samples. Furthermore, this approach allows for complete or near-complete sequencing of the viral genome, even in samples with low viral load or poor nucleic acid integrity. In addition, this strategy is highly efficient at discovering new variants by facilitating downstream investigations, such as phylogenetics, epidemiology, and evolution. Commercial kits, as well as in-house protocols, have been developed for enrichment of viral sequences. However, these kits have multiple variations in procedure, with differences in performance. This review compiles and describes studies in which hybridization capture has been applied to SARS-CoV-2 variant genomes.

Synergistic horizontal transfer of antibiotic resistance genes and transposons in the infant gut microbial genome

Transposons, plasmids, bacteriophages, and other mobile genetic elements facilitate horizontal gene transfer in the gut microbiota, allowing some pathogenic bacteria to acquire antibiotic resistance genes (ARGs). Currently, the relationship between specific ARGs and specific transposons in the comprehensive infant gut microbiome has not been elucidated. In this study, ARGs and transposons were annotated from the Unified Human Gastrointestinal Genome (UHGG) and the Early-Life Gut Genomes (ELGG). Association rules mining was used to explore the association between specific ARGs and specific transposons in UHGG, and the robustness of the association rules was validated using the external database in ELGG. Our results suggested that ARGs and transposons were more likely to be relevant in infant gut microbiota compared to adult gut microbiota, and nine robust association rules were identified, among which Klebsiella pneumoniae, Enterobacter hormaechei_A, and Escherichia coli_D played important roles in this association phenomenon. The emphasis of this study is to investigate the synergistic transfer of specific ARGs and specific transposons in the infant gut microbiota, which can contribute to the study of microbial pathogenesis and the ARG dissemination dynamics.IMPORTANCEThe transfer of transposons carrying antibiotic resistance genes (ARGs) among microorganisms accelerates antibiotic resistance dissemination among infant gut microbiota. Nonetheless, it is unclear what the relationship between specific ARGs and specific transposons within the infant gut microbiota. K. pneumoniae, E. hormaechei_A, and E. coli_D were identified as key players in the nine robust association rules we discovered. Meanwhile, we found that infant gut microorganisms were more susceptible to horizontal gene transfer events about specific ARGs and specific transposons than adult gut microorganisms. These discoveries could enhance the understanding of microbial pathogenesis and the ARG dissemination dynamics within the infant gut microbiota.

Chronic COVID-19 infection in an immunosuppressed patient shows changes in lineage over time: a case report

BACKGROUND

The COVID-19 pandemic, caused by the Severe Acute Respiratory Syndrome Coronavirus 2 virus, emerged in late 2019 and spready globally. Many effects of infection with this pathogen are still unknown, with both chronic and repeated COVID-19 infection producing novel pathologies.

CASE PRESENTATION

An immunocompromised patient presented with chronic COVID-19 infection. The patient had history of Hodgkin's lymphoma, treated with chemotherapy and stem cell transplant. During the course of their treatment, eleven respiratory samples from the patient were analyzed by whole-genome sequencing followed by lineage identification. Whole-genome sequencing of the virus present in the patient over time revealed that the patient at various timepoints harboured three different lineages of the virus. The patient was initially infected with the B.1.1.176 lineage before coinfection with BA.1. When the patient was coinfected with both B.1.1.176 and BA.1, the viral populations were found in approximately equal proportions within the patient based on sequencing read abundance. Upon further sampling, the lineage present within the patient during the final two timepoints was found to be BA.2.9. The patient eventually developed respiratory failure and died.

CONCLUSIONS

This case study shows an example of the changes that can happen within an immunocompromised patient who is infected with COVID-19 multiple times. Furthermore, this case demonstrates how simultaneous coinfection with two lineages of COVID-19 can lead to unclear lineage assignment by standard methods, which are resolved by further investigation. When analyzing chronic COVID-19 infection and reinfection cases, care must be taken to properly identify the lineages of the virus present.

A standardized nomenclature for resistance-modifying agents in the Comprehensive Antibiotic Resistance Database

IMPORTANCE

While increasing rates of antimicrobial resistance undermine our current arsenal of antibiotics, resistance-modifying agents (RMAs) hold promise to extend the lifetime of these important molecules. We here provide a standardized nomenclature for RMAs within the Comprehensive Antibiotic Resistance Database in aid of RMA discovery, data curation, and genome mining.

Direct prediction of carbapenem resistance in Pseudomonas aeruginosa by whole genome sequencing and metagenomic sequencing

Carbapenem resistance is a major concern in the management of antibiotic-resistant Pseudomonas aeruginosa infections. The direct prediction of carbapenem-resistant phenotype from genotype in P. aeruginosa isolates and clinical samples would promote timely antibiotic therapy. The complex carbapenem resistance mechanism and the high prevalence of variant-driven carbapenem resistance in P. aeruginosa make it challenging to predict the carbapenem-resistant phenotype through the genotype. In this study, using whole genome sequencing (WGS) data of 1,622 P. aeruginosa isolates followed by machine learning, we screened 16 and 31 key gene features associated with imipenem (IPM) and meropenem (MEM) resistance in P. aeruginosa, including oprD(HIGH), and constructed the resistance prediction models. The areas under the curves of the IPM and MEM resistance prediction models were 0.906 and 0.925, respectively. For the direct prediction of carbapenem resistance in P. aeruginosa from clinical samples by the key gene features selected and prediction models constructed, 72 P. aeruginosa-positive sputum samples were collected and sequenced by metagenomic sequencing (MGS) based on next-generation sequencing (NGS) or Oxford Nanopore Technology (ONT). The prediction applicability of MGS based on NGS outperformed that of MGS based on ONT. In 72 P. aeruginosa-positive sputum samples, 65.0% (26/40) of IPM-insensitive and 63.2% (24/38) of MEM-insensitive P. aeruginosa were directly predicted by NGS-based MGS with positive predictive values of 0.897 and 0.889, respectively. By the direct detection of the key gene features associated with carbapenem resistance of P. aeruginosa, the carbapenem resistance of P. aeruginosa could be directly predicted from cultured isolates by WGS or from clinical samples by NGS-based MGS, which could assist the timely treatment and surveillance of carbapenem-resistant P. aeruginosa.

Genomic and transcriptomic characterization of delta SARS-CoV-2 infection in free-ranging white-tailed deer (Odocoileus virginianus)

White-tailed deer (WTD) are susceptible to SARS-CoV-2 and represent an important species for surveillance. Samples from WTD (n = 258) collected in November 2021 from Québec, Canada were analyzed for SARS-CoV-2 RNA. We employed viral genomics and host transcriptomics to further characterize infection and investigate host response. We detected Delta SARS-CoV-2 (B.1.617.2) in WTD from the Estrie region; sequences clustered with human sequences from October 2021 from Vermont, USA, which borders this region. Mutations in the S-gene and a deletion in ORF8 were detected. Host expression patterns in SARS-CoV-2 infected WTD were associated with the innate immune response, including signaling pathways related to anti-viral, pro- and anti-inflammatory signaling, and host damage. We found limited correlation between genes associated with innate immune response from human and WTD nasal samples, suggesting differences in responses to SARS-CoV-2 infection. Our findings provide preliminary insights into host response to SARS-CoV-2 infection in naturally infected WTD.

Characterisation of neonatal Staphylococcus capitis NRCS-A isolates compared with non NRCS-A Staphylococcus capitis from neonates and adults

Staphylococcus capitis is a frequent cause of late-onset sepsis in neonates admitted to Neonatal Intensive Care Units (NICU). One clone of S. capitis, NRCS-A has been isolated from NICUs globally although the reasons for the global success of this clone are not well understood.We analysed a collection of S. capitis colonising babies admitted to two NICUs, one in the UK and one in Germany as well as corresponding pathological clinical isolates. Genome analysis identified a population structure of three groups; non-NRCS-A isolates, NRCS-A isolates, and a group of 'proto NRCS-A' - isolates closely related to NRCS-A but not associated with neonatal infection. All bloodstream isolates belonged to the NRCS-A group and were indistinguishable from strains carried on the skin or in the gut. NRCS-A isolates showed increased tolerance to chlorhexidine and antibiotics relative to the other S. capitis as well as enhanced ability to grow at higher pH values. Analysis of the pangenome of 138 isolates identified characteristic nsr and tarJ genes in both the NRCS-A and proto groups. A CRISPR-cas system was only seen in NRCS-A isolates which also showed enrichment of genes for metal acquisition and transport.We found evidence for transmission of S. capitis NRCS-A within NICU, with related isolates shared between babies and multiple acquisitions by some babies. Our data show NRCS-A strains commonly colonise uninfected babies in NICU representing a potential reservoir for potential infection. This work provides more evidence that adaptation to survive in the gut and on skin facilitates spread of NRCS-A, and that metal acquisition and tolerance may be important to the biology of NRCS-A. Understanding how NRCS-A survives in NICUs can help develop infection control procedures against this clone.

Use of whole genome sequencing to identify low-frequency mutations in SARS-CoV-2 patients treated with remdesivir

BACKGROUND

Remdesivir (RDV) has been shown to reduce hospitalization and mortality in COVID-19 patients. Resistance mutations caused by RDV are rare and have been predominantly reported in patients who are on prolonged therapy and immunocompromised. We investigate the effects of RDV treatment on intra-host SARS-CoV-2 diversity and low-frequency mutations in moderately ill hospitalized COVID-19 patients and compare them to patients without RDV treatment.

METHODS

From March 2020 to April 2022, sequential collections of nasopharyngeal and mid-turbinate swabs were obtained from 14 patients with and 30 patients without RDV treatment. Demographic and clinical data on all patients were reviewed. A total of 109 samples were sequenced and mutation analyses were performed.

RESULTS

Previously reported drug resistant mutations in nsp12 were not identified during short courses of RDV therapy. In genes encoding and surrounding the replication complex (nsp6-nsp14), low-frequency minority variants were detected in 7/14 (50%) and 18/30 (60%) patients with and without RDV treatment, respectively. We did not detect significant differences in within-host diversity and positive selection between the RDV-treated and untreated groups.

CONCLUSIONS

Minimal intra-host variability and stochastic low-frequency variants detected in moderately ill patients suggests little selective pressure in patients receiving short courses of RDV. The barrier to RDV resistance is high in patients with moderate disease. Patients undergoing short regimens of RDV therapy should continue to be monitored.

SARS-CoV-2 Outbreak Investigation Using Contact Tracing and Whole-Genome Sequencing in an Ontario Tertiary Care Hospital

Genomic epidemiology can facilitate an understanding of evolutionary history and transmission dynamics of a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreak. We used next-generation sequencing techniques to study SARS-CoV-2 genomes isolated from patients and health care workers (HCWs) across five wards of a Canadian hospital with an ongoing SARS-CoV-2 outbreak. Using traditional contact tracing methods, we show transmission events between patients and HCWs, which were also supported by the SARS-CoV-2 lineage assignments. The outbreak predominantly involved SARS-CoV-2 B.1.564.1 across all five wards, but we also show evidence of community introductions of lineages B.1, B.1.1.32, and B.1.231, falsely assumed to be outbreak related. Altogether, our study exemplifies the value of using contact tracing in combination with genomic epidemiology to understand the transmission dynamics and genetic underpinnings of a SARS-CoV-2 outbreak. IMPORTANCE Our manuscript describes a SARS-CoV-2 outbreak investigation in an Ontario tertiary care hospital. We use traditional contract tracing paired with whole-genome sequencing to facilitate an understanding of the evolutionary history and transmission dynamics of this SARS-CoV-2 outbreak in a clinical setting. These advancements have enabled the incorporation of phylogenetics and genomic epidemiology into the understanding of clinical outbreaks. We show that genomic epidemiology can help to explore the genetic evolution of a pathogen in real time, enabling the identification of the index case and helping understand its transmission dynamics to develop better strategies to prevent future spread of SARS-CoV-2 in congregate, clinical settings such as hospitals.

Diagnostics and analysis of SARS-CoV-2: current status, recent advances, challenges and perspectives

The disastrous spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has induced severe public healthcare issues and weakened the global economy significantly. Although SARS-CoV-2 infection is not as fatal as the initial outbreak, many infected victims suffer from long COVID. Therefore, rapid and large-scale testing is critical in managing patients and alleviating its transmission. Herein, we review the recent advances in techniques to detect SARS-CoV-2. The sensing principles are detailed together with their application domains and analytical performances. In addition, the advantages and limits of each method are discussed and analyzed. Besides molecular diagnostics and antigen and antibody tests, we also review neutralizing antibodies and emerging SARS-CoV-2 variants. Further, the characteristics of the mutational locations in the different variants with epidemiological features are summarized. Finally, the challenges and possible strategies are prospected to develop new assays to meet different diagnostic needs. Thus, this comprehensive and systematic review of SARS-CoV-2 detection technologies may provide insightful guidance and direction for developing tools for the diagnosis and analysis of SARS-CoV-2 to support public healthcare and effective long-term pandemic management and control.

CBPH assay for the highest sensitive detection of SARS-COV-2 in the semen

BACKGROUND

Great concerns have been raised on SARS-CoV-2 impact on men's andrological well-being, and many studies have attempted to determine whether SARS-CoV-2 is present in the semen and till now the data are unclear and somehow ambiguous. However, these studies used quantitative real-time (qRT) PCR, which is not sufficiently sensitive to detect nucleic acids in clinical samples with a low viral load.

METHODS

The clinical performance of various nucleic acid detection methods (qRT-PCR, OSN-qRT-PCR, cd-PCR, and CBPH) was assessed for SARS-CoV-2 using 236 clinical samples from laboratory-confirmed COVID-19 cases. Then, the presence of SARS-CoV-2 in the semen of 12 recovering patients was investigated using qRT-PCR, OSN-qRT-PCR, cd-PCR, and CBPH in parallel using 24 paired semen, blood, throat swab, and urine samples.

RESULTS

The sensitivity and specificity along with AUC of CBPH was markedly higher than the other 3methods. Although qRT-PCR, OSN-qRT-PCR and cdPCR detected no SARS-CoV-2 RNA in throat swab, blood, urine, and semen samples of the 12 patients, CBPH detected the presence of SARS-CoV-2 genome fragments in semen samples, but not in paired urine samples, of 3 of 12 patients. The existing SARS-CoV-2 genome fragments were metabolized over time.

CONCLUSIONS

Both OSN-qRT-PCR and cdPCR had better performance than qRT-PCR, and CBPH had the highest diagnostic performance in detecting SARS-CoV-2, which contributed the most improvement to the determination of the critical value in gray area samples with low vrial load, which then provides a rational screening strategy for studying the clearance of coronavirus in the semen over time in patients recovering from COVID-19. Although the presence of SARS-CoV-2 fragments in the semen was demonstrated by CBPH, COVID-19 is unlikely to be sexually transmitted from male partners for at least 3 months after hospital discharge.

The Quantification of Spike Proteins in the Inactivated SARS-CoV-2 Vaccines of the Prototype, Delta, and Omicron Variants by LC-MS

Developing variant vaccines or multivalent vaccines is a feasible way to address the epidemic as the SARS-CoV-2 variants of concern (VOCs) posed an increased risk to global public health. The spike protein of the SARS-CoV-2 virus was usually used as the main antigen in many types of vaccines to produce neutralizing antibodies against the virus. However, the spike (S) proteins of different variants were only differentiated by a few amino acids, making it difficult to obtain specific antibodies that can distinguish different VOCs, thereby challenging the accurate distinction and quantification of the variants using immunological methods such as ELISA. Here, we established a method based on LC-MS to quantify the S proteins in inactivated monovalent vaccines or trivalent vaccines (prototype, Delta, and Omicron strains). By analyzing the S protein sequences of the prototype, Delta, and Omicron strains, we identified peptides that were different and specific among the three strains and synthesized them as references. The synthetic peptides were isotopically labeled as internal targets. Quantitative analysis was performed by calculating the ratio between the reference and internal target. The verification results have shown that the method we established had good specificity, accuracy, and precision. This method can not only accurately quantify the inactivated monovalent vaccine but also could be applied to each strain in inactivated trivalent SARS-CoV-2 vaccines. Hence, the LC-MS method established in this study can be applied to the quality control of monovalent and multivalent SARS-CoV-2 variation vaccines. By enabling more accurate quantification, it will help to improve the protection of the vaccine to some extent.

A review of current effective COVID-19 testing methods and quality control

COVID-19 is a highly infectious disease caused by the SARS-CoV-2 virus, which primarily affects the respiratory system and can lead to severe illness. The virus is extremely contagious, early and accurate diagnosis of SARS-CoV-2 is crucial to contain its spread, to provide prompt treatment, and to prevent complications. Currently, the reverse transcriptase polymerase chain reaction (RT-PCR) is considered to be the gold standard for detecting COVID-19 in its early stages. In addition, loop-mediated isothermal amplification (LMAP), clustering rule interval short palindromic repeats (CRISPR), colloidal gold immunochromatographic assay (GICA), computed tomography (CT), and electrochemical sensors are also common tests. However, these different methods vary greatly in terms of their detection efficiency, specificity, accuracy, sensitivity, cost, and throughput. Besides, most of the current detection methods are conducted in central hospitals and laboratories, which is a great challenge for remote and underdeveloped areas. Therefore, it is essential to review the advantages and disadvantages of different COVID-19 detection methods, as well as the technology that can enhance detection efficiency and improve detection quality in greater details.

A comparison of five Illumina, Ion Torrent, and nanopore sequencing technology-based approaches for whole genome sequencing of SARS-CoV-2

Rapid identification of the rise and spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern remains critical for monitoring of the efficacy of diagnostics, therapeutics, vaccines, and control strategies. A wide range of SARS-CoV-2 next-generation sequencing (NGS) methods have been developed over the last years, but cross-sequence technology benchmarking studies have been scarce. In the current study, 26 clinical samples were sequenced using five protocols: AmpliSeq SARS-CoV-2 (Illumina), EasySeq RC-PCR SARS-CoV-2 (Illumina/NimaGen), Ion AmpliSeq SARS-CoV-2 (Thermo Fisher), custom primer sets (Oxford Nanopore Technologies (ONT)), and capture probe-based viral metagenomics (Roche/Illumina). Studied parameters included genome coverage, depth of coverage, amplicon distribution, and variant calling. The median SARS-CoV-2 genome coverage of samples with cycle threshold (Ct) values of 30 and lower ranged from 81.6 to 99.8% for, respectively, the ONT protocol and Illumina AmpliSeq protocol. Correlation of coverage with PCR Ct values varied per protocol. Amplicon distribution signatures differed across the methods, with peak differences of up to 4 log₁₀ at disbalanced positions in samples with high viral loads (Ct values ≤ 23). Phylogenetic analyses of consensus sequences showed clustering independent of the workflow used. The proportion of SARS-CoV-2 reads in relation to background sequences, as a (cost-)efficiency metric, was the highest for the EasySeq protocol. The hands-on time was the lowest when using EasySeq and ONT protocols, with the latter additionally having the shortest sequence runtime. In conclusion, the studied protocols differed on a variety of the studied metrics. This study provides data that assist laboratories when selecting protocols for their specific setting.

The P323L substitution in the SARS-CoV-2 polymerase (NSP12) confers a selective advantage during infection

BACKGROUND

The mutational landscape of SARS-CoV-2 varies at the dominant viral genome sequence and minor genomic variant population. During the COVID-19 pandemic, an early substitution in the genome was the D614G change in the spike protein, associated with an increase in transmissibility. Genomes with D614G are accompanied by a P323L substitution in the viral polymerase (NSP12). However, P323L is not thought to be under strong selective pressure.

RESULTS

Investigation of P323L/D614G substitutions in the population shows rapid emergence during the containment phase and early surge phase during the first wave. These substitutions emerge from minor genomic variants which become dominant viral genome sequence. This is investigated in vivo and in vitro using SARS-CoV-2 with P323 and D614 in the dominant genome sequence and L323 and G614 in the minor variant population. During infection, there is rapid selection of L323 into the dominant viral genome sequence but not G614. Reverse genetics is used to create two viruses (either P323 or L323) with the same genetic background. L323 shows greater abundance of viral RNA and proteins and a smaller plaque morphology than P323.

CONCLUSIONS

These data suggest that P323L is an important contribution in the emergence of variants with transmission advantages. Sequence analysis of viral populations suggests it may be possible to predict the emergence of a new variant based on tracking the frequency of minor variant genomes. The ability to predict an emerging variant of SARS-CoV-2 in the global landscape may aid in the evaluation of medical countermeasures and non-pharmaceutical interventions.

The clinical characteristics of pediatric patients infected by SARS-CoV-2 Omicron variant and whole viral genome sequencing analysis

Pediatric population was generally less affected clinically by SARS-CoV-2 infection. Few pediatric cases of COVID-19 have been reported compared to those reported in infected adults. However, a rapid increase in the hospitalization rate of SARS-CoV-2 infected pediatric patients was observed during Omicron variant dominated COVID-19 outbreak. In this study, we analyzed the B.1.1.529 (Omicron) genome sequences collected from pediatric patients by whole viral genome amplicon sequencing using Illumina next generation sequencing platform, followed by phylogenetic analysis. The demographic, epidemiologic and clinical data of these pediatric patients are also reported in this study. Fever, cough, running nose, sore throat and vomiting were the more commonly reported symptoms in children infected by Omicron variant. A novel frameshift mutation was found in the ORF1b region (NSP12) of the genome of Omicron variant. Seven mutations were identified in the target regions of the WHO listed SARS-CoV-2 primers and probes. On protein level, eighty-three amino acid substitutions and fifteen amino acid deletions were identified. Our results indicate that asymptomatic infection and transmission among children infected by Omicron subvariants BA.2.2 and BA.2.10.1 are not common. Omicron may have different pathogenesis in pediatric population.

Monitoring of SARS-CoV-2 Infection in Ragusa Area: Next Generation Sequencing and Serological Analysis

The coronavirus disease 19 (COVID-19) post pandemic evolution is correlated to the development of new variants. Viral genomic and immune response monitoring are fundamental to the surveillance of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Since 1 January to 31 July 2022, we monitored the SARS-CoV-2 variants trend in Ragusa area sequencing n.600 samples by next generation sequencing (NGS) technology: n.300 were healthcare workers (HCWs) of ASP Ragusa. The evaluation of anti-Nucleocapside (N), receptor-binding domain (RBD), the two subunit of S protein (S1 and S2) IgG levels in 300 exposed vs. 300 unexposed HCWs to SARS-CoV-2 was performed. Differences in immune response and clinical symptoms related to the different variants were investigated. The SARS-CoV-2 variants trend in Ragusa area and in Sicily region were comparable. BA.1 and BA.2 were the most representative variants, whereas the diffusion of BA.3 and BA.4 affected some places of the region. Although no correlation was found between variants and clinical manifestations, anti-N and anti-S2 levels were positively correlated with an increase in the symptoms number. SARS-CoV-2 infection induced a statistically significant enhancement in antibody titers compared to that produced by SARS-CoV-2 vaccine administration. In post-pandemic period, the evaluation of anti-N IgG could be used as an early marker to identify asymptomatic subjects.

A case for investment in clinical metagenomics in low-income and middle-income countries

Clinical metagenomics is the diagnostic approach with the broadest capacity to detect both known and novel pathogens. Clinical metagenomics is costly to run and requires infrastructure, but the use of next-generation sequencing for SARS-CoV-2 molecular epidemiology in low-income and middle-income countries (LMICs) offers an opportunity to direct this infrastructure to the establishment of clinical metagenomics programmes. Local implementation of clinical metagenomics is important to create relevant systems and evaluate cost-effective methodologies for its use, as well as to ensure that reference databases and result interpretation tools are appropriate to local epidemiology. Rational implementation, based on the needs of LMICs and the available resources, could ultimately improve individual patient care in instances in which available diagnostics are inadequate and supplement emerging infectious disease surveillance systems to ensure the next pandemic pathogen is quickly identified.

Combined Use of RT-qPCR and NGS for Identification and Surveillance of SARS-CoV-2 Variants of Concern in Residual Clinical Laboratory Samples in Miami-Dade County, Florida

Over the course of the COVID-19 pandemic, SARS-CoV-2 variants of concern (VOCs) with increased transmissibility and immune escape capabilities, such as Delta and Omicron, have triggered waves of new COVID-19 infections worldwide, and Omicron subvariants continue to represent a global health concern. Tracking the prevalence and dynamics of VOCs has clinical and epidemiological significance and is essential for modeling the progression and evolution of the COVID-19 pandemic. Next generation sequencing (NGS) is recognized as the gold standard for genomic characterization of SARS-CoV-2 variants, but it is labor and cost intensive and not amenable to rapid lineage identification. Here we describe a two-pronged approach for rapid, cost-effective surveillance of SARS-CoV-2 VOCs by combining reverse-transcriptase quantitative polymerase chain reaction (RT-qPCR) and periodic NGS with the ARTIC sequencing method. Variant surveillance by RT-qPCR included the commercially available TaqPath COVID-19 Combo Kit to track S-gene target failure (SGTF) associated with the spike protein deletion H69-V70, as well as two internally designed and validated RT-qPCR assays targeting two N-terminal-domain (NTD) spike gene deletions, NTD156-7 and NTD25-7. The NTD156-7 RT-qPCR assay facilitated tracking of the Delta variant, while the NTD25-7 RT-qPCR assay was used for tracking Omicron variants, including the BA.2, BA.4, and BA.5 lineages. In silico validation of the NTD156-7 and NTD25-7 primers and probes compared with publicly available SARS-CoV-2 genome databases showed low variability in regions corresponding to oligonucleotide binding sites. Similarly, in vitro validation with NGS-confirmed samples showed excellent correlation. RT-qPCR assays allow for near-real-time monitoring of circulating and emerging variants allowing for ongoing surveillance of variant dynamics in a local population. By performing periodic sequencing of variant surveillance by RT-qPCR methods, we were able to provide ongoing validation of the results obtained by RT-qPCR screening. Rapid SARS-CoV-2 variant identification and surveillance by this combined approach served to inform clinical decisions in a timely manner and permitted better utilization of sequencing resources.

Application of CRISPR-Based Human and Bacterial Ribosomal RNA Depletion for SARS-CoV-2 Shotgun Metagenomic Sequencing

OBJECTIVES

The aim of this study is to evaluate the effectiveness of a CRISPR-based human and bacterial ribosomal RNA (rRNA) depletion kit (JUMPCODE Genomics) on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) shotgun metagenomic sequencing in weakly positive respiratory samples.

METHODS

Shotgun metagenomics was performed on 40 respiratory specimens collected from solid organ transplant patients and deceased intensive care unit patients at UCLA Medical Center in late 2020 to early 2021. Human and bacterial rRNA depletion was performed on remnant library pools prior to sequencing by Illumina MiSeq. Data quality was analyzed using Geneious Prime, whereas the identification of SARS-CoV-2 variants and lineages was determined by Pangolin.

RESULTS

The average genome coverage of the rRNA-depleted respiratory specimens increased from 72.55% to 93.71% in overall samples and from 29.3% to 83.3% in 15 samples that failed to achieve sufficient genome coverage using the standard method. Moreover, rRNA depletion enhanced genome coverage to over 85% in 11 (73.3%) of 15 low viral load samples with cycle threshold values up to 35, resulting in the identification of genotypes.

CONCLUSION

The CRISPR-based human and bacterial rRNA depletion enhanced the sensitivity of SARS-CoV-2 shotgun metagenomic sequencing, especially in low viral load samples.

Whole-genome single molecule real-time sequencing of SARS-CoV-2 Omicron

New variants and genetic mutations of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome can only be identified using accurate sequencing methods. Single molecule real-time (SMRT) sequencing has been used to characterize Alpha and Delta variants, but not Omicron variants harboring numerous mutations in the SARS-CoV-2 genome. This study assesses the performance of a target capture SMRT sequencing protocol for whole genome sequencing (WGS) of SARS-CoV-2 Omicron variants and compared it to that of an amplicon SMRT sequencing protocol optimized for Omicron variants. The failure rate of the target capture protocol (6%) was lower than that of the amplicon protocol (34%, p < 0.001) on our data set, and the median genome coverage with the target capture protocol (98.6% [interquartile range (IQR): 86-99.4]) was greater than that with the amplicon protocol (76.6% [IQR: 66-89.6], [p < 0.001]). The percentages of samples with >95% whole genome coverage were 64% with the target capture protocol and 19% with the amplicon protocol (p < 0.05). The clades of 96 samples determined with both protocols were 93% concordant and the lineages of 59 samples were 100% concordant. Thus, target capture SMRT sequencing appears to be an efficient method for WGS, genotyping and detecting mutations of SARS-CoV-2 Omicron variants.

Whole-genome sequencing of SARS-CoV-2: Comparison of target capture and amplicon single molecule real-time sequencing protocols

Fast, accurate sequencing methods are needed to identify new variants and genetic mutations of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome. Single-molecule real-time (SMRT) Pacific Biosciences (PacBio) provides long, highly accurate sequences by circular consensus reads. This study compares the performance of a target capture SMRT PacBio protocol for whole-genome sequencing (WGS) of SARS-CoV-2 to that of an amplicon PacBio SMRT sequencing protocol. The median genome coverage was higher (p < 0.05) with the target capture protocol (99.3% [interquartile range, IQR: 96.3-99.5]) than with the amplicon protocol (99.3% [IQR: 69.9-99.3]). The clades of 65 samples determined with both protocols were 100% concordant. After adjusting for C_t values, S gene coverage was higher with the target capture protocol than with the amplicon protocol. After stratification on C_t values, higher S gene coverage with the target capture protocol was observed only for samples with C_t  > 17 (p < 0.01). PacBio SMRT sequencing protocols appear to be suitable for WGS, genotyping, and detecting mutations of SARS-CoV-2.

Two Years into the COVID-19 Pandemic: Lessons Learned

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a highly transmissible and virulent human-infecting coronavirus that emerged in late December 2019 in Wuhan, China, causing a respiratory disease called coronavirus disease 2019 (COVID-19), which has massively impacted global public health and caused widespread disruption to daily life. The crisis caused by COVID-19 has mobilized scientists and public health authorities across the world to rapidly improve our knowledge about this devastating disease, shedding light on its management and control, and spawned the development of new countermeasures. Here we provide an overview of the state of the art of knowledge gained in the last 2 years about the virus and COVID-19, including its origin and natural reservoir hosts, viral etiology, epidemiology, modes of transmission, clinical manifestations, pathophysiology, diagnosis, treatment, prevention, emerging variants, and vaccines, highlighting important differences from previously known highly pathogenic coronaviruses. We also discuss selected key discoveries from each topic and underline the gaps of knowledge for future investigations.

Phylogenetic and phylodynamic approaches to understanding and combating the early SARS-CoV-2 pandemic

Determining the transmissibility, prevalence and patterns of movement of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is central to our understanding of the impact of the pandemic and to the design of effective control strategies. Phylogenies (evolutionary trees) have provided key insights into the international spread of SARS-CoV-2 and enabled investigation of individual outbreaks and transmission chains in specific settings. Phylodynamic approaches combine evolutionary, demographic and epidemiological concepts and have helped track virus genetic changes, identify emerging variants and inform public health strategy. Here, we review and synthesize studies that illustrate how phylogenetic and phylodynamic techniques were applied during the first year of the pandemic, and summarize their contributions to our understanding of SARS-CoV-2 transmission and control.

Cross-reactive immunity against SARS-CoV-2 N protein in Central and West Africa precedes the COVID-19 pandemic

Early predictions forecasted large numbers of severe acute respiratory syndrome coronavirus (SARS-CoV-2) cases and associated deaths in Africa. To date, Africa has been relatively spared. Various hypotheses were postulated to explain the lower than anticipated impact on public health in Africa. However, the contribution of pre-existing immunity is yet to be investigated. In this study, the presence of antibodies against SARS-CoV-2 spike (S) and nucleocapsid (N) proteins in pre-pandemic samples from Africa, Europe, South and North America was examined by ELISA. The protective efficacy of N specific antibodies isolated from Central African donors was tested by in vitro neutralization and in a mouse model of SARS-CoV-2 infection. Antibodies against SARS-CoV-2 S and N proteins were rare in all populations except in Gabon and Senegal where N specific antibodies were prevalent. However, these antibodies failed to neutralize the virus either in vitro or in vivo. Overall, this study indicates that cross-reactive immunity against SARS-CoV-2 N protein was present in Africa prior to the pandemic. However, this pre-existing humoral immunity does not impact viral fitness in rodents suggesting that other human immune defense mechanisms could be involved. In Africa, seroprevalence studies using the N protein are over-estimating SARS-CoV-2 circulation.

Analysis of SARS-CoV-2 known and novel subgenomic mRNAs in cell culture, animal model, and clinical samples using LeTRS, a bioinformatic tool to identify unique sequence identifiers

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a complex strategy for the transcription of viral subgenomic mRNAs (sgmRNAs), which are targets for nucleic acid diagnostics. Each of these sgmRNAs has a unique 5' sequence, the leader-transcriptional regulatory sequence gene junction (leader-TRS junction), that can be identified using sequencing. High-resolution sequencing has been used to investigate the biology of SARS-CoV-2 and the host response in cell culture and animal models and from clinical samples. LeTRS, a bioinformatics tool, was developed to identify leader-TRS junctions and can be used as a proxy to quantify sgmRNAs for understanding virus biology. LeTRS is readily adaptable for other coronaviruses such as Middle East respiratory syndrome coronavirus or a future newly discovered coronavirus. LeTRS was tested on published data sets and novel clinical samples from patients and longitudinal samples from animal models with coronavirus disease 2019. LeTRS identified known leader-TRS junctions and identified putative novel sgmRNAs that were common across different mammalian species. This may be indicative of an evolutionary mechanism where plasticity in transcription generates novel open reading frames, which can then subject to selection pressure. The data indicated multiphasic abundance of sgmRNAs in two different animal models. This recapitulates the relative sgmRNA abundance observed in cells at early points in infection but not at late points. This pattern is reflected in some human nasopharyngeal samples and therefore has implications for transmission models and nucleic acid-based diagnostics. LeTRS provides a quantitative measure of sgmRNA abundance from sequencing data. This can be used to assess the biology of SARS-CoV-2 (or other coronaviruses) in clinical and nonclinical samples, especially to evaluate different variants and medical countermeasures that may influence viral RNA synthesis.

Rapid and Accurate Detection of SARS Coronavirus 2 by Nanopore Amplicon Sequencing

Objective

We aimed to evaluate the performance of nanopore amplicon sequencing detection for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in clinical samples.

Method

We carried out a single-center, prospective cohort study in a Wuhan hospital and collected a total of 86 clinical samples, including 54 pharyngeal swabs, 31 sputum samples, and 1 fecal sample, from 86 patients with coronavirus disease 2019 (COVID-19) from Feb 20 to May 15, 2020. We performed parallel detection with nanopore-based genome amplification and sequencing (NAS) on the Oxford Nanopore Technologies (ONT) minION platform and routine reverse transcription quantitative polymerase chain reaction (RT-qPCR). In addition, 27 negative control samples were detected using the two methods. The sensitivity and specificity of NAS were evaluated and compared with those of RT-qPCR.

Results

The viral read number and reference genome coverage were both significantly different between the two groups of samples, and the latter was a better indicator for SARS-CoV-2 detection. Based on the reference genome coverage, NAS revealed both high sensitivity (96.5%) and specificity (100%) compared with RT-qPCR (80.2 and 96.3%, respectively), although the samples had been stored for half a year before the detection. The total time cost was less than 15 h, which was acceptable compared with that of RT-qPCR (∼2.5 h). In addition, the reference genome coverage of the viral reads was in line with the cycle threshold value of RT-qPCR, indicating that this number could also be used as an indicator of the viral load in a sample. The viral load in sputum might be related to the severity of the infection, particularly in patients within 4 weeks after onset of clinical manifestations, which could be used to evaluate the infection.

Conclusion

Our results showed the high sensitivity and specificity of the NAS method for SARS-CoV-2 detection compared with RT-qPCR. The sequencing results were also used as an indicator of the viral load to display the viral dynamics during infection. This study proved the wide application prospect of nanopore sequencing detection for SARS-CoV-2 and may more knowledge about the clinical characteristics of COVID-19.

Two DNA vaccines protect against severe disease and pathology due to SARS-CoV-2 in Syrian hamsters

The SARS-CoV-2 pandemic is an ongoing threat to global health, and wide-scale vaccination is an efficient method to reduce morbidity and mortality. We designed and evaluated two DNA plasmid vaccines, based on the pIDV-II system, expressing the SARS-CoV-2 spike gene, with or without an immunogenic peptide, in mice, and in a Syrian hamster model of infection. Both vaccines demonstrated robust immunogenicity in BALB/c and C57BL/6 mice. Additionally, the shedding of infectious virus and the viral burden in the lungs was reduced in immunized hamsters. Moreover, high-titers of neutralizing antibodies with activity against multiple SARS-CoV-2 variants were generated in immunized animals. Vaccination also protected animals from weight loss during infection. Additionally, both vaccines were effective at reducing both pulmonary and extrapulmonary pathology in vaccinated animals. These data show the potential of a DNA vaccine for SARS-CoV-2 and suggest further investigation in large animal and human studies could be pursued.

Identification of SARS-CoV-2 variants using viral sequencing for the Centers for Disease Control and Prevention genomic surveillance program

BACKGROUND

The Centers for Disease Control and Prevention contracted with laboratories to sequence the SARS-CoV-2 genome from positive samples across the United States to enable public health officials to investigate the impact of variants on disease severity as well as the effectiveness of vaccines and treatment. Herein we present the initial results correlating RT-PCR quality control metrics with sample collection and sequencing methods from full SARS-CoV-2 viral genomic sequencing of 24,441 positive patient samples between April and June 2021.

METHODS

RT-PCR confirmed (N Gene Ct value < 30) positive patient samples, with nucleic acid extracted from saliva, nasopharyngeal and oropharyngeal swabs were selected for viral whole genome SARS-CoV-2 sequencing. Sequencing was performed using Illumina COVIDSeq™ protocol on either the NextSeq550 or NovaSeq6000 systems. Informatic variant calling, and lineage analysis were performed using DRAGEN COVID Lineage applications on Illumina's Basespace cloud analytical system. All sequence data and variant calls were uploaded to NCBI and GISAID.

RESULTS

An association was observed between higher sequencing coverage, quality, and samples with a lower Ct value, with < 27 being optimal, across both sequencing platforms and sample collection methods. Both nasopharyngeal swabs and saliva samples were found to be optimal samples of choice for SARS-CoV-2 surveillance sequencing studies, both in terms of strain identification and sequencing depth of coverage, with NovaSeq 6000 providing higher coverage than the NextSeq 550. The most frequent variants identified were the B.1.617.2 Delta (India) and P.1 Gamma (Brazil) variants in the samples sequenced between April 2021 and June 2021. At the time of submission, the most common variant > 99% of positives sequenced was Omicron.

CONCLUSION

These initial analyses highlight the importance of sequencing platform, sample collection methods, and RT-PCR Ct values in guiding surveillance efforts. These surveillance studies evaluating genetic changes of SARS-CoV-2 have been identified as critical by the CDC that can affect many aspects of public health including transmission, disease severity, diagnostics, therapeutics, and vaccines.

PipeCoV: a pipeline for SARS-CoV-2 genome assembly, annotation and variant identification

Motivation

Since the identification of the novel coronavirus (SARS-CoV-2), the scientific community has made a huge effort to understand the virus biology and to develop vaccines. Next-generation sequencing strategies have been successful in understanding the evolution of infectious diseases as well as facilitating the development of molecular diagnostics and treatments. Thousands of genomes are being generated weekly to understand the genetic characteristics of this virus. Efficient pipelines are needed to analyze the vast amount of data generated. Here we present a new pipeline designed for genomic analysis and variant identification of the SARS-CoV-2 virus.

Results

PipeCoV shows better performance when compared to well-established SARS-CoV-2 pipelines, with a lower content of Ns and higher genome coverage when compared to the Wuhan reference. It also provides a variant report not offered by other tested pipelines.

Availability

https://github.com/alvesrco/pipecov.

Immunogenicity of convalescent and vaccinated sera against clinical isolates of ancestral SARS-CoV-2, beta, delta, and omicron variants

Background

SARS-CoV-2 Omicron variant of concern (VOC) has evolved multiple mutations within the spike protein, raising concerns of increased antibody evasion. In this study, we assessed the neutralization potential of COVID-19 convalescent sera and sera from vaccinated individuals against ancestral SARS-CoV-2 and VOCs.

Methods

The neutralizing activity of sera from 65 coronavirus disease (COVID-19) vaccine recipients and convalescent individuals against clinical isolates of ancestral SARS-CoV-2 and Beta, Delta, and Omicron VOCs was assessed using a micro-neutralization assay.

Findings

Convalescent sera from unvaccinated individuals infected by the ancestral virus demonstrated reduced neutralization against Beta and Omicron VOCs. Sera from individuals that received three doses of the Pfizer or Moderna vaccines demonstrated reduced neutralization of the Omicron variant relative to ancestral SARS-CoV-2. Sera from individuals that were naturally infected with ancestral SARS-CoV-2 and subsequently received two doses of the Pfizer vaccine induced significantly higher neutralizing antibody levels against ancestral virus and all VOCs. Infection alone, either with ancestral SARS-CoV-2 or the Delta variant, was not sufficient to induce high neutralizing antibody titers against Omicron.

Conclusions

In summary, we demonstrate that convalescent and vaccinated sera display varying levels of SARS-CoV-2 VOC neutralization. Data from this study will inform booster vaccination strategies against SARS-CoV-2 VOCs.

Funding

This research was funded by the Canadian Institutes of Health Research (CIHR). VIDO receives operational funding from the Government of Saskatchewan through Innovation Saskatchewan and the Ministry of Agriculture and from the Canada Foundation for Innovation through the Major Science Initiatives for its CL3 facility.

SARS-CoV-2 Omicron variant of concern (VOC) has evolved multiple mutations within the spike protein, raising concerns of increased antibody evasion. In this study, we quantified neutralizing antibody levels in convalescent and vaccinated sera against ancestral SARS-CoV-2 and VOCs. Convalescent sera had lower neutralizing antibody levels against the Omicron VOC. Two doses of an mRNA vaccine following infection induced high levels of neutralizing antibodies against all VOCs, including Omicron. Three doses of authorized mRNA vaccines induced detectable but lower levels of neutralizing antibodies against VOCs in long-term care residents. Data from our study, along with other published studies, support the utility of third vaccine doses and will help inform future booster vaccination strategies to tackle the ongoing COVID-19 pandemic.

Multiplex PCR Assays for Identifying all Major Severe Acute Respiratory Syndrome Coronavirus 2 Variants

Variants of concern (VOC) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), including alpha, beta, gamma, delta, and omicron, threaten to prolong the pandemic, leading to more global morbidity and mortality. Genome sequencing is the mainstay of tracking the evolution of the virus, but is costly, slow, and not easily accessible. Multiplex quantitative RT-PCR assays for SARS-CoV-2 have been developed that identify all VOCs as well as other mutations of interest in the viral genome, nine mutations in total, using single-nucleotide discriminating molecular beacons. The presented variant molecular beacon assays showed a limit of detection of 50 copies of viral RNA, with 100% specificity. Twenty-six SARS-CoV-2-positive patient samples were blinded and tested using a two-tube assay. When testing patient samples, the assay was in full agreement with results from deep sequencing with a sensitivity and specificity of 100% (26 of 26). We have used our design methodology to rapidly design an assay that detects the new omicron variant. This omicron assay was used to accurately identify this variant in 17 of 33 additional patient samples. These quantitative RT-PCR assays identify all currently circulating VOCs of SARS-CoV-2, as well as other important mutations in the spike protein coding sequence. These assays can be easily implemented on broadly available five-color thermal cyclers and will help track the spread of these variants.

Advances in laboratory detection methods and technology application of SARS-CoV-2

At present, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is raging worldwide, and the coronavirus disease 2019 outbreak caused by SARS-CoV-2 seriously threatens the life and health of all humankind. There is no specific medicine for novel coronavirus yet. So, laboratory diagnoses of novel coronavirus as soon as possible and isolation of the source of infection play a vital role in preventing and controlling the epidemic. Therefore, selecting appropriate detection techniques and methods is particularly important to improve the efficiency of disease diagnosis and treatment and to curb the outbreak of infectious diseases. In this paper, virus nucleic acid, protein, and serum immunology were reviewed to provide a reference for further developing virus detection technology to provide better prevention and treatment strategies and research ideas for clinicians and researchers.

The dynamics of circulating SARS-CoV-2 lineages in Bogor and surrounding areas reflect variant shifting during the first and second waves of COVID-19 in Indonesia

Background

Indonesia is one of the Southeast Asian countries with high case numbers of COVID-19 with up to 4.2 million confirmed cases by 29 October 2021. Understanding the genome of SARS-CoV-2 is crucial for delivering public health intervention as certain variants may have different attributes that can potentially affect their transmissibility, as well as the performance of diagnostics, vaccines, and therapeutics.

Objectives

We aimed to investigate the dynamics of circulating SARS-CoV-2 variants over a 15-month period in Bogor and its surrounding areas in correlation with the first and second wave of COVID-19 in Indonesia.

Methods

Nasopharyngeal and oropharyngeal swab samples collected from suspected patients from Bogor, Jakarta and Tangerang were confirmed for SARS-CoV-2 infection with RT-PCR. RNA samples of those confirmed patients were subjected to whole genome sequencing using the ARTIC Network protocol and sequencer platform from Oxford Nanopore Technologies (ONT).

Results

We successfully identified 16 lineages and six clades out of 202 samples (male n = 116, female n = 86). Genome analysis revealed that Indonesian lineage B.1.466.2 dominated during the first wave (n = 48, 23.8%) while Delta variants (AY.23, AY.24, AY.39, AY.42, AY.43 dan AY.79) were dominant during the second wave (n = 53, 26.2%) following the highest number of confirmed cases in Indonesia. In the spike protein gene, S_D614G and S_P681R changes were dominant in both B.1.466.2 and Delta variants, while N439K was only observed in B.1.466.2 (n = 44) and B.1.470 (n = 1). Additionally, the S_T19R, S_E156G, S_F157del, S_R158del, S_L452R, S_T478K, S_D950N and S_V1264L changes were only detected in Delta variants, consistent with those changes being characteristic of Delta variants in general.

Conclusions

We demonstrated a shift in SARS-CoV-2 variants from the first wave of COVID-19 to Delta variants in the second wave, during which the number of confirmed cases surpassed those in the first wave of COVID-19 pandemic. Higher proportion of unique mutations detected in Delta variants compared to the first wave variants indicated potential mutational effects on viral transmissibility that correlated with a higher incidence of confirmed cases. Genomic surveillance of circulating variants, especially those with higher transmissibility, should be continuously conducted to rapidly inform decision making and support outbreak preparedness, prevention, and public health response.

High-Integrity Sequencing of Spike Gene for SARS-CoV-2 Variant Determination

For tiling of the SARS-CoV-2 genome, the ARTIC Network provided a V4 protocol using 99 pairs of primers for amplicon production and is currently the widely used amplicon-based approach. However, this technique has regions of low sequence coverage and is labour-, time-, and cost-intensive. Moreover, it requires 14 pairs of primers in two separate PCRs to obtain spike gene sequences. To overcome these disadvantages, we proposed a single PCR to efficiently detect spike gene mutations. We proposed a bioinformatic protocol that can process FASTQ reads into spike gene consensus sequences to accurately call spike protein variants from sequenced samples or to fairly express the cases of missing amplicons. We evaluated the in silico detection rate of primer sets that yield amplicon sizes of 400, 1200, and 2500 bp for spike gene sequencing of SARS-CoV-2 to be 59.49, 76.19, and 92.20%, respectively. The in silico detection rate of our proposed single PCR primers was 97.07%. We demonstrated the robustness of our analytical protocol against 3000 Oxford Nanopore sequencing runs of distinct datasets, thus ensuring high-integrity sequencing of spike genes for variant SARS-CoV-2 determination. Our protocol works well with the data yielded from versatile primer designs, making it easy to determine spike protein variants.

Surface and Air Contamination With Severe Acute Respiratory Syndrome Coronavirus 2 From Hospitalized Coronavirus Disease 2019 Patients in Toronto, Canada, March-May 2020

BACKGROUND

We determined the burden of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in air and on surfaces in rooms of patients hospitalized with coronavirus disease 2019 (COVID-19) and investigated patient characteristics associated with SARS-CoV-2 environmental contamination.

METHODS

Nasopharyngeal swabs, surface, and air samples were collected from the rooms of 78 inpatients with COVID-19 at 6 acute care hospitals in Toronto from March to May 2020. Samples were tested for SARS-CoV-2 ribonucleic acid (RNA), cultured to determine potential infectivity, and whole viral genomes were sequenced. Association between patient factors and detection of SARS-CoV-2 RNA in surface samples were investigated.

RESULTS

Severe acute respiratory syndrome coronavirus 2 RNA was detected from surfaces (125 of 474 samples; 42 of 78 patients) and air (3 of 146 samples; 3 of 45 patients); 17% (6 of 36) of surface samples from 3 patients yielded viable virus. Viral sequences from nasopharyngeal and surface samples clustered by patient. Multivariable analysis indicated hypoxia at admission, polymerase chain reaction-positive nasopharyngeal swab (cycle threshold of ≤30) on or after surface sampling date, higher Charlson comorbidity score, and shorter time from onset of illness to sampling date were significantly associated with detection of SARS-CoV-2 RNA in surface samples.

CONCLUSIONS

The infrequent recovery of infectious SARS-CoV-2 virus from the environment suggests that the risk to healthcare workers from air and near-patient surfaces in acute care hospital wards is likely limited.

SARS-CoV-2 in hospital air as revealed by comprehensive respiratory viral panel sequencing

BACKGROUND

Nosocomially acquired severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) infection has become the most significant pandemic of our lifetime. Though its transmission was essentially attributed to droplets from an infected person, with recent advancements in knowledge, aerosol transmission seems to be a viable pathway, as well. Because of the lower biological load in ambient aerosol, detection of SARS-CoV-2 is challenging. A few recent attempts of sampling large aerosol volumes and using next-generation sequencing (NGS) to detect the presence of SARS-CoV-2 in the air at very low levels gave positive results. These results suggest the potential of using this technique to detect the presence of SARS-CoV-2 and use it as an early warning signal for possible outbreak or recurrence of coronavirus disease 2019 (COVID-19).

AIM

To assess efficacy of comprehensive respiratory viral panel (CRVP) sequencing and RT-PCR for low-level identification of SARS-CoV-2 and other respiratory viruses in indoor air.

METHODS

A large volume of indoor aerosol samples from three major hospitals involved in COVID-19 care in Kuwait was collected. Viral RNA was isolated and subjected to comprehensive respiratory viral panel sequencing (CRVP) as per the standard protocol to detect the SARS-CoV-2 and other respiratory viruses in the hospital aerosol and monitor variations within the sequences. RT-PCR was also employed to estimate the viral load of SARS-CoV-2.

FINDINGS

13 of 15 (86.7%) samples exhibited SARS-CoV-2 with a relative abundance of 0.2-33.3%. The co-occurrence of human adenoviruses (type C1, C2, C5, C4), respiratory syncytial virus (RSV), influenza B, and non-SARS-CoV-229E were also recorded. Alignment of SARS-CoV-2 sequences against the reference strain of Wuhan China revealed variations in the form of single nucleotide polymorphisms (SNPs-17), insertions and deletions (indels-1). These variations were predicted to create missense (16), synonymous (15), frameshift (1) and stop-gained (1) mutations with a high (2), low (15), and moderate (16) impact.

CONCLUSIONS

Our results suggest that using CRVP on a large volume aerosol sample was a valuable tool for detecting SARS-CoV-2 in indoor aerosols of health care settings. Owing to its higher sensitivity, it can be employed as a surveillance strategy in the post COVID times to act as an early warning system to possibly control future outbreaks.

Targeted Sanger sequencing to recover key mutations in SARS-CoV-2 variant genome assemblies produced by next-generation sequencing

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is adaptively evolving to ensure its persistence within human hosts. It is therefore necessary to continuously monitor the emergence and prevalence of novel variants that arise. Importantly, some mutations have been associated with both molecular diagnostic failures and reduced or abrogated next-generation sequencing (NGS) read coverage in some genomic regions. Such impacts are particularly problematic when they occur in genomic regions such as those that encode the spike (S) protein, which are crucial for identifying and tracking the prevalence and dissemination dynamics of concerning viral variants. Targeted Sanger sequencing presents a fast and cost-effective means to accurately extend the coverage of whole-genome sequences. We designed a custom set of primers to amplify a 401 bp segment of the receptor-binding domain (RBD) (between positions 22698 and 23098 relative to the Wuhan-Hu-1 reference). We then designed a Sanger sequencing wet-laboratory protocol. We applied the primer set and wet-laboratory protocol to sequence 222 samples that were missing positions with key mutations K417N, E484K, and N501Y due to poor coverage after NGS sequencing. Finally, we developed SeqPatcher, a Python-based computational tool to analyse the trace files yielded by Sanger sequencing to generate consensus sequences, or take preanalysed consensus sequences in fasta format, and merge them with their corresponding whole-genome assemblies. We successfully sequenced 153 samples of 222 (69 %) using Sanger sequencing and confirmed the occurrence of key beta variant mutations (K417N, E484K, N501Y) in the S genes of 142 of 153 (93 %) samples. Additionally, one sample had the Y508F mutation and four samples the S477N. Samples with RT-PCR Ct scores ranging from 13.85 to 37.47 (mean=25.70) could be Sanger sequenced efficiently. These results show that our method and pipeline can be used to improve the quality of whole-genome assemblies produced using NGS and can be used with any pairs of the most used NGS and Sanger sequencing platforms.

Development and validation of a high throughput SARS-CoV-2 whole genome sequencing workflow in a clinical laboratory

Monitoring new mutations in SARS-CoV-2 provides crucial information for identifying diagnostic and therapeutic targets and important insights to achieve a more effective COVID-19 control strategy. Next generation sequencing (NGS) technologies have been widely used for whole genome sequencing (WGS) of SARS-CoV-2. While various NGS methods have been reported, one chief limitation has been the complexity of the workflow, limiting the scalability. Here, we overcome this limitation by designing a laboratory workflow optimized for high-throughput studies. The workflow utilizes modified ARTIC network v3 primers for SARS-CoV-2 whole genome amplification. NGS libraries were prepared by a 2-step PCR method, similar to a previously reported tailed PCR method, with further optimizations to improve amplicon balance, to minimize amplicon dropout for viral genomes harboring primer-binding site mutation(s), and to integrate robotic liquid handlers. Validation studies demonstrated that the optimized workflow can process up to 2688 samples in a single sequencing run without compromising sensitivity and accuracy and with fewer amplicon dropout events compared to the standard ARTIC protocol. We additionally report results for over 65,000 SARS-CoV-2 whole genome sequences from clinical specimens collected in the United States between January and September of 2021, as part of an ongoing national genomics surveillance effort.