Identification of potential conserved RNA secondary structure throughout influenza A coding regions

A Conserved Secondary Structural Element in the Coding Region of the Influenza A Virus Nucleoprotein (NP) mRNA Is Important for the Regulation of Viral Proliferation

Influenza A virus is a threat to humans due to seasonal epidemics and infrequent, but dangerous, pandemics that lead to widespread infection and death. Eight segments of RNA constitute the genome of this virus and they encode greater than eight proteins via alternative splicing of coding (+)RNAs generated from the genomic (-)RNA template strand. RNA is essential in its life cycle. A bioinformatics analysis of segment 5, which encodes nucleoprotein, revealed a conserved structural motif in the (+)RNA. The secondary structure proposed by energy minimization and comparative analysis agrees with structure predicted based on experimental data using a 121 nucleotide in vitro RNA construct comprising an influenza A virus consensus sequence and also an entire segment 5 (+)RNA (strain A/VietNam/1203/2004 (H5N1)). The conserved motif consists of three hairpins with one being especially thermodynamically stable. The biological importance of this conserved secondary structure is supported in experiments using antisense oligonucleotides in cell line, which found that disruption of this motif led to inhibition of viral fitness. These results suggest that this conserved motif in the segment 5 (+)RNA might be a candidate for oligonucleotide-based antiviral therapy.

Generation of Live Attenuated Influenza Virus by Using Codon Usage Bias

Seasonal influenza epidemics and occasional pandemics threaten public health worldwide. New alternative strategies for generating recombinant viruses with vaccine potential are needed. Interestingly, influenza viruses circulating in different hosts have been found to have distinct codon usage patterns, which may reflect host adaptation. We therefore hypothesized that it is possible to make a human seasonal influenza virus that is specifically attenuated in human cells but not in eggs by converting its codon usage so that it is similar to that observed from avian influenza viruses. This approach might help to generate human live attenuated viruses without affecting their yield in eggs. To test this hypothesis, over 300 silent mutations were introduced into the genome of a seasonal H1N1 influenza virus. The resultant mutant was significantly attenuated in mammalian cells and mice, yet it grew well in embryonated eggs. A single dose of intranasal vaccination induced potent innate, humoral, and cellular immune responses, and the mutant could protect mice against homologous and heterologous viral challenges. The attenuated mutant could also be used as a vaccine master donor strain by introducing hemagglutinin and neuraminidase genes derived from other strains. Thus, our approach is a successful strategy to generate attenuated viruses for future application as vaccines.

IMPORTANCE

Vaccination has been one of the best protective measures in combating influenza virus infection. Current licensed influenza vaccines and their production have various limitations. Our virus attenuation strategy makes use of the codon usage biases of human and avian influenza viruses to generate a human-derived influenza virus that is attenuated in mammalian hosts. This method, however, does not affect virus replication in eggs. This makes the resultant mutants highly compatible with existing egg-based vaccine production pipelines. The viral proteins generated from the codon bias mutants are identical to the wild-type viral proteins. In addition, our massive genome-wide mutational approach further minimizes the concern over reverse mutations. The potential use of this kind of codon bias mutant as a master donor strain to generate other live attenuated viruses is also demonstrated. These findings put forward a promising live attenuated influenza vaccine generation strategy to control influenza.

In silico discovery and modeling of non-coding RNA structure in viruses

This review covers several computational methods for discovering structured non-coding RNAs in viruses and modeling their putative secondary structures. Here we will use examples from two target viruses to highlight these approaches: influenza A virus-a relatively small, segmented RNA virus; and Epstein-Barr virus-a relatively large DNA virus with a complex transcriptome. Each system has unique challenges to overcome and unique characteristics to exploit. From these particular cases, generically useful approaches can be derived for the study of additional viral targets.

The Influenza A PB1-F2 and N40 Start Codons Are Contained within an RNA Pseudoknot

Influenza A is a negative-sense RNA virus with an eight-segment genome. Some segments encode more than one polypeptide product, but how the virus accesses alternate internal open reading frames (ORFs) is not completely understood. In segment 2, ribosomal scanning produces two internal ORFs, PB1-F2 and N40. Here, chemical mapping reveals a Mg(2+)-dependent pseudoknot structure that includes the PB1-F2 and N40 start codons. The results suggest that interactions of the ribosome with the pseudoknot may affect the level of translation for PB1-F2 and N40.

Structural features of a 3' splice site in influenza a

