Double-stranded RNA as a not-self alarm signal: to evade, most viruses purine-load their RNAs, but some (HTLV-1, Epstein-Barr) pyrimidine-load

Hot springs viruses at Yellowstone National Park have ancient origins and are adapted to thermophilic hosts

Geothermal springs house unicellular red algae in the class Cyanidiophyceae that dominate the microbial biomass at these sites. Little is known about host-virus interactions in these environments. We analyzed the virus community associated with red algal mats in three neighboring habitats (creek, endolithic, soil) at Lemonade Creek, Yellowstone National Park (YNP), USA. We find that despite proximity, each habitat houses a unique collection of viruses, with the giant viruses, Megaviricetes, dominant in all three. The early branching phylogenetic position of genes encoded on metagenome assembled virus genomes (vMAGs) suggests that the YNP lineages are of ancient origin and not due to multiple invasions from mesophilic habitats. The existence of genomic footprints of adaptation to thermophily in the vMAGs is consistent with this idea. The Cyanidiophyceae at geothermal sites originated ca. 1.5 Bya and are therefore relevant to understanding biotic interactions on the early Earth.

The hide-and-seek game of the oncogenic Epstein-Barr virus-encoded EBNA1 protein with the immune system: An RNA G-quadruplex tale

The Epstein-Barr virus (EBV) is the first oncogenic virus described in human. EBV infects more than 90% of the human population worldwide, but most EBV infections are asymptomatic. After the primary infection, the virus persists lifelong in the memory B cells of the infected individuals. Under certain conditions the virus can cause several human cancers, that include lymphoproliferative disorders such as Burkitt and Hodgkin lymphomas and non-lymphoid malignancies such as 100% of nasopharyngeal carcinoma and 10% of gastric cancers. Each year, about 200,000 EBV-related cancers emerge, hence accounting for at least 1% of worldwide cancers. Like all gammaherpesviruses, EBV has evolved a strategy to escape the host immune system. This strategy is mainly based on the tight control of the expression of its Epstein-Barr nuclear antigen-1 (EBNA1) protein, the EBV-encoded genome maintenance protein. Indeed, EBNA1 is essential for viral genome replication and maintenance but, at the same time, is also highly antigenic and T cells raised against EBNA1 exist in infected individuals. For this reason, EBNA1 is considered as the Achilles heel of EBV and the virus has seemingly evolved a strategy that employs the binding of nucleolin, a host cell factor, to RNA G-quadruplex (rG4) within EBNA1 mRNA to limit its expression to the minimal level required for function while minimizing immune recognition. This review recapitulates in a historical way the knowledge accumulated on EBNA1 immune evasion and discusses how this rG4-dependent mechanism can be exploited as an intervention point to unveil EBV-related cancers to the immune system.

Consequences of asymmetric mutational pressure for the dynamic of linear B-cell epitopes repertoire of influenza a virus neuraminidase rearrangement

Full-length nucleotide sequences of avian influenza A virus neuraminidase coding region (20,631 sequences) were analyzed and compared with those isolated from viruses infecting human and swine (63,750 sequences). If in fourfold degenerate sites there is asymmetric A-bias that may be more or less asymmetric depending on the type of neuraminidase and the host, than in twofold degenerate sites from third codon positions there is a strong asymmetric U-bias in coding regions of N4, N5, and N8 isolated from viruses infecting birds, as well as in those of N1 and N2 isolated from viruses infecting human, swine, and birds, while in coding regions of N9 isolated from birds, there is surprisingly strong C-bias, and in sequences of N3, N6, and N7 the usage of C is quite close to the level of U. Revealed stabilization of both U and C in twofold degenerate sites is the evidence of frequent changes in mutational pressure direction. Asymmetric mutational pressure was one of the sources of amino acid replacements that resulted in an equal percentage of sites with appeared and disappeared linear B-cell epitopes in N1, N2, N4, and N5 (33.62-35.33% vs. 32.41-36.45%, respectively), and controlled by the immune pressure it resulted in a stronger tendency to disappear for B-cell epitopes of N3, N6, N7, N8, and N9 of avian viruses (8.74-28.77% vs. 28.96-38.89%). The lack of correlation between nucleotide usages in fourfold and twofold degenerate sites for three nucleotides, except U, is a strong evidence of mutational pressure theory.

Germline mutations directions are different between introns of the same gene: case study of the gene coding for amyloid-beta precursor protein

Amyloid-beta precursor protein (APP) is highly conserved in mammals. This feature allowed us to compare nucleotide usage biases in fourfold degenerated sites along the length of its coding region for 146 species of mammals and birds in search of fragments with significant deviations. Even though cytosine usage has the highest value in fourfold degenerated sites in APP coding region from all tested placental mammals, in contrast to marsupial mammals with the bias toward thymine usage, the most frequent germline and somatic mutations in human APP coding region are C to T and G to A transitions. The same mutational AT-pressure is characteristic for germline mutations in introns of human APP gene. However, surprisingly, there are several exceptional introns with deviations in germline mutations rates. The most of those introns surround exons with exceptional biases in nucleotide usage in fourfold degenerated sites. Existence of such fragments in exons 4 and 5, as well as in exon 14, can be connected with the presence of lncRNA genes in complementary strand of DNA. Exceptional nucleotide usage bias in exons 16 and 17 that contain a sequence encoding amyloid-beta peptides can be explained either by the presence of yet unmapped lncRNA(s), or by the autonomous expression of a short mRNA that encodes just C-terminal part of the APP providing an alternative source of amyloid-beta peptides. This hypothesis is supported by the increased rate of T to C transitions in introns 16-17 and 17-18 of Human APP gene relatively to other introns.

