Whole genome deep sequencing of HIV-1 reveals the impact of early minor variants upon immune recognition during acute infection

Addressing pandemic-wide systematic errors in the SARS-CoV-2 phylogeny

The SARS-CoV-2 genome occupies a unique place in infection biology - it is the most highly sequenced genome on earth (making up over 20% of public sequencing datasets) with fine scale information on sampling date and geography, and has been subject to unprecedented intense analysis. As a result, these phylogenetic data are an incredibly valuable resource for science and public health. However, the vast majority of the data was sequenced by tiling amplicons across the full genome, with amplicon schemes that changed over the pandemic as mutations in the viral genome interacted with primer binding sites. In combination with the disparate set of genome assembly workflows and lack of consistent quality control (QC) processes, the current genomes have many systematic errors that have evolved with the virus and amplicon schemes. These errors have significant impacts on the phylogeny, and therefore over the last few years, many thousands of hours of researchers time has been spent in "eyeballing" trees, looking for artefacts, and then patching the tree. Given the huge value of this dataset, we therefore set out to reprocess the complete set of public raw sequence data in a rigorous amplicon-aware manner, and build a cleaner phylogeny. Here we provide a global tree of 3,960,704 samples, built from a consistently assembled set of high quality consensus sequences from all available public data as of March 2023, viewable at https://viridian.taxonium.org. Each genome was constructed using a novel assembly tool called Viridian (https://github.com/iqbal-lab-org/viridian), developed specifically to process amplicon sequence data, eliminating artefactual errors and mask the genome at low quality positions. We provide simulation and empirical validation of the methodology, and quantify the improvement in the phylogeny. Phase 2 of our project will address the fact that the data in the public archives is heavily geographically biased towards the Global North. We therefore have contributed new raw data to ENA/SRA from many countries including Ghana, Thailand, Laos, Sri Lanka, India, Argentina and Singapore. We will incorporate these, along with all public raw data submitted between March 2023 and the current day, into an updated set of assemblies, and phylogeny. We hope the tree, consensus sequences and Viridian will be a valuable resource for researchers.

Properties and Mechanisms of Deletions, Insertions, and Substitutions in the Evolutionary History of SARS-CoV-2

SARS-CoV-2 has accumulated many mutations since its emergence in late 2019. Nucleotide substitutions leading to amino acid replacements constitute the primary material for natural selection. Insertions, deletions, and substitutions appear to be critical for coronavirus's macro- and microevolution. Understanding the molecular mechanisms of mutations in the mutational hotspots (positions, loci with recurrent mutations, and nucleotide context) is important for disentangling roles of mutagenesis and selection. In the SARS-CoV-2 genome, deletions and insertions are frequently associated with repetitive sequences, whereas C>U substitutions are often surrounded by nucleotides resembling the APOBEC mutable motifs. We describe various approaches to mutation spectra analyses, including the context features of RNAs that are likely to be involved in the generation of recurrent mutations. We also discuss the interplay between mutations and natural selection as a complex evolutionary trend. The substantial variability and complexity of pipelines for the reconstruction of mutations and the huge number of genomic sequences are major problems for the analyses of mutations in the SARS-CoV-2 genome. As a solution, we advocate for the development of a centralized database of predicted mutations, which needs to be updated on a regular basis.

A VRC13-like bNAb response is associated with complex escape pathways in HIV-1 envelope

The rational design of HIV-1 immunogens to trigger the development of broadly neutralizing antibodies (bNAbs) requires understanding the viral evolutionary pathways influencing this process. An acute HIV-1-infected individual exhibiting >50% plasma neutralization breadth developed neutralizing antibody specificities against the CD4-binding site (CD4bs) and V1V2 regions of Env gp120. Comparison of pseudoviruses derived from early and late autologous env sequences demonstrated the development of >2 log resistance to VRC13 but not to other CD4bs-specific bNAbs. Mapping studies indicated that the V3 and CD4-binding loops of Env gp120 contributed significantly to developing resistance to the autologous neutralizing response and that the CD4-binding loop (CD4BL) specifically was responsible for the developing resistance to VRC13. Tracking viral evolution during the development of this cross-neutralizing CD4bs response identified amino acid substitutions arising at only 4 of 11 known VRC13 contact sites (K282, T283, K421, and V471). However, each of these mutations was external to the V3 and CD4BL regions conferring resistance to VRC13 and was transient in nature. Rather, complete resistance to VRC13 was achieved through the cooperative expression of a cluster of single amino acid changes within and immediately adjacent to the CD4BL, including a T359I substitution, exchange of a potential N-linked glycosylation (PNLG) site to residue S362 from N363, and a P369L substitution. Collectively, our data characterize complex HIV-1 env evolution in an individual developing resistance to a VRC13-like neutralizing antibody response and identify novel VRC13-associated escape mutations that may be important to inducing VRC13-like bNAbs for lineage-based immunogens.IMPORTANCEThe pursuit of eliciting broadly neutralizing antibodies (bNAbs) through vaccination and their use as therapeutics remains a significant focus in the effort to eradicate HIV-1. Key to our understanding of this approach is a more extensive understanding of bNAb contact sites and susceptible escape mutations in HIV-1 envelope (env). We identified a broad neutralizer exhibiting VRC13-like responses, a non-germline restricted class of CD4-binding site antibody distinct from the well-studied VRC01-class. Through longitudinal envelope sequencing and Env-pseudotyped neutralization assays, we characterized a complex escape pathway requiring the cooperative evolution of four amino acid changes to confer complete resistance to VRC13. This suggests that VRC13-class bNAbs may be refractory to rapid escape and attractive for therapeutic applications. Furthermore, the identification of longitudinal viral changes concomitant with the development of neutralization breadth may help identify the viral intermediates needed for the maturation of VRC13-like responses and the design of lineage-based immunogens.

HIV reservoirs are dominated by genetically younger and clonally enriched proviruses

IMPORTANCE

Characterizing the human immunodeficiency virus (HIV) reservoir that endures despite antiretroviral therapy (ART) is critical to cure efforts. We observed that the oldest proviruses persisting during ART were exclusively defective, while intact proviruses (and rebound HIV) dated to nearer ART initiation. This helps explain why studies that sampled sub-genomic proviruses on-ART (which are largely defective) routinely found sequences dating to early infection, whereas those that sampled replication-competent HIV found almost none. Together with our findings that intact proviruses were more likely to be clonal, and that on-ART low-level/isolated viremia originated from proviruses of varying ages (including possibly defective ones), our observations indicate that (i) on-ART and rebound viremia can have distinct within-host origins, (ii) intact proviruses have shorter lifespans than grossly defective ones and thus depend more heavily on clonal expansion for persistence, and (iii) an HIV reservoir predominantly "dating" to near ART initiation will be substantially adapted to within-host pressures, complicating immune-based cure strategies.

Towards a post-pandemic future for global pathogen genome sequencing

Pathogen genome sequencing has become a routine part of our response to active outbreaks of infectious disease and should be an important part of our preparations for future epidemics. In this Essay, we discuss the innovations that have enabled routine pathogen genome sequencing, as well as how genome sequences can be used to understand and control the spread of infectious disease. We also explore the impact of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) pandemic on the field of pathogen genomics and outline the challenges we must address to further improve the utility of pathogen genome sequencing in the future.

Inspecting the interaction between human immunodeficiency virus and the immune system through genetic turnover

