Phase-dependent immune evasion of herpesviruses

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.

Engineering Tolerance toward Allogeneic CAR-T Cells by Regulation of MHC Surface Expression with Human Herpes Virus-8 Proteins

Allogeneic, off-the-shelf (OTS) chimeric antigen receptor (CAR) cell therapies have the potential to reduce manufacturing costs and variability while providing broader accessibility to cancer patients and those with other diseases. However, host-versus-graft reactivity can limit the durability and efficacy of OTS cell therapies requiring new strategies to evade adaptive and innate-immune responses. Human herpes virus-8 (HHV8) maintains infection, in part, by evading host T and natural killer (NK) cell attack. The viral K3 gene encodes a membrane-tethered E3 ubiquitin ligase that discretely targets major histocompatibility complex (MHC) class I components, whereas K5 encodes a similar E3 ligase with broader specificity, including MHC-II and the MHC-like MHC class I polypeptide-related sequence A (MIC-A)- and sequence B (MIC-B)-activating ligands of NK cells. We created γ-retroviruses encoding K3 and/or K5 transgenes that efficiently transduce primary human T cells. Expression of K3 or K5 resulted in dramatic downregulation of MHC-IA (human leukocyte antigen [HLA]-A, -B, and -C) and MHC class II (HLA-DR) cell-surface expression. K3 expression was sufficient for T cells to resist exogenously loaded peptide-MHC-specific cytotoxicity, as well as recognition in one-way allogeneic mixed lymphocyte reactions. Further, in immunodeficient mice engrafted with allogeneic T cells, K3-transduced T cells selectively expanded in vivo. Ectopic K5 expression in MHC class I^-, MIC-A⁺/B⁺ K562 cells also reduced targeting by primary NK cells. Coexpression of K3 in prostate stem cell antigen (PSCA)-directed, inducible MyD88/CD40 (iMC)-enhanced CAR-T cells did not impact cytotoxicity, T cell growth, or cytokine production against HPAC pancreatic tumor target cells, whereas K5-expressing cells showed a modest reduction in interleukin (IL)-2 production without effect on cytotoxicity. Together, these results support application of these E3 ligases to advance development of OTS CAR-T cell products.

Battle between Host Immune Cellular Responses and HCMV Immune Evasion

Human cytomegalovirus (HCMV) is ubiquitously prevalent. HCMV infection is typically asymptomatic and controlled by the immune system in healthy individuals, yet HCMV can be severely pathogenic for the fetus during pregnancy and in immunocompromised persons, such as transplant recipients or HIV infected patients. HCMV has co-evolved with the hosts, developed strategies to hide from immune effector cells and to successfully survive in the human organism. One strategy for evading or delaying the immune response is maintenance of the viral genome to establish the phase of latency. Furthermore, HCMV immune evasion involves the downregulation of human leukocyte antigens (HLA)-Ia molecules to hide infected cells from T-cell recognition. HCMV expresses several proteins that are described for downregulation of the HLA class I pathway via various mechanisms. Here, we review the wide range of immune evasion mechanisms of HCMV. Understanding the mechanisms of HCMV immune evasion will contribute to the development of new customized therapeutic strategies against the virus.

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.

The limited benefit of "product-by-process" claim

A product-by-process claim was initially developed when the invention could not otherwise be adequately defined. In the US, a patent applicant can use a product-by-process claim for a new product, no matter whether the product can also be defined by using structure features. However, the applicant mainly bears the burden to establish the novelty between the claimed product and the prior art product. Moreover, in terms of infringement, the very recent CAFC cases indicate that the product must be made by a process recited in the claim to infringe a product-by-process claim. Thus, an assignee almost can not take advantage from product-by-process claims during the patent enforcement by comparing with process claims. From the points of the burden during patent prosecution and confined scope in determining infringement, the benefit of using a product-by-process claim would be very limited.

Revelation of Influencing Factors in Overall Codon Usage Bias of Equine Influenza Viruses

