Optimization of the SARS-CoV-2 ARTIC Network V4 Primers and Whole Genome Sequencing Protocol

Wastewater multiplex PCR amplicon sequencing revealed community transmission of SARS-CoV-2 lineages during the outbreak of infection in Chinese Mainland

During the COVID-19, wastewater-based epidemiology (WBE) has become a powerful epidemic surveillance tool widely used worldwide. However, the development and application of this technology in Chinese Mainland are relatively lagging. Herein, we for the first time monitored the community circulation of SARS-CoV-2 lineages using WBE methods in Chinese Mainland. During the peak period of infection outbreak at the end of 2022, six precious sewage samples were collected from the manhole in the student dormitory area on Wangjiang Campus of Sichuan University. RT-qPCR revealed that the six sewage samples were all positive for SARS-CoV-2 RNA. Multiplex PCR amplicon sequencing of the sewage samples reflected the local transmission of SARS-CoV-2 variants. The results of two deconvolution methods indicate that the main virus lineages have clear evolutionary genetic correlations. Furthermore, the sampling time is consistent with the timeline of concern for these virus lineages, as well as the timeline of uploading the nucleic acid sequences from the corresponding lineages in Sichuan to the database. These results demonstrate the reliability of the sewage sequencing results. Multiplex PCR amplicon sequencing is by far the most powerful analytical tool of WBE, enabling quantitative detection of virus lineages transmission and evolution at the community level.

Tracking the prevalence and emergence of SARS-CoV-2 variants of concern using a regional genomic surveillance program

SARS-CoV-2 molecular testing coupled with whole-genome sequencing is instrumental for real-time genomic surveillance. Genomic surveillance is critical for monitoring the spread of variants of concern (VOCs) as well as discovery of novel variants. Since the beginning of the pandemic, millions of SARS-CoV-2 genomes have been deposited into public sequence databases. This is the result of efforts of both national and regional diagnostic laboratories. In this study, we describe the results of SARS-CoV-2 genomic surveillance from February 2021 to June 2022 at a metropolitan hospital in the United States. We demonstrate that consistent daily sampling is sufficient to track the regional prevalence and emergence of VOCs and recapitulate national trends. Similar sampling efforts should be considered a viable option for local SARS-CoV-2 genomic surveillance at other regional laboratories.

IMPORTANCE

In our manuscript, we describe the results of SARS-CoV-2 genomic surveillance from February 2021 to June 2022 at a metropolitan hospital in the United States. We demonstrate that consistent daily sampling is sufficient to track the regional prevalence and emergence of variants of concern (VOCs). Similar sampling efforts should be considered a viable option for local SARS-CoV-2 genomic surveillance at other regional laboratories. While the SARS-CoV-2 pandemic has evolved into a more endemic form, we still believe that additional real-world information about sampling, procedures, and data interpretation is valuable for ongoing as well as future genomic surveillance efforts. Our study should be of substantial interest to clinical virologists.

Clinical Evaluation of the Accuracy of the Panbio™ COVID-19/Flu A&B Rapid Panel: A Combination Antigen Rapid Diagnostic Test for the Omicron Variant and Influenza A Virus

It is difficult to differentiate between coronavirus disease 2019 (COVID-19) and influenza based on the symptoms. In the present study, a newly developed antigen rapid diagnostic test (Ag-RDT) called Panbio™ COVID-19/Flu A&B that can simultaneously detect severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and influenza A/B virus was evaluated. Its accuracy was evaluated using 235 pairs of nasopharyngeal samples collected from patients with respiratory symptoms and fever (>37.5°C). Reverse transcription polymerase chain reaction was used as a reference method to evaluate the accuracy of the SARS-CoV-2 detection. We confirmed the accuracy of the developed Ag-RDT against the Omicron variant where the sensitivity and specificity were 94.8% and 100%, respectively. In addition, to identify the influenza A virus, a noninferiority test was conducted using a commercial Ag-RDT, which has a sensitivity and specificity in comparison with viral culture of 94.8% and 98.4%, respectively. The positive and negative predictive values for influenza A virus were 98.5% and 98.1%, respectively, for the Panbio COVID-19/Flu A&B test. The evaluation of this newly developed Ag-RDT using clinical samples suggests that it has a high efficacy in clinical settings.

Dengue Virus Surveillance in Nepal Yields the First On-Site Whole Genome Sequences of Isolates from the 2022 Outbreak

Background

The 4 serotypes of dengue virus (DENV1-4) can each cause potentially deadly dengue disease, and are spreading globally from tropical and subtropical areas to more temperate ones. Nepal provides a microcosm of this global phenomenon, having met each of these grim benchmarks. To better understand DENV transmission dynamics and spread into new areas, we chose to study dengue in Nepal and, in so doing, to build the onsite infrastructure needed to manage future, larger studies.

Methods and Results

During the 2022 dengue season, we enrolled 384 patients presenting at a hospital in Kathmandu with dengue-like symptoms; 79% of the study participants had active or recent DENV infection (NS1 antigen and IgM). To identify circulating serotypes, we screened serum from 50 of the NS1 + participants by RT-PCR and identified DENV1, 2, and 3 - with DENV1 and 3 codominant. We also performed whole-genome sequencing of DENV, for the first time in Nepal, using our new on-site capacity. Sequencing analysis demonstrated the DENV1 and 3 genomes clustered with sequences reported from India in 2019, and the DENV2 genome clustered with a sequence reported from China in 2018.

Conclusion

These findings highlight DENV's geographic expansion from neighboring countries, identify China and India as the likely origin of the 2022 DENV cases in Nepal, and demonstrate the feasibility of building onsite capacity for more rapid genomic surveillance of circulating DENV. These ongoing efforts promise to protect populations in Nepal and beyond by informing the development and deployment of DENV drugs and vaccines in real time.

A solution to achieve sequencing from SARS-CoV-2 specimens with low viral loads: concatenation of reads from independent reactions

BACKGROUND

During the pandemic, whole genome sequencing was critical to characterize SARS-CoV-2 for surveillance, clinical and therapeutical purposes. However, low viral loads in specimens often led to suboptimal sequencing, making lineage assignment and phylogenetic analysis difficult. We propose an alternative approach to sequencing these specimens that involves sequencing in triplicate and concatenation of the reads obtained using bioinformatics. This proposal is based on the hypothesis that the uncovered regions in each replicate differ and that concatenation would compensate for these gaps and recover a larger percentage of the sequenced genome.

RESULTS

