Rare HIV-1 transmitted/founder lineages identified by deep viral sequencing contribute to rapid shifts in dominant quasispecies during acute and early infection

Reconciling founder variant multiplicity of HIV-1 infection with the rate of CD4+ decline

HIV-1 transmission precipitates a stringent genetic bottleneck, with 75% of new infections initiated by a single genetic variant. Where multiple variants initiate infection, recipient set point viral load (SpVL) and the rate of CD4+ T cell decline may be elevated, but these findings remain inconsistent. Here, we summarised the evidence for this phenomenon, then tested whether previous studies possessed sufficient statistical power to reliably identify a true effect of multiple variant infection leading to higher SpVL. Next, we combined models of HIV-1 transmission, heritability and disease progression to understand whether available data suggest a faster CD4+ T cell decline would be expected to associated with multiple variant infection, without an explicit dependency between the two. First, we found that most studies had insufficient power to identify a true significant difference, prompting an explanation for previous inconsistencies. Next, our model framework revealed we would not expect to observe a positive association between multiple variant infections and faster CD4+ T cell decline, in the absence of an explicit dependency. Consequently, while empirical evidence may be consistent with a positive association between multiple variant infection and faster CD4+ T cell decline, further investigation is required to establish a causal basis.

HIV-1 diversity in viral reservoirs obtained from circulating T-cell subsets during early ART and beyond

Even during extended periods of effective immunological control, a substantial dynamic of the viral genome can be observed in different cellular compartments in HIV-1 positive individuals, indicating the persistence of active viral reservoirs. To obtain further insights, we studied changes in the proviral as well as in the viral HIV-1 envelope (Env) sequence along with transcriptional, translational and viral outgrowth activity as indicators for viral dynamics and genomic intactness. Our study identified distinct reservoir patterns that either represented highly sequence-diverse HIV-1 populations or only a single / few persisting virus variants. The single dominating variants were more often found in individuals starting ART during early infection phases, indicating that early treatment might limit reservoir diversification. At the same time, more sequence-diverse HIV reservoirs correlated with a poorer immune status, indicated by lower CD4 count, a higher number of regimen changes and more co-morbidities. Furthermore, we noted that in T-cell populations in the peripheral blood, replication-competent HIV-1 is predominantly present in Lymph node homing TN (naïve) and TCM (central memory) T cells. Provirus genomes archived in TTM (transitional memory) and TEM (effector memory) T cells more frequently tended to carry inactivating mutations and, population-wise, possess changes in the genetic diversity. These discriminating properties of the viral reservoir in T-cell subsets may have important implications for new early therapy strategies, underscoring the critical role of early therapy in preserving robust immune surveillance and constraining the viral reservoir.

Evolution of Antibody Responses in HIV-1 CRF01_AE Acute Infection: Founder Envelope V1V2 Impacts the Timing and Magnitude of Autologous Neutralizing Antibodies

Understanding the dynamics of early immune responses to HIV-1 infection, including the evolution of initial neutralizing and antibody-dependent cellular cytotoxicity (ADCC)-mediating antibodies, will inform HIV vaccine design. In this study, we assess the development of autologous neutralizing antibodies (ANAbs) against founder envelopes (Envs) from 18 participants with HIV-1 CRF01_AE acute infection. The timing of ANAb development directly associated with the magnitude of the longitudinal ANAb response. Participants that developed ANAbs within 6 months of infection had significantly higher ANAb responses at 1 year (50% inhibitory concentration [IC50] geometric mean titer [GMT] = 2,010 versus 184; P = 0.001) and 2 years (GMT = 3,479 versus 340; P = 0.015), compared to participants that developed ANAb responses after 6 months. Participants with later development of ANAb tended to develop an earlier, potent heterologous tier 1 (92TH023) neutralizing antibody (NAb) response (P = 0.049). CRF01_AE founder Env V1V2 loop lengths correlated indirectly with the timing (P = 0.002, r = -0.675) and directly with magnitude (P = 0.005, r = 0.635) of ANAb responses; Envs with longer V1V2 loop lengths elicited earlier and more potent ANAb responses. While ANAb responses did not associate with viral load, the viral load set point correlated directly with neutralization of the heterologous 92TH023 strain (P = 0.007, r = 0.638). In contrast, a striking inverse correlation was observed between viral load set point and peak ADCC against heterologous 92TH023 Env strain (P = 0.0005, r = -0.738). These data indicate that specific antibody functions can be differentially related to viral load set point and may affect HIV-1 pathogenesis. Exploiting Env properties, such as V1V2 length, could facilitate development of subtype-specific vaccines that elicit more effective immune responses and improved protection. IMPORTANCE Development of an effective HIV-1 vaccine will be facilitated by better understanding the dynamics between the founder virus and the early humoral responses. Variations between subtypes may influence the evolution of immune responses and should be considered as we strive to understand these dynamics. In this study, autologous founder envelope neutralization and heterologous functional humoral responses were evaluated after acute infection by HIV-1 CRF01_AE, a subtype that has not been thoroughly characterized. The evolution of these humoral responses was assessed in relation to envelope characteristics, magnitude of elicited immune responses, and viral load. Understanding immune parameters in natural infection will improve our understanding of protective responses and aid in the development of immunogens that elicit protective functional antibodies. Advancing our knowledge of correlates of positive clinical outcomes should lead to the design of more efficacious vaccines.

A longitudinal analysis of immune escapes from HLA-B*13-restricted T-cell responses at early stage of CRF01_AE subtype HIV-1 infection and implications for vaccine design

BACKGROUND

Identifying immunogens which can elicit effective T cell responses against human immunodeficiency virus type 1 (HIV-1) is important for developing a T-cell based vaccine. It has been reported that human leukocyte antigen (HLA)-B*13-restricted T-cell responses contributed to HIV control in subtype B' and C infected individuals. However, the kinetics of B*13-restricted T-cell responses, viral evolution within epitopes, and the impact on disease progression in CRF01_AE subtype HIV-1-infected men who have sex with men (MSM) are not known.

RESULTS