Influenza A is an RNA virus with a genome of eight negative sense segments. Segment 7 mRNA contains a 3′ splice site for alternative splicing to encode the essential M2 protein. On the basis of sequence alignment and chemical mapping experiments, the secondary structure surrounding the 3′ splice site has an internal loop, adenine bulge, and hairpin loop when it is in the hairpin conformation that exposes the 3′ splice site. We report structural features of a three-dimensional model of the hairpin derived from nuclear magnetic resonance spectra and simulated annealing with restrained molecular dynamics. Additional insight was provided by modeling based on ¹H chemical shifts. The internal loop containing the 3′ splice site has a dynamic guanosine and a stable imino (cis Watson–Crick/Watson–Crick) GA pair. The adenine bulge also appears to be dynamic with the A either stacked in the stem or forming a base triple with a Watson–Crick GC pair. The hairpin loop is a GAAA tetraloop closed by an AC pair.

RNase L targets distinct sites in influenza A virus RNAs

Influenza A virus (IAV) infections are influenced by type 1 interferon-mediated antiviral defenses and by viral countermeasures to these defenses. When IAV NS1 protein is disabled, RNase L restricts virus replication; however, the RNAs targeted for cleavage by RNase L under these conditions have not been defined. In this study, we used deep-sequencing methods to identify RNase L cleavage sites within host and viral RNAs from IAV PR8ΔNS1-infected A549 cells. Short hairpin RNA knockdown of RNase L allowed us to distinguish between RNase L-dependent and RNase L-independent cleavage sites. RNase L-dependent cleavage sites were evident at discrete locations in IAV RNA segments (both positive and negative strands). Cleavage in PB2, PB1, and PA genomic RNAs suggests that viral RNPs are susceptible to cleavage by RNase L. Prominent amounts of cleavage mapped to specific regions within IAV RNAs, including some areas of increased synonymous-site conservation. Among cellular RNAs, RNase L-dependent cleavage was most frequent at precise locations in rRNAs. Our data show that RNase L targets specific sites in both host and viral RNAs to restrict influenza virus replication when NS1 protein is disabled.

IMPORTANCE

RNase L is a critical component of interferon-regulated and double-stranded-RNA-activated antiviral host responses. We sought to determine how RNase L exerts its antiviral activity during influenza virus infection. We enhanced the antiviral activity of RNase L by disabling a viral protein, NS1, that inhibits the activation of RNase L. Then, using deep-sequencing methods, we identified the host and viral RNAs targeted by RNase L. We found that RNase L cleaved viral RNAs and rRNAs at very precise locations. The direct cleavage of IAV RNAs by RNase L highlights an intimate battle between viral RNAs and an antiviral endonuclease.

Microarrays for identifying binding sites and probing structure of RNAs

Oligonucleotide microarrays are widely used in various biological studies. In this review, application of oligonucleotide microarrays for identifying binding sites and probing structure of RNAs is described. Deep sequencing allows fast determination of DNA and RNA sequence. High-throughput methods for determination of secondary structures of RNAs have also been developed. Those methods, however, do not reveal binding sites for oligonucleotides. In contrast, microarrays directly determine binding sites while also providing structural insights. Microarray mapping can be used over a wide range of experimental conditions, including temperature, pH, various cations at different concentrations and the presence of other molecules. Moreover, it is possible to make universal microarrays suitable for investigations of many different RNAs, and readout of results is rapid. Thus, microarrays are used to provide insight into oligonucleotide sequences potentially able to interfere with biological function. Better understanding of structure-function relationships of RNA can be facilitated by using microarrays to find RNA regions capable to bind oligonucleotides. That information is extremely important to design optimal sequences for antisense oligonucleotides and siRNA because both bind to single-stranded regions of target RNAs.

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

Influenza virus non-structural protein NS1: interferon antagonism and beyond

Most viruses express one or several proteins that counter the antiviral defences of the host cell. This is the task of non-structural protein NS1 in influenza viruses. Absent in the viral particle, but highly expressed in the infected cell, NS1 dramatically inhibits cellular gene expression and prevents the activation of key players in the IFN system. In addition, NS1 selectively enhances the translation of viral mRNAs and may regulate the synthesis of viral RNAs. Our knowledge of the virus and of NS1 has increased dramatically during the last 15 years. The atomic structure of NS1 has been determined, many cellular partners have been identified and its multiple activities have been studied in depth. This review presents our current knowledge, and attempts to establish relationships between the RNA sequence, the structure of the protein, its ligands, its activities and the pathogenicity of the virus. A better understanding of NS1 could help in elaborating novel antiviral strategies, based on either live vaccines with altered NS1 or on small-compound inhibitors of NS1.

