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.

Conserved structured domains in plant non-coding RNA enod40, their evolution and recruitment of sequences from transposable elements.

Plant long noncoding RNA enod40 is involved in the regulation of symbiotic associations with bacteria, in particular, in nitrogen-fixing root nodules of legumes, and with fungi in phosphate-acquiring arbuscular mycorrhizae formed by various plants. The presence of enod40 genes in plants that do not form such symbioses indicates its other roles in cell physiology. The molecular mechanisms of enod40 RNA function are poorly understood. Enod40 RNAs form several structured domains, conserved to different extents. Due to relatively low sequence similarity, identification of enod40 sequences in plant genomes is not straightforward, and many enod40 genes remain unannotated even in complete genomes. Here, we used comparative structure analysis and sequence similarity searches in order to locate enod40 genes and determine enod40 RNA structures in nitrogen-fixing clade plants and in grasses. The structures combine conserved features with considerable diversity of structural elements, including insertions of structured domain modules originating from transposable elements. Remarkably, these insertions contain sequences similar to tandem repeats and several stem-loops are homologous to microRNA precursors.

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.

All genera of Flaviviridae host a conserved Xrn1-resistant RNA motif.

After infection by flaviviruses like Zika and West Nile virus, eukaryotic hosts employ the well-conserved endoribonuclease Xrn1 to degrade the viral genomic RNA. Within the 3' untranslated regions, this enzyme encounters intricate Xrn1-resistant structures. This results in the accumulation of subgenomic flaviviral RNAs, an event that improves viral growth and aggravates viral pathogenicity. Xrn1-resistant RNAs have been established throughout the flaviviral genus, but not yet throughout the entire Flaviviridae family. In this work, we use previously determined characteristics of these structures to identify homologous sequences in many members of the genera pegivirus, hepacivirus and pestivirus. We used structural alignment and mutational analyses to establish that these sequences indeed represent Xrn1-resistant RNA and that they employ the general features of the flaviviral xrRNAs, consisting of a double pseudoknot formed by five base-paired regions stitched together by a crucial triple base interaction. Furthermore, we demonstrate that the pestivirus Bungowannah virus produces subgenomic RNA in vivo. Altogether, these results indicate that viruses make use of a universal Xrn1-resistant RNA throughout the Flaviviridae family.

A database of flavivirus RNA structures with a search algorithm for pseudoknots and triple base interactions.

MOTIVATION: The Flavivirus genus includes several important pathogens, such as Zika, dengue and yellow fever virus. Flavivirus RNA genomes contain a number of functionally important structures in their 3' untranslated regions (3'UTRs). Due to the diversity of sequences and topologies of these structures, their identification is often difficult. In contrast, predictions of such structures are important for understanding of flavivirus replication cycles and development of antiviral strategies. RESULTS: We have developed an algorithm for structured pattern search in RNA sequences, including secondary structures, pseudoknots and triple base interactions. Using the data on known conserved flavivirus 3'UTR structures, we constructed structural descriptors which covered the diversity of patterns in these motifs. The descriptors and the search algorithm were used for the construction of a database of flavivirus 3'UTR structures. Validating this approach, we identified a number of domains matching a general pattern of exoribonuclease Xrn1-resistant RNAs in the growing group of insect-specific flaviviruses. AVAILABILITY AND IMPLEMENTATION: The Leiden Flavivirus RNA Structure Database is available at https://rna.liacs.nl. The search algorithm is available at https://github.com/LeidenRNA/SRHS. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Xrn1-resistant RNA structures are well-conserved within the genus flavivirus.

Subgenomic RNAs are produced by several RNA viruses through incomplete degradation of their genomic RNA by the exoribonuclease Xrn1, and have been shown to be essential for viral growth and pathogenicity. Within the flavivirus genus of the Flaviviridae family, two distinct classes of Xrn1-resistant RNA motifs have been proposed; one for mosquito-borne and insect-specific flaviviruses, and one for tick-borne flaviviruses and no-known-vector flaviviruses. We investigated tick-borne and no-known-vector flavivirus Xrn1-resistant RNA motifs through systematic in vitro mutational analysis and showed that both classes actually possess very similar structural configurations, including a double pseudoknot and a base-triple at identical, conserved locations. For the no-known-vector flavivirus Modoc virus, we show that in vivo generation of subgenomic flaviviral RNA was affected by mutations targeted at nucleotides involved in the structural features of flaviviral Xrn1-resistant RNA motifs that were defined in this work. Our results suggest that throughout the genus flavivirus Xrn1-resistant RNA motifs adopt the same topologically conserved structure.

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.

