Subtype-specific structural constraints in the evolution of influenza A virus hemagglutinin genes

Locations and structures of influenza A virus packaging-associated signals and other functional elements via an in silico pipeline for predicting constrained features in RNA viruses

Influenza A virus contains regions of its segmented genome associated with ability to package the segments into virions, but many such regions are poorly characterised. We provide detailed predictions of the key locations within these packaging-associated regions, and their structures, by applying a recently-improved pipeline for delineating constrained regions in RNA viruses and applying structural prediction algorithms. We find and characterise other known constrained regions within influenza A genomes, including the region associated with the PA-X frameshift, regions associated with alternative splicing, and constraint around the initiation motif for a truncated PB1 protein, PB1-N92, associated with avian viruses. We further predict the presence of constrained regions that have not previously been described. The extra characterisation our work provides allows investigation of these key regions for drug target potential, and points towards determinants of packaging compatibility between segments.

Transient RNA structures underlie highly pathogenic avian influenza virus genesis

Highly pathogenic avian influenza viruses (HPAIVs) cause severe disease and high fatality in poultry1. They emerge exclusively from H5 and H7 low pathogenic avian influenza viruses (LPAIVs)2. Although insertion of a furin-cleavable multibasic cleavage site (MBCS) in the hemagglutinin gene was identified decades ago as the genetic basis for LPAIV-to-HPAIV transition3,4, the exact mechanisms underlying said insertion have remained unknown. Here we used an innovative combination of bioinformatic models to predict RNA structures forming around the influenza virus RNA polymerase during replication, and circular sequencing5 to reliably detect nucleotide insertions. We show that transient H5 hemagglutinin RNA structures predicted to trap the polymerase on purine-rich sequences drive nucleotide insertions characteristic of MBCSs, providing the first strong empirical evidence of RNA structure involvement in MBCS acquisition. Insertion frequencies at the H5 cleavage site were strongly affected by substitutions in flanking genomic regions altering predicted transient RNA structures. Introduction of H5-like cleavage site sequences and structures into an H6 hemagglutinin resulted in MBCS-yielding insertions never observed before in H6 viruses. Our results demonstrate that nucleotide insertions that underlie H5 HPAIV emergence result from a previously unknown RNA-structure-driven diversity-generating mechanism, which could be shared with other RNA viruses.

Discovery of RNA secondary structural motifs using sequence-ordered thermodynamic stability and comparative sequence analysis

Major advances in RNA secondary structural motif prediction have been achieved in the last few years; however, few methods harness the predictive power of multiple approaches to deliver in-depth characterizations of local RNA motifs and their potential functionality. Additionally, most available methods do not predict RNA pseudoknots. This work combines complementary bioinformatic systems into one robust discovery pipeline where: •RNA sequences are folded to search for thermodynamically favorable motifs utilizing ScanFold.•Motifs are expanded and refolded into alternate pseudoknot conformations by Knotty/Iterative HFold.•All conformations are evaluated for covariance via the cm-builder pipeline (Infernal and R-scape).

Comprehensive in virio structure probing analysis of the influenza A virus identifies functional RNA structures involved in viral genome replication

The influenza A virus genome is segmented into eight viral RNAs (vRNA). Secondary structures of vRNA are known to be involved in the viral proliferation process. Comprehensive vRNA structures in vitro, in virio, and in cellulo have been analyzed. However, the resolution of the structure map can be improved by comparative analysis and statistical modeling. Construction of a more high-resolution and reliable RNA structure map can identify uncharacterized functional structure motifs on vRNA in virion. Here, we establish the global map of the vRNA secondary structure in virion using the combination of dimethyl sulfate (DMS)-seq and selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE)-seq with a robust statistical analysis. Our high-resolution analysis identified a stem-loop structure at nucleotide positions 39 - 60 of segment 6 and further validated the structure at nucleotide positions 87 - 130 of segment 5 that was previously predicted to form a pseudoknot structure in silico. Notably, when the cells were infected with recombinant viruses which possess the mutations to disrupt the structure, the replication and packaging of the viral genome were drastically decreased. Our results provide comprehensive and high-resolution information on the influenza A virus genome structures in virion and evidence that the functional RNA structure motifs on the influenza A virus genome are associated with appropriate replication and packaging of the viral genome.

