Numerous experimental evolution studies have suggested that adaptation rate of microbial populations evolving in stable environments decline over time. Despite the generality of this phenomenon across different domains of life, the timing and magnitude of decline in adaptation can vary greatly based on the idiosyncrasies of the biological system. To investigate the characteristics of adaptation deceleration in a fast-evolving virus, we propagated HIV-1 in two human T-cell lines (MT-2 and MT-4) for approximately 4.8 years and tracked its genome evolution through next-generation sequencing. The curated sequencing data covering the whole-genome can be accessed and explored via LTEEviz, an interactive web application. Time-resolved sequencing data uncovered that despite constant fixation rate of 0.085 (MT-2) and 0.042 (MT-4) mutations per generation, the fixation kinetics of adaptive mutations changed considerably over time. The rate of fixation of adaptive parallel mutations decreased by 44% per 300 generations, while their conferred fitness gain decreased by 27% (MT-2) and 18% (MT-4) per every added adaptive mutation in their genetic background. The early and substantial deceleration of adaptation in our HIV-1 populations can, at least in part, be explained by diminishing gains of adaptive mutations. Furthermore, we identified genomic patterns consistent with a hard selective sweep that occurred in one population later in the experiment. Together, our results confirm that HIV-1 genomic evolution is characterized by a swift and substantial deceleration of adaptation, while also revealing that episodes of positive selection can occur beyond the initial adaptive phase.

Background:People with human immunodeficiency virus type 1 (HIV-1; PWH) on antiretroviral therapy (ART) and suppressed viremia usually experience a decrease in HIV DNA over time, but about 25% experience an increase. Some also experience intermittent viremia. The reasons for and potential clinical implications of increases in the HIV reservoir remain unclear and a major concern. Methods:In this study, we longitudinally characterized the proviral landscape in 4 distinct groups of PWH (n = 40) successfully treated with ART over 10.4 years without any viral failure, presenting either an increase or decrease of total HIV DNA levels and experiencing or not experiencing intermittent viremia, by intact proviral DNA assay and near-full-length HIV proviral sequencing in bulk and on the single proviral level. Results:A decrease in intact proviruses was observed in all groups, independent of total HIV DNA level dynamics and viral load kinetics by both intact proviral DNA assay and single proviral sequencing. Genetic distances and diversities of individuals' proviral sequences did not increase over time in any group. Although new drug resistance mutations were occasionally observed in proviral DNA, numbers did not differ significantly between the groups. Conclusions:Our results show that the increase in HIV DNA levels is driven by an increase in defective proviruses, also in PWH experiencing intermittent viremia. Furthermore, we did not see evidence of evolution of the HIV reservoir, regardless of HIV reservoir size dynamics and viral load kinetics over a follow-up period of 10 years with ART. Nevertheless, PWH with intermittent viremia should be monitored frequently.

Numerous HIV related outcomes can be determined on the viral genome, for example, resistance associated mutations, population transmission dynamics, viral heritability traits, or time since infection. Viral sequences of people with HIV (PWH) are therefore essential for therapeutic and research purposes. While in the first three decades of the HIV pandemic viral genomes were mainly sequenced using Sanger sequencing, the last decade has seen a shift towards next-generation sequencing (NGS) as the preferred method. NGS can achieve near full length genome sequence coverage and simultaneously, it accurately encapsulates the within-host diversity by characterizing HIV subpopulations. NGS opens new avenues for HIV research, but it also presents challenges concerning data management and analysis. We therefore set up the Swiss HIV Cohort Study Viral NGS Database (SHCND) to address key issues in the handling of NGS data including high loads of raw- and processed NGS data, data storage solutions, downstream application of sophisticated bioinformatic tools, high-performance computing resources, and reproducibility. The database is nested within the Swiss HIV Cohort Study (SHCS) and the Zurich Primary HIV Infection Cohort Study (ZPHI), which together enrolled 21,876 PWH since 1988 and include a biobank dating back to the early nineties. Since its initiation in 2018, the SHCND accumulated NGS sequences (plasma and proviral origin) of 5,178 unique PWH. We here describe the design, set-up, and use of this NGS database. Overall, the SHCND has contributed to several research projects on HIV pathogenesis, treatment, drug resistance, and molecular epidemiology, and has thereby become a central part of HIV-genomics research in Switzerland.

