Founder effects in the assessment of HIV polymorphisms and HLA allele associations

Accurate Detection of Convergent Mutations in Large Protein Alignments With ConDor

Evolutionary convergences are observed at all levels, from phenotype to DNA and protein sequences, and changes at these different levels tend to be correlated. Notably, convergent mutations can lead to convergent changes in phenotype, such as changes in metabolism, drug resistance, and other adaptations to changing environments. We propose a two-component approach to detect mutations subject to convergent evolution in protein alignments. The "Emergence" component selects mutations that emerge more often than expected, while the "Correlation" component selects mutations that correlate with the convergent phenotype under study. With regard to Emergence, a phylogeny deduced from the alignment is provided by the user and is used to simulate the evolution of each alignment position. These simulations allow us to estimate the expected number of mutations in a neutral model, which is compared to the observed number of mutations in the data studied. In Correlation, a comparative phylogenetic approach, is used to measure whether the presence of each of the observed mutations is correlated with the convergent phenotype. Each component can be used on its own, for example Emergence when no phenotype is available. Our method is implemented in a standalone workflow and a webserver, called ConDor. We evaluate the properties of ConDor using simulated data, and we apply it to three real datasets: sedge PEPC proteins, HIV reverse transcriptase, and fish rhodopsin. The results show that the two components of ConDor complement each other, with an overall accuracy that compares favorably to other available tools, especially on large datasets.

Quantifying how single dose Ad26.COV2.S vaccine efficacy depends on Spike sequence features

In the ENSEMBLE randomized, placebo-controlled phase 3 trial (NCT04505722), estimated single-dose Ad26.COV2.S vaccine efficacy (VE) was 56% against moderate to severe-critical COVID-19. SARS-CoV-2 Spike sequences were determined from 484 vaccine and 1,067 placebo recipients who acquired COVID-19. In this set of prespecified analyses, we show that in Latin America, VE was significantly lower against Lambda vs. Reference and against Lambda vs. non-Lambda [family-wise error rate (FWER) p < 0.05]. VE differed by residue match vs. mismatch to the vaccine-insert at 16 amino acid positions (4 FWER p < 0.05; 12 q-value ≤ 0.20); significantly decreased with physicochemical-weighted Hamming distance to the vaccine-strain sequence for Spike, receptor-binding domain, N-terminal domain, and S1 (FWER p < 0.001); differed (FWER ≤ 0.05) by distance to the vaccine strain measured by 9 antibody-epitope escape scores and 4 NTD neutralization-impacting features; and decreased (p = 0.011) with neutralization resistance level to vaccinee sera. VE against severe-critical COVID-19 was stable across most sequence features but lower against the most distant viruses.

Prevention efficacy of the broadly neutralizing antibody VRC01 depends on HIV-1 envelope sequence features

In the Antibody Mediated Prevention (AMP) trials (HVTN 704/HPTN 085 and HVTN 703/HPTN 081), prevention efficacy (PE) of the monoclonal broadly neutralizing antibody (bnAb) VRC01 (vs. placebo) against HIV-1 acquisition diagnosis varied according to the HIV-1 Envelope (Env) neutralization sensitivity to VRC01, as measured by 80% inhibitory concentration (IC80). Here, we performed a genotypic sieve analysis, a complementary approach to gaining insight into correlates of protection that assesses how PE varies with HIV-1 sequence features. We analyzed HIV-1 Env amino acid (AA) sequences from the earliest available HIV-1 RNA-positive plasma samples from AMP participants diagnosed with HIV-1 and identified Env sequence features that associated with PE. The strongest Env AA sequence correlate in both trials was VRC01 epitope distance that quantifies the divergence of the VRC01 epitope in an acquired HIV-1 isolate from the VRC01 epitope of reference HIV-1 strains that were most sensitive to VRC01-mediated neutralization. In HVTN 704/HPTN 085, the Env sequence-based predicted probability that VRC01 IC80 against the acquired isolate exceeded 1 µg/mL also significantly associated with PE. In HVTN 703/HPTN 081, a physicochemical-weighted Hamming distance across 50 VRC01 binding-associated Env AA positions of the acquired isolate from the most VRC01-sensitive HIV-1 strain significantly associated with PE. These results suggest that incorporating mutation scoring by BLOSUM62 and weighting by the strength of interactions at AA positions in the epitope:VRC01 interface can optimize performance of an Env sequence-based biomarker of VRC01 prevention efficacy. Future work could determine whether these results extend to other bnAbs and bnAb combinations.

Insights to HIV-1 coreceptor usage by estimating HLA adaptation with Bayesian generalized linear mixed models

The mechanisms triggering the human immunodeficiency virus type I (HIV-1) to switch the coreceptor usage from CCR5 to CXCR4 during the course of infection are not entirely understood. While low CD4+ T cell counts are associated with CXCR4 usage, a predominance of CXCR4 usage with still high CD4+ T cell counts remains puzzling. Here, we explore the hypothesis that viral adaptation to the human leukocyte antigen (HLA) complex, especially to the HLA class II alleles, contributes to the coreceptor switch. To this end, we sequence the viral gag and env protein with corresponding HLA class I and II alleles of a new cohort of 312 treatment-naive, subtype C, chronically-infected HIV-1 patients from South Africa. To estimate HLA adaptation, we develop a novel computational approach using Bayesian generalized linear mixed models (GLMMs). Our model allows to consider the entire HLA repertoire without restricting the model to pre-learned HLA-polymorphisms. In addition, we correct for phylogenetic relatedness of the viruses within the model itself to account for founder effects. Using our model, we observe that CXCR4-using variants are more adapted than CCR5-using variants (p-value = 1.34e-2). Additionally, adapted CCR5-using variants have a significantly lower predicted false positive rate (FPR) by the geno2pheno[coreceptor] tool compared to the non-adapted CCR5-using variants (p-value = 2.21e-2), where a low FPR is associated with CXCR4 usage. Consequently, estimating HLA adaptation can be an asset in predicting not only coreceptor usage, but also an approaching coreceptor switch in CCR5-using variants. We propose the usage of Bayesian GLMMs for modeling virus-host adaptation in general.

The applications of machine learning in HIV neutralizing antibodies research-A systematic review

Machine learning algorithms play an essential role in bioinformatics and allow exploring the vast and noisy biological data in unrivaled ways. This paper is a systematic review of the applications of machine learning in the study of HIV neutralizing antibodies. This significant and vast research domain can pave the way to novel treatments and to a vaccine. We selected the relevant papers by investigating the available literature from the Web of Science and PubMed databases in the last decade. The computational methods are applied in neutralization potency prediction, neutralization span prediction against multiple viral strains, antibody-virus binding sites detection, enhanced antibodies design, and the study of the antibody-induced immune response. These methods are viewed from multiple angles spanning data processing, model description, feature selection, evaluation, and sometimes paper comparisons. The algorithms are diverse and include supervised, unsupervised, and generative types. Both classical machine learning and modern deep learning were taken into account. The review ends with our ideas regarding future research directions and challenges.

SARS-CoV-2 intra-host single-nucleotide variants associated with disease severity

Variants of severe acute respiratory syndrome coronavirus 2 frequently arise within infected individuals. Here, we explored the level and pattern of intra-host viral diversity in association with disease severity. Then, we analyzed information underlying these nucleotide changes to infer the impetus including mutational signatures and immune selection from neutralizing antibody or T-cell recognition. From 23 January to 31 March 2020, a set of cross-sectional samples were collected from individuals with homogeneous founder virus regardless of disease severity. Intra-host single-nucleotide variants (iSNVs) were enumerated using deep sequencing. Human leukocyte antigen (HLA) alleles were genotyped by Sanger sequencing. Medical records were collected and reviewed by attending physicians. A total of 836 iSNVs (3-106 per sample) were identified and distributed in a highly individualized pattern. The number of iSNVs paced with infection duration peaked within days and declined thereafter. These iSNVs did not stochastically arise due to a strong bias toward C > U/G > A and U > C/A > G substitutions in reciprocal proportion with escalating disease severity. Eight nonsynonymous iSNVs in the receptor-binding domain could escape from neutralization, and eighteen iSNVs were significantly associated with specific HLA alleles. The level and pattern of iSNVs reflect the in vivo viral-host interaction and the disease pathogenesis.