When acting as a reproductive barrier for sympatric speciation, hybrid sterility can only be primary

Animal gametes unite to form a zygote that develops into an adult with gonads that, in turn, produce gametes. Interruption of this germinal cycle by prezygotic or postzygotic reproductive barriers can result in two cycles, each with the potential to evolve into a new species. When the speciation process is complete, members of each species are fully reproductively isolated from those of the other. During speciation a primary barrier may be supported and eventually superceded by a later-appearing secondary barrier. For those holding certain cases of prezygotic isolation to be primary (e.g. elephant cannot copulate with mouse), the onus is to show that they had not been preceded over evolutionary time by periods of postzygotic hybrid inviability (genically determined) or sterility (genically or chromosomally determined). Likewise, the onus is upon those holding cases of hybrid inviability to be primary (e.g. Dobzhansky–Muller epistatic incompatibilities) to show that they had not been preceded by periods, however brief, of hybrid sterility. The latter, when acting as a sympatric barrier causing reproductive isolation, can only be primary. In many cases, hybrid sterility may result from incompatibilities between parental chromosomes that attempt to pair during meiosis in the gonad of their offspring (Winge-Crowther-Bateson incompatibilities). While such incompatibilities have long been observed on a microscopic scale, there is growing evidence for a role of dispersed finer DNA sequence differences (i.e. in base k-mers).

Mutational Pressure in Zika Virus: Local ADAR-Editing Areas Associated with Pauses in Translation and Replication

Zika virus (ZIKV) spread led to the recent medical health emergency of international concern. Understanding the variations in virus system is of utmost need. Using available complete sequences of ZIKV we estimated directions of mutational pressure along the length of consensus sequences of three lineages of the virus. Results showed that guanine usage is growing in ZIKV RNA plus strand due to adenine to guanine transitions, while adenine usage is growing due to cytosine to adenine transversions. Especially high levels of guanine have been found in two-fold degenerated sites of certain areas of RNA plus strand with high amount of secondary structure. The usage of cytosine in two-fold degenerated sites shows direct dependence on the amount of secondary structure in 52% (consensus sequence of East African ZIKV lineage)—32% (consensus sequence of epidemic strains) of the length of RNA minus strand. These facts are the evidences of ADAR-editing of both strands of ZIKV genome during pauses in replication. RNA plus strand can also be edited by ADAR during pauses in translation caused by the appearance of groups of rare codons. According to our results, RNA minus strand of epidemic ZIKV strain has lower number of points in which polymerase can be stalled (allowing ADAR-editing) compared to other strains. The data on preferable directions of mutational pressure in epidemic ZIKV strain is useful for future vaccine development and understanding the evolution of new strains.

Local Mutational Pressures in Genomes of Zaire Ebolavirus and Marburg Virus

Heterogeneities in nucleotide content distribution along the length of Zaire ebolavirus and Marburg virus genomes have been analyzed. Results showed that there is asymmetric mutational A-pressure in the majority of Zaire ebolavirus genes; there is mutational AC-pressure in the coding region of the matrix protein VP40, probably, caused by its high expression at the end of the infection process; there is also AC-pressure in the 3'-part of the nucleoprotein (NP) coding gene associated with low amount of secondary structure formed by the 3'-part of its mRNA; in the middle of the glycoprotein (GP) coding gene that kind of mutational bias is linked with the high amount of secondary structure formed by the corresponding fragment of RNA negative (-) strand; there is relatively symmetric mutational AU-pressure in the polymerase (Pol) coding gene caused by its low expression level. In Marburg virus all genes, including C-rich fragment of GP coding region, demonstrate asymmetric mutational A-bias, while the last gene (Pol) demonstrates more symmetric mutational AU-pressure. The hypothesis of a newly synthesized RNA negative (-) strand shielding by complementary fragments of mRNAs has been described in this work: shielded fragments of RNA negative (-) strand should be better protected from oxidative damage and prone to ADAR-editing.

Lymphocyte repertoire selection and intracellular self/non-self-discrimination: historical overview

Immunological self/non-self-discrimination is conventionally seen as an extracellular event, involving interactions been receptors on T cells pre-educated to discriminate and peptides bound to major histocompatibility complex proteins (pMHCs). Mechanisms by which non-self peptides might first be sorted intracellularly to distinguish them from the vast excess of self-peptides have long been called for. Recent demonstrations of endogenous peptide-specific clustering of pMHCs on membrane rafts are indicative of intracellular enrichment before surface display. The clustering could follow the specific aggregation of a foreign protein that exceeded its solubility limit in the crowded intracellular environment. Predominantly entropy-driven, this homoaggregation would colocalize identical peptides, thus facilitating their collective presentation. Concentrations of self-proteins are fine-tuned over evolutionary time to avoid this. Disparate observations, such as pyrexia and female susceptibility to autoimmune disease, can be explained in terms of the need to cosegregate cognate pMHC complexes internally before extracellular display.

Effects of messenger RNA structure and other translational control mechanisms on major histocompatibility complex-I mediated antigen presentation

