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

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.

TFOFinder: Python program for identifying purine-only double-stranded stretches in the predicted secondary structure(s) of RNA targets

Nucleic acid probes are valuable tools in biology and chemistry and are indispensable for PCR amplification of DNA, RNA quantification and visualization, and downregulation of gene expression. Recently, triplex-forming oligonucleotides (TFO) have received increased attention due to their improved selectivity and sensitivity in recognizing purine-rich double-stranded RNA regions at physiological pH by incorporating backbone and base modifications. For example, triplex-forming peptide nucleic acid (PNA) oligomers have been used for imaging a structured RNA in cells and inhibiting influenza A replication. Although a handful of programs are available to identify triplex target sites (TTS) in DNA, none are available that find such regions in structured RNAs. Here, we describe TFOFinder, a Python program that facilitates the identification of intramolecular purine-only RNA duplexes that are amenable to forming parallel triple helices (pyrimidine/purine/pyrimidine) and the design of the corresponding TFO(s). We performed genome- and transcriptome-wide analyses of TTS in Drosophila melanogaster and found that only 0.3% (123) of total unique transcripts (35,642) show the potential of forming 12-purine long triplex forming sites that contain at least one guanine. Using minimization algorithms, we predicted the secondary structure(s) of these transcripts, and using TFOFinder, we found that 97 (79%) of the identified 123 transcripts are predicted to fold to form at least one TTS for parallel triple helix formation. The number of transcripts with potential purine TTS increases when the strict search conditions are relaxed by decreasing the length of the probe or by allowing up to two pyrimidine inversions or 1-nucleotide bulge in the target site. These results are encouraging for the use of modified triplex forming probes for live imaging of endogenous structured RNA targets, such as pre-miRNAs, and inhibition of target-specific translation and viral replication.

In vivo secondary structural analysis of Influenza A virus genomic RNA

Influenza A virus (IAV) is a respiratory virus that causes epidemics and pandemics. Knowledge of IAV RNA secondary structure in vivo is crucial for a better understanding of virus biology. Moreover, it is a fundament for the development of new RNA-targeting antivirals. Chemical RNA mapping using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) coupled with Mutational Profiling (MaP) allows for the thorough examination of secondary structures in low-abundance RNAs in their biological context. So far, the method has been used for analyzing the RNA secondary structures of several viruses including SARS-CoV-2 in virio and in cellulo. Here, we used SHAPE-MaP and dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) for genome-wide secondary structure analysis of viral RNA (vRNA) of the pandemic influenza A/California/04/2009 (H1N1) strain in both in virio and in cellulo environments. Experimental data allowed the prediction of the secondary structures of all eight vRNA segments in virio and, for the first time, the structures of vRNA5, 7, and 8 in cellulo. We conducted a comprehensive structural analysis of the proposed vRNA structures to reveal the motifs predicted with the highest accuracy. We also performed a base-pairs conservation analysis of the predicted vRNA structures and revealed many highly conserved vRNA motifs among the IAVs. The structural motifs presented herein are potential candidates for new IAV antiviral strategies.

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.

Targeting RNA structures with small molecules

RNA adopts 3D structures that confer varied functional roles in human biology and dysfunction in disease. Approaches to therapeutically target RNA structures with small molecules are being actively pursued, aided by key advances in the field including the development of computational tools that predict evolutionarily conserved RNA structures, as well as strategies that expand mode of action and facilitate interactions with cellular machinery. Existing RNA-targeted small molecules use a range of mechanisms including directing splicing - by acting as molecular glues with cellular proteins (such as branaplam and the FDA-approved risdiplam), inhibition of translation of undruggable proteins and deactivation of functional structures in noncoding RNAs. Here, we describe strategies to identify, validate and optimize small molecules that target the functional transcriptome, laying out a roadmap to advance these agents into the next decade.

Programmable antivirals targeting critical conserved viral RNA secondary structures from influenza A virus and SARS-CoV-2

Influenza A virus's (IAV's) frequent genetic changes challenge vaccine strategies and engender resistance to current drugs. We sought to identify conserved and essential RNA secondary structures within IAV's genome that are predicted to have greater constraints on mutation in response to therapeutic targeting. We identified and genetically validated an RNA structure (packaging stem-loop 2 (PSL2)) that mediates in vitro packaging and in vivo disease and is conserved across all known IAV isolates. A PSL2-targeting locked nucleic acid (LNA), administered 3 d after, or 14 d before, a lethal IAV inoculum provided 100% survival in mice, led to the development of strong immunity to rechallenge with a tenfold lethal inoculum, evaded attempts to select for resistance and retained full potency against neuraminidase inhibitor-resistant virus. Use of an analogous approach to target SARS-CoV-2, prophylactic administration of LNAs specific for highly conserved RNA structures in the viral genome, protected hamsters from efficient transmission of the SARS-CoV-2 USA_WA1/2020 variant. These findings highlight the potential applicability of this approach to any virus of interest via a process we term 'programmable antivirals', with implications for antiviral prophylaxis and post-exposure therapy.

Nuclear magnetic resonance reveals a two hairpin equilibrium near the 3'-splice site of influenza A segment 7 mRNA that can be shifted by oligonucleotides

Influenza A kills hundreds of thousands of people globally every year and has the potential to generate more severe pandemics. Influenza A's RNA genome and transcriptome provide many potential therapeutic targets. Here, nuclear magnetic resonance (NMR) experiments suggest that one such target could be a hairpin loop of 8 nucleotides in a pseudoknot that sequesters a 3' splice site in canonical pairs until a conformational change releases it into a dynamic 2 × 2-nt internal loop. NMR experiments reveal that the hairpin loop is dynamic and able to bind oligonucleotides as short as pentamers. A 3D NMR structure of the complex contains 4 and likely 5 bp between pentamer and loop. Moreover, a hairpin sequence was discovered that mimics the equilibrium of the influenza hairpin between its structure in the pseudoknot and upon release of the splice site. Oligonucleotide binding shifts the equilibrium completely to the hairpin secondary structure required for pseudoknot folding. The results suggest this hairpin can be used to screen for compounds that stabilize the pseudoknot and potentially reduce splicing.