Chronic infections of the human immunodeficiency virus (HIV) create a very complex coevolutionary process, where the virus tries to escape the continuously adapting host immune system. Quantitative details of this process are largely unknown and could help in disease treatment and vaccine development. Here we study a longitudinal dataset of ten HIV-infected people, where both the B-cell receptors and the virus are deeply sequenced. We focus on simple measures of turnover, which quantify how much the composition of the viral strains and the immune repertoire change between time points. At the single-patient level, the viral-host turnover rates do not show any statistically significant correlation, however, they correlate if one increases the amount of statistics by aggregating the information across patients. We identify an anti-correlation: large changes in the viral pool composition come with small changes in the B-cell receptor repertoire. This result seems to contradict the naïve expectation that when the virus mutates quickly, the immune repertoire needs to change to keep up. However, a simple model of antagonistically evolving populations can explain this signal. If it is sampled at intervals comparable with the sweep time, one population has had time to sweep while the second cannot start a counter-sweep, leading to the observed anti-correlation. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

A review of current effective COVID-19 testing methods and quality control

COVID-19 is a highly infectious disease caused by the SARS-CoV-2 virus, which primarily affects the respiratory system and can lead to severe illness. The virus is extremely contagious, early and accurate diagnosis of SARS-CoV-2 is crucial to contain its spread, to provide prompt treatment, and to prevent complications. Currently, the reverse transcriptase polymerase chain reaction (RT-PCR) is considered to be the gold standard for detecting COVID-19 in its early stages. In addition, loop-mediated isothermal amplification (LMAP), clustering rule interval short palindromic repeats (CRISPR), colloidal gold immunochromatographic assay (GICA), computed tomography (CT), and electrochemical sensors are also common tests. However, these different methods vary greatly in terms of their detection efficiency, specificity, accuracy, sensitivity, cost, and throughput. Besides, most of the current detection methods are conducted in central hospitals and laboratories, which is a great challenge for remote and underdeveloped areas. Therefore, it is essential to review the advantages and disadvantages of different COVID-19 detection methods, as well as the technology that can enhance detection efficiency and improve detection quality in greater details.

HIV reservoirs are dominated by genetically younger and clonally enriched proviruses

In order to cure HIV, we need to better understand the within-host evolutionary origins of the small reservoir of genome-intact proviruses that persists within infected cells during antiretroviral therapy (ART). Most prior studies on reservoir evolutionary dynamics however did not discriminate genome-intact proviruses from the vast background of defective ones. We reconstructed within-host pre-ART HIV evolutionary histories in six individuals and leveraged this information to infer the ages of intact and defective proviruses sampled after an average >9 years on ART, along with the ages of rebound and low-level/isolated viremia occurring during this time. We observed that the longest-lived proviruses persisting on ART were exclusively defective, usually due to large deletions. In contrast, intact proviruses and rebound HIV exclusively dated to the years immediately preceding ART. These observations are consistent with genome-intact proviruses having shorter lifespans, likely due to the cumulative risk of elimination following viral reactivation and protein production. Consistent with this, intact proviruses (and those with packaging signal defects) were three times more likely to be genetically identical compared to other proviral types, highlighting clonal expansion as particularly important in ensuring their survival. By contrast, low-level/isolated viremia sequences were genetically heterogeneous and sometimes ancestral, where viremia may have originated from defective proviruses. Results reveal that the HIV reservoir is dominated by clonally-enriched and genetically younger sequences that date to the untreated infection period when viral populations had been under within-host selection pressures for the longest duration. Knowledge of these qualities may help focus strategies for reservoir elimination.

Inferring the multiplicity of founder variants initiating HIV-1 infection: a systematic review and individual patient data meta-analysis

BACKGROUND

HIV-1 infections initiated by multiple founder variants are characterised by a higher viral load and a worse clinical prognosis than those initiated with single founder variants, yet little is known about the routes of exposure through which transmission of multiple founder variants is most probable. Here we used individual patient data to calculate the probability of multiple founders stratified by route of HIV exposure and study methodology.

METHODS

We conducted a systematic review and meta-analysis of studies that estimated founder variant multiplicity in HIV-1 infection, searching MEDLINE, Embase, and Global Health databases for papers published between Jan 1, 1990, and Sept 14, 2020. Eligible studies must have reported original estimates of founder variant multiplicity in people with acute or early HIV-1 infections, have clearly detailed the methods used, and reported the route of exposure. Studies were excluded if they reported data concerning people living with HIV-1 who had known or suspected superinfection, who were documented as having received pre-exposure prophylaxis, or if the transmitting partner was known to be receiving antiretroviral treatment. Individual patient data were collated from all studies, with authors contacted if these data were not publicly available. We applied logistic meta-regression to these data to estimate the probability that an HIV infection is initiated by multiple founder variants. We calculated a pooled estimate using a random effects model, subsequently stratifying this estimate across exposure routes in a univariable analysis. We then extended our model to adjust for different study methods in a multivariable analysis, recalculating estimates across the exposure routes. This study is registered with PROSPERO, CRD42020202672.

FINDINGS

We included 70 publications in our analysis, comprising 1657 individual patients. Our pooled estimate of the probability that an infection is initiated by multiple founder variants was 0·25 (95% CI 0·21-0·29), with moderate heterogeneity (Q=132·3, p<0·0001, I²=64·2%). Our multivariable analysis uncovered differences in the probability of multiple variant infection by exposure route. Relative to a baseline of male-to-female transmission, the predicted probability for female-to-male multiple variant transmission was significantly lower at 0·13 (95% CI 0·08-0·20), and the probabilities were significantly higher for transmissions in people who inject drugs (0·37 [0·24-0·53]) and men who have sex with men (0·30 [0·33-0·40]). There was no significant difference in the probability of multiple variant transmission between male-to-female transmission (0·21 [0·14-0·31]), post-partum transmission (0·18 [0·03-0·57]), pre-partum transmission (0·17 [0·08-0·33]), and intra-partum transmission (0·27 [0·14-0·45]).

INTERPRETATION

We identified that transmissions in people who inject drugs and men who have sex with men are significantly more likely to result in an infection initiated by multiple founder variants, and female-to-male infections are significantly less probable. Quantifying how the routes of HIV infection affect the transmission of multiple variants allows us to better understand how the evolution and epidemiology of HIV-1 determine clinical outcomes.

FUNDING

Medical Research Council Precision Medicine Doctoral Training Programme and a European Research Council Starting Grant.

Advances in Molecular Genetics Enabling Studies of Highly Pathogenic RNA Viruses

Experimental work with viruses that are highly pathogenic for humans and animals requires specialized Biosafety Level 3 or 4 facilities. Such pathogens include some spectacular but also rather seldomly studied examples such as Ebola virus (requiring BSL-4), more wide-spread and commonly studied viruses such as HIV, and the most recent example, SARS-CoV-2, which causes COVID-19. A common characteristic of these virus examples is that their genomes consist of single-stranded RNA, which requires the conversion of their genomes into a DNA copy for easy manipulation; this can be performed to study the viral life cycle in detail, develop novel therapies and vaccines, and monitor the disease course over time for chronic virus infections. We summarize the recent advances in such new genetic applications for RNA viruses in Switzerland over the last 25 years, from the early days of the HIV/AIDS epidemic to the most recent developments in research on the SARS-CoV-2 coronavirus. We highlight game-changing collaborative efforts between clinical and molecular disciplines in HIV research on the path to optimal clinical disease management. Moreover, we summarize how the modern technical evolution enabled the molecular studies of emerging RNA viruses, confirming that Switzerland is at the forefront of SARS-CoV-2 research and potentially other newly emerging viruses.

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.

Immunoescape of HIV-1 in Env-EL9 CD8 + T cell response restricted by HLA-B*14:02 in a Non progressor who lost twenty-seven years of HIV-1 control

Background

Long-Term Non-Progressors (LTNPs) are untreated Human Immunodeficiency virus type 1 (HIV-1) infected individuals able to control disease progression for prolonged periods. However, the LTNPs status is temporary, as viral load increases followed by decreases in CD4 + T-cell counts. Control of HIV-1 infection in LTNPs viremic controllers, have been associated with effective immunodominant HIV-1 Gag-CD8 + T-cell responses restricted by protective HLA-B alleles. Individuals carrying HLA-B*14:02 control HIV-1 infection is related to an immunodominant Env-CD8 + T-cell response. Limited data are available on the contribution of HLA-B*14:02 CD8 + T -cells in LTNPs.