Structural and Functional RNA Motifs of SARS-CoV-2 and Influenza A Virus as a Target of Viral Inhibitors

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the COVID-19 pandemic, whereas the influenza A virus (IAV) causes seasonal epidemics and occasional pandemics. Both viruses lead to widespread infection and death. SARS-CoV-2 and the influenza virus are RNA viruses. The SARS-CoV-2 genome is an approximately 30 kb, positive sense, 5' capped single-stranded RNA molecule. The influenza A virus genome possesses eight single-stranded negative-sense segments. The RNA secondary structure in the untranslated and coding regions is crucial in the viral replication cycle. The secondary structure within the RNA of SARS-CoV-2 and the influenza virus has been intensively studied. Because the whole of the SARS-CoV-2 and influenza virus replication cycles are dependent on RNA with no DNA intermediate, the RNA is a natural and promising target for the development of inhibitors. There are a lot of RNA-targeting strategies for regulating pathogenic RNA, such as small interfering RNA for RNA interference, antisense oligonucleotides, catalytic nucleic acids, and small molecules. In this review, we summarized the knowledge about the inhibition of SARS-CoV-2 and influenza A virus propagation by targeting their RNA secondary structure.

The influenza A virus genome packaging network - complex, flexible and yet unsolved

The genome of influenza A virus (IAV) consists of eight unique viral RNA segments. This genome organization allows genetic reassortment between co-infecting IAV strains, whereby new IAVs with altered genome segment compositions emerge. While it is known that reassortment events can create pandemic IAVs, it remains impossible to anticipate reassortment outcomes with pandemic prospects. Recent research indicates that reassortment is promoted by a viral genome packaging mechanism that delivers the eight genome segments as a supramolecular complex into the virus particle. This finding holds promise of predicting pandemic IAVs by understanding the intermolecular interactions governing this genome packaging mechanism. Here, we critically review the prevailing mechanistic model postulating that IAV genome packaging is orchestrated by a network of intersegmental RNA-RNA interactions. Although we find supporting evidence, including segment-specific packaging signals and experimentally proposed RNA-RNA interaction networks, this mechanistic model remains debatable due to a current shortage of functionally validated intersegmental RNA-RNA interactions. We speculate that identifying such functional intersegmental RNA-RNA contacts might be hampered by limitations of the utilized probing techniques and the inherent complexity of the genome packaging mechanism. Nevertheless, we anticipate that improved probing strategies combined with a mutagenesis-based validation could facilitate their discovery.

Hemagglutinin Subtype Specificity and Mechanisms of Highly Pathogenic Avian Influenza Virus Genesis

Highly Pathogenic Avian Influenza Viruses (HPAIVs) arise from low pathogenic precursors following spillover from wild waterfowl into poultry populations. The main virulence determinant of HPAIVs is the presence of a multi-basic cleavage site (MBCS) in the hemagglutinin (HA) glycoprotein. The MBCS allows for HA cleavage and, consequently, activation by ubiquitous proteases, which results in systemic dissemination in terrestrial poultry. Since 1959, 51 independent MBCS acquisition events have been documented, virtually all in HA from the H5 and H7 subtypes. In the present article, data from natural LPAIV to HPAIV conversions and experimental in vitro and in vivo studies were reviewed in order to compile recent advances in understanding HA cleavage efficiency, protease usage, and MBCS acquisition mechanisms. Finally, recent hypotheses that might explain the unique predisposition of the H5 and H7 HA sequences to obtain an MBCS in nature are discussed.

In Silico Analyses of the Role of Codon Usage at the Hemagglutinin Cleavage Site in Highly Pathogenic Avian Influenza Genesis