Interferon-γ ELISPOT assays and deep sequencing of viral RNAs were done in 14 early HLA-B*13-positive CRF01_AE subtype HIV-1-infected MSM. We found that responses to RQEILDLWV (Nef_106-114, RV9), GQMREPRGSDI (Gag_226-236, GI11), GQDQWTYQI (Pol_487-498, GI9), and VQNAQGQMV (Gag_135-143, VV9) were dominant. A higher relative magnitude of Gag-specific T-cell responses, contributed to viral control, whereas Nef-specific T-cell responses were associated with rapid disease progression. GI11 (Gag) was conserved and strong GI11 (Gag)-specific T-cell responses showed cross-reactivity with a dominant variant, M228I, found in 3/12 patients; GI11 (Gag)-specific T-cell responses were positively associated with CD4 T-cell counts (R = 0.716, P = 0.046). Interestingly, the GI9 (Pol) epitope was also conserved, but GI9 (Pol)-specific T-cell responses did not influence disease progression (P > 0.05), while a D490G variant identified in one patient did not affect CD4 T-cell counts. All the other epitopes studied [VV9 (Gag), RQYDQILIEI (Pol_113-122, RI10), HQSLSPRTL (Gag_144-152, HL9), and RQANFLGRL (Gag_429-437, RL9)] developed escape mutations within 1 year of infection, which may have contributed to overall disease progression. Intriguingly, we found early RV9 (Nef)-specific T-cell responses were associated with rapid disease progression, likely due to escape mutations.

CONCLUSIONS

Our study strongly suggested the inclusion of GI11 (Gag) and exclusion of RV9 (Nef) for T-cell-based vaccine design for B*13-positive CRF01_AE subtype HIV-1-infected MSM and high-risk individuals.

The mucosal barrier and anti-viral immune responses can eliminate portions of the viral population during transmission and early viral growth

Little is known about how specific individual viral lineages replicating systemically during acute Human Immunodeficiency Virus or Simian Immunodeficiency Virus (HIV/SIV) infection persist into chronic infection. In this study, we use molecularly barcoded SIV (SIVmac239M) to track distinct viral lineages for 12 weeks after intravenous (IV) or intrarectal (IR) challenge in macaques. Two Mafa-A1*063+ cynomolgus macaques (Macaca fascicularis, CM) were challenged IV, and two Mamu-A1*001+ rhesus macaques (Macaca mulatta, RM) were challenged IR with 200,000 Infectious Units (IU) of SIVmac239M. We sequenced the molecular barcode of SIVmac239M from all animals over the 12 weeks of the study to characterize the diversity and persistence of virus lineages. During the first three weeks post-infection, we found ~70-560 times more unique viral lineages circulating in the animals challenged IV compared to those challenged IR, which is consistent with the hypothesis that the challenge route is the primary driver restricting the transmission of individual viral lineages. We also characterized the sequences of T cell epitopes targeted during acute SIV infection, and found that the emergence of escape variants in acutely targeted epitopes can occur on multiple virus templates simultaneously, but that elimination of some of these templates is likely a consequence of additional host factors. These data imply that virus lineages present during acute infection can still be eliminated from the systemic virus population even after initial selection.

Quasispecies of SARS-CoV-2 revealed by single nucleotide polymorphisms (SNPs) analysis

New SARS-CoV-2 mutants have been continuously indentified with enhanced transmission ever since its outbreak in early 2020. As an RNA virus, SARS-CoV-2 has a high mutation rate due to the low fidelity of RNA polymerase. To study the single nucleotide polymorphisms (SNPs) dynamics of SARS-CoV-2, 158 SNPs with high confidence were identified by deep meta-transcriptomic sequencing, and the most common SNP type was C > T. Analyses of intra-host population diversity revealed that intra-host quasispecies' composition varies with time during the early onset of symptoms, which implicates viral evolution during infection. Network analysis of co-occurring SNPs revealed the most abundant non-synonymous SNP 22,638 in the S glycoprotein RBD region and 28,144 in the ORF8 region. Furthermore, SARS-CoV-2 variations differ in an individual's respiratory tissue (nose, throat, BALF, or sputum), suggesting independent compartmentalization of SARS-CoV-2 populations in patients. The positive selection analysis of the SARS-CoV-2 genome uncovered the positive selected amino acid G251V on ORF3a. Alternative allele frequency spectrum (AAFS) of all variants revealed that ORF8 could bear alternate alleles with high frequency. Overall, the results show the quasispecies' profile of SARS-CoV-2 in the respiratory tract in the first two months after the outbreak.

Predator-Prey Dynamics of Intra-Host Simian Immunodeficiency Virus Evolution Within the Untreated Host

The dynamic nature of the SIV population during disease progression in the SIV/macaque model of AIDS and the factors responsible for its behavior have not been documented, largely owing to the lack of sufficient spatial and temporal sampling of both viral and host data from SIV-infected animals. In this study, we detail Bayesian coalescent inference of the changing collective intra-host viral effective population size (Ne ) from various tissues over the course of infection and its relationship with what we demonstrate is a continuously changing immune cell repertoire within the blood. Although the relative contribution of these factors varied among hosts and time points, the adaptive immune response best explained the overall periodic dynamic behavior of the effective virus population. Data exposing the nature of the relationship between the virus and immune cell populations revealed the plausibility of an eco-evolutionary mathematical model, which was able to mimic the large-scale oscillations in Ne through virus escape from relatively few, early immunodominant responses, followed by slower escape from several subdominant and weakened immune populations. The results of this study suggest that SIV diversity within the untreated host is governed by a predator-prey relationship, wherein differing phases of infection are the result of adaptation in response to varying immune responses. Previous investigations into viral population dynamics using sequence data have focused on single estimates of the effective viral population size (Ne ) or point estimates over sparse sampling data to provide insight into the precise impact of immune selection on virus adaptive behavior. Herein, we describe the use of the coalescent phylogenetic frame- work to estimate the relative changes in Ne over time in order to quantify the relationship with empirical data on the dynamic immune composition of the host. This relationship has allowed us to expand on earlier simulations to build a predator-prey model that explains the deterministic behavior of the virus over the course of disease progression. We show that sequential viral adaptation can occur in response to phases of varying immune pressure, providing a broader picture of the viral response throughout the entire course of progression to AIDS.

New SHIVs and Improved Design Strategy for Modeling HIV-1 Transmission, Immunopathogenesis, Prevention and Cure