Results

In this study, we performed a virological and immunological detailed analysis of an HLA-B*14:02 LNTP individual that lost viral control (LVC) 27 years after HIV-1 diagnosis. We analysed viral evolution and immune escape in HLA-B*14:02 restricted CD8 + T -cell epitopes and identified viral evolution at the Env-EL9 epitope selecting the L592R mutation. By IFN-γ ELISpot and immune phenotype, we characterized HLA- B*14:02 HIV-1 CD8 + T cell responses targeting, Gag-DA9 and Env-EL9 epitopes before and after LVC. We observed an immunodominant response against the Env-EL9 epitope and a decreased of the CD8 T + cell response over time with LVC. Loss of Env-EL9 responses was concomitant with selecting K588R + L592R mutations at Env-EL9. Finally, we evaluated the impact of Env-EL9 escape mutations on HIV-1 infectivity and Env protein structure. The K588R + L592R escape variant was directly related to HIV-1 increase replicative capacity and stability of Env at the LVC.

Conclusions

These findings support the contribution of immunodominant Env-EL9 CD8 + T-cell responses and the imposition of immune escape variants with higher replicative capacity associated with LVC in this LNTP. These data highlight the importance of Env-EL9 specific-CD8 + T-cell responses restricted by the HLA-B*14:02 and brings new insights into understanding long-term HIV-1 control mediated by Env mediated CD8 + T-cell responses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12977-022-00591-7.

Evolution during primary HIV infection does not require adaptive immune selection

Modern HIV research depends crucially on both viral sequencing and population measurements. To directly link mechanistic biological processes and evolutionary dynamics during HIV infection, we developed multiple within-host phylodynamic models of HIV primary infection for comparative validation against viral load and evolutionary dynamics data. The optimal model of primary infection required no positive selection, suggesting that the host adaptive immune system reduces viral load but surprisingly does not drive observed viral evolution. Rather, the fitness (infectivity) of mutant variants is drawn from an exponential distribution in which most variants are slightly less infectious than their parents (nearly neutral evolution). This distribution was not largely different from either in vivo fitness distributions recorded beyond primary infection or in vitro distributions that are observed without adaptive immunity, suggesting the intrinsic viral fitness distribution may drive evolution. Simulated phylogenetic trees also agree with independent data and illuminate how phylogenetic inference must consider viral and immune-cell population dynamics to gain accurate mechanistic insights.

Analyzing rare mutations in metagenomes assembled using long and accurate reads

The advent of long and accurate "HiFi" reads has greatly improved our ability to generate complete metagenome-assembled genomes (MAGs), enabling "complete metagenomics" studies that were nearly impossible to conduct with short reads. In particular, HiFi reads simplify the identification and phasing of mutations in MAGs: It is increasingly feasible to distinguish between positions that are prone to mutations and positions that rarely ever mutate, and to identify co-occurring groups of mutations. However, the problems of identifying rare mutations in MAGs, estimating the false-discovery rate (FDR) of these identifications, and phasing identified mutations remain open in the context of HiFi data. We present strainFlye, a pipeline for the FDR-controlled identification and analysis of rare mutations in MAGs assembled using HiFi reads. We show that deep HiFi sequencing has the potential to reveal and phase tens of thousands of rare mutations in a single MAG, identify hotspots and coldspots of these mutations, and detail MAGs' growth dynamics.

OUP accepted manuscript

A detailed understanding of how and when severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission occurs is crucial for designing effective prevention measures. Other than contact tracing, genome sequencing provides information to help infer who infected whom. However, the effectiveness of the genomic approach in this context depends on both (high enough) mutation and (low enough) transmission rates. Today, the level of resolution that we can obtain when describing SARS-CoV-2 outbreaks using just genomic information alone remains unclear. In order to answer this question, we sequenced forty-nine SARS-CoV-2 patient samples from ten local clusters in NW Spain for which partial epidemiological information was available and inferred transmission history using genomic variants. Importantly, we obtained high-quality genomic data, sequencing each sample twice and using unique barcodes to exclude cross-sample contamination. Phylogenetic and cluster analyses showed that consensus genomes were generally sufficient to discriminate among independent transmission clusters. However, levels of intrahost variation were low, which prevented in most cases the unambiguous identification of direct transmission events. After filtering out recurrent variants across clusters, the genomic data were generally compatible with the epidemiological information but did not support specific transmission events over possible alternatives. We estimated the effective transmission bottleneck size to be one to two viral particles for sample pairs whose donor–recipient relationship was likely. Our analyses suggest that intrahost genomic variation in SARS-CoV-2 might be generally limited and that homoplasy and recurrent errors complicate identifying shared intrahost variants. Reliable reconstruction of direct SARS-CoV-2 transmission based solely on genomic data seems hindered by a slow mutation rate, potential convergent events, and technical artifacts. Detailed contact tracing seems essential in most cases to study SARS-CoV-2 transmission at high resolution.

Space and Genotype-Dependent Virus Distribution during Infection Progression

The paper is devoted to a nonlocal reaction-diffusion equation describing the development of viral infection in tissue, taking into account virus distribution in the space of genotypes, the antiviral immune response, and natural genotype-dependent virus death. It is shown that infection propagates as a reaction-diffusion wave. In some particular cases, the 2D problem can be reduced to a 1D problem by separation of variables, allowing for proof of wave existence and stability. In general, this reduction provides an approximation of the 2D problem by a 1D problem. The analysis of the reduced problem allows us to determine how viral load and virulence depend on genotype distribution, the strength of the immune response, and the level of immunity.

HIV Proviral Burden, Genetic Diversity, and Dynamics in Viremic Controllers Who Subsequently Initiated Suppressive Antiretroviral Therapy

Curing HIV will require eliminating the reservoir of integrated, replication-competent proviruses that persist despite antiretroviral therapy (ART). Understanding the burden, genetic diversity, and longevity of persisting proviruses in diverse individuals with HIV is critical to this goal, but these characteristics remain understudied in some groups. Among them are viremic controllers—individuals who naturally suppress HIV to low levels but for whom therapy is nevertheless recommended. We reconstructed within-host HIV evolutionary histories from longitudinal single-genome amplified viral sequences in four viremic controllers who eventually initiated ART and used this information to characterize the age and diversity of proviruses persisting on therapy. We further leveraged these within-host proviral age distributions to estimate rates of proviral turnover prior to ART. This is an important yet understudied metric, since pre-ART proviral turnover dictates reservoir composition at ART initiation (and thereafter), which is when curative interventions, once developed, would be administered. Despite natural viremic control, all participants displayed significant within-host HIV evolution pretherapy, where overall on-ART proviral burden and diversity broadly reflected the extent of viral replication and diversity pre-ART. Consistent with recent studies of noncontrollers, the proviral pools of two participants were skewed toward sequences that integrated near ART initiation, suggesting dynamic proviral turnover during untreated infection. In contrast, proviruses recovered from the other two participants dated to time points that were more evenly spread throughout infection, suggesting slow or negligible proviral decay following deposition. HIV cure strategies will need to overcome within-host proviral diversity, even in individuals who naturally controlled HIV replication before therapy.

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 (N_e ) 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 N_e 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 (N_e ) 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 N_e 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.

Functional impairment of HIV-specific CD8+ T cells precedes aborted spontaneous control of viremia

