Prevalence of persistent SARS-CoV-2 in a large community surveillance study.

Persistent SARS-CoV-2 infections may act as viral reservoirs that could seed future outbreaks1-5, give rise to highly divergent lineages6-8 and contribute to cases with post-acute COVID-19 sequelae (long COVID)9,10. However, the population prevalence of persistent infections, their viral load kinetics and evolutionary dynamics over the course of infections remain largely unknown. Here, using viral sequence data collected as part of a national infection survey, we identified 381 individuals with SARS-CoV-2 RNA at high titre persisting for at least 30 days, of which 54 had viral RNA persisting at least 60 days. We refer to these as 'persistent infections' as available evidence suggests that they represent ongoing viral replication, although the persistence of non-replicating RNA cannot be ruled out in all. Individuals with persistent infection had more than 50% higher odds of self-reporting long COVID than individuals with non-persistent infection. We estimate that 0.1-0.5% of infections may become persistent with typically rebounding high viral loads and last for at least 60 days. In some individuals, we identified many viral amino acid substitutions, indicating periods of strong positive selection, whereas others had no consensus change in the sequences for prolonged periods, consistent with weak selection. Substitutions included mutations that are lineage defining for SARS-CoV-2 variants, at target sites for monoclonal antibodies and/or are commonly found in immunocompromised people11-14. This work has profound implications for understanding and characterizing SARS-CoV-2 infection, epidemiology and evolution.

The epidemiological impact of the NHS COVID-19 app.

The COVID-19 pandemic has seen the emergence of digital contact tracing to help to prevent the spread of the disease. A mobile phone app records proximity events between app users, and when a user tests positive for COVID-19, their recent contacts can be notified instantly. Theoretical evidence has supported this new public health intervention1-6, but its epidemiological impact has remained uncertain7. Here we investigate the impact of the National Health Service (NHS) COVID-19 app for England and Wales, from its launch on 24 September 2020 to the end of December 2020. It was used regularly by approximately 16.5 million users (28% of the total population), and sent approximately 1.7 million exposure notifications: 4.2 per index case consenting to contact tracing. We estimated that the fraction of individuals notified by the app who subsequently showed symptoms and tested positive (the secondary attack rate (SAR)) was 6%, similar to the SAR for manually traced close contacts. We estimated the number of cases averted by the app using two complementary approaches: modelling based on the notifications and SAR gave an estimate of 284,000 (central 95% range of sensitivity analyses 108,000-450,000), and statistical comparison of matched neighbouring local authorities gave an estimate of 594,000 (95% confidence interval 317,000-914,000). Approximately one case was averted for each case consenting to notification of their contacts. We estimated that for every percentage point increase in app uptake, the number of cases could be reduced by 0.8% (using modelling) or 2.3% (using statistical analysis). These findings support the continued development and deployment of such apps in populations that are awaiting full protection from vaccines.

Viral burden is associated with age, vaccination, and viral variant in a population-representative study of SARS-CoV-2 that accounts for time-since-infection-related sampling bias.

In this study, we evaluated the impact of viral variant, in addition to other variables, on within-host viral burden, by analysing cycle threshold (Ct) values derived from nose and throat swabs, collected as part of the UK COVID-19 Infection Survey. Because viral burden distributions determined from community survey data can be biased due to the impact of variant epidemiology on the time-since-infection of samples, we developed a method to explicitly adjust observed Ct value distributions to account for the expected bias. By analysing the adjusted Ct values using partial least squares regression, we found that among unvaccinated individuals with no known prior exposure, viral burden was 44% lower among Alpha variant infections, compared to those with the predecessor strain, B.1.177. Vaccination reduced viral burden by 67%, and among vaccinated individuals, viral burden was 286% higher among Delta variant, compared to Alpha variant, infections. In addition, viral burden increased by 17% for every 10-year age increment of the infected individual. In summary, within-host viral burden increases with age, is reduced by vaccination, and is influenced by the interplay of vaccination status and viral variant.

