An Oxford Nanopore Technology-Based Hepatitis B Virus Sequencing Protocol Suitable For Genomic Surveillance Within Clinical Diagnostic Settings

Chronic hepatitis B virus (HBV) infection remains a significant public health concern, particularly in Africa, where there is a substantial burden. HBV is an enveloped virus, with isolates being classified into ten phylogenetically distinct genotypes (A - J) determined based on full-genome sequence data or reverse hybridization-based diagnostic tests. In practice, limitations are noted in that diagnostic sequencing, generally using Sanger sequencing, tends to focus only on the S-gene, yielding little or no information on intra-patient HBV genetic diversity with very low-frequency variants and reverse hybridization detects only known genotype-specific mutations. To resolve these limitations, we developed an Oxford Nanopore Technology (ONT)-based HBV genotyping protocol suitable for clinical virology, yielding complete HBV genome sequences and extensive data on intra-patient HBV diversity. Specifically, the protocol involves tiling-based PCR amplification of HBV sequences, library preparation using the ONT Rapid Barcoding Kit, ONT GridION sequencing, genotyping using Genome Detective software, recombination analysis using jpHMM and RDP5 software, and drug resistance profiling using Geno2pheno software. We prove the utility of our protocol by efficiently generating and characterizing high-quality near full-length HBV genomes from 148 left-over diagnostic Hepatitis B patient samples obtained in the Western Cape province of South Africa, providing valuable insights into the genetic diversity and epidemiology of HBV in this region of the world.

Evaluation of miniaturized Illumina DNA preparation protocols for SARS-CoV-2 whole genome sequencing

The global pandemic caused by SARS-CoV-2 has increased the demand for scalable sequencing and diagnostic methods, especially for genomic surveillance. Although next-generation sequencing has enabled large-scale genomic surveillance, the ability to sequence SARS-CoV-2 in some settings has been limited by the cost of sequencing kits and the time-consuming preparations of sequencing libraries. We compared the sequencing outcomes, cost and turn-around times obtained using the standard Illumina DNA Prep kit protocol to three modified protocols with fewer clean-up steps and different reagent volumes (full volume, half volume, one-tenth volume). We processed a single run of 47 samples under each protocol and compared the yield and mean sequence coverage. The sequencing success rate and quality for the four different reactions were as follows: the full reaction was 98.2%, the one-tenth reaction was 98.0%, the full rapid reaction was 97.5% and the half-reaction, was 97.1%. As a result, uniformity of sequence quality indicated that libraries were not affected by the change in protocol. The cost of sequencing was reduced approximately seven-fold and the time taken to prepare the library was reduced from 6.5 hours to 3 hours. The sequencing results obtained using the miniaturised volumes showed comparability to the results obtained using full volumes as described by the manufacturer. The adaptation of the protocol represents a lower-cost, streamlined approach for SARS-CoV-2 sequencing, which can be used to produce genomic data quickly and more affordably, especially in resource-constrained settings.

Genomic epidemiology of SARS-CoV-2 during the first four waves in Mozambique

Mozambique reported the first case of coronavirus disease 2019 (COVID-19) in March 2020 and it has since spread to all provinces in the country. To investigate the introductions and spread of SARS-CoV-2 in Mozambique, 1 142 whole genome sequences sampled within Mozambique were phylogenetically analyzed against a globally representative set, reflecting the first 25 months of the epidemic. The epidemic in the country was marked by four waves of infection, the first associated with B.1 ancestral lineages, while the Beta, Delta, and Omicron Variants of Concern (VOCs) were responsible for most infections and deaths during the second, third, and fourth waves. Large-scale viral exchanges occurred during the latter three waves and were largely attributed to southern African origins. Not only did the country remain vulnerable to the introductions of new variants but these variants continued to evolve within the borders of the country. Due to the Mozambican health system already under constraint, and paucity of data in Mozambique, there is a need to continue to strengthen and support genomic surveillance in the country as VOCs and Variants of interests (VOIs) are often reported from the southern African region.

The evolving SARS-CoV-2 epidemic in Africa: Insights from rapidly expanding genomic surveillance

Investment in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequencing in Africa over the past year has led to a major increase in the number of sequences that have been generated and used to track the pandemic on the continent, a number that now exceeds 100,000 genomes. Our results show an increase in the number of African countries that are able to sequence domestically and highlight that local sequencing enables faster turnaround times and more-regular routine surveillance. Despite limitations of low testing proportions, findings from this genomic surveillance study underscore the heterogeneous nature of the pandemic and illuminate the distinct dispersal dynamics of variants of concern-particularly Alpha, Beta, Delta, and Omicron-on the continent. Sustained investment for diagnostics and genomic surveillance in Africa is needed as the virus continues to evolve while the continent faces many emerging and reemerging infectious disease threats. These investments are crucial for pandemic preparedness and response and will serve the health of the continent well into the 21st century.

Emergence of SARS-CoV-2 Omicron lineages BA.4 and BA.5 in South Africa

Three lineages (BA.1, BA.2 and BA.3) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant of concern predominantly drove South Africa's fourth Coronavirus Disease 2019 (COVID-19) wave. We have now identified two new lineages, BA.4 and BA.5, responsible for a fifth wave of infections. The spike proteins of BA.4 and BA.5 are identical, and similar to BA.2 except for the addition of 69-70 deletion (present in the Alpha variant and the BA.1 lineage), L452R (present in the Delta variant), F486V and the wild-type amino acid at Q493. The two lineages differ only outside of the spike region. The 69-70 deletion in spike allows these lineages to be identified by the proxy marker of S-gene target failure, on the background of variants not possessing this feature. BA.4 and BA.5 have rapidly replaced BA.2, reaching more than 50% of sequenced cases in South Africa by the first week of April 2022. Using a multinomial logistic regression model, we estimated growth advantages for BA.4 and BA.5 of 0.08 (95% confidence interval (CI): 0.08-0.09) and 0.10 (95% CI: 0.09-0.11) per day, respectively, over BA.2 in South Africa. The continued discovery of genetically diverse Omicron lineages points to the hypothesis that a discrete reservoir, such as human chronic infections and/or animal hosts, is potentially contributing to further evolution and dispersal of the virus.

Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa

