An algorithm for mapping positively selected members of quasispecies-type viruses

Genetic variability of avian influenza virus subtype H5N8 in Egypt in 2017 and 2018

Since the incursion of avian influenza virus subtype H5N8 in Egypt in late 2016, it has spread rapidly, causing severe losses in poultry production. Multiple introductions of different reassorted strains were observed in 2017. In this study, a genetic characterization of the HA gene was carried out with 31 isolates selected from different governorates and sectors. Fifteen isolates were selected for NA gene sequence analysis. The HA and NA genes were divided into two subgroups (I and II) with positive selection pressure identified at positions 174 and 29, respectively. The HA gene contained two novel mutations in the antigenic sites, A and E. The HA nucleotide sequence identity ranged from 77 to 90% with different vaccine seeds. Full-genome sequence analysis was carried out for eight viruses, representing different governorates and sectors, to identify the predominant reassorted strain in Egypt. All viruses were similar to a reassorted strain of clade 2.3.4.4b that has been identified in Germany, among other countries. Analysis of these viruses revealed mutations specific to Egyptian strains and not the original virus characterized in 2017 (A/duck/Egypt/F446/2017), with a novel antiviral resistance marker, V27A, indicating resistance to amantadine in the M2 protein of two strains. The results indicate increased variability of circulating H5N8 viruses compared to earlier viruses sequenced in 2016 and 2017. The predominant reassorted virus circulating in 2017 and 2018 originated from an early 2017 strain. It is important to continue this surveillance of avian influenza viruses to monitor the evolution of circulating viruses.

Broad and Potent Neutralizing Antibodies Recognize the Silent Face of the HIV Envelope

Broadly neutralizing antibodies (bNAbs) against HIV-1 envelope (Env) inform vaccine design and are potential therapeutic agents. We identified SF12 and related bNAbs with up to 62% neutralization breadth from an HIV-infected donor. SF12 recognized a glycan-dominated epitope on Env's silent face and was potent against clade AE viruses, which are poorly covered by V3-glycan bNAbs. A 3.3Å cryo-EM structure of a SF12-Env trimer complex showed additional contacts to Env protein residues by SF12 compared with VRC-PG05, the only other known donor-derived silentface antibody, explaining SF12's increased neutralization breadth, potency, and resistance to Env mutation routes. Asymmetric binding of SF12 was associated with distinct N-glycan conformations across Env protomers, demonstrating intra-Env glycan heterogeneity. Administrating SF12 to HIV-1-infected humanized mice suppressed viremia and selected for viruses lacking the N448gp120 glycan. Effective bNAbs can therefore be raised against HIV-1 Env's silent face, suggesting their potential for HIV-1 prevention, therapy, and vaccine development.

Identification and Characterization of Positively Selected Mutations in Nef of Four HIV-1 Major Subtypes from Los Alamos National Laboratory

BACKGROUND

Human immunodeficiency virus-1 (HIV-1) mutates rapidly to escape host immune pressure. This results in the generation of positively selected mutations (PSM) throughout the viral genome. Escape mutations in Nef, one of the accessory proteins of HIV-1, which plays an important role in viral pathogenicity have previously been identified in several large cohort studies, but the evolution of PSMs overtime in various HIV-1 subtypes remains unknown.

METHODS

161 clade A1, 3093 clade B, 647 clade C and 115 clade D HIV-1 nef sequences were obtained from the HIV Database of Los Alamos National Laboratory and aligned using MEGA 6.0. The sequences from each clade were grouped based on the year of collection. Quasi analysis was used to identify PSMs and the number and locations of PSMs were compared among different subtypes.

RESULTS

PSMs for all four subtypes were distributed across the sequence of Nef, and conserved residues F90, W113, PxxPxR (a.a 72-77) remain unaltered overtime. The frequency of PSMs was stable among subtype B sequences but increased overtime for other subtypes. Phylogenetic analysis shows that sequences containing PSMs tend to cluster together at both inter and intra- subtype levels.

CONCLUSION

Identification of PSMs and their changes overtime within various subtypes of HIV-1 is important in defining global viral evolutionary patterns that can provide insights for designing therapeutic strategies.

A Neutralizing Antibody Recognizing Primarily N-Linked Glycan Targets the Silent Face of the HIV Envelope

Virtually the entire surface of the HIV-1-envelope trimer is recognized by neutralizing antibodies, except for a highly glycosylated region at the center of the "silent face" on the gp120 subunit. From an HIV-1-infected donor, #74, we identified antibody VRC-PG05, which neutralized 27% of HIV-1 strains. The crystal structure of the antigen-binding fragment of VRC-PG05 in complex with gp120 revealed an epitope comprised primarily of N-linked glycans from N262, N295, and N448 at the silent face center. Somatic hypermutation occurred preferentially at antibody residues that interacted with these glycans, suggesting somatic development of glycan recognition. Resistance to VRC-PG05 in donor #74 involved shifting of glycan-N448 to N446 or mutation of glycan-proximal residue E293. HIV-1 neutralization can thus be achieved at the silent face center by glycan-recognizing antibody; along with other known epitopes, the VRC-PG05 epitope completes coverage by neutralizing antibody of all major exposed regions of the prefusion closed trimer.