Whole genome sequencing was performed in triplicate on 30 samples with Ct > 32 and the benefit of replicate read concatenation was assessed. After concatenation: i) 28% of samples reached the standard quality coverage threshold (> 90% genome covered > 30x); ii) 39% of samples did not reach the coverage quality thresholds but coverage improved by more than 40%; and iii) SARS-CoV-2 lineage assignment was possible in 68.7% of samples where it had been impaired.

CONCLUSIONS

Concatenation of reads from replicate sequencing reactions provides a simple way to access hidden information in the large proportion of SARS-CoV-2-positive specimens eliminated from analysis in standard sequencing schemes. This approach will enhance our potential to rule out involvement in outbreaks, to characterize reinfections and to identify lineages of concern for surveillance or therapeutical purposes.

Genomic surveillance during the first two years of the COVID-19 pandemic - country experience and lessons learned from Türkiye

Background

Türkiye confirmed its first case of SARS-CoV-2 on March 11, 2020, coinciding with the declaration of the global COVID-19 pandemic. Subsequently, Türkiye swiftly increased testing capacity and implemented genomic sequencing in 2020. This paper describes Türkiye's journey of establishing genomic surveillance as a middle-income country with limited prior sequencing capacity and analyses sequencing data from the first two years of the pandemic. We highlight the achievements and challenges experienced and distill globally relevant lessons.

Methods

We tracked the evolution of the COVID-19 pandemic in Türkiye from December 2020 to February 2022 through a timeline and analysed epidemiological, vaccination, and testing data. To investigate the phylodynamic and phylogeographic aspects of SARS-CoV-2, we used Nextstrain to analyze 31,629 high-quality genomes sampled from seven regions nationwide.

Results

Türkiye's epidemiological curve, mirroring global trends, featured four distinct waves, each coinciding with the emergence and spread of variants of concern (VOCs). Utilizing locally manufactured kits to expand testing capacity and introducing variant-specific quantitative reverse transcription polymerase chain reaction (RT-qPCR) tests developed in partnership with a private company was a strategic advantage in Türkiye, given the scarcity and fragmented global supply chain early in the pandemic. Türkiye contributed more than 86,000 genomic sequences to global databases by February 2022, ensuring that Turkish data was reflected globally. The synergy of variant-specific RT-qPCR kits and genomic sequencing enabled cost-effective monitoring of VOCs. However, data analysis was constrained by a weak sequencing sampling strategy and fragmented data management systems, limiting the application of sequencing data to guide the public health response. Phylodynamic analysis indicated that Türkiye's geographical position as an international travel hub influenced both national and global transmission of each VOC despite travel restrictions.

Conclusion

This paper provides valuable insights into the testing and genomic surveillance systems adopted by Türkiye during the COVID-19 pandemic, proposing important lessons for countries developing national systems. The findings underscore the need for robust testing and sampling strategies, streamlined sample referral, and integrated data management with metadata linkage and data quality crucial for impactful epidemiological analysis. We recommend developing national genomic surveillance strategies to guide sustainable and integrated expansion of capacities built for COVID-19 and to optimize the effective utilization of sequencing data for public health action.

Considerations and Opportunities for Probe Capture Enrichment Sequencing of Emerging Viruses from Wastewater

Until recently, wastewater-based monitoring for pathogens of public health concern primarily used PCR-based quantification methods and targeted sequencing for specific pathogens (e.g., SARS-CoV-2). In the past three years, researchers have expanded sequencing to monitor a broad range of pathogens, applying probe capture enrichment to wastewater. The goals of those studies included (1) monitoring and expanding fundamental knowledge of disease dynamics for known pathogens and (2) evaluating the potential for early detection of emerging diseases resulting from zoonotic spillover or novel viral variants. Several studies using off-the-shelf probe panels designed for clinical and environmental surveillance reported that enrichment increased virus relative abundance but did not recover complete genomes for most nonenteric viruses. Based on our experience and recent results reported by others using these panels for wastewater, clinical, and synthetic samples, we discuss challenges and technical factors that affect the rates of false positive and false negative results. We identify trade-offs and opportunities throughout the workflow, including in wastewater sample processing, probe panel design, and bioinformatic analysis. We suggest tailored methods of virus concentration and background removal, carefully designed probe panels, and multithresholded bioinformatics analysis.

SARS-CoV-2 shedding, infectivity, and evolution in an immunocompromised adult patient

This study aimed to provide further insight into the evolutionary dynamics of SARS-CoV-2 by analyzing the case of a 40-year-old man who had previously undergone autologous hematopoietic stem cell transplantation due to a diffuse large B-cell lymphoma. He developed a persistent SARS-CoV-2 infection lasting at least 218 days and did not manifest a humoral immune response to the virus during this follow-up period. Whole-genome sequencing and viral cultures confirmed a persistent infection with a replication-positive virus that had undergone genetic variation for at least 196 days after symptom onset.

Wastewater sequencing as a powerful tool to reveal SARS-CoV-2 variant introduction and spread in French Guiana, South America

The origin of introduction of a new pathogen in a country, the evolutionary dynamics of an epidemic within a country, and the role of cross-border areas on pathogen dynamics remain complex to disentangle and are often poorly understood. For instance, cross-border areas represent the ideal location for the sharing of viral variants between countries, with international air travel, land travel and waterways playing an important role in the cross-border spread of infectious diseases. Unfortunately, monitoring the point of entry and the evolutionary dynamics of viruses in space and time within local populations remain challenging. Here we tested the efficiency of wastewater-based epidemiology and genotyping in monitoring Covid-19 epidemiology and SARS-CoV-2 variant dynamics in French Guiana, a tropical country located in South America. Our results suggest that wastewater-based epidemiology and genotyping are powerful tools to monitor variant introduction and disease evolution within a tropical country but the inclusion of both clinical and wastewater samples could still improve our understanding of genetic diversity co-circulating. Wastewater sequencing also revealed the cryptic transmission of SARS-CoV-2 variants within the country. Interestingly, we found some amino acid changes specific to the variants co-circulating in French Guiana, suggesting a local evolution of the SARS-CoV-2 variants after their introduction. More importantly, our results showed that the proximity to bordering countries was not the origin of the emergence of the French Guianese B.1.160.25 variant, but rather that this variant emerged from an ancestor B.1.160 variant introduced by European air plane travelers, suggesting thus that air travel remains a significant risk for cross-border spread of infectious diseases. Overall, we suggest that wastewater-based epidemiology and genotyping provides a cost effective and non-invasive approach for pathogen monitoring and an early-warning tool for disease emergence and spread within a tropical country.