Effective T‐cell surveillance of antigen‐presenting cells is dependent on the expression of an array of antigenic peptides bound to major histocompatibility complex (MHC) class I (MHC‐I) or class II (MHC‐II) molecules. Pathogens co‐evolving with their hosts exploit crucial translational regulatory mechanisms in order to evade host immune recognition and thereby sustain their infection. Evasion strategies that downregulate viral protein synthesis and thereby restrict antigen presentation to cytotoxic T‐cells through the endogenous MHC‐I pathway have been implicated in the pathogenesis of viral‐associated malignancies. An understanding of the mechanisms by which messenger RNA (mRNA) structure modulates both viral mRNA translation and the antigen processing machinery to escape immune surveillance, will stimulate the development of alternative therapeutic strategies focused on RNA‐directed drugs designed to enhance immune responses against infected cells. In this review, we discuss regulatory aspects of the MHC‐I pathway and summarize current knowledge of the role attributed by mRNA structure and other translational regulatory mechanisms in immune evasion. In particular we highlight the impact of recently identified G‐quadruplex structures within virally encoded transcripts as unique regulatory signals for translational control and antigen presentation. WIREs RNA 2015, 6:157–171. doi: 10.1002/wrna.1262

This article is categorized under: 1

RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems

2

Translation > Translation Regulation

G-quadruplexes regulate Epstein-Barr virus-encoded nuclear antigen 1 mRNA translation

Viruses that establish latent infections have evolved unique mechanisms to avoid host immune recognition. Maintenance proteins of these viruses regulate their synthesis to levels sufficient for maintaining persistent infection but below threshold levels for host immune detection. The mechanisms governing this finely tuned regulation of viral latency are unknown. Here we show that mRNAs encoding gammaherpesviral maintenance proteins contain within their open reading frames clusters of unusual structural elements, G-quadruplexes, which are responsible for the cis-acting regulation of viral mRNA translation. By studying the Epstein-Barr virus-encoded nuclear antigen 1 (EBNA1) mRNA, we demonstrate that destabilization of G-quadruplexes using antisense oligonucleotides increases EBNA1 mRNA translation. In contrast, pretreatment with a G-quadruplex-stabilizing small molecule, pyridostatin, decreases EBNA1 synthesis, highlighting the importance of G-quadruplexes within virally encoded transcripts as unique regulatory signals for translational control and immune evasion. Furthermore, these findings suggest alternative therapeutic strategies focused on targeting RNA structure within viral ORFs.

A yeast-based assay identifies drugs that interfere with immune evasion of the Epstein-Barr virus

Epstein-Barr virus (EBV) is tightly associated with certain human cancers, but there is as yet no specific treatment against EBV-related diseases. The EBV-encoded EBNA1 protein is essential to maintain viral episomes and for viral persistence. As such, EBNA1 is expressed in all EBV-infected cells, and is highly antigenic. All infected individuals, including individuals with cancer, have CD8(+) T cells directed towards EBNA1 epitopes, yet the immune system fails to detect and destroy cells harboring the virus. EBV immune evasion depends on the capacity of the Gly-Ala repeat (GAr) domain of EBNA1 to inhibit the translation of its own mRNA in cis, thereby limiting the production of EBNA1-derived antigenic peptides presented by the major histocompatibility complex (MHC) class I pathway. Here we establish a yeast-based assay for monitoring GAr-dependent inhibition of translation. Using this assay we identify doxorubicin (DXR) as a compound that specifically interferes with the GAr effect on translation in yeast. DXR targets the topoisomerase-II-DNA complexes and thereby causes genomic damage. We show, however, that the genotoxic effect of DXR and various analogs thereof is uncoupled from the effect on GAr-mediated translation control. This is further supported by the observation that etoposide and teniposide, representing another class of topoisomerase-II-DNA targeting drugs, have no effect on GAr-mediated translation control. DXR and active analogs stimulate, in a GAr-dependent manner, EBNA1 expression in mammalian cells and overcome GAr-dependent restriction of MHC class I antigen presentation. These results validate our approach as an effective high-throughput screening assay to identify drugs that interfere with EBV immune evasion and, thus, constitute candidates for treating EBV-related diseases, in particular EBV-associated cancers.

Implications of HIV RNA structure for recombination, speciation, and the neutralism-selectionism controversy

The conflict between the needs to encode both a protein (impaired by non-synonymous mutation), and nucleic acid structure (impaired by synonymous or non-synonymous mutation), can sometimes be resolved in favour of the nucleic acid because its structure is critical for a selectively advantageous genome-wide activity--recombination. However, above a sequence difference threshold, recombination is impaired. It may then be advantageous for new species to arise. Building on the work of Grantham and others critical of the neutralist viewpoint, heuristic support for this hypothesis emerged from studies of the base composition and structure of retroviral genomes. The extreme enrichment in the purine A of the RNA of human immunodeficiency virus (HIV-1), parallels the mild purine-loading of the RNAs of most organisms, for which there is an adaptive explanation--immune evasion. However, human T cell leukaemia virus (HTLV-1), with the potential to invade the same host cell, shows extreme enrichment in the pyrimidine C. Assuming the low GC% HIV and the high GC% HTLV-1 to share a common ancestor, it was postulated that differences in GC% had arisen to prevent homologous recombination between these emerging lentiviral species. Sympatrically isolated by this intracellular reproductive barrier, prototypic HIV-1 seized the AU-rich (low GC%) high ground (thus committing to purine A rather than purine G). Prototypic HTLV-1 forwent this advantage and evolved an independent evolutionary strategy--similar to that of the GC%-rich Epstein-Barr virus--profound latency maintained by transcription of one purine-rich mRNA. The evidence supporting these interpretations is reviewed.

