HIV transmission dynamics and population-wide drug resistance in rural South Africa.

Despite expanded antiretroviral therapy (ART) in South Africa, HIV-1 transmission persists. Integrase strand transfer inhibitors (INSTI) and long-acting injectables offer potential for superior viral suppression, but pre-existing drug resistance could threaten their effectiveness. In a community-based study in rural KwaZulu-Natal, prior to widespread INSTI usage, we enroled 18,025 individuals to characterise HIV-1 drug resistance and transmission networks to inform public health strategies. HIV testing and reflex viral load quantification were performed, with deep sequencing (20% variant threshold) used to detect resistance mutations. Phylogenetic and geospatial analyses characterised transmission clusters. One-third of participants were HIV-positive, with 21.7% having detectable viral loads; 62.1% of those with detectable viral loads were ART-naïve. Resistance to older reverse transcriptase (RT)-targeting drugs was found, but INSTI resistance remained low (<1%). Non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance, particularly to rilpivirine (RPV) even in ART-naïve individuals, was concerning. Twenty percent of sequenced individuals belonged to transmission clusters, with geographic analysis highlighting higher clustering in peripheral and rural areas. Our findings suggest promise for INSTI-based strategies in this setting but underscore the need for RPV resistance screening before implementing long-acting cabotegravir (CAB) + RPV. The significant clustering emphasises the importance of geographically targeted interventions to effectively curb HIV-1 transmission.

HIV transmission dynamics and population-wide drug resistance in rural South Africa.

Despite the scale-up of antiretroviral therapy (ART) in South Africa, HIV-1 incidence remains high. The anticipated use of potent integrase strand transfer inhibitors and long-acting injectables aims to enhance viral suppression at the population level and diminish transmission. Nevertheless, pre-existing drug resistance could impede the efficacy of long-acting injectable ART combinations, such as rilpivirine (an NNRTI) and cabotegravir (an INSTI). Consequently, a thorough understanding of transmission networks and geospatial distributions is vital for tailored interventions, including pre-exposure prophylaxis with long-acting injectables. However, empirical data on background resistance and transmission networks remain limited. In a community-based study in rural KwaZulu-Natal (2018-2019), prior to the widespread use of integrase inhibitor-based first-line ART, we performed HIV testing with reflex HIV-1 RNA viral load quantification on 18,025 participants. From this cohort, 6,096 (33.9%) tested positive for HIV via ELISA, with 1,323 (21.7%) exhibiting detectable viral loads (> 40 copies/mL). Of those with detectable viral loads, 62.1% were ART-naïve, and the majority of the treated were on an efavirenz + cytosine analogue + tenofovir regimen. Deep sequencing analysis, with a variant abundance threshold of 20%, revealed NRTI resistance mutations such as M184V in 2% of ART-naïve and 32% of treated individuals. Tenofovir resistance mutations K65R and K70E were found in 12% and 5% of ART-experienced individuals, respectively, and in less than 1% of ART-naïve individuals. Integrase inhibitor resistance mutations were notably infrequent (< 1%). Prevalence of pre-treatment drug resistance to NNRTIs was 10%, predominantly consisting of the K103N mutation. Among those with viraemic ART, NNRTI resistance was 50%, with rilpivirine-associated mutations observed in 9% of treated and 6% of untreated individuals. Cluster analysis revealed that 20% (205/1,050) of those sequenced were part of a cluster. We identified 171 groups with at least two linked participants; three quarters of clusters had only two individuals, and a quarter had 3-6 individuals. Integrating phylogenetic with geospatial analyses, we revealed a complex transmission network with significant clustering in specific regions, notably peripheral and rural areas. These findings derived from population scale genomic analyses are encouraging in terms of the limited resistance to DTG, but indicate that transitioning to long-acting cabotegravir + rilpivirine for transmission reduction should be accompanied by prior screening for rilpivirine resistance. Whole HIV-1 genome sequencing allowed identification of significant proportions of clusters with multiple individuals, and geospatial analyses suggesting decentralised networks can inform targeting public health interventions to effectively curb HIV-1 transmission.

Clinical severity of Omicron subvariants BA.1, BA.2, and BA.5 in a population-based cohort study in British Columbia, Canada.