The intact proviral DNA assay (IPDA) has become a gold standard for HIV-1 reservoir quantification in HIV-1 latency research, as well as in the evaluation of HIV-1 cure strategies. In this work, we adjusted the IPDA assay to a chip-based digital PCR (dPCR) system, established the use of CCR5 as a different reference gene, and evaluated the performance on cell lines, clinical samples from people with HIV (PWH) off and on antiretroviral therapy (ART), and different HIV-1 subtypes. Our adapted IPDA performs well on the chip-based dPCR system with no false positive intact HIV-1 in negative controls and with less hands-on time compared to droplet-based dPCR systems. Undetectable intact HIV-1 DNA is common in clinical samples on ART and correlated with total HIV-1 reservoir size and sample input concentration for HIV-1 subtype B samples. Performance on non-B subtypes varies depending on the subtype and should be improved with subtype specific primer/probe combinations.

Background:Bacillus Calmette-Guérin (BCG) vaccination, used against tuberculosis, is recognized for its immunomodulatory properties, a phenomenon referred to as "trained immunity." Given these effects, there is increasing interest in evaluating its safety and impact on immune function in people living with HIV-1 (PWH). Historically, BCG was contraindicated in PWH due to safety concerns in immunocompromised individuals. This study aims to assess both the safety of BCG in PWH and its effects on the HIV-1 latent reservoir size. Methods:This Phase IIA randomized, double-blind, placebo-controlled, single-center trial enrolled 60 PWH with a suppressed viral load and CD4 T-cell count >350/μL. Participants were randomized in a stepped-wedge design into equal groups for early or late BCG vaccination. Each participant received a single intradermal dose of BCG vaccine followed by a placebo 3 months later, or vice versa. The HIV-1 latent reservoir was quantified at 3-month intervals to day 270. The primary endpoint was the HIV-1 reservoir size 6 months postvaccination, with secondary endpoints including safety outcomes. Results:No significant differences were found in intact proviral HIV-1 DNA levels at 6 months compared to baseline. Local reactions occurred in 96% of participants, leading to scarring in 73%. No systemic infections or serious BCG-related adverse events were observed. Conclusions:BCG vaccination is safe in PWH, but local skin reactions including scarring are common. There was no significant effect on the HIV-1 reservoir. These findings provide valuable insights into the safety profile of BCG vaccination in PWH, emphasizing its potential for broader immunological studies.

HIV-1 RNA genetic diversity predicts time since infection, which is important for clinical care and research. It is unclear, however, whether proviral DNA genetic diversity sampled under suppressive antiretroviral therapy can be used for this purpose. We tested whether proviral genetic diversity from next-generation sequencing predicts time since infection and recency in 221 people with HIV-1 with known infection time. Proviral diversity was significantly associated with time since infection (P < 5x10-7, R2 up to 25%) and predictive of treatment initiation during recent infection (area under the curve-receiver operating characteristic up to 0.85). This shows the utility of proviral genetic diversity as a proxy for time since infection.

The pathogenesis of HIV-1 infection is governed by a highly dynamic, time-dependent interaction between the host and the viral genome. In this study, we developed a novel systematic approach to assess the host-virus interaction, using average pairwise viral diversity as a proxy for time since infection, and applied this method to nearly whole viral genome sequences (n = 4,464), human leukocyte antigen (HLA) genotyping data (n = 1,044), and viral RNA load (VL) measurements during the untreated chronic phase (n = 829) of Swiss HIV Cohort Study participants. Our systematic genome-wide screen revealed for 98 HLA/viral-variant pairs a signature of immune-driven selection in the form of an HLA-dependent effect of infection time on the presence of HIV amino acid variants. Of these pairs, 12 were found to have an effect on VL. Furthermore, 28/58 pairs were validated by time-to-event analyses and 48/92 by computational HLA-epitope predictions. Our diversity-based approach allows a powerful and systematic investigation of the interaction between the virus and cellular immunity, revealing a notable subset of such interaction effects. From an evolutionary perspective, these observations underscore the complexity of HLA-mediated selection pressures on the virus that shape viral evolution and pathogenesis.