Genomic surveillance of SARS-CoV-2 using long-range PCR primers

Introduction

Whole Genome Sequencing (WGS) of the SARS-CoV-2 virus is crucial in the surveillance of the COVID-19 pandemic. Several primer schemes have been developed to sequence nearly all of the ~30,000 nucleotide SARS-CoV-2 genome, using a multiplex PCR approach to amplify cDNA copies of the viral genomic RNA. Midnight primers and ARTIC V4.1 primers are the most popular primer schemes that can amplify segments of SARS-CoV-2 (400 bp and 1200 bp, respectively) tiled across the viral RNA genome. Mutations within primer binding sites and primer-primer interactions can result in amplicon dropouts and coverage bias, yielding low-quality genomes with 'Ns' inserted in the missing amplicon regions, causing inaccurate lineage assignments, and making it challenging to monitor lineage-specific mutations in Variants of Concern (VoCs).

Methods

In this study we used a set of seven long-range PCR primer pairs to sequence clinical isolates of SARS-CoV-2 on Oxford Nanopore sequencer. These long-range primers generate seven amplicons approximately 4500 bp that covered whole genome of SARS-CoV-2. One of these regions includes the full-length S-gene by using a set of flanking primers. We also evaluated the performance of these long-range primers with Midnight primers by sequencing 94 clinical isolates in a Nanopore flow cell.

Results and discussion

Using a small set of long-range primers to sequence SARS-CoV-2 genomes reduces the possibility of amplicon dropout and coverage bias. The key finding of this study is that long range primers can be used in single-molecule sequencing of RNA viruses in surveillance of emerging variants. We also show that by designing primers flanking the S-gene, we can obtain reliable identification of SARS-CoV-2 variants.

Targeted amplification-based whole genome sequencing of Monkeypox virus in clinical specimens

IMPORTANCE

We present a protocol to efficiently sequence genomes of the MPXV-causing mpox. This enables researchers and public health agencies to acquire high-quality genomic data using a rapid and cost-effective approach. Genomic data can be used to conduct surveillance and investigate mpox outbreaks. We present 91 mpox genomes that show the diversity of the 2022 mpox outbreak in Ontario, Canada.

NASCarD (Nanopore Adaptive Sampling with Carrier DNA): A Rapid, PCR-Free Method for SARS-CoV-2 Whole-Genome Sequencing in Clinical Samples

Whole-genome sequencing (WGS) represents the main technology for SARS-CoV-2 lineage characterization in diagnostic laboratories worldwide. The rapid, near-full-length sequencing of the viral genome is commonly enabled by high-throughput sequencing of PCR amplicons derived from cDNA molecules. Here, we present a new approach called NASCarD (Nanopore Adaptive Sampling with Carrier DNA), which allows a low amount of nucleic acids to be sequenced while selectively enriching for sequences of interest, hence limiting the production of non-target sequences. Using COVID-19 positive samples available during the omicron wave, we demonstrate how the method may lead to >99% genome completeness of the SARS-CoV-2 genome sequences within 7 h of sequencing at a competitive cost. The new approach may have applications beyond SARS-CoV-2 sequencing for other DNA or RNA pathogens in clinical samples.

Genome stability assessment of PRRS vaccine strain with new ARTIC-style sequencing protocol

A tiling amplicon sequencing protocol was developed to analyse the genome sequence stability of the modified live PRRSV vaccine strain, Porcilis MLV. The backbone of the ARTIC-style protocol was formed by 34 individual primer pairs, which were divided into two primer pools. Primer pairs were designed to amplify 532 to 588 bp fragments of the corresponding genomic region. The amplicons are suitable for sequencing on Illumina DNA sequencers with available 600-cycle sequencing kits. The concentration of primer pairs in the pools was optimized to obtain a balanced sequencing depth along the genome. Deep sequencing data of three vaccine batches were also analysed. All three vaccine batches were very similar to each other, although they also showed single nucleotide variations (SNVs) affecting less than 1 % of the genome. In the three vaccine strains, 113 to 122 SNV sites were identified; at these sites, the minority variants represented a frequency range of 1 to 48.7 percent. Additionally, the strains within the batches contained well-known length polymorphisms; the genomes of these minority deletion mutants were 135 to 222 bp shorter than the variant with the complete genome. Our results show the usefulness of ARTIC-style protocols in the evaluation of the genomic stability of PRRS MLV strains.

A rapid and affordable amplicon-based method for next-generation genome sequencing of the infectious bursal disease virus

The infectious bursal disease virus (IBDV) causes a severe immunosuppressive disorder in young chickens. IBDV evolution resulted in the emergence of strains with divergent genetic, antigenic, and pathogenic characteristics. Genetic classification is typically performed by sequencing the coding region of the most immunogenic region of the viral protein 2 (VP2). Sequencing both double-stranded RNA genome segments is essential to achieve a more comprehensive IBDV classification that can detect recombinants and reassortments. Here, we report the development and standardization of a tiled PCR amplicon protocol for the direct and cost-effective genome sequencing of global IBDV strains using next-generation technology. Primers for tiled PCR were designed with adapters to bypass expensive and time-consuming library preparation steps. Sequencing was performed on Illumina MiniSeq equipment, and fourteen complete genomes of field strains were assembled using reference sequences. The PCR-enrichment step was used to obtain genomes from low-titer biological samples that were difficult to amplify using traditional sequencing. Phylogenetic analyses of the obtained genomes confirmed previous strain classification. By combining the enrichment methodology with massive sequencing, it is possible to obtain IBDV genomic sequences in a fast and affordable manner. This procedure can be a valuable tool to better understand virus epidemiology.

SARS-CoV-2 genomics and impact on clinical care for COVID-19

The emergence and worldwide spread of SARS-CoV-2 during the COVID-19 pandemic necessitated the adaptation and rapid deployment of viral WGS and analysis techniques that had been previously applied on a more limited basis to other viral pathogens, such as HIV and influenza viruses. The need for WGS was driven in part by the low mutation rate of SARS-CoV-2, which necessitated measuring variation along the entire genome sequence to effectively differentiate lineages and characterize viral evolution. Several WGS approaches designed to maximize throughput and accuracy were quickly adopted by surveillance labs around the world. These broad-based SARS-CoV-2 genomic sequencing efforts revealed ongoing evolution of the virus, highlighted by the successive emergence of new viral variants throughout the course of the pandemic. These genomic insights were instrumental in characterizing the effects of viral mutations on transmissibility, immune escape and viral tropism, which in turn helped guide public health policy, the use of monoclonal antibody therapeutics and vaccine development strategies. As the use of direct-acting antivirals for the treatment of COVID-19 became more widespread, the potential for emergence of antiviral resistance has driven ongoing efforts to delineate resistance mutations and to monitor global sequence databases for their emergence. Given the critical role of viral genomics in the international effort to combat the COVID-19 pandemic, coordinated efforts should be made to expand global genomic surveillance capacity and infrastructure towards the anticipation and prevention of future pandemics.