Contribution of RNA-RNA Interactions Mediated by the Genome Packaging Signals for the Selective Genome Packaging of Influenza A Virus

The influenza A virus genome is composed of eight single-stranded negative-sense viral RNA segments (vRNAs). The eight vRNAs are selectively packaged into each progeny virion. This process likely involves specific interactions between the vRNAs via segment-specific packaging signals located in both the 3'- and 5'-terminal regions of the respective vRNAs. To assess the importance of vRNA-vRNA interactions via packaging signals for selective genome packaging, we generated mutant viruses possessing silent mutations in the packaging signal region of the hemagglutinin (HA) vRNA. A mutant virus possessing silent mutations in nucleotides (nt) 1664 to 1676 resulted in defects in HA vRNA incorporation and showed a reduction in viral growth. After serial passage, the mutant virus acquired additional mutations in the 5'-terminal packaging signal regions of both the HA and polymerase basic 2 (PB2) vRNAs. These mutations contributed to the recovery of viral growth and HA vRNA packaging efficiency. In addition, an RNA-RNA interaction between the 5' ends of HA and PB2 vRNAs was confirmed in vitro, and this interaction was disrupted following the introduction of silent mutations in the HA vRNA. Thus, our results demonstrated that RNA-RNA interactions between the packaging signal regions of HA vRNA and PB2 vRNA are important for selective genome packaging. IMPORTANCE While numerous viral genomes comprise a single genome segment, the influenza A virus possesses eight segmented genomes. Influenza A virus can benefit from having a segmented genome because the segments can reassort with other strains of the influenza virus to create new genetically distinct strains. The influenza A virus efficiently incorporates one copy of each of its eight genomic segments per viral particle. However, the mechanism by which each segment is specifically selected is poorly understood. The genome segments contain RNA signals that facilitate the incorporation of segments into virus particles. These regions may facilitate specific interactions between the genome segments, creating an eight-segment complex, which can then be packaged into individual particles. In this study, we provide evidence that RNA signals contribute to specific interactions between two of the influenza virus genome segments.

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.

RNA-Targeting Splicing Modifiers: Drug Development and Screening Assays

RNA splicing is an essential step in producing mature messenger RNA (mRNA) and other RNA species. Harnessing RNA splicing modifiers as a new pharmacological modality is promising for the treatment of diseases caused by aberrant splicing. This drug modality can be used for infectious diseases by disrupting the splicing of essential pathogenic genes. Several antisense oligonucleotide splicing modifiers were approved by the U.S. Food and Drug Administration (FDA) for the treatment of spinal muscular atrophy (SMA) and Duchenne muscular dystrophy (DMD). Recently, a small-molecule splicing modifier, risdiplam, was also approved for the treatment of SMA, highlighting small molecules as important warheads in the arsenal for regulating RNA splicing. The cellular targets of these approved drugs are all mRNA precursors (pre-mRNAs) in human cells. The development of novel RNA-targeting splicing modifiers can not only expand the scope of drug targets to include many previously considered “undruggable” genes but also enrich the chemical-genetic toolbox for basic biomedical research. In this review, we summarized known splicing modifiers, screening methods for novel splicing modifiers, and the chemical space occupied by the small-molecule splicing modifiers.

Mutations designed to modify the NS gene mRNA secondary structure affect influenza A pathogenicity <em>in vivo</em>

The influenza A virus genome consists of eight segments of negative-sense RNA that encode up to 18 proteins. During the process of viral replication, positive-sense (+)RNA (cRNA) or messenger RNA (mRNA) is synthesized. Today, there is only a partial understanding of the function of several secondary structures within vRNA and cRNA promoters, and splice sites in the M and NS genes. The most precise secondary structure of (+)RNA has been determined for the NS segment of influenza A virus.  The influenza A virus NS gene features two regions with a conserved mRNA secondary structure located near splice sites. Here, we compared 4 variants of the A/Puerto Rico/8/1934 strain featuring different combinations of secondary structures at the NS segment (+)RNA regions 82-148 and 497-564. We found that RNA structures did not affect viral replication in cell culture. However, one of the viruses demonstrated lower NS1 and NEP expression levels during early stage cell infection as well as reduced pathogenicity in mice compared to other variants. In particular, this virus is characterized by an RNA hairpin in the 82-148 region and a stable hairpin in the 497-564 region.

RNA Secondary Structure as a First Step for Rational Design of the Oligonucleotides towards Inhibition of Influenza A Virus Replication

Influenza is an important research subject around the world because of its threat to humanity. Influenza A virus (IAV) causes seasonal epidemics and sporadic, but dangerous pandemics. A rapid antigen changes and recombination of the viral RNA genome contribute to the reduced effectiveness of vaccination and anti-influenza drugs. Hence, there is a necessity to develop new antiviral drugs and strategies to limit the influenza spread. IAV is a single-stranded negative sense RNA virus with a genome (viral RNA—vRNA) consisting of eight segments. Segments within influenza virion are assembled into viral ribonucleoprotein (vRNP) complexes that are independent transcription-replication units. Each step in the influenza life cycle is regulated by the RNA and is dependent on its interplay and dynamics. Therefore, viral RNA can be a proper target to design novel therapeutics. Here, we briefly described examples of anti-influenza strategies based on the antisense oligonucleotide (ASO), small interfering RNA (siRNA), microRNA (miRNA) and catalytic nucleic acids. In particular we focused on the vRNA structure-function relationship as well as presented the advantages of using secondary structure information in predicting therapeutic targets and the potential future of this field.

Attenuation Methods for Live Vaccines

Vaccination was developed by Edward Jenner in 1796. Since then, vaccination and vaccine development research has been a hotspot of research in the scientific community. Various ways of vaccine development are successfully employed in mass production of vaccines. One of the most successful ways to generate vaccines is the method of virulence attenuation in pathogens. The attenuated strains of viruses, bacteria, and parasites are used as vaccines which elicit robust immune response and confers protection against virulent pathogens. This chapter brings together the most common and efficient ways of generating live attenuated vaccine strains in viruses, bacteria, and parasites.

RNA Secondary Structure Motifs of the Influenza A Virus as Targets for siRNA-Mediated RNA Interference