A vast diversity of 16 influenza hemagglutinin (HA) subtypes are found in birds. Interestingly, viruses from only two subtypes, H5 and H7, have so far evolved into highly pathogenic avian influenza viruses (HPAIVs) following insertions or substitutions at the HA cleavage site by the viral polymerase. The mechanisms underlying this striking subtype specificity are still unknown. Here, we compiled a comprehensive dataset of 20,488 avian influenza virus HA sequences to investigate differences in nucleotide and amino acid usage at the HA cleavage site between subtypes and how these might impact the genesis of HPAIVs by polymerase stuttering and realignment. We found that sequences of the H5 and H7 subtypes stand out by their high purine content at the HA cleavage site. In addition, fewer substitutions were necessary in H5 and H7 HAs than in HAs from other subtypes to acquire an insertion-prone HA cleavage site sequence, as defined based on in vitro and in vivo data from the literature. Codon usage was more favorable for HPAIV genesis in sequences of viruses isolated from species or geographical regions in which HPAIV genesis is more frequently observed in nature. The results of the present analyses suggest that the subtype restriction of HPAIV genesis to H5 and H7 influenza viruses might be due to the particular codon usage at the HA cleavage site in these subtypes.

In silico analysis of local RNA secondary structure in influenza virus A, B and C finds evidence of widespread ordered stability but little evidence of significant covariation

Influenza virus is a persistent threat to human health; indeed, the deadliest modern pandemic was in 1918 when an H1N1 virus killed an estimated 50 million people globally. The intent of this work is to better understand influenza from an RNA-centric perspective to provide local, structural motifs with likely significance to the influenza infectious cycle for therapeutic targeting. To accomplish this, we analyzed over four hundred thousand RNA sequences spanning three major clades: influenza A, B and C. We scanned influenza segments for local secondary structure, identified/modeled motifs of likely functionality, and coupled the results to an analysis of evolutionary conservation. We discovered 185 significant regions of predicted ordered stability, yet evidence of sequence covariation was limited to 7 motifs, where 3-found in influenza C-had higher than expected amounts of sequence covariation.

Structural Requirements in the Hemagglutinin Cleavage Site-Coding RNA Region for the Generation of Highly Pathogenic Avian Influenza Virus

Highly pathogenic avian influenza viruses (HPAIVs) with H5 and H7 hemagglutinin (HA) subtypes are derived from their low pathogenic counterparts following the acquisition of multiple basic amino acids in their HA cleavage site. It has been suggested that consecutive adenine residues and a stem-loop structure in the viral RNA region that encodes the cleavage site are essential for the acquisition of the polybasic cleavage site. By using a reporter assay to detect non-templated nucleotide insertions, we found that insertions more frequently occurred in the RNA region (29 nucleotide-length) encoding the cleavage site of an H5 HA gene that was predicted to have a stem-loop structure containing consecutive adenines than in a mutated corresponding RNA region that had a disrupted loop structure with fewer adenines. In virus particles generated by using reverse genetics, nucleotide insertions that created additional codons for basic amino acids were found in the RNA region encoding the cleavage site of an H5 HA gene but not in the mutated RNA region. We confirmed the presence of virus clones with the ability to replicate without trypsin in a plaque assay and to cause lethal infection in chicks. These results demonstrate that the stem-loop structure containing consecutive adenines in HA genes is a key molecular determinant for the emergence of H5 HPAIVs.

Universal and strain specific structure features of segment 8 genomic RNA of influenza A virus-application of 4-thiouridine photocrosslinking

