Although combination antiretroviral therapy (cART) initiated in the acute phase of HIV-1 infection may prevent expansion of the latent reservoir, its benefits remain controversial. In the current study, HIV-1 RNA transcription patterns in peripheral blood mononuclear cells (PBMC) were monitored during acute cART to assess the effect of early treatment on cellular viral reservoirs.

The HIV-1 envelope trimer adopts a quaternary conformation that effectively shields neutralization-sensitive domains and thus represents a major obstacle for natural and vaccine-elicited antibody responses. By using a structure-function analysis based on a specifically devised mathematical model, we demonstrate in this study that protection from neutralization is enforced by intersubunit contact between the variable loops 1 and 2 (V1V2) and domains of neighboring gp120 subunits in the trimer encompassing the V3 loop. Our data are consistent with an interaction of the V1V2 and V3 loop at the spike apex as proposed by cryoelectron tomography experiments. By defining the orientation of the V1V2 loop within the trimer toward the neighboring gp120 subunit's V3 loop, our data close an important gap in the understanding of the architecture of the trimeric spike. Knowledge on how the V1V2 barrier functions in the context of the trimer to mask conserved epitopes on gp120 may aid future vaccine design.

Interference with virus entry is known to be the principle mechanism of HIV neutralization by antibodies, including 2F5 and 4E10, which bind to the membrane-proximal external region (MPER) of the gp41 envelope protein. However, to date, the precise molecular events underlying neutralization by MPER-specific antibodies remain incompletely understood. In this study, we investigated the capacity of these antibodies to irrevocably sterilize HIV virions. Long-term effects of antibodies on virions can differ, rendering neutralization either reversible or irreversible. MPER-specific antibodies irreversibly neutralize virions, and this capacity is associated with induction of gp120 shedding. Both processes have similar thermodynamic properties and slow kinetics requiring several hours. Antibodies directed to the CD4 binding site, V3 loop, and the MPER can induce gp120 shedding, and shedding activity is detected with high frequency in plasma from patients infected with divergent genetic HIV-1 subtypes. Importantly, as we show in this study, induction of gp120 shedding is closely associated with MPER antibody inhibition, constituting either a primary event leading to virion neutralization or representing an immediate consequence thereof, and thus needs to be factored into the mechanistic processes underlying their activity.

The time passed since the infection of a human immunodeficiency virus (HIV)-infected individual (the age of infection) is an important but often only poorly known quantity. We assessed whether the fraction of ambiguous nucleotides obtained from bulk sequencing as done for genotypic resistance testing can serve as a proxy of this parameter.

Mucosal-associated invariant T (MAIT) cells are an abundant antibacterial innate-like lymphocyte population. There are conflicting reports as to their fate in HIV infection. The objective of this study was to determine whether MAIT cells are truly depleted in HIV infection. In this case-control study of HIV-positive patients and healthy controls, quantitative real-time polymerase chain reaction was used to assess the abundance of messenger RNA (mRNA) and genomic DNA (gDNA) encoding the canonical MAIT cell T cell receptor (Vα7.2-Jα33). Comparison was made with flow cytometry. Significant depletion of both Vα7.2-Jα33 mRNA and gDNA was seen in HIV infection. Depletion of Vα7.2+CD161++ T cells was confirmed by flow cytometry. In HIV infection, the abundance of Vα7.2-Jα33 mRNA correlated most strongly with the frequency of Vα7.2+CD161++ cells. No increase was observed in the frequency of Vα7.2+CD161- cells among CD3+CD4- lymphocytes. MAIT cells are depleted from blood in HIV infection as confirmed by independent assays. Significant accumulation of a CD161- MAIT cell population is unlikely. Molecular approaches represent a suitable alternative to flow cytometry-based assays for tracking of MAIT cells in HIV and other settings.