The SARS-CoV-2 variant Omicron emerged in late 2021. In British Columbia (BC), Canada, and globally, three genetically distinct subvariants of Omicron, BA.1, BA.2, and BA.5, emerged and became dominant successively within an 8-month period. SARS-CoV-2 subvariants continue to circulate in the population, acquiring new mutations that have the potential to alter infectivity, immunity, and disease severity. Here, we report a propensity-matched severity analysis from residents of BC over the course of the Omicron wave, including 39,237 individuals infected with BA.1, BA.2, or BA.5 based on paired high-quality sequence data and linked to comprehensive clinical outcomes data between December 23, 2021 and August 31, 2022. Relative to BA.1, BA.2 cases were associated with a 15% and 28% lower risk of hospitalization and intensive care unit (ICU) admission (aHRhospital = 1.17; 95% confidence interval [CI] = 1.096-1.252; aHRICU = 1.368; 95% CI = 1.152-1.624), whereas BA.5 infections were associated with an 18% higher risk of hospitalization (aHRhospital = 1.18; 95% CI = 1.133-1.224) after accounting for age, sex, comorbidities, vaccination status, geography, and social determinants of health. Phylogenetic analysis revealed no specific subclades associated with more severe clinical outcomes for any Omicron subvariant. In summary, BA.1, BA.2, and BA.5 subvariants were associated with differences in clinical severity, emphasizing how variant-specific monitoring programs remain critical components of patient and population-level public health responses as the pandemic continues.

Cov2clusters: genomic clustering of SARS-CoV-2 sequences.

BACKGROUND: The COVID-19 pandemic remains a global public health concern. Advances in sequencing technologies has allowed for high numbers of SARS-CoV-2 whole genome sequence (WGS) data and rapid sharing of sequences through global repositories to enable almost real-time genomic analysis of the pathogen. WGS data has been used previously to group genetically similar viral pathogens to reveal evidence of transmission, including methods that identify distinct clusters on a phylogenetic tree. Identifying clusters of linked cases can aid in the regional surveillance and management of the disease. In this study, we present a novel method for producing stable genomic clusters of SARS-CoV-2 cases, cov2clusters, and compare the accuracy and stability of our approach to previous methods used for phylogenetic clustering using real-world SARS-CoV-2 sequence data obtained from British Columbia, Canada. RESULTS: We found that cov2clusters produced more stable clusters than previously used phylogenetic clustering methods when adding sequence data through time, mimicking an increase in sequence data through the pandemic. Our method also showed high accuracy when predicting epidemiologically informed clusters from sequence data. CONCLUSIONS: Our new approach allows for the identification of stable clusters of SARS-CoV-2 from WGS data. Producing high-resolution SARS-CoV-2 clusters from sequence data alone can a challenge and, where possible, both genomic and epidemiological data should be used in combination.

Transmission of B.1.617.2 Delta variant between vaccinated healthcare workers.

Breakthrough infections with SARS-CoV-2 Delta variant have been reported in doubly-vaccinated recipients and as re-infections. Studies of viral spread within hospital settings have highlighted the potential for transmission between doubly-vaccinated patients and health care workers and have highlighted the benefits of high-grade respiratory protection for health care workers. However the extent to which vaccination is preventative of viral spread in health care settings is less well studied. Here, we analysed data from 118 vaccinated health care workers (HCW) across two hospitals in India, constructing two probable transmission networks involving six HCWs in Hospital A and eight HCWs in Hospital B from epidemiological and virus genome sequence data, using a suite of computational approaches. A maximum likelihood reconstruction of transmission involving known cases of infection suggests a high probability that doubly vaccinated HCWs transmitted SARS-CoV-2 between each other and highlights potential cases of virus transmission between individuals who had received two doses of vaccine. Our findings show firstly that vaccination may reduce rates of transmission, supporting the need for ongoing infection control measures even in highly vaccinated populations, and secondly we have described a novel approach to identifying transmissions that is scalable and rapid, without the need for an infection control infrastructure.

Epidemiological analysis of the emergence and disappearance of the SARS-CoV-2 Kappa variant within a region of British Columbia, Canada.