RNA structure in the influenza A virus (IAV) has been the focus of several studies that have shown connections between conserved secondary structure motifs and their biological function in the virus replication cycle. Questions have arisen on how to best recognize and understand the pandemic properties of IAV strains from an RNA perspective, but determination of the RNA secondary structure has been challenging. Herein, we used chemical mapping to determine the secondary structure of segment 8 viral RNA (vRNA) of the pandemic A/California/04/2009 (H1N1) strain of IAV. Additionally, this long, naturally occurring RNA served as a model to evaluate RNA mapping with 4-thiouridine (4sU) crosslinking. We explored 4-thiouridine as a probe of nucleotides in close proximity, through its incorporation into newly transcribed RNA and subsequent photoactivation. RNA secondary structural features both universal to type A strains and unique to the A/California/04/2009 (H1N1) strain were recognized. 4sU mapping confirmed and facilitated RNA structure prediction, according to several rules: 4sU photocross-linking forms efficiently in the double-stranded region of RNA with some flexibility, in the ends of helices, and across bulges and loops when their structural mobility is permitted. This method highlighted three-dimensional properties of segment 8 vRNA secondary structure motifs and allowed to propose several long-range three-dimensional interactions. 4sU mapping combined with chemical mapping and bioinformatic analysis could be used to enhance the RNA structure determination as well as recognition of target regions for antisense strategies or viral RNA detection.

Phylodynamic Study of the Conserved RNA Structure Encompassing the Hemagglutinin Cleavage Site Encoding Region of H5 and H7 Low Pathogenic Avian Influenza Viruses

Highly Pathogenic Avian Influenza Viruses (HPAIV) evolve from Low Pathogenic Avian Influenza Viruses (LPAIV) of the H5 and H7 subtypes. This evolution is characterized by the acquisition of a multi-basic cleavage site (MBCS) motif in the hemagglutinin (HA) that leads to an extended viral tropism and severe disease in poultry. One key unanswered question is whether the risk of transition to HPAIV is similar for all LPAIV H5 or H7 strains, or whether specific determinants in the HA sequence of some H5 or H7 LPAIV strains correlate with a higher risk of transition to HPAIV. Here we determined if specific features of the conserved RNA stem loop located at the hemagglutinin cleavage site-encoding region could be detected along the LPAIV to HPAIV evolutionary pathway. Analysis of the thermodynamic stability of the predicted RNA structures showed no specific patterns common to HA sequences leading to HPAIV and distinct from those remaining LPAIV. However, RNA structure clustering analysis revealed that most of the American lineage ancestors leading to H7 emergences via recombination shared the same vRNA structure topology at the HA1/HA2 boundary region. Our study thus identified predicted secondary RNA structures present in the HA of H7 viruses, which could promote genetic recombination and acquisition of a MBCS.

Conserved Structural Motifs of Two Distant IAV Subtypes in Genomic Segment 5 RNA

The functionality of RNA is fully dependent on its structure. For the influenza A virus (IAV), there are confirmed structural motifs mediating processes which are important for the viral replication cycle, including genome assembly and viral packaging. Although the RNA of strains originating from distant IAV subtypes might fold differently, some structural motifs are conserved, and thus, are functionally important. Nowadays, NGS-based structure modeling is a source of new in vivo data helping to understand RNA biology. However, for accurate modeling of in vivo RNA structures, these high-throughput methods should be supported with other analyses facilitating data interpretation. In vitro RNA structural models complement such approaches and offer RNA structures based on experimental data obtained in a simplified environment, which are needed for proper optimization and analysis. Herein, we present the secondary structure of the influenza A virus segment 5 vRNA of A/California/04/2009 (H1N1) strain, based on experimental data from DMS chemical mapping and SHAPE using NMIA, supported by base-pairing probability calculations and bioinformatic analyses. A comparison of the available vRNA5 structures among distant IAV strains revealed that a number of motifs present in the A/California/04/2009 (H1N1) vRNA5 model are highly conserved despite sequence differences, located within previously identified packaging signals, and the formation of which in in virio conditions has been confirmed. These results support functional roles of the RNA secondary structure motifs, which may serve as candidates for universal RNA-targeting inhibitory methods.

Insertions of codons encoding basic amino acids in H7 hemagglutinins of influenza A viruses occur by recombination with RNA at hotspots near snoRNA binding sites