The SARS-CoV-2 epidemic in southern Africa has been characterized by three distinct waves. The first was associated with a mix of SARS-CoV-2 lineages, while the second and third waves were driven by the Beta (B.1.351) and Delta (B.1.617.2) variants, respectively1-3. In November 2021, genomic surveillance teams in South Africa and Botswana detected a new SARS-CoV-2 variant associated with a rapid resurgence of infections in Gauteng province, South Africa. Within three days of the first genome being uploaded, it was designated a variant of concern (Omicron, B.1.1.529) by the World Health Organization and, within three weeks, had been identified in 87 countries. The Omicron variant is exceptional for carrying over 30 mutations in the spike glycoprotein, which are predicted to influence antibody neutralization and spike function4. Here we describe the genomic profile and early transmission dynamics of Omicron, highlighting the rapid spread in regions with high levels of population immunity.

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.

A year of genomic surveillance reveals how the SARS-CoV-2 pandemic unfolded in Africa

The progression of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic in Africa has so far been heterogeneous, and the full impact is not yet well understood. In this study, we describe the genomic epidemiology using a dataset of 8746 genomes from 33 African countries and two overseas territories. We show that the epidemics in most countries were initiated by importations predominantly from Europe, which diminished after the early introduction of international travel restrictions. As the pandemic progressed, ongoing transmission in many countries and increasing mobility led to the emergence and spread within the continent of many variants of concern and interest, such as B.1.351, B.1.525, A.23.1, and C.1.1. Although distorted by low sampling numbers and blind spots, the findings highlight that Africa must not be left behind in the global pandemic response, otherwise it could become a source for new variants.

Implementation of an efficient SARS-CoV-2 specimen pooling strategy for high throughput diagnostic testing

The rapid identification and isolation of infected individuals remains a key strategy for controlling the spread of SARS-CoV-2. Frequent testing of populations to detect infection early in asymptomatic or presymptomatic individuals can be a powerful tool for intercepting transmission, especially when the viral prevalence is low. However, RT-PCR testing-the gold standard of SARS-CoV-2 diagnosis-is expensive, making regular testing of every individual unfeasible. Sample pooling is one approach to lowering costs. By combining samples and testing them in groups the number of tests required is reduced, substantially lowering costs. Here we report on the implementation of pooling strategies using 3-d and 4-d hypercubes to test a professional sports team in South Africa. We have shown that infected samples can be reliably detected in groups of 27 and 81, with minimal loss of assay sensitivity for samples with individual Ct values of up to 32. We report on the automation of sample pooling, using a liquid-handling robot and an automated web interface to identify positive samples. We conclude that hypercube pooling allows for the reliable RT-PCR detection of SARS-CoV-2 infection, at significantly lower costs than lateral flow antigen (LFA) tests.

Multiplex qPCR discriminates variants of concern to enhance global surveillance of SARS-CoV-2

With the emergence of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants that may increase transmissibility and/or cause escape from immune responses, there is an urgent need for the targeted surveillance of circulating lineages. It was found that the B.1.1.7 (also 501Y.V1) variant, first detected in the United Kingdom, could be serendipitously detected by the Thermo Fisher TaqPath COVID-19 PCR assay because a key deletion in these viruses, spike Δ69-70, would cause a "spike gene target failure" (SGTF) result. However, a SGTF result is not definitive for B.1.1.7, and this assay cannot detect other variants of concern (VOC) that lack spike Δ69-70, such as B.1.351 (also 501Y.V2), detected in South Africa, and P.1 (also 501Y.V3), recently detected in Brazil. We identified a deletion in the ORF1a gene (ORF1a Δ3675-3677) in all 3 variants, which has not yet been widely detected in other SARS-CoV-2 lineages. Using ORF1a Δ3675-3677 as the primary target and spike Δ69-70 to differentiate, we designed and validated an open-source PCR assay to detect SARS-CoV-2 VOC. Our assay can be rapidly deployed in laboratories around the world to enhance surveillance for the local emergence and spread of B.1.1.7, B.1.351, and P.1.

PCR assay to enhance global surveillance for SARS-CoV-2 variants of concern

With the emergence of SARS-CoV-2 variants that may increase transmissibility and/or cause escape from immune responses 1-3 , there is an urgent need for the targeted surveillance of circulating lineages. It was found that the B.1.1.7 (also 501Y.V1) variant first detected in the UK 4,5 could be serendipitously detected by the ThermoFisher TaqPath COVID-19 PCR assay because a key deletion in these viruses, spike Δ69-70, would cause a "spike gene target failure" (SGTF) result. However, a SGTF result is not definitive for B.1.1.7, and this assay cannot detect other variants of concern that lack spike Δ69-70, such as B.1.351 (also 501Y.V2) detected in South Africa 6 and P.1 (also 501Y.V3) recently detected in Brazil 7 . We identified a deletion in the ORF1a gene (ORF1a Δ3675-3677) in all three variants, which has not yet been widely detected in other SARS-CoV-2 lineages. Using ORF1a Δ3675-3677 as the primary target and spike Δ69-70 to differentiate, we designed and validated an open source PCR assay to detect SARS-CoV-2 variants of concern 8 . Our assay can be rapidly deployed in laboratories around the world to enhance surveillance for the local emergence spread of B.1.1.7, B.1.351, and P.1.