In people with HIV-1 (PWH), Mycobacterium tuberculosis (MTB) infection poses a significant threat. While active tuberculosis (TB) accelerates immunodeficiency, the interaction between MTB and HIV-1 during asymptomatic phases remains unclear. Analysis of peripheral blood mononuclear cells (PBMC) transcriptomic profiles in PWH, with and without controlled viral loads, revealed distinct clustering in MTB-infected individuals. Functional annotation identified alterations in IL-6, TNF, and KRAS pathways. Notably, MTB-related genes displayed an inverse correlation with HIV-1 viremia, at both individual and signature score levels. These findings suggest that MTB infection in PWH induces a shift in immune system activation, inversely related to HIV-1 viral load. These results may explain the observed enhanced antiretroviral control in MTB-infected PWH. This study highlights the complex interplay between MTB and HIV-1, emphasizing the importance of understanding their interaction for managing co-infections in this population.

People with HIV may report neurocognitive complaints, with or without associated neurocognitive impairment, varying between individuals and populations. While the HIV genome could play a major role, large systematic viral genome-wide screens to date are lacking. The Swiss HIV Cohort Study biannually enquires neurocognitive complaints. We quantified broad-sense heritability estimates using partial 'pol' sequences from the Swiss HIV Cohort Study resistance database and performed a viral near full-length genome-wide association study for the longitudinal area under the curve of neurocognitive complaints. We performed all analysis (i) restricted to HIV Subtype B and (ii) including all HIV subtypes. From 8547 people with HIV with neurocognitive complaints, we obtained 6966 partial 'pol' sequences and 2334 near full-length HIV sequences. Broad-sense heritability estimates for presence of memory loss complaints ranged between 1% and 17% (Subtype B restricted 1-22%) and increased with the stringency of the phylogenetic distance thresholds. The genome-wide association study revealed one amino acid (Env L641E), after adjusting for multiple testing, positively associated with memory loss complaints (P = 4.3 * 10-6). Other identified mutations, while insignificant after adjusting for multiple testing, were reported in other smaller studies (Tat T64N, Env *291S). We present the first HIV genome-wide association study analysis of neurocognitive complaints and report a first estimate for the heritability of neurocognitive complaints through HIV. Moreover, we could identify one mutation significantly associated with the presence of memory loss complaints. Our findings indicate that neurocognitive complaints are polygenetic and highlight advantages of a whole genome approach for pathogenicity determination. Zeeb et al. report that self-reported neurocognitive complaints in people with HIV are in part heritable through the HIV genome, but in a viral genome-wide associations study, they found only one significant amino acid substitution. They conclude that self-reported neurocognitive complaints are a highly polygenic phenotype. Graphical Abstract

Abstract Background Starting combination antiretroviral therapy (cART) during primary human immunodeficiency virus type 1 (HIV-1) infection results in a smaller HIV-1 latent reservoir, reduced immune activation, and less viral diversity compared to starting cART during chronic infection. We report results of a 4-year study designed to determine whether these properties would allow sustained virological suppression after simplification of cART to dolutegravir (DTG) monotherapy. Methods EARLY-SIMPLIFIED is a randomized, open-label, noninferiority trial. People with HIV (PWH) who started cART <180 days after a documented primary HIV-1 infection with suppressed viral load were randomized (2:1) to DTG monotherapy with 50 mg daily or continuation of cART. The primary endpoints were the proportion of PWH with viral failure at 48, 96, 144, and 192 weeks; noninferiority margin was 10%. After 96 weeks, randomization was lifted and patients were permitted to switch treatment groups as desired. Results Of 101 PWH randomized, 68 were assigned to DTG monotherapy and 33 to cART. At week 96 in the per-protocol population, 64/64 (100%) showed virological response in the DTG monotherapy group versus 30/30 (100%) in the cART group (difference, 0.00%; upper bound of 95% confidence interval 6.22%). This demonstrated noninferiority of DTG monotherapy at the prespecified level. At week 192, the study end, no virological failure occurred in either group during 13 308 and 4897 person weeks of follow-up for the DTG monotherapy (n = 80) and cART groups, respectively. Conclusions This trial suggests that early cART initiation during primary HIV infection allows sustained virological suppression after switching to DTG monotherapy.

