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

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.

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

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

The RNA pseudoknots in foot-and-mouth disease virus are dispensable for genome replication, but essential for the production of infectious virus

Non-coding regions of viral RNA (vRNA) genomes are critically important in the regulation of gene expression. In particular, pseudoknot (PK) structures, which are present in a wide range of RNA molecules, have a variety of roles. The 5' untranslated region (5' UTR) of foot-and-mouth disease virus (FMDV) vRNA is considerably longer than in other viruses from the picornavirus family and consists of a number of distinctive structural motifs that includes multiple (2, 3 or 4 depending on the virus strain) putative PKs linked in tandem. The role(s) of the PKs in the FMDV infection are not fully understood. Here, using bioinformatics, sub-genomic replicons and recombinant viruses we have investigated the structural conservation and importance of the PKs in the FMDV lifecycle. Our results show that despite the conservation of two or more PKs across all FMDVs, a replicon lacking PKs was replication competent, albeit at reduced levels. Furthermore, in competition experiments, GFP FMDV replicons with less than two (0 or 1) PK structures were outcompeted by a mCherry FMDV wt replicon that had 4 PKs, whereas GFP replicons with 2 or 4 PKs were not. This apparent replicative advantage offered by the additional PKs correlates with the maintenance of at least two PKs in the genomes of FMDV field isolates. Despite a replicon lacking any PKs retaining the ability to replicate, viruses completely lacking PK were not viable and at least one PK was essential for recovery of infections virus, suggesting a role for the PKs in virion assembly. Thus, our study points to roles for the PKs in both vRNA replication and virion assembly, thereby improving understanding the molecular biology of FMDV replication and the wider roles of PK in RNA functions.

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.

Influenza A virus segments five and six can harbor artificial introns allowing expanded coding capacity

Influenza A viruses encode their genomes across eight, negative sense RNA segments. The six largest segments produce mRNA transcripts that do not generally splice; however, the two smallest segments are actively spliced to produce the essential viral proteins NEP and M2. Thus, viral utilization of RNA splicing effectively expands the viral coding capacity without increasing the number of genomic segments. As a first step towards understanding why splicing is not more broadly utilized across genomic segments, we designed and inserted an artificial intron into the normally nonsplicing NA segment. This insertion was tolerated and, although viral mRNAs were incompletely spliced, we observed only minor effects on viral fitness. To take advantage of the unspliced viral RNAs, we encoded a reporter luciferase gene in frame with the viral ORF such that when the intron was not removed the reporter protein would be produced. This approach, which we also show can be applied to the NP encoding segment and in different viral genetic backgrounds, led to high levels of reporter protein expression with minimal effects on the kinetics of viral replication or the ability to cause disease in experimentally infected animals. These data together show that the influenza viral genome is more tolerant of splicing than previously appreciated and this knowledge can be leveraged to develop viral genetic platforms with utility for biotechnology applications.

Unlike most host mRNAs, some viral mRNAs encode multiple discrete, functional proteins. One method influenza A viruses use to increase the protein products from two of their eight RNA genome segments is splicing. Splicing requires host machinery to remove part of the viral mRNA, the intron, to generate a different mRNA product. Although only certain influenza viral segments naturally splice, we were interested in whether additional segments could splice to produce multiple proteins. We inserted artificial introns harboring reporter genes into otherwise nonsplicing genomic segments of an H1N1 influenza A virus and found that this modification was well tolerated by the virus. We further demonstrated that an unrelated H3N2 influenza A virus could similarly support splicing and express a reporter protein from an artificial intron. These findings have implications for our understanding of how viruses expand their coding capacity with a limited genome. Additionally, encoding reporter proteins in spliced intronic sequences also represents a new method of generating reporter viruses requiring limited manipulation of the viral RNA.

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.

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.

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.

A new method for detecting signal regions in ordered sequences of real numbers, and application to viral genomic data