Correction: Viral burden is associated with age, vaccination, and viral variant in a population-representative study of SARS-CoV-2 that accounts for time-since-infection-related sampling bias.

[This corrects the article DOI: 10.1371/journal.ppat.1011461.].

Lineage replacement and evolution captured by 3 years of the United Kingdom Coronavirus (COVID-19) Infection Survey.

The Office for National Statistics Coronavirus (COVID-19) Infection Survey (ONS-CIS) is the largest surveillance study of SARS-CoV-2 positivity in the community, and collected data on the United Kingdom (UK) epidemic from April 2020 until March 2023 before being paused. Here, we report on the epidemiological and evolutionary dynamics of SARS-CoV-2 determined by analysing the sequenced samples collected by the ONS-CIS during this period. We observed a series of sweeps or partial sweeps, with each sweeping lineage having a distinct growth advantage compared to their predecessors, although this was also accompanied by a gradual fall in average viral burdens from June 2021 to March 2023. The sweeps also generated an alternating pattern in which most samples had either S-gene target failure (SGTF) or non-SGTF over time. Evolution was characterized by steadily increasing divergence and diversity within lineages, but with step increases in divergence associated with each sweeping major lineage. This led to a faster overall rate of evolution when measured at the between-lineage level compared to within lineages, and fluctuating levels of diversity. These observations highlight the value of viral sequencing integrated into community surveillance studies to monitor the viral epidemiology and evolution of SARS-CoV-2, and potentially other pathogens.

Limited emergence of resistance to integrase strand transfer inhibitors (INSTIs) in ART-experienced participants failing dolutegravir-based antiretroviral therapy: a cross-sectional analysis of a Northeast Nigerian cohort.

BACKGROUND: Due to the high prevalence of resistance to NNRTI-based ART since 2018, consolidated recommendations from the WHO have indicated dolutegravir as the preferred drug of choice for HIV treatment globally. There is a paucity of resistance outcome data from HIV-1 non-B subtypes circulating across West Africa. AIMS: We characterized the mutational profiles of persons living with HIV from a cross-sectional cohort in North-East Nigeria failing a dolutegravir-based ART regimen. METHODS: WGS of plasma samples collected from 61 HIV-1-infected participants following virological failure of dolutegravir-based ART were sequenced using the Illumina platform. Sequencing was successfully completed for samples from 55 participants. Following quality control, 33 full genomes were analysed from participants with a median age of 40 years and median time on ART of 9 years. HIV-1 subtyping was performed using SNAPPy. RESULTS: Most participants had mutational profiles reflective of exposure to previous first- and second-line ART regimens comprised NRTIs and NNRTIs. More than half of participants had one or more drug resistance-associated mutations (DRMs) affecting susceptibility to NRTIs (17/33; 52%) and NNRTIs (24/33; 73%). Almost a quarter of participants (8/33; 24.4%) had one or more DRMs affecting tenofovir susceptibility. Only one participant, infected with HIV-1 subtype G, had evidence of DRMs affecting dolutegravir susceptibility-this was characterized by the T66A, G118R, E138K and R263K mutations. CONCLUSIONS: This study found a low prevalence of resistance to dolutegravir; the data are therefore supportive of the continual rollout of dolutegravir as the primary first-line regimen for ART-naive participants and the preferred switch to second-line ART across the region. However, population-level, longer-term data collection on dolutegravir outcomes are required to further guide implementation and policy action across the region.

Two Distinct Mechanisms Leading to Loss of Virological Control in the Rare Group of Antiretroviral Therapy-Naive, Transiently Aviremic Children Living with HIV.