The humoral immune response plays a critical role in controlling infection, and the rapid adaptation to a broad range of pathogens depends on a highly diverse antibody repertoire. The advent of high-throughput sequencing technologies in the past decade has enabled insights into this immense diversity. However, not only the variable, but also the constant region of antibodies determines their in vivo activity. Antibody isotypes differ in effector functions and are thought to play a defining role in elicitation of immune responses, both in natural infection and in vaccination. We have developed an Illumina MiSeq high-throughput sequencing protocol that allows determination of the human IgG subtype alongside sequencing full-length antibody variable heavy chain regions. We thereby took advantage of the Illumina procedure containing two additional short reads as identifiers. By performing paired-end sequencing of the variable regions and customizing one of the identifier sequences to distinguish IgG subtypes, IgG transcripts with linked information of variable regions and IgG subtype can be retrieved. We applied our new method to the analysis of the IgG variable region repertoire from PBMC of an HIV-1 infected individual confirmed to have serum antibody reactivity to the Membrane Proximal External Region (MPER) of gp41. We found that IgG3 subtype frequencies in the memory B cell compartment increased after halted treatment and coincided with increased plasma antibody reactivity against the MPER domain. The sequencing strategy we developed is not restricted to analysis of IgG. It can be adopted for any Ig subtyping and beyond that for any research question where phasing of distant regions on the same amplicon is needed.

Phylogenetic clustering approaches can elucidate HIV transmission dynamics. Comparisons across countries are essential for evaluating public health policies. Here, we used a standardised approach to compare the UK HIV Drug Resistance Database and the Swiss HIV Cohort Study while maintaining data-protection requirements. Clusters were identified in subtype A1, B and C pol phylogenies. We generated degree distributions for each risk group and compared distributions between countries using Kolmogorov-Smirnov (KS) tests, Degree Distribution Quantification and Comparison (DDQC) and bootstrapping. We used logistic regression to predict cluster membership based on country, sampling date, risk group, ethnicity and sex. We analysed >8,000 Swiss and >30,000 UK subtype B sequences. At 4.5% genetic distance, the UK was more clustered and MSM and heterosexual degree distributions differed significantly by the KS test. The KS test is sensitive to variation in network scale, and jackknifing the UK MSM dataset to the size of the Swiss dataset removed the difference. Only heterosexuals varied based on the DDQC, due to UK male heterosexuals who clustered exclusively with MSM. Their removal eliminated this difference. In conclusion, the UK and Swiss HIV epidemics have similar underlying dynamics and observed differences in clustering are mainly due to different population sizes.

Long-term benefits of combination antiretroviral therapy (cART) initiation during primary HIV-1 infection are debated.

Next-generation sequencing (NGS) is a valuable tool for the detection and quantification of HIV-1 variants in vivo. However, these technologies require detailed characterization and control of artificially induced errors to be applicable for accurate haplotype reconstruction. To investigate the occurrence of substitutions, insertions, and deletions at the individual steps of RT-PCR and NGS, 454 pyrosequencing was performed on amplified and non-amplified HIV-1 genomes. Artificial recombination was explored by mixing five different HIV-1 clonal strains (5-virus-mix) and applying different RT-PCR conditions followed by 454 pyrosequencing. Error rates ranged from 0.04-0.66% and were similar in amplified and non-amplified samples. Discrepancies were observed between forward and reverse reads, indicating that most errors were introduced during the pyrosequencing step. Using the 5-virus-mix, non-optimized, standard RT-PCR conditions introduced artificial recombinants in a fraction of at least 30% of the reads that subsequently led to an underestimation of true haplotype frequencies. We minimized the fraction of recombinants down to 0.9-2.6% by optimized, artifact-reducing RT-PCR conditions. This approach enabled correct haplotype reconstruction and frequency estimations consistent with reference data obtained by single genome amplification. RT-PCR conditions are crucial for correct frequency estimation and analysis of haplotypes in heterogeneous virus populations. We developed an RT-PCR procedure to generate NGS data useful for reliable haplotype reconstruction and quantification.

The implementation of genomic-based medicine is hindered by unresolved questions regarding data privacy and delivery of interpreted results to health-care practitioners. We used DNA-based prediction of HIV-related outcomes as a model to explore critical issues in clinical genomics.