We present a fast, robust and parsimonious approach to detecting signals in an ordered sequence of numbers. Our motivation is in seeking a suitable method to take a sequence of scores corresponding to properties of positions in virus genomes, and find outlying regions of low scores. Suitable statistical methods without using complex models or making many assumptions are surprisingly lacking. We resolve this by developing a method that detects regions of low score within sequences of real numbers. The method makes no assumptions a priori about the length of such a region; it gives the explicit location of the region and scores it statistically. It does not use detailed mechanistic models so the method is fast and will be useful in a wide range of applications. We present our approach in detail, and test it on simulated sequences. We show that it is robust to a wide range of signal morphologies, and that it is able to capture multiple signals in the same sequence. Finally we apply it to viral genomic data to identify regions of evolutionary conservation within influenza and rotavirus.

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/.

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

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

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

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

Secondary structure model of the naked segment 7 influenza A virus genomic RNA

The influenza A virus (IAV) genome comprises eight negative-sense viral (v)RNA segments. The seventh segment of the genome encodes two essential viral proteins and is specifically packaged alongside the other seven vRNAs. To gain insights into the possible roles of RNA structure both within and without virions, a secondary structure model of a naked (protein-free) segment 7 vRNA (vRNA7) has been determined using chemical mapping and thermodynamic energy minimization. The proposed structure model was validated using microarray mapping, RNase H cleavage and comparative sequence analysis. Additionally, the detailed structures of three vRNA7 fragment constructs - comprising independently folded subdomains - were determined. Much of the proposed vRNA7 structure is preserved between IAV strains, suggesting their importance in the influenza replication cycle. Possible structure rearrangements, which allow or preclude long-range RNA interactions, are also proposed.

Seasonality and selective trends in viral acute respiratory tract infections

Influenza A and B, and many unrelated viruses including rhinovirus, RSV, adenovirus, metapneumovirus and coronavirus share the same seasonality, since these viral acute respiratory tract infections (vARIs) are much more common in winter than summer. Unfortunately, early investigations that used recycled “pedigree” virus strains seem to have led microbiologists to dismiss the common folk belief that vARIs often follow chilling. Today, incontrovertible evidence shows that ambient temperature dips and host chilling increase the incidence and severity of vARIs. This review considers four possible mechanisms, M1 - 4, that can explain this link: (M1) increased crowding in winter may enhance viral transmission; (M2) lower temperatures may increase the stability of virions outside the body; (M3) chilling may increase host susceptibility; (M4) lower temperatures or host chilling may activate dormant virions. There is little evidence for M1 or M2, which are incompatible with tropical observations. Epidemiological anomalies such as the repeated simultaneous arrival of vARIs over wide geographical areas, the rapid cessation of influenza epidemics, and the low attack rate of influenza within families are compatible with M4, but not M3 (in its simple form). M4 seems to be the main driver of seasonality, but M3 may also play an important role.

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

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

Microarrays for identifying binding sites and probing structure of RNAs

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

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

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

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

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

Influenza viruses and mRNA splicing: doing more with less

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

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

BACKGROUND

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

FINDINGS

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

CONCLUSIONS

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

Evidence of pervasive biologically functional secondary structures within the genomes of eukaryotic single-stranded DNA viruses

Single-stranded DNA (ssDNA) viruses have genomes that are potentially capable of forming complex secondary structures through Watson-Crick base pairing between their constituent nucleotides. A few of the structural elements formed by such base pairings are, in fact, known to have important functions during the replication of many ssDNA viruses. Unknown, however, are (i) whether numerous additional ssDNA virus genomic structural elements predicted to exist by computational DNA folding methods actually exist and (ii) whether those structures that do exist have any biological relevance. We therefore computationally inferred lists of the most evolutionarily conserved structures within a diverse selection of animal- and plant-infecting ssDNA viruses drawn from the families Circoviridae, Anelloviridae, Parvoviridae, Nanoviridae, and Geminiviridae and analyzed these for evidence of natural selection favoring the maintenance of these structures. While we find evidence that is consistent with purifying selection being stronger at nucleotide sites that are predicted to be base paired than at sites predicted to be unpaired, we also find strong associations between sites that are predicted to pair with one another and site pairs that are apparently coevolving in a complementary fashion. Collectively, these results indicate that natural selection actively preserves much of the pervasive secondary structure that is evident within eukaryote-infecting ssDNA virus genomes and, therefore, that much of this structure is biologically functional. Lastly, we provide examples of various highly conserved but completely uncharacterized structural elements that likely have important functions within some of the ssDNA virus genomes analyzed here.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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