Secondary structure of a conserved domain in an intron of influenza A M1 mRNA

Influenza A virus utilizes RNA throughout infection. Little is known, however, about the roles of RNA structures. A previous bioinformatics survey predicted multiple regions of influenza A virus that are likely to generate evolutionarily conserved and stable RNA structures. One predicted conserved structure is in the pre-mRNA coding for essential proteins, M1 and M2. This structure starts 79 nucleotides downstream of the M2 mRNA 5' splice site. Here, a combination of biochemical structural mapping, mutagenesis, and NMR confirms the predicted three-way multibranch structure of this RNA. Imino proton NMR spectra reveal no change in secondary structure when 80 mM KCl is supplemented with 4 mM MgCl2. Optical melting curves in 1 M NaCl and in 100 mM KCl with 10 mM MgCl2 are very similar, with melting temperatures ∼14 °C higher than that for 100 mM KCl alone. These results provide a firm basis for designing experiments and potential therapeutics to test for function in cell culture.

RNA structural constraints in the evolution of the influenza A virus genome NP segment

Conserved RNA secondary structures were predicted in the nucleoprotein (NP) segment of the influenza A virus genome using comparative sequence and structure analysis. A number of structural elements exhibiting nucleotide covariations were identified over the whole segment length, including protein-coding regions. Calculations of mutual information values at the paired nucleotide positions demonstrate that these structures impose considerable constraints on the virus genome evolution. Functional importance of a pseudoknot structure, predicted in the NP packaging signal region, was confirmed by plaque assays of the mutant viruses with disrupted structure and those with restored folding using compensatory substitutions. Possible functions of the conserved RNA folding patterns in the influenza A virus genome are discussed.

Influenza viruses and mRNA splicing: doing more with less

During their nuclear replication stage, influenza viruses hijack the host splicing machinery to process some of their RNA segments, the M and NS segments. In this review, we provide an overview of the current knowledge gathered on this interplay between influenza viruses and the cellular spliceosome, with a particular focus on influenza A viruses (IAV). These viruses have developed accurate regulation mechanisms to reassign the host spliceosome to alter host cellular expression and enable an optimal expression of specific spliced viral products throughout infection. Moreover, IAV segments undergoing splicing display high levels of similarity with human consensus splice sites and their viral transcripts show noteworthy secondary structures. Sequence alignments and consensus analyses, along with recently published studies, suggest both conservation and evolution of viral splice site sequences and structure for improved adaptation to the host. Altogether, these results emphasize the ability of IAV to be well adapted to the host's splicing machinery, and further investigations may contribute to a better understanding of splicing regulation with regard to viral replication, host range, and pathogenesis.

Identification of conserved RNA secondary structures at influenza B and C splice sites reveals similarities and differences between influenza A, B, and C

BACKGROUND

Influenza B and C are single-stranded RNA viruses that cause yearly epidemics and infections. Knowledge of RNA secondary structure generated by influenza B and C will be helpful in further understanding the role of RNA structure in the progression of influenza infection.

FINDINGS

All available protein-coding sequences for influenza B and C were analyzed for regions with high potential for functional RNA secondary structure. On the basis of conserved RNA secondary structure with predicted high thermodynamic stability, putative structures were identified that contain splice sites in segment 8 of influenza B and segments 6 and 7 of influenza C. The sequence in segment 6 also contains three unused AUG start codon sites that are sequestered within a hairpin structure.

CONCLUSIONS

When added to previous studies on influenza A, the results suggest that influenza splicing may share common structural strategies for regulation of splicing. In particular, influenza 3' splice sites are predicted to form secondary structures that can switch conformation to regulate splicing. Thus, these RNA structures present attractive targets for therapeutics aimed at targeting one or the other conformation.

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.

Unmasking the information encoded as structural motifs of viral RNA genomes: a potential antiviral target

