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.

Virologic Nonsuppression and HIV Drug Resistance Among People Who Inject Drugs and Their Sexual and Injecting Partners in Kenya

We evaluated the prevalence and correlates of HIV viral nonsuppression and HIV drug resistance (HIV-DR) in a cohort of people who inject drugs living with HIV (PWID-LH) and their sexual and injecting partners living with HIV in Kenya. HIV-DR testing was performed on participants with viral nonsuppression. Of 859 PWID-LH and their partners, 623 (72.5%) were on antiretroviral therapy (ART) ≥4 months and 148/623 (23.8%) were not virally suppressed. Viral nonsuppression was more common among younger participants and those on ART for a shorter duration. Among 122/148 (82.4%) successfully sequenced samples, 55 (45.1%) had detectable major HIV-DR mutations, mainly to non-nucleoside and nucleotide reverse transcriptase inhibitors (NNRTI and NRTI). High levels of HIV-DR among those with viral nonsuppression suggests need for viral load monitoring, adherence counseling, and timely switching to alternate ART regimens in this key population.

An emerging clade of Chikungunya West African genotype discovered in real-time during 2023 outbreak in Senegal

Chikungunya (CHIKV) is a re-emerging endemic arbovirus in West Africa. Since July 2023, Senegal and Burkina Faso have been experiencing an ongoing outbreak, with over 300 confirmed cases detected so far in the regions of Kédougou and Tambacounda in Senegal, the largest recorded outbreak yet. CHIKV is typically maintained in a sylvatic cycle in Senegal but its evolution and factors contributing to re-emergence are so far unknown in West Africa, leaving a gap in understanding and responding to recurrent epidemics. We produced, in real-time, the first locally-generated and publicly available CHIKV whole genomes in West Africa, to characterize the genetic diversity of circulating strains, along with phylodynamic analysis to estimate time of emergence and population growth dynamics. A novel strain of the West African genotype, phylogenetically distinct from strains circulating in previous outbreaks, was identified. This suggests a likely new spillover from sylvatic cycles in rural Senegal and potential of seeding larger epidemics in urban settings in Senegal and elsewhere.

Insights into SARS-CoV-2 in Angola during the COVID-19 peak: Molecular epidemiology and genome surveillance

Background

In Angola, COVID-19 cases have been reported in all provinces, resulting in >105,000 cases and >1900 deaths. However, no detailed genomic surveillance into the introduction and spread of the SARS-CoV-2 virus has been conducted in Angola. We aimed to investigate the emergence and epidemic progression during the peak of the COVID-19 pandemic in Angola.

Methods

We generated 1210 whole-genome SARS-CoV-2 sequences, contributing West African data to the global context, that were phylogenetically compared against global strains. Virus movement events were inferred using ancestral state reconstruction.

Results

The epidemic in Angola was marked by four distinct waves of infection, dominated by 12 virus lineages, including VOCs, VOIs, and the VUM C.16, which was unique to South-Western Africa and circulated for an extended period within the region. Virus exchanges occurred between Angola and its neighboring countries, and strong links with Brazil and Portugal reflected the historical and cultural ties shared between these countries. The first case likely originated from southern Africa.

Conclusion

A lack of a robust genome surveillance network and strong dependence on out-of-country sequencing limit real-time data generation to achieve timely disease outbreak responses, which remains of the utmost importance to mitigate future disease outbreaks in Angola.

Correction: Subramoney et al. Molecular Epidemiology of SARS-CoV-2 during Five COVID-19 Waves and the Significance of Low-Frequency Lineages. Viruses 2023, 15, 1194

In the original publication [...].

Performance of rapid antigen tests in identifying Omicron BA.4 and BA.5 infections in South Africa

BACKGROUND

Concerns around accuracy and performance of rapid antigen tests continue to be raised with the emergence of new SARS-CoV-2 variants.

OBJECTIVE

To evaluate the performance of two widely used SARS-CoV-2 rapid antigen tests during BA.4/BA.5 SARS-CoV-2 wave in South Africa (May - June 2022).

STUDY DESIGN

A prospective field evaluation compared the SARS-CoV-2 Antigen Rapid test from Hangzhou AllTest Biotech (nasal swab) and the Standard Q COVID-19 Rapid Antigen test from SD Biosensor (nasopharyngeal swab) to the Abbott RealTime SARS-CoV-2 assay (nasopharyngeal swab) on samples collected from 540 study participants.

RESULTS

Overall 28.52% (154/540) were SARS-CoV-2 RT-PCR positive with median cycle number value of 12.30 (IQR 9.30-19.40). Out of the 99 successfully sequenced SARS-CoV-2 positive samples, 18 were classified as BA.4 and 56 were classified as BA.5. The overall sensitivities of the AllTest SARS-CoV-2 Ag test and Standard Q COVID-19 Ag test were 73.38% (95% CI 65.89-79.73) and 74.03% (95% CI 66.58-80.31) and their specificities were 97.41% (95% CI 95.30-98.59) and 99.22% (95% CI 97.74-99.74) respectively. Sensitivity was >90% when the cycle number value was <20. The sensitivity of both rapid tests was >90% in samples infected with Omicron sub-lineage BA.4 and BA.5.

CONCLUSION

Accuracy of tested rapid antigen tests that target the nucleocapsid SARS-CoV-2 protein, were not adversely affected by BA.4 and BA.5 Omicron sub-variants.

SARS-CoV-2 spike protein diversity at an intra-host level, among SARS-CoV-2 infected individuals in South Africa, 2020 to 2022

Intra-host diversity studies are used to characterise the mutational heterogeneity of SARS-CoV-2 infections in order to understand the impact of virus-host adaptations. This study investigated the frequency and diversity of the spike (S) protein mutations within SARS-CoV-2 infected South African individuals. The study included SARS-CoV-2 respiratory samples, from individuals of all ages, received at the National Health Laboratory Service at Charlotte Maxeke Johannesburg Academic hospital, Gauteng, South Africa, from June 2020 to May 2022. Single nucleotide polymorphism (SNP) assays and whole genome sequencing were performed on a random selection of SARS-CoV-2 positive samples. The allele frequency (AF) was determined using TaqMan Genotyper software for SNP PCR analysis and galaxy.eu for analysis of FASTQ reads from sequencing. The SNP assays identified 5.3% (50/948) of Delta cases with heterogeneity at delY144 (4%; 2/50), E484Q (6%; 3/50), N501Y (2%; 1/50) and P681H (88%; 44/50), however only heterogeneity for E484Q and delY144 were confirmed by sequencing. From sequencing we identified 9% (210/2381) of cases with Beta, Delta, Omicron BA.1, BA.2.15, and BA.4 lineages that had heterogeneity in the S protein. Heterogeneity was primarily identified at positions 19 (1.4%) with T19IR (AF 0.2-0.7), 371 (92.3%) with S371FP (AF 0.1-1.0), and 484 (1.9%) with E484AK (0.2-0.7), E484AQ (AF 0.4-0.5) and E484KQ (AF 0.1-0.4). Mutations at heterozygous amino acid positions 19, 371 and 484 are known antibody escape mutations, however the impact of the combination of multiple substitutions identified at the same position is unknown. Therefore, we hypothesise that intra-host SARS-CoV-2 quasispecies with heterogeneity in the S protein facilitate competitive advantage of variants that can completely/partially evade host's natural and vaccine-induced immune responses.