The influenza A virus is a human pathogen that poses a serious public health threat due to rapid antigen changes and emergence of new, highly pathogenic strains with the potential to become easily transmitted in the human population. The viral genome is encoded by eight RNA segments, and all stages of the replication cycle are dependent on RNA. In this study, we designed small interfering RNA (siRNA) targeting influenza segment 5 nucleoprotein (NP) mRNA structural motifs that encode important functions. The new criterion for choosing the siRNA target was the prediction of accessible regions based on the secondary structure of segment 5 (+)RNA. This design led to siRNAs that significantly inhibit influenza virus type A replication in Madin-Darby canine kidney (MDCK) cells. Additionally, chemical modifications with the potential to improve siRNA properties were introduced and systematically validated in MDCK cells against the virus. A substantial and maximum inhibitory effect was achieved at concentrations as low as 8 nM. The inhibition of viral replication reached approximately 90% for the best siRNA variants. Additionally, selected siRNAs were compared with antisense oligonucleotides targeting the same regions; this revealed that effectiveness depends on both the target accessibility and oligonucleotide antiviral strategy. Our new approach of target-site preselection based on segment 5 (+)RNA secondary structure led to effective viral inhibition and a better understanding of the impact of RNA structural motifs on the influenza replication cycle.

Changes in RNA secondary structure affect NS1 protein expression during early stage influenza virus infection

RNA secondary structures play a key role in splicing, gene expression, microRNA biogenesis, RNA editing, and other biological processes. The importance of RNA structures has been demonstrated in the life cycle of RNA-containing viruses, including the influenza virus. At least two regions of conserved secondary structure in NS segment (+) RNA are predicted to vary among influenza virus strains with respect to thermodynamic stability; both fall in the NS1 open reading frame. The NS1 protein is involved in multiple virus-host interaction processes, and its main function is to inhibit the cellular immune response to viral infection. Using a reverse genetics approach, four influenza virus strains were constructed featuring mutations that have different effects on RNA secondary structure. Growth curve experiments and ELISA data show that, at least in the first viral replication cycle, mutations G123A and A132G affecting RNA structure in the (82-148) NS RNA region influence NS1 protein expression.

Mapping the RNA structural landscape of viral genomes