RNA viruses show enormous capacity to evolve and adapt to new cellular and molecular contexts, a consequence of mutations arising from errors made by viral RNA-dependent RNA polymerase during replication. Sequence variation must occur, however, without compromising functions essential for the completion of the viral cycle. RNA viruses are safeguarded in this respect by their genome carrying conserved information that does not code only for proteins but also for the formation of structurally conserved RNA domains that directly perform these critical functions. Functional RNA domains can interact with other regions of the viral genome and/or proteins to direct viral translation, replication and encapsidation. They are therefore potential targets for novel therapeutic strategies. This review summarises our knowledge of the functional RNA domains of human RNA viruses and examines the achievements made in the design of antiviral compounds that interfere with their folding and therefore their function.

Genome-wide analysis of selective constraints on high stability regions of mRNA reveals multiple compensatory mutations in Escherichia coli

Message RNA (mRNA) carries a large number of local secondary structures, with structural stability to participate in the regulations of gene expression. A worthy question is how the local structural stability is maintained under the constraint that multiple selective pressures are imposed on mRNA local regions. Here, we performed the first genome-wide study of natural selection operating on high structural stability regions (HSRs) of mRNAs in Escherichia coli. We found that HSR tends to adjust the folded conformation to reduce the harm of mutations, showing a high level of mutational robustness. Moreover, guanine preference in HSR was observed, supporting the hypothesis that the selective constraint for high structural stability may partly account for the high percentage of G content in Escherichia coli genome. Notably, we found a substantially reduced synonymous substitution rate in HSRs compared with that in their adjacent regions. Surprisingly and interestingly, the non-key sites in HSRs, which have slight effect on structural stability, have synonymous substitution rate equivalent to background regions. To explain this result, we identified compensatory mutations in HSRs based on structural stability, and found that a considerable number of synonymous mutations occur to restore the structural stability decreased heavily by the mutations on key sites. Overall, these results suggest a significant role of local structural stability as a selective force operating on mRNA, which furthers our understanding of the constraints imposed on protein-coding RNAs.

Secondary structure of a conserved domain in the intron of influenza A NS1 mRNA

Influenza A virus is a segmented single-stranded (−)RNA virus that causes substantial annual morbidity and mortality. The transcriptome of influenza A is predicted to have extensive RNA secondary structure. The smallest genome segment, segment 8, encodes two proteins, NS1 and NEP, via alternative splicing. A conserved RNA domain in the intron of segment 8 may be important for regulating production of NS1. Two different multi-branch loop structures have been proposed for this region. A combination of in vitro chemical mapping and isoenergetic microarray techniques demonstrate that the consensus sequence for this region folds into a hairpin conformation. These results provide an alternative folding for this region and a foundation for designing experiments to probe its functional role in the influenza life cycle.

Influenza B virus has global ordered RNA structure in (+) and (-) strands but relatively less stable predicted RNA folding free energy than allowed by the encoded protein sequence

Background

Influenza A virus contributes to seasonal epidemics and pandemics and contains Global Ordered RNA structure (GORS) in the nucleoprotein (NP), non-structural (NS), PB2, and M segments. A related virus, influenza B, is also a major annual public health threat, but unlike influenza A is very selective to human hosts. This study extends the search for GORS to influenza B.

Findings

A survey of all available influenza B sequences reveals GORS in the (+) and (−)RNAs of the NP, NS, PB2, and PB1 gene segments. The results are similar to influenza A, except GORS is observed for the M1 segment of influenza A but not for PB1. In general, the folding free energies of human-specific influenza B RNA segments are less stable than allowable by the encoded amino acid sequence. This is consistent with findings in influenza A, where human-specific influenza RNA folds are less stable than avian and swine strains.

Conclusions

These results reveal fundamental molecular similarities and differences between Influenza A and B and suggest a rational basis for choosing segments to target with therapeutics and for viral attenuation for live vaccines by altering RNA folding stability.

Number variation of high stability regions is correlated with gene functions

Various regulatory elements in messenger RNAs (mRNAs) carrying the secondary structure play important roles in a wide range of expression processes. Numerous recent works have focused on the discovery of these functional elements that contain the conserved mRNA structures. However, to date, regions with high structural stability have been largely overlooked. In this study, we defined high stability regions (HSRs) in the coding sequences (CDSs) in bacteria based on the normalized folding free energy. We found that CDSs had high number of HSRs, and these HSRs showed high structural context robustness compared with random sequences, indicating a direct selective constraint imposed on HSRs. A reduced ribosome speed was detected near the start position of HSR, implying a possibility that HSR acted as obstacle to drive translational pausing that coordinated protein synthesis. Interestingly, we found that genes with high HSR density were enriched in the processes of translation, protein folding, and cell division. In addition, essential genes exhibited higher HSR density than nonessential genes. Overall, our study presented the previously unappreciated correlation between the number variation of HSRs and cellular processes.