Infection with Mycobacterium tuberculosis (MTB) remains one of the most important opportunistic infections in people with HIV-1 (PWH). While active Tuberculosis (TB) leads to rapid progression of immunodeficiency in PWH, the interaction between MTB and HIV-1 during the asymptomatic phase of both infections remains poorly understood. In a cohort of individuals with HIV (PWH) with and without suppressed HIV-1 viral load, the transcriptomic profiles of peripheral blood mononuclear cells (PBMC) clustered in individuals infected with Mycobacterium tuberculosis (MTB) compared to carefully matched controls. Subsequent functional annotation analysis disclosed alterations in the IL-6, TNF, and KRAS pathways. Notably, MTB-associated genes demonstrated an inverse correlation with HIV-1 viremia, evident at both on individual gene level and when employed as a gene score. In sum, our data show that MTB infection in PWH is associated with a shift in the activation state of the immune system, displaying an inverse relationship with HIV-1 viral load. These results could provide an explanation for the observed increased antiretroviral control associated with MTB infection in PWH.

Background Genotypic resistance testing (GRT) is routinely performed upon diagnosis of HIV-1 infection or during virological failure using plasma viral RNA. An alternative source for GRT could be cellular HIV-1 DNA. Objectives A substantial number of participants in the Swiss HIV Cohort Study (SHCS) never received GRT. We applied a method that enables access to the near full-length proviral HIV-1 genome without requiring detectable viraemia. Methods Nine hundred and sixty-two PBMC specimens were received. Our two-step nested PCR protocol was applied to generate two overlapping long-range amplicons of the HIV-1 genome, sequenced by next-generation sequencing (NGS) and analysed by MinVar, a pipeline to detect drug resistance mutations (DRMs). Results Six hundred and eighty-one (70.8%) of the samples were successfully amplified, sequenced and analysed by MinVar. Only partial information of the pol gene was contained in 82/681 (12%), probably due to naturally occurring deletions in the proviral sequence. All common HIV-1 subtypes were successfully sequenced. We detected at least one major DRM at high frequency (& GE;15%) in 331/599 (55.3%) individuals. Excluding APOBEC-signature (G-to-A mutation) DRMs, 145/599 (24.2%) individuals carried at least one major DRM. RT-inhibitor DRMs were most prevalent. The experienced time on ART was significantly longer in DRM carriers (P = 0.001) independent of inclusion or exclusion of APOBEC-signature DRMs. Conclusions We successfully applied a reliable and efficient method to analyse near full-length HIV-1 proviral DNA and investigated DRMs in individuals with undetectable or low viraemia. Additionally, our data underscore the need for new computational tools to exclude APOBEC-related hypermutated NGS sequence reads for reporting DRMs.

Human immunodeficiency virus type 1 (HIV-1) infection is treated with antiretroviral therapy (ART), usually consisting of 2-3 different drugs, referred to as combination ART (cART). Our recent randomized clinical trial comparing a switch to dolutegravir monotherapy with continuation of cART in early-treated individuals demonstrated sustained virological suppression over 48 weeks. Here, we characterize the longitudinal landscape of the HIV-1 reservoir in these participants, with particular attention to potential differences between treatment groups regarding evidence of evolution as a proxy for low-level replication. Near full-length HIV-1 proviral polymerase chain reaction and next-generation sequencing was applied to longitudinal peripheral blood mononuclear cell samples to assess proviral evolution and the potential emergence of drug resistance mutations (DRMs). Neither an increase in genetic distance nor diversity over time was detected in participants of both treatment groups. Single proviral analysis showed high proportions of defective proviruses and low DRM numbers. No evidence for evolution during dolutegravir monotherapy was found in these early-treated individuals.