Molecular Epidemiology of SARS-CoV-2 during Five COVID-19 Waves and the Significance of Low-Frequency Lineages

SARS-CoV-2 lineages and variants of concern (VOC) have gained more efficient transmission and immune evasion properties with time. We describe the circulation of VOCs in South Africa and the potential role of low-frequency lineages on the emergence of future lineages. Whole genome sequencing was performed on SARS-CoV-2 samples from South Africa. Sequences were analysed with Nextstrain pangolin tools and Stanford University Coronavirus Antiviral & Resistance Database. In 2020, 24 lineages were detected, with B.1 (3%; 8/278), B.1.1 (16%; 45/278), B.1.1.348 (3%; 8/278), B.1.1.52 (5%; 13/278), C.1 (13%; 37/278) and C.2 (2%; 6/278) circulating during the first wave. Beta emerged late in 2020, dominating the second wave of infection. B.1 and B.1.1 continued to circulate at low frequencies in 2021 and B.1.1 re-emerged in 2022. Beta was outcompeted by Delta in 2021, which was thereafter outcompeted by Omicron sub-lineages during the 4th and 5th waves in 2022. Several significant mutations identified in VOCs were also detected in low-frequency lineages, including S68F (E protein); I82T (M protein); P13L, R203K and G204R/K (N protein); R126S (ORF3a); P323L (RdRp); and N501Y, E484K, D614G, H655Y and N679K (S protein). Low-frequency variants, together with VOCs circulating, may lead to convergence and the emergence of future lineages that may increase transmissibility, infectivity and escape vaccine-induced or natural host immunity.

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.

Slow progression of pediatric HIV associates with early CD8+ T cell PD-1 expression and a stem-like phenotype

HIV nonprogression despite persistent viremia is rare among adults who are naive to antiretroviral therapy (ART) but relatively common among ART-naive children. Previous studies indicate that ART-naive pediatric slow progressors (PSPs) adopt immune evasion strategies similar to those described in natural hosts of SIV. However, the mechanisms underlying this immunophenotype are not well understood. In a cohort of early-treated infants who underwent analytical treatment interruption (ATI) after 12 months of ART, expression of PD-1 on CD8+ T cells immediately before ATI was the main predictor of slow progression during ATI. PD-1+CD8+ T cell frequency was also negatively correlated with CCR5 and HLA-DR expression on CD4+ T cells and predicted stronger HIV-specific T lymphocyte responses. In the CD8+ T cell compartment of PSPs, we identified an enrichment of stem-like TCF-1+PD-1+ memory cells, whereas pediatric progressors and viremic adults had a terminally exhausted PD-1+CD39+ population. TCF-1+PD-1+ expression on CD8+ T cells was associated with higher proliferative activity and stronger Gag-specific effector functionality. These data prompted the hypothesis that the proliferative burst potential of stem-like HIV-specific cytotoxic cells could be exploited in therapeutic strategies to boost the antiviral response and facilitate remission in infants who received early ART with a preserved and nonexhausted T cell compartment.

Consequences of rpoB mutations missed by the GenoType MTBDRplus assay in a programmatic setting in South Africa

Background

Rifampicin resistance missed by commercial rapid molecular assays but detected by phenotypic assays may lead to discordant susceptibility results and affect patient management.

Objective

This study was conducted to evaluate the causes of rifampicin resistance missed by the GenoType MTBDRplus and its impact on the programmatic management of tuberculosis in KwaZulu-Natal, South Africa.

Methods

We analysed routine tuberculosis programme data from January 2014 to December 2014 on isolates showing rifampicin susceptibility on the GenoType MTBDRplus assay but resistance on the phenotypic agar proportion method. Whole-genome sequencing was performed on a subset of these isolates.

Results

Out of 505 patients with isoniazid mono-resistant tuberculosis on the MTBDRplus, 145 (28.7%) isolates showed both isoniazid and rifampicin resistance on the phenotypic assay. The mean time from MTBDRplus results to initiation of drug-resistant tuberculosis therapy was 93.7 days. 65.7% of the patients had received previous tuberculosis treatment. The most common mutations detected in the 36 sequenced isolates were I491F (16; 44.4%) and L452P (12; 33.3%). Among the 36 isolates, resistance to other anti-tuberculosis drugs was 69.4% for pyrazinamide, 83.3% for ethambutol, 69.4% for streptomycin, and 50% for ethionamide.

Conclusion

Missed rifampicin resistance was mostly due to the I491F mutation located outside the MTBDRplus detection area and the L452P mutation, which was not included in the initial version 2 of the MTBDRplus. This led to substantial delays in the initiation of appropriate therapy. The previous tuberculosis treatment history and the high level of resistance to other anti-tuberculosis drugs suggest an accumulation of resistance.

HIV specific CD8+ TRM-like cells in tonsils express exhaustive signatures in the absence of natural HIV control

Lymphoid tissues are an important HIV reservoir site that persists in the face of antiretroviral therapy and natural immunity. Targeting these reservoirs by harnessing the antiviral activity of local tissue-resident memory (TRM) CD8+ T-cells is of great interest, but limited data exist on TRM-like cells within lymph nodes of people living with HIV (PLWH). Here, we studied tonsil CD8+ T-cells obtained from PLWH and uninfected controls from South Africa. We show that these cells are preferentially located outside the germinal centers (GCs), the main reservoir site for HIV, and display a low cytolytic and a transcriptionally TRM-like profile distinct from blood CD8+ T-cells. In PLWH, CD8+ TRM-like cells are expanded and adopt a more cytolytic, activated, and exhausted phenotype not reversed by antiretroviral therapy (ART). This phenotype was enhanced in HIV-specific CD8+ T-cells from tonsils compared to matched blood suggesting a higher antigen burden in tonsils. Single-cell transcriptional and clonotype resolution showed that these HIV-specific CD8+ T-cells in the tonsils express heterogeneous signatures of T-cell activation, clonal expansion, and exhaustion ex-vivo. Interestingly, this signature was absent in a natural HIV controller, who expressed lower PD-1 and CXCR5 levels and reduced transcriptional evidence of T-cell activation, exhaustion, and cytolytic activity. These data provide important insights into lymphoid tissue-derived HIV-specific CD8+ TRM-like phenotypes in settings of HIV remission and highlight their potential for immunotherapy and targeting of the HIV reservoirs.

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.