Genetic background matters: a plant-virus gene-for-gene interaction is strongly influenced by genetic contexts

Evolutionary processes responsible for parasite adaptation to their hosts determine our capacity to manage sustainably resistant plant crops. Most plant-parasite interactions studied so far correspond to gene-for-gene models in which the nature of the alleles present at a plant resistance locus and at a pathogen pathogenicity locus determine entirely the outcome of their confrontation. The interaction between the pepper pvr2 resistance locus and Potato virus Y (PVY) genome-linked protein VPg locus obeys this kind of model. Using synthetic chimeras between two parental PVY cDNA clones, we showed that the viral genetic background surrounding the VPg pathogenicity locus had a strong impact on the resistance breakdown capacity of the virus. Indeed, recombination of the cylindrical inclusion (CI) coding region between two PVY cDNA clones multiplied by six the virus capacity to break down the pvr2(3) -mediated resistance. High-throughput sequencing allowed the exploration of the diversity of PVY populations in response to the selection pressure of the pvr2(3) resistance. The CI chimera, which possessed an increased resistance breakdown capacity, did not show an increased mutation accumulation rate. Instead, selection of the most frequent resistance-breaking mutation seemed to be more efficient for the CI chimera than for the parental virus clone. These results echoed previous observations, which showed that the plant genetic background in which the pvr2(3) resistance gene was introduced modified strongly the efficiency of selection of resistance-breaking mutations by PVY. In a broader context, the PVY CI coding region is one of the first identified genetic factors to determine the evolvability of a plant virus.

Expression, characterization and structural modelling of a feruloyl esterase from the thermophilic fungus Myceliophthora thermophila

A ferulic acid esterase (FAE) from the thermophilic fungus Myceliophthora thermophila (synonym Sporotrichum thermophile), belonging to the carbohydrate esterase family 1 (CE-1), was functionally expressed in methylotrophic yeast Pichia pastoris. The putative FAE from the genomic DNA was successfully cloned in P. pastoris X-33 to confirm that the enzyme exhibits FAE activity. The recombinant FAE was purified to its homogeneity (39 kDa) and subsequently characterized using a series of model substrates including methyl esters of hydroxycinnamates, alkyl ferulates and monoferuloylated 4-nitrophenyl glycosides. The substrate specificity profiling reveals that the enzyme shows a preference for the hydrolysis of methyl caffeate and p-coumarate and a strong preference for the hydrolysis of n-butyl and iso-butyl ferulate. The enzyme was active on substrates containing ferulic acid ester linked to the C-5 and C-2 linkages of arabinofuranose, whilst it was found capable of de-esterifying acetylated glucuronoxylans. Ferulic acid (FA) was efficiently released from destarched wheat bran when the esterase was incubated together with an M3 xylanase from Trichoderma longibrachiatum (a maximum of 41% total FA released after 1 h incubation). Prediction of the secondary structure of MtFae1a was performed in the PSIPRED server whilst modelling the 3D structure was accomplished by the use of the HH 3D structure prediction server.

QUASI analysis of host immune responses to Gag polyproteins of human immunodeficiency virus type 1 by a systematic bioinformatics approach

There is a consensus that Gag-specific cytotoxic T lymphocyte (CTL) response plays a key role in the immune control of human immunodeficiency virus type 1 (HIV-1) infection. In this study, we analyzed all currently available gag sequences in the Los Alamos HIV sequence database and identified positive selection (PS) sites likely restricted by the host immune responses. We found that between 23.4% and 47.4% of PS sites were shared by clades A, B, and C of Gag, indicating similar positive selection pressure on Gag in different subtypes of HIV-1. Furthermore, a significant correlation was observed between the combined CTL and antibody responses and PS sites. The Gag regions of free from PS contained 9 CTL epitopes restricted by 11 HLA class I alleles associated with disease progression to acquired immune deficiency syndrome (AIDS). These analyses provide information important for the identification of cross-clade epitopes and development of a global HIV-1 vaccine.

Role of positive selection pressure on the evolution of H5N1 hemagglutinin