Spontaneous control of HIV infection has been repeatedly linked to antiviral CD8+ T cells but is not always permanent. To address mechanisms of durable and aborted control of viremia, we evaluated immunologic and virologic parameters longitudinally among 34 HIV-infected subjects with differential outcomes. Despite sustained recognition of autologous virus, HIV-specific proliferative and cytolytic T cell effector functions became selectively and intrinsically impaired prior to aborted control. Longitudinal transcriptomic profiling of functionally impaired HIV-specific CD8+ T cells revealed altered expression of genes related to activation, cytokine-mediated signaling, and cell cycle regulation, including increased expression of the antiproliferative transcription factor KLF2 but not of genes associated with canonical exhaustion. Lymphoid HIV-specific CD8+ T cells also exhibited poor functionality during aborted control relative to durable control. Our results identify selective functional impairment of HIV-specific CD8+ T cells as prognostic of impending aborted HIV control, with implications for clinical monitoring and immunotherapeutic strategies.

Identification of Genetically Related HCV Infections Among Self-Described Injecting Partnerships

BACKGROUND

The current opioid epidemic across the United States has fueled a surge in the rate of new hepatitis C virus (HCV) infections among young persons who inject drugs (PWIDs). Paramount to interrupting transmission is targeting these high-risk populations and understanding the underlying network structures facilitating transmission within these communities.

METHODS

Deep sequencing data were obtained for 52 participants from 32 injecting partnerships enrolled in the U-Find-Out (UFO) Partner Study, which is a prospective study of self-described injecting dyad partnerships from a large community-based study of HCV infection in young adult PWIDs from San Francisco. Phylogenetically linked transmission events were identified using traditional genetic-distance measures and viral deep sequence phylogenies reconstructed to determine the statistical support of inferences and the direction of transmission within partnerships.

RESULTS

Using deep sequencing data, we found that 12 of 32 partnerships were genetically similar and clustered. Three additional phylogenetic clusters were found describing novel putative transmission links outside of the injecting relationship. Transmission direction was inferred correctly for 5 partnerships with the incorrect transmission direction inferred in more than 50% of cases. Notably, we observed that phylogenetic linkage was most often associated with a lower number of network partners and involvement in a sexual relationship.

CONCLUSIONS

Deep sequencing of HCV among self-described injecting partnerships demonstrates that the majority of transmission events originate from outside of the injecting partnership. Furthermore, these findings caution that phylogenetic methods may be unable to routinely infer the direction of transmission among PWIDs especially when transmission events occur in rapid succession within high-risk networks.

Differentiation of exhausted CD8+ T cells after termination of chronic antigen stimulation stops short of achieving functional T cell memory

T cell exhaustion is associated with failure to clear chronic infections and malignant cells. Defining the molecular mechanisms of T cell exhaustion and reinvigoration is essential to improving immunotherapeutic modalities. Here we confirmed pervasive phenotypic, functional, and transcriptional differences between memory and exhausted antigen-specific CD8+ T cells in human hepatitis C virus (HCV) infection before and after treatment. After viral cure, phenotypic changes in clonally stable exhausted T cell populations suggested differentiation towards a memory-like profile. However, functionally, the cells showed little improvement and critical transcriptional regulators remained in the exhaustion state. Notably, T cells from chronic HCV infection that were exposed to antigen for less time because of viral escape mutations were functionally and transcriptionally more similar to memory T cells from spontaneously resolved HCV infection. Thus, T cell stimulation duration impacts exhaustion recovery, with antigen removal after long-term exhaustion being insufficient for development of functional T cell memory.

Hepatitis C Virus Transmission Clusters in Public Health and Correctional Settings, Wisconsin, USA, 2016-20171

Ending the hepatitis C virus (HCV) epidemic requires stopping transmission among networks of persons who inject drugs. Identifying transmission networks by using genomic epidemiology may inform community responses that can quickly interrupt transmission. We retrospectively identified HCV RNA–positive specimens corresponding to 459 persons in settings that use the state laboratory, including correctional facilities and syringe services programs, in Wisconsin, USA, during 2016–2017. We conducted next-generation sequencing of HCV and analyzed it for phylogenetic linkage by using the Centers for Disease Control and Prevention Global Hepatitis Outbreak Surveillance Technology platform. Analysis showed that 126 persons were linked across 42 clusters. Phylogenetic clustering was higher in rural communities and associated with female sex and younger age among rural residents. These data highlight that HCV transmission could be reduced by expanding molecular-based surveillance strategies to rural communities affected by the opioid crisis.

Cytotoxic CD8+ T Cells Expressing CXCR5 Are Detectable in HIV-1 Elite Controllers After Prolonged In Vitro Peptide Stimulation

Antiretroviral therapy (ART) is not curative as HIV-1 persists in long-lived viral reservoirs. Consequently, patients are dependent on life-long drug adherence with possible side effects. To overcome these limitations strategies of a functional cure aim at ART free viral remission. In this study, we sought to identify detailed subsets of anti-viral CD8⁺ T cell immunity linked to natural long-term control of HIV-1 infection. Here, we analyzed HIV controllers and ART suppressed progressors for in vitro viral suppressive capacity (VSC) at baseline and after peptide stimulation. Functional properties and phenotypes of CD8⁺ T cells were assessed by IFN-γ ELISPOT and 18 color flow cytometry. HIV controllers showed significantly increased suppression at baseline as well as after peptide stimulation. IFN-γ secretion and the proliferation marker Ki67 positively correlated with VSC. Moreover, the detailed phenotype of three distinct multifunctional memory CD8⁺ T cell subsets were specific traits of HIV controllers of which two correlated convincingly with VSC. Our results underline the importance of multifunctional CD8⁺ T cell responses during natural control. Especially the role of CXCR5 expressing cytotoxic subsets emphasizes potential surveillance in sites of reservoir persistence and demand further study.

Epidemiological data analysis of viral quasispecies in the next-generation sequencing era

The unprecedented coverage offered by next-generation sequencing (NGS) technology has facilitated the assessment of the population complexity of intra-host RNA viral populations at an unprecedented level of detail. Consequently, analysis of NGS datasets could be used to extract and infer crucial epidemiological and biomedical information on the levels of both infected individuals and susceptible populations, thus enabling the development of more effective prevention strategies and antiviral therapeutics. Such information includes drug resistance, infection stage, transmission clusters and structures of transmission networks. However, NGS data require sophisticated analysis dealing with millions of error-prone short reads per patient. Prior to the NGS era, epidemiological and phylogenetic analyses were geared toward Sanger sequencing technology; now, they must be redesigned to handle the large-scale NGS datasets and properly model the evolution of heterogeneous rapidly mutating viral populations. Additionally, dedicated epidemiological surveillance systems require big data analytics to handle millions of reads obtained from thousands of patients for rapid outbreak investigation and management. We survey bioinformatics tools analyzing NGS data for (i) characterization of intra-host viral population complexity including single nucleotide variant and haplotype calling; (ii) downstream epidemiological analysis and inference of drug-resistant mutations, age of infection and linkage between patients; and (iii) data collection and analytics in surveillance systems for fast response and control of outbreaks.

Validation of multiplex PCR sequencing assay of SIV

Background

The generation of accurate and reproducible viral sequence data is necessary to understand the diversity present in populations of RNA viruses isolated from clinical samples. While various sequencing methods are available, they often require high quality templates and high viral titer to ensure reliable data.

Methods

We modified a multiplex PCR and sequencing approach to characterize populations of simian immunodeficiency virus (SIV) isolated from nonhuman primates. We chose this approach with the aim of reducing the number of required input templates while maintaining fidelity and sensitivity. We conducted replicate sequencing experiments using different numbers of quantified viral RNA (vRNA) or viral cDNA as input material. We performed assays with clonal SIVmac239 to detect false positives, and we mixed SIVmac239 and a variant with 24 point mutations (SIVmac239-24X) to measure variant detection sensitivity.

Results