Evolution of the North American Lineage H7 Avian Influenza Viruses in Association with H7 Virus's Introduction to Poultry

The incursions of H7 subtype low-pathogenicity avian influenza virus (LPAIV) from wild birds into poultry and its mutations to highly pathogenic avian influenza virus (HPAIV) have been an ongoing concern in North America. Since 2000, 10 phylogenetically distinct H7 virus outbreaks from wild birds have been detected in poultry, six of which mutated to HPAIV. To study the molecular evolution of the H7 viruses that occurs when changing hosts from wild birds to poultry, we performed analyses of the North American H7 hemagglutinin (HA) genes to identify amino acid changes as the virus circulated in wild birds from 2000 to 2019. Then, we analyzed recurring HA amino acid changes and gene constellations of the viruses that spread from wild birds to poultry. We found six HA amino acid changes occurring during wild bird circulation and 10 recurring changes after the spread to poultry. Eight of the changes were in and around the HA antigenic sites, three of which were supported by positive selection. Viruses from each H7 outbreak had a unique genotype, with no specific genetic group associated with poultry outbreaks or mutation to HPAIV. However, the genotypes of the H7 viruses in poultry outbreaks tended to contain minor genetic groups less observed in wild bird H7 viruses, suggesting either a biased sampling of wild bird AIVs or a tendency of having reassortment with minor genetic groups prior to the virus's introduction to poultry. IMPORTANCE Wild bird-origin H7 subtype avian influenza viruses are a constant threat to commercial poultry, both directly by the disease they cause and indirectly through trade restrictions that can be imposed when the virus is detected in poultry. It is important to understand the genetic basis of why the North American lineage H7 viruses have repeatedly crossed the species barrier from wild birds to poultry. We examined the amino acid changes in the H7 viruses associated with poultry outbreaks and tried to determine gene reassortment related to poultry adaptation and mutations to HPAIV. The findings in this study increase the understanding of the evolutionary pathways of wild bird AIV before infecting poultry and the HA changes associated with adaptation of the virus in poultry.

Functional Profile of CD8+ T-Cells in Response to HLA-A*02:01-Restricted Mutated Epitopes Derived from the Gag Protein of Circulating HIV-1 Strains from Medellín, Colombia

CD8⁺ T-cells play a crucial role in the control of HIV replication. HIV-specific CD8⁺ T-cell responses rapidly expand since the acute phase of the infection, and it has been observed that HIV controllers harbor CD8⁺ T-cells with potent anti-HIV capacity. The development of CD8⁺ T-cell-based vaccine against HIV-1 has focused on searching for immunodominant epitopes. However, the strong immune pressure of CD8⁺ T-cells causes the selection of viral variants with mutations in immunodominant epitopes. Since HIV-1 mutations are selected under the context of a specific HLA-I, the circulation of viral variants with these mutations is highly predictable based on the most prevalent HLA-I within a population. We previously demonstrated the adaptation of circulating strains of HIV-1 to the HLA-A*02 molecule by identifying mutations under positive selection located in GC9 and SL9 epitopes derived from the Gag protein. Also, we used an in silico prediction approach and evaluated whether the mutations found had a higher or lower affinity to the HLA-A*02. Although this strategy allowed predicting the interaction between mutated peptides and HLA-I, the functional response of CD8⁺ T-cells that these peptides induce is unknown. In the present work, peripheral blood mononuclear cells from 12 HIV-1⁺ HLA-A*02:01⁺ individuals were stimulated with the mutated and wild-type peptides derived from the GC9 and SL9 epitopes. The functional profile of CD8⁺ T-cells was evaluated using flow cytometry, and the frequency of subpopulations was determined according to their number of functions and the polyfunctionality index. The results suggest that the quality of the response (polyfunctionality) could be associated with the binding affinity of the peptide to the HLA molecule, and the functional profile of specific CD8⁺ T-cells to mutated epitopes in individuals under cART is maintained.

HAMdetector: A Bayesian regression model that integrates information to detect HLA-associated mutations

MOTIVATION

A key process in anti-viral adaptive immunity is that the Human Leukocyte Antigen system (HLA) presents epitopes as Major Histocompatibility Complex I (MHC I) protein-peptide complexes on cell surfaces and in this way alerts CD8+ cytotoxic T-Lymphocytes (CTLs). This pathway exerts strong selection pressure on viruses, favoring viral mutants that escape recognition by the HLA/CTL system. Naturally, such immune escape mutations often emerge in highly variable viruses, e.g. HIV or HBV, as HLA-associated mutations (HAMs), specific to the hosts MHC I proteins. The reliable identification of HAMs is not only important for understanding viral genomes and their evolution, but it also impacts the development of broadly effective anti-viral treatments and vaccines against variable viruses. By their very nature, HAMs are amenable to detection by statistical methods in paired sequence/HLA data. However, HLA alleles are very polymorphic in the human host population which makes the available data relatively sparse and noisy. Under these circumstances, one way to optimize HAM detection is to integrate all relevant information in a coherent model. Bayesian inference offers a principled approach to achieve this.

RESULTS

We present a new Bayesian regression model for the detection of HAMs that integrates a sparsity-inducing prior, epitope predictions, and phylogenetic bias assessment, and that yields easily interpretable quantitative information on HAM candidates. The model predicts experimentally confirmed HAMs as having high posterior probabilities, and it performs well in comparison to state-of-the-art models for several data sets from individuals infected with HBV, HDV, and HIV.

AVAILABILITY

The source code of this software is available at https://github.com/HAMdetector/Escape.jl under a permissive MIT license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genomic Characterization of SARS-CoV-2 Isolated from Patients with Distinct Disease Outcomes in Mexico

The coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), has shown a wide spectrum of clinical manifestations ranging from asymptomatic infections to severe disease and death. Pre-existing medical conditions and age have been mainly linked to the development of severe disease; however, the potential association of viral genetic characteristics with different clinical conditions remains unclear. SARS-CoV-2 variants with increased transmissibility were detected early in the pandemics, and several variants with potential relevance for public health are currently circulating around the world. In this study, we characterized 57 complete SARS-CoV-2 genomes during the exponential growth phase of the early epidemiological curve in Mexico, in April 2020. Patients were categorized under distinct disease severity outcomes: mild disease or ambulatory care, severe disease or hospitalized, and deceased. To reduce bias related to risk factors, the patients were less than 60 years old and with no diagnosed comorbidities A trait-association phylogenomic approach was used to explore genotype–phenotype associations, represented by the co-occurrence of mutations, disease severity outcome categories, and clusters of Mexican sequences. Phylogenetic results revealed a higher genomic diversity compared to the initial viruses detected during the early stage of the local epidemic. We identified a total of 90 single nucleotide variants compared to the Wuhan-Hu-1 genome, including 54 nonsynonymous mutations. We did not find evidence for the co-occurrence of mutations associated with specific disease outcomes. Therefore, in the group of patients studied, disease severity was likely mainly driven by the host genetic background and other demographic factors.

IMPORTANCE The genetic association of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) with different clinical conditions remains unclear and needs further investigation. In this study, we characterized 57 complete SARS-CoV-2 genomes from patients in Mexico with distinct disease severity outcomes: mild disease or ambulatory care, severe disease or hospitalized, and deceased. To reduce bias related to risk factors the patients were less than 60 years old and with no diagnosed comorbidities. We did not find evidence for the co-occurrence of mutations associated with specific disease outcomes. Therefore, in the group of patients studied, disease severity was likely mainly driven by the host genetic background and other demographic factors.

A comprehensive analysis of genetic diversity of EBV reveals potential high-risk subtypes associated with nasopharyngeal carcinoma in China