Background HIV-1 replication capacity (RC) of transmitted/founder viruses may influence the further course of HIV-1 infection. Methods RCs of 355 whole-genome primary HIV-1 isolates derived from samples acquired during acute and recent primary HIV-1 infection (PHI) were determined using a novel high-throughput infection assay in primary cells. The RCs were used to elucidate potential factors that could be associated with RC during PHI. Results Increased RC was found to be associated with increased set point viral load (VL), and significant differences in RCs among 13 different HIV-1 subtypes were discerned. Notably, we observed an increase in RCs for primary HIV-1 isolates of HIV-1 subtype B over a 17-year period. Associations were not observed between RC and CD4 count at sample date of RC measurement, CD4 recovery after initiation of antiretroviral treatment, CD4 decline in untreated individuals, and acute retroviral syndrome severity scores. Conclusions These findings highlight that RCs of primary HIV-1 isolates acquired during the acute and recent phase of infection are more associated with viral factors, that is set point VL, than with host factors. Furthermore, we observed a temporal increase in RC for HIV-1 subtype B viruses over a period of 17 years. Clinical Trials Registration NCT00537966.

HIV-1 reservoir size and dynamics are promising parameters to ensure the safe prescription of simplified maintenance antiretroviral therapy in chronically HIV-1 infected patients. In the SIMPL’HIV trial, HIV-1 DNA was quantified in peripheral blood mononuclear cells obtained at baseline and week 48 to investigate changes over time and evidence of a predictive relationship to maintain HIV-1 RNA <20 copies/ml. Measurements were available for 175 patients, with no differences observed between treatment strategies. Findings showed that baseline HIV-1 DNA was lower in those with durable HIV-1 RNA <20 copies/ml compared with patients with incomplete viral suppression over 48 weeks.

HIV-1 replication capacity is an important characteristic to understand the replication competence of single variants or virus populations. It can further aid in the understanding of HIV-1 pathogenicity, disease progression, and drug resistance mutations. To effectively study RC, many assays have been established. However, there is still demand for a high throughput replication capacity assay using primary cells which is robust and reproducible. In this study, we established such an assay and validated it using 346 primary HIV-1 isolates from patients enrolled in the Zurich Primary HIV Infection study (ZPHI) and two control viruses, HIV-1 JR-CSFWT and HIV-1 JR-CSFK65R_M184V. Replication capacity was determined by measuring the viral growth on PBMCs over 10 days by longitudinally transferring cell culture supernatant to TZM-bl reporter cells. By utilizing the TZM-bl luciferase reporter assay, we determined replication capacity by measuring viral infectivity. The simplicity of the experimental setup allowed for all 346 primary HIV-1 isolates to be replicated at one time. Although the infectious input dose for each virus was normalized, a broad range of replication capacity values over 4 logs was observed. The approach was confirmed by two repeated experiments and we demonstrated that the reproducibility of the replication capacity values is statistically comparable between the two separate experiments. In summary, these results endorse our high throughput replication capacity assay as reproducible and robust and can be utilized for large scale HIV-1 replication capacity experiments in primary cells.

HIV-1 is capable of integrating its genome into that of its host cell. We examined the influence of the activation state of CD4+ T cells, the effect of antiretroviral therapy (ART), and the clinical stage of HIV-1 infection on HIV-1 integration site features and selection. HIV-1 integration sites were sequenced from longitudinally sampled resting and activated CD4+ T cells from 12 HIV-1–infected individuals. In total, 589 unique HIV-1 integration sites were analyzed: 147, 391, and 51 during primary, chronic, and late presentation of HIV-1 infection, respectively. As early as during primary HIV-1 infection and independent of the activation state of CD4+ T cells collected on and off ART, HIV-1 integration sites were preferentially detected in recurrent integration genes, genes associated with clonal expansion of latently HIV-1–infected CD4+ T cells, cancer-related genes, and highly expressed genes. The preference for cancer-related genes was more pronounced at late stages of HIV-1 infection. Host genomic features of HIV-1 integration site selection remained stable during HIV-1 infection in both resting and activated CD4+ T cells. In summary, characteristic HIV-1 integration site features are preestablished as early as during primary HIV-1 infection and are found in both resting and activated CD4+ T cells.