Previously, we showed that substitution of HIV-1 Env residue 375-Ser by bulky aromatic residues enhances binding to rhesus CD4 and enables primary HIV-1 Envs to support efficient replication as simian-human immunodeficiency virus (SHIV) chimeras in rhesus macaques (RMs). Here, we test this design strategy more broadly by constructing SHIVs containing ten primary Envs corresponding to HIV-1 subtypes A, B, C, AE and AG. All ten SHIVs bearing wildtype Env375 residues replicated efficiently in human CD4+ T cells, but only one replicated efficiently in primary rhesus cells. This was a subtype AE SHIV that naturally contained His at Env375. Replacement of wildtype Env375 residues by Trp, Tyr, Phe or His in the other nine SHIVs led to efficient replication in rhesus CD4+ T cells in vitro and in vivo Nine SHIVs containing optimized Env375 alleles were grown large-scale in primary rhesus CD4+ T cells to serve as challenge stocks in preclinical prevention trials. These virus stocks were genetically homogeneous, native-like in Env antigenicity and tier-2 neutralization sensitivity, and transmissible by rectal, vaginal, penile, oral or intravenous routes. To facilitate future SHIV constructions, we engineered a simplified second-generation design scheme and validated it in RMs. Overall, our findings demonstrate that SHIVs bearing primary Envs with bulky aromatic substitutions at Env375 consistently replicate in RMs, recapitulating many features of HIV-1 infection in humans. Such SHIVs are efficiently transmitted by mucosal routes common to HIV-1 infection and can be used to test vaccine efficacy in preclinical monkey trials.ImportanceSHIV infection of Indian rhesus macaques is an important animal model for studying HIV-1 transmission, prevention, immunopathogenesis and cure. Such research is timely, given recent progress with active and passive immunization and novel approaches to HIV-1 cure. Given the multifaceted roles of HIV-1 Env in cell tropism and virus entry, and as a target for neutralizing and non-neutralizing antibodies, Envs selected for SHIV construction are of paramount importance. Until recently, it has been impossible to strategically design SHIVs bearing clinically relevant Envs that replicate consistently in monkeys. This changed with the discovery that bulky aromatic substitutions at residue Env375 confer enhanced affinity to rhesus CD4. Here, we show that 10 new SHIVs bearing primary HIV-1 Envs with residue 375 substitutions replicated efficiently in RMs and could be transmitted efficiently across rectal, vaginal, penile and oral mucosa. These findings suggest an expanded role for SHIVs as a model of HIV-1 infection.

B cell engagement with HIV-1 founder virus envelope predicts development of broadly neutralizing antibodies

Determining which immunological mechanisms contribute to the development of broad neutralizing antibodies (bNAbs) during HIV-1 infection is a major goal to inform vaccine design. Using samples from a longitudinal HIV-1 acute infection cohort, we found key B cell determinants within the first 14-43 days of viremia that predict the development of bNAbs years later. Individuals who develop neutralization breadth had significantly higher B cell engagement with the autologous founder HIV envelope (Env) within 1 month of initial viremia. A higher frequency of founder-Env-specific naive B cells was associated with increased B cell activation and differentiation and predictive of bNAb development. These data demonstrate that the initial B cell interaction with the founder HIV Env is important for the development of broadly neutralizing antibodies and provide evidence that events within HIV acute infection lead to downstream functional outcomes.

Factors influencing estimates of HIV-1 infection timing using BEAST

While large datasets of HIV-1 sequences are increasingly being generated, many studies rely on a single gene or fragment of the genome and few comparative studies across genes have been done. We performed genome-based and gene-specific Bayesian phylogenetic analyses to investigate how certain factors impact estimates of the infection dates in an acute HIV-1 infection cohort, RV217. In this cohort, HIV-1 diagnosis corresponded to the first RNA positive test and occurred a median of four days after the last negative test, allowing us to compare timing estimates using BEAST to a narrow window of infection. We analyzed HIV-1 sequences sampled one week, one month and six months after HIV-1 diagnosis in 39 individuals. We found that shared diversity and temporal signal was limited in acute infection, and insufficient to allow timing inferences in the shortest HIV-1 genes, thus dated phylogenies were primarily analyzed for env, gag, pol and near full-length genomes. There was no one best-fitting model across participants and genes, though relaxed molecular clocks (73% of best-fitting models) and the Bayesian skyline (49%) tended to be favored. For infections with single founders, the infection date was estimated to be around one week pre-diagnosis for env (IQR: 3–9 days) and gag (IQR: 5–9 days), whilst the genome placed it at a median of 10 days (IQR: 4–19). Multiply-founded infections proved problematic to date. Our ability to compare timing inferences to precise estimates of HIV-1 infection (within a week) highlights that molecular dating methods can be applied to within-host datasets from early infection. Nonetheless, our results also suggest caution when using uniform clock and population models or short genes with limited information content.

Molecular dating using phylogenetics allows us to estimate the date of an infection from time-stamped within-host sequences alone. There are large datasets of HIV-1 sequences, but genome and gene analyses are not often performed in parallel and rarely with the possibility to compare results against a known narrow window of infection. We showed that all but the longest genes are near-clonal in acute infection, with little information for dating purposes. For infections with single founders, we estimated the eclipse phase—the time between HIV-1 exposure and the first positive diagnostic test—to last between one and two weeks using env, gag, pol and near full-length genomes. This approach could be used to narrow the date of suspected infection in ongoing clinical trials for the prevention of HIV-1 infection.

Cross-scale dynamics and the evolutionary emergence of infectious diseases

When emerging pathogens encounter new host species for which they are poorly adapted, they must evolve to escape extinction. Pathogens experience selection on traits at multiple scales, including replication rates within host individuals and transmissibility between hosts. We analyze a stochastic model linking pathogen growth and competition within individuals to transmission between individuals. Our analysis reveals a new factor, the cross-scale reproductive number of a mutant virion, that quantifies how quickly mutant strains increase in frequency when they initially appear in the infected host population. This cross-scale reproductive number combines with viral mutation rates, single-strain reproductive numbers, and transmission bottleneck width to determine the likelihood of evolutionary emergence, and whether evolution occurs swiftly or gradually within chains of transmission. We find that wider transmission bottlenecks facilitate emergence of pathogens with short-term infections, but hinder emergence of pathogens exhibiting cross-scale selective conflict and long-term infections. Our results provide a framework to advance the integration of laboratory, clinical, and field data in the context of evolutionary theory, laying the foundation for a new generation of evidence-based risk assessment of emergence threats.