Antiretroviral therapy adherence patterns, virological suppression, and emergence of drug resistance: A nested case–control study from Uganda and South Africa

Background

Relationships between distinct antiretroviral therapy (ART) adherence patterns and risk of drug resistance are not well understood.

Methods

We conducted a nested case–control analysis within a longitudinal cohort study of individuals initiating efavirenz-based ART. Primary outcomes of interest, measured at 6 and 12 months after treatment initiation, were: 1) virologic suppression, 2) virologic failure with resistance, and 3) virologic failure without resistance. Our primary exposure of interest was ART adherence, measured over the 6 months before each visit with electronic pill monitors, and categorized in three ways: 1) 6 months average adherence; 2) running adherence, defined as the proportion of days with average adherence over 9 days of less than or equal to 10%, 20%, and 30%; and 3) number of 3-, 7-, and 28-day treatment gaps in the prior 6 months

Results

We analyzed data from 166 individuals (107 had virologic failure during observation and 59 had virologic suppression at 6 and 12 months). Average adherence was higher among those with virologic suppression (median 83%, IQR 58–96%) versus those with virologic failure with resistance (median 35%, IQR 20–77%, pairwise P < 0.01) and those with virologic failure without resistance (median 21%, IQR 2–54%, pairwise P < 0.01). Although treatment gaps generally predicted virologic failure ( P < 0.01), they did not differentiate failure with and without drug resistance ( P > 0.6).

Conclusions

Average adherence patterns, but not the assessed frequency of treatment gaps, differentiated failure with versus without drug resistance among individuals initiating efavirenz-based ART. Future work should explore adherence-resistance relationships for integrase inhibitor-based regimens.

Molecular Epidemiology and Trends in HIV-1 Transmitted Drug Resistance in Mozambique 1999–2018

HIV drug resistance (HIVDR) can become a public health concern, especially in low- and middle-income countries where genotypic testing for people initiating antiretroviral therapy (ART) is not available. For first-line regimens to remain effective, levels of transmitted drug resistance (TDR) need to be monitored over time. To determine the temporal trends of TDR in Mozambique, a search for studies in PubMed and sequences in GenBank was performed. Only studies covering the pol region that described HIVDR and genetic diversity from treatment naïve patients were included. A dataset from seven published studies and one novel unpublished study conducted between 1999 and 2018 were included. The Calibrated Population Resistance tool (CPR) and REGA HIV-1 Subtyping Tool version 3 for sequences pooled by sampling year were used to determine resistance mutations and subtypes, respectively. The prevalence of HIVDR amongst treatment-naïve individuals increased over time, reaching 14.4% in 2018. The increase was most prominent for non-nucleoside reverse transcriptase inhibitors (NNRTIs), reaching 12.7% in 2018. Subtype C was predominant in all regions, but a higher genetic variability (19% non-subtype C) was observed in the north region of Mozambique. These findings confirm a higher diversity of HIV in the north of the country and an increased prevalence of NNRTI resistance among treatment naïve individuals over time.

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.

Clinical evaluation of SARS-CoV-2 rapid antigen tests during the Omicron wave in South Africa

We evaluated the performance of nasal and nasopharyngeal Standard Q COVID-19 Ag tests (SD Biosensor) and the Panbio COVID-19 Ag Rapid Test Device (nasal, Abbott) against the Abbott RealTime SARS-CoV-2 assay during the Omicron (21 M, 21K and 21L) wave in South Africa. Overall, all evaluated tests performed well with high level of sensitivity (ranging from 77.78-81.42%) and excellent specificity values (>99%). The sensitivity of rapid antigen tests increased above 90% in samples with Ct <20, and all three tests performed best within first week of symptom onset. The Panbio Ag test performed best in the context of Omicron 21L infections, that were also characterised by significantly lower Ct values compared to infections with 21K and 21 M.

Low-frequency HIV-1 drug resistance mutations in antiretroviral naïve individuals in Botswana

BACKGROUND

Individuals living with human immunodeficiency virus (HIV) who experience virological failure (VF) after combination antiretroviral therapy (cART) initiation may have had low-frequency drug resistance mutations (DRMs) at cART initiation. There are no data on low-frequency DRMs among cART-naïve HIV-positive individuals in Botswana.

METHODS

We evaluated the prevalence of low-frequency DRMs among cART-naïve individuals previously sequenced using Sanger sequencing. The generated pol amplicons were sequenced by next-generation sequencing.

RESULTS