Genome-wide analyses of Epstein-Barr virus reveal conserved RNA structures and a novel stable intronic sequence RNA

BACKGROUND

Epstein-Barr virus (EBV) is a human herpesvirus implicated in cancer and autoimmune disorders. Little is known concerning the roles of RNA structure in this important human pathogen. This study provides the first comprehensive genome-wide survey of RNA and RNA structure in EBV.

RESULTS

Novel EBV RNAs and RNA structures were identified by computational modeling and RNA-Seq analyses of EBV. Scans of the genomic sequences of four EBV strains (EBV-1, EBV-2, GD1, and GD2) and of the closely related Macacine herpesvirus 4 using the RNAz program discovered 265 regions with high probability of forming conserved RNA structures. Secondary structure models are proposed for these regions based on a combination of free energy minimization and comparative sequence analysis. The analysis of RNA-Seq data uncovered the first observation of a stable intronic sequence RNA (sisRNA) in EBV. The abundance of this sisRNA rivals that of the well-known and highly expressed EBV-encoded non-coding RNAs (EBERs).

CONCLUSION

This work identifies regions of the EBV genome likely to generate functional RNAs and RNA structures, provides structural models for these regions, and discusses potential functions suggested by the modeled structures. Enhanced understanding of the EBV transcriptome will guide future experimental analyses of the discovered RNAs and RNA structures.

Influenza A virus PB1-F2 protein expression is regulated in a strain-specific manner by sequences located downstream of the PB1-F2 initiation codon

Translation of influenza A virus PB1-F2 occurs in a second open reading frame (ORF) of the PB1 gene segment. PB1-F2 has been implicated in regulation of polymerase activity, immunopathology, susceptibility to secondary bacterial infection, and induction of apoptosis. Experimental evidence of PB1-F2 molecular function during infection has been collected primarily from human and avian viral isolates. As the 2009 H1N1 (H1N1pdm09) strain highlighted, some swine-derived influenza viruses have the capacity to infect human hosts and emerge as a pandemic. Understanding the impact that virulence factors from swine isolates have on both human and swine health could aid in early identification of viruses with pandemic potential. Studies examining PB1-F2 from swine isolates have focused primarily on H1N1pdm09, which does not encode PB1-F2 but was engineered to carry a full-length PB1-F2 ORF to assess the impact on viral replication and pathogenicity. However, experimental evidence of PB1-F2 protein expression from swine lineage viruses has not been demonstrated. Here, we reveal that during infection, PB1-F2 expression levels are substantially different in swine and human influenza viruses. We provide evidence that PB1-F2 expression is regulated at the translational level, with very low levels of PB1-F2 expression from swine lineage viruses relative to a human isolate PB1-F2. Translational regulation of PB1-F2 expression was partially mapped to two independent regions within the PB1 mRNA, located downstream of the PB1-F2 start site. Our data suggest that carrying a full-length PB1-F2 ORF may not be predictive of PB1-F2 expression in infected cells for all influenza A viruses.

Conservation of mRNA secondary structures may filter out mutations in Escherichia coli evolution

Recent reports indicate that mutations in viral genomes tend to preserve RNA secondary structure, and those mutations that disrupt secondary structural elements may reduce gene expression levels, thereby serving as a functional knockout. In this article, we explore the conservation of secondary structures of mRNA coding regions, a previously unknown factor in bacterial evolution, by comparing the structural consequences of mutations in essential and nonessential Escherichia coli genes accumulated over 40 000 generations in the course of the ‘long-term evolution experiment’. We monitored the extent to which mutations influence minimum free energy (MFE) values, assuming that a substantial change in MFE is indicative of structural perturbation. Our principal finding is that purifying selection tends to eliminate those mutations in essential genes that lead to greater changes of MFE values and, therefore, may be more disruptive for the corresponding mRNA secondary structures. This effect implies that synonymous mutations disrupting mRNA secondary structures may directly affect the fitness of the organism. These results demonstrate that the need to maintain intact mRNA structures imposes additional evolutionary constraints on bacterial genomes, which go beyond preservation of structure and function of the encoded proteins.

The effect of gene overlapping on the rate of RNA virus evolution