The presence of multiple basic amino acids in the protease cleavage site of the hemagglutinin (HA) protein is the main molecular determinant of virulence of highly pathogenic avian influenza (HPAI) viruses. Recombination of HA RNA with other RNA molecules of host or virus origin is a dominant mechanism of multibasic cleavage site (MBCS) acquisition for H7 subtype HA. Using alignments of HA RNA sequences from documented cases of MBCS insertion due to recombination, we show that such recombination with host RNAs is most likely to occur at particular hotspots in ribosomal RNAs (rRNAs), transfer RNAs (tRNAs), and viral RNAs. The locations of these hotspots in highly abundant RNAs indicate that RNA recombination is facilitated by the binding of small nucleolar RNA (snoRNA) near the recombination points.

A risk marker of tribasic hemagglutinin cleavage site in influenza A (H9N2) virus

Low pathogenic avian influenza A(H9N2) virus is endemic worldwide and continually recruit internal genes to generate human-infecting H5N1, H5N6, H7N9, and H10N8 influenza variants. Here we show that hemagglutinin cleavage sites (HACS) of H9N2 viruses tended to mutate towards hydrophilic via evolutionary transition, and the tribasic HACS were found at high prevalence in Asia and the Middle East. Our finding suggested that the tribasic H9N2 viruses increased the viral replication, stability, pathogenicity and transmission in chickens and the virulence of mice compared to the monobasic H9N2 viruses. Notably, the enlarged stem-loop structures of HACS in the RNA region were found in the increasing tribasic H9N2 viruses. The enlarged HACS RNA secondary structures of H9N2 viruses did not influence the viral replication but accelerated the frequency of nucleotide insertion in HACS. With the prevailing tendency of the tribasic H9N2 viruses, the tribasic HACS in H9N2 viruses should be paid more attention.

Emergence and Selection of a Highly Pathogenic Avian Influenza H7N3 Virus

Low-pathogenicity avian influenza (LPAI) viruses of subtypes H5 and H7 have the ability to spontaneously mutate to highly pathogenic (HPAI) virus variants, causing high mortality in poultry. The highly pathogenic phenotype is caused by mutation of the hemagglutinin (HA) cleavage site, but additional mutations may play a role. Evidence from the field for the switch to high pathogenicity remains scarce. This study provides direct evidence for LPAI-to-HPAI virus mutation during H7N3 infection of a turkey farm in the Netherlands. No severe clinical symptoms were reported at the farm, but deep sequencing of isolates from the infected turkeys revealed a minority of HPAI virus sequences (0.06%) in the virus population. The HPAI virus contained a 12-nucleotide insertion in the HA cleavage site that was likely introduced by a single event as no intermediates with shorter inserts were identified. This suggests nonhomologous recombination as the mechanism of insertion. Analysis of different organs of the infected turkeys showed the largest amount of HPAI virus in the lung (4.4%). The HPAI virus was rapidly selected in experimentally infected chickens after both intravenous and intranasal/intratracheal inoculation with a mixed virus preparation. Full-genome sequencing revealed that both pathotypes contained a deletion in the stalk region of the neuraminidase protein. We identified additional mutations in HA and polymerase basic protein 1 (PB1) in the HPAI virus, which were already present as minority variants in the LPAI virus population. Our findings provide more insight into the molecular changes and mechanisms involved in the emergence and selection of HPAI viruses.IMPORTANCE Low-pathogenicity avian influenza (LPAI) viruses circulate in wild birds and can be transmitted to poultry. LPAI viruses can mutate to become highly pathogenic avian influenza (HPAI) viruses causing severe disease and death in poultry. Little is known about this switch to high pathogenicity. We isolated an LPAI H7N3 virus from an infected turkey farm and showed that this contains small amounts of HPAI virus. The HPAI virus rapidly outcompeted the LPAI virus in chickens that were experimentally infected with this mixture of viruses. We analyzed the genome sequences of the LPAI and HPAI viruses and identified several changes that may be important for a virus to become highly pathogenic. This knowledge may be used for timely identification of LPAI viruses that pose a risk of becoming highly pathogenic in the field.