We observed low-frequency DRMs (detected at <20% in 33/103 (32%) of the successfully sequenced individuals, of whom four also had mutations detected at >20%. K65R was the most common low-frequency DRM detected in 8 individuals. Eighty-two of the 103 individuals had follow-up viral load data while on cART. Twenty-seven of the 82 individuals harbored low-frequency DRMs. Only 12 of 82 individuals experienced VF. The following low-frequency DRMs were observed in four individuals experiencing VF: K65R, K103N, V108I, and Y188C. No statistically significant difference was observed in the prevalence of low-frequency DRMs between individuals experiencing VF (4/12) and those not experiencing VF (23/70) (P = .97). However, individuals with non-nucleoside reverse transcriptase inhibitors-associated low-frequency DRMs were 2.68 times more likely to experience VF (odds ratio, 2.68; 95% confidential interval, 0.4-13.9) compared with those without (P = .22).

CONCLUSION

Next-generation sequencing was able to detect low-frequency DRMs in this cohort in Botswana, but these DRMs did not contribute significantly to VF.

Omicron infection enhances Delta antibody immunity in vaccinated persons

The extent to which Omicron infection1-9, with or without previous vaccination, elicits protection against the previously dominant Delta (B.1.617.2) variant is unclear. Here we measured the neutralization capacity against variants of severe acute respiratory syndrome coronavirus 2 in 39 individuals in South Africa infected with the Omicron sublineage BA.1 starting at a median of 6 (interquartile range 3-9) days post symptom onset and continuing until last follow-up sample available, a median of 23 (interquartile range 19-27) days post symptoms to allow BA.1-elicited neutralizing immunity time to develop. Fifteen participants were vaccinated with Pfizer's BNT162b2 or Johnson & Johnson's Ad26.CoV2.S and had BA.1 breakthrough infections, and 24 were unvaccinated. BA.1 neutralization increased from a geometric mean 50% focus reduction neutralization test titre of 42 at enrolment to 575 at the last follow-up time point (13.6-fold) in vaccinated participants and from 46 to 272 (6.0-fold) in unvaccinated participants. Delta virus neutralization also increased, from 192 to 1,091 (5.7-fold) in vaccinated participants and from 28 to 91 (3.0-fold) in unvaccinated participants. At the last time point, unvaccinated individuals infected with BA.1 had low absolute levels of neutralization for the non-BA.1 viruses and 2.2-fold lower BA.1 neutralization, 12.0-fold lower Delta neutralization, 9.6-fold lower Beta variant neutralization, 17.9-fold lower ancestral virus neutralization and 4.8-fold lower Omicron sublineage BA.2 neutralization relative to vaccinated individuals infected with BA.1. These results indicate that hybrid immunity formed by vaccination and Omicron BA.1 infection should be protective against Delta and other variants. By contrast, infection with Omicron BA.1 alone offers limited cross-protection despite moderate enhancement.

Rapid Replacement of SARS-CoV-2 Variants by Delta and Subsequent Arrival of Omicron, Uganda, 2021

Genomic surveillance in Uganda showed rapid replacement of severe acute respiratory syndrome coronavirus 2 over time by variants, dominated by Delta. However, detection of the more transmissible Omicron variant among travelers and increasing community transmission highlight the need for near-real-time genomic surveillance and adherence to infection control measures to prevent future pandemic waves.

Comparison of SARS-CoV-2 sequencing using the ONT GridION and the Illumina MiSeq

Background

Over 4 million SARS-CoV-2 genomes have been sequenced globally in the past 2 years. This has been crucial in elucidating transmission chains within communities, the development of new diagnostic methods, vaccines, and antivirals. Although several sequencing technologies have been employed, Illumina and Oxford Nanopore remain the two most commonly used platforms. The sequence quality between these two platforms warrants a comparison of the genomes produced by the two technologies. Here, we compared the SARS-CoV-2 consensus genomes obtained from the Oxford Nanopore Technology GridION and the Illumina MiSeq for 28 sequencing runs.

Results

Our results show that the MiSeq had a significantly higher number of consensus genomes classified by Nextclade as good and mediocre compared to the GridION. The MiSeq also had a significantly higher genome coverage and mutation counts than the GridION.

Conclusion

Due to the low genome coverage, high number of indels, and sensitivity to SARS-CoV-2 viral load noted with the GridION when compared to MiSeq, we can conclude that the MiSeq is more favourable for SARS-CoV-2 genomic surveillance, as successful genomic surveillance is dependent on high quality, near-whole consensus genomes.

Identification of SARS-CoV-2 Omicron variant using spike gene target failure and genotyping assays, Gauteng, South Africa, 2021

The circulation of Omicron BA.1 led to the rapid increase in severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) cases in South Africa in November 2021, which warranted the use of more rapid detection methods. We, therefore, assessed the ability to detect Omicron BA.1 using genotyping assays to identify specific mutations in SARS‐CoV‐2 positive samples, Gauteng province, South Africa. The TaqPath™ COVID‐19 real‐time polymerase chain reaction assay was performed on all samples selected to identify spike gene target failure (SGTF). SARS‐CoV‐2 genotyping assays were used for the detection of del69/70 and K417N mutation. Whole‐genome sequencing was performed on a subset of genotyped samples to confirm these findings. Of the positive samples received, 11.0% (175/1589) were randomly selected to assess if SGTF and genotyping assays, that detect del69/70 and K417N mutations, could identify Omicron BA.1. We identified SGTF in 98.9% (173/175) of samples, of which 88.0% (154/175) had both the del69/70 and K417N mutation. The genotyped samples (45.7%; 80/175) that were sequenced confirmed Omicron BA.1 (97.5%; 78/80). Our data show that genotyping for the detection of the del69/70 and K417N coupled with SGTF is efficient to exclude Alpha and Beta variants and rapidly detect Omicron BA.1. However, we still require assays for the detection of unique mutations that will allow for the differentiation between other Omicron sublineages. Therefore, the use of genotyping assays to detect new dominant or emerging lineages of SARS‐CoV‐2 will be beneficial in limited‐resource settings.

Selection Analysis Identifies Clusters of Unusual Mutational Changes in Omicron Lineage BA.1 That Likely Impact Spike Function

Among the 30 nonsynonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (1) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (2) interactions of Spike with ACE2 receptors, and (3) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any virus within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and, in combination with other mutations, adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron overall previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.

Persistent Hepatitis B Viraemia with Polymerase Mutations among HIV/HBV Co-Infected Patients on HBV-Active ART in KwaZulu-Natal, South Africa

To understand the problem of persistent Hepatitis B virus (HBV) viraemia in HIV/HBV co-infected patients on HBV-active antiretroviral therapy (ART), we assessed the rate of HBV virological response in patients on HBV-active ART in KwaZulu-Natal, South Africa, and analysed factors associated with persistent HBV viraemia. One hundred and fifty eligible participants with a chronic HBV diagnosis, with or without HIV coinfection, were enrolled and followed up after 6 months. The HBV pol gene was sequenced by next-generation sequencing and mutations were determined using the Stanford HBVseq database. Logistic regression analysis was used to assess factors associated with HBV viraemia at 6-month follow-up. The mean duration of HBV-active ART was 24 months. Thirty-seven of one hundred and six (35%) participants receiving HBV-active ART for longer than 6 months had virological failure. Advanced immunosuppression with CD4+ cell counts <200 cells/μL was independently associated with persistent HBV viraemia, aOR 5.276 (95% CI 1.575−17.670) p = 0.007. A high proportion of patients on HBV-active ART are unsuppressed, which will ultimately have an impact on global elimination goals. Better monitoring should be implemented, especially in HIV-coinfected patients with low CD4+ cell counts and followed by early HBV drug-resistance testing.

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.

T cell responses to SARS-CoV-2 spike cross-recognize Omicron