Streamlined Subpopulation, Subtype, and Recombination Analysis of HIV-1 Half-Genome Sequences Generated by High-Throughput Sequencing

The highly recombinogenic nature of human immunodeficiency virus type 1 (HIV-1) leads to recombination and emergence of quasispecies. It is important to reliably identify subpopulations to understand the complexity of a viral population for drug resistance surveillance and vaccine development. High-throughput sequencing (HTS) provides improved resolution over Sanger sequencing for the analysis of heterogeneous viral subpopulations. However, current methods of analysis of HTS reads are unable to fully address accurate population reconstruction. Hence, there is a dire need for a more sensitive, accurate, user-friendly, and cost-effective method to analyze viral quasispecies. For this purpose, we have improved the HIVE-hexahedron algorithm that we previously developed with in silico short sequences to analyze raw HTS short reads. The significance of this study is that our standalone algorithm enables a streamlined analysis of quasispecies, subtype, and recombination patterns from long HIV-1 genome regions without the need of additional sequence analysis tools. Distinct viral populations and recombination patterns identified by HIVE-hexahedron are further validated by comparison with sequences obtained by single genome sequencing (SGS).

High-throughput sequencing (HTS) has been widely used to characterize HIV-1 genome sequences. There are no algorithms currently that can directly determine genotype and quasispecies population using short HTS reads generated from long genome sequences without additional software. To establish a robust subpopulation, subtype, and recombination analysis workflow, we amplified the HIV-1 3′-half genome from plasma samples of 65 HIV-1-infected individuals and sequenced the entire amplicon (∼4,500 bp) by HTS. With direct analysis of raw reads using HIVE-hexahedron, we showed that 48% of samples harbored 2 to 13 subpopulations. We identified various subtypes (17 A1s, 4 Bs, 27 Cs, 6 CRF02_AGs, and 11 unique recombinant forms) and defined recombinant breakpoints of 10 recombinants. These results were validated with viral genome sequences generated by single genome sequencing (SGS) or the analysis of consensus sequence of the HTS reads. The HIVE-hexahedron workflow is more sensitive and accurate than just evaluating the consensus sequence and also more cost-effective than SGS.

IMPORTANCE The highly recombinogenic nature of human immunodeficiency virus type 1 (HIV-1) leads to recombination and emergence of quasispecies. It is important to reliably identify subpopulations to understand the complexity of a viral population for drug resistance surveillance and vaccine development. High-throughput sequencing (HTS) provides improved resolution over Sanger sequencing for the analysis of heterogeneous viral subpopulations. However, current methods of analysis of HTS reads are unable to fully address accurate population reconstruction. Hence, there is a dire need for a more sensitive, accurate, user-friendly, and cost-effective method to analyze viral quasispecies. For this purpose, we have improved the HIVE-hexahedron algorithm that we previously developed with in silico short sequences to analyze raw HTS short reads. The significance of this study is that our standalone algorithm enables a streamlined analysis of quasispecies, subtype, and recombination patterns from long HIV-1 genome regions without the need of additional sequence analysis tools. Distinct viral populations and recombination patterns identified by HIVE-hexahedron are further validated by comparison with sequences obtained by single genome sequencing (SGS).

A Department of Defense Laboratory Consortium Approach to Next Generation Sequencing and Bioinformatics Training for Infectious Disease Surveillance in Kenya

Epidemics of emerging and re-emerging infectious diseases are a danger to civilian and military populations worldwide. Health security and mitigation of infectious disease threats is a priority of the United States Government and the Department of Defense (DoD). Next generation sequencing (NGS) and Bioinformatics (BI) enhances traditional biosurveillance by providing additional data to understand transmission, identify resistance and virulence factors, make predictions, and update risk assessments. As more and more laboratories adopt NGS and BI technologies they encounter challenges in building local capacity. In addition to choosing the right sequencing platform and approach, considerations must also be made for the complexity of bioinformatics analyses, data storage, as well as personnel and computational requirements. To address these needs, a comprehensive training program was developed covering wet lab and bioinformatics approaches to NGS. The program is meant to be modular and adaptive to meet both common and individualized needs of medical research and public health laboratories across the DoD. The training program was first deployed internationally to the Basic Science Laboratory of the US Army Medical Research Directorate-Africa in Kisumu, Kenya, which is an overseas Lab of the Walter Reed Army Institute of Research (WRAIR). A week-long workshop with intensive focus on targeted sequencing and the bioinformatics of genome assembly (n = 24 participants) was held. Post-workshop self-assessment (completed by 21 participants) noted significant median gains in knowledge domains related to NGS targeted sequencing, bioinformatics for genome assembly, and sequence quality assessment. The participants also reported that the information on study design, sample preparation, sequencing quality control, data quality assessment, reporting, and basic and advanced bioinformatics analysis were the most useful information presented in the training. While longer-term evaluations are planned, the training resulted in significant short-term improvement of a laboratory’s self-reported wet lab and bioinformatics capabilities. This framework can be used for future DoD laboratory development in the area of NGS and BI for infectious disease surveillance, ultimately enhancing this global DoD capability.

Dynamics of Human Immunodeficiency Virus-1 Genetic Diversification During Acute Infection

We analyzed human immunodeficiency virus envelope diversity in 98 acute infections. The within-host genetic diversity, divergence from transmitted/founder (T/F) strain, and the observed frequency of multiple T/F infections increased with Fiebig stage. These data identify rapid viral dynamics during acute infection with implications for clinical trials conducted in this setting.

Infection with multiple HIV-1 founder variants is associated with lower viral replicative capacity, faster CD4+ T cell decline and increased immune activation during acute infection