MiRNAs and other small noncoding RNAs (sncRNAs) are key players in post-transcriptional gene regulation. HIV-1 derived small noncoding RNAs (sncRNAs) have been described in HIV-1 infected cells, but their biological functions still remain to be elucidated. Here, we approached the question whether viral sncRNAs may play a role in the RNA interference (RNAi) pathway or whether viral mRNAs are targeted by cellular miRNAs in human monocyte derived macrophages (MDM).

Candida esophagitis belongs to the most common AIDS-defining diseases; however, a comprehensive immune pathogenic concept is lacking.

HIV-1 infects CD4+ T cells and completes its replication cycle in approximately 24 hours. We employed repeated measurements in a standardized cell system and rigorous mathematical modeling to characterize the emergence of the viral replication intermediates and their impact on the cellular transcriptional response with high temporal resolution. We observed 7,991 (73%) of the 10,958 expressed genes to be modulated in concordance with key steps of viral replication. Fifty-two percent of the overall variability in the host transcriptome was explained by linear regression on the viral life cycle. This profound perturbation of cellular physiology was investigated in the light of several regulatory mechanisms, including transcription factors, miRNAs, host-pathogen interaction, and proviral integration. Key features were validated in primary CD4+ T cells, and with viral constructs using alternative entry strategies. We propose a model of early massive cellular shutdown and progressive upregulation of the cellular machinery to complete the viral life cycle.

In order to facilitate and improve the use of antiretroviral therapy (ART), international recommendations are released and updated regularly. We aimed to study if adherence to the recommendations is associated with better treatment outcomes in the Swiss HIV Cohort Study (SHCS).

Many viruses, including the clinically relevant RNA viruses HIV (human immunodeficiency virus) and HCV (hepatitis C virus), exist in large populations and display high genetic heterogeneity within and between infected hosts. Assessing intra-patient viral genetic diversity is essential for understanding the evolutionary dynamics of viruses, for designing effective vaccines, and for the success of antiviral therapy. Next-generation sequencing (NGS) technologies allow the rapid and cost-effective acquisition of thousands to millions of short DNA sequences from a single sample. However, this approach entails several challenges in experimental design and computational data analysis. Here, we review the entire process of inferring viral diversity from sample collection to computing measures of genetic diversity. We discuss sample preparation, including reverse transcription and amplification, and the effect of experimental conditions on diversity estimates due to in vitro base substitutions, insertions, deletions, and recombination. The use of different NGS platforms and their sequencing error profiles are compared in the context of various applications of diversity estimation, ranging from the detection of single nucleotide variants (SNVs) to the reconstruction of whole-genome haplotypes. We describe the statistical and computational challenges arising from these technical artifacts, and we review existing approaches, including available software, for their solution. Finally, we discuss open problems, and highlight successful biomedical applications and potential future clinical use of NGS to estimate viral diversity.

Understanding pathways that promote HIV-1 broadly neutralizing antibody (bnAb) induction is crucial to advance bnAb-based vaccines. We recently demarcated host, viral, and disease parameters associated with bnAb development in a large HIV-1 cohort screen. By establishing comprehensive antibody signatures based on IgG1, IgG2, and IgG3 activity to 13 HIV-1 antigens in 4,281 individuals in the same cohort, we now show that the same four parameters that are significantly linked with neutralization breadth, namely viral load, infection length, viral diversity, and ethnicity, also strongly influence HIV-1-binding antibody responses. However, the effects proved selective, shaping binding antibody responses in an antigen and IgG subclass-dependent manner. IgG response landscapes in bnAb inducers indicated a differentially regulated, IgG1-driven HIV-1 antigen response, and IgG1 binding of the BG505 SOSIP trimer proved the best predictor of HIV-1 neutralization breadth in plasma. Our findings emphasize the need to unravel immune modulators that underlie the differentially regulated IgG response in bnAb inducers to guide vaccine development.