Equine influenza viruses (EIVs) of H3N8 subtype are culprits of severe acute respiratory infections in horses, and are still responsible for significant outbreaks worldwide. Adaptability of influenza viruses to a particular host is significantly influenced by their codon usage preference, due to an absolute dependence on the host cellular machinery for their replication. In the present study, we analyzed genome-wide codon usage patterns in 92 EIV strains, including both H3N8 and H7N7 subtypes by computing several codon usage indices and applying multivariate statistical methods. Relative synonymous codon usage (RSCU) analysis disclosed bias of preferred synonymous codons towards A/U-ended codons. The overall codon usage bias in EIVs was slightly lower, and mainly affected by the nucleotide compositional constraints as inferred from the RSCU and effective number of codon (ENc) analysis. Our data suggested that codon usage pattern in EIVs is governed by the interplay of mutation pressure, natural selection from its hosts and undefined factors. The H7N7 subtype was found less fit to its host (horse) in comparison to H3N8, by possessing higher codon bias, lower mutation pressure and much less adaptation to tRNA pool of equine cells. To the best of our knowledge, this is the first report describing the codon usage analysis of the complete genomes of EIVs. The outcome of our study is likely to enhance our understanding of factors involved in viral adaptation, evolution, and fitness towards their hosts.

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.

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.

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

Kaposi's sarcoma-associated herpesvirus K3 and K5 ubiquitin E3 ligases have stage-specific immune evasion roles during lytic replication

The downregulation of immune synapse components such as major histocompatibility complex class I (MHC-I) and ICAM-1 is a common viral immune evasion strategy that protects infected cells from targeted elimination by cytolytic effector functions of the immune system. Kaposi's sarcoma-associated herpesvirus (KSHV) encodes two membrane-bound ubiquitin E3 ligases, called K3 and K5, which share the ability to induce internalization and degradation of MHC-I molecules. Although individual functions of K3 and K5 outside the viral genome are well characterized, their roles during the KSHV life cycle are still unclear. In this study, we individually introduced the amino acid-coding sequences of K3 or K5 into a ΔK3 ΔK5 recombinant virus, at either original or interchanged genomic positions. Recombinants harboring coding sequences within the K5 locus showed higher K3 and K5 protein expression levels and more rapid surface receptor downregulation than cognate recombinants in which coding sequences were introduced into the K3 locus. To identify infected cells undergoing K3-mediated downregulation of MHC-I, we employed a novel reporter virus, called red-green-blue-BAC16 (RGB-BAC16), which was engineered to harbor three fluorescent protein expression cassettes: EF1α-monomeric red fluorescent protein 1 (mRFP1), polyadenylated nuclear RNA promoter (pPAN)-enhanced green fluorescent protein (EGFP), and pK8.1-monomeric blue fluorescent protein (tagBFP), marking latent, immediate early, and late viral gene expression, respectively. Analysis of RGB-derived K3 and K5 deletion mutants showed that while the K5-mediated downregulation of MHC-I was concomitant with pPAN induction, the reduction of MHC-I surface expression by K3 was evident in cells that were enriched for pPAN-driven EGFP(high) and pK8.1-driven blue fluorescent protein-positive (BFP(+)) populations. These data support the notion that immunoreceptor downregulation occurs by a sequential process wherein K5 is critical during the immediately early phase and K3 plays a significant role during later stages.

IMPORTANCE

Although the roles of K3 and K5 outside the viral genome are well characterized, the function of these proteins in the context of the KSHV life cycle has remained unclear, particularly in the case of K3. This study examined the relative contributions of K3 and K5 to the downregulation of MHC-I during the lytic replication of KSHV. We show that while K5 acts immediately upon entry into the lytic phase, K3-mediated downregulation of MHC-I was evident during later stages of lytic replication. The identification of distinctly timed K3 and K5 activities significantly advances our understanding of KSHV-mediated immune evasion. Crucial to this study was the development of a novel recombinant KSHV, called RGB-BAC16, which facilitated the delineation of stage-specific phenotypes.