A proportion of EBV+ B cells are infectible by CXCR4-tropic HIV-1 in vitro and in vivo, and these HIV-infected B cells can transmit infection to previously virus-naïve humanized mice. HIV and EBV are human pathogens that cause a considerable burden to worldwide health. In combination, these viruses are linked to AIDS-associated lymphomas. We found that EBV, which transforms B cells, renders them susceptible to HIV-1 infection in a CXCR4 and CD4-dependent manner in vitro and that CXCR4-tropic HIV-1 integrates into the genome of these B cells with the same molecular profile as in autologous CD4+ T cells. In addition, we established a humanized mouse model to investigate the in vivo interactions of EBV and HIV-1 upon coinfection. The respective mice that reconstitute human immune system components upon transplantation with CD34+ human hematopoietic progenitor cells could recapitulate aspects of EBV and HIV immunobiology observed in dual-infected patients. Upon coinfection of humanized mice, EBV/HIV dual-infected B cells could be detected, but were susceptible to CD8+ T-cell–mediated immune control.

Background: The primary hurdle for the eradication of HIV-1 is the establishment of a latent viral reservoir early after primary infection. Here, we investigated the potential influence of human genetic variation on the HIV-1 reservoir size and its decay rate during suppressive antiretroviral treatment. Setting: Genome-wide association study and exome sequencing study to look for host genetic determinants of HIV-1 reservoir measurements in patients enrolled in the Swiss HIV Cohort Study, a nation-wide prospective observational study. Methods: We measured total HIV-1 DNA in peripheral blood mononuclear cells from study participants, as a proxy for the reservoir size at 3 time points over a median of 5.4 years, and searched for associations between human genetic variation and 2 phenotypic readouts: the reservoir size at the first time point and its decay rate over the study period. We assessed the contribution of common genetic variants using genome-wide genotyping data from 797 patients with European ancestry enrolled in the Swiss HIV Cohort Study and searched for a potential impact of rare variants and exonic copy number variants using exome sequencing data generated in a subset of 194 study participants. Results: Genome-wide and exome-wide analyses did not reveal any significant association with the size of the HIV-1 reservoir or its decay rate on suppressive antiretroviral treatment. Conclusions: Our results point to a limited influence of human genetics on the size of the HIV-1 reservoir and its long-term dynamics in successfully treated individuals.

Background The persistence of the latent HIV-1 reservoir is a major obstacle to curing HIV-1 infection. HIV-1 integrates into the cellular genome and some targeted genomic loci are frequently detected in clonally expanded latently HIV-1 infected cells, for instance, the gene BTB domain and CNC homology 2 (BACH2) . Methods We investigated HIV-1 promoter activity after integration into specific sites in BACH2 in Jurkat T-cells. The HIV-1-based vector LTatCL[M] contains two fluorophores: (1) Cerulean, which reports the activity of the HIV-1 promoter and (2) mCherry driven by a constitutive promotor and flanked by genetic insulators. This vector was inserted into introns 2 and 5 of BACH2 of Jurkat T-cells via CRISPR/Cas9 technology in the same and convergent transcriptional orientation of BACH2 , and into the genomic safe harbour AAVS1. Single cell clones representing active (Cerulean + /mCherry + ) and inactive (Cerulean – /mCherry + ) HIV-1 promoters were characterised. Results Upon targeted integration of the 5.3 kb vector LTatCL[M] into BACH2 , the HIV-1 promoter was gradually silenced as reflected by the decrease in Cerulean expression over a period of 162 days. Silenced HIV-1 promoters could be reactivated by TNF-α and Romidepsin. This observation was independent of the targeted intron and the transcriptional orientation. BACH2 mRNA and protein expression was not impaired by mono-allelic integration of LTatCL[M]. Conclusion Successful targeted integration of the HIV-1-based vector LTatCL[M] allows longitudinal analyses of HIV-1 promoter activity.