A multiplexed RT-PCR assay for nanopore whole genome sequencing of Tilapia lake virus (TiLV)

Tilapia lake virus (TiLV) is a highly contagious viral pathogen that affects tilapia, a globally significant and affordable source of fish protein. To prevent the introduction and spread of TiLV and its impact, there is an urgent need for increased surveillance, improved biosecurity measures, and continuous development of effective diagnostic and rapid sequencing methods. In this study, we have developed a multiplexed RT-PCR assay that can amplify all ten complete genomic segments of TiLV from various sources of isolation. The amplicons generated using this approach were immediately subjected to real-time sequencing on the Nanopore system. By using this approach, we have recovered and assembled 10 TiLV genomes from total RNA extracted from naturally TiLV-infected tilapia fish, concentrated tilapia rearing water, and cell culture. Our phylogenetic analysis, consisting of more than 36 TiLV genomes from both newly sequenced and publicly available TiLV genomes, provides new insights into the high genetic diversity of TiLV. This work is an essential steppingstone towards integrating rapid and real-time Nanopore-based amplicon sequencing into routine genomic surveillance of TiLV, as well as future vaccine development.

Bioinformatic investigation of discordant sequence data for SARS-CoV-2: insights for robust genomic analysis during pandemic surveillance

The COVID-19 pandemic has necessitated the rapid development and implementation of whole-genome sequencing (WGS) and bioinformatic methods for managing the pandemic. However, variability in methods and capabilities between laboratories has posed challenges in ensuring data accuracy. A national working group comprising 18 laboratory scientists and bioinformaticians from Australia and New Zealand was formed to improve data concordance across public health laboratories (PHLs). One effort, presented in this study, sought to understand the impact of the methodology on consensus genome concordance and interpretation. SARS-CoV-2 WGS proficiency testing programme (PTP) data were retrospectively obtained from the 2021 Royal College of Pathologists of Australasia Quality Assurance Programmes (RCPAQAP), which included 11 participating Australian laboratories. The submitted consensus genomes and reads from eight contrived specimens were investigated, focusing on discordant sequence data and findings were presented to the working group to inform best practices. Despite using a variety of laboratory and bioinformatic methods for SARS-CoV-2 WGS, participants largely produced concordant genomes. Two participants returned five discordant sites in a high-Cτ replicate, which could be resolved with reasonable bioinformatic quality thresholds. We noted ten discrepancies in genome assessment that arose from nucleotide heterogeneity at three different sites in three cell-culture-derived control specimens. While these sites were ultimately accurate after considering the participants' bioinformatic parameters, it presented an interesting challenge for developing standards to account for intrahost single nucleotide variation (iSNV). Observed differences had little to no impact on key surveillance metrics, lineage assignment and phylogenetic clustering, while genome coverage <90 % affected both. We recommend PHLs bioinformatically generate two consensus genomes with and without ambiguity thresholds for quality control and downstream analysis, respectively, and adhere to a minimum 90 % genome coverage threshold for inclusion in surveillance interpretations. We also suggest additional PTP assessment criteria, including primer efficiency, detection of iSNVs and minimum genome coverage of 90 %. This study underscores the importance of multidisciplinary national working groups in informing guidelines in real time for bioinformatic quality acceptance criteria. It demonstrates the potential for enhancing public health responses through improved data concordance and quality control in SARS-CoV-2 genomic analysis during pandemic surveillance.

"Nano COVID-19": Nanopore sequencing of spike gene to identify SARS-CoV-2 variants of concern

Coronavirus disease 2019 (COVID-19) is a worldwide pandemic infectious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). World Health Organization (WHO) has defined the viral variants of concern (VOC) which cause more severe disease, higher transmissibility, and reduced vaccine efficacy. In this study, the "Nano COVID-19" workflow based on Oxford nanopore sequencing of the full-length spike gene combined with flexible data analysis options was developed to identify SARS-CoV-2 VOCs. The primers were designed to cover the full-length spike gene and can amplify all VOC strains. The results of VOC identification based on phylogenetic analysis of the full-length spike gene were comparable to the whole genome sequencing (WGS). Compared to the standard VOC identification pipeline, the fast analysis based on Read Assignment, Mapping, and Phylogenetic Analysis in Real Time (RAMPART) and the user-friendly method based on EPI2ME yielded 89.3% and 97.3% accuracy, respectively. The EPI2ME pipeline is recommended for researchers without bioinformatic skills, whereas RAMPART is more suitable for bioinformaticians. This workflow provides a cost-effective, simplified pipeline with a rapid turnaround time. Furthermore, it is portable to point-of-care SARS-CoV-2 VOC identification and compatible with large-scale analysis. Therefore, "Nano COVID-19" is an alternative viral epidemic screening and transmission tracking workflow.

Nanopore sequencing technology and its applications

Since the development of Sanger sequencing in 1977, sequencing technology has played a pivotal role in molecular biology research by enabling the interpretation of biological genetic codes. Today, nanopore sequencing is one of the leading third-generation sequencing technologies. With its long reads, portability, and low cost, nanopore sequencing is widely used in various scientific fields including epidemic prevention and control, disease diagnosis, and animal and plant breeding. Despite initial concerns about high error rates, continuous innovation in sequencing platforms and algorithm analysis technology has effectively addressed its accuracy. During the coronavirus disease (COVID-19) pandemic, nanopore sequencing played a critical role in detecting the severe acute respiratory syndrome coronavirus-2 virus genome and containing the pandemic. However, a lack of understanding of this technology may limit its popularization and application. Nanopore sequencing is poised to become the mainstream choice for preventing and controlling COVID-19 and future epidemics while creating value in other fields such as oncology and botany. This work introduces the contributions of nanopore sequencing during the COVID-19 pandemic to promote public understanding and its use in emerging outbreaks worldwide. We discuss its application in microbial detection, cancer genomes, and plant genomes and summarize strategies to improve its accuracy.