In vitro transfection of bone marrow-derived dendritic cells with TATp-liposomes

Dendritic cells (DC) are antigen-presenting cells uniquely capable of priming naïve T cells and cross-presenting antigens, and they determine the type of immune response elicited against an antigen. TAT peptide (TATp), is an amphipathic, arginine-rich, cationic peptide that promotes penetration and translocation of various molecules and nanoparticles into cells. TATp-liposomes (TATp-L) used for DC transfection were prepared using TATp derivatized with a lipid-terminated polymer capable of anchoring in the liposomal membrane. Here, we show that the addition of TATp to DNA-loaded liposomes increased the uptake of DNA in DC. DNA-loaded TATp-L increased the in vitro transfection efficiency in DC cultures as evidenced by a higher expression of the enhanced green fluorescent protein and bovine herpes virus type 1 glycoprotein D (gD). The de novo synthesized gD protein was immunologically stimulating when transfections were performed with TATp-L, as indicated by the secretion of interleukin 6.

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.

Evasion of adaptive and innate immune response mechanisms by γ-herpesviruses

γ-Herpesviral immune evasion mechanisms are optimized to support the acute, lytic and the longterm, latent phase of infection. During acute infection, specific immune modulatory proteins limit, but also exploit, the antiviral activities of cell intrinsic innate immune responses as well as those of innate and adaptive immune cells. During latent infection, a restricted gene expression program limits immune targeting and cis-acting mechanisms to reduce the antigen presentation as well as antigenicity of latency-associated proteins. Here, we will review recent progress in our understanding of γ-herpesviral immune evasion strategies.

Predictor for the effect of amino acid composition on CD4+ T cell epitopes preprocessing

Predictive tools for all levels of CD8+ T cell epitopes processing have reached a maturation level. Good prediction algorithms have been developed for proteasomal cleavage, TAP and MHC class I peptide binding. The same cannot be said of CD4+ T cell epitopes. While multiple algorithms of varying accuracy have been proposed for MHC class II peptide binding, the preprocessing of CD4+ T cell epitopes is still lacking a good prediction algorithm. CD4+ T cell epitopes generation includes several stages, not all which are well-defined. We here group these stages to produce a generic preprocessing stage predictor for the cleavage processes preceding the presentation of epitopes to CD4+ T cell. The predictor is learnt using a combination of in vitro cleavage experiments and observed naturally processed MHC class II binding peptides. The properties of the predictor highlight the effect of different factors on CD4+ T cell epitopes preprocessing. The most important factor emerging from the predictor is the secondary structure of the cleaved region in the protein. The effect of the secondary structure is expected since CD4+ T cell epitopes are not denatured before cleavage. A website developed based on this predictor is available at: http://peptibase.cs.biu.ac.il/PepCleave_cd4/.

Professional antigen presenting cells in human herpesvirus 8 infection

Professional antigen presenting cells (APC), i.e., dendritic cells (DC), monocytes/macrophages, and B lymphocytes, are critically important in the recognition of an invading pathogen and presentation of antigens to the T cell-mediated arm of immunity. Human herpesvirus 8 (HHV-8) is one of the few human viruses that primarily targets these APC for infection, altering their cytokine profiles, manipulating their surface expression of MHC molecules, and altering their ability to activate HHV-8-specific T cells. This could be why T cell responses to HHV-8 antigens are not very robust. Of these APC, only B cells support complete, lytic HHV-8 infection. However, both complete and abortive virus replication cycles in APC could directly affect viral pathogenesis and progression to Kaposi's sarcoma (KS) and HHV-8-associated B cell cancers. In this review, we discuss the effects of HHV-8 infection on professional APC and their relationship to the development of KS and B cell lymphomas.

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.

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.

Loss of HIV-1-derived cytotoxic T lymphocyte epitopes restricted by protective HLA-B alleles during the HIV-1 epidemic

OBJECTIVE AND DESIGN