Low complexity regions (LCRs) contribute to the hypervariability of the HIV-1 gp120 protein

Low complexity regions (LCRs) are sequences of nucleic acids or proteins defined by a compositional bias. Their occurrence has been confirmed in sequences of the three cellular lineages (Bacteria, Archaea and Eucarya), and has also been reported in viral genomes. We present here the results of a detailed computer analysis of the LCRs present in the HIV-1 glycoprotein 120 (gp120) encoded by the viral gene env. The analysis was performed using a sample of 3637 Env polyprotein sequences derived from 4117 completely sequenced and translated HIV-1 genomes available in public databases as of December 2012. We have identified 1229 LCRs located in four different regions of the gp120 protein that correspond to four of the five regions that have been identified as hypervariable (V1, V2, V4 and V5). The remaining 29 LCRs are found in the signal peptide and in the conserved regions C2, C3, C4 and C5. No LCR has been identified in the hypervariable region V3. The LCRs detected in the V1, V2, V4, and V5 hypervariable regions exhibit a high Asn content in their amino acid composition, which very likely correspond to glycosylation sites, which may contribute to the retroviral ability to avoid the immune system. In sharp contrast with what is observed in gp120 proteins lacking LCRs, the glycosylation sites present in LCRs tend to be clustered towards the center of the region forming well-defined islands. The results presented here suggest that LCRs represent a hitherto undescribed source of genomic variability in lentivirus, and that these repeats may represent an important source of antigenic variation in HIV-1 populations. The results reported here may exemplify the evolutionary processes that may have increased the size of primitive cellular RNA genomes and the role of LCRs as a source of raw material during the processes of evolutionary acquisition of new functions.

Messenger RNA sequence rather than protein sequence determines the level of self-synthesis and antigen presentation of the EBV-encoded antigen, EBNA1

Unique purine-rich mRNA sequences embedded in the coding sequences of a distinct group of gammaherpesvirus maintenance proteins underlie the ability of the latently infected cell to minimize immune recognition. The Epstein-Barr virus nuclear antigen, EBNA1, a well characterized lymphocryptovirus maintenance protein has been shown to inhibit in cis antigen presentation, due in part to a large internal repeat domain encoding glycine and alanine residues (GAr) encoded by a purine-rich mRNA sequence. Recent studies have suggested that it is the purine-rich mRNA sequence of this repeat region rather than the encoded GAr polypeptide that directly inhibits EBNA1 self-synthesis and contributes to immune evasion. To test this hypothesis, we generated a series of EBNA1 internal repeat frameshift constructs and assessed their effects on cis-translation and endogenous antigen presentation. Diverse peptide sequences resulting from alternative repeat reading frames did not alleviate the translational inhibition characteristic of EBNA1 self-synthesis or the ensuing reduced surface presentation of EBNA1-specific peptide-MHC class I complexes. Human cells expressing the EBNA1 frameshift variants were also poorly recognized by antigen-specific T-cells. Furthermore, a comparative analysis of the mRNA sequences of the corresponding repeat regions of different viral maintenance homologues highlights the high degree of identity between the nucleotide sequences despite very little homology in the encoded amino acid sequences. Based on these combined observations, we propose that the cis-translational inhibitory effect of the EBNA1 internal repeat sequence operates mechanistically at the nucleotide level, potentially through RNA secondary structural elements, and is unlikely to be mediated through the GAr polypeptide. The demonstration that the EBNA1 repeat mRNA sequence and not the encoded protein sequence underlies immune evasion in this class of virus suggests a novel approach to therapeutic development through the use of anti-sense strategies or small molecules targeting EBNA1 mRNA structure.

Viruses establishing persistent latent infections have evolved various mechanisms to avoid immune surveillance. The Epstein-Barr virus-encoded nuclear antigen, EBNA1, expressed in all EBV-associated malignancies, modulates its own protein levels at quantities sufficient to maintain viral infection but low enough so as to minimize an immune response by the infected host cell. This evasion mechanism is regulated through an internal purine-rich mRNA repeat sequence encoding glycine and alanine residues. In this study we assess the impact of the repeat's nucleotide versus peptide sequence on inhibiting EBNA1 self-synthesis and antigen presentation. We demonstrate that altered peptide sequences resulting from frameshift mutations within the repeat do not alleviate the immune-evasive function of EBNA1, suggesting that the repetitive purine-rich mRNA sequence itself is responsible for inhibiting EBNA1 synthesis and subsequent poor immunogenicity. Our comparative analysis of the mRNA sequences of the corresponding repeat regions of different gammaherpesvirus maintenance homologues to EBNA1 highlights the high degree of identity between the nucleotide sequences despite very little homology in the encoded amino acid sequences. These studies demonstrate the importance of gammaherpesvirus purine-rich mRNA repeat sequences on antigenic epitope generation and evasion from T-cell mediated immune control, suggesting novel approaches to prevention and treatment of latent infection by this class of virus.