HIV-1 transmission is associated with a severe bottleneck in which a limited number of variants from a pool of genetically diverse quasispecies establishes infection. The IAVI protocol C cohort of discordant couples, female sex workers, other heterosexuals and men who have sex with men (MSM) present varying risks of HIV infection, diverse HIV-1 subtypes and represent a unique opportunity to characterize transmitted/founder viruses (TF) where disease outcome is known. To identify the TF, the HIV-1 repertoire of 38 MSM participants' samples was sequenced close to transmission (median 21 days post infection, IQR 18-41) and assessment of multivariant infection done. Patient derived gag genes were cloned into an NL4.3 provirus to generate chimeric viruses which were characterized for replicative capacity (RC). Finally, an evaluation of how the TF virus predicted disease progression and modified the immune response at both acute and chronic HIV-1 infection was done. There was higher prevalence of multivariant infection compared with previously described heterosexual cohorts. A link was identified between multivariant infection and replicative capacity conferred by gag, whereby TF gag tended to be of lower replicative capacity in multivariant infection (p = 0.02) suggesting an overall lowering of fitness requirements during infection with multiple variants. Notwithstanding, multivariant infection was associated with rapid CD4+ T cell decline and perturbances in the CD4+ T cell and B cell compartments compared to single variant infection, which were reversible upon control of viremia. Strategies aimed at identifying and mitigating multivariant infection could contribute toward improving HIV-1 prognosis and this may involve strategies that tighten the stringency of the transmission bottleneck such as treatment of STI. Furthermore, the sequences and chimeric viruses help with TF based experimental vaccine immunogen design and can be used in functional assays to probe effective immune responses against TF.

The HIV-1 Env gp120 Inner Domain Shapes the Phe43 Cavity and the CD4 Binding Site

The HIV-1 envelope glycoproteins (Env) undergo conformational changes upon interaction of the gp120 exterior glycoprotein with the CD4 receptor. The gp120 inner domain topological layers facilitate the transition of Env to the CD4-bound conformation. CD4 engages gp120 by introducing its phenylalanine 43 (Phe43) in a cavity ("the Phe43 cavity") located at the interface between the inner and outer gp120 domains. Small CD4-mimetic compounds (CD4mc) can bind within the Phe43 cavity and trigger conformational changes similar to those induced by CD4. Crystal structures of CD4mc in complex with a modified CRF01_AE gp120 core revealed the importance of these gp120 inner domain layers in stabilizing the Phe43 cavity and shaping the CD4 binding site. Our studies reveal a complex interplay between the gp120 inner domain and the Phe43 cavity and generate useful information for the development of more-potent CD4mc.IMPORTANCE The Phe43 cavity of HIV-1 envelope glycoproteins (Env) is an attractive druggable target. New promising compounds, including small CD4 mimetics (CD4mc), were shown to insert deeply into this cavity. Here, we identify a new network of residues that helps to shape this highly conserved CD4 binding pocket and characterize the structural determinants responsible for Env sensitivity to small CD4 mimetics.

Into the Deep (Sequence) of the Foot-and-Mouth Disease Virus Gene Pool: Bottlenecks and Adaptation during Infection in Naïve and Vaccinated Cattle

Foot-and-mouth disease virus (FMDV) infects hosts as a population of closely related viruses referred to as a quasispecies. The behavior of this quasispecies has not been described in detail in natural host species. In this study, virus samples collected from vaccinated and non-vaccinated cattle up to 35 days post-experimental infection with FMDV A24-Cruzeiro were analyzed by deep-sequencing. Vaccination induced significant differences compared to viruses from non-vaccinated cattle in substitution rates, entropy, and evidence for adaptation. Genomic variation detected during early infection reflected the diversity inherited from the source virus (inoculum), whereas by 12 days post infection, dominant viruses were defined by newly acquired mutations. Mutations conferring recognized fitness gain occurred and were associated with selective sweeps. Persistent infections always included multiple FMDV subpopulations, suggesting distinct foci of infection within the nasopharyngeal mucosa. Subclinical infection in vaccinated cattle included very early bottlenecks associated with reduced diversity within virus populations. Viruses from both animal cohorts contained putative antigenic escape mutations. However, these mutations occurred during later stages of infection, at which time transmission is less likely to occur. This study improves upon previously published work by analyzing deep sequences of samples, allowing for detailed characterization of FMDV populations over time within multiple hosts.

In vitro replicative fitness of early Transmitted founder HIV-1 variants and sensitivity to Interferon alpha

Type I interferons, particularly interferon-alpha (IFN-α), play a vital role in the host's anti-viral defenses by interfering with viral replication. However, the virus rapidly evolves to exploit the IFN-α response for its replication, spread, and pathogenic function. In this study, we attempted to determine IFN-α susceptibility and productivity of infectious transmitted/founder (TF) (n = 8) and non-transmitted (NT) viruses (n = 8) derived from HIV-1 infected infants. Independent experiments were carried out to determine IFN-α resistance, replication fitness, and viral productivity in CD4⁺ T cells over a short period. In vitro studies showed that TF viruses were resistant to IFN-α during the very near moment of transmission, but in the subsequent time points, they became susceptible to IFN-α. We did not observe much difference in replicative fitness of the TF viruses in cultures treated with and without IFN-α, but the difference was significant in the case of NT viruses obtained from the same individual. Despite increased susceptibility to IFN-α, NT viruses produced more viral particles than TF viruses. Similar results were also obtained in cultures treated with maraviroc (MVC). The study identified unique characteristics of TF viruses thus prompting further investigation into virus-host interaction occurring during the early stages of HIV infection.

Genetically barcoded SIV reveals the emergence of escape mutations in multiple viral lineages during immune escape

The rapidity of replication coupled with a high mutation rate enables HIV to evade selective pressures imposed by host immune responses. Investigating the ability of HIV to escape different selection forces has generally relied on population-level measures, such as the time to detectable escape mutations in plasma and the rate these mutations subsequently take over the virus population. Here we employed a barcoded synthetic swarm of simian immunodeficiency virus (SIV) in rhesus macaques to investigate the generation and selection of escape mutations within individual viral lineages at the Mamu-A*01-restricted Tat-SL8 epitope. We observed the persistence of more than 1,000 different barcode lineages following selection after acquiring escape mutations. Furthermore, the increased resolution into the virus population afforded by barcode analysis revealed changes in the population structure of the viral quasispecies as it adapted to immune pressure. The high frequency of emergence of escape mutations in parallel viral lineages at the Tat-SL8 epitope highlights the challenge posed by viral escape for the development of T cell-based vaccines. Importantly, the level of viral replication required for generating escape mutations in individual lineages can be directly estimated using the barcoded virus, thereby identifying the level of efficacy required for a successful vaccine to limit escape. Overall, assessing the survival of barcoded viral lineages during selection provides a direct and quantitative measure of the stringency of the underlying genetic bottleneck, making it possible to predict the ability of the virus to escape selective forces induced by host immune responses as well as during therapeutic interventions.