Fungal metabarcoding data integration framework for the MycoDiversity DataBase (MDDB).

Fungi have crucial roles in ecosystems, and are important associates for many organisms. They are adapted to a wide variety of habitats, however their global distribution and diversity remains poorly documented. The exponential growth of DNA barcode information retrieved from the environment is assisting considerably the traditional ways for unraveling fungal diversity and detection. The raw DNA data in association to environmental descriptors of metabarcoding studies are made available in public sequence read archives. While this is potentially a valuable source of information for the investigation of Fungi across diverse environmental conditions, the annotation used to describe environment is heterogenous. Moreover, a uniform processing pipeline still needs to be applied to the available raw DNA data. Hence, a comprehensive framework to analyses these data in a large context is still lacking. We introduce the MycoDiversity DataBase, a database which includes public fungal metabarcoding data of environmental samples for the study of biodiversity patterns of Fungi. The framework we propose will contribute to our understanding of fungal biodiversity and aims to become a valuable source for large-scale analyses of patterns in space and time, in addition to assisting evolutionary and ecological research on Fungi.

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.

Structural features of an Xrn1-resistant plant virus RNA.

Xrn1 is a major 5'-3' exoribonuclease involved in the RNA metabolism of many eukaryotic species. RNA viruses have evolved ways to thwart Xrn1 in order to produce subgenomic non-coding RNA that affects the hosts RNA metabolism. The 3' untranslated region of several beny- and cucumovirus RNAs harbors a so-called 'coremin' motif that is required for Xrn1 stalling. The structural features of this motif have not been studied in detail yet. Here, by using in vitro Xrn1 degradation assays, we tested over 50 different RNA constructs based on the Beet necrotic yellow vein virus sequence to deduce putative structural features responsible for Xrn1 stalling. We demonstrated that the minimal benyvirus stalling site consists of two hairpins of 3 and 4 base pairs respectively. The 5' proximal hairpin requires a YGAD (Y = U/C, D = G/A/U) consensus loop sequence, whereas the 3' proximal hairpin loop sequence is variable. The sequence of the 10-nucleotide spacer that separates the hairpins is highly conserved and potentially involved in tertiary interactions. Similar coremin motifs were identified in plant virus isolates from other families including Betaflexiviridae, Virgaviridae, Potyviridae and Secoviridae (order of the Picornavirales). We conclude that Xrn1-stalling motifs are more widespread among RNA viruses than previously realized.

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.

Optimisations and Challenges Involved in the Creation of Various Bioluminescent and Fluorescent Influenza A Virus Strains for In Vitro and In Vivo Applications.

Bioluminescent and fluorescent influenza A viruses offer new opportunities to study influenza virus replication, tropism and pathogenesis. To date, several influenza A reporter viruses have been described. These strategies typically focused on a single reporter gene (either bioluminescent or fluorescent) in a single virus backbone. However, whilst bioluminescence is suited to in vivo imaging, fluorescent viruses are more appropriate for microscopy. Therefore, the idea l reporter virus varies depending on the experiment in question, and it is important that any reporter virus strategy can be adapted accordingly. Herein, a strategy was developed to create five different reporter viruses in a single virus backbone. Specifically, enhanced green fluorescent protein (eGFP), far-red fluorescent protein (fRFP), near-infrared fluorescent protein (iRFP), Gaussia luciferase (gLUC) and firefly luciferase (fLUC) were inserted into the PA gene segment of A/PR/8/34 (H1N1). This study provides a comprehensive characterisation of the effects of different reporter genes on influenza virus replication and reporter activity. In vivo reporter gene expression, in lung tissues, was only detected for eGFP, fRFP and gLUC expressing viruses. In vitro, the eGFP-expressing virus displayed the best reporter stability and could be used for correlative light electron microscopy (CLEM). This strategy was then used to create eGFP-expressing viruses consisting entirely of pandemic H1N1, highly pathogenic avian influenza (HPAI) H5N1 and H7N9. The HPAI H5N1 eGFP-expressing virus infected mice and reporter gene expression was detected, in lung tissues, in vivo. Thus, this study provides new tools and insights for the creation of bioluminescent and fluorescent influenza A reporter viruses.

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

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

Tospovirus ambisense genomic RNA segments use almost complete repertoire of stable tetraloops in the intergenic region.