The surface glycoprotein hemagglutinin (HA) helps the influenza A virus to evade the host immune system by antigenic variation and is a major driving force for viral evolution. In this study, the selection pressure on HA of H5N1 influenza A virus was analyzed using bioinformatics algorithms. Most of the identified positive selection (PS) sites were found to be within or adjacent to epitope sites. Some of the identified PS sites are consistent with previous experimental studies, providing further support to the biological significance of our findings. The highest frequency of PS sites was observed in recent strains isolated during 2005-2007. Phylogenetic analysis was also conducted on HA sequences from various hosts. Viral drift is almost similar in both avian and human species with a progressive trend over the years. Our study reports new mutations in functional regions of HA that might provide markers for vaccine design or can be used to predict isolates of pandemic potential.

Identification and characterization of HLA-A*0301 epitopes in HIV-1 gag proteins using a novel approach

Identification of CTL epitopes correlated to immune protection is important for the development of vaccines that enhance T-cell mediated immune responses. The correlation of positively selected amino acids (PS) of HIV-1 with host HLA alleles can identify regions containing potential T-cell epitopes. However, the specific epitopes have to be identified and characterized using overlapping peptides through T-cell functional assays. In this study we used a new approach to identify and characterize potential epitopes in the gag region containing PS mutations that significantly correlated with HLA-A*0301. The iTopia Epitope Discovery System was used to rapidly screen a panel of peptides overlapping the regions containing PS mutations and the peptides identified were assessed for relative affinity and complex stability. The potential epitopes were then validated by interferon gamma (IFN-gamma) ELISpot assays with patient PBMCs. Using this approach we identified/confirmed the predicted HLA-A*0301 epitopes in two regions of gag containing PS mutations V7I and K403R, one previously reported and the other novel. Five of the seven peptides that bound to A*0301 contained the K403R mutation and corresponded to the documented LARNCRAPRK-A3 supertype epitope. Two epitope variants, RASVLSGGK and RASILSGGK containing the V7I mutation, were identified using the iTopia Epitope Discovery System, however only the consensus variant (RAK9C) was confirmed using the ELISpot assay and it represents a novel A*0301 epitope.

Multiple T-cell epitopes overlap positively-selected residues in the p1 spacer protein of HIV-1 gag

OBJECTIVES

The p1 region of HIV-1 gag contains the frameshift stem-loop, gag-pol transframe and a protease cleavage site that are crucial for viral assembly, replication and infectivity. Identifying and characterizing CD8+ epitopes that are under host immune selection in this region will help in designing effective vaccines for HIV-1.

DESIGN

An approach combining bioinformatical analysis and interferon gamma enzyme-linked immunosorbent spot (ELISPOT) assays is used to identify and characterize the epitopes. Potential human leukocyte antigen (HLA)-restricted epitopes were identified by correlating the positively-selected mutations with host HLA alleles.

METHODS

ELISPOT analysis with overlapping peptides was used to confirm and characterize the epitopes.

RESULTS

Four positively-selected residues were significantly associated with HLA class I alleles, including HLA B*1302 (K4R, P = 0.0008 and I5L, P = 0.0108), A*7401 (S9N, P = 0.0002) and A*30 genotypes (P7S, P = 0.009), suggesting epitopes restricted by these alleles are present in this region. ELISPOT analysis with patient peripheral blood mononuclear cells (PBMCs) identified seven novel epitopes restricted by the 3 alleles. Two types of epitopes were observed in this region based on the ELISPOT responses, Type I: the positively-selected variation does not affect CD8+ T-cell responses; and Type II: the CD8+ T-cell responses are determined by the epitope variants.

CONCLUSION

We identified and characterized seven novel CD8+ epitopes in the p1 spacer protein region. Classifying the effects of positively-selected variants on CD8+ T-cell responses will help in designing effective vaccines for HIV-1.

Prediction of protein structural class using a complexity-based distance measure

Knowledge of structural class plays an important role in understanding protein folding patterns. So it is necessary to develop effective and reliable computational methods for prediction of protein structural class. To this end, we present a new method called NN-CDM, a nearest neighbor classifier with a complexity-based distance measure. Instead of extracting features from protein sequences as done previously, distance between each pair of protein sequences is directly evaluated by a complexity measure of symbol sequences. Then the nearest neighbor classifier is adopted as the predictive engine. To verify the performance of this method, jackknife cross-validation tests are performed on several benchmark datasets. Results show that our approach achieves a high prediction accuracy over some classical methods.

An integrative bioinformatic approach for studying escape mutations in human immunodeficiency virus type 1 gag in the Pumwani Sex Worker Cohort