Epstein-Barr virus (EBV), a widespread oncovirus, is associated with multiple cancers including nasopharyngeal carcinoma (NPC), gastric cancer and diverse lymphoid malignancies. Recent studies reveal that specific EBV strains or subtypes are associated with NPC development in endemic regions. However, these NPC specific subtypes were only identified in a portion of infected individuals due possibly to the limited samples size studied or the complicated population structures of the virus. To identify additional high-risk EBV subtypes, we conducted a comprehensive genetic analysis of 22 critical viral proteins by using the largest dataset of 628 EBV genomes and 792 sequences of single target genes/proteins from GenBank. The phylogenetic, principal component and genetic structure analyses of these viral proteins were performed through worldwide populations. In addition to the general Asia-Western/Africa geographic segregation, population structure analysis showed a 'Chinese-unique' cluster (96.57% isolates from China) was highly enriched in the NPC patients, compared to the healthy individuals (89.6% vs. 44.5%, P < 0.001). The newly identified EBV subtypes, which contains four Chinese-specific NPC-associated amino acid substitutions (BALF2 V317M, BNRF1 G696R, V1222I and RPMS1 D51E), showed a robust positive association with the risk of NPC in China (Odds Ratio = 4.80, 20.00, 18.24 and 32.00 for 1, 2, 3 and 4 substitutions, respectively, P _trend <0.001). Interestingly, the coincidence of positively selected sites with NPC-associated substitutions suggests that adaptive nonsynonymous mutation on critical proteins, such as BNRF1, may interact with host immune system and contribute to the carcinogenesis of NPC. Our findings provide a comprehensive overview of EBV genetic structure for worldwide populations and offer novel clues to EBV carcinogenesis from the aspect of evolution.

Mutations that confer resistance to broadly-neutralizing antibodies define HIV-1 variants of transmitting mothers from that of non-transmitting mothers

Despite considerable reduction of mother-to-child transmission (MTCT) of HIV through use of maternal and infant antiretroviral therapy (ART), over 150,000 infants continue to become infected with HIV annually, falling far short of the World Health Organization goal of reaching <20,000 annual pediatric HIV cases worldwide by 2020. Prior to the widespread use of ART in the setting of pregnancy, over half of infants born to HIV-infected mothers were protected against HIV acquisition. Yet, the role of maternal immune factors in this protection against vertical transmission is still unclear, hampering the development of synergistic strategies to further reduce MTCT. It has been established that infant transmitted/founder (T/F) viruses are often resistant to maternal plasma, yet it is unknown if the neutralization resistance profile of circulating viruses predicts the maternal risk of transmission to her infant. In this study, we amplified HIV-1 envelope genes (env) by single genome amplification and produced representative Env variants from plasma of 19 non-transmitting mothers from the U.S. Women Infant Transmission Study (WITS), enrolled in the pre-ART era. Maternal HIV Env variants from non-transmitting mothers had similar sensitivity to autologous plasma as observed for non-transmitting variants from transmitting mothers. In contrast, infant variants were on average 30% less sensitive to paired plasma neutralization compared to non-transmitted maternal variants from both transmitting and non-transmitting mothers (p = 0.015). Importantly, a signature sequence analysis revealed that motifs enriched in env sequences from transmitting mothers were associated with broadly neutralizing antibody (bnAb) resistance. Altogether, our findings suggest that circulating maternal virus resistance to bnAb-mediated neutralization, but not autologous plasma neutralization, near the time of delivery, predicts increased MTCT risk. These results caution that enhancement of maternal plasma neutralization through passive or active vaccination during pregnancy may potentially drive the evolution of variants fit for vertical transmission.

Despite widespread, effective use of ART among HIV infected pregnant women, new pediatric HIV infections increase by about 150,000 every year. Thus, alternative strategies will be required to reduce MTCT and eliminate pediatric HIV infections. Interestingly, in the absence of ART, less than half of HIV-infected pregnant women will transmit HIV, suggesting natural immune protection of infants from virus acquisition. To understand the impact of maternal plasma autologous virus neutralization responses on MTCT, we compared the plasma and bnAb neutralization sensitivity of the circulating viral population present at the time of delivery in untreated, HIV-infected transmitting and non-transmitting mothers. While there was no significant difference in the ability of transmitting and non-transmitting women to neutralize their own circulating virus strains, specific genetic motifs enriched in variants from transmitting mothers were associated with resistance to bnAbs, suggesting that acquired bnAb resistance is a common feature of vertically-transmitted variants. This work suggests that enhancement of plasma neutralization responses in HIV-infected mothers through passive or active vaccination could further drive selection of variants that could be vertically transmitted, and cautions the use of passive bnAbs for HIV-1 prophylaxis or therapy during pregnancy.

Mucosal IL-4R antagonist HIV vaccination with SOSIP-gp140 booster can induce high-quality cytotoxic CD4+/CD8+ T cells and humoral responses in macaques

Inducing humoral, cellular and mucosal immunity is likely to improve the effectiveness of HIV-1 vaccine strategies. Here, we tested a vaccine regimen in pigtail macaques using an intranasal (i.n.) recombinant Fowl Pox Virus (FPV)-gag pol env-IL-4R antagonist prime, intramuscular (i.m.) recombinant Modified Vaccinia Ankara Virus (MVA)-gag pol-IL-4R antagonist boost followed by an i.m SOSIP-gp140 boost. The viral vector-expressed IL-4R antagonist transiently inhibited IL-4/IL-13 signalling at the vaccination site. The SOSIP booster not only induced gp140-specific IgG, ADCC (antibody-dependent cellular cytotoxicity) and some neutralisation activity, but also bolstered the HIV-specific cellular and humoral responses. Specifically, superior sustained systemic and mucosal HIV Gag-specific poly-functional/cytotoxic CD4⁺ and CD8⁺ T cells were detected with the IL-4R antagonist adjuvanted strategy compared to the unadjuvanted control. In the systemic compartment elevated Granzyme K expression was linked to CD4⁺ T cells, whilst Granzyme B/TIA-1 to CD8⁺ T cells. In contrast, the cytotoxic marker expression by mucosal CD4⁺ and CD8⁺ T cells differed according to the mucosal compartment. This vector-based mucosal IL-4R antagonist/SOSIP booster strategy, which promotes cytotoxic mucosal CD4⁺ T cells at the first line of defence, and cytotoxic CD4⁺ and CD8⁺ T cells plus functional antibodies in the blood, may prove valuable in combating mucosal infection with HIV-1 and warrants further investigation.

HIV-1 p24Gag adaptation to modern and archaic HLA-allele frequency differences in ethnic groups contributes to viral subtype diversification

Pathogen-driven selection and past interbreeding with archaic human lineages have resulted in differences in human leukocyte antigen (HLA)-allele frequencies between modern human populations. Whether or not this variation affects pathogen subtype diversification is unknown. Here we show a strong positive correlation between ethnic diversity in African countries and both human immunodeficiency virus (HIV)-1 p24gag and subtype diversity. We demonstrate that ethnic HLA-allele differences between populations have influenced HIV-1 subtype diversification as the virus adapted to escape common antiviral immune responses. The evolution of HIV Subtype B (HIV-B), which does not appear to be indigenous to Africa, is strongly affected by immune responses associated with Eurasian HLA variants acquired through adaptive introgression from Neanderthals and Denisovans. Furthermore, we show that the increasing and disproportionate number of HIV-infections among African Americans in the USA drive HIV-B evolution towards an Africa-centric HIV-1 state. Similar adaptation of other pathogens to HLA variants common in affected populations is likely.

Parallel evolution in the emergence of highly pathogenic avian influenza A viruses

Parallel molecular evolution and adaptation are important phenomena commonly observed in viruses. Here, we exploit parallel molecular evolution to understand virulence evolution in avian influenza viruses (AIV). Highly-pathogenic AIVs evolve independently from low-pathogenic ancestors via acquisition of polybasic cleavage sites. Why some AIV lineages but not others evolve in this way is unknown. We hypothesise that the parallel emergence of highly-pathogenic AIV may be facilitated by permissive or compensatory mutations occurring across the viral genome. We combine phylogenetic, statistical and structural approaches to discover parallel mutations in AIV genomes associated with the highly-pathogenic phenotype. Parallel mutations were screened using a statistical test of mutation-phenotype association and further evaluated in the contexts of positive selection and protein structure. Our resulting mutational panel may help to reveal new links between virulence evolution and other traits, and raises the possibility of predicting aspects of AIV evolution.

SARS-CoV-2 ORF8 and SARS-CoV ORF8ab: Genomic Divergence and Functional Convergence