The intergenic regions of the ambisense RNA segments of viruses from the Tospovirus genus form large extended RNA structures that regulate virus replication. Using comparative structure analysis, we show the presence of conserved alternative conformations at the apical parts of these structures. In one conformation, a branched Y-shape, the 5'-proximal hairpin arms are mostly capped by exceptionally stable tetraloop motifs. The tetraloop hairpins are folded in both virus and virus-complementary sense RNAs, and different tetraloops can functionally replace each other. Folding simulations show that the branched Y-shape structures can undergo a conformational transition to alternative extended rod-like conformations. Functional importance of both alternatives is supported by nucleotide covariations. The balanced equilibrium between alternative structures is evidenced by native gel electrophoresis of mutant RNA transcripts with shifted equilibria. The tetraloops play a role in the stability and dynamics of structures but may also be recognized by proteins involved in translation and/or replication.

Upstream start codon in segment 4 of North American H2 avian influenza A viruses.

H2N2 influenza A virus was the cause of the 1957 pandemic. Due to its constant presence in birds, the H2 subtype remains a topic of interest. In this work, comparison of H2 leader sequences of influenza A segment 4 revealed the presence of an upstream in-frame start codon in a majority of North American avian strains. This AUG is located seven codons upstream of the conventional start codon and is in a good Kozak context. In vivo experiments, using a luciferase reporter gene fused to leader sequences derived from North American avian H2 strains, support the efficient use of the upstream start codon. These results were corroborated by in vitro translation data using full-length segment 4 mRNA. Phylogenic analyses indicate that the upstream AUG, first detected in 1976, is stably nested in the North American avian lineage of H2 strains nowadays. The possible consequences of the upstream AUG are discussed.

Influenza virus RNA structure: unique and common features.

The influenza A virus genome consists of eight negative-sense RNA segments. Here we review the currently available data on structure-function relationships in influenza virus RNAs. Various ideas and hypotheses about the roles of influenza virus RNA folding in the virus replication are also discussed in relation to other viruses.

A family of non-classical pseudoknots in influenza A and B viruses.

A very conserved pseudoknot structure has been shown to fold in influenza virus RNA. The pseudoknot encompasses the 3' splice site of segment 8 RNA in both influenza A and B viruses. By sequence comparison of influenza virus strains, we derive a consensus motif that defines a novel RNA pseudoknot family. The orientation of the coaxially stacked stems in the influenza pseudoknot differs from that in classical H-pseudoknots. Apart from the size of the loops, the topology of the influenza pseudoknot resembles that of some long-range pseudoknotted conformations. A seed alignment of the influenza pseudoknot family, containing representative strain sequences together with a consensus structure description, has been submitted to the RNA families (Rfam) database.

PseudoBase++: an extension of PseudoBase for easy searching, formatting and visualization of pseudoknots.

Pseudoknots have been recognized to be an important type of RNA secondary structures responsible for many biological functions. PseudoBase, a widely used database of pseudoknot secondary structures developed at Leiden University, contains over 250 records of pseudoknots obtained in the past 25 years through crystallography, NMR, mutational experiments and sequence comparisons. To promptly address the growing analysis requests of the researchers on RNA structures and bring together information from multiple sources across the Internet to a single platform, we designed and implemented PseudoBase++, an extension of PseudoBase for easy searching, formatting and visualization of pseudoknots. PseudoBase++ (http://pseudobaseplusplus.utep.edu) maps the PseudoBase dataset into a searchable relational database including additional functionalities such as pseudoknot type. PseudoBase++ links each pseudoknot in PseudoBase to the GenBank record of the corresponding nucleotide sequence and allows scientists to automatically visualize RNA secondary structures with PseudoViewer. It also includes the capabilities of fine-grained reference searching and collecting new pseudoknot information.

Identification of conserved secondary structures and expansion segments in enod40 RNAs reveals new enod40 homologues in plants.

enod40 is a plant gene that participates in the regulation of symbiotic interaction between leguminous plants and bacteria or fungi. Furthermore, it has been suggested to play a general role in non-symbiotic plant development. Although enod40 seems to have multiple functions, being present in many land plants, the molecular mechanisms of its activity are unclear; they may be determined though, by short peptides and/or RNA structures encoded in the enod40 genes. We utilized conserved RNA structures in enod40 sequences to search nucleotide sequence databases and identified a number of new enod40 homologues in plant species that belong to known, but also, to yet unknown enod40-containing plant families. RNA secondary structure predictions and comparative sequence analysis of enod40 RNAs allowed us to determine the most conserved structural features, present in all known enod40 genes. Remarkably, the topology and evolution of one of the conserved structural domains are similar to those of the expansion segments found in structural RNAs such as rRNAs, RNase P and SRP RNAs. Surprisingly, the enod40 RNA structural elements are much more stronger conserved than the encoded peptides. This finding suggests that some general functions of enod40 gene could be determined by the encoded RNA structure, whereas short peptides may be responsible for more diverse functions found only in certain plant families.