Terminal Effector CD8 T Cells Defined by an IKZF2+IL-7R- Transcriptional Signature Express FcγRIIIA, Expand in HIV Infection, and Mediate Potent HIV-Specific Antibody-Dependent Cellular Cytotoxicity

Chronic HIV-1 is associated with increased levels of FcγRIIIA⁺ CD8 T cells.

FcγRIIIA⁺ CD8 T cells display an innate transcriptomic profile akin to NK cells.

ADCC is mediated by FcγRIIIA⁺ CD8 T cells at levels comparable with NK cells.

HIV-1 infection expands large populations of late-stage differentiated CD8 T cells that may persist long after viral escape from TCR recognition. In this study, we investigated whether such CD8 T cell populations can perform unconventional innate-like antiviral effector functions. Chronic untreated HIV-1 infection was associated with elevated numbers of CD45RA⁺CD57⁺ terminal effector CD8 T cells expressing FcγRIIIA (CD16). The FcγRIIIA⁺ CD8 T cells displayed a distinctive transcriptional profile between conventional CD8 T cells and NK cells, characterized by high levels of IKZF2 and low expression of IL7R. This transcriptional profile translated into a distinct NKp80⁺ IL-7Rα⁻ surface phenotype with high expression of the Helios transcription factor. Interestingly, the FcγRIIIA⁺ CD8 T cells mediated HIV-specific Ab-dependent cellular cytotoxicity (ADCC) activity at levels comparable with NK cells on a per cell basis. The FcγRIIIA⁺ CD8 T cells were highly activated in a manner that correlated positively with expansion of the CD8 T cell compartment and with plasma levels of soluble mediators of antiviral immunity and inflammation such as IP-10, TNF, IL-6, and TNFRII. The frequency of FcγRIIIA⁺ CD8 T cells persisted as patients initiated suppressive antiretroviral therapy, although their activation levels declined. These data indicate that terminally differentiated effector CD8 T cells acquire enhanced innate cell-like characteristics during chronic viral infection and suggest that HIV-specific ADCC is a function CD8 T cells use to target HIV-infected cells. Furthermore, as the FcγRIIIA⁺ CD8 T cells persist in treatment, they contribute significantly to the ADCC-capable effector cell pool in patients on antiretroviral therapy.

Deep Sequencing Reveals Central Nervous System Compartmentalization in Multiple Transmitted/Founder Virus Acute HIV-1 Infection

HIV-1 disseminates to a broad range of tissue compartments during acute HIV-1 infection (AHI). The central nervous system (CNS) can serve as an early and persistent site of viral replication, which poses a potential challenge for HIV-1 remission strategies that target the HIV reservoir. CNS compartmentalization is a key feature of HIV-1 neuropathogenesis. Thus far, the timing of how early CNS compartmentalization develops after infection is unknown. We examined whether HIV-1 transmitted/founder (T/F) viruses differ between CNS and blood during AHI using single-genome sequencing of envelope gene and further examined subregions in pol and env using next-generation sequencing in paired plasma and cerebrospinal fluid (CSF) from 18 individuals. Different proportions of mostly minor variants were found in six of the eight multiple T/F-infected individuals, indicating enrichment of some variants in CSF that may lead to significant compartmentalization in the later stages of infection. This study provides evidence for the first time that HIV-1 compartmentalization in the CNS can occur within days of HIV-1 exposure in multiple T/F infections. Further understanding of factors that determine enrichment of T/F variants in the CNS, as well as potential long-term implications of these findings for persistence of HIV-1 reservoirs and neurological impairment in HIV, is needed.

Dolutegravir reshapes the genetic diversity of HIV-1 reservoirs

Objectives

Better understanding of the dynamics of HIV reservoirs under ART is a critical step to achieve a functional HIV cure. Our objective was to assess the genetic diversity of archived HIV-1 DNA over 48 weeks in blood cells of individuals starting treatment with a dolutegravir-based regimen.

Methods

Eighty blood samples were prospectively and longitudinally collected from 20 individuals (NCT02557997) including: acutely (n = 5) and chronically (n = 5) infected treatment-naive individuals, as well as treatment-experienced individuals who switched to a dolutegravir-based regimen and were either virologically suppressed (n = 5) or had experienced treatment failure (n = 5). The integrase and V3 loop regions of HIV-1 DNA isolated from PBMCs were analysed by pyrosequencing at baseline and weeks 4, 24 and 48. HIV-1 genetic diversity was calculated using Shannon entropy.

Results

All individuals achieved or maintained viral suppression throughout the study. A low and stable genetic diversity of archived HIV quasispecies was observed in individuals starting treatment during acute infection. A dramatic reduction of the genetic diversity was observed at week 4 of treatment in the other individuals. In these patients and despite virological suppression, a recovery of the genetic diversity of the reservoirs was observed up to 48 weeks. Viral variants bearing dolutegravir resistance-associated substitutions at integrase position 50, 124, 230 or 263 were detected in five individuals (n = 5/20, 25%) from all groups except those who were ART-failing at baseline. None of these substitutions led to virological failure.

Conclusions

These data demonstrate that the genetic diversity of the HIV-1 reservoir is reshaped following the initiation of a dolutegravir-based regimen and strongly suggest that HIV-1 can continue to replicate despite successful treatment.

Humoral Response to the HIV-1 Envelope V2 Region in a Thai Early Acute Infection Cohort