Gene overlapping is widely employed by RNA viruses to generate genetic novelty while retaining a small genome size. However, gene overlapping also increases the deleterious effect of mutations as they affect more than one gene, thereby reducing the evolutionary rate of RNA viruses and hence their adaptive capacity. Although there is general agreement on the benefits of gene overlapping as a mechanism of genomic compression for rapidly evolving organisms, its effect on the pace of RNA virus evolution remains a source of debate. To address this issue, we collected sequence data from 117 instances of gene overlapping across 19 families, 30 genera, and 55 species of RNA viruses. On these data, we analyzed how genetic distances, selective pressures, and the distribution of RNA secondary structures and conserved protein functional domains vary between overlapping (OV) and nonoverlapping (NOV) regions. We show that gene overlapping generally results in a decrease in the rate of RNA virus evolution through a reduction in the frequency of synonymous mutations. However, this effect is less pronounced in genes with a terminal rather than an internal gene overlap, which might result from a greater proportion of protein functional conserved domains in NOV than in OV regions, in turn reducing the number of nonsynonymous mutations in the former. Overall, our analyses clarify the role of gene overlapping as a modulator of the evolutionary rates exhibited by RNA viruses and shed light on the factors that shape the genetic diversity of this important group of pathogens.

Sounds of silence: synonymous nucleotides as a key to biological regulation and complexity

Messenger RNA is a key component of an intricate regulatory network of its own. It accommodates numerous nucleotide signals that overlap protein coding sequences and are responsible for multiple levels of regulation and generation of biological complexity. A wealth of structural and regulatory information, which mRNA carries in addition to the encoded amino acid sequence, raises the question of how these signals and overlapping codes are delineated along non-synonymous and synonymous positions in protein coding regions, especially in eukaryotes. Silent or synonymous codon positions, which do not determine amino acid sequences of the encoded proteins, define mRNA secondary structure and stability and affect the rate of translation, folding and post-translational modifications of nascent polypeptides. The RNA level selection is acting on synonymous sites in both prokaryotes and eukaryotes and is more common than previously thought. Selection pressure on the coding gene regions follows three-nucleotide periodic pattern of nucleotide base-pairing in mRNA, which is imposed by the genetic code. Synonymous positions of the coding regions have a higher level of hybridization potential relative to non-synonymous positions, and are multifunctional in their regulatory and structural roles. Recent experimental evidence and analysis of mRNA structure and interspecies conservation suggest that there is an evolutionary tradeoff between selective pressure acting at the RNA and protein levels. Here we provide a comprehensive overview of the studies that define the role of silent positions in regulating RNA structure and processing that exert downstream effects on proteins and their functions.

Components of selection in the evolution of the influenza virus: linkage effects beat inherent selection

The influenza virus is an important human pathogen, with a rapid rate of evolution in the human population. The rate of homologous recombination within genes of influenza is essentially zero. As such, where two alleles within the same gene are in linkage disequilibrium, interference between alleles will occur, whereby selection acting upon one allele has an influence upon the frequency of the other. We here measured the relative importance of selection and interference effects upon the evolution of influenza. We considered time-resolved allele frequency data from the global evolutionary history of the haemagglutinin gene of human influenza A/H3N2, conducting an in-depth analysis of sequences collected since 1996. Using a model that accounts for selection-caused interference between alleles in linkage disequilibrium, we estimated the inherent selective benefit of individual polymorphisms in the viral population. These inherent selection coefficients were in turn used to calculate the total selective effect of interference acting upon each polymorphism, considering the effect of the initial background upon which a mutation arose, and the subsequent effect of interference from other alleles that were under selection. Viewing events in retrospect, we estimated the influence of each of these components in determining whether a mutant allele eventually fixed or died in the global viral population. Our inherent selection coefficients, when combined across different regions of the protein, were consistent with previous measurements of dN/dS for the same system. Alleles going on to fix in the global population tended to be under more positive selection, to arise on more beneficial backgrounds, and to avoid strong negative interference from other alleles under selection. However, on average, the fate of a polymorphism was determined more by the combined influence of interference effects than by its inherent selection coefficient.