HIV-specific CD8+ T cells play a central role in immune control of adult HIV, but their contribution in pediatric infection is less well characterized. Previously, we identified a group of ART-naive children with persistently undetectable plasma viremia, termed "elite controllers," and a second group who achieved aviremia only transiently. To investigate the mechanisms of failure to maintain aviremia, we characterized in three transient aviremic individuals (TAs), each of whom expressed the disease-protective HLA-B*81:01, longitudinal HIV-specific T-cell activity, and viral sequences. In two TAs, a CD8+ T-cell response targeting the immunodominant epitope TPQDLNTML (Gag-TL9) was associated with viral control, followed by viral rebound and the emergence of escape variants with lower replicative capacity. Both TAs mounted variant-specific responses, but only at low functional avidity, resulting in immunological progression. In contrast, in TA-3, intermittent viremic episodes followed aviremia without virus escape or a diminished CD4+ T-cell count. High quality and magnitude of the CD8+ T-cell response were associated with aviremia. We therefore identify two distinct mechanisms of loss of viral control. In one scenario, CD8+ T-cell responses initially cornered low-replicative-capacity escape variants, but with insufficient avidity to prevent viremia and disease progression. In the other, loss of viral control was associated with neither virus escape nor progression but with a decrease in the quality of the CD8+ T-cell response, followed by recovery of viral control in association with improved antiviral response. These data suggest the potential for a consistently strong and polyfunctional antiviral response to achieve long-term viral control without escape. IMPORTANCE Very early initiation of antiretroviral therapy (ART) in pediatric HIV infection offers a unique opportunity to limit the size and diversity of the viral reservoir. However, only rarely is ART alone sufficient to achieve remission. Additional interventions that likely include contributions from host immunity are therefore required. The HIV-specific T-cell response plays a central role in immune control of adult HIV, often mediated through protective alleles such as HLA-B*57/58:01/81:01. However, due to the tolerogenic and type 2 biased immune response in early life, HLA-I-mediated immune suppression of viremia is seldom observed in children. We assessed a rare group of HLA-B*81:01-positive, ART-naive children who achieved aviremia, albeit only transiently, and investigated the role of the CD8+ T-cell response in the establishment and loss of viral control. We identified a mechanism by which the HIV-specific response can achieve viremic control without viral escape that can be explored in strategies to achieve remission.

Validation of an HIV whole genome sequencing method for HIV drug resistance testing in an Australian clinical microbiology laboratory.

Detection of HIV drug resistance (HIVDR) is vital to successful anti-retroviral therapy (ART). HIVDR testing to determine drug-resistance mutations is routinely performed in Australia to guide ART choice in newly diagnosed people living with HIV or in cases of treatment failure. In 2022, our clinical microbiology laboratory sought to validate a next-generation sequencing (NGS)-based HIVDR assay to replace the previous Sanger-sequencing (SS)-based ViroSeq. NGS solutions for HIVDR offer higher throughput, lower costs and higher sensitivity for variant detection. We sought to validate the previously described low-cost probe-based NGS method (veSEQ-HIV) for whole-genome recovery and HIVDR-testing in a diagnostic setting. veSEQ-HIV displayed 100% and 98% accuracy in major and minor mutation detection, respectively, and 100% accuracy of subtyping (provided > 1000 mapped reads were obtained). Pairwise comparison exhibited low inter-and intrarun variability across the whole-genome (Jaccard index [J] = 0.993; J = 0.972) and the Pol gene (J = 0.999; J = 0.999), respectively. veSEQ-HIV met all our pre-set criteria based on WHO recommendations and successfully replaced ViroSeq in our laboratory. Scaling-down veSEQ-HIV to a limited batch size and sequencing on Illumina iSeq. 100, allowed easy implementation of the assay into the workflow of a small sequencing laboratory with minimal staff and equipment and the ability to meet clinically relevant test turn-around times. As HIVDR-testing moves from SS- to NGS-based methods and new ART drugs come to market (particularly those with targets outside the Pol region), whole-genome recovery using veSEQ-HIV provides a robust, cost-effective and "future-proof" NGS method for HIVDR-testing.

Evaluation of Intraperitoneal [18F]-FDOPA Administration for Micro-PET Imaging in Mice and Assessment of the Effect of Subchronic Ketamine Dosing on Dopamine Synthesis Capacity.