We found that utilizing a starting material of quantified viral cDNA templates had a lower rate of false positives and increased reproducibility when compared to that of quantified vRNA templates. This study identifies the importance of rigorously validating deep sequencing methods and including replicate samples when using a new method to characterize low frequency variants in a population with a small number of templates.

Conclusions

Because the need to generate reproducible and accurate sequencing data from diverse viruses from low titer samples, we modified a multiplex PCR and sequencing approach to characterize SIV from populations from non-human primates. We found that increasing starting template numbers increased the reproducibility and decreased the number of false positives identified, and this was further seen when cDNA was used as a starting material. Ultimately, we highlight the importance of vigorously validating methods to prevent overinterpretation of low frequency variants in a sample.

Physical 'strength' of the multi-protein chain connecting immune cells: Does the weakest link limit antibody affinity maturation?: The weakest link in the multi-protein chain facilitating antigen acquisition by B cells in germinal centres limits antibody affinity maturation

The affinities of antibodies (Abs) for their target antigens (Ags) gradually increase in vivo following an infection or vaccination, but reach saturation at values well below those realisable in vitro. This 'affinity ceiling' could in many cases restrict our ability to fight infections and compromise vaccines. What determines the affinity ceiling has been an unresolved question for decades. Here, we argue that it arises from the strength of the chain of protein complexes that is pulled by B cells during the process of Ag acquisition. The affinity ceiling is determined by the strength of the weakest link in the chain. We identify the weakest link and show that the resulting affinity ceiling can explain the Ab affinities realized in vivo, providing a conceptual understanding of Ab affinity maturation. We explore plausible evolutionary underpinnings of the affinity ceiling, examine supporting evidence and alternative hypotheses and discuss implications for vaccination strategies.

Molecular epidemiology in the HIV and SARS-CoV-2 pandemics

PURPOSE OF REVIEW

The aim of this review was to compare and contrast the application of molecular epidemiology approaches for the improved management and understanding of the HIV versus SARS-CoV-2 epidemics.

RECENT FINDINGS

Molecular biology approaches, including PCR and whole genome sequencing (WGS), have become powerful tools for epidemiological investigation. PCR approaches form the basis for many high-sensitivity diagnostic tests and can supplement traditional contact tracing and surveillance strategies to define risk networks and transmission patterns. WGS approaches can further define the causative agents of disease, trace the origins of the pathogen, and clarify routes of transmission. When coupled with clinical datasets, such as electronic medical record data, these approaches can investigate co-correlates of disease and pathogenesis. In the ongoing HIV epidemic, these approaches have been effectively deployed to identify treatment gaps, transmission clusters and risk factors, though significant barriers to rapid or real-time implementation remain critical to overcome. Likewise, these approaches have been successful in addressing some questions of SARS-CoV-2 transmission and pathogenesis, but the nature and rapid spread of the virus have posed additional challenges.

SUMMARY

Overall, molecular epidemiology approaches offer unique advantages and challenges that complement traditional epidemiological tools for the improved understanding and management of epidemics.

Analytical validity of nanopore sequencing for rapid SARS-CoV-2 genome analysis

Viral whole-genome sequencing (WGS) provides critical insight into the transmission and evolution of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Long-read sequencing devices from Oxford Nanopore Technologies (ONT) promise significant improvements in turnaround time, portability and cost, compared to established short-read sequencing platforms for viral WGS (e.g., Illumina). However, adoption of ONT sequencing for SARS-CoV-2 surveillance has been limited due to common concerns around sequencing accuracy. To address this, here we perform viral WGS with ONT and Illumina platforms on 157 matched SARS-CoV-2-positive patient specimens and synthetic RNA controls, enabling rigorous evaluation of analytical performance. We report that, despite the elevated error rates observed in ONT sequencing reads, highly accurate consensus-level sequence determination was achieved, with single nucleotide variants (SNVs) detected at >99% sensitivity and >99% precision above a minimum ~60-fold coverage depth, thereby ensuring suitability for SARS-CoV-2 genome analysis. ONT sequencing also identified a surprising diversity of structural variation within SARS-CoV-2 specimens that were supported by evidence from short-read sequencing on matched samples. However, ONT sequencing failed to accurately detect short indels and variants at low read-count frequencies. This systematic evaluation of analytical performance for SARS-CoV-2 WGS will facilitate widespread adoption of ONT sequencing within local, national and international COVID-19 public health initiatives.

Quasispecies dynamics of a hepatitis E virus 4 from the feces and liver biopsy of an acute hepatitis E patient during virus clearance

The presence of a quasispecies in hepatitis E virus (HEV) infection has been documented, however, the implications of a quasispecies in HEV-host interaction are poorly understood. Here, we analyzed the whole genome sequences of a HEV 4d from the feces and liver biopsy of a patient during the icteric and convalescent phases in an acute hepatitis E infection. Viral RNAs were extracted, reversely transcribed and seven fragments encompassing the entire viral genome were amplified and cloned. By sequencing multiple colonies of each cloned viral genome amplicon with Sanger sequencing, we verified the existence of the HEV quasispecies or intra-host genetic variations within the fecal and liver biopsy samples. There were broader genetic variations in the HEV ORF1 region including the PCP, HPX, and RdRp regions during the convalescent phase whereas more genetic variations in the ORF2 P domain during the icteric phase. The quasispecies dynamics might reflect host immune pressure during viral clearance.

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.

Evaluation of haplotype callers for next-generation sequencing of viruses

Currently, the standard practice for assembling next-generation sequencing (NGS) reads of viral genomes is to summarize thousands of individual short reads into a single consensus sequence, thus confounding useful intra-host diversity information for molecular phylodynamic inference. It is hypothesized that a few viral strains may dominate the intra-host genetic diversity with a variety of lower frequency strains comprising the rest of the population. Several software tools currently exist to convert NGS sequence variants into haplotypes. Previous benchmarks of viral haplotype reconstruction programs used simulation scenarios that are useful from a mathematical perspective but do not reflect viral evolution and epidemiology. Here, we tested twelve NGS haplotype reconstruction methods using viral populations simulated under realistic evolutionary dynamics. We simulated coalescent-based populations that spanned known levels of viral genetic diversity, including mutation rates, sample size and effective population size, to test the limits of the haplotype reconstruction methods and to ensure coverage of predicted intra-host viral diversity levels (especially HIV-1). All twelve investigated haplotype callers showed variable performance and produced drastically different results that were mainly driven by differences in mutation rate and, to a lesser extent, in effective population size. Most methods were able to accurately reconstruct haplotypes when genetic diversity was low. However, under higher levels of diversity (e.g., those seen intra-host HIV-1 infections), haplotype reconstruction quality was highly variable and, on average, poor. All haplotype reconstruction tools, except QuasiRecomb and ShoRAH, greatly underestimated intra-host diversity and the true number of haplotypes. PredictHaplo outperformed, in regard to highest precision, recall, and lowest UniFrac distance values, the other haplotype reconstruction tools followed by CliqueSNV, which, given more computational time, may have outperformed PredictHaplo. Here, we present an extensive comparison of available viral haplotype reconstruction tools and provide insights for future improvements in haplotype reconstruction tools using both short-read and long-read technologies.

Validation of Variant Assembly Using HAPHPIPE with Next-Generation Sequence Data from Viruses