Genomic Surveillance of SARS-CoV-2 Using Long-Range PCR Primers

Whole Genome Sequencing (WGS) of the SARS-CoV-2 virus is crucial in the surveillance of the COVID-19 pandemic. Several primer schemes have been developed to sequence the ~30,000 nucleotide SARS-CoV-2 genome that use a multiplex PCR approach to amplify cDNA copies of the viral genomic RNA. Midnight primers and ARTIC V4.1 primers are the most popular primer schemes that can amplify segments of SARS-CoV-2 (400 bp and 1200 bp, respectively) tiled across the viral RNA genome. Mutations within primer binding sites and primer-primer interactions can result in amplicon dropouts and coverage bias, yielding low-quality genomes with 'Ns' inserted in the missing amplicon regions, causing inaccurate lineage assignments, and making it challenging to monitor lineage-specific mutations in Variants of Concern (VoCs). This study uses seven long-range PCR primers with an amplicon size of ~4500 bp to tile across the complete SARS-CoV-2 genome. One of these regions includes the full-length S-gene by using a set of flanking primers. Using a small set of long-range primers to sequence SARS-CoV-2 genomes reduces the possibility of amplicon dropout and coverage bias.

Toward a global virus genomic surveillance network

The COVID-19 pandemic galvanized the field of virus genomic surveillance, demonstrating its utility for public health. Now, we must harness the momentum that led to increased infrastructure, training, and political will to build a sustainable global genomic surveillance network for other epidemic and endemic viruses. We suggest a generalizable modular sequencing framework wherein users can easily switch between virus targets to maximize cost-effectiveness and maintain readiness for new threats. We also highlight challenges associated with genomic surveillance and when global inequalities persist. We propose solutions to mitigate some of these issues, including training and multilateral partnerships. Exploring alternatives to clinical sequencing can also reduce the cost of surveillance programs. Finally, we discuss how establishing genomic surveillance would aid control programs and potentially provide a warning system for outbreaks, using a global respiratory virus (RSV), an arbovirus (dengue virus), and a regional zoonotic virus (Lassa virus) as examples.

Novel Amplicon-Based Sequencing Approach to West Nile Virus

West Nile virus is a re-emerging arbovirus whose impact on public health is increasingly important as more and more epidemics and epizootics occur, particularly in America and Europe, with evidence of active circulation in Africa. Because birds constitute the main reservoirs, migratory movements allow the diffusion of various lineages in the world. It is therefore crucial to properly control the dispersion of these lineages, especially because some have a greater health impact on public health than others. This work describes the development and validation of a novel whole-genome amplicon-based sequencing approach to West Nile virus. This study was carried out on different strains from lineage 1 and 2 from Senegal and Italy. The presented protocol/approach showed good coverage using samples derived from several vertebrate hosts and may be valuable for West Nile genomic surveillance.

Tracking the prevalence and emergence of SARS CoV2 variants of concern using a regional genomic surveillance program

SARS-CoV-2 molecular testing coupled with whole genome sequencing is instrumental for real-time genomic surveillance. Genomic surveillance is critical for monitoring the spread of variants of concern (VOC) as well as novel variant discovery. Since the beginning of the pandemic millions of SARS-CoV-2 genomes have been deposited into public sequence databases. This is the result of efforts of both national and regional diagnostic laboratories. Here we describe the results of SARS-CoV-2 genomic surveillance from February 2021 to June 2022 at a metropolitan hospital in the USA. We demonstrate that consistent daily sampling is sufficient to track the regional prevalence and emergence of VOC. Similar sampling efforts should be considered a viable option for local SARS-CoV-2 genomic surveillance at other regional laboratories.

Comparison of SARS-CoV-2 whole genome sequencing using tiled amplicon enrichment and bait hybridization

The severe acute respiratory syndrome coronavirus 2 (SARS‑CoV‑2) pandemic has led to extensive virological monitoring by whole genome sequencing (WGS). Investigating the advantages and limitations of different protocols is key when conducting population-level WGS. SARS-CoV-2 positive samples with Ct values of 14-30 were run using three different protocols: the Twist Bioscience SARS‑CoV‑2 protocol with bait hybridization enrichment sequenced with Illumina, and two tiled amplicon enrichment protocols, ARTIC V3 and Midnight, sequenced with Illumina and Oxford Nanopore Technologies, respectively. Twist resulted in better coverage uniformity and coverage of the entire genome, but has several drawbacks: high human contamination, laborious workflow, high cost, and variation between batches. The ARTIC and Midnight protocol produced an even coverage across samples, and almost all reads were mapped to the SARS-CoV-2 reference. ARTIC and Midnight represent robust, cost-effective, and highly scalable methods that are appropriate in a clinical environment. Lineage designations were uniform across methods, representing the dominant lineages in Sweden during the period of collection. This study provides insights into methodological differences in SARS‑CoV‑2 sequencing and guidance in selecting suitable methods for various purposes.

A Multiplexed, Tiled PCR Method for Rapid Whole-Genome Sequencing of Infectious Spleen and Kidney Necrosis Virus (ISKNV) in Tilapia

Tilapia farming is one of the most important sectors in aquaculture worldwide and of major importance to global food security. Infectious spleen and kidney necrosis virus (ISKNV) has been identified as an agent of high morbidity and mortality, threatening tilapia aquaculture. ISKNV was detected in Lake Volta, Ghana, in September 2018 and spread rapidly, with mortality rates between 60 and 90% and losses of more than 10 tonnes of fish per day. Understanding the spread and evolution of viral pathogens is important for control strategies. Here, we developed a tiled-PCR sequencing approach for the whole-genome sequencing of ISKNV, using long read sequencing to enable field-based, real-time genomic surveillance. This work represents the first use of tiled-PCR for whole genome recovery of viruses in aquaculture, with the longest genome target (>110 kb dsDNA) to date. Our protocol was applied to field samples collected from the ISKNV outbreaks from four intensive tilapia cage culture systems across Lake Volta, between October 2018 and May 2022. Despite the low mutation rate of dsDNA viruses, 20 single nucleotide polymorphisms accumulated during the sampling period. Droplet digital PCR identified a minimum requirement of template in a sample to recover 50% of an ISKNV genome at 275 femtograms (2410 viral templates per 5 µL sequencing reaction). Overall, tiled-PCR sequencing of ISKNV provides an informative tool to assist in disease control in aquaculture.

Robust Performance of SARS-CoV-2 Whole-Genome Sequencing from Wastewater with a Nonselective Virus Concentration Method