Host-Virus Interaction: How Host Cells Defend against Influenza A Virus Infection

Influenza A viruses (IAVs) are highly contagious pathogens infecting human and numerous animals. The viruses cause millions of infection cases and thousands of deaths every year, thus making IAVs a continual threat to global health. Upon IAV infection, host innate immune system is triggered and activated to restrict virus replication and clear pathogens. Subsequently, host adaptive immunity is involved in specific virus clearance. On the other hand, to achieve a successful infection, IAVs also apply multiple strategies to avoid be detected and eliminated by the host immunity. In the current review, we present a general description on recent work regarding different host cells and molecules facilitating antiviral defenses against IAV infection and how IAVs antagonize host immune responses.

Conserved structural RNA domains in regions coding for cleavage site motifs in hemagglutinin genes of influenza viruses

The acquisition of a multibasic cleavage site (MBCS) in the hemagglutinin (HA) glycoprotein is the main determinant of the conversion of low pathogenic avian influenza viruses into highly pathogenic strains, facilitating HA cleavage and virus replication in a broader range of host cells. In nature, substitutions or insertions in HA RNA genomic segments that code for multiple basic amino acids have been observed only in the HA genes of two out of sixteen subtypes circulating in birds, H5 and H7. Given the compatibility of MBCS motifs with HA proteins of numerous subtypes, this selectivity was hypothesized to be determined by the existence of specific motifs in HA RNA, in particular structured domains. In H5 and H7 HA RNAs, predictions of such domains have yielded alternative conserved stem-loop structures with the cleavage site codons in the hairpin loops. Here, potential RNA secondary structures were analyzed in the cleavage site regions of HA segments of influenza viruses of different types and subtypes. H5- and H7-like stem-loop structures were found in all known influenza A virus subtypes and in influenza B and C viruses with homology modeling. Nucleotide covariations supported this conservation to be determined by RNA structural constraints that are stronger in the domain-closing bottom stems as compared to apical parts. The structured character of this region in (sub-)types other than H5 and H7 indicates its functional importance beyond the ability to evolve toward an MBCS responsible for a highly pathogenic phenotype.

RNA Sequence Features Are at the Core of Influenza A Virus Genome Packaging

The influenza A virus (IAV), a respiratory pathogen for humans, poses serious medical and economic challenges to global healthcare systems. The IAV genome, consisting of eight single-stranded viral RNA segments, is incorporated into virions by a complex process known as genome packaging. Specific RNA sequences within the viral RNA segments serve as signals that are necessary for genome packaging. Although efficient packaging is a prerequisite for viral infectivity, many of the mechanistic details about this process are still missing. In this review, we discuss the recent advances toward the understanding of IAV genome packaging and focus on the RNA features that play a role in this process.

Structural and Functional Motifs in Influenza Virus RNAs

Influenza A viruses (IAV) are responsible for recurrent influenza epidemics and occasional devastating pandemics in humans and animals. They belong to the Orthomyxoviridae family and their genome consists of eight (-) sense viral RNA (vRNA) segments of different lengths coding for at least 11 viral proteins. A heterotrimeric polymerase complex is bound to the promoter consisting of the 13 5′-terminal and 12 3′-terminal nucleotides of each vRNA, while internal parts of the vRNAs are associated with multiple copies of the viral nucleoprotein (NP), thus forming ribonucleoproteins (vRNP). Transcription and replication of vRNAs result in viral mRNAs (vmRNAs) and complementary RNAs (cRNAs), respectively. Complementary RNAs are the exact positive copies of vRNAs; they also form ribonucleoproteins (cRNPs) and are intermediate templates in the vRNA amplification process. On the contrary, vmRNAs have a 5′ cap snatched from cellular mRNAs and a 3′ polyA tail, both gained by the viral polymerase complex. Hence, unlike vRNAs and cRNAs, vmRNAs do not have a terminal promoter able to recruit the viral polymerase. Furthermore, synthesis of at least two viral proteins requires vmRNA splicing. Except for extensive analysis of the viral promoter structure and function and a few, mostly bioinformatics, studies addressing the vRNA and vmRNA structure, structural studies of the influenza A vRNAs, cRNAs, and vmRNAs are still in their infancy. The recent crystal structures of the influenza polymerase heterotrimeric complex drastically improved our understanding of the replication and transcription processes. The vRNA structure has been mainly studied in vitro using RNA probing, but its structure has been very recently studied within native vRNPs using crosslinking and RNA probing coupled to next generation RNA sequencing. Concerning vmRNAs, most studies focused on the segment M and NS splice sites and several structures initially predicted by bioinformatics analysis have now been validated experimentally and their role in the viral life cycle demonstrated. This review aims to compile the structural motifs found in the different RNA classes (vRNA, cRNA, and vmRNA) of influenza viruses and their function in the viral replication cycle.