Next-generation sequencing (NGS) offers a powerful opportunity to identify low-abundance, intra-host viral sequence variants, yet the focus of many bioinformatic tools on consensus sequence construction has precluded a thorough analysis of intra-host diversity. To take full advantage of the resolution of NGS data, we developed HAplotype PHylodynamics PIPEline (HAPHPIPE), an open-source tool for the de novo and reference-based assembly of viral NGS data, with both consensus sequence assembly and a focus on the quantification of intra-host variation through haplotype reconstruction. We validate and compare the consensus sequence assembly methods of HAPHPIPE to those of two alternative software packages, HyDRA and Geneious, using simulated HIV and empirical HIV, HCV, and SARS-CoV-2 datasets. Our validation methods included read mapping, genetic distance, and genetic diversity metrics. In simulated NGS data, HAPHPIPE generated pol consensus sequences significantly closer to the true consensus sequence than those produced by HyDRA and Geneious and performed comparably to Geneious for HIV gp120 sequences. Furthermore, using empirical data from multiple viruses, we demonstrate that HAPHPIPE can analyze larger sequence datasets due to its greater computational speed. Therefore, we contend that HAPHPIPE provides a more user-friendly platform for users with and without bioinformatics experience to implement current best practices for viral NGS assembly than other currently available options.

The effect of variant interference on de novo assembly for viral deep sequencing

BACKGROUND

Viruses have high mutation rates and generally exist as a mixture of variants in biological samples. Next-generation sequencing (NGS) approaches have surpassed Sanger for generating long viral sequences, yet how variants affect NGS de novo assembly remains largely unexplored.

RESULTS

Our results from > 15,000 simulated experiments showed that presence of variants can turn an assembly of one genome into tens to thousands of contigs. This "variant interference" (VI) is highly consistent and reproducible by ten commonly-used de novo assemblers, and occurs over a range of genome length, read length, and GC content. The main driver of VI is pairwise identities between viral variants. These findings were further supported by in silico simulations, where selective removal of minor variant reads from clinical datasets allow the "rescue" of full viral genomes from fragmented contigs.

CONCLUSIONS

These results call for careful interpretation of contigs and contig numbers from de novo assembly in viral deep sequencing.

Analysis of sequence diversity and selection pressure in HIV-1 clade C gp41 from India

Evaluation of viral diversity is critical for the rational design of treatment modalities against Human immunodeficiency virus (HIV). Predominated by HIV-1 clade C (HIV-1C), the epidemic in India represents the third largest population infected with HIV-1 globally. Glycoprotein 41 (gp41) is critical for viral replication and is a target for the design of therapeutic strategies. However, documentation of viral diversity of gp41 gene in infected individuals from India remains limited. Present study employed high throughput sequencing to examine variation in gp41 amplicons generated from blood derived viruses in 24 HIV-1C infected individuals from Mumbai, India. Sequence diversity profiles were documented in different functional domains of gp41. Furthermore, through a meta-analysis approach, all reported gp41 sequences from India (N = 70) were compared with those from South Africa (N = 126), country with the largest HIV epidemic globally, also predominated by HIV-1C. A total of 44 positions displayed statistically significant differential (p < 0.05) Shannon entropy in the two regions. This comparison also identified 11 codon sites undergoing distinct selection, 8 of which remained differentially selected in an extended comparison of data from Asia (N = 137) and Africa(N = 383). Assessment of correlated mutation networks associated with differentially selected residues revealed common as well as distinct interaction networks. Furthermore, codon usage analysis revealed 17 differentially selected codons (Mann-Whitney test, p < 0.001) in Asia and Africa. Dissimilar trends in GC content across codon positions were also observed. In depth understanding of these divergent evolutionary signatures through extended analysis with larger data-sets would assist development of effective interventions being considered for HIV-1C.

Gp41-targeted antibodies restore infectivity of a fusion-deficient HIV-1 envelope glycoprotein

The HIV-1 envelope glycoprotein (Env) mediates viral entry via conformational changes associated with binding the cell surface receptor (CD4) and coreceptor (CCR5/CXCR4), resulting in subsequent fusion of the viral and cellular membranes. While the gp120 Env surface subunit has been extensively studied for its role in viral entry and evasion of the host immune response, the gp41 transmembrane glycoprotein and its role in natural infection are less well characterized. Here, we identified a primary HIV-1 Env variant that consistently supports >300% increased viral infectivity in the presence of autologous or heterologous HIV-positive plasma. However, in the absence of HIV-positive plasma, viruses with this Env exhibited reduced infectivity that was not due to decreased CD4 binding. Using Env chimeras and sequence analysis, we mapped this phenotype to a change Q563R, in the gp41 heptad repeat 1 (HR1) region. We demonstrate that Q563R reduces viral infection by disrupting formation of the gp41 six-helix bundle required for virus-cell membrane fusion. Intriguingly, antibodies that bind cluster I epitopes on gp41 overcome this inhibitory effect, restoring infectivity to wild-type levels. We further demonstrate that the Q563R change increases HIV-1 sensitivity to broadly neutralizing antibodies (bNAbs) targeting the gp41 membrane-proximal external region (MPER). In summary, we identify an HIV-1 Env variant with impaired infectivity whose Env functionality is restored through the binding of host antibodies. These data contribute to our understanding of gp41 residues involved in membrane fusion and identify a mechanism by which host factors can alleviate a viral defect.

Clinical and evolutionary consequences of HIV adaptation to HLA: implications for vaccine and cure

PURPOSE OF REVIEW

The purpose of this review is to summarize recent advances in our understanding of HIV adaptation to human leukocyte antigen (HLA)-associated immune pressures and its relevance to HIV prevention and cure research.

RECENT FINDINGS

Recent research has confirmed that HLA is a major driver of individual and population-level HIV evolution, that HIV strains are adapting to the immunogenetic profiles of the different human ethnic groups in which they circulate, and that HIV adaptation has substantial clinical and immunologic consequences. As such, adaptation represents a major challenge to HIV prevention and cure. At the same time, there are opportunities: Studies of HIV adaptation are revealing why certain HLA alleles are protective in some populations and not others; they are identifying immunogenic viral epitopes that harbor high mutational barriers to escape, and they may help illuminate novel, vaccine-relevant HIV epitopes in regions where circulating adaptation is extensive. Elucidation of HLA-driven adapted and nonadapted viral forms in different human populations and HIV subtypes also renders 'personalized' immunogen selection, as a component of HIV cure strategies, conceptually feasible.

SUMMARY

Though adaptation represents a major challenge to HIV prevention and cure, achieving an in-depth understanding of this phenomenon can help move the design of such strategies forward.

Structural topology defines protective CD8+ T cell epitopes in the HIV proteome

Mutationally constrained epitopes of variable pathogens represent promising targets for vaccine design but are not reliably identified by sequence conservation. In this study, we employed structure-based network analysis, which applies network theory to HIV protein structure data to quantitate the topological importance of individual amino acid residues. Mutation of residues at important network positions disproportionately impaired viral replication and occurred with high frequency in epitopes presented by protective human leukocyte antigen (HLA) class I alleles. Moreover, CD8⁺ T cell targeting of highly networked epitopes distinguished individuals who naturally control HIV, even in the absence of protective HLA alleles. This approach thereby provides a mechanistic basis for immune control and a means to identify CD8⁺ T cell epitopes of topological importance for rational immunogen design, including a T cell-based HIV vaccine.

Virus Adaptation and Selection Following Challenge of Animals Vaccinated against Classical Swine Fever Virus

Vaccines against classical swine fever have proven very effective in protecting pigs from this deadly disease. However, little is known about how vaccination impacts the selective pressures acting on the classical swine fever virus (CSFV). Here we use high-throughput sequencing of viral genomes to investigate evolutionary changes in virus populations following the challenge of naïve and vaccinated pigs with the highly virulent CSFV strain “Koslov”. The challenge inoculum contained an ensemble of closely related viral sequences, with three major haplotypes being present, termed A, B, and C. After the challenge, the viral haplotype A was preferentially located within the tonsils of naïve animals but was highly prevalent in the sera of all vaccinated animals. We find that the viral population structure in naïve pigs after infection is very similar to that in the original inoculum. In contrast, the viral population in vaccinated pigs, which only underwent transient low-level viremia, displayed several distinct changes including the emergence of 16 unique non-synonymous single nucleotide polymorphisms (SNPs) that were not detectable in the challenge inoculum. Further analysis showed a significant loss of heterogeneity and an increasing positive selection acting on the virus populations in the vaccinated pigs. We conclude that vaccination imposes a strong selective pressure on viruses that subsequently replicate within the vaccinated animal.