CHROMOSOMES AS INTERDEPENDENT ACCOUNTING UNITS: THE ASSIGNED ORIENTATION OFC. ELEGANSCHROMOSOMES MINIMIZES THE TOTAL W-BASE CHARGAFF DIFFERENCE

DNAs of individual chromosomes violate, albeit perhaps by only one in a thousand bases, Chargaff's second parity rule, which is that Chargaff's first parity rule for duplex DNA (A = T, G = C) applies, to a close approximation, to single stranded DNA. If the "top" strand of one chromosome has A > T and the "top" strand of another has T > A, can they complement to approach even parity (A = T)? Assignment of orientation to the six chromosomes of Caenorhabditis elegans is said to have been arbitrary and, of 26(= 64) possible combinations of top (T) and bottom (B) strands, the GenBank orientation (designated "TTTTTT") is but one. Yet, for the W bases (A and T) the chromosomes in the GenBank orientation complement to reduce the Chargaff difference (A–T) to only 200 bases (i.e. only one in 323,658 bases does not have a potential Watson-Crick pairing partner). This suggests that the assignment was not arbitrary. However, the GenBank orientation for the S bases (G and C) allows an approach to even parity less well than many other orientations, the best of which is BBBBTT (indicating a disparity between the GenBank orientations of the first four autosomes and those of chromosomes V and X). Although only the euchromatic regions of Drosophila melanogaster chromosomes have been sequenced, there are orientations that allow an approach to even parity. We conclude that, with respect to their Chargaff differences, the chromosomes of C. elegans have the potential to engage in interdependent base accounting. Since this might also apply to D. melanogaster, even when heterochromatin-associated DNA rich in tandem repeats (microsatellite DNA) is excluded, then heterochromatic DNA might not normally participate in the hypothetical accounting process.

POSITIVE DARWINIAN SELECTION: DOES THE COMPARATIVE METHOD RULE?

To detect positive Darwinian selection it is thought essential to compare two sequences. Despite its defects, "the comparative method rules." However, genes evolving rapidly under positive selection conflict more with internal forces (the genome phenotype) than genes evolving slowly under negative selection. In particular, there is conflict with stem-loop potential. The conflict between protein-encoding potential (primary information) and stem-loop potential (secondary information) permits detection of positive selection in a single sequence. The degree to which secondary information is compromised provides a measure of the speed of transmission of primary information. Thus, the sovereignty of the comparative method is challenged not only by its own defects, but also by the availability of a single-sequence method. However, while of limited utility for positive selection, the comparative method casts new light on Darwin's great question — the origin of species. Comparison of rates of synonymous and non-synonymous mutation suggests that branching into new species begins with synonymous mutations.

Percent of highly immunogenic amino acid residues forming B-cell epitopes is higher in homologous proteins encoded by GC-rich genes