The SARS-CoV-2 Omicron variant (B.1.1.529) has multiple spike protein mutations1,2 that contribute to viral escape from antibody neutralization3-6 and reduce vaccine protection from infection7,8. The extent to which other components of the adaptive response such as T cells may still target Omicron and contribute to protection from severe outcomes is unknown. Here we assessed the ability of T cells to react to Omicron spike protein in participants who were vaccinated with Ad26.CoV2.S or BNT162b2, or unvaccinated convalescent COVID-19 patients (n = 70). Between 70% and 80% of the CD4+ and CD8+ T cell response to spike was maintained across study groups. Moreover, the magnitude of Omicron cross-reactive T cells was similar for Beta (B.1.351) and Delta (B.1.617.2) variants, despite Omicron harbouring considerably more mutations. In patients who were hospitalized with Omicron infections (n = 19), there were comparable T cell responses to ancestral spike, nucleocapsid and membrane proteins to those in patients hospitalized in previous waves dominated by the ancestral, Beta or Delta variants (n = 49). Thus, despite extensive mutations and reduced susceptibility to neutralizing antibodies of Omicron, the majority of T cell responses induced by vaccination or infection cross-recognize the variant. It remains to be determined whether well-preserved T cell immunity to Omicron contributes to protection from severe COVID-19 and is linked to early clinical observations from South Africa and elsewhere9-12.

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.

SARS-CoV-2 prolonged infection during advanced HIV disease evolves extensive immune escape

Characterizing SARS-CoV-2 evolution in specific geographies may help predict properties of the variants that come from these regions. We mapped neutralization of a SARS-CoV-2 strain that evolved over 6 months from ancestral virus in a person with advanced HIV disease in South Africa; this person was infected prior to emergence of the Beta and Delta variants. We longitudinally tracked the evolved virus and tested it against self-plasma and convalescent plasma from ancestral, Beta, and Delta infections. Early virus was similar to ancestral, but it evolved a multitude of mutations found in Omicron and other variants. It showed substantial but incomplete Pfizer BNT162b2 escape, weak neutralization by self-plasma, and despite pre-dating Delta, it also showed extensive escape of Delta infection-elicited neutralization. This example is consistent with the notion that SARS-CoV-2 evolving in individual immune-compromised hosts, including those with advanced HIV disease, may gain immune escape of vaccines and enhanced escape of Delta immunity, and this has implications for vaccine breakthrough and reinfections.

Omicron extensively but incompletely escapes Pfizer BNT162b2 neutralization

The emergence of the SARS-CoV-2 variant of concern Omicron (Pango lineage B.1.1.529), first identified in Botswana and South Africa, may compromise vaccine effectiveness and lead to re-infections1. Here we investigated Omicron escape from neutralization by antibodies from South African individuals vaccinated with Pfizer BNT162b2. We used blood samples taken soon after vaccination from individuals who were vaccinated and previously infected with SARS-CoV-2 or vaccinated with no evidence of previous infection. We isolated and sequence-confirmed live Omicron virus from an infected person and observed that Omicron requires the angiotensin-converting enzyme 2 (ACE2) receptor to infect cells. We compared plasma neutralization of Omicron relative to an ancestral SARS-CoV-2 strain and found that neutralization of ancestral virus was much higher in infected and vaccinated individuals compared with the vaccinated-only participants. However, both groups showed a 22-fold reduction in vaccine-elicited neutralization by the Omicron variant. Participants who were vaccinated and had previously been infected exhibited residual neutralization of Omicron similar to the level of neutralization of the ancestral virus observed in the vaccination-only group. These data support the notion that reasonable protection against Omicron may be maintained using vaccination approaches.

Omicron infection of vaccinated individuals enhances neutralizing immunity against the Delta variant

Omicron variant (B.1.1.529) infections are rapidly expanding worldwide, often in settings where the Delta variant (B.1.617.2) was dominant. We investigated whether neutralizing immunity elicited by Omicron infection would also neutralize the Delta variant and the role of prior vaccination. We enrolled 23 South African participants infected with Omicron a median of 5 days post-symptoms onset (study baseline) with a last follow-up sample taken a median of 23 days post-symptoms onset. Ten participants were breakthrough cases vaccinated with Pfizer BNT162b2 or Johnson and Johnson Ad26.CoV2.S. In vaccinated participants, neutralization of Omicron increased from a geometric mean titer (GMT) FRNT50 of 28 to 378 (13.7-fold). Unvaccinated participants had similar Omicron neutralization at baseline but increased from 26 to only 113 (4.4-fold) at follow-up. Delta virus neutralization increased from 129 to 790, (6.1-fold) in vaccinated but only 18 to 46 (2.5-fold, not statistically significant) in unvaccinated participants. Therefore, in Omicron infected vaccinated individuals, Delta neutralization was 2.1-fold higher at follow-up relative to Omicron. In a separate group previously infected with Delta, neutralization of Delta was 22.5-fold higher than Omicron. Based on relative neutralization levels, Omicron re-infection would be expected to be more likely than Delta in Delta infected individuals, and in Omicron infected individuals who are vaccinated. This may give Omicron an advantage over Delta which may lead to decreasing Delta infections in regions with high infection frequencies and high vaccine coverage.

Omicron infection of vaccinated individuals enhances neutralizing immunity against the Delta variant

Omicron variant (B.1.1.529) infections are rapidly expanding worldwide, often in settings where the Delta variant (B.1.617.2) was dominant. We investigated whether neutralizing immunity elicited by Omicron infection would also neutralize the Delta variant and the role of prior vaccination. We enrolled 23 South African participants infected with Omicron a median of 5 days post-symptoms onset (study baseline) with a last follow-up sample taken a median of 23 days post-symptoms onset. Ten participants were breakthrough cases vaccinated with Pfizer BNT162b2 or Johnson and Johnson Ad26.CoV2.S. In vaccinated participants, neutralization of Omicron increased from a geometric mean titer (GMT) FRNT50 of 28 to 378 (13.7-fold). Unvaccinated participants had similar Omicron neutralization at baseline but increased from 26 to only 113 (4.4-fold) at follow-up. Delta virus neutralization increased from 129 to 790, (6.1-fold) in vaccinated but only 18 to 46 (2.5-fold, not statistically significant) in unvaccinated participants. Therefore, in Omicron infected vaccinated individuals, Delta neutralization was 2.1-fold higher at follow-up relative to Omicron. In a separate group previously infected with Delta, neutralization of Delta was 22.5-fold higher than Omicron. Based on relative neutralization levels, Omicron re-infection would be expected to be more likely than Delta in Delta infected individuals, and in Omicron infected individuals who are vaccinated. This may give Omicron an advantage over Delta which may lead to decreasing Delta infections in regions with high infection frequencies and high vaccine coverage.

Selection analysis identifies unusual clustered mutational changes in Omicron lineage BA.1 that likely impact Spike function