The HIV-1 reservoir is the major hurdle to a cure. We here evaluate viral and host characteristics associated with reservoir size and long-term dynamics in 1,057 individuals on suppressive antiretroviral therapy for a median of 5.4 years. At the population level, the reservoir decreases with diminishing differences over time, but increases in 26.6% of individuals. Viral blips and low-level viremia are significantly associated with slower reservoir decay. Initiation of ART within the first year of infection, pretreatment viral load, and ethnicity affect reservoir size, but less so long-term dynamics. Viral blips and low-level viremia are thus relevant for reservoir and cure studies.

Human immunodeficiency virus (HIV) infection is associated with an increased risk of non-Hodgkin lymphoma (NHL). Even in the era of suppressive antiretroviral treatment, HIV-infected individuals remain at higher risk of developing NHL compared to the general population. To identify potential genetic risk loci, we performed case-control genome-wide association studies and a meta-analysis across three cohorts of HIV+ patients of European ancestry, including a total of 278 cases and 1924 matched controls. We observed a significant association with NHL susceptibility in the C-X-C motif chemokine ligand 12 (CXCL12) region on chromosome 10. A fine mapping analysis identified rs7919208 as the most likely causal variant (P = 4.77e-11), with the G>A polymorphism creating a new transcription factor binding site for BATF and JUND. These results suggest a modulatory role of CXCL12 regulation in the increased susceptibility to NHL observed in the HIV-infected population.

The potential of broadly neutralizing antibodies targeting the HIV-1 envelope trimer to prevent HIV-1 transmission has opened new avenues for therapies and vaccines. However, their implementation remains challenging and would profit from a deepened mechanistic understanding of HIV-antibody interactions and the mucosal transmission process. In this study we experimentally determined stoichiometric parameters of the HIV-1 trimer-antibody interaction, confirming that binding of one antibody is sufficient for trimer neutralization. This defines numerical requirements for HIV-1 virion neutralization and thereby enables mathematical modelling of in vitro and in vivo antibody neutralization efficacy. The model we developed accurately predicts antibody efficacy in animal passive immunization studies and provides estimates for protective mucosal antibody concentrations. Furthermore, we derive estimates of the probability for a single virion to start host infection and the risks of male-to-female HIV-1 transmission per sexual intercourse. Our work thereby delivers comprehensive quantitative insights into both the molecular principles governing HIV-antibody interactions and the initial steps of mucosal HIV-1 transmission. These insights, alongside the underlying, adaptable modelling framework presented here, will be valuable for supporting in silico pre-trial planning and post-hoc evaluation of HIV-1 vaccination or antibody treatment trials.

Drug resistant HIV is a major threat to the long-term efficacy of antiretroviral treatment. Around 10% of ART-naïve patients in Europe are infected with drug-resistant HIV type 1. Hence it is important to understand the dynamics of transmitted drug resistance evolution. Thanks to routinely performed drug resistance tests, HIV sequence data is increasingly available and can be used to reconstruct the phylogenetic relationship among viral lineages. In this study we employ a phylodynamic approach to quantify the fitness costs of major resistance mutations in the Swiss HIV cohort. The viral phylogeny reflects the transmission tree, which we model using stochastic birth-death-sampling processes with two types: hosts infected by a sensitive or resistant strain. This allows quantification of fitness cost as the ratio between transmission rates of hosts infected by drug resistant strains and transmission rates of hosts infected by drug sensitive strains. The resistance mutations 41L, 67N, 70R, 184V, 210W, 215D, 215S and 219Q (nRTI-related) and 103N, 108I, 138A, 181C, 190A (NNRTI-related) in the reverse trancriptase and the 90M mutation in the protease gene are included in this study. Among the considered resistance mutations, only the 90M mutation in the protease gene was found to have significantly higher fitness than the drug sensitive strains. The following mutations associated with resistance to reverse transcriptase inhibitors were found to be less fit than the sensitive strains: 67N, 70R, 184V, 219Q. The highest posterior density intervals of the transmission ratios for the remaining resistance mutations included in this study all included 1, suggesting that these mutations do not have a significant effect on viral transmissibility within the Swiss HIV cohort. These patterns are consistent with alternative measures of the fitness cost of resistance mutations. Overall, we have developed and validated a novel phylodynamic approach to estimate the transmission fitness cost of drug resistance mutations.