The COVID-19 pandemic, in the first seven months, has led to more than 15 million confirmed infected cases and 600,000 deaths. SARS-CoV-2, the causative agent for COVID-19, has proved to be a great challenge for its ability to spread in asymptomatic stages and the diverse disease spectrum it has generated. This has created a challenge of unimaginable magnitude, not only affecting human health and life but also potentially generating a long-lasting socioeconomic impact. Both medical sciences and biomedical research have also been challenged, consequently leading to a large number of clinical trials and vaccine initiatives. While known proteins of pathobiological importance are targets for these therapeutic approaches, it is imperative to explore other factors of viral significance. Accessory proteins are one such trait that have diverse roles in coronavirus pathobiology. Here, we analyze certain genomic characteristics of SARS-CoV-2 accessory protein ORF8 and predict its protein features. We have further reviewed current available literature regarding its function and comparatively evaluated these and other features of ORF8 and ORF8ab, its homolog from SARS-CoV. Because coronaviruses have been infecting humans repeatedly and might continue to do so, we therefore expect this study to aid in the development of holistic understanding of these proteins. Despite low nucleotide and protein identity and differentiating genome level characteristics, there appears to be significant structural integrity and functional proximity between these proteins pointing towards their high significance. There is further need for comprehensive genomics and structural-functional studies to lead towards definitive conclusions regarding their criticality and that can eventually define their relevance to therapeutics development.

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.

Evolutionary Analyses of Base-Pairing Interactions in DNA and RNA Secondary Structures

Pairs of nucleotides within functional nucleic acid secondary structures often display evidence of coevolution that is consistent with the maintenance of base-pairing. Here, we introduce a sequence evolution model, MESSI (Modeling the Evolution of Secondary Structure Interactions), that infers coevolution associated with base-paired sites in DNA or RNA sequence alignments. MESSI can estimate coevolution while accounting for an unknown secondary structure. MESSI can also use graphics processing unit parallelism to increase computational speed. We used MESSI to infer coevolution associated with GC, AU (AT in DNA), GU (GT in DNA) pairs in noncoding RNA alignments, and in single-stranded RNA and DNA virus alignments. Estimates of GU pair coevolution were found to be higher at base-paired sites in single-stranded RNA viruses and noncoding RNAs than estimates of GT pair coevolution in single-stranded DNA viruses. A potential biophysical explanation is that GT pairs do not stabilize DNA secondary structures to the same extent that GU pairs do in RNA. Additionally, MESSI estimates the degrees of coevolution at individual base-paired sites in an alignment. These estimates were computed for a SHAPE-MaP-determined HIV-1 NL4-3 RNA secondary structure. We found that estimates of coevolution were more strongly correlated with experimentally determined SHAPE-MaP pairing scores than three nonevolutionary measures of base-pairing covariation. To assist researchers in prioritizing substructures with potential functionality, MESSI automatically ranks substructures by degrees of coevolution at base-paired sites within them. Such a ranking was created for an HIV-1 subtype B alignment, revealing an excess of top-ranking substructures that have been previously identified as having structure-related functional importance, among several uncharacterized top-ranking substructures.

No evidence for distinct types in the evolution of SARS-CoV-2

A recent study by Tang et al. claimed that two major types of severe acute respiratory syndrome-coronavirus-2 (CoV-2) had evolved in the ongoing CoV disease-2019 pandemic and that one of these types was more 'aggressive' than the other. Given the repercussions of these claims and the intense media coverage of these types of articles, we have examined in detail the data presented by Tang et al., and show that the major conclusions of that paper cannot be substantiated. Using examples from other viral outbreaks, we discuss the difficulty in demonstrating the existence or nature of a functional effect of a viral mutation, and we advise against overinterpretation of genomic data during the pandemic.

Coevolution analysis of amino-acids reveals diversified drug-resistance solutions in viral sequences: a case study of hepatitis B virus

The study of mutational landscapes of viral proteins is fundamental for the understanding of the mechanisms of cross-resistance to drugs and the design of effective therapeutic strategies based on several drugs. Antiviral therapy with nucleos(t)ide analogues targeting the hepatitis B virus (HBV) polymerase protein (Pol) can inhibit disease progression by suppression of HBV replication and makes it an important case study. In HBV, treatment may fail due to the emergence of drug-resistant mutants. Primary and compensatory mutations have been associated with lamivudine resistance, whereas more complex mutational patterns are responsible for resistance to other HBV antiviral drugs. So far, all known drug-resistance mutations are located in one of the four Pol domains, called reverse transcriptase. We demonstrate that sequence covariation identifies drug-resistance mutations in viral sequences. A new algorithmic strategy, BIS2TreeAnalyzer, is designed to apply the coevolution analysis method BIS2, successfully used in the past on small sets of conserved sequences, to large sets of evolutionary related sequences. When applied to HBV, BIS2TreeAnalyzer highlights diversified viral solutions by discovering thirty-seven positions coevolving with residues known to be associated with drug resistance and located on the four Pol domains. These results suggest a sequential mechanism of emergence for some mutational patterns. They reveal complex combinations of positions involved in HBV drug resistance and contribute with new information to the landscape of HBV evolutionary solutions. The computational approach is general and can be applied to other viral sequences when compensatory mutations are presumed.

Exploring the interactions between the human and viral genomes

Over the last decade, genome-wide association studies led to major advances in identifying human genetic variants associated with infectious disease susceptibility. On the pathogen side, comparable methods are now applied to identify disease-modulating pathogen variants. As host and pathogen variants jointly determine disease outcomes, the most recent development has been to explore simultaneously host and pathogen genomes, through so-called genome-to-genome studies. In this review, we provide some background on the development of genome-to-genome analysis and we detail the first wave of studies in this emerging field, which focused on patients chronically infected with HIV and hepatitis C virus. We also discuss the need for novel statistical methods to better tackle the issues of population stratification and multiple testing. Finally, we speculate on future research areas where genome-to-genome analysis may prove to be particularly effective.

Aiming for protective T-cell responses: a focus on the first generation conserved-region HIVconsv vaccines in preventive and therapeutic clinical trials

Introduction: Despite life-saving antiretroviral drugs, an effective HIV-1 vaccine is the best solution and likely a necessary component of any strategy for halting the AIDS epidemic. The currently prevailing aim is to pursue antibody-mediated vaccine protection. With ample evidence for the ability of T cells to control HIV-1 replication, their protective potential should be also harnessed by vaccination. The challenge is to elicit not just any, but protective T cells.Areas covered: This article reviews the clinical experience with the first-generation conserved-region immunogen HIVconsv delivered by combinations of plasmid DNA, simian adenovirus, and poxvirus MVA. The aim of our strategy is to induce strong and broad T cells targeting functionally important parts of HIV-1 proteins common to global variants. These vaccines were tested in eight phase 1/2 preventive and therapeutic clinical trials in Europe and Africa, and induced high frequencies of broadly specific CD8⁺ T cells capable of in vitro inhibition of four major HIV-1 clades A, B, C and D, and in combination with latency-reactivating agent provided a signal of drug-free virological control in early treated patients.Expert opinion: A number of critical T-cell traits have to come together at the same time to achieve control over HIV-1.

CD8 T cells targeting adapted epitopes in chronic HIV infection promote dendritic cell maturation and CD4 T cell trans-infection

HIV-1 frequently escapes from CD8 T cell responses via HLA-I restricted adaptation, leading to the accumulation of adapted epitopes (AE). We previously demonstrated that AE compromise CD8 T cell responses during acute infection and are associated with poor clinical outcomes. Here, we examined the impact of AE on CD8 T cell responses and their biological relevance in chronic HIV infection (CHI). In contrast to acute infection, the majority of AE are immunogenic in CHI. Longitudinal analyses from acute to CHI showed an increased frequency and magnitude of AE-specific IFNγ responses compared to NAE-specific ones. These AE-specific CD8 T cells also were more cytotoxic to CD4 T cells. In addition, AE-specific CD8 T cells expressed lower levels of PD1 and CD57, as well as higher levels of CD28, suggesting a more activated and less exhausted phenotype. During CHI, viral sequencing identified AE-encoding strains as the dominant quasispecies. Despite increased CD4 T cell cytotoxicity, CD8 T cells responding to AE promoted dendritic cell (DC) maturation and CD4 T cell trans-infection perhaps explaining why AE are predominant in CHI. Taken together, our data suggests that the emergence of AE-specific CD8 T cell responses in CHI confers a selective advantage to the virus by promoting DC-mediated CD4 T cell trans-infection.