HIV-1 is known to adapt to the human immune system, leading to accumulation of escape mutations during the course of infection within an individual. Cross-sectional studies have shown an inverse correlation between the prevalence of human leukocyte antigen (HLA) alleles in a population and the number of cytotoxic T lymphocyte (CTL) escape mutations in epitopes restricted by those HLA alleles. Recently, it was demonstrated that at a population level HIV-1 is adapting to the humoral immune response, which is reflected in an increase in resistance to neutralizing antibodies over time. Here we investigated whether adaptations to cellular immunity have also accumulated during the epidemic.

METHODS

We compared the number of CTL epitopes in HIV-1 strains isolated from individuals who seroconverted at the beginning of the HIV-1 epidemic and from individuals who seroconverted in recent calendar time.

RESULTS

The number of CTL epitopes in HIV-1 variants restricted by the most common HLA alleles in the population did not change significantly during the epidemic. In contrast, we found a significant loss of CTL epitopes restricted by HLA-B alleles associated with a low relative hazard of HIV-1 disease progression during the epidemic. Such a loss was not observed for CTL epitopes restricted by HLA-A alleles.

CONCLUSION

Despite the large degree of HLA polymorphism, HIV-1 has accumulated adaptations to CTL responses within 20 years of the epidemic. The fact that such adaptations are driven by the HLA-B molecules that provide best protection against HIV-1 disease progression has important implications for our understanding of HIV evolution.

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.

Mapping the landscape of host-pathogen coevolution: HLA class I binding and its relationship with evolutionary conservation in human and viral proteins

The high diversity of HLA binding preferences has been driven by the sequence diversity of short segments of relevant pathogenic proteins presented by HLA molecules to the immune system. To identify possible commonalities in HLA binding preferences, we quantify these using a novel measure termed "targeting efficiency," which captures the correlation between HLA-peptide binding affinities and the conservation of the targeted proteomic regions. Analysis of targeting efficiencies for 95 HLA class I alleles over thousands of human proteins and 52 human viruses indicates that HLA molecules preferentially target conserved regions in these proteomes, although the arboviral Flaviviridae are a notable exception where nonconserved regions are preferentially targeted by most alleles. HLA-A alleles and several HLA-B alleles that have maintained close sequence identity with chimpanzee homologues target conserved human proteins and DNA viruses such as Herpesviridae and Adenoviridae most efficiently, while all HLA-B alleles studied efficiently target RNA viruses. These patterns of host and pathogen specialization are both consistent with coevolutionary selection and functionally relevant in specific cases; for example, preferential HLA targeting of conserved proteomic regions is associated with improved outcomes in HIV infection and with protection against dengue hemorrhagic fever. Efficiency analysis provides a novel perspective on the coevolutionary relationship between HLA class I molecular diversity, self-derived peptides that shape T-cell immunity through ontogeny, and the broad range of viruses that subsequently engage with the adaptive immune response.

Thyroid hormone controls the gene expression of HSV-1 LAT and ICP0 in neuronal cells

Various factors/pathways including hormonal regulation have been suggested to control HSV-1 latency and reactivation. Our computer analysis identified a DNA repeat containing thyroid hormone response elements (TRE) in the regulatory region of HSV-1 LAT. Thyroid hormone (T₃) exerts its function via its receptor (TR), a transcriptional factor. Present study investigated the roles of TR and T₃ on HSV-1 gene expression using cultured neuoroblastoma cell lines. We demonstrated that liganded TR activated LAT transcription but repressed ICP0 transcription in the presence of LAT TRE. The chromatin immunoprecipitation (ChIP) assays showed that TRs were recruited to LAT TREs independently of T₃ and hyperacetylated H4 was associated with promoters that were transcriptionally active. In addition, ChIP results showed that a chromatin insulator protein CTCF was enriched at the LAT TREs in the presence of TR and T₃. In addition, chromatin remodeling factor BRG1 complex is found to participate in the T₃/TR-mediated LAT activation since overexpression of BRG1 enhanced the LAT transcription and the dominant negative mutant K785R abolished the activation. This is the first report revealing that TR exerted epigenetic regulation on HSV-1 ICP0 expression in neuronal cells and could have a role in the complex processes of HSV-1 latency/reactivation.