Among the 30 non-synonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (i) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (ii) interactions of Spike with ACE2 receptors, and (iii) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any genomes within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron over all previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.

SARS-CoV-2 Omicron has extensive but incomplete escape of Pfizer BNT162b2 elicited neutralization and requires ACE2 for infection

The emergence of SARS-CoV-2 Omicron, first identified in Botswana and South Africa, may compromise vaccine effectiveness and the ability of antibodies triggered by previous infection to protect against re-infection (1). Here we investigated whether Omicron escapes antibody neutralization in South Africans, either previously SARS-CoV-2 infected or uninfected, who were vaccinated with Pfizer BNT162b2. We also investigated if Omicron requires the ACE2 receptor to infect cells. We isolated and sequence confirmed live Omicron virus from an infected person in South Africa and compared plasma neutralization of this virus relative to an ancestral SARS-CoV-2 strain with the D614G mutation, observing that Omicron still required ACE2 to infect. For neutralization, blood samples were taken soon after vaccination, so that vaccine elicited neutralization was close to peak. Neutralization capacity of the D614G virus was much higher in infected and vaccinated versus vaccinated only participants but both groups had 22-fold Omicron escape from vaccine elicited neutralization. Previously infected and vaccinated individuals had residual neutralization predicted to confer 73% protection from symptomatic Omicron infection, while those without previous infection were predicted to retain only about 35%. Both groups were predicted to have substantial protection from severe disease. These data support the notion that high neutralization capacity elicited by a combination of infection and vaccination, and possibly boosting, could maintain reasonable effectiveness against Omicron. A waning neutralization response is likely to decrease vaccine effectiveness below these estimates. However, since protection from severe disease requires lower neutralization levels and involves T cell immunity, such protection may be maintained.

SARS-CoV-2 evolved during advanced HIV disease immunosuppression has Beta-like escape of vaccine and Delta infection elicited immunity

Characterizing SARS-CoV-2 evolution in specific geographies may help predict the properties of variants coming from these regions. We mapped neutralization of a SARS-CoV-2 strain that evolved over 6 months from the ancestral virus in a person with advanced HIV disease. Infection was before the emergence of the Beta variant first identified in South Africa, and the Delta variant. We compared early and late evolved virus to the ancestral, Beta, Alpha, and Delta viruses and tested against convalescent plasma from ancestral, Beta, and Delta infections. Early virus was similar to ancestral, whereas late virus was similar to Beta, exhibiting vaccine escape and, despite pre-dating Delta, strong escape of Delta-elicited neutralization. This example is consistent with the notion that variants arising in immune-compromised hosts, including those with advanced HIV disease, may evolve immune escape of vaccines and enhanced escape of Delta immunity, with implications for vaccine breakthrough and reinfections.

HIGHLIGHTS

A prolonged ancestral SARS-CoV-2 infection pre-dating the emergence of Beta and Delta resulted in evolution of a Beta-like serological phenotypeSerological phenotype includes strong escape from Delta infection elicited immunity, intermediate escape from ancestral virus immunity, and weak escape from Beta immunityEvolved virus showed substantial but incomplete escape from antibodies elicited by BNT162b2 vaccination.

GRAPHICAL ABSTRACT

HIV-1 drug resistance in adults and adolescents on protease inhibitor-based antiretroviral treatment in KwaZulu-Natal Province, South Africa

BACKGROUND

In low- and middle-income countries, increasing levels of HIV drug resistance (HIVDR) on second-line protease inhibitor (PI)-based regimens are a cause for concern, given limited drug options for third-line antiretroviral therapy (ART).

OBJECTIVES

We conducted a retrospective analysis of routine HIV-1 genotyping laboratory data from KwaZulu-Natal, in South Africa, to describe the frequency and patterns of HIVDR mutations and their consequent impact on standardized third-line regimens.

METHODS

This was a cross-sectional analysis of all HIV-1 genotypic resistance tests conducted by the National Health Laboratory Service in KwaZulu-Natal, South Africa (Jan 2015 - Dec 2016), for adults and adolescents (age ≥10 years) on second-line PI-based ART with virological failure. We assigned a third-line regimen to each record, based on a national treatment algorithm and calculated the genotypic susceptibility score (GSS) for that regimen.

RESULTS

Of 348 samples analyzed, 287 (83%) had at least one drug resistance mutation (DRM) and 114 (33%) had at least one major PI DRM. Major PI resistance was associated with longer duration on second-line ART (aOR per 6-months, 1.11, 95% CI 1.04-1.19) and older age (aOR 1.03, 95% CI 1.01-1.05). Of 112 patients requiring third-line ART, 12 (11%) had a GSS of <2 for the algorithm-assigned third-line regimen.

CONCLUSIONS

One in three people failing second-line ART had significant PI DRMs. A subgroup of these individuals had extensive HIVDR, where the predicted activity of third-line ART was suboptimal, highlighting the need for continuous evaluation of outcomes on third-line regimens and close monitoring for emergent HIV-1 integrase-inhibitor resistance.

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.

Acquired HIV drug resistance and virologic monitoring in a HIV hyper-endemic setting in KwaZulu-Natal Province, South Africa

Background

Introduction of tenofovir (TDF) plus lamivudine (3TC) and dolutegravir (DTG) in first- and second-line HIV treatment regimens in South Africa warrants characterization of acquired HIV-1 drug resistance (ADR) mutations that could impact DTG-based antiretroviral therapy (ART). In this study, we sought to determine prevalence of ADR mutations and their potential impact on susceptibility to drugs used in combination with DTG among HIV-positive adults (≥ 18 years) accessing routine care at a selected ART facility in KwaZulu-Natal, South Africa.

Methods

We enrolled adult participants in a cross-sectional study between May and September 2019. Eligible participants had a most recent documented viral load (VL) ≥ 1000 copies/mL after at least 6 months on ART. We genotyped HIV-1 reverse transcriptase and protease genes by Sanger sequencing and assessed ADR. We characterized the effect of ADR mutations on the predicted susceptibility to drugs used in combination with DTG.

Results

From 143 participants enrolled, we obtained sequence data for 115 (80%), and 92.2% (95% CI 85.7–96.4) had ADR. The proportion with ADR was similar for participants on first-line ART (65/70, 92.9%, 95% CI 84.1–97.6) and those on second-line ART (40/44, 90.9%, 95% CI 78.3–97.5), and was present for the single participant on third-line ART. Approximately 89% (62/70) of those on first-line ART had dual class NRTI and NNRTI resistance and only six (13.6%) of those on second-line ART had major PI mutations. Most participants (82%) with first-line viraemia maintained susceptibility to Zidovudine (AZT), and the majority of them had lost susceptibility to TDF (71%) and 3TC (84%). Approximately two in every five TDF-treated individuals had thymidine analogue mutations (TAMs).

