Viruses selectively mutate their CD8+ T-cell epitopes--a large-scale immunomic analysis

A Structural In Silico Analysis of the Immunogenicity of L-Asparaginase from Penicillium cerradense

L-asparaginase is an essential drug used to treat acute lymphoid leukemia (ALL), a cancer of high prevalence in children. Several adverse reactions associated with L-asparaginase have been observed, mainly caused by immunogenicity and allergenicity. Some strategies have been adopted, such as searching for new microorganisms that produce the enzyme and applying protein engineering. Therefore, this work aimed to elucidate the molecular structure and predict the immunogenic profile of L-asparaginase from Penicillium cerradense, recently revealed as a new fungus of the genus Penicillium and producer of the enzyme, as a motivation to search for alternatives to bacterial L-asparaginase. In the evolutionary relationship, L-asparaginase from P. cerradense closely matches Aspergillus species. Using in silico tools, we characterized the enzyme as a protein fragment of 378 amino acids (39 kDa), including a signal peptide containing 17 amino acids, and the isoelectric point at 5.13. The oligomeric state was predicted to be a homotetramer. Also, this L-asparaginase presented a similar immunogenicity response (T- and B-cell epitopes) compared to Escherichia coli and Dickeya chrysanthemi enzymes. These results suggest a potentially useful L-asparaginase, with insights that can drive strategies to improve enzyme production.

The impact of HLA polymorphism on herpesvirus infection and disease

Human Leukocyte Antigens (HLA) are cell surface molecules, central in coordinating innate and adaptive immune responses, that are targets of strong diversifying natural selection by pathogens. Of these pathogens, human herpesviruses have a uniquely ancient relationship with our species, where coevolution likely has reciprocating impact on HLA and viral genomic diversity. Consistent with this notion, genetic variation at multiple HLA loci is strongly associated with modulating immunity to herpesvirus infection. Here, we synthesize published genetic associations of HLA with herpesvirus infection and disease, both from case/control and genome-wide association studies. We analyze genetic associations across the eight human herpesviruses and identify HLA alleles that are associated with diverse herpesvirus-related phenotypes. We find that whereas most HLA genetic associations are virus- or disease-specific, HLA-A*01 and HLA-A*02 allotypes may be more generally associated with immune susceptibility and control, respectively, across multiple herpesviruses. Connecting genetic association data with functional corroboration, we discuss mechanisms by which diverse HLA and cognate receptor allotypes direct variable immune responses during herpesvirus infection and pathogenesis. Together, this review examines the complexity of HLA-herpesvirus interactions driven by differential T cell and Natural Killer cell immune responses.

Predicting CD4 T-cell epitopes based on antigen cleavage, MHCII presentation, and TCR recognition