HIV-1 infection remains a critical public health threat across the world. Over the past two decades, CD8 T cells have been clearly shown to exert immune pressure on HIV and drive viral adaptation. Previously, our group reported that such HLA-I associated adaptations can predict clinical outcomes and are beneficial to HIV-1 as CD8 T cells are unable to recognize epitopes with adaptation in acute HIV infection. However, it is still unclear how HIV-1 adaptation impacts CD8 T cells during chronic HIV infection. In this study, we observed an enhancement of CD8 T cell responses targeting adapted epitopes in chronic infection. Although these responses were cytotoxic, they also exhibited a “helper” effect by promoting viral infection of CD4 T cells via interaction with dendritic cells. This phenomenon may contribute to the persistence of adapted viruses. In summary, these findings present a novel mechanism of CD8 T cell driven HIV-1 adaptation.

Mapping the drivers of within-host pathogen evolution using massive data sets

Differences among hosts, resulting from genetic variation in the immune system or heterogeneity in drug treatment, can impact within-host pathogen evolution. Genetic association studies can potentially identify such interactions. However, extensive and correlated genetic population structure in hosts and pathogens presents a substantial risk of confounding analyses. Moreover, the multiple testing burden of interaction scanning can potentially limit power. We present a Bayesian approach for detecting host influences on pathogen evolution that exploits vast existing data sets of pathogen diversity to improve power and control for stratification. The approach models key processes, including recombination and selection, and identifies regions of the pathogen genome affected by host factors. Our simulations and empirical analysis of drug-induced selection on the HIV-1 genome show that the method recovers known associations and has superior precision-recall characteristics compared to other approaches. We build a high-resolution map of HLA-induced selection in the HIV-1 genome, identifying novel epitope-allele combinations.

Difficult-to-neutralize global HIV-1 isolates are neutralized by antibodies targeting open envelope conformations

The HIV-1 envelope (Env) is the target for neutralizing antibodies and exists on the surface of virions in open or closed conformations. Difficult-to-neutralize viruses (tier 2) express Env in a closed conformation antigenic for broadly neutralizing antibodies (bnAbs) but not for third variable region (V3) antibodies. Here we show that select V3 macaque antibodies elicited by Env vaccination can neutralize 26% of otherwise tier 2 HIV-1 isolates in standardized virus panels. The V3 antibodies only bound to Env in its open conformation. Thus, Envs on tier 2 viruses sample a state where the V3 loop is not in its closed conformation position. Envelope second variable region length, glycosylation sites and V3 amino acids were signatures of neutralization sensitivity. This study determined that open conformations of Env with V3 exposed are present on a subset of otherwise neutralization-resistant virions, therefore neutralization of tier 2 HIV-1 does not always indicate bnAb induction.

Parallel evolution of HIV-1 in a long-term experiment

One of the most intriguing puzzles in biology is the degree to which evolution is repeatable. The repeatability of evolution, or parallel evolution, has been studied in a variety of model systems, but has rarely been investigated with clinically relevant viruses. To investigate parallel evolution of HIV-1, we passaged two replicate HIV-1 populations for almost 1 year in each of two human T-cell lines. For each of the four evolution lines, we determined the genetic composition of the viral population at nine time points by deep sequencing the entire genome. Mutations that were carried by the majority of the viral population accumulated continuously over 1 year in each evolution line. Many majority mutations appeared in more than one evolution line, that is, our experiments showed an extreme degree of parallel evolution. In one of the evolution lines, 62% of the majority mutations also occur in another line. The parallelism impairs our ability to reconstruct the evolutionary history by phylogenetic methods. We show that one can infer the correct phylogenetic topology by including minority mutations in our analysis. We also find that mutation diversity at the beginning of the experiment is predictive of the frequency of majority mutations at the end of the experiment.

HIV-1 Neutralizing Antibody Signatures and Application to Epitope-Targeted Vaccine Design

Eliciting HIV-1-specific broadly neutralizing antibodies (bNAbs) remains a challenge for vaccine development, and the potential of passively delivered bNAbs for prophylaxis and therapeutics is being explored. We used neutralization data from four large virus panels to comprehensively map viral signatures associated with bNAb sensitivity, including amino acids, hypervariable region characteristics, and clade effects across four different classes of bNAbs. The bNAb signatures defined for the variable loop 2 (V2) epitope region of HIV-1 Env were then employed to inform immunogen design in a proof-of-concept exploration of signature-based epitope targeted (SET) vaccines. V2 bNAb signature-guided mutations were introduced into Env 459C to create a trivalent vaccine, and immunization of guinea pigs with V2-SET vaccines resulted in increased breadth of NAb responses compared with Env 459C alone. These data demonstrate that bNAb signatures can be utilized to engineer HIV-1 Env vaccine immunogens capable of eliciting antibody responses with greater neutralization breadth.

•

HIV-1 bNAb sensitivity signatures from 4 large virus panels mapped across 4 Ab classes

•

Non-contact hypervariable region characteristics are critical for bNAb sensitivity

•

HIV-1 Env 459C used alone as a vaccine can elicit modest tier 2 NAbs in guinea pigs

•

V2 bNAb signature-guided modifications in 459C enhanced neutralization breadth

Parallel molecular evolution and adaptation in viruses

Parallel molecular evolution is the independent evolution of the same genotype or phenotype from distinct ancestors. The simple genomes and rapid evolution of many viruses mean they are useful model systems for studying parallel evolution by natural selection. Parallel adaptation occurs in the context of several viral behaviours, including cross-species transmission, drug resistance, and host immune escape, and its existence suggests that at least some aspects of virus evolution and emergence are repeatable and predictable. We introduce examples of virus parallel evolution and summarise key concepts. We outline the difficulties in detecting parallel adaptation using virus genomes, with a particular focus on phylogenetic and structural approaches, and we discuss future approaches that may improve our understanding of the phenomenon.

Differential escape of HCV from CD8+ T cell selection pressure between China and Germany depends on the presenting HLA class I molecule

Adaptation of hepatitis C virus (HCV) to CD8⁺ T cell selection pressure is well described; however, it is unclear if HCV differentially adapts in different populations. Here, we studied HLA class I-associated viral sequence polymorphisms in HCV 1b isolates in a Chinese population and compared viral substitution patterns between Chinese and German populations. We identified three HLA class I-restricted epitopes in HCV NS3 with statistical support for selection pressure and found evidence for differential escape pathways between isolates from China and Germany depending on the HLA class I molecule. The substitution patterns particularly differed in the epitope VTLTHPITK_1635-1643 , which was presented by HLA-A*03 as well as HLA-A*11, two alleles with highly different frequencies in the two populations. In Germany, a substitution in position seven of the epitope was the most frequent substitution in the presence of HLA-A*03, functionally associated with immune escape and nearly absent in Chinese isolates. In contrast, the most frequent substitution in China was located at position two of the epitope and became the predominant consensus residue. Moreover, substitutions in position one of the epitope were significantly enriched in HLA-A*11-positive individuals in China and associated with different patterns of CD8⁺ T cell reactivity. Our study confirms the differential escape pathways selected by HCV that depended on different HLA class I alleles in Chinese and German populations, indicating that HCV differentially adapts to distinct HLA class I alleles in these populations. This result has important implications for vaccine design against highly variable and globally distributed pathogens, which may require matching antigen sequences to geographic regions for T cell-based vaccine strategies.

Genotypic and Mechanistic Characterization of Subtype-Specific HIV Adaptation to Host Cellular Immunity