Human tumor-associated viruses and new insights into the molecular mechanisms of cancer

The study of acute-transforming retroviruses and their oncogenes and of the multiple mechanisms deployed by DNA viruses to circumvent the growth-suppressive and proapoptotic function of tumor suppressor genes has provided the foundation of our current understanding of cancer biology. Unlike acute-transforming animal viruses, however, human tumor-associated viruses lead to malignancies with a prolonged latency and in conjunction with other environmental and host-related cooperating events. The relevance of viral infection to human cancer development has often been debated. We now know that at least six human viruses, Epstein-Barr virus (EBV), hepatitis B virus (HBV), hepatitis C virus (HCV), human papilloma virus (HPV), human T-cell lymphotropic virus (HTLV-1) and Kaposi's associated sarcoma virus (KSHV) contribute to 10-15% of the cancers worldwide. Hence, the opportunity exists to fight cancer at the global scale by preventing the spread of these viruses, by the development and distribution of effective and safe antiviral vaccines, and by identifying their oncogenic mechanism. Here, we discuss the molecular events underlying the neoplastic potential of the human tumor-associated viruses, with emphasis on the enigmatic KSHV and its numerous virally hijacked proangiogenic, immune-evasive and tumor-promoting genes. The emerging information may facilitate the development of new molecular-targeted approaches to prevent and treat virally associated human malignancies.

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

Motivation: Viruses employ various means to evade immune detection. One common evasion strategy is the removal of CD8+cytotoxic T-lymphocyte epitopes. We here use a combination of multiple bioinformatic tools and large amount of genomic data to compute the epitope repertoire presented by over 1300 viruses in many HLA alleles. We define the ‘Size of Immune Repertoire score’, which represents the ratio between the epitope density within a protein and the expected density. This score is used to study viral immune evasion.

Results: We show that viral proteins in general have a higher epitope density than human proteins. This difference is due to a good fit of the human MHC molecules to the typical amino-acid usage of viruses. Among different viruses, viruses infecting humans present less epitopes than non-human viruses. This selection is not at the amino-acid usage level, but through the removal of specific epitopes. Within a single virus, not all proteins express the same epitopes density. Proteins expressed early in the viral life cycle have a lower epitope density than late proteins. Such a difference is not observed in non-human viruses. The removal of early epitopes and the targeting of the cellular immune response to late viral proteins, allow the virus a time interval to propagate before its host cells are destroyed by T cells.

Contact:louzouy@math.biu.ac.il

The HIV hide and seek game: an immunogenomic analysis of the HIV epitope repertoire

BACKGROUND

Viruses employ various means to evade immune detection. One common evasion strategy is the removal of CD8 cytotoxic T-lymphocyte (CTL) epitopes.

METHOD

Here, we use bioinformatic tools to compute the HIV CTL epitope repertoire presented by over 8000 HIV sequences in multiple Human Leukocyte Antigen alleles. We define the 'Size of Immune Repertoire' (SIR) score, which represents the ratio between the number of the predicted epitopes within a protein and their expected number within a scrambled version of the same protein.

RESULTS

We show that HIV proteins present less epitopes than expected and that the number of epitopes gradually decreases from SIV to recent HIV sequences. The decrease of the SIR score of HIV is accompanied by a high frequency of replacement mutations within epitopes. The SIR score of the different HIV proteins is not uniform. The regulatory proteins, Tat and Rev, expressed early during cellular infection have a low SIR score, whereas virion-associated genes that are expressed later, such as Env, Pol and Gag, have a higher SIR score. Actually, the SIR score of Gag keeps increasing over time.

CONCLUSION

We hypothesize that our results reflect an HIV immune evasion strategy. This involves the targeting of the CTL immune response to viral structural and enzyme proteins, allowing the virus a time interval to propagate before its host cells are destroyed by CTLs. An efficient anti-HIV CTL response against HIV should thus also target the regulatory genes that HIV seeks to hide from the immune system.