Reduced risk of HIV-1 infection correlated with antibody responses to the envelope variable 1 and 2 regions in the RV144 vaccine trial. To understand the relationship between antibody responses, V2 sequence, and structure, plasma samples (n = 16) from an early acute HIV-1 infection cohort from Thailand infected with CRF01_AE strain were analyzed for binding to V2 peptides by surface plasmon resonance. Five participants with a range of V2 binding responses at week 24 post-infection were further analyzed against a set of four overlapping V2 peptides that were designed based on envelope single-genome amplification. Antibody responses that were relatively consistent over the four segments of the V2 region or a focused response to the C-strand (residues 165–186) of the V2 region were observed. Viral escape in the V2 region resulted in significantly reduced antibody binding. Structural modeling indicated that the C-strand and the sites of viral variation were highly accessible in the open conformation of the HIV-1 Env trimer. V2 residues, 165–186 are preferentially targeted during acute infection. Residues 169–184 were also preferentially targeted by the protective immune response in the RV144 trial, thus emphasizing the importance of these residues for vaccine design.

Caution is needed in interpreting HIV transmission chains by ultradeep sequencing

OBJECTIVES

Molecular epidemiology is applied to various aspects of HIV transmission analyses. With ultradeep sequencing (UDS), in-depth characterization of transmission episodes involving minority variants is permitted. We explored HIV-1 epidemiological linkage and evaluated characteristics of transmission dynamics and transmitted drug resistance (TDR) detection through the added value of UDS.

DESIGN

HIV pol gene fragments were sequenced by UDS and Sanger sequencing on samples of 70 HIV-1-infected, treatment-naive recently diagnosed MSM.

METHODS

Pairwise genetic distances and maximum likelihood phylogenies were computed. Transmission events were identified as clades with branch support at least 70% and intraclade genetic difference less than 4.5%. TDR mutations were recognized from the TDR consensus list. Transmission directionality, directness and inoculum size were inferred from tree topologies.

RESULTS

Both datasets concurred in the identification of seven transmission pairs and one cluster of three patients. With UDS, direction of transmission was inferred in four out of eight chains. Evidence for multiple founder viruses was found in two out of eight chains. No transmission of minority-resistant variants was evidenced. TDR mutations prevalence in protease and reverse transcriptase fragments was 4.3% with Sanger sequencing and 18.6% with UDS.

CONCLUSION

Although Sanger sequencing and UDS identified the same transmission chains, UDS provided additional information on founder viruses, direction of transmission and levels of TDR. Nevertheless, topology of clusters was not always consistent across gene fragments, calling for a cautious interpretation of the data. Moreover, unobserved intermediary links cannot be excluded. Phylogenetic analysis use as a forensic technique for HIV transmission investigations is risky.

Next-generation sequencing of HIV-1 single genome amplicons

The analysis of HIV-1 sequences has helped understand the viral molecular epidemiology, monitor the development of antiretroviral drug resistance, and design candidate vaccines. The introduction of single genome amplification (SGA) has been a major advancement in the field, allowing for the characterization of multiple sequences per patient while preserving linkage among polymorphisms in the same viral genome copy. Sequencing of SGA amplicons is performed by capillary Sanger sequencing, which presents low throughput, requires a high amount of template, and is highly sensitive to template/primer mismatching. In order to meet the increasing demand for HIV-1 SGA amplicon sequencing, we have developed a platform based on benchtop next-generation sequencing (NGS) (IonTorrent) accompanied by a bioinformatics pipeline capable of running on computer resources commonly available at research laboratories. During assay validation, the NGS-based sequencing of 10 HIV-1 env SGA amplicons was fully concordant with Sanger sequencing. The field test was conducted on plasma samples from 10 US Navy and Marine service members with recent HIV-1 infection (sampling interval: 2005–2010; plasma viral load: 5,884–194,984 copies/ml). The NGS analysis of 101 SGA amplicons (median: 10 amplicons/individual) showed within-individual viral sequence profiles expected in individuals at this disease stage, including individuals with highly homogeneous quasispecies, individuals with two highly homogeneous viral lineages, and individuals with heterogeneous viral populations. In a scalability assessment using the Ion Chef automated system, 41/43 tested env SGA amplicons (95%) multiplexed on a single Ion 318 chip showed consistent gene-wide coverage >50×. With lower sample requirements and higher throughput, this approach is suitable to support the increasing demand for high-quality and cost-effective HIV-1 sequences in fields such as molecular epidemiology, and development of preventive and therapeutic strategies.

Evolution of HIV-1 within untreated individuals and at the population scale in Uganda

HIV-1 undergoes multiple rounds of error-prone replication between transmission events, resulting in diverse viral populations within and among individuals. In addition, the virus experiences different selective pressures at multiple levels: during the course of infection, at transmission, and among individuals. Disentangling how these evolutionary forces shape the evolution of the virus at the population scale is important for understanding pathogenesis, how drug- and immune-escape variants are likely to spread in populations, and the development of preventive vaccines. To address this, we deep-sequenced two regions of the HIV-1 genome (p24 and gp41) from 34 longitudinally-sampled untreated individuals from Rakai District in Uganda, infected with subtypes A, D, and inter-subtype recombinants. This dataset substantially increases the availability of HIV-1 sequence data that spans multiple years of untreated infection, in particular for different geographical regions and viral subtypes. In line with previous studies, we estimated an approximately five-fold faster rate of evolution at the within-host compared to the population scale for both synonymous and nonsynonymous substitutions, and for all subtypes. We determined the extent to which this mismatch in evolutionary rates can be explained by the evolution of the virus towards population-level consensus, or the transmission of viruses similar to those that establish infection within individuals. Our findings indicate that both processes are likely to be important.

Identification of Unequally Represented Founder Viruses Among Tissues in Very Early SIV Rectal Transmission

Characterizing the transmitted/founder (T/F) viruses of multi-variant SIV infection may shed new light on the understanding of mucosal transmission. We intrarectally inoculated six Chinese rhesus macaques with a single high dose of SIVmac251 (3.1 × 10⁴ TCID₅₀) and obtained 985 full-length env sequences from multiple tissues at 6 and 10 days post-infection by single genome amplification (SGA). All 6 monkeys were infected with a range of 2 to 8 T/F viruses and the dominant variants from the inoculum were still dominant in different tissues from each monkey. Interestingly, our data showed that a cluster of rare T/F viruses was unequally represented in different tissues. This cluster of rare T/F viruses phylogenetically related to the non-dominant SIV variants in the inoculum and was not detected in any rectum tissues, but could be identified in the descending colon, jejunum, spleen, or plasma. In 2 out of 6 macaques, identical SIVmac251 variants belonging to this cluster were detected simultaneously in descending colon/jejunum and the inoculum. We also demonstrated that the average CG dinucleotide frequency of these rare T/F viruses found in tissues, as well as non-dominant variants in the inoculum, was significantly higher than the dominant T/F viruses in tissues and the inoculum. Collectively, these findings suggest that descending colon/jejunum might be more susceptible than rectum to SIV in the very early phase of infection. And host CG suppression, which was previously shown to inhibit HIV replication in vitro, may also contribute to the bottleneck selection during in vivo transmission.