Human immunodeficiency virus type 1 (HIV-1) is able to evade the host cytotoxic T-lymphocyte (CTL) response through a variety of escape avenues. Epitopes that are presented to CTLs are first processed in the presenting cell in several steps, including proteasomal cleavage, transport to the endoplasmic reticulum, binding by the HLA molecule, and finally presentation to the T-cell receptor. An understanding of the potential of the virus to escape CTL responses can aid in designing an effective vaccine. To investigate such a potential, we analyzed HIV-1 gag from 468 HIV-1-positive Kenyan women by using several bioinformatic approaches that allowed the identification of positively selected amino acids in the HIV-1 gag region and study of the effects that these mutations could have on the various stages of antigen processing. Correlations between positively selected residues and mean CD4 counts also allowed study of the effect of mutation on HIV disease progression. A number of mutations that could create or destroy proteasomal cleavage sites or reduce binding affinity of the transport antigen processing protein, effectively hindering epitope presentation, were identified. Many mutations correlated with the presence of specific HLA alleles and with lower or higher CD4 counts. For instance, the mutation V190I in subtype A1-infected individuals is associated with HLA-B*5802 (P = 4.73 x 10(-4)), a rapid-progression allele according to other studies, and also to a decreased mean CD4 count (P = 0.019). Thus, V190I is a possible HLA escape mutant. This method classifies many positively selected mutations across the entire gag region according to their potential for immune escape and their effect on disease progression.

QUASI analysis of the HIV-1 envelope sequences in the Los Alamos National Laboratory HIV sequence database: pattern and distribution of positive selection sites and their frequencies over years

The envelope (env) protein of human immunodeficiency virus type 1 (HIV-1) plays a crucial role in virus entry and is a central target for HIV vaccine design. Using the QUASI program, we analyzed the conserved regions of all currently available env sequences in the Los Alamos National Laboratory HIV Sequence Database and identified positive selection (PS) sites that are likely to be restricted by host immune responses. We found that PS sites are dispersed across conserved regions of env sequence, and that the C3, C4, and C5 regions were the most targeted. Several regions were identified as being PS free and were mainly distributed in the C1 and C2 regions. When comparing individual QUASI PS site frequencies across clades and geographical regions with the overall frequency of the entire env database, the env sequences from North America showed significantly lower PS site frequency, while those from Asia were significantly higher using Student's t test. The QUASI PS site frequency of env proteins from viruses isolated from different years showed that the PS site frequencies of the env population increased over time. Our study provides an overview of PS sites across the conserved regions of HIV-1 env sequences.

Problems in using statistical analysis of replacement and silent mutations in antibody genes for determining antigen-driven affinity selection

The analysis of molecular signatures of antigen-driven affinity selection of B cells is of immense use in studies on normal and abnormal B cell development. Most of the published literature compares the expected and observed frequencies of replacement (R) and silent (S) mutations in the complementarity-determining regions (CDRs) and the framework regions (FRs) of antibody genes to identify the signature of antigenic selection. The basic assumption of this statistical method is that antigenic selection creates a bias for R mutations in the CDRs and for S mutations in the FRs. However, it has been argued that the differences in intrinsic mutability among different regions of an antibody gene can generate a statistically significant bias even in the absence of any antigenic selection. We have modified the existing statistical method to include the effects of intrinsic mutability of different regions of an antibody gene. We used this method to analyse sequences of several B cell-derived monoclonals against T-dependent antigens, T-independent antigens, clones derived from lymphoma and amyloidogenic clones. Our sequence analysis indicates that even after correcting for the intrinsic mutability of antibody genes, statistical parameters fail to reflect the role of antigen-driven affinity selection in maturation of many clones. We suggest that, contrary to the basic assumption of such statistical methods, selection can act both for and against R mutations in the CDR as well as in the FR regions. In addition we have identified different methodological difficulties in the current uses of such statistical analysis of antibody genes.

Presence of Intact vpu and nef genes in nonpathogenic SHIV is essential for acquisition of pathogenicity of this virus by serial passage in macaques

Use of the macaque model of human immunodeficiency virus (HIV) pathogenesis has shown that the accessory genes nef and vpu are important in the pathogenicity of simian immunodeficiency virus (SIV) and simian-human immunodeficiency virus (SHIV). We examined the ability of two nonpathogenic SHIVs, SHIV(PPC) and DeltavpuDeltanefSHIV(PPC), to gain pathogenicity by rapid serial passage in macaques. In this study, each virus was passaged by blood intravenously four times at 4-week intervals in macaques. Animals were monitored for 40 weeks for levels of CD4 T cells and quantitative measures of virus infection. DeltavpuDeltanefSHIV(PPC) maintained a limited phase of productive replication in the four animals, with no loss of CD4(+) T cells, whereas SHIV(PPC) became more pathogenic in later passages, judging by plasma viral load and viral mRNA in lymph nodes, infectious peripheral blood mononuclear cells and CD4(+) T cell loss. The nef, LTR, and env of the SHIV(PPC) viruses underwent numerous mutations, compared to DeltavpuDeltanefSHIV(PPC). This study confirms the seminal role that nef, LTR, and vpu could play in regulation of pathogenesis of HIV infection.