The extent to which viral genetic context influences HIV adaptation to human leukocyte antigen (HLA) class I-restricted immune pressures remains incompletely understood. The Ugandan HIV epidemic, where major pandemic group M subtypes A1 and D cocirculate in a single host population, provides an opportunity to investigate this question. We characterized plasma HIV RNA gag, pol, and nef sequences, along with host HLA genotypes, in 464 antiretroviral-naive individuals chronically infected with HIV subtype A1 or D. Using phylogenetically informed statistical approaches, we identified HLA-associated polymorphisms and formally compared their strengths of selection between viral subtypes. A substantial number (32%) of HLA-associated polymorphisms identified in subtype A1 and/or D had previously been reported in subtype B, C, and/or circulating recombinant form 01_AE (CRF01_AE), confirming the shared nature of many HLA-driven escape pathways regardless of viral genetic context. Nevertheless, 34% of the identified HLA-associated polymorphisms were significantly differentially selected between subtypes A1 and D. Experimental investigation of select examples of subtype-specific escape revealed distinct underlying mechanisms with important implications for vaccine design: whereas some were attributable to subtype-specific sequence variation that influenced epitope-HLA binding, others were attributable to differential mutational barriers to immune escape. Overall, our results confirm that HIV genetic context is a key modulator of viral adaptation to host cellular immunity and highlight the power of combined bioinformatic and mechanistic studies, paired with knowledge of epitope immunogenicity, to identify appropriate viral regions for inclusion in subtype-specific and universal HIV vaccine strategies.IMPORTANCE The identification of HIV polymorphisms reproducibly selected under pressure by specific HLA alleles and the elucidation of their impact on viral function can help identify immunogenic viral regions where immune escape incurs a fitness cost. However, our knowledge of HLA-driven escape pathways and their functional costs is largely limited to HIV subtype B and, to a lesser extent, subtype C. Our study represents the first characterization of HLA-driven adaptation pathways in HIV subtypes A1 and D, which dominate in East Africa, and the first statistically rigorous characterization of differential HLA-driven escape across viral subtypes. The results support a considerable impact of viral genetic context on HIV adaptation to host HLA, where HIV subtype-specific sequence variation influences both epitope-HLA binding and the fitness costs of escape. Integrated bioinformatic and mechanistic characterization of these and other instances of differential escape could aid rational cytotoxic T-lymphocyte-based vaccine immunogen selection for both subtype-specific and universal HIV vaccines.

Predicting Promiscuous T Cell Epitopes for Designing a Vaccine Against Streptococcus pyogenes

One of the most dangerous human pathogens with high prevalence worldwide is Streptococcus pyogenes, which has major impacts on global morbidity and mortality. A major challenge for S. pyogenes vaccine development is the detection of epitopes that confer protection from infection by multiple S. pyogenes types. Our aim was to identify the most conserved and immunogenic antigens of S. pyogenes, which can be a potential candidate for vaccine design in the future. Eight important surface proteins were analyzed. Using different prediction servers, strongest epitopes were selected. They had the ability to stimulate the humoral and cell-mediated immune system. Molecular docking was performed for measuring free-binding energy of selected epitopes. Seven epitopes from three surface proteins were selected as potential candidates for vaccine development. Conservation of selected epitopes among different Streptococcus types was checked. Further in vitro and in vivo tests are required to validate the suitability of the epitopes for vaccine design.

Protection against a mixed SHIV challenge by a broadly neutralizing antibody cocktail

HIV-1 sequence diversity presents a major challenge for the clinical development of broadly neutralizing antibodies (bNAbs) for both therapy and prevention. Sequence variation in critical bNAb epitopes has been observed in most HIV-1-infected individuals and can lead to viral escape after bNAb monotherapy in humans. We show that viral sequence diversity can limit both the therapeutic and prophylactic efficacy of bNAbs in rhesus monkeys. We first demonstrate that monotherapy with the V3 glycan-dependent antibody 10-1074, but not PGT121, results in rapid selection of preexisting viral variants containing N332/S334 escape mutations and loss of therapeutic efficacy in simian-HIV (SHIV)-SF162P3-infected rhesus monkeys. We then show that the V3 glycan-dependent antibody PGT121 alone and the V2 glycan-dependent antibody PGDM1400 alone both fail to protect against a mixed challenge with SHIV-SF162P3 and SHIV-325c. In contrast, the combination of both bNAbs provides 100% protection against this mixed SHIV challenge. These data reveal that single bNAbs efficiently select resistant viruses from a diverse challenge swarm to establish infection, demonstrating the importance of bNAb cocktails for HIV-1 prevention.

Staged induction of HIV-1 glycan-dependent broadly neutralizing antibodies

A preventive HIV-1 vaccine should induce HIV-1-specific broadly neutralizing antibodies (bnAbs). However, bnAbs generally require high levels of somatic hypermutation (SHM) to acquire breadth, and current vaccine strategies have not been successful in inducing bnAbs. Because bnAbs directed against a glycosylated site adjacent to the third variable loop (V3) of the HIV-1 envelope protein require limited SHM, the V3-glycan epitope is an attractive vaccine target. By studying the cooperation among multiple V3-glycan B cell lineages and their coevolution with autologous virus throughout 5 years of infection, we identify key events in the ontogeny of a V3-glycan bnAb. Two autologous neutralizing antibody lineages selected for virus escape mutations and consequently allowed initiation and affinity maturation of a V3-glycan bnAb lineage. The nucleotide substitution required to initiate the bnAb lineage occurred at a low-probability site for activation-induced cytidine deaminase activity. Cooperation of B cell lineages and an improbable mutation critical for bnAb activity defined the necessary events leading to breadth in this V3-glycan bnAb lineage. These findings may, in part, explain why initiation of V3-glycan bnAbs is rare, and suggest an immunization strategy for inducing similar V3-glycan bnAbs.

HIV-1 Vpu Mediates HLA-C Downregulation

Many pathogens evade cytotoxic T lymphocytes (CTLs) by downregulating HLA molecules on infected cells, but the loss of HLA can trigger NK cell-mediated lysis. HIV-1 is thought to subvert CTLs while preserving NK cell inhibition by Nef-mediated downregulation of HLA-A and -B but not HLA-C molecules. We find that HLA-C is downregulated by most primary HIV-1 clones, including transmitted founder viruses, in contrast to the laboratory-adapted NL4-3 virus. HLA-C reduction is mediated by viral Vpu and reduces the ability of HLA-C restricted CTLs to suppress viral replication in CD4+ cells in vitro. HLA-A/B are unaffected by Vpu, and primary HIV-1 clones vary in their ability to downregulate HLA-C, possibly in response to whether CTLs or NK cells dominate immune pressure through HLA-C. HIV-2 also suppresses HLA-C expression through distinct mechanisms, underscoring the immune pressure HLA-C exerts on HIV. This viral immune evasion casts new light on the roles of CTLs and NK cells in immune responses against HIV.

Dynamic Co-evolution and Interaction of Avian Leukosis Virus Genetic Variants and Host Immune Responses

Subgroup J avian leukosis virus (ALV-J), a typical retrovirus, is characterized of existence of a cloud of diverse variants and considerable genetic diversity. Previous studies describing the evolutionary dynamics of ALV-J genetic variants mainly focused on the early infection period or few randomly selected clones. Here, we inoculated 30 specific-pathogen-free chickens with the same founder ALV-J stock of known genetic background. Six (three antibody positive and three antibody negative) chickens were selected among 15 chickens with viremia. Viruses were serially isolated in 36 weeks and then sequenced using MiSeq high-throughput sequencing platform. This produced the largest ALV-J dataset to date, composed of more than three million clean reads. Our results showed that host humoral immunity could greatly enhance the genetic diversity of ALV-J genetic variants. In particular, selection pressures promoted a dynamic proportional changes in ALV-J genetic variants frequency. Cross-neutralization experiment showed that along with the change of the dominant variant, the antibody titers specific to infectious clones corresponding to the most dominant variants in weeks 12 and 28 have also changed significantly in sera collected in weeks 16 and 32. In contrast, no shift of dominant variant was observed in antibody-negative chickens. Moreover, we identified a novel hypervariable region in the gp85 gene. Our study reveals the interaction between ALV-J and the host, which could facilitate the development of vaccines and antiviral drugs.

Genome-to-genome analysis highlights the effect of the human innate and adaptive immune systems on the hepatitis C virus