Conclusions

Susceptibility to AZT among most participants with first-line viraemia suggests that a new second-line regimen of AZT + 3TC + DTG could be effective. However, atypical occurrence of TAMs in TDF-treated individuals suggests a less effective AZT + 3TC + DTG regimen in a subpopulation of patients. As most patients with first-line viraemia had at least low-level resistance to TDF and 3TC, identifying viraemia before switch to TDF + 3TC + DTG is important to avoid DTG functional monotherapy. These findings highlight a need for close monitoring of outcomes on new standardized treatment regimens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12981-021-00393-5.

Temporal Changes in Vaginal Microbiota and Genital Tract Cytokines Among South African Women Treated for Bacterial Vaginosis

The standard treatment for bacterial vaginosis (BV) with oral metronidazole is often ineffective, and recurrence rates are high among African women. BV-associated anaerobes are closely associated with genital inflammation and HIV risk, which underscores the importance of understanding the interplay between vaginal microbiota and genital inflammation in response to treatment. In this cohort study, we therefore investigated the effects of metronidazole treatment on the vaginal microbiota and genital cytokines among symptomatic South African women with BV [defined as Nugent score (NS) ≥4] using 16S rRNA gene sequencing and multiplex bead arrays. Among 56 BV-positive women, we observed short-term BV clearance (NS <4) in a proportion of women six weeks after metronidazole treatment, with more than half of these experiencing recurrence by 12 weeks post-treatment. BV treatment temporarily reduced the relative abundance of BV-associated anaerobes (particularly Gardnerella vaginalis and Atopobium vaginae) and increased lactobacilli species (mainly L. iners), resulting in significantly altered mucosal immune milieu over time. In a linear mixed model, the median concentrations of pro-inflammatory cytokines and chemokines were significantly reduced in women who cleared BV compared to pre-treatment. BV persistence and recurrence were strongly associated with mucosal cytokine profiles that may increase the risk of HIV acquisition. Concentrations of these cytokines were differentially regulated by changes in the relative abundance of BVAB1 and G. vaginalis. We conclude that metronidazole for the treatment of BV induced short-term shifts in the vaginal microbiota and mucosal cytokines, while treatment failures promoted persistent elevation of pro-inflammatory cytokine concentrations in the genital tract. These data suggest the need to improve clinical management of BV to minimize BV related reproductive risk factors.

The emergence and ongoing convergent evolution of the SARS-CoV-2 N501Y lineages

The independent emergence late in 2020 of the B.1.1.7, B.1.351, and P.1 lineages of SARS-CoV-2 prompted renewed concerns about the evolutionary capacity of this virus to overcome public health interventions and rising population immunity. Here, by examining patterns of synonymous and non-synonymous mutations that have accumulated in SARS-CoV-2 genomes since the pandemic began, we find that the emergence of these three "501Y lineages" coincided with a major global shift in the selective forces acting on various SARS-CoV-2 genes. Following their emergence, the adaptive evolution of 501Y lineage viruses has involved repeated selectively favored convergent mutations at 35 genome sites, mutations we refer to as the 501Y meta-signature. The ongoing convergence of viruses in many other lineages on this meta-signature suggests that it includes multiple mutation combinations capable of promoting the persistence of diverse SARS-CoV-2 lineages in the face of mounting host immune recognition.

Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2 with grinch

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected.

Strong correlation between urine and vaginal swab samples for bacterial vaginosis

Background

Vaginal swabs have been traditionally used for the diagnosis of bacterial vaginosis (BV). Currently, there are limited studies that have investigated the use of other sample types other than vaginal swabs for the detection of BV from South African populations. This study investigated whether urine can be used for the detection of BV-associated microorganisms in South African pregnant women.

Methods

One-hundred self-collected vaginal swabs and urine samples were obtained from women presenting for antenatal care at King Edward VIII Hospital in Durban. The BD MAX™ vaginal panel assay was used for diagnosing BV and droplet digital polymerase chain reaction was used to quantify Gardnerella vaginalis, Prevotella bivia, Atopobium vaginae and Lactobacillus crispatus. The absolute counts were determined on the QX200 Droplet Reader (Bio-Rad) using the QuantaSoft Software. Data analysis was performed with statistical computing software called R, version 3.6.1.

Results

Median copy numbers obtained for G. vaginalis and P. bivia across urine and swabs in BV-positive samples were not significantly different (p = 0.134 and p = 0.652, respectively). This was confirmed by the correlation analysis that showed a good correlation between the two sample types (G. vaginalis [r = 0.63] and P. bivia [r = 0.50]). However, the data obtained for A. vaginae differed, and a weak correlation between urine and swabs was observed (r = 0.21). Bacterial vaginosis-negative samples had no significant difference in median copy numbers for L. crispatus across the urine and swabs (p = 0.062), and a good correlation between the sample types was noted (r = 0.71).

Conclusion

This study highlights the appropriateness of urine for the detection of microorganisms associated with BV.

The emergence and ongoing convergent evolution of the N501Y lineages coincides with a major global shift in the SARS-CoV-2 selective landscape

The emergence and rapid rise in prevalence of three independent SARS-CoV-2 "501Y lineages", B.1.1.7, B.1.351 and P.1, in the last three months of 2020 prompted renewed concerns about the evolutionary capacity of SARS-CoV-2 to adapt to both rising population immunity, and public health interventions such as vaccines and social distancing. Viruses giving rise to the different 501Y lineages have, presumably under intense natural selection following a shift in host environment, independently acquired multiple unique and convergent mutations. As a consequence, all have gained epidemiological and immunological properties that will likely complicate the control of COVID-19. Here, by examining patterns of mutations that arose in SARSCoV-2 genomes during the pandemic we find evidence of a major change in the selective forces acting on various SARS-CoV-2 genes and gene segments (such as S, nsp2 and nsp6), that likely coincided with the emergence of the 501Y lineages. In addition to involving continuing sequence diversification, we find evidence that a significant portion of the ongoing adaptive evolution of the 501Y lineages also involves further convergence between the lineages. Our findings highlight the importance of monitoring how members of these known 501Y lineages, and others still undiscovered, are convergently evolving similar strategies to ensure their persistence in the face of mounting infection and vaccine induced host immune recognition.