RNA2DMut: a web tool for the design and analysis of RNA structure mutations

With the widespread application of high-throughput sequencing, novel RNA sequences are being discovered at an astonishing rate. The analysis of function, however, lags behind. In both the cis- and trans-regulatory functions of RNA, secondary structure (2D base-pairing) plays essential regulatory roles. In order to test RNA function, it is essential to be able to design and analyze mutations that can affect structure. This was the motivation for the creation of the RNA2DMut web tool. With RNA2DMut, users can enter in RNA sequences to analyze, constrain mutations to specific residues, or limit changes to purines/pyrimidines. The sequence is analyzed at each base to determine the effect of every possible point mutation on 2D structure. The metrics used in RNA2DMut rely on the calculation of the Boltzmann structure ensemble and do not require a robust 2D model of RNA structure for designing mutations. This tool can facilitate a wide array of uses involving RNA: for example, in designing and evaluating mutants for biological assays, interrogating RNA-protein interactions, identifying key regions to alter in SELEX experiments, and improving RNA folding and crystallization properties for structural biology. Additional tools are available to help users introduce other mutations (e.g., indels and substitutions) and evaluate their effects on RNA structure. Example calculations are shown for five RNAs that require 2D structure for their function: the MALAT1 mascRNA, an influenza virus splicing regulatory motif, the EBER2 viral noncoding RNA, the Xist lncRNA repA region, and human Y RNA 5. RNA2DMut can be accessed at https://rna2dmut.bb.iastate.edu/.

Synonymous Mutations at the Beginning of the Influenza A Virus Hemagglutinin Gene Impact Experimental Fitness

The fitness effects of synonymous mutations can provide insights into biological and evolutionary mechanisms. We analyzed the experimental fitness effects of all single-nucleotide mutations, including synonymous substitutions, at the beginning of the influenza A virus hemagglutinin (HA) gene. Many synonymous substitutions were deleterious both in bulk competition and for individually isolated clones. Investigating protein and RNA levels of a subset of individually expressed HA variants revealed that multiple biochemical properties contribute to the observed experimental fitness effects. Our results indicate that a structural element in the HA segment viral RNA may influence fitness. Examination of naturally evolved sequences in human hosts indicates a preference for the unfolded state of this structural element compared to that found in swine hosts. Our overall results reveal that synonymous mutations may have greater fitness consequences than indicated by simple models of sequence conservation, and we discuss the implications of this finding for commonly used evolutionary tests and analyses.

Know Your Enemy: Successful Bioinformatic Approaches to Predict Functional RNA Structures in Viral RNAs

Structured RNA elements may control virus replication, transcription and translation, and their distinct features are being exploited by novel antiviral strategies. Viral RNA elements continue to be discovered using combinations of experimental and computational analyses. However, the wealth of sequence data, notably from deep viral RNA sequencing, viromes, and metagenomes, necessitates computational approaches being used as an essential discovery tool. In this review, we describe practical approaches being used to discover functional RNA elements in viral genomes. In addition to success stories in new and emerging viruses, these approaches have revealed some surprising new features of well-studied viruses e.g., human immunodeficiency virus, hepatitis C virus, influenza, and dengue viruses. Some notable discoveries were facilitated by new comparative analyses of diverse viral genome alignments. Importantly, comparative approaches for finding RNA elements embedded in coding and non-coding regions differ. With the exponential growth of computer power we have progressed from stem-loop prediction on single sequences to cutting edge 3D prediction, and from command line to user friendly web interfaces. Despite these advances, many powerful, user friendly prediction tools and resources are underutilized by the virology community.