Positron emission tomography (PET) using the radiotracer [18F]-FDOPA provides a tool for studying brain dopamine synthesis capacity in animals and humans. We have previously standardised a micro-PET methodology in mice by intravenously administering [18F]-FDOPA via jugular vein cannulation and assessment of striatal dopamine synthesis capacity, indexed as the influx rate constant K i Mod of [18F]-FDOPA, using an extended graphical Patlak analysis with the cerebellum as a reference region. This enables a direct comparison between preclinical and clinical output values. However, chronic intravenous catheters are technically difficult to maintain for longitudinal studies. Hence, in this study, intraperitoneal administration of [18F]-FDOPA was evaluated as a less-invasive alternative that facilitates longitudinal imaging. Our experiments comprised the following assessments: (i) comparison of [18F]-FDOPA uptake between intravenous and intraperitoneal radiotracer administration and optimisation of the time window used for extended Patlak analysis, (ii) comparison of Ki Mod in a within-subject design of both administration routes, (iii) test-retest evaluation of Ki Mod in a within-subject design of intraperitoneal radiotracer administration, and (iv) validation of Ki Mod estimates by comparing the two administration routes in a mouse model of hyperdopaminergia induced by subchronic ketamine. Our results demonstrate that intraperitoneal [18F]-FDOPA administration resulted in good brain uptake, with no significant effect of administration route on Ki Mod estimates (intraperitoneal: 0.024 ± 0.0047 min-1, intravenous: 0.022 ± 0.0041 min-1, p = 0.42) and similar coefficient of variation (intraperitoneal: 19.6%; intravenous: 18.4%). The technique had a moderate test-retest validity (intraclass correlation coefficient (ICC) = 0.52, N = 6) and thus supports longitudinal studies. Following subchronic ketamine administration, elevated K i Mod as compared to control condition was measured with a large effect size for both methods (intraperitoneal: Cohen's d = 1.3; intravenous: Cohen's d = 0.9), providing further evidence that ketamine has lasting effects on the dopamine system, which could contribute to its therapeutic actions and/or abuse liability.

Efficacy of ChAdOx1 nCoV-19 (AZD1222) vaccine against SARS-CoV-2 variant of concern 202012/01 (B.1.1.7): an exploratory analysis of a randomised controlled trial.

BACKGROUND: A new variant of SARS-CoV-2, B.1.1.7, emerged as the dominant cause of COVID-19 disease in the UK from November, 2020. We report a post-hoc analysis of the efficacy of the adenoviral vector vaccine, ChAdOx1 nCoV-19 (AZD1222), against this variant. METHODS: Volunteers (aged ≥18 years) who were enrolled in phase 2/3 vaccine efficacy studies in the UK, and who were randomly assigned (1:1) to receive ChAdOx1 nCoV-19 or a meningococcal conjugate control (MenACWY) vaccine, provided upper airway swabs on a weekly basis and also if they developed symptoms of COVID-19 disease (a cough, a fever of 37·8°C or higher, shortness of breath, anosmia, or ageusia). Swabs were tested by nucleic acid amplification test (NAAT) for SARS-CoV-2 and positive samples were sequenced through the COVID-19 Genomics UK consortium. Neutralising antibody responses were measured using a live-virus microneutralisation assay against the B.1.1.7 lineage and a canonical non-B.1.1.7 lineage (Victoria). The efficacy analysis included symptomatic COVID-19 in seronegative participants with a NAAT positive swab more than 14 days after a second dose of vaccine. Participants were analysed according to vaccine received. Vaccine efficacy was calculated as 1 - relative risk (ChAdOx1 nCoV-19 vs MenACWY groups) derived from a robust Poisson regression model. This study is continuing and is registered with ClinicalTrials.gov, NCT04400838, and ISRCTN, 15281137. FINDINGS: Participants in efficacy cohorts were recruited between May 31 and Nov 13, 2020, and received booster doses between Aug 3 and Dec 30, 2020. Of 8534 participants in the primary efficacy cohort, 6636 (78%) were aged 18-55 years and 5065 (59%) were female. Between Oct 1, 2020, and Jan 14, 2021, 520 participants developed SARS-CoV-2 infection. 1466 NAAT positive nose and throat swabs were collected from these participants during the trial. Of these, 401 swabs from 311 participants were successfully sequenced. Laboratory virus neutralisation activity by vaccine-induced antibodies was lower against the B.1.1.7 variant than against the Victoria lineage (geometric mean ratio 8·9, 95% CI 7·2-11·0). Clinical vaccine efficacy against symptomatic NAAT positive infection was 70·4% (95% CI 43·6-84·5) for B.1.1.7 and 81·5% (67·9-89·4) for non-B.1.1.7 lineages. INTERPRETATION: ChAdOx1 nCoV-19 showed reduced neutralisation activity against the B.1.1.7 variant compared with a non-B.1.1.7 variant in vitro, but the vaccine showed efficacy against the B.1.1.7 variant of SARS-CoV-2. FUNDING: UK Research and Innovation, National Institute for Health Research (NIHR), Coalition for Epidemic Preparedness Innovations, NIHR Oxford Biomedical Research Centre, Thames Valley and South Midlands NIHR Clinical Research Network, and AstraZeneca.