Accurate predictions of T-cell epitopes would be useful for designing vaccines, immunotherapies for cancer and autoimmune diseases, and improved protein therapies. The humoral immune response involves uptake of antigens by antigen presenting cells (APCs), APC processing and presentation of peptides on MHC class II (pMHCII), and T-cell receptor (TCR) recognition of pMHCII complexes. Most in silico methods predict only peptide-MHCII binding, resulting in significant over-prediction of CD4 T-cell epitopes. We present a method, ITCell, for prediction of T-cell epitopes within an input protein antigen sequence for given MHCII and TCR sequences. The method integrates information about three stages of the immune response pathway: antigen cleavage, MHCII presentation, and TCR recognition. First, antigen cleavage sites are predicted based on the cleavage profiles of cathepsins S, B, and H. Second, for each 12-mer peptide in the antigen sequence we predict whether it will bind to a given MHCII, based on the scores of modeled peptide-MHCII complexes. Third, we predict whether or not any of the top scoring peptide-MHCII complexes can bind to a given TCR, based on the scores of modeled ternary peptide-MHCII-TCR complexes and the distribution of predicted cleavage sites. Our benchmarks consist of epitope predictions generated by this algorithm, checked against 20 peptide-MHCII-TCR crystal structures, as well as epitope predictions for four peptide-MHCII-TCR complexes with known epitopes and TCR sequences but without crystal structures. ITCell successfully identified the correct epitopes as one of the 20 top scoring peptides for 22 of 24 benchmark cases. To validate the method using a clinically relevant application, we utilized five factor VIII-specific TCR sequences from hemophilia A subjects who developed an immune response to factor VIII replacement therapy. The known HLA-DR1-restricted factor VIII epitope was among the six top-scoring factor VIII peptides predicted by ITCall to bind HLA-DR1 and all five TCRs. Our integrative approach is more accurate than current single-stage epitope prediction algorithms applied to the same benchmarks. It is freely available as a web server (http://salilab.org/itcell).

HLA class I haplotype diversity is consistent with selection for frequent existing haplotypes

The major histocompatibility complex (MHC) contains the most polymorphic genetic system in humans, the human leukocyte antigen (HLA) genes of the adaptive immune system. High allelic diversity in HLA is argued to be maintained by balancing selection, such as negative frequency-dependent selection or heterozygote advantage. Selective pressure against immune escape by pathogens can maintain appreciable frequencies of many different HLA alleles. The selection pressures operating on combinations of HLA alleles across loci, or haplotypes, have not been extensively evaluated since the high HLA polymorphism necessitates very large sample sizes, which have not been available until recently. We aimed to evaluate the effect of selection operating at the HLA haplotype level by analyzing HLA A~C~B~DRB1~DQB1 haplotype frequencies derived from over six million individuals genotyped by the National Marrow Donor Program registry. In contrast with alleles, HLA haplotype diversity patterns suggest purifying selection, as certain HLA allele combinations co-occur in high linkage disequilibrium. Linkage disequilibrium is positive (D_ij'>0) among frequent haplotypes and negative (D_ij'<0) among rare haplotypes. Fitting the haplotype frequency distribution to several population dynamics models, we found that the best fit was obtained when significant positive frequency-dependent selection (FDS) was incorporated. Finally, the Ewens-Watterson test of homozygosity showed excess homozygosity for 5-locus haplotypes within 23 US populations studied, with an average Fnd of 28.43. Haplotype diversity is most consistent with purifying selection for HLA Class I haplotypes (HLA-A, -B, -C), and was not inferred for HLA Class II haplotypes (-DRB1 and—DQB1). We discuss our empirical results in the context of evolutionary theory, exploring potential mechanisms of selection that maintain high linkage disequilibrium in MHC haplotype blocks.

The adaptive immune system presents antigens derived from pathogenic and normal self proteins on the cell surface using human leukocyte antigen (HLA) molecules. The HLA loci coding for these molecules are found in major histocompatibility complex (MHC) region, the most polymorphic region in the human genome, with over 15,000 HLA alleles observed so far in the world population. A high frequency of many different HLA alleles is thought be sustained by balancing selection. New HLA alleles may have an advantage over existing frequent alleles since immune escape mutations in pathogens within a population are maintained primarily in epitopes presented on frequent HLA alleles. Host immune function is not determined by single HLA alleles, but by both copies of autosomal HLA genes together (genotypes). Complementarity in function across the two potentially-variant copies of HLA at each locus can result in overdominance and heterozygote advantage at the genotype level. Less explored are selection mechanisms that may be operating across combinations of HLA alleles across loci (haplotypes). Indeed, in addition to high allelic diversity, HLA also has distinctive patterns of haplotype diversity, as certain HLA alleles co-occur in high linkage disequilibrium across five classical HLA loci (HLA-A, -B, -C, -DRB1, -DQB1). We applied multiple population genetic models to a dataset of HLA haplotype frequencies derived from over six million individuals with the goal of determining what type of selection may impact HLA haplotype diversity. We found frequent haplotypes were preferentially maintained in the population across 23 US populations studied. Thus, balancing selection at the allele level and purifying selection at the haplotype level may together affect HLA diversity in human populations.

Converging evolution leads to near maximal junction diversity through parallel mechanisms in B and T cell receptors

T and B cell receptor (TCR and BCR) complementarity determining region 3 (CDR3) genetic diversity is produced through multiple diversification and selection stages. Potential holes in the CDR3 repertoire were argued to be linked to immunodeficiencies and diseases. In contrast with BCRs, TCRs have practically no Dβ germline genetic diversity, and the question emerges as to whether they can produce a diverse CDR3 repertoire. In order to address the genetic diversity of the adaptive immune system, appropriate quantitative measures for diversity and large-scale sequencing are required. Such a diversity method should incorporate the complex diversification mechanisms of the adaptive immune response and the BCR and TCR loci structure. We combined large-scale sequencing and diversity measures to show that TCRs have a near maximal CDR3 genetic diversity. Specifically, TCR have a larger junctional and V germline diversity, which starts more 5' in Vβ than BCRs. Selection decreases the TCR repertoire diversity, but does not affect BCR repertoire. As a result, TCR is as diverse as BCR repertoire, with a biased CDR3 length toward short TCRs and long BCRs. These differences suggest parallel converging evolutionary tracks to reach the required diversity to avoid holes in the CDR3 repertoire.

Immunological evasion of immediate-early varicella zoster virus proteins

The varicella zoster virus (VZV) causes the childhood disease commonly known as chickenpox and can later in life reactivate as herpes zoster. The adaptive immune system is known to play an important role in suppressing VZV reactivation. A central aspect of this system is the presentation of VZV-derived peptides by the major histocompatibility complex (MHC) proteins. Here, we investigate if key VZV proteins have evolved their amino acid sequence to avoid presentation by MHC based on predictive models of MHC-peptide affinity. This study shows that the immediate-early proteins of all characterized VZV strains are profoundly depleted for high-affinity MHC-I-restricted epitopes. The same depletion can be found in its closest animal analog, the simian varicella virus. Further orthology analysis towards other herpes viruses suggests that the protein amino acid frequency is one of the primary drivers of targeted epitope depletion.

T cell epitope redundancy: cross-conservation of the TCR face between pathogens and self and its implications for vaccines and autoimmunity

T cells are extensively trained on 'self' in the thymus and then move to the periphery, where they seek out and destroy infections and regulate immune response to self-antigens. T cell receptors (TCRs) on T cells' surface recognize T cell epitopes, short linear strings of amino acids presented by antigen-presenting cells. Some of these epitopes activate T effectors, while others activate regulatory T cells. It was recently discovered that T cell epitopes that are highly conserved on their TCR face with human genome sequences are often associated with T cells that regulate immune response. These TCR-cross-conserved or 'redundant epitopes' are more common in proteins found in pathogens that have co-evolved with humans than in other non-commensal pathogens. Epitope redundancy might be the link between pathogens and autoimmune disease. This article reviews recently published data and addresses epitope redundancy, the "elephant in the room" for vaccine developers and T cell immunologists.

Estimate of within population incremental selection through branch imbalance in lineage trees

Incremental selection within a population, defined as limited fitness changes following mutation, is an important aspect of many evolutionary processes. Strongly advantageous or deleterious mutations are detected using the synonymous to non-synonymous mutations ratio. However, there are currently no precise methods to estimate incremental selection. We here provide for the first time such a detailed method and show its precision in multiple cases of micro-evolution. The proposed method is a novel mixed lineage tree/sequence based method to detect within population selection as defined by the effect of mutations on the average number of offspring. Specifically, we propose to measure the log of the ratio between the number of leaves in lineage trees branches following synonymous and non-synonymous mutations. The method requires a high enough number of sequences, and a large enough number of independent mutations. It assumes that all mutations are independent events. It does not require of a baseline model and is practically not affected by sampling biases. We show the method's wide applicability by testing it on multiple cases of micro-evolution. We show that it can detect genes and inter-genic regions using the selection rate and detect selection pressures in viral proteins and in the immune response to pathogens.

Immune camouflage: relevance to vaccines and human immunology

High strain sequence variability, interference with innate immune mechanisms, and epitope deletion are all examples of strategies that pathogens have evolved to subvert host defenses. To this list we would add another strategy: immune camouflage. Pathogens whose epitope sequences are cross-conserved with multiple human proteins at the TCR-facing residues may be exploiting “ignorance and tolerance," which are mechanisms by which mature T cells avoid immune responses to self-antigens. By adopting amino acid configurations that may be recognized by autologous regulatory T cells, pathogens may be actively suppressing protective immunity. Using the new JanusMatrix TCR-homology-mapping tool, we have identified several such ‘camouflaged’ tolerizing epitopes that are present in the viral genomes of pathogens such as emerging H7N9 influenza. Thus in addition to the overall low number of T helper epitopes that is present in H7 hemaglutinin (as described previously, see http://dx.doi.org/10.4161/hv.24939), the presence of such tolerizing epitopes in H7N9 could explain why, in recent vaccine trials, whole H7N9-HA was poorly immunogenic and associated with low seroconversion rates (see http://dx.doi.org/10.4161/hv.28135). In this commentary, we provide an overview of the immunoinformatics process leading to the discovery of tolerizing epitopes in pathogen genomic sequences, provide a brief summary of laboratory data that validates the discovery, and point the way forward. Removal of viral, bacterial and parasite tolerizing epitopes may permit researchers to develop more effective vaccines and immunotherapeutics in the future.

Hit-and-run, hit-and-stay, and commensal bacteria present different peptide content when viewed from the perspective of the T cell

Different types of bacteria face different pressures from the immune system, with those that persist ("hit-and-stay") potentially having to adapt more in order to escape than those prone to short-lived infection ("hit-and-run"), and with commensal bacteria potentially different from both due to additional physical mechanisms for avoiding immune detection. The Janus Immunogenicity Score (JIS) was recently developed to assess the likelihood of T cell recognition of an antigen, using an analysis that considers both binding of a peptide within the antigen by major histocompatability complex (MHC) and recognition of the peptide:MHC complex by cognate T cell receptor (TCR). This score was shown to be predictive of T effector vs. T regulatory or null responses in experimental data, as well as to distinguish viruses representative of the hit-and-stay vs. hit-and-run phenotypes. Here, JIS-based analyses were conducted in order to characterize the extent to which the pressure to avoid T cell recognition is manifested in genomic differences among representative hit-and-run, hit-and-stay, and commensal bacteria. Overall, extracellular proteins were found to have different JIS profiles from cytoplasmic ones. Contrasting the bacterial groups, extracellular proteins were shown to be quite different across the groups, much more so than intracellular proteins. The differences were evident even at the level of corresponding peptides in homologous protein pairs from hit-and-run and hit-and-stay bacteria. The multi-level analysis of patterns of immunogenicity across different groups of bacteria provides a new way to approach questions of bacterial immune camouflage or escape, as well as to approach the selection and optimization of candidates for vaccine design.

H7N9 T-cell epitopes that mimic human sequences are less immunogenic and may induce Treg-mediated tolerance

Avian-origin H7N9 influenza is a novel influenza A virus (IAV) that emerged in humans in China in 2013. Using immunoinformatics tools, we identified several H7N9 T cell epitopes with T cell receptor (TCR)-facing residues identical to those of multiple epitopes from human proteins. We hypothesized that host tolerance to these peptides may impair T helper response and contribute to the low titer, weak hemagglutination inhibiting (HI) antibody responses and diminished seroconversion rates that have been observed in human H7N9 infections and vaccine trials. We found that the magnitude of human T effector responses to individual H7N9 peptides was inversely correlated with the peptide's resemblance to self. Furthermore, a promiscuous T cell epitope from the hemagglutinin (HA) protein suppressed responses to other H7N9 peptides when co-administered in vitro. Along with other highly ‘human-like’ peptides from H7N9, this peptide was also shown to expand FoxP3⁺ regulatory T cells (Tregs). Thus, H7N9 may be camouflaged from effective human immune response by T cell epitope sequences that avert or regulate effector T cell responses through host tolerance.

Viral CD8 T cell epitope nucleotide composition shows evidence of short- and long-term evolutionary strategies

Viral epitopes have a distinct codon usage that reflects their dual role in infection and immunity. On the one hand, epitopes are part of proteins important to viral function; on the other hand, they are targets of the immune response. Studies of selection are most commonly based on changes of amino acid and seen through the accumulation of non-synonymous mutations. An independent measure of selection is the codon usage and underlying changeability of the nucleotide sequences. We here use multiple tools and a large-scale analysis of viral genomes to demonstrate that viral epitopes have a distinct codon usage and that this codon usage reflects distinct short- and long-term types of selection during viral evolution. We show that CD8(+) T cell epitopes are encoded by codons more distant from stop codons and more changeable than codons outside epitopes. This biased codon usage reflects the viral population toggling back and forth from a wild-type sequence to an escape mode, which enable them to avoid immune detection when needed, and go back to the functionally favorable form when the threat is removed (i.e., in a new host).

Integrated assessment of predicted MHC binding and cross-conservation with self reveals patterns of viral camouflage

Background

Immune recognition of foreign proteins by T cells hinges on the formation of a ternary complex sandwiching a constituent peptide of the protein between a major histocompatibility complex (MHC) molecule and a T cell receptor (TCR). Viruses have evolved means of "camouflaging" themselves, avoiding immune recognition by reducing the MHC and/or TCR binding of their constituent peptides. Computer-driven T cell epitope mapping tools have been used to evaluate the degree to which particular viruses have used this means of avoiding immune response, but most such analyses focus on MHC-facing 'agretopes'. Here we set out a new means of evaluating the TCR faces of viral peptides in addition to their agretopes, integrating evaluations of both sides of the ternary complex in a single analysis.

Methods

This paper develops what we call the Janus Immunogenicity Score (JIS), bringing together a well-established method for predicting MHC binding, with a novel assessment of the potential for TCR binding based on similarity with self. Intuitively, both good MHC binding and poor self-similarity are required for high immunogenicity (i.e., a robust T effector response).

Results

Focusing on the class II antigen-processing pathway, we show that the JIS of T effector epitopes and null or regulatory epitopes deposited in a large database of epitopes (Immune Epitope Database) are significantly different. We then show that different types of viruses display significantly different patterns of scores over their constituent peptides, with viruses causing chronic infection (Epstein-Barr and cytomegalovirus) strongly shifted to lower scores relative to those causing acute infection (Ebola and Marburg). Similarly we find distinct patterns among influenza proteins in H1N1 (a strain against which human populations rapidly developed immunity) and H5N1 and H7N9 (highly pathogenic avian flu strains, with significantly greater case mortality rates).

Conclusion

The Janus Immunogenicity Score, which integrates MHC binding and TCR cross-reactivity, provides a new tool for studying immunogenicity of pathogens and may improve the selection and optimization of antigenic elements for vaccine design.

Viral proteome size and CD8+ T cell epitope density are correlated: the effect of complexity on selection

•

We analyze the relation between viral complexity and their adaptation to the host immune system.

•

Viruses with few proteins and low number of nucleotides remove more CD8+ T cell epitopes.

•

Within a virus, short proteins (with fewer amino acids) adapt better than long ones.

•

The relation between total size and adaptation is host specific.

•

Complexity limits genetic adaptation in the high-mutation rate strong selection regime.

The relation between the complexity of organisms and proteins and their evolution rates has been discussed in the context of multiple generic models. The main robust claim from most such models is the negative relation between complexity and the accumulation rate of mutations.

Viruses accumulate escape mutations in their epitopes to avoid detection and destruction of their host cell by CD8+ T cells. The extreme regime of immune escape, namely, strong selection and high mutation rate, provide an opportunity to extend and validate the existing models of relation between complexity and evolution rate as proposed by Fisher and Kimura.

Using epitope prediction algorithms to compute the epitopes presented on the most frequent human HLA alleles in over 100 fully sequenced human viruses, and over 900 non-human viruses, we here study the correlation between viruses/proteins complexity (as measured by the number of proteins in the virus and the length of each protein, respectively) and the rate of accumulation of escape mutation. The latter is evaluated by measuring the normalized epitope density of viral proteins.

If the virus/protein complexity prevents the accumulation of escape mutations, the epitope density is expected to be positively correlated with both the number of proteins in the virus and the length of proteins. We show that such correlations are indeed observed for most human viruses. For non-human viruses the correlations were much less significant, indicating that the correlation is indeed induced by human HLA molecules.

Positional bias of MHC class I restricted T-cell epitopes in viral antigens is likely due to a bias in conservation

The immune system rapidly responds to intracellular infections by detecting MHC class I restricted T-cell epitopes presented on infected cells. It was originally thought that viral peptides are liberated during constitutive protein turnover, but this conflicts with the observation that viral epitopes are detected within minutes of their synthesis even when their source proteins exhibit half-lives of days. The DRiPs hypothesis proposes that epitopes derive from Defective Ribosomal Products (DRiPs), rather than degradation of mature protein products. One potential source of DRiPs is premature translation termination. If this is a major source of DRiPs, this should be reflected in positional bias towards the N-terminus. By contrast, if downstream initiation is a major source of DRiPs, there should be positional bias towards the C-terminus. Here, we systematically assessed positional bias of epitopes in viral antigens, exploiting the large set of data available in the Immune Epitope Database and Analysis Resource. We show a statistically significant degree of positional skewing among epitopes; epitopes from both ends of antigens tend to be under-represented. Centric-skewing correlates with a bias towards class I binding peptides being over-represented in the middle, in parallel with a higher degree of evolutionary conservation.

Evolution of viral life-cycle in response to cytotoxic T lymphocyte-mediated immunity

Viruses in mammals are constantly faced with the problem of elimination by the host immunity. Cytotoxic T lymphocyte (CTL) responses are thought to play a major role in the control and clearance of several viral infections in mice and humans. It is therefore expected that over evolutionary time, viruses would be forced to evolve to avoid recognition by CTLs. Indeed, a number of studies have documented the accumulation of viral variants with escape mutations. These mutations allow viruses to hide from CTL responses common in the host population. CTLs recognize viruses by short protein sequences, named epitopes, derived from viral proteins. The efficiency of viral recognition by epitope-specific CTL responses depends on the expression pattern of the proteins carrying these epitopes, and the total amount of that protein (and thus epitopes) in the cell. When a virus replicates in a cell, some viral genes are expressed early in the life cycle of the virus, while other proteins are expressed late. For example, HIV infected cells first express Rev and Tat proteins, and the Gag proteins are expressed late. Here we propose a dynamical model of the viral life cycle to study how expression level of early vs. late genes may affect viral dynamics within the host and virus transmission over the course of infection. We find that for acute and chronic viral infections lower expression of early genes than that of the late genes is expected to give selective advantage and higher transmission to viruses.

Comparing the similarity of different groups of bacteria to the human proteome

Numerous aspects of the relationship between bacteria and human have been investigated. One aspect that has recently received attention is sequence overlap at the proteomic level. However, there has not yet been a study that comprehensively characterizes the level of sequence overlap between bacteria and human, especially as it relates to bacterial characteristics like pathogenicity, G-C content, and proteome size. In this study, we began by performing a general characterization of the range of bacteria-human similarity at the proteomic level, and identified characteristics of the most- and least-similar bacterial species. We then examined the relationship between proteomic similarity and numerous other variables. While pathogens and nonpathogens had comparable similarity to the human proteome, pathogens causing chronic infections were found to be more similar to the human proteome than those causing acute infections. Although no general correspondence between a bacterium’s proteome size and its similarity to the human proteome was noted, no bacteria with small proteomes had high similarity to the human proteome. Finally, we discovered an interesting relationship between similarity and a bacterium’s G-C content. While the relationship between bacteria and human has been studied from many angles, their proteomic similarity still needs to be examined in more detail. This paper sheds further light on this relationship, particularly with respect to immunity and pathogenicity.

Bacteria modulate the CD8+ T cell epitope repertoire of host cytosol-exposed proteins to manipulate the host immune response

The main adaptive immune response to bacteria is mediated by B cells and CD4+ T-cells. However, some bacterial proteins reach the cytosol of host cells and are exposed to the host CD8+ T-cells response. Both gram-negative and gram-positive bacteria can translocate proteins to the cytosol through type III and IV secretion and ESX-1 systems, respectively. The translocated proteins are often essential for the bacterium survival. Once injected, these proteins can be degraded and presented on MHC-I molecules to CD8+ T-cells. The CD8+ T-cells, in turn, can induce cell death and destroy the bacteria's habitat. In viruses, escape mutations arise to avoid this detection. The accumulation of escape mutations in bacteria has never been systematically studied. We show for the first time that such mutations are systematically present in most bacteria tested. We combine multiple bioinformatic algorithms to compute CD8+ T-cell epitope libraries of bacteria with secretion systems that translocate proteins to the host cytosol. In all bacteria tested, proteins not translocated to the cytosol show no escape mutations in their CD8+ T-cell epitopes. However, proteins translocated to the cytosol show clear escape mutations and have low epitope densities for most tested HLA alleles. The low epitope densities suggest that bacteria, like viruses, are evolutionarily selected to ensure their survival in the presence of CD8+ T-cells. In contrast with most other translocated proteins examined, Pseudomonas aeruginosa's ExoU, which ultimately induces host cell death, was found to have high epitope density. This finding suggests a novel mechanism for the manipulation of CD8+ T-cells by pathogens. The ExoU effector may have evolved to maintain high epitope density enabling it to efficiently induce CD8+ T-cell mediated cell death. These results were tested using multiple epitope prediction algorithms, and were found to be consistent for most proteins tested.

Optimal viral immune surveillance evasion strategies

Following cell entry, viruses can be detected by cytotoxic T lymphocytes. These cytotoxic T lymphocytes can induce host cell apoptosis and prevent the propagation of the virus. Viruses with fewer epitopes have a higher survival probability, and are selected through evolution. However, mutations have a fitness cost and on evolutionary periods viruses maintain some epitopes. The number of epitopes in each viral protein is a balance between the selective advantage of having fewer epitopes and the reduced fitness following the epitope removing mutations. We discuss a bioinformatic analysis of the number of epitopes in various viral proteins and propose an optimization framework to explain these numbers. We show, using a genomic analysis and a theoretical optimization framework, that a critical factor affecting the number of presented epitopes is the expression stage in the viral life cycle of the gene coding for the protein. The early expression of epitopes can lead to the destruction of the host cell before budding can take place. We show that a lower number of epitopes is expected in early proteins even if late proteins have a much higher copy number.

Universal peptide vaccines - optimal peptide vaccine design based on viral sequence conservation

Rapidly mutating viruses such as the hepatitis C virus (HCV), the human immunodeficiency virus (HIV), or influenza viruses (Flu) call for highly effective universal peptide vaccines, i.e. vaccines that do not only yield broad population coverage but also broad coverage of various viral strains. The efficacy of such vaccines is determined by multiple properties of the epitopes they comprise. Beyond the specific properties of each epitope, properties of the corresponding source antigens are of great importance. If a response is mounted against viral proteins with a low copy number within the cell or against proteins expressed very late, this response may fail to induce lysis of the infected cells before budding can take place. We here propose a novel methodology to optimize the epitope composition and assembly in order to induce maximum protection. In order for a peptide vaccine to yield the best possible universal protection, several conditions should be met: (a) an optimal choice of target antigens, (b) an optimal choice of highly conserved epitopes, (c) maximum coverage of the target population, and (d) the proper ordering of the epitopes in the final vaccine to ensure favorable cleavage. We propose a mathematical formalism for epitope selection and ordering that balances the constraints imposed by these different conditions. Focusing on HCV, HIV, and Flu, we show that not all of the conditions can be satisfied for all viruses. Depending on the virus, different constraints are harder to fulfill: for Flu, the conservation constraint is violated first, while for HIV, it is difficult to focus the response at the optimal target antigens. The proposed methodology can be applied to any virus to assess the feasibility of optimally combining the above-mentioned constraints.

Signal peptides and trans-membrane regions are broadly immunogenic and have high CD8+ T cell epitope densities: Implications for vaccine development

Cell mediated immune response has a major role in controlling the elimination of infectious agents. The rational design of sub-unit peptide vaccines against intracellular pathogens or cancer requires the use of antigenic sequence/s that can induce highly potent, long lasting and antigen-specific responses in the majority of the population. A promising peptide selection strategy is the detection of multi-epitope peptide sequences with an ability to bind multiple MHC alleles. While past research sought the best epitopes based on their specific antigenicity, we ask whether specific defined domains have high epitope densities. Signal peptides and trans-membrane domains were found to have exceptionally high epitope densities. The improved MHC binding of these domains relies on their hydrophobic nature and, in signal peptides, also on their specific sequence. The high epitope density of SP was computed using in-silico methods and corroborated by the high percentage of identified SP epitope in the IEDB (immune epitope database). The enhanced immunogenicity of SP was then experimentally confirmed using a panel of nine peptides derived from Mycobacterium tuberculosis (MTb) proteins used in human PBMC proliferation assays and T cell lines functional assays. Our results show the exceptionally high antigen specific response rates and population coverage to SP sequences compared with non-SP peptide antigens derived from the same proteins. The results suggest a novel scheme for the rational design of T cell vaccines using a domain based rather than an epitope based approach.

Immune-induced evolutionary selection focused on a single reading frame in overlapping hepatitis B virus proteins

Viruses employ various means to evade immune detection. Reduction of CD8(+) T cell epitopes is one of the common strategies used for this purpose. Hepatitis B virus (HBV), a member of the Hepadnaviridae family, has four open reading frames, with about 50% overlap between the genes they encode. We computed the CD8(+) T cell epitope density within HBV proteins and the mutations within the epitopes. Our results suggest that HBV accumulates escape mutations that reduce the number of epitopes. These mutations are not equally distributed among genes and reading frames. While the highly expressed core and X proteins are selected to have low epitope density, polymerase, which is expressed at low levels, does not undergo the same selection. In overlapping regions, mutations in one protein-coding sequence also affect the other protein-coding sequence. We show that mutations lead to the removal of epitopes in X and surface proteins even at the expense of the addition of epitopes in polymerase. The total escape mutation rate for overlapping regions is lower than that for nonoverlapping regions. The lower epitope replacement rate for overlapping regions slows the evolutionary escape rate of these regions but leads to the accumulation of mutations more robust in the transfer between hosts, such as mutations preventing proteasomal cleavage into epitopes.

HLA class I binding of HBZ determines outcome in HTLV-1 infection

CD8(+) T cells can exert both protective and harmful effects on the virus-infected host. However, there is no systematic method to identify the attributes of a protective CD8(+) T cell response. Here, we combine theory and experiment to identify and quantify the contribution of all HLA class I alleles to host protection against infection with a given pathogen. In 432 HTLV-1-infected individuals we show that individuals with HLA class I alleles that strongly bind the HTLV-1 protein HBZ had a lower proviral load and were more likely to be asymptomatic. We also show that in general, across all HTLV-1 proteins, CD8(+) T cell effectiveness is strongly determined by protein specificity and produce a ranked list of the proteins targeted by the most effective CD8(+) T cell response through to the least effective CD8(+) T cell response. We conclude that CD8(+) T cells play an important role in the control of HTLV-1 and that CD8(+) cells specific to HBZ, not the immunodominant protein Tax, are the most effective. We suggest that HBZ plays a central role in HTLV-1 persistence. This approach is applicable to all pathogens, even where data are sparse, to identify simultaneously the HLA Class I alleles and the epitopes responsible for a protective CD8(+) T cell response.

Clinical relevance of genetic heterogeneity in HCV

Infection by HCV affects an estimated 170 million people worldwide and it represents one of the major causes of liver transplantation and a heavy burden to healthcare systems. As with many other RNA viruses, HCV is characterized by very high levels of genetic variation, which have been associated to differences in disease progression and efficiency of antiviral treatment. Studies show many contradictory results and little consensus on such associations. Nevertheless, some general guidelines translating research results to clinical practice have been postulated. Here, we review the main research results obtained on HCV variation so far and explore the reasons for their lack of congruence under a population genetics framework. Understanding the factors responsible for the variable dynamics of HCV diversity in human populations and variation within infected individuals is even more necessary in face of the soon-to-arrive new HCV therapies.