Human self-protein CD8+ T-cell epitopes are both positively and negatively selected

The cellular immune system recognizes self-epitopes in the context of MHC-I molecules. The immunological general view presumes that these self-epitopes are just a background, both positively and negatively selecting T cells. We here estimate the number of epitopes in each human protein for many frequent HLA alleles, and a score representing over or under presentation of epitopes on these proteins. We further show that there is a clear selection for the presentation of specific self-protein types. Proteins presenting many epitopes include, for example, autoimmune regulator (AIRE) upregulated tissue-specific antigens, immune system receptors and proteins with a high expression level. On the other hand, proteins that may be considered less "useful" for the immune system, such as low expression level proteins, are under-presented. We combine our epitope estimate with single nucleotide polymorphism (SNP) measures to show that this selection can be directly observed through the fraction of non-synonymous SNP (replacement fraction), which is significantly higher inside epitopes than outside.

Patterns of evolution and host gene mimicry in influenza and other RNA viruses

It is well known that the dinucleotide CpG is under-represented in the genomic DNA of many vertebrates. This is commonly thought to be due to the methylation of cytosine residues in this dinucleotide and the corresponding high rate of deamination of 5-methycytosine, which lowers the frequency of this dinucleotide in DNA. Surprisingly, many single-stranded RNA viruses that replicate in these vertebrate hosts also have a very low presence of CpG dinucleotides in their genomes. Viruses are obligate intracellular parasites and the evolution of a virus is inexorably linked to the nature and fate of its host. One therefore expects that virus and host genomes should have common features. In this work, we compare evolutionary patterns in the genomes of ssRNA viruses and their hosts. In particular, we have analyzed dinucleotide patterns and found that the same patterns are pervasively over- or under-represented in many RNA viruses and their hosts suggesting that many RNA viruses evolve by mimicking some of the features of their host's genes (DNA) and likely also their corresponding mRNAs. When a virus crosses a species barrier into a different host, the pressure to replicate, survive and adapt, leaves a footprint in dinucleotide frequencies. For instance, since human genes seem to be under higher pressure to eliminate CpG dinucleotide motifs than avian genes, this pressure might be reflected in the genomes of human viruses (DNA and RNA viruses) when compared to those of the same viruses replicating in avian hosts. To test this idea we have analyzed the evolution of the influenza virus since 1918. We find that the influenza A virus, which originated from an avian reservoir and has been replicating in humans over many generations, evolves in a direction strongly selected to reduce the frequency of CpG dinucleotides in its genome. Consistent with this observation, we find that the influenza B virus, which has spent much more time in the human population, has adapted to its human host and exhibits an extremely low CpG dinucleotide content. We believe that these observations directly show that the evolution of RNA viral genomes can be shaped by pressures observed in the host genome. As a possible explanation, we suggest that the strong selection pressures acting on these RNA viruses are most likely related to the innate immune response and to nucleotide motifs in the host DNA and RNAs.

Viruses are obligate intracellular parasites that use different strategies to sequester host cell machinery and avoid the host immune system. In this paper we explore the genomes of viruses that encode their genetic information in single-stranded RNA, a different material than the one used by their hosts (double-stranded DNA). It is interesting to observe that these viruses share some of the host's characteristics. For instance, one of the most underrepresented motifs in the DNA of vertebrates is the dinucleotide CpG. This is commonly thought to be due to methylation and deamination of cytosine residues in this dinucleotide. Surprisingly, the same CpG suppression is observed in vertebrate RNA viruses but not in RNA phages. We show that RNA viruses present similar dinucleotide pressures as their host genes. We find that the influenza A virus, which originated from an avian reservoir and replicated in humans over many generations, evolves to reduce the frequency of CpG dinucleotides mimicking the human genes. Influenza B, which has been in humans longer, exhibits an extremely low CpG dinucleotide content. These observations suggest that the evolution of RNA viruses is shaped by pressures observed in the host genome.