Impact of pretreatment low-abundance HIV-1 drug-resistant variants on virological failure among HIV-1/TB-co-infected individuals

Objectives

To determine the impact of pretreatment low-abundance HIV-1 drug-resistant variants (LA-DRVs) on virological failure (VF) among HIV-1/TB-co-infected individuals treated with NNRTI first-line ART.

Methods

We conducted a case–control study of 170 adults with HIV-1/TB co-infection. Cases had at least one viral load (VL) ≥1000 RNA copies/mL after ≥6 months on NNRTI-based ART, and controls had sustained VLs <1000 copies/mL. We sequenced plasma viruses by Sanger and MiSeq next-generation sequencing (NGS). We assessed drug resistance mutations (DRMs) using the Stanford drug resistance database, and analysed NGS data for DRMs at ≥20%, 10%, 5% and 2% thresholds. We assessed the effect of pretreatment drug resistance (PDR) on VF.

Results

We analysed sequences from 45 cases and 125 controls. Overall prevalence of PDR detected at a ≥20% threshold was 4.7% (8/170) and was higher in cases than in controls (8.9% versus 3.2%), P = 0.210. Participants with PDR at ≥20% had almost 4-fold higher odds of VF (adjusted OR 3.7, 95% CI 0.8–18.3) compared with those without, P = 0.104. PDR prevalence increased to 18.2% (31/170) when LA-DRVs at ≥2% were included. Participants with pretreatment LA-DRVs only had 1.6-fold higher odds of VF (adjusted OR 1.6, 95% CI 0.6–4.3) compared with those without, P = 0.398.

Conclusions

Pretreatment DRMs and LA-DRVs increased the odds of developing VF on NNRTI-based ART, although without statistical significance. NGS increased detection of DRMs but provided no additional benefit in identifying participants at risk of VF at lower thresholds. More studies assessing mutation thresholds predictive of VF are required to inform use of NGS in treatment decisions.

Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected.

Escape of SARS-CoV-2 501Y.V2 from neutralization by convalescent plasma

SARS-CoV-2 variants of concern (VOC) have arisen independently at multiple locations^1,2 and may reduce the efficacy of current vaccines that target the spike glycoprotein of SARS-CoV-2³. Here, using a live-virus neutralization assay, we compared the neutralization of a non-VOC variant with the 501Y.V2 VOC (also known as B.1.351) using plasma collected from adults who were hospitalized with COVID-19 during the two waves of infection in South Africa, the second wave of which was dominated by infections with the 501Y.V2 variant. Sequencing demonstrated that infections of plasma donors from the first wave were with viruses that did not contain the mutations associated with 501Y.V2, except for one infection that contained the E484K substitution in the receptor-binding domain. The 501Y.V2 virus variant was effectively neutralized by plasma from individuals who were infected during the second wave. The first-wave virus variant was effectively neutralized by plasma from first-wave infections. However, the 501Y.V2 variant was poorly cross-neutralized by plasma from individuals with first-wave infections; the efficacy was reduced by 15.1-fold relative to neutralization of 501Y.V2 by plasma from individuals infected in the second wave. By contrast, cross-neutralization of first-wave virus variants using plasma from individuals with second-wave infections was more effective, showing only a 2.3-fold decrease relative to neutralization of first-wave virus variants by plasma from individuals infected in the first wave. Although we tested only one plasma sample from an individual infected with a SARS-CoV-2 variant with only the E484K substitution, this plasma sample potently neutralized both variants. The observed effective neutralization of first-wave virus by plasma from individuals infected with 501Y.V2 provides preliminary evidence that vaccines based on VOC sequences could retain activity against other circulating SARS-CoV-2 lineages.

Early transmission of SARS-CoV-2 in South Africa: An epidemiological and phylogenetic report

OBJECTIVES

The Network for Genomic Surveillance in South Africa (NGS-SA) was formed to investigate the introduction and understand the early transmission dynamics of the SARS-CoV-2 epidemic in South-Africa.

DESIGN

This paper presents the first results from this group, which is a molecular epidemiological study of the first 21 SARS-CoV-2 whole genomes sampled in the first port of entry - KwaZulu-Natal (KZN) - during the first month of the epidemic. By combining this with calculations of the effective reproduction number (R), it aimed to shed light on the patterns of infections in South Africa.

RESULTS

Two of the largest provinces - Gauteng and KZN - had a slow growth rate for the number of detected cases, while the epidemic spread faster in the Western Cape and Eastern Cape. The estimates of transmission potential suggested a decrease towards R = 1 since the first cases and deaths, but a subsequent estimated R average of 1.39 between 6-18 May 2020. It was also demonstrated that early transmission in KZN was associated with multiple international introductions and dominated by lineages B1 and B. Evidence for locally acquired infections in a hospital in Durban within the first month of the epidemic was also provided.

CONCLUSION

The COVID-19 pandemic in South Africa was very heterogeneous in its spatial dimension, with many distinct introductions of SARS-CoV2 in KZN and evidence of nosocomial transmission, which inflated early mortality in KZN. The epidemic at the local level was still developing and NGS-SA aimed to clarify the dynamics in South Africa and devise the most effective measures as the outbreak evolved.

Transmission dynamics of SARS-CoV-2 within-host diversity in two major hospital outbreaks in South Africa

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes acute, highly transmissible respiratory infection in humans and a wide range of animal species. Its rapid global spread has resulted in a major public health emergency, necessitating commensurately rapid research to improve control strategies. In particular, the ability to effectively retrace transmission chains in outbreaks remains a major challenge, partly due to our limited understanding of the virus' underlying evolutionary dynamics within and between hosts. We used high-throughput sequencing whole-genome data coupled with bottleneck analysis to retrace the pathways of viral transmission in two nosocomial outbreaks that were previously characterised by epidemiological and phylogenetic methods. Additionally, we assessed the mutational landscape, selection pressures, and diversity at the within-host level for both outbreaks. Our findings show evidence of within-host selection and transmission of variants between samples. Both bottleneck and diversity analyses highlight within-host and consensus-level variants shared by putative source-recipient pairs in both outbreaks, suggesting that certain within-host variants in these outbreaks may have been transmitted upon infection rather than arising de novo independently within multiple hosts. Overall, our findings demonstrate the utility of combining within-host diversity and bottleneck estimations for elucidating transmission events in SARS-CoV-2 outbreaks, provide insight into the maintenance of viral genetic diversity, provide a list of candidate targets of positive selection for further investigation, and demonstrate that within-host variants can be transferred between patients. Together these results will help in developing strategies to understand the nature of transmission events and curtail the spread of SARS-CoV-2.