Highly Sensitive Lineage Discrimination of SARS-CoV-2 Variants through Allele-Specific Probe PCR.

Tools to detect SARS-CoV-2 variants of concern and track the ongoing evolution of the virus are necessary to support public health efforts and the design and evaluation of novel COVID-19 therapeutics and vaccines. Although next-generation sequencing (NGS) has been adopted as the gold standard method for discriminating SARS-CoV-2 lineages, alternative methods may be required when processing samples with low viral loads or low RNA quality. To this aim, an allele-specific probe PCR (ASP-PCR) targeting lineage-specific single nucleotide polymorphisms (SNPs) was developed and used to screen 1,082 samples from two clinical trials in the United Kingdom and Brazil. Probit regression models were developed to compare ASP-PCR performance against 1,771 NGS results for the same cohorts. Individual SNPs were shown to readily identify specific variants of concern. ASP-PCR was shown to discriminate SARS-CoV-2 lineages with a higher likelihood than NGS over a wide range of viral loads. The comparative advantage for ASP-PCR over NGS was most pronounced in samples with cycle threshold (CT) values between 26 and 30 and in samples that showed evidence of degradation. Results for samples screened by ASP-PCR and NGS showed 99% concordant results. ASP-PCR is well suited to augment but not replace NGS. The method can differentiate SARS-CoV-2 lineages with high accuracy and would be best deployed to screen samples with lower viral loads or that may suffer from degradation. Future work should investigate further destabilization from primer-target base mismatch through altered oligonucleotide chemistry or chemical additives.

The mutational spectrum of SARS-CoV-2 genomic and antigenomic RNA.

The raw material for viral evolution is provided by intra-host mutations occurring during replication, transcription or post-transcription. Replication and transcription of Coronaviridae proceed through the synthesis of negative-sense 'antigenomes' acting as templates for positive-sense genomic and subgenomic RNA. Hence, mutations in the genomes of SARS-CoV-2 and other coronaviruses can occur during (and after) the synthesis of either negative-sense or positive-sense RNA, with potentially distinct patterns and consequences. We explored for the first time the mutational spectrum of SARS-CoV-2 (sub)genomic and anti(sub)genomic RNA. We use a high-quality deep sequencing dataset produced using a quantitative strand-aware sequencing method, controlled for artefacts and sequencing errors, and scrutinized for accurate detection of within-host diversity. The nucleotide differences between negative- and positive-sense strand consensus vary between patients and do not show dependence on age or sex. Similarities and differences in mutational patterns between within-host minor variants on the two RNA strands suggested strand-specific mutations or editing by host deaminases and oxidative damage. We observe generally neutral and slight negative selection on the negative strand, contrasting with purifying selection in ORF1a, ORF1b and S genes of the positive strand of the genome.

Reproducing the dopamine pathophysiology of schizophrenia and approaches to ameliorate it: a translational imaging study with ketamine.