Functional RNA structures are prevalent in viral genomes, and have been shown to play roles in almost every aspect of their biology. However, the majority of viral RNA remains structurally uncharacterized. This is likely to remain true as the cost of sequencing decreases much faster than the cost of structural characterizations. Because of this, there is a need for rapid, inexpensive methods to highlight regions of viral RNA which are ideal candidates for structure-function analyses. The ScanFold method was developed as a single sequence alternative to traditional RNA structural motif discovery pipelines, which rely heavily on well curated sequence alignments to identify conserved RNA structures. ScanFold focuses on identifying (based on their more stable than expected folding energies) the most likely functional structures encoded within a single large RNA sequence, while allowing predicted motifs to be tested for evidence of structural conservation later. Decoupling these processes can be a benefit to researchers studying viruses lacking the ideal phylogenetic depth to yield evidence of structural conservation. Here, we demonstrate how the most significant ScanFold predicted structures correspond to higher base pairing probabilities, SHAPE reactivities, and predict known functional structures within the ZIKV and HIV-1 genomes with accuracy. Best practices and examples are also shown to aid users in utilizing ScanFold for their own systems of interest. ScanFold is available as a Webserver (https://mosslabtools.bb.iastate.edu/scanfold) or can be downloaded (https://github.com/moss-lab/ScanFold) and run locally.

Allele-specific nonstationarity in evolution of influenza A virus surface proteins

Influenza A virus (IAV) is a major public health problem and a pandemic threat. Its evolution is largely driven by diversifying positive selection so that relative fitness of different amino acid variants changes with time due to changes in herd immunity or genomic context, and novel amino acid variants attain fitness advantage. Here, we hypothesize that diversifying selection also has another manifestation: the fitness associated with a particular amino acid variant should decline with time since its origin, as the herd immunity adapts to it. By tracing the evolution of antigenic sites at IAV surface proteins, we show that an amino acid variant becomes progressively more likely to become replaced by another variant with time since its origin-a phenomenon we call "senescence." Senescence is particularly pronounced at experimentally validated antigenic sites, implying that it is largely driven by host immunity. By contrast, at internal sites, existing variants become more favorable with time, probably due to arising contingent mutations at other epistatically interacting sites. Our findings reveal a previously undescribed facet of adaptive evolution and suggest approaches for prediction of evolutionary dynamics of pathogens.

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.

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.

In vivo analysis of influenza A mRNA secondary structures identifies critical regulatory motifs

The influenza A virus (IAV) is a continuous health threat to humans as well as animals due to its recurring epidemics and pandemics. The IAV genome is segmented and the eight negative-sense viral RNAs (vRNAs) are transcribed into positive sense complementary RNAs (cRNAs) and viral messenger RNAs (mRNAs) inside infected host cells. A role for the secondary structure of IAV mRNAs has been hypothesized and debated for many years, but knowledge on the structure mRNAs adopt in vivo is currently missing. Here we solve, for the first time, the in vivo secondary structure of IAV mRNAs in living infected cells. We demonstrate that, compared to the in vitro refolded structure, in vivo IAV mRNAs are less structured but exhibit specific locally stable elements. Moreover, we show that the targeted disruption of these high-confidence structured domains results in an extraordinary attenuation of IAV replicative capacity. Collectively, our data provide the first comprehensive map of the in vivo structural landscape of IAV mRNAs, hence providing the means for the development of new RNA-targeted antivirals.

The RNA encoding the microtubule-associated protein tau has extensive structure that affects its biology

Tauopathies are neurodegenerative diseases that affect millions of people worldwide including those with Alzheimer’s disease. While many efforts have focused on understanding the role of tau protein in neurodegeneration, there has been little done to systematically analyze and study the structures within tau’s encoding RNA and their connection to disease pathology. Knowledge of RNA structure can provide insights into disease mechanisms and how to affect protein production for therapeutic benefit. Using computational methods based on thermodynamic stability and evolutionary conservation, we identified structures throughout the tau pre-mRNA, especially at exon-intron junctions and within the 5′ and 3′ untranslated regions (UTRs). In particular, structures were identified at twenty exon-intron junctions. The 5′ UTR contains one structured region, which lies within a known internal ribosome entry site. The 3′ UTR contains eight structured regions, including one that contains a polyadenylation signal. A series of functional experiments were carried out to assess the effects of mutations associated with mis-regulation of alternative splicing of exon 10 and to identify regions of the 3′ UTR that contain cis-regulatory elements. These studies defined novel structural regions within the mRNA that affect stability and pre-mRNA splicing and may lead to new therapeutic targets for treating tau-associated diseases.

Endless Forms: Within-Host Variation in the Structure of the West Nile Virus RNA Genome during Serial Passage in Bird Hosts

The enzymes that copy RNA genomes lack proofreading, and viruses that possess RNA genomes, such as West Nile virus, rapidly diversify into swarms of mutant lineages within a host. Intrahost variation of the primary genomic sequence of RNA viruses has been studied extensively because the extent of this variation shapes key virus phenotypes. However, RNA genomes also form complex secondary structures based on within-genome nucleotide complementarity, which are critical regulators of the cyclization of the virus genome that is necessary for efficient replication and translation. We sought to characterize variation in these secondary structures within populations of West Nile virus during serial passage in three bird species. Our study indicates that the intrahost population of West Nile virus is a diverse assortment of RNA secondary structures that should be considered in future analyses of intrahost viral diversity, but some regions that are critical for genome cyclization are conserved within hosts. Besides potential impacts on viral replication, structural diversity can influence the efficacy of small RNA antiviral therapies.

RNA viruses are infamous for their high rates of mutation, which produce swarms of genetic variants within individual hosts. To date, analyses of intrahost genetic diversity have focused on the primary genome sequence. However, virus phenotypes are shaped not only by primary sequence but also by the secondary structures into which this sequence folds. Such structures enable viral replication, translation, and binding of small RNAs, yet within-host variation at the structural level has not been adequately explored. We characterized the structural diversity of the 5′ untranslated region (UTR) of populations of West Nile virus (WNV) that had been subject to five serial passages in triplicate in each of three bird species. Viral genomes were sampled from host serum samples at each passage (n = 45 populations) and subjected to next-generation sequencing. For populations derived from passages 1, 3, and 5 (n = 9 populations), we predicted the impact of each mutation occurring at a frequency of ≥1% on the secondary structure of the 5′ UTR. As expected, mutations in double-stranded (DS) regions of the 5′ UTR stem structures caused structural changes of significantly greater magnitude than did mutations in single-stranded (SS) regions. Despite the greater impact of mutations in DS regions, mutations in DS and SS regions occurred at similar frequencies, with no evidence of enhanced selection against mutation in DS regions. In contrast, mutations in two regions that mediate genome cyclization and thereby regulate replication and translation, the 5′ cyclization sequence and the UAR flanking stem (UFS), were suppressed in all three hosts.

IMPORTANCE The enzymes that copy RNA genomes lack proofreading, and viruses that possess RNA genomes, such as West Nile virus, rapidly diversify into swarms of mutant lineages within a host. Intrahost variation of the primary genomic sequence of RNA viruses has been studied extensively because the extent of this variation shapes key virus phenotypes. However, RNA genomes also form complex secondary structures based on within-genome nucleotide complementarity, which are critical regulators of the cyclization of the virus genome that is necessary for efficient replication and translation. We sought to characterize variation in these secondary structures within populations of West Nile virus during serial passage in three bird species. Our study indicates that the intrahost population of West Nile virus is a diverse assortment of RNA secondary structures that should be considered in future analyses of intrahost viral diversity, but some regions that are critical for genome cyclization are conserved within hosts. Besides potential impacts on viral replication, structural diversity can influence the efficacy of small RNA antiviral therapies.

Regulation of influenza A virus mRNA splicing by CLK1

Influenza A virus carries eight negative single-stranded RNAs and uses spliced mRNAs to increase the number of proteins produced from them. Several genome-wide screens for essential host factors for influenza A virus replication revealed a necessity for splicing and splicing-related factors, including Cdc-like kinase 1 (CLK1). This CLK family kinase plays a role in alternative splicing regulation through phosphorylation of serine-arginine rich (SR) proteins. To examine the influence that modulation of splicing regulation has on influenza infection, we analyzed the effect of CLK1 knockdown and inhibition. CLK1 knockdown in A549 cells reduced influenza A/WSN/33 virus replication and increased the level of splicing of segment 7, which encodes the viral M1 and M2 proteins. CLK1-/- mice infected with influenza A/England/195/2009 (H1N1pdm09) virus supported lower levels of virus replication than wild-type mice. Screening of newly developed CLK inhibitors revealed several compounds that have an effect on the level of splicing of influenza A gene segment M in different models and decrease influenza A/WSN/33 virus replication in A549 cells. The promising inhibitor KH-CB19, an indole-based enaminonitrile with unique binding mode for CLK1, and its even more selective analogue NIH39 showed high specificity towards CLK1 and had a similar effect on influenza mRNA splicing regulation. Taken together, our findings indicate that targeting host factors that regulate splicing of influenza mRNAs may represent a novel therapeutic approach.

Conserved Secondary Structures in Viral mRNAs

RNA secondary structure in untranslated and protein coding regions has been shown to play an important role in regulatory processes and the viral replication cycle. While structures in non-coding regions have been investigated extensively, a thorough overview of the structural repertoire of protein coding mRNAs, especially for viruses, is lacking. Secondary structure prediction of large molecules, such as long mRNAs remains a challenging task, as the contingent of structures a sequence can theoretically fold into grows exponentially with sequence length. We applied a structure prediction pipeline to Viral Orthologous Groups that first identifies the local boundaries of potentially structured regions and subsequently predicts their functional importance. Using this procedure, the orthologous groups were split into structurally homogenous subgroups, which we call subVOGs. This is the first compilation of potentially functional conserved RNA structures in viral coding regions, covering the complete RefSeq viral database. We were able to recover structural elements from previous studies and discovered a variety of novel structured regions. The subVOGs are available through our web resource RNASIV (RNA structure in viruses; http://rnasiv.bio.wzw.tum.de).

Computational approaches for the discovery of splicing regulatory RNA structures

Global RNA structure and local functional motifs mediate interactions important in determining the rates and patterns of mRNA splicing. In this review, we overview approaches for the computational prediction of RNA secondary structure with a special emphasis on the discovery of motifs important to RNA splicing. The process of identifying and modeling potential splicing regulatory structures is illustrated using a recently-developed approach for RNA structural motif discovery, the ScanFold pipeline, which is applied to the identification of a known splicing regulatory structure in influenza virus.

Precise small-molecule cleavage of an r(CUG) repeat expansion in a myotonic dystrophy mouse model

Development of small-molecule lead medicines that potently and specifically modulate RNA function is challenging. We designed a small molecule that cleaves r(CUG)^exp, the RNA repeat expansion that causes myotonic dystrophy type 1. In cells and in an animal model, the small-molecule cleaver specifically recognizes the 3-dimensional structure of r(CUG)^exp, cleaving it more selectively among transcripts containing short, nonpathogenic r(CUG) repeats than an oligonucleotide that recognizes RNA sequence via Watson-Crick base pairing. The small molecule broadly relieves disease-associated phenotype in a mouse model. Thus, small molecules that recognize and cleave RNA structures should be considered a therapeutic strategy for targeting RNA in vivo.

Myotonic dystrophy type 1 (DM1) is an incurable neuromuscular disorder caused by an expanded CTG repeat that is transcribed into r(CUG)^exp. The RNA repeat expansion sequesters regulatory proteins such as Muscleblind-like protein 1 (MBNL1), which causes pre-mRNA splicing defects. The disease-causing r(CUG)^exp has been targeted by antisense oligonucleotides, CRISPR-based approaches, and RNA-targeting small molecules. Herein, we describe a designer small molecule, Cugamycin, that recognizes the structure of r(CUG)^exp and cleaves it in both DM1 patient-derived myotubes and a DM1 mouse model, leaving short repeats of r(CUG) untouched. In contrast, oligonucleotides that recognize r(CUG) sequence rather than structure cleave both long and short r(CUG)-containing transcripts. Transcriptomic, histological, and phenotypic studies demonstrate that Cugamycin broadly and specifically relieves DM1-associated defects in vivo without detectable off-targets. Thus, small molecules that bind and cleave RNA have utility as lead chemical probes and medicines and can selectively target disease-causing RNA structures to broadly improve defects in preclinical animal models.

Sequence- And Structure-Specific Probing of RNAs by Short Nucleobase-Modified dsRNA-Binding PNAs Incorporating a Fluorescent Light-up Uracil Analog

RNAs are emerging as important biomarkers and therapeutic targets. The strategy of directly targeting double-stranded RNA (dsRNA) by triplex-formation is relatively underexplored mainly due to the weak binding at physiological conditions for the traditional triplex-forming oligonucleotides (TFOs). Compared to DNA and RNA, peptide nucleic acids (PNAs) are chemically stable and have a neutral peptide-like backbone, and thus, they show significantly enhanced binding to natural nucleic acids. We have successfully developed nucleobase-modified dsRNA-binding PNAs (dbPNAs) to facilitate structure-specific and selective recognition of dsRNA over single-stranded RNA (ssRNA) and dsDNA regions at near-physiological conditions. The triplex formation strategy facilitates the targeting of not only the sequence but also the secondary structure of RNA. Here, we report the development of novel dbPNA-based fluorescent light-up probes through the incorporation of A-U pair-recognizing 5-benzothiophene uracil (^btU). The incorporation of ^btU into dbPNAs does not affect the binding affinity toward dsRNAs significantly, in most cases, as evidenced by our nondenaturing gel shift assay data. The blue fluorescence emission intensity of ^btU-modified dbPNAs is sequence- and structure-specifically enhanced by dsRNAs, including the influenza viral RNA panhandle duplex and HIV-1-1 ribosomal frameshift-inducing RNA hairpin, but not ssRNAs or DNAs, at 200 mM NaCl, pH 7.5. Thus, dbPNAs incorporating ^btU-modified and other further modified fluorescent nucleobases will be useful biochemical tools for probing and detecting RNA structures, interactions, and functions.

Secondary structure of the segment 5 genomic RNA of influenza A virus and its application for designing antisense oligonucleotides

Influenza virus causes seasonal epidemics and dangerous pandemic outbreaks. It is a single stranded (-)RNA virus with a segmented genome. Eight segments of genomic viral RNA (vRNA) form the virion, which are then transcribed and replicated in host cells. The secondary structure of vRNA is an important regulator of virus biology and can be a target for finding new therapeutics. In this paper, the secondary structure of segment 5 vRNA is determined based on chemical mapping data, free energy minimization and structure-sequence conservation analysis for type A influenza. The revealed secondary structure has circular folding with a previously reported panhandle motif and distinct novel domains. Conservations of base pairs is 87% on average with many structural motifs that are highly conserved. Isoenergetic microarray mapping was used to additionally validate secondary structure and to discover regions that easy bind short oligonucleotides. Antisense oligonucleotides, which were designed based on modeled secondary structure and microarray mapping, inhibit influenza A virus proliferation in MDCK cells. The most potent oligonucleotides lowered virus titer by ~90%. These results define universal for type A structured regions that could be important for virus function, as well as new targets for antisense therapeutics.

ScanFold: an approach for genome-wide discovery of local RNA structural elements-applications to Zika virus and HIV

In addition to encoding RNA primary structures, genomes also encode RNA secondary and tertiary structures that play roles in gene regulation and, in the case of RNA viruses, genome replication. Methods for the identification of functional RNA structures in genomes typically rely on scanning analysis windows, where multiple partially-overlapping windows are used to predict RNA structures and folding metrics to deduce regions likely to form functional structure. Separate structural models are produced for each window, where the step size can greatly affect the returned model. This makes deducing unique local structures challenging, as the same nucleotides in each window can be alternatively base paired. We are presenting here a new approach where all base pairs from analysis windows are considered and weighted by favorable folding. This results in unique base pairing throughout the genome and the generation of local regions/structures that can be ranked by their propensity to form unusually thermodynamically stable folds. We applied this approach to the Zika virus (ZIKV) and HIV-1 genomes. ZIKV is linked to a variety of neurological ailments including microcephaly and Guillain-Barré syndrome and its (+)-sense RNA genome encodes two, previously described, functionally essential structured RNA regions. HIV, the cause of AIDS, contains multiple functional RNA motifs in its genome, which have been extensively studied. Our approach is able to successfully identify and model the structures of known functional motifs in both viruses, while also finding additional regions likely to form functional structures. All data have been archived at the RNAStructuromeDB (www.structurome.bb.iastate.edu), a repository of RNA folding data for humans and their pathogens.

The ensemble diversity of non-coding RNA structure is lower than random sequence

In addition to energetically optimal structures, RNAs can fold into near energy suboptimal conformations that may be populated and play functional roles. The diversity of this structural ensemble can be estimated using a metric derived from the calculated RNA partition function: the ensemble diversity. In this report, 10 classes of functional RNAs were analyzed: the 5.8S and 5S rRNAs, ribozyme, RNase P, snoRNA, snRNA, SRP RNA, tmRNA, Vault RNA and Y RNA. Representative sequences from each class were mutagenized in two ways: firstly, all possible point mutations were generated and secondly, wild type sequences were randomized to generate multiple scrambled mutants. Compared to the mutants, the native RNA ensemble diversity was predicted to be lower. This finding held true when all available sequences (378,455 sequences) for each RNA class (archived in the RNAcentral database) were analyzed. This suggests that a compact structural ensemble is an evolved characteristic of functional RNAs.

Directed Evolution of an Influenza Reporter Virus To Restore Replication and Virulence and Enhance Noninvasive Bioluminescence Imaging in Mice

Reporter viruses provide a powerful tool to study infection, yet incorporating a nonessential gene often results in virus attenuation and genetic instability. Here, we used directed evolution of a luciferase-expressing pandemic H1N1 (pH1N1) 2009 influenza A virus in mice to restore replication kinetics and virulence, increase the bioluminescence signal, and maintain reporter gene expression. An unadapted pH1N1 virus with NanoLuc luciferase inserted into the 5' end of the PA gene segment grew to titers 10-fold less than those of the wild type in MDCK cells and in DBA/2 mice and was less virulent. For 12 rounds, we propagated DBA/2 lung samples with the highest bioluminescence-to-titer ratios. Every three rounds, we compared in vivo replication, weight loss, mortality, and bioluminescence. Mouse-adapted virus after 9 rounds (MA-9) had the highest relative bioluminescence signal and had wild-type-like fitness and virulence in DBA/2 mice. Using reverse genetics, we discovered fitness was restored in virus rPB2-MA9/PA-D479N by a combination of PA-D479N and PB2-E158G amino acid mutations and PB2 noncoding mutations C1161T and C1977T. rPB2-MA9/PA-D479N has increased mRNA transcription, which helps restore wild-type-like phenotypes in DBA/2 and BALB/c mice. Overall, the results demonstrate that directed evolution that maximizes foreign-gene expression while maintaining genetic stability is an effective method to restore wild-type-like in vivo fitness of a reporter virus. Virus rPB2-MA9/PA-D479N is expected to be a useful tool for noninvasive imaging of pH1N1 influenza virus infection and clearance while analyzing virus-host interactions and developing new therapeutics and vaccines.IMPORTANCE Influenza viruses contribute to 290,000 to 650,000 deaths globally each year. Infection is studied in mice to learn how the virus causes sickness and to develop new drugs and vaccines. During experiments, scientists have needed to euthanize groups of mice at different times to measure the amount of infectious virus in mouse tissues. By inserting a foreign gene that causes infected cells to light up, scientists could see infection spread in living mice. Unfortunately, adding an extra gene not needed by the virus slowed it down and made it weaker. Here, we used a new strategy to restore the fitness and lethality of an influenza reporter virus; we adapted it to mouse lungs and selected for variants that had the greatest light signal. The adapted virus can be used to study influenza virus infection, immunology, and disease in living mice. The strategy can also be used to adapt other viruses.

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.

Nucleotide resolution mapping of influenza A virus nucleoprotein-RNA interactions reveals RNA features required for replication

Influenza A virus nucleoprotein (NP) association with viral RNA (vRNA) is essential for packaging, but the pattern of NP binding to vRNA is unclear. Here we applied photoactivatable ribonucleoside enhanced cross-linking and immunoprecipitation (PAR-CLIP) to assess the native-state of NP-vRNA interactions in infected human cells. NP binds short fragments of RNA (~12 nucleotides) non-uniformly and without apparent sequence specificity. Moreover, NP binding is reduced at specific locations within the viral genome, including regions previously identified as required for viral genome segment packaging. Synonymous mutations designed to alter the predicted RNA structures in these low-NP-binding regions impact genome packaging and result in virus attenuation, whereas control mutations or mutagenesis of NP-bound regions have no effect. Finally, we demonstrate that the sequence conservation of low-NP-binding regions is required in multiple genome segments for propagation of diverse mammalian and avian IAV in host cells.

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.

RNAStructuromeDB: A genome-wide database for RNA structural inference

RNA plays important roles in almost every aspect of biology, and every aspect of RNA biology is influenced by its folding. This is a particularly important consideration in the era of high-throughput sequencing, when the discovery of novel transcripts far outpaces our knowledge of their functions. To gain a comprehensive picture of biology requires a structural framework for making functional inferences on RNA. To this end we have developed the RNA Structurome Database ( https://structurome.bb.iastate.edu ), a comprehensive repository of RNA secondary structural information that spans the entire human genome. Here, we compile folding information for every base pair of the genome that may be transcribed: coding, noncoding, and intergenic regions, as well as repetitive elements, telomeres, etc. This was done by fragmenting the GRCh38 reference genome into 154,414,320 overlapping sequence fragments and, for each fragment, calculating a set of metrics based on the sequence's folding properties. These data will facilitate a wide array of investigations: e.g. discovery of structured regulatory elements in differential gene expression data or noncoding RNA discovery, as well as allow genome-scale analyses of RNA folding.

Influenza virus segment 5 (+)RNA - secondary structure and new targets for antiviral strategies

Influenza A virus is a threat for humans due to seasonal epidemics and occasional pandemics. This virus can generate new strains that are dangerous through nucleotide/amino acid changes or through segmental recombination of the viral RNA genome. It is important to gain wider knowledge about influenza virus RNA to create new strategies for drugs that will inhibit its spread. Here, we present the experimentally determined secondary structure of the influenza segment 5 (+)RNA. Two RNAs were studied: the full-length segment 5 (+)RNA and a shorter construct containing only the coding region. Chemical mapping data combined with thermodynamic energy minimization were used in secondary structure prediction. Sequence/structure analysis showed that the determined secondary structure of segment 5 (+)RNA is mostly conserved between influenza virus type A strains. Microarray mapping and RNase H cleavage identified accessible sites for oligonucleotides in the revealed secondary structure of segment 5 (+)RNA. Antisense oligonucleotides were designed based on the secondary structure model and tested against influenza virus in cell culture. Inhibition of influenza virus proliferation was noticed, identifying good targets for antisense strategies. Effective target sites fall within two domains, which are conserved in sequence/structure indicating their importance to the virus.

Visualization of conformational variability in the domains of long single-stranded RNA molecules

We demonstrate an application of atomic force microscopy (AFM) for the structural analysis of long single-stranded RNA (>1 kb), focusing on 28S ribosomal RNA (rRNA). Generally, optimization of the conditions required to obtain three-dimensional (3D) structures of long RNA molecules is a challenging or nearly impossible process. In this study, we overcome these limitations by developing a method using AFM imaging combined with automated, MATLAB-based image analysis algorithms for extracting information about the domain organization of single RNA molecules. We examined the 5 kb human 28S rRNA since it is the largest RNA molecule for which a 3D structure is available. As a proof of concept, we determined a domain structure that is in accordance with previously described secondary structural models. Importantly, we identified four additional small (200–300 nt), previously unreported domains present in these molecules. Moreover, the single-molecule nature of our method enabled us to report on the relative conformational variability of each domain structure identified, and inter-domain associations within subsets of molecules leading to molecular compaction, which may shed light on the process of how these molecules fold into the final tertiary structure.

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.

Nuclear Magnetic Resonance Structure of an 8 × 8 Nucleotide RNA Internal Loop Flanked on Each Side by Three Watson-Crick Pairs and Comparison to Three-Dimensional Predictions

The prediction of RNA three-dimensional structure from sequence alone has been a long-standing goal. High-resolution, experimentally determined structures of simple noncanonical pairings and motifs are critical to the development of prediction programs. Here, we present the nuclear magnetic resonance structure of the (5'CCAGAAACGGAUGGA)₂ duplex, which contains an 8 × 8 nucleotide internal loop flanked by three Watson-Crick pairs on each side. The loop is comprised of a central 5'AC/3'CA nearest neighbor flanked by two 3RRs motifs, a known stable motif consisting of three consecutive sheared GA pairs. Hydrogen bonding patterns between base pairs in the loop, the all-atom root-mean-square deviation for the loop, and the deformation index were used to compare the structure to automated predictions by MC-sym, RNA FARFAR, and RNAComposer.

Structure and Dynamics of RNA Repeat Expansions That Cause Huntington's Disease and Myotonic Dystrophy Type 1

RNA repeat expansions cause a host of incurable, genetically defined diseases. The most common class of RNA repeats consists of trinucleotide repeats. These long, repeating transcripts fold into hairpins containing 1 × 1 internal loops that can mediate disease via a variety of mechanism(s) in which RNA is the central player. Two of these disorders are Huntington's disease and myotonic dystrophy type 1, which are caused by r(CAG) and r(CUG) repeats, respectively. We report the structures of two RNA constructs containing three copies of a r(CAG) [r(3×CAG)] or r(CUG) [r(3×CUG)] motif that were modeled with nuclear magnetic resonance spectroscopy and simulated annealing with restrained molecular dynamics. The 1 × 1 internal loops of r(3×CAG) are stabilized by one-hydrogen bond (cis Watson-Crick/Watson-Crick) AA pairs, while those of r(3×CUG) prefer one- or two-hydrogen bond (cis Watson-Crick/Watson-Crick) UU pairs. Assigned chemical shifts for the residues depended on the identity of neighbors or next nearest neighbors. Additional insights into the dynamics of these RNA constructs were gained by molecular dynamics simulations and a discrete path sampling method. Results indicate that the global structures of the RNA are A-form and that the loop regions are dynamic. The results will be useful for understanding the dynamic trajectory of these RNA repeats but also may aid in the development of therapeutics.

The effective rate of influenza reassortment is limited during human infection

We characterise the evolutionary dynamics of influenza infection described by viral sequence data collected from two challenge studies conducted in human hosts. Viral sequence data were collected at regular intervals from infected hosts. Changes in the sequence data observed across time show that the within-host evolution of the virus was driven by the reversion of variants acquired during previous passaging of the virus. Treatment of some patients with oseltamivir on the first day of infection did not lead to the emergence of drug resistance variants in patients. Using an evolutionary model, we inferred the effective rate of reassortment between viral segments, measuring the extent to which randomly chosen viruses within the host exchange genetic material. We find strong evidence that the rate of effective reassortment is low, such that genetic associations between polymorphic loci in different segments are preserved during the course of an infection in a manner not compatible with epistasis. Combining our evidence with that of previous studies we suggest that spatial heterogeneity in the viral population may reduce the extent to which reassortment is observed. Our results do not contradict previous findings of high rates of viral reassortment in vitro and in small animal studies, but indicate that in human hosts the effective rate of reassortment may be substantially more limited.

The influenza virus is an important cause of disease in the human population. During the course of an infection the virus can evolve rapidly. An important mechanism of viral evolution is reassortment, whereby different segments of the influenza genome are shuffled with other segments, producing new viral combinations. Here we study natural selection and reassortment during the course of infections occurring in human hosts. Examining viral genome sequence data from these infections, we note that genetic variants that were acquired during the growth of viruses in culture are selected against in the human host. In addition, we find evidence that the effective rate of reassortment is low. We suggest that the spatial separation between viruses in different parts of the host airway may limit the extent to which genetically distinct segments reassort with one another. Within the global population of influenza viruses, reassortment remains an important factor. However, reassortment is not so rapid as to exclude the possibility of interactions between genome segments affecting the course of influenza evolution during a single infection.

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

The influenza A virus genome consists of eight RNA segments. RNA structures within these segments and complementary (cRNA) and protein-coding mRNAs may play a role in virus replication. Here, conserved putative secondary structures that impose significant evolutionary constraints on the gene segment encoding the surface glycoprotein hemagglutinin (HA) were investigated using available sequence data on tens of thousands of virus strains. Structural constraints were identified by analysis of covariations of nucleotides suggested to be paired by structure prediction algorithms. The significance of covariations was estimated by mutual information calculations and tracing multiple covariation events during virus evolution. Covariation patterns demonstrated that structured domains in HA RNAs were mostly subtype-specific, whereas some structures were conserved in several subtypes. The influence of RNA folding on virus replication was studied by plaque assays of mutant viruses with disrupted structures. The results suggest that over the whole length of the HA segment there are local structured domains which contribute to the virus fitness but individually are not essential for the virus. Existence of subtype-specific structured regions in the segments of the influenza A virus genome is apparently an important factor in virus evolution and reassortment of its genes.

Computational and molecular analysis of conserved influenza A virus RNA secondary structures involved in infectious virion production

As well as encoding viral proteins, genomes of RNA viruses harbor secondary and tertiary RNA structures that have been associated with functions essential for successful replication and propagation. Here, we identified stem-loop structures that are extremely conserved among 1,884 M segment sequences of influenza A virus (IAV) strains from various subtypes and host species using computational and evolutionary methods. These structures were predicted within the 3' and 5' ends of the coding regions of M1 and M2, respectively, where packaging signals have been previously proposed to exist. These signals are thought to be required for the incorporation of a single copy of 8 different negative-strand RNA segments (vRNAs) into an IAV particle. To directly test the functionality of conserved stem-loop structures, we undertook reverse genetic experiments to introduce synonymous mutations designed to disrupt secondary structures predicted at 3 locations and found them to attenuate infectivity of recombinant virus. In one mutant, predicted to disrupt stem loop structure at nucleotide positions 219-240, attenuation was more evident at increased temperature and was accompanied by an increase in the production of defective virus particles. Our results suggest that the conserved secondary structures predicted in the M segment are involved in the production of infectious viral particles during IAV replication.

Mutations Designed by Ensemble Defect to Misfold Conserved RNA Structures of Influenza A Segments 7 and 8 Affect Splicing and Attenuate Viral Replication in Cell Culture

Influenza A virus is a significant public health threat, but little is understood about the viral RNA structure and function. Current vaccines and therapeutic options to control influenza A virus infections are mostly protein-centric and of limited effectiveness. Here, we report using an ensemble defect approach to design mutations to misfold regions of conserved mRNA structures in influenza A virus segments 7 and 8. Influenza A mutant viruses inhibit pre-mRNA splicing and attenuate viral replication in cell culture, thus providing evidence for functions of the targeted regions. Targeting these influenza A viral RNA regions provides new possibilities for designing vaccines and therapeutics against this important human respiratory pathogen. The results also demonstrate that the ensemble defect approach is an efficient way to test for function of RNA sequences.

The influenza A virus NS genome segment displays lineage-specific patterns in predicted RNA secondary structure

BACKGROUND

Influenza A virus (IAV) is a segmented negative-sense RNA virus that causes seasonal epidemics and periodic pandemics in humans. Two regions (nucleotide positions 82-148 and 497-564) in the positive-sense RNA of the NS segment fold into a multi-branch loop or hairpin structures.

RESULTS

We studied 25,384 NS segment positive-sense RNA unique sequences of human and non-human IAVs in order to predict secondary RNA structures of the 82-148 and 497-564 regions using RNAfold software, and determined their host- and lineage-specific distributions. Hairpins prevailed in avian and avian-origin human IAVs, including H1N1pdm1918 and H5N1. In human and swine IAV hairpins distribution varied between evolutionary lineages.

CONCLUSIONS

These results suggest a possible functional role for these RNA secondary structures and the need for experimental evaluation of these structures in the influenza life cycle.

Elevation of CpG frequencies in influenza A genome attenuates pathogenicity but enhances host response to infection

Previously, we demonstrated that frequencies of CpG and UpA dinucleotides profoundly influence the replication ability of echovirus 7 (Tulloch et al., 2014). Here, we show that that influenza A virus (IAV) with maximised frequencies of these dinucleotides in segment 5 showed comparable attenuation in cell culture compared to unmodified virus and a permuted control (CDLR). Attenuation was also manifested in vivo, with 10-100 fold reduced viral loads in lungs of mice infected with 200PFU of CpG-high and UpA-high mutants. However, both induced powerful inflammatory cytokine and adaptive (T cell and neutralising antibody) responses disproportionate to their replication. CpG-high infected mice also showed markedly reduced clinical severity, minimal weight loss and reduced immmunopathology in lung, yet sterilising immunity to lethal dose WT challenge was achieved after low dose (20PFU) pre-immunisation with this mutant. Increasing CpG dinucleotide frequencies represents a generic and potentially highly effective method for generating safe, highly immunoreactive vaccines.

Self-Folding of Naked Segment 8 Genomic RNA of Influenza A Virus

Influenza A is a negative sense RNA virus that kills hundreds of thousands of humans each year. Base pairing in RNA is very favorable, but possibilities for RNA secondary structure of the influenza genomic RNA have not been investigated. This work presents the first experimentally-derived exploration of potential secondary structure in an influenza A naked (protein-free) genomic segment. Favorable folding regions are revealed by in vitro chemical structure mapping, thermodynamics, bioinformatics, and binding to isoenergetic microarrays of an entire natural sequence of the 875 nt segment 8 vRNA and of a smaller fragment. Segment 8 has thermodynamically stable and evolutionarily conserved RNA structure and encodes essential viral proteins NEP and NS1. This suggests that vRNA self-folding may generate helixes and loops that are important at one or more stages of the influenza life cycle.