Success in life is the product of many factors. Inherent ability often underlies great achievement. But other factors may play their part. The circumstances a child is born into may help or hinder his or her progress. Later events also have their effect; a life may be influenced by a lucky break, or an unforeseen disaster. In this work, we examine the factors underlying success for mutations in the HA gene of human influenza virus A/H3N2, defining success as the attainment of a high frequency in the global population. We examined the history of the gene from 1968 until 2010. For each observed mutation, a mathematical model was used to estimate the inherent benefit or disadvantage it conferred to the virus. We calculated the advantageousness or otherwise of the background upon which it arose, and the subsequent effect of interference from other mutations under selection. We found that successful mutations tended to have an advantageous background, and were subsequently fortunate in avoiding negative events throughout their lifetime. Beneficial mutations were more likely to be successful. But a mutation's chances of success were influenced more by circumstances of birth and subsequent events, than by its inherent effect on the virus.

The influenza A segment 7 mRNA 3' splice site pseudoknot/hairpin family

The 3′ splice site of the influenza A segment 7 transcript is utilized to produce mRNA for the critical M2 ion-channel protein. In solution a 63 nt fragment that includes this region can adopt two conformations: a pseudoknot and a hairpin. In each conformation, the splice site, a binding site for the SF2/ASF exonic splicing enhancer and a polypyrimidine tract, each exists in a different structural context. The most dramatic difference occurs for the splice site. In the hairpin the splice site is between two residues that are involved in a 2 by 2 nucleotide internal loop. In the pseudoknot, however, these bases are canonically paired within one of the pseudoknotted helices. The conformational switching observed in this region has implications for the regulation of splicing of the segment 7 mRNA. A measure of stability of the structures also shows interesting trends with respect to host specificity: avian strains tend to be the most stable, followed by swine and then human.

Specific temperature-induced perturbations of secondary mRNA structures are associated with the cold-adapted temperature-sensitive phenotype of influenza A virus

For decades, cold-adapted, temperature-sensitive (ca/ts) strains of influenza A virus have been used as live attenuated vaccines. Due to their great public health importance it is crucial to understand the molecular mechanism(s) of cold adaptation and temperature sensitivity that are currently unknown. For instance, secondary RNA structures play important roles in influenza biology. Thus, we hypothesized that a relatively minor change in temperature (32-39°C) can lead to perturbations in influenza RNA structures and, that these structural perturbations may be different for mRNAs of the wild type (wt) and ca/ts strains. To test this hypothesis, we developed a novel in silico method that enables assessing whether two related RNA molecules would undergo (dis)similar structural perturbations upon temperature change. The proposed method allows identifying those areas within an RNA chain where dissimilarities of RNA secondary structures at two different temperatures are particularly pronounced, without knowing particular RNA shapes at either temperature. We identified such areas in the NS2, PA, PB2 and NP mRNAs. However, these areas are not identical for the wt and ca/ts mutants. Differences in temperature-induced structural changes of wt and ca/ts mRNA structures may constitute a yet unappreciated molecular mechanism of the cold adaptation/temperature sensitivity phenomena.

The 3' splice site of influenza A segment 7 mRNA can exist in two conformations: a pseudoknot and a hairpin

The 3′ splice site of influenza A segment 7 is used to produce mRNA for the M2 ion-channel protein, which is critical to the formation of viable influenza virions. Native gel analysis, enzymatic/chemical structure probing, and oligonucleotide binding studies of a 63 nt fragment, containing the 3′ splice site, key residues of an SF2/ASF splicing factor binding site, and a polypyrimidine tract, provide evidence for an equilibrium between pseudoknot and hairpin structures. This equilibrium is sensitive to multivalent cations, and can be forced towards the pseudoknot by addition of 5 mM cobalt hexammine. In the two conformations, the splice site and other functional elements exist in very different structural environments. In particular, the splice site is sequestered in the middle of a double helix in the pseudoknot conformation, while in the hairpin it resides in a two-by-two nucleotide internal loop. The results suggest that segment 7 mRNA splicing can be controlled by a conformational switch that exposes or hides the splice site.

Influenza A virus coding regions exhibit host-specific global ordered RNA structure

Influenza A is a significant public health threat, partially because of its capacity to readily exchange gene segments between different host species to form novel pandemic strains. An understanding of the fundamental factors providing species barriers between different influenza hosts would facilitate identification of strains capable of leading to pandemic outbreaks and could also inform vaccine development. Here, we describe the difference in predicted RNA secondary structure stability that exists between avian, swine and human coding regions. The results predict that global ordered RNA structure exists in influenza A segments 1, 5, 7 and 8, and that ranges of free energies for secondary structure formation differ between host strains. The predicted free energy distributions for strains from avian, swine, and human species suggest criteria for segment reassortment and strains that might be ideal candidates for viral attenuation and vaccine development.