Patients with schizophrenia show increased striatal dopamine synthesis capacity in imaging studies. The mechanism underlying this is unclear but may be due to N-methyl-D-aspartate receptor (NMDAR) hypofunction and parvalbumin (PV) neuronal dysfunction leading to disinhibition of mesostriatal dopamine neurons. Here, we develop a translational mouse model of the dopamine pathophysiology seen in schizophrenia and test approaches to reverse the dopamine changes. Mice were treated with sub-chronic ketamine (30 mg/kg) or saline and then received in vivo positron emission tomography of striatal dopamine synthesis capacity, analogous to measures used in patients. Locomotor activity was measured using the open-field test. In vivo cell-type-specific chemogenetic approaches and pharmacological interventions were used to manipulate neuronal excitability. Immunohistochemistry and RNA sequencing were used to investigate molecular mechanisms. Sub-chronic ketamine increased striatal dopamine synthesis capacity (Cohen's d = 2.5) and locomotor activity. These effects were countered by inhibition of midbrain dopamine neurons, and by activation of PV interneurons in pre-limbic cortex and ventral subiculum of the hippocampus. Sub-chronic ketamine reduced PV expression in these cortical and hippocampal regions. Pharmacological intervention with SEP-363856, a novel psychotropic agent with agonism at trace amine receptor 1 (TAAR1) and 5-HT1A receptors but no appreciable action at dopamine D2 receptors, significantly reduced the ketamine-induced increase in dopamine synthesis capacity. These results show that sub-chronic ketamine treatment in mice mimics the dopaminergic alterations in patients with psychosis, that this requires activation of midbrain dopamine neurons, and can be ameliorated by activating PV interneurons and by a TAAR1/5-HT1A agonist. This identifies novel therapeutic approaches for targeting presynaptic dopamine dysfunction in patients with schizophrenia and effects of ketamine relevant to its therapeutic use for  treating major depression.

OpenABM-Covid19-An agent-based model for non-pharmaceutical interventions against COVID-19 including contact tracing.

SARS-CoV-2 has spread across the world, causing high mortality and unprecedented restrictions on social and economic activity. Policymakers are assessing how best to navigate through the ongoing epidemic, with computational models being used to predict the spread of infection and assess the impact of public health measures. Here, we present OpenABM-Covid19: an agent-based simulation of the epidemic including detailed age-stratification and realistic social networks. By default the model is parameterised to UK demographics and calibrated to the UK epidemic, however, it can easily be re-parameterised for other countries. OpenABM-Covid19 can evaluate non-pharmaceutical interventions, including both manual and digital contact tracing, and vaccination programmes. It can simulate a population of 1 million people in seconds per day, allowing parameter sweeps and formal statistical model-based inference. The code is open-source and has been developed by teams both inside and outside academia, with an emphasis on formal testing, documentation, modularity and transparency. A key feature of OpenABM-Covid19 are its Python and R interfaces, which has allowed scientists and policymakers to simulate dynamic packages of interventions and help compare options to suppress the COVID-19 epidemic.

Virological Characterization of Critically Ill Patients With COVID-19 in the United Kingdom: Interactions of Viral Load, Antibody Status, and B.1.1.7 Infection.

BACKGROUND: Convalescent plasma containing neutralizing antibody to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is under investigation for coronavirus disease 2019 (COVID-19) treatment. We report diverse virological characteristics of UK intensive care patients enrolled in the Immunoglobulin Domain of the REMAP-CAP randomized controlled trial that potentially influence treatment outcomes. METHODS: SARS-CoV-2 RNA in nasopharyngeal swabs collected pretreatment was quantified by PCR. Antibody status was determined by spike-protein ELISA. B.1.1.7 was differentiated from other SARS-CoV-2 strains using allele-specific probes or restriction site polymorphism (SfcI) targeting D1118H. RESULTS: Of 1274 subjects, 90% were PCR positive with viral loads 118-1.7 × 1011IU/mL. Median viral loads were 40-fold higher in those IgG seronegative (n = 354; 28%) compared to seropositives (n = 939; 72%). Frequencies of B.1.1.7 increased from <1% in November 2020 to 82% of subjects in January 2021. Seronegative individuals with wild-type SARS-CoV-2 had significantly higher viral loads than seropositives (medians 5.8 × 106 and 2.0 × 105 IU/mL, respectively; P = 2 × 10-15). CONCLUSIONS: High viral loads in seropositive B.1.1.7-infected subjects and resistance to seroconversion indicate less effective clearance by innate and adaptive immune responses. SARS-CoV-2 strain, viral loads, and antibody status define subgroups for analysis of treatment efficacy.

Simultaneous Viral Whole-Genome Sequencing and Differential Expression Profiling in Respiratory Syncytial Virus Infection of Infants.

Targeted metagenomics using strand-specific libraries with target enrichment is a sensitive, generalized approach to pathogen sequencing and transcriptome profiling. Using this method, we recovered 13 (76%) complete human respiratory syncytial virus (RSV) genomes from 17 clinical respiratory samples, reconstructed the phylogeny of the infecting viruses, and detected differential gene expression between 2 RSV subgroups, specifically, a lower expression of the P gene and a higher expression of the M2 gene in RSV-A than in RSV-B. This methodology can help to relate viral genetics to clinical phenotype and facilitate ongoing population-level RSV surveillance and vaccine development. Clinical Trials Registration. NCT03627572 and NCT03756766.

A Comprehensive Genomics Solution for HIV Surveillance and Clinical Monitoring in Low-Income Settings.

Viral genetic sequencing can be used to monitor the spread of HIV drug resistance, identify appropriate antiretroviral regimes, and characterize transmission dynamics. Despite decreasing costs, next-generation sequencing (NGS) is still prohibitively costly for routine use in generalized HIV epidemics in low- and middle-income countries. Here, we present veSEQ-HIV, a high-throughput, cost-effective NGS sequencing method and computational pipeline tailored specifically to HIV, which can be performed using leftover blood drawn for routine CD4 cell count testing. This method overcomes several major technical challenges that have prevented HIV sequencing from being used routinely in public health efforts; it is fast, robust, and cost-efficient, and generates full genomic sequences of diverse strains of HIV without bias. The complete veSEQ-HIV pipeline provides viral load estimates and quantitative summaries of drug resistance mutations; it also exploits information on within-host viral diversity to construct directed transmission networks. We evaluated the method's performance using 1,620 plasma samples collected from individuals attending 10 large urban clinics in Zambia as part of the HPTN 071-2 study (PopART Phylogenetics). Whole HIV genomes were recovered from 91% of samples with a viral load of >1,000 copies/ml. The cost of the assay (30 GBP per sample) compares favorably with existing VL and HIV genotyping tests, proving an affordable option for combining HIV clinical monitoring with molecular epidemiology and drug resistance surveillance in low-income settings.

Performance of a high-throughput next-generation sequencing method for analysis of HIV drug resistance and viral load.

OBJECTIVES: To evaluate the performance of a high-throughput research assay for HIV drug resistance testing based on whole genome next-generation sequencing (NGS) that also quantifies HIV viral load. METHODS: Plasma samples (n = 145) were obtained from HIV-positive MSM (HPTN 078). Samples were analysed using clinical assays (the ViroSeq HIV-1 Genotyping System and the Abbott RealTime HIV-1 Viral Load assay) and a research assay based on whole-genome NGS (veSEQ-HIV). RESULTS: HIV protease and reverse transcriptase sequences (n = 142) and integrase sequences (n = 138) were obtained using ViroSeq. Sequences from all three regions were obtained for 100 (70.4%) of the 142 samples using veSEQ-HIV; results were obtained more frequently for samples with higher viral loads (93.5% for 93 samples with >5000 copies/mL; 50.0% for 26 samples with 1000-5000 copies/mL; 0% for 23 samples with <1000 copies/mL). For samples with results from both methods, drug resistance mutations (DRMs) were detected in 33 samples using ViroSeq and 42 samples using veSEQ-HIV (detection threshold: 5.0%). Overall, 146 major DRMs were detected; 107 were detected by both methods, 37 were detected by veSEQ-HIV only (frequency range: 5.0%-30.6%) and two were detected by ViroSeq only. HIV viral loads estimated by veSEQ-HIV strongly correlated with results from the Abbott RealTime Viral Load assay (R2 = 0.85; n = 142). CONCLUSIONS: The NGS-based veSEQ-HIV method provided results for most samples with higher viral loads, was accurate for detecting major DRMs, and detected mutations at lower levels compared with a method based on population sequencing. The veSEQ-HIV method also provided HIV viral load data.

High prevalence of integrase mutation L74I in West African HIV-1 subtypes prior to integrase inhibitor treatment.

OBJECTIVES: HIV-1 integrase inhibitors are recommended as first-line therapy by WHO, though efficacy and resistance data for non-B subtypes are limited. Two recent trials have identified the integrase L74I mutation to be associated with integrase inhibitor treatment failure in HIV-1 non-B subtypes. We sought to define the prevalence of integrase resistance mutations, including L74I, in West Africa. METHODS: We studied a Nigerian cohort of recipients prior to and during receipt of second-line PI-based therapy, who were integrase inhibitor-naive. Illumina next-generation sequencing with target enrichment was used on stored plasma samples. Drug resistance was interpreted using the Stanford Resistance Database and the IAS-USA 2019 mutation lists. RESULTS: Of 115 individuals, 59.1% harboured CRF02_AG HIV-1 and 40.9% harboured subtype G HIV-1. Four participants had major IAS-USA integrase resistance-associated mutations detected at low levels (2%-5% frequency). Two had Q148K minority variants and two had R263K (one of whom also had L74I). L74I was detected in plasma samples at over 2% frequency in 40% (46/115). Twelve (26.1%) had low-level minority variants of between 2% and 20% of the viral population sampled. The remaining 34 (73.9%) had L74I present at >20% frequency. L74I was more common among those with subtype G infection (55.3%, 26/47) than those with CRF02_AG infection (29.4%, 20/68) (P = 0.005). CONCLUSIONS: HIV-1 subtypes circulating in West Africa appear to have very low prevalence of major integrase mutations, but significant prevalence of L74I. A combination of in vitro and clinical studies is warranted to understand the potential implications.

Resistance analysis of genotype 3 hepatitis C virus indicates subtypes inherently resistant to nonstructural protein 5A inhibitors.

Hepatitis C virus (HCV) genotype (gt) 3 is highly prevalent globally, with non-gt3a subtypes common in Southeast Asia. Resistance-associated substitutions (RASs) have been shown to play a role in treatment failure. However, the role of RASs in gt3 is not well understood. We report the prevalence of RASs in a cohort of direct-acting antiviral treatment-naive, gt3-infected patients, including those with rarer subtypes, and evaluate the effect of these RASs on direct-acting antivirals in vitro. Baseline samples from 496 gt3 patients enrolled in the BOSON clinical trial were analyzed by next-generation sequencing after probe-based enrichment for HCV. Whole viral genomes were analyzed for the presence of RASs to approved direct-acting antivirals. The resistance phenotype of RASs in combination with daclatasvir, velpatasvir, pibrentasvir, elbasvir, and sofosbuvir was measured using the S52 ΔN gt3a replicon model. The nonstructural protein 5A A30K and Y93H substitutions were the most common at 8.9% (n = 44) and 12.3% (n = 61), respectively, and showed a 10-fold and 11-fold increase in 50% effect concentration for daclatasvir compared to the unmodified replicon. Paired RASs (A30K + L31M and A30K + Y93H) were identified in 18 patients (9 of each pair); these combinations were shown to be highly resistant to daclatasvir, velpatasvir, elbasvir, and pibrentasvir. The A30K + L31M combination was found in all gt3b and gt3g samples. Conclusion: Our study reveals high frequencies of RASs to nonstructural protein 5A inhibitors in gt3 HCV; the paired A30K + L31M substitutions occur in all patients with gt3b and gt3g virus, and in vitro analysis suggests that these subtypes may be inherently resistant to all approved nonstructural protein 5A inhibitors for gt3 HCV. (Hepatology 2018).