A compensatory mutagenesis study of a conserved hairpin in the M gene segment of influenza A virus shows its role in virus replication

RNA structures are increasingly recognized to be of importance during influenza A virus replication. Here, we investigated a predicted conserved hairpin in the M gene segment (nt 967-994) within the region of the vRNA 5′ packaging signal. The existence of this RNA structure and its possible role in virus replication was investigated using a compensatory mutagenesis approach. Mutations were introduced in the hairpin stem, based on natural variation. Virus replication properties were studied for the mutant viruses with disrupted and restored RNA structures. Viruses with structure-disrupting mutations had lower virus titers and a significantly reduced median plaque size when compared with the wild-type (WT) virus, while viruses with structure restoring-mutations replicated comparable to WT. Moreover, virus replication was also reduced when mutations were introduced in the hairpin loop, suggesting its involvement in RNA interactions. Northern blot and FACS experiments were performed to study differences in RNA levels as well as production of M1 and M2 proteins, expressed via alternative splicing. Stem-disruptive mutants caused lower vRNA and M2 mRNA levels and reduced M2 protein production at early time-points. When the RNA structure was restored, vRNA, M2 mRNA and M2 protein levels were increased, demonstrating a compensatory effect. Thus, this study provides evidence for functional importance of the predicted M RNA structure and suggests its role in splicing regulation.

Heterologous Packaging Signals on Segment 4, but Not Segment 6 or Segment 8, Limit Influenza A Virus Reassortment

Influenza A virus (IAV) RNA packaging signals serve to direct the incorporation of IAV gene segments into virus particles, and this process is thought to be mediated by segment-segment interactions. These packaging signals are segment and strain specific, and as such, they have the potential to impact reassortment outcomes between different IAV strains. Our study aimed to quantify the impact of packaging signal mismatch on IAV reassortment using the human seasonal influenza A/Panama/2007/99 (H3N2) and pandemic influenza A/Netherlands/602/2009 (H1N1) viruses. Focusing on the three most divergent segments, we constructed pairs of viruses that encoded identical proteins but differed in the packaging signal regions on a single segment. We then evaluated the frequency with which segments carrying homologous versus heterologous packaging signals were incorporated into reassortant progeny viruses. We found that, when segment 4 (HA) of coinfecting parental viruses was modified, there was a significant preference for the segment containing matched packaging signals relative to the background of the virus. This preference was apparent even when the homologous HA constituted a minority of the HA segment population available in the cell for packaging. Conversely, when segment 6 (NA) or segment 8 (NS) carried modified packaging signals, there was no significant preference for homologous packaging signals. These data suggest that movement of NA and NS segments between the human H3N2 and H1N1 lineages is unlikely to be restricted by packaging signal mismatch, while movement of the HA segment would be more constrained. Our results indicate that the importance of packaging signals in IAV reassortment is segment dependent.IMPORTANCE Influenza A viruses (IAVs) can exchange genes through reassortment. This process contributes to both the highly diverse population of IAVs found in nature and the formation of novel epidemic and pandemic IAV strains. Our study sought to determine the extent to which IAV packaging signal divergence impacts reassortment between seasonal IAVs. Our knowledge in this area is lacking, and insight into the factors that influence IAV reassortment will inform and strengthen ongoing public health efforts to anticipate the emergence of new viruses. We found that the packaging signals on the HA segment, but not the NA or NS segments, restricted IAV reassortment. Thus, the packaging signals of the HA segment could be an important factor in determining the likelihood that two IAV strains of public health interest will undergo reassortment.