Outcomes of hepatitis C virus (HCV) infection and treatment depend on viral and host genetic factors. We use human genome-wide genotyping arrays and new whole-genome HCV viral sequencing technologies to perform a systematic genome-to-genome study of 542 individuals chronically infected with HCV, predominately genotype 3. We show that both HLA alleles and interferon lambda innate immune system genes drive viral genome polymorphism, and that IFNL4 genotypes determine HCV viral load through a mechanism that is dependent on a specific polymorphism in the HCV polyprotein. We highlight the interplay between innate immune responses and the viral genome in HCV control.

Histidine 375 Modulates CD4 Binding in HIV-1 CRF01_AE Envelope Glycoproteins

The envelope glycoproteins (Envs) from human immunodeficiency virus type 1 (HIV-1) mediate viral entry. The binding of the HIV-1 gp120 glycoprotein to CD4 triggers conformational changes in gp120 that allow high-affinity binding to its coreceptors. In contrast to all other Envs from the same phylogenetic group, M, which possess a serine (S) at position 375, those from CRF01_AE strains possess a histidine (H) at this location. This residue is part of the Phe43 cavity, where residue 43 of CD4 (a phenylalanine) engages with gp120. Here we evaluated the functional consequences of replacing this residue in two CRF01_AE Envs (CM244 and 92TH023) by a serine. We observed that reversion of amino acid 375 to a serine (H375S) resulted in a loss of functionality of both CRF01_AE Envs as measured by a dramatic loss in infectivity and ability to mediate cell-to-cell fusion. While no effects on processing or trimer stability of these variants were observed, decreased functionality could be linked to a major defect in CD4 binding induced by the replacement of H375 by a serine. Importantly, mutations of residues 61 (layer 1), 105 and 108 (layer 2), and 474 to 476 (layer 3) of the CRF01_AE gp120 inner domain layers to the consensus residues present in group M restored CD4 binding and wild-type levels of infectivity and cell-to-cell fusion. These results suggest a functional coevolution between the Phe43 cavity and the gp120 inner domain layers. Altogether, our observations describe the functional importance of amino acid 375H in CRF01_AE envelopes.

IMPORTANCE

A highly conserved serine located at position 375 in group M is replaced by a histidine in CRF01_AE Envs. Here we show that H375 is required for efficient CRF01_AE Env binding to CD4. Moreover, this work suggests that specific residues of the gp120 inner domain layers have coevolved with H375 in order to maintain its ability to mediate viral entry.

Potent and broad HIV-neutralizing antibodies in memory B cells and plasma

Induction of broadly neutralizing antibodies (bnAbs) is a goal of HIV-1 vaccine development. Antibody 10E8, reactive with the distal portion of the membrane-proximal external region (MPER) of HIV-1 gp41, is broadly neutralizing. However, the ontogeny of distal MPER antibodies and the relationship of memory B cell to plasma bnAbs are poorly understood. HIV-1-specific memory B cell flow sorting and proteomic identification of anti-MPER plasma antibodies from an HIV-1-infected individual were used to isolate broadly neutralizing distal MPER bnAbs of the same B cell clonal lineage. Structural analysis demonstrated that antibodies from memory B cells and plasma recognized the envelope gp41 bnAb epitope in a distinct orientation compared with other distal MPER bnAbs. The unmutated common ancestor of this distal MPER bnAb was autoreactive, suggesting lineage immune tolerance control. Construction of chimeric antibodies of memory B cell and plasma antibodies yielded a bnAb that potently neutralized most HIV-1 strains.

Virus-Host Gene Interactions Define HIV-1 Disease Progression

In this chapter, we will review recent research on the virology of HIV-1 transmission and the impact of the transmitted virus genotype on subsequent disease progression. In most instances of HIV-1 sexual transmission, a single genetic variant, or a very limited number of variants from the diverse viral quasi-species present in the transmitting partner establishes systemic infection. Transmission involves both stochastic and selective processes, such that in general a minority variant in the donor is transmitted. While there is clear evidence for selection, the biological properties that mediate transmission remain incompletely defined. Nevertheless, the genotype of the transmitted founder virus, which reflects prior exposure to and escape from host immune responses, clearly influences disease progression. Some escape mutations impact replicative capacity, while others effectively cloak the virus from the newly infected host's immune response by preventing recognition. It is the balance between the impact of escape mutations on viral fitness and susceptibility to the host immunogenetics that defines HIV-1 disease progression.

HIV and host immunogenetics: unraveling the role of HLA-C

Host genetic factors play a major role in determining the outcome of many infections including human immunodeficiency virus (HIV). Multiple host factors have been studied till date showing their varied role in susceptibility or resistance to HIV infection. HLA-C, however, has been recently started gaining interest in researchers mind revealing its polymorphisms to have an important effect on viral load set-points, disease progression as well as transmission. In this review report, we have compiled these significant findings of HLA-C in HIV infection, in an attempt to highlight the need for further research in the area in different ethnic population to establish its role in the infection.

Balance between transmitted HLA preadapted and nonassociated polymorphisms is a major determinant of HIV-1 disease progression

In a cohort of Zambian heterosexual transmission pairs, the authors show that HIV-1–transmitted variants already exhibit a significant degree of preadaptation to the new host's HLA alleles, which, modulated by polymorphisms that decrease viral fitness, determines early set-point VL and the rate of disease progression in the newly infected individual.

HIV-1 adapts to a new host through mutations that facilitate immune escape. Here, we evaluate the impact on viral control and disease progression of transmitted polymorphisms that were either preadapted to or nonassociated with the new host’s HLA. In a cohort of 169 Zambian heterosexual transmission pairs, we found that almost one-third of possible HLA-linked target sites in the transmitted virus Gag protein are already adapted, and that this transmitted preadaptation significantly reduced early immune recognition of epitopes. Transmitted preadapted and nonassociated polymorphisms showed opposing effects on set-point VL and the balance between the two was significantly associated with higher set-point VLs in a multivariable model including other risk factors. Transmitted preadaptation was also significantly associated with faster CD4 decline (<350 cells/µl) and this association was stronger after accounting for nonassociated polymorphisms, which were linked with slower CD4 decline. Overall, the relative ratio of the two classes of polymorphisms was found to be the major determinant of CD4 decline in a multivariable model including other risk factors. This study reveals that, even before an immune response is mounted in the new host, the balance of these opposing factors can significantly influence the outcome of HIV-1 infection.

Bayesian Inference of the Evolution of a Phenotype Distribution on a Phylogenetic Tree

The distribution of a phenotype on a phylogenetic tree is often a quantity of interest. Many phenotypes have imperfect heritability, so that a measurement of the phenotype for an individual can be thought of as a single realization from the phenotype distribution of that individual. If all individuals in a phylogeny had the same phenotype distribution, measured phenotypes would be randomly distributed on the tree leaves. This is, however, often not the case, implying that the phenotype distribution evolves over time. Here we propose a new model based on this principle of evolving phenotype distribution on the branches of a phylogeny, which is different from ancestral state reconstruction where the phenotype itself is assumed to evolve. We develop an efficient Bayesian inference method to estimate the parameters of our model and to test the evidence for changes in the phenotype distribution. We use multiple simulated data sets to show that our algorithm has good sensitivity and specificity properties. Since our method identifies branches on the tree on which the phenotype distribution has changed, it is able to break down a tree into components for which this distribution is unique and constant. We present two applications of our method, one investigating the association between HIV genetic variation and human leukocyte antigen and the other studying host range distribution in a lineage of Salmonella enterica, and we discuss many other potential applications.

Impact of pre-adapted HIV transmission

Human leukocyte antigen class I (HLA)-restricted CD8(+) T lymphocyte (CTL) responses are crucial to HIV-1 control. Although HIV can evade these responses, the longer-term impact of viral escape mutants remains unclear, as these variants can also reduce intrinsic viral fitness. To address this, we here developed a metric to determine the degree of HIV adaptation to an HLA profile. We demonstrate that transmission of viruses that are pre-adapted to the HLA molecules expressed in the recipient is associated with impaired immunogenicity, elevated viral load and accelerated CD4(+) T cell decline. Furthermore, the extent of pre-adaptation among circulating viruses explains much of the variation in outcomes attributed to the expression of certain HLA alleles. Thus, viral pre-adaptation exploits 'holes' in the immune response. Accounting for these holes may be key for vaccine strategies seeking to elicit functional responses from viral variants, and to HIV cure strategies that require broad CTL responses to achieve successful eradication of HIV reservoirs.