Time Intervals in Sequence Sampling, Not Data Modifications, Have a Major Impact on Estimates of HIV Escape Rates

The ability of human immunodeficiency virus (HIV) to avoid recognition by humoral and cellular immunity (viral escape) is well-documented, but the strength of the immune response needed to cause such a viral escape remains poorly quantified. Several previous studies observed a more rapid escape of HIV from CD8 T cell responses in the acute phase of infection compared to chronic infection. The rate of HIV escape was estimated with the help of simple mathematical models, and results were interpreted to suggest that CD8 T cell responses causing escape in acute HIV infection may be more efficient at killing virus-infected cells than responses that cause escape in chronic infection, or alternatively, that early escapes occur in epitopes mutations in which there is minimal fitness cost to the virus. However, these conclusions were challenged on several grounds, including linkage and interference of multiple escape mutations due to a low population size and because of potential issues associated with modifying the data to estimate escape rates. Here we use a sampling method which does not require data modification to show that previous results on the decline of the viral escape rate with time since infection remain unchanged. However, using this method we also show that estimates of the escape rate are highly sensitive to the time interval between measurements, with longer intervals biasing estimates of the escape rate downwards. Our results thus suggest that data modifications for early and late escapes were not the primary reason for the observed decline in the escape rate with time since infection. However, longer sampling periods for escapes in chronic infection strongly influence estimates of the escape rate. More frequent sampling of viral sequences in chronic infection may improve our understanding of factors influencing the rate of HIV escape from CD8 T cell responses.

Kinetics of HIV-Specific CTL Responses Plays a Minimal Role in Determining HIV Escape Dynamics

Cytotoxic T lymphocytes (CTLs) have been suggested to play an important role in controlling human immunodeficiency virus (HIV-1 or simply HIV) infection. HIV, due to its high mutation rate, can evade recognition of T cell responses by generating escape variants that cannot be recognized by HIV-specific CTLs. Although HIV escape from CTL responses has been well documented, factors contributing to the timing and the rate of viral escape from T cells have not been fully elucidated. Fitness costs associated with escape and magnitude of the epitope-specific T cell response are generally considered to be the key in determining timing of HIV escape. Several previous analyses generally ignored the kinetics of T cell responses in predicting viral escape by either considering constant or maximal T cell response; several studies also considered escape from different T cell responses to be independent. Here, we focus our analysis on data from two patients from a recent study with relatively frequent measurements of both virus sequences and HIV-specific T cell response to determine impact of CTL kinetics on viral escape. In contrast with our expectation, we found that including temporal dynamics of epitope-specific T cell response did not improve the quality of fit of different models to escape data. We also found that for well-sampled escape data, the estimates of the model parameters including T cell killing efficacy did not strongly depend on the underlying model for escapes: models assuming independent, sequential, or concurrent escapes from multiple CTL responses gave similar estimates for CTL killing efficacy. Interestingly, the model assuming sequential escapes (i.e., escapes occurring along a defined pathway) was unable to accurately describe data on escapes occurring rapidly within a short-time window, suggesting that some of model assumptions must be violated for such escapes. Our results thus suggest that the current sparse measurements of temporal CTL dynamics in blood bear little quantitative information to improve predictions of HIV escape kinetics. More frequent measurements using more sensitive techniques and sampling in secondary lymphoid tissues may allow to better understand whether and how CTL kinetics impacts viral escape.

HIV-1 and hepatitis C virus selection bottleneck in Chinese people who inject drugs

OBJECTIVES

For both HIV-1 and hepatitis C virus (HCV), assessing the stringency of the transmission process is a scientific priority. Enumerations of transmitted/founder (TF) viruses have shown a strict transmission bottleneck in sexual transmission of HIV-1 and a wide range in the multiplicity of infection in HCV. Here, we aim to determine the stringency of parenteral transmission for HIV-1 and HCV in people who inject drugs (PWID).

DESIGN

We used molecular sequencing and several complementary analyses to enumerate the TF HIV-1 and HCV variants in a well described cohort of PWID in Xinjiang, China.

METHODS

We performed single genome sequencing of HIV-1 env and 5' half HCV genomes, then applied phylogenetic analysis and validated models of early virus diversification to enumerate TF viruses in 60 PWID. We used multivariate analysis to determine correlates of multivariant transmission (MVT).

RESULTS

We generated 1070 env region sequences from 33 HIV-1 early infected individuals and 773 5' half region sequences from 27 HCV early infected individuals. We found rates of MVT of 39 and 54%, respectively, for HIV-1 and HCV, with a limited range in the number of TF viruses in both infections. Behavioural characteristics suggested high-risk injection practices and lower risk sexual practices; we did not find an association between any specific behaviours and MVT.

CONCLUSION

MVT is frequent in parenteral transmission of both HIV-1 and HCV in Xinjiang PWID, indicating a less stringent transmission process than sexual transmission.

The impact of HIV-1 within-host evolution on transmission dynamics

The adaptive potential of HIV-1 is a vital mechanism to evade host immune responses and antiviral treatment. However, high evolutionary rates during persistent infection can impair transmission efficiency and alter disease progression in the new host, resulting in a delicate trade-off between within-host virulence and between-host infectiousness. This trade-off is visible in the disparity in evolutionary rates at within-host and between-host levels, and preferential transmission of ancestral donor viruses. Understanding the impact of within-host evolution for epidemiological studies is essential for the design of preventive and therapeutic measures. Herein, we review recent theoretical and experimental work that generated new insights into the complex link between within-host evolution and between-host fitness, revealing temporal and selective processes underlying the structure and dynamics of HIV-1 transmission.

Correction: Rare HIV-1 transmitted/founder lineages identified by deep viral sequencing contribute to rapid shifts in dominant quasispecies during acute and early infection

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