Background: The Kappa variant is designated as a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant of interest (VOI). We identified 195 Kappa variant cases in a region of British Columbia, Canada-the largest published cluster in North America. Objectives: To describe the epidemiology of the Kappa variant in relation to other circulating SARS-CoV-2 variants of concern (VOC) in the region to determine if the epidemiology of the Kappa variant supports a VOI or VOC status. Methods: Clinical specimens testing positive for SARS-CoV-2 collected between March 10 and May 2, 2021, were screened for the detection of known circulating VOCs; approximately 50% of specimens were subsequently selected for whole genome sequencing (WGS). Epidemiological analysis was performed comparing the characteristics of Kappa cases to the main circulating variants in the region (Alpha and Gamma) and to non-VOC/VOI cases. Results: A total of 2,079 coronavirus disease 2019 (COVID-19) cases were reported in the region during the study period, of which 54% were selected for WGS. The 1,131 sequenced cases were categorized into Kappa, Alpha, Gamma and non-VOC/VOI. While Alpha and Gamma cases were found to have a significantly higher attack rate among household contacts compared to non-VOI/VOC cases, Kappa was not. Conclusion: Epidemiological analysis supports the designation of Kappa as a VOI and not a VOC. The Alpha and Gamma variants were found to be more transmissible, explaining their subsequent dominance in the region and the rapid disappearance of the Kappa variant. Variant surveillance strategies should focus on both detection of established VOCs and detection of potential new VOCs.

Pre-treatment integrase inhibitor resistance is uncommon in antiretroviral therapy-naive individuals with HIV-1 subtype A1 and D infections in Uganda.

OBJECTIVE: Dolutegravir (DTG) is now a preferred component of first-line antiretroviral therapy (ART). However, prevalence data on natural resistance to integrase inhibitors [integrase strand transfer inhibitors (INSTIs)] in circulating non-subtype B HIV-1 in sub-Saharan Africa is scarce. Our objective is to report prevalence of pre-treatment integrase polymorphisms associated with resistance to INSTIs in an ART-naive cohort with diverse HIV-1 subtypes. DESIGN: We retrospectively examined HIV-1 integrase sequences from Uganda. METHODS: Plasma samples were derived from the Uganda AIDS Rural Treatment Outcomes (UARTO) cohort, reflecting enrollment from 2002 to 2010, prior to initiation of ART. HIV-1 integrase was amplified using nested-PCR and Sanger-sequenced (HXB2 4230-5093). Stanford HIVdb v8.8 was used to infer clinically significant INSTI-associated mutations. Human leukocyte antigen (HLA) typing was performed for all study participants. RESULTS: Plasma samples from 511 ART-naive individuals (subtype: 48% A1, 39% D) yielded HIV-1 integrase genotyping results. Six out of 511 participants (1.2%) had any major INSTI-associated mutations. Of these, two had E138T (subtype A1), three had E138E/K (subtype D), and one had T66T/I (subtype D). No participants had mutations traditionally associated with high levels of INSTI resistance. HLA genotypes A∗02:01/05/14, B∗44:15, and C∗04:07 predicted the presence of L74I, a mutation recently observed in association with long-acting INSTI cabotegravir virologic failure. CONCLUSION: We detected no HIV-1 polymorphisms associated with high levels of DTG resistance in Uganda in the pre-DTG era. Our results support widespread implementation of DTG but careful monitoring of patients on INSTI with virologic failure is warranted to determine if unique mutations predict failure for non-B subtypes of HIV-1.

Improvements to the ARTIC multiplex PCR method for SARS-CoV-2 genome sequencing using nanopore.