The sequencing of human virus genomes from wastewater samples is an efficient method for tracking viral transmission and evolution at the community level. However, this requires the recovery of viral nucleic acids of high quality. We developed a reusable tangential-flow filtration system to concentrate and purify viruses from wastewater for genome sequencing. A pilot study was conducted with 94 wastewater samples from four local sewersheds, from which viral nucleic acids were extracted, and the whole genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was sequenced using the ARTIC V4.0 primers. Our method yielded a high probability (0.9) of recovering complete or near-complete SARS-CoV-2 genomes (>90% coverage at 10× depth) from wastewater when the COVID-19 incidence rate exceeded 33 cases per 100 000 people. The relative abundances of sequenced SARS-CoV-2 variants followed the trends observed from patient-derived samples. We also identified SARS-CoV-2 lineages in wastewater that were underrepresented or not present in the clinical whole-genome sequencing data. The developed tangential-flow filtration system can be easily adopted for the sequencing of other viruses in wastewater, particularly those at low concentrations.

Maximum likelihood pandemic-scale phylogenetics

Phylogenetics has a crucial role in genomic epidemiology. Enabled by unparalleled volumes of genome sequence data generated to study and help contain the COVID-19 pandemic, phylogenetic analyses of SARS-CoV-2 genomes have shed light on the virus's origins, spread, and the emergence and reproductive success of new variants. However, most phylogenetic approaches, including maximum likelihood and Bayesian methods, cannot scale to the size of the datasets from the current pandemic. We present 'MAximum Parsimonious Likelihood Estimation' (MAPLE), an approach for likelihood-based phylogenetic analysis of epidemiological genomic datasets at unprecedented scales. MAPLE infers SARS-CoV-2 phylogenies more accurately than existing maximum likelihood approaches while running up to thousands of times faster, and requiring at least 100 times less memory on large datasets. This extends the reach of genomic epidemiology, allowing the continued use of accurate phylogenetic, phylogeographic and phylodynamic analyses on datasets of millions of genomes.

High-throughput sequencing approaches applied to SARS-CoV-2

High-throughput sequencing is crucial for surveillance and control of viral outbreaks. During the ongoing coronavirus disease 2019 (COVID-19) pandemic, advances in the high-throughput sequencing technology resources have enhanced diagnosis, surveillance, and vaccine discovery. From the onset of the pandemic in December 2019, several genome-sequencing approaches have been developed and supported across the major sequencing platforms such as Illumina, Oxford Nanopore, PacBio, MGI DNBSEQTM and Ion Torrent. Here, we share insights from the sequencing approaches developed for sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) between December 2019 and October 2022.

Mudhune G, M. Morobe J, Said Mohammed K, O. Makori T, Ndwiga L, W. Mburu M, O. Moraa E, Musyoki J, Murunga N, N. Waliaula I, K. Mumelo A, Bejon P, Isabella Ochola-Oyier L, Githinji G, Nokes J, Agoti C. Temporal distribution and clinical characteristics of the Alpha, Delta and Omicron SARS-CoV-2 variants of concern in Laikipia, Kenya: institutional and community-based genomic surveillance

Background: Understanding the molecular epidemiology and clinical presentation of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOC) in rural-urban populations in Kenya is important for informing future public health responses and clinical care. Methods: We undertook a retrospective analysis of the clinical presentation and phylogenetic relatedness of specimens from 97 SARS-CoV-2 cases collected between 24th April and 31st December 2021 in Laikipia county, Kenya. VOC were related to observed symptoms. Phylogenetic analyses included contemporaneous sequences from across Kenya and the globe, to contextualise local transmission dynamics. Results: These sequences fell into three VOC; Alpha (n=8), Delta (n=52) and Omicron (n=37). We estimated 75 independent SARS-CoV-2 introductions into the county. The Alpha and Delta VOC were commonly detected in persons aged 31 to 45 years, 50.0% and 30.8%, respectively. The Omicron VOC was mostly detected in 16 to 30-year-olds (51.4%). Whereas relative to the other VOCs, Omicron was associated with mild upper-respiratory tract symptoms (cough, OR 3.78; 95% CI 1.1 – 16.74, p= 0.026) and sore throat, OR 22.42; 95% CI 7.11 – 81.40, p<0.001), Delta was associated with moderate to severe lower-respiratory tract symptoms (shortness of breath, OR 26.8; 95% CI 3.89 – 1158.14, p<0.001) and fever (OR 6.11; 95% CI 1.57 – 35.35, p= 0.004). Post-acute phase neurological complications were suspected in four Delta infected cases (neuralgia, neuritis, peripheral neuropathy, numbness of hand and tinnitus). Conclusion: We highlight the distinctive clinical characteristics of SARS-CoV-2 VOCs, as observed in Laikipia, Kenya, to support evidence-based clinical decisions. Multiple introductions of the VOCs were recorded despite the public health measures that were in place questioning their effectiveness during the study period.

Accurate Detection of SARS-CoV-2 by Next-Generation Sequencing in Low Viral Load Specimens

As new SARS-CoV-2 variants emerge, there is an urgent need to increase the efficiency and availability of viral genome sequencing, notably to detect the lineage in samples with a low viral load. SARS-CoV-2 genome next-generation sequencing (NGS) was performed retrospectively in a single center on 175 positive samples from individuals. An automated workflow used the Ion AmpliSeq SARS-CoV-2 Insight Research Assay on the Genexus Sequencer. All samples were collected in the metropolitan area of the city of Nice (France) over a period of 32 weeks (from 19 July 2021 to 11 February 2022). In total, 76% of cases were identified with a low viral load (Ct ≥ 32, and ≤200 copies/µL). The NGS analysis was successful in 91% of cases, among which 57% of cases harbored the Delta variant, and 34% the Omicron BA.1.1 variant. Only 9% of cases had unreadable sequences. There was no significant difference in the viral load in patients infected with the Omicron variant compared to the Delta variant (Ct values, p = 0.0507; copy number, p = 0.252). We show that the NGS analysis of the SARS-CoV-2 genome provides reliable detection of the Delta and Omicron SARS-CoV-2 variants 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.

Rapid System to Detect Variants of SARS-CoV-2 in Nasopharyngeal Swabs