HLA-I Associated Adaptation Dampens CD8 T-Cell Responses in HIV Ad5-Vectored Vaccine Recipients

HLA-I-associated human immunodeficiency virus (HIV) adaptation is known to negatively affect disease progression and CD8 T-cell responses. We aimed to assess how HLA-I-associated adaptation affects HIV vaccine-induced CD8 T-cell responses in 2 past vaccine efficacy trials. We found that vaccine-encoded adapted epitopes were less immunogenic than vaccine-encoded nonadapted epitopes, and adapted epitope-specific responses were less polyfunctional than nonadapted epitope-specific responses. Along those lines, vaccine recipients with higher HLA-I adaptation to the Gag vaccine insert mounted less polyfunctional CD8 T-cell responses at the protein level. Breadth of response, which correlated with viral control in recipients who became infected, is also dampened by HLA-I adaptation. These findings suggest that HLA-I-associated adaptation is an important consideration for strategies aiming to induce robust CD8 T-cell responses.

Immunization of BLT Humanized Mice Redirects T Cell Responses to Gag and Reduces Acute HIV-1 Viremia

BLT (bone marrow-liver-thymus) humanized mice, which reconstitute a functional human immune system, develop prototypic human virus-specific CD8⁺ T cell responses following infection with human immunodeficiency virus type 1 (HIV-1). We explored the utility of the BLT model for HIV-1 vaccine development by immunizing BLT mice against the conserved viral Gag protein, utilizing a rapid prime-boost protocol of poly(lactic-co-glycolic) acid microparticles and a replication-defective herpes simplex virus (HSV) recombinant vector. After HIV-1 challenge, the mice developed broad, proteome-wide gamma interferon-positive (IFN-γ⁺) T cell responses against HIV-1 that reached magnitudes equivalent to what is observed in HIV-1-infected individuals. The functionality of these responses was underscored by the consistent emergence of escape mutations in multiple CD8⁺ T cell epitopes during the course of infection. Although prechallenge vaccine-induced responses were largely undetectable, the Gag immunization increased both the magnitude and the kinetics of anamnestic Gag-specific T cell responses following HIV-1 infection, and the magnitude of these postchallenge Gag-specific responses was inversely correlated with acute HIV-1 viremia. Indeed, Gag immunization was associated with a modest but significant 0.5-log reduction in HIV-1 viral load when analyzed across four experimental groups of BLT mice. Notably, the HSV vector induced elevated plasma concentrations of polarizing cytokines and chemotactic factors, including interleukin-12p70 (IL-12p70) and MIP-1α, which were positively correlated with the magnitude of Gag-specific responses. Overall, these results support the ability of BLT mice to recapitulate human pathogen-specific T cell responses and to respond to immunization; however, additional improvements to the model are required to develop a robust system for testing HIV-1 vaccine efficacy.IMPORTANCE Advances in the development of humanized mice have raised the possibility of a small-animal model for preclinical testing of an HIV-1 vaccine. Here, we describe the capacity of BLT humanized mice to mount broadly directed HIV-1-specific human T cell responses that are functionally active, as indicated by the rapid emergence of viral escape mutations. Although immunization of BLT mice with the conserved viral Gag protein did not result in detectable prechallenge responses, it did increase the magnitude and kinetics of postchallenge Gag-specific T cell responses, which was associated with a modest but significant reduction in acute HIV-1 viremia. Additionally, the BLT model revealed immunization-associated increases in the plasma concentrations of immunomodulatory cytokines and chemokines that correlated with more robust T cell responses. These data support the potential utility of the BLT humanized mouse for HIV-1 vaccine development but suggest that additional improvements to the model are warranted.

Effect of diversity in gp41 membrane proximal external region of primary HIV-1 Indian subtype C sequences on interaction with broadly neutralizing antibodies 4E10 and 10E8

Human Immunodeficiency Virus-1 Clade C (HIV-1C) dominates the AIDS epidemic in India, afflicting 2.1 million individuals within the country and more than 15 million people worldwide. Membrane proximal external region (MPER) is an attractive target for broadly neutralizing antibody (bNAb) based therapies. However, information on MPER sequence diversity from India is meagre due to limited sampling of primary viral sequences. In the present study, we examined the variation in MPER of HIV-1C from 24 individuals in Mumbai, India by high throughput sequencing of uncultured viral sequences. Deep sequencing of MPER (662-683; HXB2 envelope amino acid numbering) allowed quantification of intra-individual variation up to 65% at positions 662, 665, 668, 674 and 677 within this region. These variable positions included contact sites targeted by bNAbs 2F5, Z13e1, 4E10 as well as 10E8. Both major and minor epitope variants i.e. 'haplotypes' were generated for each sample dataset. A total of 23, 34 and 25 unique epitope haplotypes could be identified for bNAbs 2F5, Z13e1 and 4E10/10E8 respectively. Further analysis of 4E10 and 10E8 epitopes from our dataset and meta-analysis of previously reported HIV-1 sequences from India revealed 26 epitopes (7 India-specific), heretofore untested for neutralization sensitivity. Peptide-Ab docking predicted 13 of these to be non-binding to 10E8. ELISA, Surface Plasmon Resonance and peptide inhibition of HIV-1 neutralization assays were then performed which validated predicted weak/non-binding interactions for peptides corresponding to six of these epitopes. These results highlight the under-representation of 10E8 non-binding HIV-1C MPER sequences from India. Our study thus underscores the need for increased surveillance of primary circulating envelope sequences for development of efficacious bNAb-based interventions in India.

Clonal interference can cause wavelet-like oscillations of multilocus linkage disequilibrium

Within-host adaptation of pathogens such as human immunodeficiency virus (HIV) often occurs at more than two loci. Multiple beneficial mutations may arise simultaneously on different genetic backgrounds and interfere, affecting each other's fixation trajectories. Here, we explore how these evolutionary dynamics are mirrored in multilocus linkage disequilibrium (MLD), a measure of multi-way associations between alleles. In the parameter regime corresponding to HIV, we show that deterministic early infection models induce MLD to oscillate over time in a wavelet-like fashion. We find that the frequency of these oscillations is proportional to the rate of adaptation. This signature is robust to drift, but can be eroded by high variation in fitness effects of beneficial mutations. Our findings suggest that MLD oscillations could be used as a signature of interference among multiple equally advantageous mutations and may aid the interpretation of MLD in data.

There is no Difference Between Sequences of HIV-1 Infected Patients with Stable Clinical Status and HIV-1 Reference Sequence

BACKGROUND

The rate of human immunodeficiency virus type 1 (HIV-1) infection in Iran has increased dramatically in the past few years. HIV-1 genome sequences are pivotal for large-scale studies of inter- and intra-host evolution. To understand the molecular difference between reference HIV-1 isolate and two HIV-1 infected patients in Iran, we conducted this study to analyze some genome segments of Iranian HIV-1 isolates.

METHODS

Two HIV-1-infected individuals who were under antiretroviral therapy (ARV) for 8 years with stable clinical status were enrolled. The patient's plasma samples were used for the Gag-Pol genome sequences (4500 nt). The phylogenetic tree and similarity plotty were obtained based on Gag-Pol sequences.

RESULTS

Both HIV-1-infected isolates belonged to CRF35_AD subtype even though one of them had drug resistance. The HIV genome and protein sequences showed no clear difference between genome and protein sequences of our samples and the reference sequence.

CONCLUSIONS

Our patient's stable clinical status had no connection to genome sequence; which could be owing to immunological factors or other patient's mode which are still unknown.

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.

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

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