Genome sequencing has been widely deployed to study the evolution of SARS-CoV-2 with more than 90,000 genome sequences uploaded to the GISAID database. We published a method for SARS-CoV-2 genome sequencing (https://www.protocols.io/view/ncov-2019-sequencing-protocol-bbmuik6w) online on January 22, 2020. This approach has rapidly become the most popular method for sequencing SARS-CoV-2 due to its simplicity and cost-effectiveness. Here we present improvements to the original protocol: i) an updated primer scheme with 22 additional primers to improve genome coverage, ii) a streamlined library preparation workflow which improves demultiplexing rate for up to 96 samples and reduces hands-on time by several hours and iii) cost savings which bring the reagent cost down to £10 per sample making it practical for individual labs to sequence thousands of SARS-CoV-2 genomes to support national and international genomic epidemiology efforts.

Phylogenetic surveillance of travel-related Zika virus infections through whole-genome sequencing methods.

In 2018, the World Health Organization identified the Zika virus (ZIKV) as a pathogen that should be prioritized for public health research due to its epidemic potential. In this study, whole-genome sequencing (WGS) of travel-acquired ZIKV infections was used to examine the limitations of phylogenetic analysis. WGS and phylogenetic analysis were performed to investigate geographic clustering of samples from five Canadians with travel-acquired ZIKV infections and to assess the limitations of phylogenetic analysis of ZIKV sequences using a phylogenetic cluster approach. Genomic variability of ZIKV samples was assessed and for context, compared with hepatitis C virus (HCV) samples. Phylogenetic analysis confirmed the suspected region of ZIKV infection for one of five samples and one sample failed to cluster with sequences from its suspected country of infection. Travel-acquired ZIKV samples depicted low genomic variability relative to HCV samples. A floating patristic distance threshold classified all pre-2000 ZIKV sequences into separate clusters, while only Cambodian, Peruvian, Malaysian, and South Korean sequences were similarly classifiable. While phylogenetic analysis of ZIKV data can identify the broad geographical region of ZIKV infection, ZIKV's low genomic variability is likely to limit precise interpretations of phylogenetic analysis of the origins of travel-related cases.

Prevalence of Human Immunodeficiency Virus-1 Integrase Strand Transfer Inhibitor Resistance in British Columbia, Canada Between 2009 and 2016: A Longitudinal Analysis.

BACKGROUND: Integrase strand transfer inhibitors (INSTIs) are highly efficacious and well tolerated antiretrovirals with fewer adverse side-effects relative to other classes of antiretrovirals. The use of INSTIs raltegravir, elvitegravir, and dolutegravir has increased dramatically over recent years. However, there is limited information about the evolution and prevalence of INSTI resistance mutations in clinical human immunodeficiency virus populations. METHODS: Human immunodeficiency virus-1-positive individuals ≥19 years were included if they received ≥1 dispensed prescription of antiretroviral therapy (ART) in British Columbia between 2009 and 2016 (N = 9358). Physician-ordered drug resistance tests were analyzed and protease inhibitor (PI), reverse-transcriptase inhibitor (RT), and INSTI resistance were defined as having ≥1 sample with a combined, cumulative score ≥30 by Stanford HIV Drug Resistance Algorithm version 7.0.1. RESULTS: Although most ART-treated individuals were tested for PI and RT resistance, INSTI resistance testing lagged behind the uptake of INSTIs among INSTI-treated individuals (11% in 2009; 34% in 2016). The prevalence of INSTI resistance was relatively low, but it increased from 1 to 7 per 1000 ART-treated individuals between 2009 and 2016 (P < .0001, R2 = 0.98). Integrase strand transfer inhibitor resistance mutations increased at integrase codons 66, 97, 140, 148, 155, and 263. CONCLUSIONS: The prevalence of INSTI resistance remains low compared with PI and RT resistance in ART-treated populations but is expanding with increased INSTI use.

Increasing Prevalence of HIV Pretreatment Drug Resistance in Women But Not Men in Rural Uganda During 2005-2013.

The prevalence of HIV pretreatment drug resistance (PDR) is increasing in sub-Saharan Africa. We sought to describe correlates of PDR and evaluate effects of PDR on clinical outcomes in rural Uganda. We analyzed data from the Uganda AIDS Rural Treatment Outcomes study, a cohort of antiretroviral therapy (ART)-naive adults with HIV (2005-2015). We performed resistance testing on pre-ART specimens. We defined PDR as any World Health Organization (WHO) 2009 surveillance drug resistance mutation and classified PDR level using the Stanford algorithm. We fit unadjusted and sex-stratified log binomial regression and Cox proportional hazard models to identify correlates of PDR and the impact of PDR on viral suppression, loss to follow-up (LTFU), and death. We analyzed data from 738 participants (median age 33 years, 69% female). Overall, prevalence of PDR was 3.5% (n = 26), owing mostly to resistance to non-nucleoside reverse transcriptase inhibitors. PDR increased over time in women (1.8% in those enrolling in clinic in 2001-2006, vs. 7.0% in 2007-2013; p = 0.006), but not in men (1.15% vs. 0.72%, p = 0.737). Lower pre-ART log10 HIV RNA was also associated with higher prevalence of PDR. We identified longer time to viral suppression among those with PDR compared with without PDR (0.5 and 0.3 years, respectively, p = 0.023), but there was no significant relationship with mortality or LTFU (p = 0.139). We observed increasing rates of PDR in women in southwestern Uganda. Implications of this trend, particularly to prevention of mother-to-child transmission programs in the region, require attention due to delayed viral suppression among those with PDR.