Currently, the reference method for identifying the presence of variants of SARS-CoV-2 is whole genome sequencing. Although it is less expensive than in the past, it is still time-consuming, and interpreting the results is difficult, requiring staff with specific skills who are not always available in diagnostic laboratories. The test presented in this study aimed to detect, using traditional real-time PCR, the presence of the main variants described for the spike protein of the SARS-CoV-2 genome. The primers and probes were designed to detect the main deletions that characterize the different variants. The amplification targets were deletions in the S gene: 25-27, 69-70, 241-243, and 157-158. In the ORF1a gene, the deletion 3675-3677 was chosen. Some of these mutations can be considered specific variants, while others can be identified by the simultaneous presence of one or more deletions. We avoided using point mutations in order to improve the speed of the test. Our test can help clinical and medical microbiologists quickly recognize the presence of variants in biological samples (particularly nasopharyngeal swabs). The test can also be used to identify variants of the virus that could potentially be more diffusive as well as not responsive to the vaccine.

Nyiro J, M. Morobe J, Makori TN, Mburu MW, Ndwiga L, Moraa EO, Musyoki J, Murunga N, Bejon P, Ochola-Oyier LI, Nokes DJ, Agoti CN, Githinji G. Detection of a SARS-CoV-2 Beta-like variant in coastal Kenya after more than a year of disappearance

Continuous genomic surveillance is necessary and important to inform emergence of new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants and resurgence of previous circulating variants. In this study, we provide genomic characterisation of a Beta variant sequence identified through out genomic surveillance platform more than a year since the last reported case of Beta. This variant contained additional mutations associated with immune escape that have been observed in other newer variants such as Delta and Omicron implying ongoing convergent evolution of this variant in the community or in an immunocompromised patient.

Identification of Genomic Variants of SARS-CoV-2 Using Nanopore Sequencing

Background and Objectives: SARS-CoV-2 is the first global threat and life-changing event of the twenty-first century. Although efficient treatments and vaccines have been developed, due to the virus's ability to mutate in key regions of the genome, whole viral genome sequencing is needed for efficient monitoring, evaluation of the spread, and even the adjustment of the molecular diagnostic assays. Materials and Methods: In this study, Nanopore and Ion Torrent sequencing technologies were used to detect the main SARS-CoV-2 circulating strains in Timis County, Romania, between February 2021 and May 2022. Results: We identified 22 virus lineages belonging to seven clades: 20A, 20I (Alpha, V1), 21B (Kappa), 21I (Delta), 21J (Delta), 21K (Omicron), and 21L (Omicron). Conclusions: Results obtained with both methods are comparable, and we confirm the utility of Nanopore sequencing in large-scale epidemiological surveillance due to the lower cost and reduced time for library preparation.

Genomic evolution of SARS-CoV-2 in Reunion Island

Island communities are interesting study sites for microbial evolution during epidemics, as their insular nature reduces the complexity of the population's connectivity. This was particularly true on Reunion Island during the first half of 2021, when international travel was restricted in order to mitigate the risk for SARS-CoV-2 introductions. Concurrently, the SARS-CoV-2 Beta variant became dominant and started to circulate at high levels for several months before being completely replaced by the Delta variant as of October 2021. Here, we explore some of the particularities of SARS-CoV-2 genomic evolution within the insular context of Reunion Island. We show that island isolation allowed the amplification and expansion of unique genetic lineages that remained uncommon across the globe. Islands are therefore potential hotspots for the emergence of new genetic variants, meaning that they will play a key role in the continued evolution and propagation of COVID-19 as the pandemic persists.

Next Generation Sequencing Approaches to Characterize the Respiratory Tract Virome

The COVID-19 pandemic and heightened perception of the risk of emerging viral infections have boosted the efforts to better understand the virome or complete repertoire of viruses in health and disease, with a focus on infectious respiratory diseases. Next-generation sequencing (NGS) is widely used to study microorganisms, allowing the elucidation of bacteria and viruses inhabiting different body systems and identifying new pathogens. However, NGS studies suffer from a lack of standardization, in particular, due to various methodological approaches and no single format for processing the results. Here, we review the main methodological approaches and key stages for studies of the human virome, with an emphasis on virome changes during acute respiratory viral infection, with applications for clinical diagnostics and epidemiologic analyses.

Challenges and Opportunities for Global Genomic Surveillance Strategies in the COVID-19 Era

Global SARS-CoV-2 genomic surveillance efforts have provided critical data on the ongoing evolution of the virus to inform best practices in clinical care and public health throughout the pandemic. Impactful genomic surveillance strategies generally follow a multi-disciplinary pipeline involving clinical sample collection, viral genotyping, metadata linkage, data reporting, and public health responses. Unfortunately, current limitations in each of these steps have compromised the overall effectiveness of these strategies. Biases from convenience-based sampling methods can obfuscate the true distribution of circulating variants. The lack of standardization in genotyping strategies and bioinformatic expertise can create bottlenecks in data processing and complicate interpretation. Limitations and inconsistencies in clinical and demographic data collection and sharing can slow the compilation and limit the utility of comprehensive datasets. This likewise can complicate data reporting, restricting the availability of timely data. Finally, gaps and delays in the implementation of genomic surveillance data in the public health sphere can prevent officials from formulating effective mitigation strategies to prevent outbreaks. In this review, we outline current SARS-CoV-2 global genomic surveillance methods and assess roadblocks at each step of the pipeline to identify potential solutions. Evaluating the current obstacles that impede effective surveillance can improve both global coordination efforts and pandemic preparedness for future outbreaks.

SARS-CoV-2 Variants Identification: Overview of Molecular Existing Methods

Since the beginning of COVID-19 pandemic the Real Time sharing of genome sequences of circulating virus supported the diagnostics and surveillance of SARS-CoV-2 and its transmission dynamics. SARS-CoV-2 straightaway showed its tendency to mutate and adapt to the host, culminating in the emergence of variants; so it immediately became of crucial importance to be able to detect them quickly but also to be able to monitor in depth the changes on the whole genome to early identify the new possibly emerging variants. In this scenario, this manuscript aims to provide an overview of the existing methods for the identification of SARS-CoV-2 variants (from rapid method based on identification of one or more specific mutations to Whole Genome sequencing approach-WGS), taking into account limitations, advantages and applications of them in the field of diagnosis and surveillance of SARS-CoV-2.

Mudhune G, M. Morobe J, Said Mohammed K, O. Makori T, Ndwiga L, W. Mburu M, O. Moraa E, Musyoki J, Murunga N, N. Waliaula I, K. Mumelo A, Bejon P, Isabella Ochola-Oyier L, Githinji G, Nokes J, Agoti C. Temporal distribution and clinical characteristics of the Alpha, Delta and Omicron SARS-CoV-2 variants of concern in Laikipia, Kenya: institutional and community-based genomic surveillance