We analyzed the dependence of the percent of highly immunogenic amino acid residues included in B-cell epitopes of homologous proteins on the GC-content (G+C) of genes coding for them in twenty-seven lineages of proteins (and subsequent genes), which belong to seven Varicello and five Simplex viruses. We found out that proteins encoded by genes of a high GC-content usually contain more targets for humoral immune response than their homologs encoded by GC-poor genes. This tendency is characteristic not only to the lineages of glycoproteins, which are the main targets for humoral immune response against Simplex and Varicello viruses, but also to the lineages of capsid proteins and even "housekeeping" enzymes. The percent of amino acids included in linear B-cell epitopes has been predicted for 324 proteins by BepiPred algorithm (www.cbs.dtu.dk/services/BepiPred), the percent of highly immunogenic amino acids included in discontinuous B-cell epitopes and the percent of exposed amino acid residues have been predicted by Epitopia algorithm (http://epitopia.tau.ac.il/). Immunological consequences of the directional mutational GC-pressure are mostly due to the decrease in the total usage of highly hydrophobic amino acids and due to the increase in proline and glycine levels of usage in proteins. The weaker the negative selection on amino acid substitutions caused by symmetric mutational pressure, the higher the slope of direct dependence of the percent of highly immunogenic amino acids included in B-cell epitopes on G+C.

The level of cytosine is usually much higher than the level of guanine in two-fold degenerated sites from third codon positions of genes from Simplex- and Varicelloviruses with G+C higher than 50%

We studied usage of cytosine and guanine in 914 genes from completely sequenced genomes of five Simplex- and seven Varicelloviruses. In genes with total GC-content higher than 50% usage of cytosine is usually higher than usage of guanine (an average difference for genes with G+C higher than 70% reaches 4.0%). This difference is caused mostly by the elevated usage of cytosine in two-fold degenerated sites situated in third codon positions relatively to the usage of guanine in two-fold degenerated sites situated in third codon positions (an average difference for genes with G+C higher than 70% is equal to 28.2%). The usage of amino acids that are encoded by codons containing cytosine in two-fold degenerated sites situated in third codon positions (AA2TC) is much higher than the usage of amino acids encoded by codons containing guanine in two-fold degenerated sites situated in third codon positions (AA2AG). The usage of AA2AG declines much more steeply with the growth of GC-content than the usage of AA2TC. This effect is the consequence of the nature of genetic code and of the negative selection. In GC-rich genes the usage of cytosine in four-fold degenerated sites is only a little (but significantly) higher than the usage of guanine (in genes with G+C higher than 70% an average difference is equal to 4.3%). This difference may be caused by transcription-associated mutational pressure.

The synthesis of truncated polypeptides for immune surveillance and viral evasion

Background

Cytotoxic T cells detect intracellular pathogens by surveying peptide loaded MHC class I molecules (pMHC I) on the cell surface. Effective immune surveillance also requires infected cells to present pMHC I promptly before viral progeny can escape. Rapid pMHC I presentation apparently occurs because infected cells can synthesize and present peptides from antigenic precursors called defective ribosomal products (DRiPs). The molecular characteristics of DRiPs are not known.

Methodology/Principal Findings

Here, using a novel method for detecting antigenic precursors and proteolytic intermediates, we tracked the synthesis and processing of Epstein-Barr Virus encoded nuclear antigen 1 (EBNA1). We find that ribosomes initiated translation appropriately, but rapidly produced DRiPs representing ∼120 amino acid truncated EBNA1 polypeptides by premature termination. Moreover, specific sequences in EBNA1 mRNA strongly inhibited the generation of truncated DRiPs and pMHC I presentation.

Significance

Our results reveal the first characterization of virus DRiPs as truncated translation products. Furthermore, production of EBNA1-derived DRiPs is down-regulated in cells, possibly limiting the antigenicity of EBNA1.

Immune evasion by gammaherpesvirus genome maintenance proteins

Viruses that establish lifelong latent infections must ensure that the viral genome is maintained within the latently infected cell throughout the life of the host, yet at the same time must also be capable of avoiding elimination by the immune surveillance system. Gammaherpesviruses, which include the human viruses Epstein-Barr virus and Kaposi's sarcoma-associated herpesvirus, establish latent infections in lymphocytes. Infection of this dynamic host-cell population requires that the viruses have appropriate strategies for enabling the viral genome to persist while these cells go through rounds of mitosis, but at the same time must avoid detection by host CD8(+) cytotoxic T lymphocytes (CTLs). The majority of gammaherpesviruses studied have been found to encode a specific protein that is critical for maintenance of the viral genome within latently infected cells. This protein is termed the genome maintenance protein (GMP). Due to its vital role in long-term latency, this offers the immune system a crucial target for detection and elimination of virus-infected cells. GMPs from different gammaherpesviruses have evolved related strategies that allow the protein to be present within latently infected cells, but to remain effectively hidden from circulating CD8(+) CTLs. In this review, I will summarize the role of the GMPs and highlight the available data describing the immune-evasion properties of these proteins.

X chromosome reactivation perturbs intracellular self/not-self discrimination

New reports indicate a chromosomal rather than hormonal basis for the susceptibility of females to autoimmune disease. It is held that if females reactivate an inactivated X chromosome, there will be overexpression of certain X-located genes affecting immune function. Hence, normal mechanisms of self/not-self discrimination might be impaired resulting in immune reaction to self antigens. However, the data are also consistent with the long-held view that the demands of intracellular self/not-self discrimination have driven the evolution of X-chromosome dosage compensation. It was proposed that, whether cells are in male or female bodies, concentrations of proteins are fine-tuned up to their aggregation thresholds. A disruption of this equilibrium, by agents originating either externally (for example, virus) or internally (for example, reactivated X chromosome), generates homoaggregates that trigger responses against the respective not-self or self antigens. Thus, female susceptibility to autoimmune disease may not be because certain immune system genes happen to be X-located, but because self/not-self discrimination was the raison d'être for X-chromosome dosage compensation in the first place.

Epstein-Barr virus evasion of CD8(+) and CD4(+) T cell immunity via concerted actions of multiple gene products

Upon primary infection, EBV establishes a latent infection in B cells, characterized by maintenance of the viral genome in the absence of viral replication. The Epstein-Barr Nuclear Antigen 1 (EBNA1) plays a crucial role in maintenance of the viral DNA episome and is consistently expressed in all EBV-associated malignancies. Compared to other EBV latent gene products, EBNA1 is poorly recognized by CD8(+) T lymphocytes. Recent studies are discussed that shed new light on the mechanisms that underlie this unusual lack of CD8(+) T cell activation. Whereas the latent phase is characterized by the expression of a limited subset of viral gene products, the full repertoire of over 80 EBV lytic gene products is expressed during the replicative phase. Despite this abundance of potential T cell antigens, which indeed give rise to a strong response of CD4(+) and CD8(+) T lymphocytes, the virus can replicate successfully. Evidence is accumulating that this paradoxical situation is the result of actions of multiple viral gene products, inhibiting discrete stages of the MHC class I and class II antigen presentation pathways. Immediately after initiation of the lytic cycle, BNLF2a prevents peptide-loading of MHC class I molecules through inhibition of the Transporter associated with Antigen Processing, TAP. This will reduce presentation of viral antigens by the large ER-resident pool of MHC class I molecules. Synthesis of new MHC class I molecules is blocked by BGLF5. Viral-IL10 causes a reduction in mRNA levels of TAP1 and bli/LMP2, a subunit of the immunoproteasome. MHC class I molecules present at the cell surface are downregulated by BILF1. Also the antigen presenting capacity of MHC class II molecules is severely compromised by multiple EBV lytic gene products, including gp42/gH/gL, BGLF5, and vIL-10. In this review, we discuss how concerted actions of these EBV lytic proteins result in highly effective interference with CD8(+) and CD4(+) T cell surveillance, thereby providing the virus with a window for undisturbed generation of viral progeny.

Reduction in RNA levels rather than retardation of translation is responsible for the inhibition of major histocompatibility complex class I antigen presentation by the glutamic acid-rich repeat of herpesvirus saimiri open reading frame 73

Herpesvirus saimiri (HVS) establishes a persistent infection in squirrel monkeys by maintaining its episome within T lymphocytes. The product of open reading frame 73 (ORF73) plays a key role in episomal maintenance and is the functional homologue of Epstein-Barr virus EBNA1 and Kaposi's sarcoma-associated herpesvirus LANA1 proteins. There is little sequence homology among these proteins, although all contain a central domain of repeating amino acids. The repeat domains of EBNA1 and LANA1 enhance the stability of these proteins and cause a retardation in self-protein synthesis, leading to poor recognition by CD8(+) cytotoxic T lymphocytes (CTL). The HVS ORF73 repeat domain is composed of a glutamic acid and glycine repeat linked to a glutamic acid and alanine repeat (EG-EA repeat). Here we show that the EG-EA repeat similarly causes a reduction in the recognition of ORF73 by CD8(+) CTL. However, deletion of the EG-EA repeat from HVS ORF73 had no affect on the stability of the protein or its rate of translation. In contrast, the presence of the EG-EA repeat was found to decrease the steady-state levels of ORF73 mRNA. The inhibitory properties of the EG-EA repeat were maintained when transferred to a heterologous protein, and manipulation of the repeat revealed that the motif EEAEEAEEE was sufficient to cause a reduction in recognition of ORF73 by CD8(+) CTL. Thus, the EG-EA repeat of HVS ORF73 plays a role in immune evasion but utilizes a mechanism distinct from that of the EBNA1 and LANA1 repeats.

Regulation of protein translation through mRNA structure influences MHC class I loading and T cell recognition

Many viruses avoid immune surveillance during latent infection through reduction in the synthesis of virally encoded proteins. Although antigen presentation critically depends on the level of viral protein synthesis, the precise mechanism used to regulate the generation of antigenic peptide precursors remains elusive. Here, we demonstrate that a purine overloaded virally encoded mRNA lacking secondary structure significantly impacts the efficiency of protein translation and prevents endogenous antigen presentation. Reducing this purine bias through the generation of constructs expressing codon-modified sequences, while maintaining the encoded protein sequence, increased the stem-loop structure of the corresponding mRNA and dramatically enhanced self-synthesis of the viral protein. As a consequence, a higher number of HLA-peptide complexes were detected on the surface of cells expressing this viral protein. Furthermore, these cells were more efficiently recognized by virus-specific T cells compared with those expressing the same antigen expressed by a purine-biased mRNA. These findings delineate a mechanism by which viruses regulate self-synthesis of proteins and offer an effective strategy to evade CD8(+) T cell-mediated immune regulation.

Prokaryotes that grow optimally in acid have purine-poor codons in long open reading frames

In nucleic acids the N-glycosyl bonds between purines and their ribose sugar moities are broken under acid conditions. If one strand of a duplex DNA segment were more vulnerable to mutation than the other, then the archaeon Picrophilus torridus, with an optimum growth pH near zero, could have adapted by decreasing the purine content of that strand. Yet, P. torridus has an optimum growth temperature near 60 degrees C, and thermophiles prefer purine-rich codons. We found that, as in other thermophiles, high growth temperature correlates with the use of purine-rich codons. The extra purines are often in third, non-amino acid determining, codon positions. However, as in other acidophiles, as open reading frame lengths increase, there is increased use of purine-poor codons, particularly those without purines in second, amino acid-determining, codon positions. Thus, P. torridus can be seen as adapting (a) to temperature by increasing its purines in all open reading frames without greatly impacting protein amino acid compositions, and (b) to pH by decreasing purines in longer open reading frames, thereby potentially impacting protein amino acid compositions. It is proposed that longer open reading frames, being larger mutational targets, have become less vulnerable to depurination by virtue of pyrimidine for purine substitutions.

Amino acids as placeholders: base-composition pressures on protein length in malaria parasites and prokaryotes

BACKGROUND

The composition and sequence of amino acids in a protein may serve the underlying needs of the nucleic acids that encode the protein (the genome phenotype). In extreme form, amino acids become mere placeholders inserted between functional segments or domains, and--apart from increasing protein length--playing no role in the specific function or structure of a protein (the conventional phenotype).

METHODS

We studied the genomes of two malarial parasites and 521 prokaryotes (144 complete) that differ widely in GC% and optimum growth temperature, comparing the base compositions of the protein coding regions and corresponding lengths (kilobases).

RESULTS

Malarial parasites show distinctive responses to base-compositional pressures that increase as protein lengths increase. A low-GC% species (Plasmodium falciparum) is likely to have more placeholder amino acids than an intermediate-GC% species (P. vivax), so that homologous proteins are longer. In prokaryotes, GC% is generally greater and AG% is generally less in open reading frames (ORFs) encoding long proteins. The increased GC% in long ORFs increases as species' GC% increases, and decreases as species' AG% increases. In low- and intermediate-GC% prokaryotic species, increases in ORF GC% as encoded proteins increase in length are largely accounted for by the base compositions of first and second (amino acid-determining) codon positions. In high-GC% prokaryotic species, first and third (non-amino acid-determining) codon positions play this role.

CONCLUSION

In low- and intermediate-GC% prokaryotes, placeholder amino acids are likely to be well defined, corresponding to codons enriched in G and/or C at first and second positions. In high-GC% prokaryotes, placeholder amino acids are likely to be less well defined. Increases in ORF GC% as encoded proteins increase in length are greater in mesophiles than in thermophiles, which are constrained from increasing protein lengths in response to base-composition pressures.

Optimum growth temperature and the base composition of open reading frames in prokaryotes

The purine-loading index (PLI) is the difference between the numbers of purines (A+G) and pyrimidines (T+C) per kilobase of single-stranded nucleic acid. By purine-loading their mRNAs organisms may minimize unnecessary RNA-RNA interactions and prevent inadvertent formation of "self" double-stranded RNA. Since RNA-RNA interactions have a strong entropy-driven component, this need to minimize should increase as temperature increases. Consistent with this, we report for 550 prokaryotic species that optimum growth temperature is related to the average PLI of open reading frames. With increasing temperature prokaryotes tend to acquire base A and lose base C, while keeping bases T and G relatively constant. Accordingly, while the PLI increases, the (G+C)% decreases. The previously observed positive correlation between (G+C)% and optimum growth temperature, which applies to RNA species whose structure is of major importance for their function (ribosomal and transfer RNAs) does not apply to mRNAs, and hence is unlikely to apply generally to genomic DNA.

Low-complexity segments in Plasmodium falciparum proteins are primarily nucleic acid level adaptations

Protein segments that contain few of the possible 20 amino acids, sometimes in tandem repeat arrays, are referred to as containing "simple" or "low-complexity" sequence. Many Plasmodium falciparum proteins are longer than their homologs in other species by virtue of their content of such low-complexity segments that have no known function; these are interspersed among segments of higher complexity to which function can often be ascribed. If there is low complexity at the protein level, there is likely to be low complexity at the corresponding nucleic acid level (departure from equifrequency of the four bases). Thus, low complexity may have been selected primarily at the nucleic acid level and low complexity at the protein level may be secondary. In this case, the amino acid composition of low-complexity segments should be more reflective than that of high complexity segments on forces operating at the nucleic acid level, which include GC-pressure and AG-pressure. Consistent with this, for amino acid determining first and second codon positions, open reading frames containing low-complexity segments show increased contributions to downward GC-pressure (revealed as decreased percentage of G+C) and to upward AG-pressure (revealed as increased percentage A+G). When not countermanded by high contributions to AG-pressure, low-complexity segments can contribute to base order-dependent fold potential; in this respect, they resemble introns. Thus, in P. falciparum, low-complexity segments appear as adaptations primarily serving nucleic acid level functions.

Immunity as a function of the unicellular state: implications of emerging genomic data

Instead of being greeted as supporting the growing corpus of immunological theory, recent advances in the bioinformatic analysis of genomes have often surprised the discoverers and failed to attract the attention of immunologists. In fact, the view that multicellular immune systems are adaptations of already highly evolved unicellular immune systems that are capable of self/not-self discrimination can assist our comprehension of phenomena, such as 'junk' DNA, genetic polymorphism and the ubiquity of repetitive elements. For instance, the 'hidden transcriptome', revealed by run-on transcription of genes or repetitive elements, contains a diverse repertoire of RNA 'immune receptors' with the potential to form double-stranded RNA with viral RNA 'antigens', thus triggering intracellular alarms.

Selective pressures that decrease synonymous mutations in Plasmodium falciparum

Rich and Ayala propose that the zero rate of non-amino-acid-changing (synonymous) mutations in some proteins of Plasmodium falciparum reflects a recent population bottleneck. Alternatively, Arnot and Saul propose sequence conservation in response to selective pressures other than the pressure to encode protein. Among these are fold pressure and purine-loading pressure. Genomes adapt to these by acquisition of introns and/or low-complexity (simple-sequence) segments in proteins. Adaptive explanations include facilitation of intragenic recombination (and hence diversification of the encoded protein) by DNA stem-loop secondary structures.

Adaptive value of polymorphism in intracellular self/not-self discrimination?

A microbial pathogen species can adapt to its host species to the extent that members of the host species are uniform. Loss of this uniformity would make it difficult for a pathogen species to transfer, from one member of the host species to another, what it had "learned" through selection of its members with advantageous mutations. The existence of major histocompatibility complex (MHC) polymorphism indicates that non-uniformity within a species is an effective host defence strategy. By virtue of this molecular discontinuity among its members the host species can "present a moving target" to the pathogen. Many proteins other than MHC proteins show polymorphism - a phenomenon which has suggested that mutations in regions of protein molecules which do not affect overt function are neutral. However, in the context of the author's differential aggregation theory of intracellular self/not-self discrimination as previously applied to the problem of the antigenicity of cancer cells, such polymorphism should serve for the recruitment of subsets of self-antigens into the antigenic repertoire of an infected cell. These would act as "intracellular antibodies" by virtue of their weak, but specific, aggregation with pathogen proteins. Peptides from the self-antigens, as well as (or instead of) those from the antigens of the pathogen, would then serve as targets for attack by cytotoxic T cells. Thus, polymorphism of intracellular proteins should be of adaptive value, serving to amplify and individualize the immune response to intracellular pathogens.