Differences in the Selection Bottleneck between Modes of Sexual Transmission Influence the Genetic Composition of the HIV-1 Founder Virus

Due to the stringent population bottleneck that occurs during sexual HIV-1 transmission, systemic infection is typically established by a limited number of founder viruses. Elucidation of the precise forces influencing the selection of founder viruses may reveal key vulnerabilities that could aid in the development of a vaccine or other clinical interventions. Here, we utilize deep sequencing data and apply a genetic distance-based method to investigate whether the mode of sexual transmission shapes the nascent founder viral genome. Analysis of 74 acute and early HIV-1 infected subjects revealed that 83% of men who have sex with men (MSM) exhibit a single founder virus, levels similar to those previously observed in heterosexual (HSX) transmission. In a metadata analysis of a total of 354 subjects, including HSX, MSM and injecting drug users (IDU), we also observed no significant differences in the frequency of single founder virus infections between HSX and MSM transmissions. However, comparison of HIV-1 envelope sequences revealed that HSX founder viruses exhibited a greater number of codon sites under positive selection, as well as stronger transmission indices possibly reflective of higher fitness variants. Moreover, specific genetic “signatures” within MSM and HSX founder viruses were identified, with single polymorphisms within gp41 enriched among HSX viruses while more complex patterns, including clustered polymorphisms surrounding the CD4 binding site, were enriched in MSM viruses. While our findings do not support an influence of the mode of sexual transmission on the number of founder viruses, they do demonstrate that there are marked differences in the selection bottleneck that can significantly shape their genetic composition. This study illustrates the complex dynamics of the transmission bottleneck and reveals that distinct genetic bottleneck processes exist dependent upon the mode of HIV-1 transmission.

While the global spread of HIV-1 has been fueled by sexual transmission the genetic determinants underlying the transmission bottleneck remains poorly understood. Here we characterized founder virus population diversity from next generation sequencing data in a cohort of 74 acute and early HIV-1 infected individuals. We observe that the risk of multi-variant infection in men-who-have-sex-with-men (MSM) is not greater than that observed for heterosexuals (HSX), contrary to reports of higher rates of multiple founder virus infections in higher-risk MSM transmissions. These findings were further supported through a metadata analysis of 354 acute and early HIV-1 subjects. We did, however, observe differences between HSM and MSM founder viruses, including a higher selection barrier in HSX transmission with founder viruses being more cohort consensus-like that may be reflective of increased replicative fitness. We also identified a number of residues within Envelope that behave in a risk-dependent manner and could be key for HIV-1 transmission. These novel insights improve our understanding of the HIV-1 transmission bottleneck and underscore the differential selective pressures that founder viruses within the two major transmission risk groups are subjected to.

A strong association of human leukocyte antigen-associated Pol and Gag mutations with clinical parameters in HIV-1 subtype A/E infection

OBJECTIVES

Identification of human leukocyte antigen-associated HIV-1 polymorphisms (HLA-APs) in different global populations furthers our understanding of HIV-1 pathogenesis and may help identify candidate immunogens for HIV vaccines targeted to these populations. Although numerous population-based studies identifying HLA-APs have been conducted in HIV-1 subtype B- and subtype C-infected cohorts, few have focused on subtype A/E.

DESIGN

We investigated HLA-APs in a cohort of chronically HIV-1 subtype A/E-infected Vietnamese individuals.

METHODS

HLA-APs in HIV-1 Gag, Pol, and Nef regions from 388 treatment-naive individuals chronically infected with HIV-1 subtype A/E were analyzed using phylogenetically informed approaches.

RESULTS

A total of 303 HLA-APs were identified. HLA-APs occurring at six positions in Gag and six positions in Pol were significantly associated with higher plasma viral load (pVL), whereas HLA-APs occurring at two positions in Gag and 13 positions in Pol were significantly associated with lower CD4 T-cell counts. Furthermore, the proportion of Pol codons harboring an HLA-AP specific to the host's HLA correlated positively with HIV-1 pVL (R = 0.22; P < 0.0001) and inversely with CD4 T-cell counts (R = -0.32; P < 0.0001). Similarly, the proportion of HLA-associated Gag codons harboring host-specific HLA-AP correlated inversely with CD4 T-cell counts (R = -0.13; P = 0.01).

CONCLUSION

These significant associations between HIV-1 amino acids adapted to Vietnamese HLA alleles and higher pVL and lower CD4 T-cell counts suggests that accumulation of cytotoxic T cells escape mutations may influence clinical outcomes in HIV-1 subtype A/E-infected Vietnamese individuals.

Role of HLA Adaptation in HIV Evolution

Killing of HIV-infected cells by CD8⁺ T-cells imposes strong selection pressure on the virus toward escape. The HLA class I molecules that are successful in mediating some degree of control over the virus are those that tend to present epitopes in conserved regions of the proteome, such as in p24 Gag, in which escape also comes at a significant cost to viral replicative capacity (VRC). In some instances, compensatory mutations can fully correct for the fitness cost of such an escape variant; in others, correction is only partial. The consequences of these events within the HIV-infected host, and at the population level following transmission of escape variants, are discussed. The accumulation of escape mutants in populations over the course of the epidemic already shows instances of protective HLA molecules losing their impact, and in certain cases, a modest decline in HIV virulence in association with population-level increase in mutants that reduce VRC.

Population-Level Immune-Mediated Adaptation in HIV-1 Polymerase during the North American Epidemic

Human leukocyte antigen (HLA) class I-associated polymorphisms in HIV-1 that persist upon transmission to HLA-mismatched hosts may spread in the population as the epidemic progresses. Transmission of HIV-1 sequences containing such adaptations may undermine cellular immune responses to the incoming virus in future hosts. Building upon previous work, we investigated the extent of HLA-associated polymorphism accumulation in HIV-1 polymerase (Pol) through comparative analysis of linked HIV-1/HLA class I genotypes sampled during historic (1979 to 1989; n = 338) and modern (2001 to 2011; n = 278) eras from across North America (Vancouver, BC, Canada; Boston, MA; New York, NY; and San Francisco, CA). Phylogenies inferred from historic and modern HIV-1 Pol sequences were star-like in shape, with an inferred most recent common ancestor (epidemic founder virus) sequence nearly identical to the modern North American subtype B consensus sequence. Nevertheless, modern HIV-1 Pol sequences exhibited roughly 2-fold-higher patristic (tip-to-tip) genetic distances than historic sequences, with HLA pressures likely driving ongoing diversification. Moreover, the frequencies of published HLA-associated polymorphisms in individuals lacking the selecting HLA class I allele was on average ∼2.5-fold higher in the modern than in the historic era, supporting their spread in circulation, though some remained stable in frequency during this time. Notably, polymorphisms restricted by protective HLA alleles appear to be spreading to a greater relative extent than others, though these increases are generally of modest absolute magnitude. However, despite evidence of polymorphism spread, North American hosts generally remain at relatively low risk of acquiring an HIV-1 polymerase sequence substantially preadapted to their HLA profiles, even in the present era.

IMPORTANCE

HLA class I-restricted cytotoxic T-lymphocyte (CTL) escape mutations in HIV-1 that persist upon transmission may accumulate in circulation over time, potentially undermining host antiviral immunity to the transmitted viral strain. We studied >600 experimentally collected HIV-1 polymerase sequences linked to host HLA information dating back to 1979, along with phylogenetically reconstructed HIV-1 sequences dating back to the virus' introduction into North America. Overall, our results support the gradual spread of many-though not all-HIV-1 polymerase immune escape mutations in circulation over time. This is consistent with recent observations from other global regions, though the extent of polymorphism accumulation in North America appears to be lower than in populations with high seroprevalence, older epidemics, and/or limited HLA diversity. Importantly, the risk of acquiring an HIV-1 polymerase sequence at transmission that is substantially preadapted to one's HLA profile remains relatively low in North America, even in the present era.