Background: Understanding the molecular epidemiology and clinical presentation of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOC) in rural-urban populations in Kenya is important for informing future public health responses and clinical care. Methods: We undertook a retrospective analysis of the clinical presentation and phylogenetic relatedness of specimens from 97 SARS-CoV-2 cases collected between 24th April and 31st December 2021 in Laikipia county, Kenya. VOC were related to observed symptoms. Phylogenetic analyses included contemporaneous sequences from across Kenya and the globe, to contextualise local transmission dynamics. Results: These sequences fell into three VOC; Alpha (n=8), Delta (n=52) and Omicron (n=37). We estimated 75 independent SARS-CoV-2 introductions into the county. The Alpha and Delta VOC were commonly detected in persons aged 31 to 45 years, 50.0% and 30.8%, respectively. The Omicron VOC was mostly detected in 16 to 30-year-olds (51.4%). Whereas relative to the other VOCs, Omicron was associated with mild upper-respiratory tract symptoms (cough, OR 3.78; 95% CI 1.1 – 16.74, p= 0.026) and sore throat, OR 22.42; 95% CI 7.11 – 81.40, p<0.001), Delta was associated with moderate to severe lower-respiratory tract symptoms (shortness of breath, OR 26.8; 95% CI 3.89 – 1158.14, p<0.001) and fever (OR 6.11; 95% CI 1.57 – 35.35, p= 0.004). Post-acute phase neurological complications were suspected in four Delta infected cases (neuralgia, neuritis, peripheral neuropathy, numbness of hand and tinnitus). Conclusion: We highlight the distinctive clinical characteristics of SARS-CoV-2 VOCs, as observed in Laikipia, Kenya, to support evidence-based clinical decisions. Multiple introductions of the VOCs were recorded despite the public health measures that were in place questioning their effectiveness during the study period.

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.

Evaluation of EPISEQ SARS-CoV-2 and a Fully Integrated Application to Identify SARS-CoV-2 Variants from Several Next-Generation Sequencing Approaches

Whole-genome sequencing has become an essential tool for real-time genomic surveillance of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) worldwide. The handling of raw next-generation sequencing (NGS) data is a major challenge for sequencing laboratories. We developed an easy-to-use web-based application (EPISEQ SARS-CoV-2) to analyse SARS-CoV-2 NGS data generated on common sequencing platforms using a variety of commercially available reagents. This application performs in one click a quality check, a reference-based genome assembly, and the analysis of the generated consensus sequence as to coverage of the reference genome, mutation screening and variant identification according to the up-to-date Nextstrain clade and Pango lineage. In this study, we validated the EPISEQ SARS-CoV-2 pipeline against a reference pipeline and compared the performance of NGS data generated by different sequencing protocols using EPISEQ SARS-CoV-2. We showed a strong agreement in SARS-CoV-2 clade and lineage identification (>99%) and in spike mutation detection (>99%) between EPISEQ SARS-CoV-2 and the reference pipeline. The comparison of several sequencing approaches using EPISEQ SARS-CoV-2 revealed 100% concordance in clade and lineage classification. It also uncovered reagent-related sequencing issues with a potential impact on SARS-CoV-2 mutation reporting. Altogether, EPISEQ SARS-CoV-2 allows an easy, rapid and reliable analysis of raw NGS data to support the sequencing efforts of laboratories with limited bioinformatics capacity and those willing to accelerate genomic surveillance of SARS-CoV-2.

Sequencing during Times of Change: Evaluating SARS-CoV-2 Clinical Samples during the Transition from the Delta to Omicron Wave

The pandemic of SARS-CoV-2 is characterized by the emergence of new variants of concern (VOCs) that supplant previous waves of infection. Here, we describe our investigation of the lineages and host-specific mutations identified in a particularly vulnerable population of predominantly older and immunosuppressed SARS-CoV-2-infected patients seen at our medical center in Chicago during the transition from the Delta to Omicron wave. We compare two primer schemes, ArticV4.1 and VarSkip2, used for short read amplicon sequencing, and describe our strategy for bioinformatics analysis that facilitates identifying lineage-associated mutations and host-specific mutations that arise during infection. This study illustrates the ongoing evolution of SARS-CoV-2 VOCs in our community and documents novel constellations of mutations that arise in individual patients. The ongoing evaluation of the evolution of SARS-CoV-2 during this pandemic is important for informing our public health strategies.

Maximum likelihood pandemic-scale phylogenetics.. [PMID: 35350209 PMCID: PMC8963701 DOI: 10.1101/2022.03.22.485312]

Genomic data plays an essential role in the study of transmissible disease, as exemplified by its current use in identifying and tracking the spread of novel SARS-CoV-2 variants. However, with the increase in size of genomic epidemiological datasets, their phylogenetic analyses become increasingly impractical due to high computational demand. In particular, while maximum likelihood methods are go-to tools for phylogenetic inference, the scale of datasets from the ongoing pandemic has made apparent the urgent need for more computationally efficient approaches. Here we propose a new likelihood-based phylogenetic framework that greatly reduces both the memory and time demand of popular maximum likelihood approaches when analysing many closely related genomes, as in the scenario of SARS-CoV-2 genome data and more generally throughout genomic epidemiology. To achieve this, we rewrite the classical Felsenstein pruning algorithm so that we can infer phylogenetic trees on at least 10 times larger datasets with higher accuracy than existing maximum likelihood methods. Our algorithms provide a powerful framework for maximum-likelihood genomic epidemiology and could facilitate similarly groundbreaking applications in Bayesian phylogenomic analyses as well.

PredictION: a predictive model to establish the performance of Oxford sequencing reads of SARS-CoV-2

The optimization of resources for research in developing countries forces us to consider strategies in the wet lab that allow the reuse of molecular biology reagents to reduce costs. In this study, we used linear regression as a method for predictive modeling of coverage depth given the number of MinION reads sequenced to define the optimum number of reads necessary to obtain >200X coverage depth with a good lineage-clade assignment of SARS-CoV-2 genomes. The research aimed to create and implement a model based on machine learning algorithms to predict different variables (e.g., coverage depth) given the number of MinION reads produced by Nanopore sequencing to maximize the yield of high-quality SARS-CoV-2 genomes, determine the best sequencing runtime, and to be able to reuse the flow cell with the remaining nanopores available for sequencing in a new run. The best accuracy was -0.98 according to the R squared performance metric of the models. A demo version is available at https://genomicdashboard.herokuapp.com/.