Role of HCV Viremia in Corroborated HCV Transmission Events Within Young Adult Injecting Partnerships

Background

Hepatitis C virus (HCV), a major cause of morbidity and mortality, is common and rising among young persons who inject drugs (PWID). Reducing the level of viremia may be an intervention, yet the impact of viremia on HCV transmission is unknown.

Methods

We conducted a prospective study of injecting partnerships (Partner Study) of young adult (age < 30 years) PWID within the UFO Study, which enrolled those at risk for HCV or with seronegative viremic infection and up to 3 HCV RNA-positive regular injecting partners. We examined the level of HCV viremia and stage of infection in the HCV-positive partner in regression analyses of HCV transmission events that were corroborated via HCV phylogenetic linkage analyses.

Results

We enrolled 69 at-risk/acutely infected PWID. There were 25 new HCV infections (incidence rate, 35.9 per 100 person-years; 95% confidence interval [CI], 24.3-53.2 per 100 person-years); 12/25 (48%) were phylogenetically linked to at least 1 partner. We found no association between the infected partner's quantitative level of HCV viremia and likely transmission in multivariate analyses (adjusted odds ratio [AOR], 0.90; 95% confidence interval [CI], 0.55-1.46); however, seronegative viremic infection in the infected partner was associated with increased transmission (AOR, 28.02; 95% CI, 5.61-139.95).

Conclusions

The HCV viremia level was not associated with increased odds of transmission, yet acute HCV infection (seronegative viremic) was. Explanations include high-risk behavior during acute infection or missed fluctuations in viremia during acute infection. Both point to the need for frequent testing to detect new infection and attempt to prevent onward transmission.

PCR Amplification Strategies Towards Full-length HIV-1 Genome Sequencing

The advent of next-generation sequencing has enabled greater resolution of viral diversity and improved feasibility of full viral genome sequencing allowing routine HIV-1 full genome sequencing in both research and diagnostic settings. Regardless of the sequencing platform selected, successful PCR amplification of the HIV-1 genome is essential for sequencing template preparation. As such, full HIV-1 genome amplification is a crucial step in dictating the successful and reliable sequencing downstream. Here we reviewed existing PCR protocols leading to HIV-1 full genome sequencing. In addition to the discussion on basic considerations on relevant PCR design, the advantages as well as the pitfalls of the published protocols were reviewed.

Development of a Versatile, Near Full Genome Amplification and Sequencing Approach for a Broad Variety of HIV-1 Group M Variants

Near full genome sequencing (NFGS) of HIV-1 is required to assess the genetic composition of HIV-1 strains comprehensively. Population-wide, it enables a determination of the heterogeneity of HIV-1 and the emergence of novel/recombinant strains, while for each individual it constitutes a diagnostic instrument to assist targeted therapeutic measures against viral components. There is still a lack of robust and adaptable techniques for efficient NFGS from miscellaneous HIV-1 subtypes. Using rational primer design, a broad primer set was developed for the amplification and sequencing of diverse HIV-1 group M variants from plasma. Using pure subtypes as well as diverse, unique recombinant forms (URF), variable amplicon approaches were developed for NFGS comprising all functional genes. Twenty-three different genomes composed of subtypes A (A1), B, F (F2), G, CRF01_AE, CRF02_AG, and CRF22_01A1 were successfully determined. The NFGS approach was robust irrespective of viral loads (≥306 copies/mL) and amplification method. Third-generation sequencing (TGS), single genome amplification (SGA), cloning, and bulk sequencing yielded similar outcomes concerning subtype composition and recombinant breakpoint patterns. The introduction of a simple and versatile near full genome amplification, sequencing, and cloning method enables broad application in phylogenetic studies of diverse HIV-1 subtypes and can contribute to personalized HIV therapy and diagnosis.

HIV Subtype and Nef-Mediated Immune Evasion Function Correlate with Viral Reservoir Size in Early-Treated Individuals

The HIV accessory protein Nef modulates key immune evasion and pathogenic functions, and its encoding gene region exhibits high sequence diversity. Given the recent identification of early HIV-specific adaptive immune responses as novel correlates of HIV reservoir size, we hypothesized that viral factors that facilitate the evasion of such responses-namely, Nef genetic and functional diversity-might also influence reservoir establishment and/or persistence. We isolated baseline plasma HIV RNA-derived nef clones from 30 acute/early-infected individuals who participated in a clinical trial of early combination antiretroviral therapy (cART) (<6 months following infection) and assessed each Nef clone's ability to downregulate CD4 and human leukocyte antigen (HLA) class I in vitro We then explored the relationships between baseline clinical, immunological, and virological characteristics and the HIV reservoir size measured 48 weeks following initiation of suppressive cART (where the reservoir size was quantified in terms of the proviral DNA loads as well as the levels of replication-competent HIV in CD4⁺ T cells). Maximal within-host Nef-mediated downregulation of HLA, but not CD4, correlated positively with post-cART proviral DNA levels (Spearman's R = 0.61, P = 0.0004) and replication-competent reservoir sizes (Spearman's R = 0.36, P = 0.056) in univariable analyses. Furthermore, the Nef-mediated HLA downregulation function was retained in final multivariable models adjusting for established clinical and immunological correlates of reservoir size. Finally, HIV subtype B-infected persons (n = 25) harbored significantly larger viral reservoirs than non-subtype B-infected persons (2 infected with subtype CRF01_AE and 3 infected with subtype G). Our results highlight a potentially important role of viral factors-in particular, HIV subtype and accessory protein function-in modulating viral reservoir establishment and persistence.IMPORTANCE While combination antiretroviral therapies (cART) have transformed HIV infection into a chronic manageable condition, they do not act upon the latent HIV reservoir and are therefore not curative. As HIV cure or remission should be more readily achievable in individuals with smaller HIV reservoirs, achieving a deeper understanding of the clinical, immunological, and virological determinants of reservoir size is critical to eradication efforts. We performed a post hoc analysis of 30 participants of a clinical trial of early cART who had previously been assessed in detail for their clinical, immunological, and reservoir size characteristics. We observed that the HIV subtype and autologous Nef-mediated HLA downregulation function correlated with the viral reservoir size measured approximately 1 year post-cART initiation. Our findings highlight virological characteristics-both genetic and functional-as possible novel determinants of HIV reservoir establishment and persistence.

Laboratory Methods in Molecular Epidemiology: Viral Infections

Viruses, which are the most abundant biological entities on the planet, have been regarded as the "dark matter" of biology in the sense that despite their ubiquity and frequent presence in large numbers, their detection and analysis are not always straightforward. The majority of them are very small (falling under the limit of 0.5 μm), and collectively, they are extraordinarily diverse. In fact, the majority of the genetic diversity on the planet is found in the so-called virosphere, or the world of viruses. Furthermore, the most frequent viral agents of disease in humans display an RNA genome, and frequently evolve very fast, due to the fact that most of their polymerases are devoid of proofreading activity. Therefore, their detection, genetic characterization, and epidemiological surveillance are rather challenging. This review (part of the Curated Collection on Advances in Molecular Epidemiology of Infectious Diseases) describes many of the methods that, throughout the last few decades, have been used for viral detection and analysis. Despite the challenge of having to deal with high genetic diversity, the majority of these methods still depend on the amplification of viral genomic sequences, using sequence-specific or sequence-independent approaches, exploring thermal profiles or a single nucleic acid amplification temperature. Furthermore, viral populations, and especially those with RNA genomes, are not usually genetically uniform but encompass swarms of genetically related, though distinct, viral genomes known as viral quasispecies. Therefore, sequence analysis of viral amplicons needs to take this fact into consideration, as it constitutes a potential analytic problem. Possible technical approaches to deal with it are also described here. *This article is part of a curated collection.