RNA structure interactions and ribonucleoprotein processes of the influenza A virus

Abstract

In one more years, we will ‘celebrate’ an exact centenary of the Spanish flu pandemic. With the rapid evolution of the influenza virus, the possibility of novel pandemic remains ever a concern. This review covers our current knowledge of the influenza A virus: on the role of RNA in translation, replication, what is known of the expressed proteins and the protein products generated from alternative splicing, and on the role of base pairing in RNA structure. We highlight the main events associated with viral entry into the cell, the transcription and replication process, an export of the viral genetic material from the nucleus and the final release of the virus. We discuss the observed potential roles of RNA secondary structure (the RNA base-pairing arrangement) and RNA/RNA interactions in this scheme.

Introduction

The influenza virus belongs to the class known as Orthomyxoviruses [1, 2] and can be divided into three distinct subtypes (influenza A, influenza B or influenza C), depending on its nucleoproteins (NPs) and the antigen determinants of its matrix proteins. The influenza virus is a negative-sense, single-stranded RNA (ssRNA) virus that is segmented into eight (Subtypes A, B) or seven (Subtype C) separate strands. Each segment codes for at least one protein, and based on the extensive studies of the influenza A virus (IAV) in recent years, it seems likely that all subtypes express at least twice as many proteins as its number of segments [3].

Influenza is an envelope virus, where viral nucleocapsid (genome + NPs) is protected by a lipid bilayer derived from the modified host cell’s membrane [2]. The outer layer of the influenza virus is studded with the glycoproteins hemagglutinin (HA) and neuraminidase (NA). On the virion, the HA is found in the form of a homotrimer (HA-trimer), while the NA forms into a homotetramer (NA-tetramer). The HA-trimers are relatively evenly distributed on the surface of the virion in the lipid bilayer, whereas the NA-tetramers appear to localize at one end of the virion, most likely at the exit point of the budding virion [4]. Recently, the distribution of the HA-trimers into larger conglomerates was observed; however, the biological consequences of these processes are not well understood [5].

The IAVs are among the fastest evolving families of viruses [6] and are notable for their ability to rapidly change from one season to the next or even within seasons [7]. Over the course of a few generations, the virus can change sufficiently to become virulent. Like many emerging infectious diseases such as Ebola virus or severe acute respiratory syndrome (SARS) [8], IAV depends on animal reservoirs, where they can reside [9–13] and transfer between species that share the same environment [14]. As one major reservoir involves migratory birds [15], the potential to rapidly spread an infection throughout the world is of major concern and IAV is under constant watch (http://www.who.int/csr/don/archive/year/en/). During the past 100 years, there have been three major epidemics with extensive fatalities [16] and one minor one (mostly because of quick action from the WHO and health practitioners) [17, 18], the Spanish flu (1918, H1N1, 75 million fatalities worldwide), the Asian flu (1957, H2N2, 2 million fatalities), the Hong Kong flu (1968, H3N2, 1 million fatalities) and more recently the swine flu (2009, H1N1, 14 000 fatalities).

Here, we briefly review topics heavily associated with RNA in the IAV, particularly our current knowledge on how the IAV genome is packaged, how the RNA interacts with proteins, what protein products are generated from alternative splicing in the IAV genome and our present knowledge of RNA structure. We also summarize our current understanding of the infection cycle with a focus on how the viral RNA (vRNA) sequences manage to enter the nucleus of the cell, carry out the transcription and replication process and export of the viral genetic material from the nucleus before the final budding and release of the virion. Finally, considering more than one virion can infect a given cell at the same time, it is also possible to shuffle segments between different viral strains in a process known as reassortment. We discuss briefly on the recent discoveries of RNA/RNA interactions between the vRNA segments that help in the packaging of the viral genome. Such interactions would have implications on the likelihood of such reassortment events.

Structure and the viral genome

This section is organized according to the structural features of the viral ribonucleoproteins (vRNP) rather than the numbering of the RNA segments in the genome. Hence, the organizational description of the genome will be with respect to how the virus is built from the individual components rather than the order in which the segments are numbered.

As previously mentioned, IAV consists of a negative-sense, ssRNA genome. The genome is multipartite and contains eight separate segments of various lengths; see Table 1 for details [3, 5, 9, 19–38]. More than 10 major proteins are encoded by these eight segments together with a variety of additional frame-shifted proteins and alternatively spliced products; a short description of what is known of each IAV protein is provided in the corresponding table. Other examples of segmented negative-sense ssRNA viruses are the bi-segmented arenavirus [39, 40], where the segments are called the ‘large segment’ and ‘small segment’, and the tri-segmented bunyavirus [41] (large, medium and small segments).

Each of the eight segments is individually packaged into a structurally distinct vRNP particle. A cartoon of a typical vRNP is shown in Figure 1 (top panel). The medium-gray/green beads represent the NP (Segment 5 in Table 1), where a polymer of NP molecules forms a scaffold consisting of a right-handed double-stranded (ds) helix of NPs in an ensiform (needle-like) structure [42–44]. The single-stranded vRNA sequence is represented by the light-gray/cyan lines (on the 5′ half of the sequence) that connects adjacent NP beads together and similarly the blackish/navy-blue lines (on the 3′ half). Thus, the vRNA sequence is left partially exposed as it winds around the NP-scaffolding; as though the NP forms a lining for the vRNA. The NP-scaffold does not play a significant role in protecting the vRNA because the NP-lined vRNA can be cut by ribonucleases [42]. The only part of the vRNA sequence where any dsRNA base pairing occurs is the highly conserved 12 and 13 nucleotide (nt) sequence on the 5′- and 3′-most ends of the vRNA segment, respectively [45]. Likewise, the NP-lined vRNA complex also appears to strongly restrict the formation of RNA secondary structure, though the last section of this review will discuss an emerging nuanced view.

An additional component of the vRNP is the viral polymerase. In Figure 1 (top panel, left-hand side), the viral polymerase is a hetrotrimer comprised the following proteins: polymerase acidic (PA; Segment 3), polymerase basic 1 (PB1; Segment 2) and polymerase basic 2 (PB2; Segment 1), with the polymerase located near the 5′-/3′-termini of the vRNA sequence [42]. It is important to note that [43] reports the location of the viral polymerase at the hairpin side (the exact opposite side of the 5′/3′-termini) of the ensiform structure, and the ds helix of NPs has a left-handed twist instead of a right-handed twist. Nevertheless, either way, the nucleocapsid proteins are PB1, PB2, PA and NP, where NP adopts a needle-shaped scaffold for the ssRNA and the viral polymerase is positioned at one end. The viral polymerase is essential for the transcription of the viral mRNA sequences and replication of the viral genome [44].

In the nucleocapsid, electron microscopy and tomography reveal eight rod-shaped vRNPs structures, which are organized inside of the influenza virus in a rolled-up daisy chain pattern called the ‘7 + 1’ structure. The vRNPs make contact with the matrix layer (lipid bilayer) at the budding tip of each of the virions [46]. This cluster also appears to be present in the virion during the fusion process (i.e. when the virus enters the host cell), during the budding process (i.e. when the new virion escapes from the cell) and even on entering the nucleus. Hence, some significant binding interactions are probably present that can be broken only in the nucleus of host cell.

Two additional proteins are also known to be lodged in the nucleocapsid: the nonstructural protein 1 (NS1) and 2 [NS2/nuclear export protein (NEP)], both encoded in Segment 8. NS1 plays a critical role in suppression of the 3′-end processing of the host mRNA and blocking the hosts innate immune response stimulated by interferon-α/β [47–49]. NS2, which is also called the NEP, is a multifunctional protein involved in regulating the accumulation of vRNA, antigenomic complementary RNA and viral mRNA synthesized by the vRNA polymerase [9, 50, 51]. The protein is also implicated in the efficient release of budding virions [51, 52]. It seems likely that all of the major expressed proteins are stored in the nucleocapsid to provide additional export, import and viral polymerase functionality, but the exact contents and proportions are unknown.

The virions of the IAV are observed in both a spherical or filamentous shape [4, 53, 54], where in both cases, only one copy of viral genomic appears to be present [54]; the capsid is neither a random collection of vRNPs nor multiple copies of the genome. The virions in which ‘7 + 1’ clusters of vRNP where observed were all spherical [46, 55]. The filamentous structure can sometimes encompass same space as eight spherical virions [46]. One might assume that as the virion can fit in a sphere, it would be the most likely shape of the virus. However, the filamentous structure appears to be the more commonly observed form. It is not presently known exactly what determines whether the virion is spherical or filamentous nor is it clear how the vRNPs combine in the filamentous virion; however, it seems reasonable that the filamentous-shaped virus also preserves this ‘7 + 1’ pattern.

The viral envelop is built from lipids derived from the host membrane and comprises the matrix protein 1 (M1; from Segment 7), which is located at the interface between the nucleocapsid and the lipid envelop. M1 is partly responsible for the curvature of the emerging viroid, supports the arrangement of the eight RNP on the inner surface of the lipid envelope and regulates the import and export of vRNPs into and out of the nucleus. Anchored to the viral envelope, with part of the peptide sequence occupying the border of the lipid and the nucleocapsid, extending through the lipid layer and jutting out from the surface, are two externally oriented glycoproteins: the hemmaglutinin (HA) and the NA. It is important to keep in mind that none of these envelop surface proteins are monomers. The HA forms a homotrimer (HA-trimer) that binds the sialylated glycans located on the surface of epithelial cells and initializes the process of fusion between viral envelope and the host’s endosome. The NA forms a homotetramer (NA-tetramer) and appears to play a major role in the final steps of the budding process. In addition to these glycoproteins, matrix protein 2 (M2, Segment 7) is a small protein that combines into a homotetramer (M2-tetramer) serves as a transmembrane ion channel (a proton pump). HA and NA are the main antigens of IAV and are used to identify various subtypes of IAV, e.g. H3N2 specifies an IAV with the HA belonging to Subtype 3 and the NA also belonging to Subtype 2.

Infection and the role of vRNPs

The IAV infection begins with the virus binding to the surface of an epithelial cell on recognizing the sialic acid glycans. The main sialic acid is N-acetylneuramic acid in mammals and birds [56], typically followed by galactose. These molecules can be connected with an α2,6 bond, which is abundant on the surface of cells in the upper respiratory tract of humans, or with an α2,3 bond, which is found in the digestive tract in birds or in the lungs or bronchial tissue in humans. Although sialylated glycans in epithelial cells are studied mostly for their role in disease, the function of such glycans is diverse, interfacing with signaling proteins and various structural and modulatory roles [57, 58], and in epithelial cells in general, defense against oxidative stress [59] and even nutrition [60, 61]. To defend against infections, epithelial cells secrete mucin (mucous) in the respiratory tract [27, 57, 62]. The major components of mucin are glycans terminated with sialic acid together with various antimicrobials compounds [62]. For the virions to traverse the thick mucosal layer covering the epithelial cells in the lungs, the NA-tetramers on the surface of the influenza A virion must cleave through the sialic acid within the mucous [27]. This allows the HA-trimers of the IAV to interact with the sialylated glycans on the apical cell surface [63], an event known as virus adsorption.

This event is followed by internalization of the virus and the fusion process, which takes advantage of the specific features of HA. Owing to the process of proteolytic cleavage of the HA0 precursor [64], the HA protein is divided into two subunits HA1 and HA2. HA1 forms the globular head that binds to specific glycans on the epithelial cells. HA2 anchors the HA to the viral membrane and facilitates membrane fusion of the virus with the lipid bilayer of the endosomes, where the virus enters an endosome and the fusion peptide (from HA2) interacts with the lipid bilayer. As a result of lowering the pH in the endosome, the HA subunit adopts a low pH conformation [65]. The fusion process also depends on whether the virion is spherical or filamentous. Binding for the spherical virion triggers the internalization into an endosome, where acidification causes the HA-trimers to aggregate or cluster [26, 66–68]. Filamentous virions are, however, too large to fit into a clathrin-coated pit; therefore, they must enter the cell through a process known as macropinocytosis [44, 69–71]; a pathway also used by the Ebola virus and SARS virus [8].

After membrane fusion, the virion remains ensconced within the endosome until the acidification process further dissociates the viral membrane from the vRNPs. Before the endosomes transport the contents to lysosomes for recycling, the acidification process opens the endosome resulting in dissociation of M1 from the vRNPs and releasing the vRNP into the cytosol. After dissociation of M1 and release of the viral RNPs, the vRNPs appears to remain in a single cluster of eight vRNPs (the abovementioned ‘7 + 1’ structure) and enters the nucleus through the nuclear pores still as a clustered unit. This behavior was observed, regardless of whether the virion entered the cell via endocytosis or macropinocytosis. The path to the nucleus appears to be diffusion driven [72].

In the nucleus, the complex of vRNPs separates. The heterotrimer (PA-PB1-PB2) of each vRNP is drawn to the nuclear matrix, where it can bind to RNA polymerase II (Pol II) at the C-terminal domain of the large subunit [73]. By hijacking the host transcriptional machinery, the virus can gain access to many of the host’s mRNA processing proteins (the spliceosome complex), the 5′-capped RNA primer sequences, poly(A) tail processing proteins and the host’s NEPs; this recruitment process is shown in cartoon style in Figure 1. The Pol II provides many of the required tools and raw materials needed to build a properly constructed mRNA for export and translation into proteins. One such process is known as cap-snatching [21], where a short primer sequence at the 5′ end of the host’s mRNA is captured and spliced onto the forming viral mRNA (first two panels of Figure 1). Such a primer sequence include the methylated GTP cap (m⁷GDP), a 13 nt primer sequence that permits binding of NEPs and recruits the associated cap-binding proteins (CBPs). By binding Pol II, the transport proteins are snatched before the host CBP can finish processing the host pre-mRNAs [36, 73]. The ability of the vRNA polymerase to inhibit Pol II may slow down or block important host antiviral responses [48, 74]. Other factors involved in packaging the mRNA and constructing the poly(A) tail are also required and provided by Pol II (last two panels of Figure 1). M1 and NP also appear to bind to the histones of the nucleosome [36], which may be necessary to anchor the viral polymerase during replication processes.

Each individual vRNP is capable of transcribing its own sequence using the heterotrimer (PA, PB1 and PB2) to accomplish this task. However, replication cannot be carried out by a single vRNA. Therefore, it is important to distinguish the difference between transcription and replication.

Transcription of the vRNA sequence involves obtaining a positive-sense copy of a subsequence of vRNA. Combined with the editing steps provided by the Pol II complex, a positive-sense messenger RNA (mRNA) sequence with full mRNA functionality is generated and subsequently transported to the cytosol and translated into a protein. Multiple splice products of mRNA are produced as a result of using the host spliceosome in alternative splicing. Some of these alternatively spliced products appear to be viable [22, 23, 25, 29, 75]. The general process is shown in cartoon form in Figure 1.

Replication means making a complete unedited copy of the original viral genome (with possible mutations introduced by replication ‘errors’). In the first step of replication, the vRNA sequence uses the original viral polymerase as an anchor and recruits a second free viral polymerase (in trans) to build a positive-sense copy of the vRNA sequence (vRNA → cRNA). This is a complete positive-sense strand copy of the original segment. The final result is a complete vRNP that contains only the cRNA sequence. This new cRNA must be again replicated to obtain a negative-sense copy of the sequence (cRNA → ccRNA = vRNA), where ccRNA is the newly created complement of cRNA (which will serve as a vRNA for next generations of virus). This freshly made vRNA sequence must also have access to newly synthesized copies of M1, PA, PB1, PB2 and NP; hence, these proteins (generated by the transcribed mRNA products) must be imported back into the nucleus. The PA + PB1 + PB2 must assemble into a new viral polymerase, which together with the nascent NPs and vRNA is used to form the new vRNP [20, 44, 69].

Currently, how exactly the selection between transcription and replication is determined is not completely clear. Nevertheless, based on the images of the RNPs [42, 43], there are some indications as to how the process works. One possible explanation is that the selection is influenced by the different proportions of M1 and NP [76]. In this respect, all processes are essentially transcription and are regulated by these ratios and other factors to establish whether the product will be mRNA or cRNA/ccRNA. Another view [50, 51] is that transcription involves a single viral polymerase, the so-called cis-acting polymerase, which binds somewhere in the vicinity of the promoter region (near the 3′ end of the RNA sequence) and copies the sequence moving in the 3′–5′ direction [42] and possibly responding to splicing factors, as illustrated in Figure 1 (middle). On the other hand, replication of vRNA requires transcription by an independent viral polymerase transcribing starting from the 3′-most end of the negative-sense viral ssRNA sequence—the so-called trans-acting polymerase [42] (Figure 2).

Transcription (viral mRNA production) is required to generate new copies of the viral proteins for creating a new virion. Several additional editing steps must be done before a viable mRNA can be exported and translated; the mRNA must have a properly constructed 5′-end on the nucleic acid sequence, a cap with a 5′ untranslated region (5′-UTR) and a poly-adenylated [poly(A)] tract inserted at the 3′-end of the sequence. The viral polymerase catalyzes not only RNA transcription but also cleaves the host cell’s mRNA and splices it onto the forming viral mRNA using cap-snatching by way of the PA and PB2 subunits [21]. To do this, the host Pol II and resident spliceosome factors become involved in this transcription process. For example, although M1 is translated from the unspliced mRNA, both M2 and M42 are the spliced variants. Similarly, NS1 is the unspliced version, while the NEP (NS2) is the spliced variant [44]. Additional known splice variants that have been identified in the viral genome are listed in Table 1. The Pol II is also captured to recruit its capping machinery, which provides the required cap and 13 nt primer sequence at the 5′-end of the viral transcript using PA and PB2 (a process shown as the second and third steps from the top in Figure 1). The mimicry must be complete, so that the viral mRNA is processed in the same way as the host’s mRNA. Midway along the vRNA sequence in Figure 1 (middle), one can see an image of various splice factors that introduce the possibility of alternative splicing in the transcription. Just like the host’s own mRNA, by the end of this process, the viral mRNA is probably coated with the host’s exon-junction-core proteins and associated interacting proteins, the cap-binding complex, poly(A) binding proteins and most likely other important proteins needed for export to the cytosol and translation into proteins by the host’s ribosomal RNA [77].

Some of the proteins generated by this transcription process must find their way back to the nucleus to complete the replication process. Nuclear import of the nascent copies of M1, NP, PA, PB1 and PB2 (produced from the viral mRNA) requires the special actin-activated proteins—importin β and various induced fit supplementary proteins—importins α1, α3, α5 and α7 [47, 78]; a list of proteins and their nuclear import proteins is given in Table 2 [20]. In essence, the importin β cannot by itself import the viral proteins and requires NS1 and, possibly, other alternatively spliced units of the viral mRNA transcripts to induce sufficient fit to allow the import of M1, NP and PB2. PA and PB1 are imported together with Ran-binding protein 5 (RanBP5)—a member of the importin β family. This PA + PB1 complex must then associate with PB2. NEP is small enough that it is able to pass freely through the nuclear pores. This complexity may seem like a wild ride ‘around Robin Hood’s barn’; yet, it shows that this simple virus (with <20 proteins) has evolved a sophisticated strategy to circumvent its host.

As the vRNPs enter the nucleus as a cluster, and yet the whole process of vRNA synthesis must begin with the viral polymerase transcribing in trans on a subsegment of the vRNP template (initially, the replication of cRNA), this means that there must be some free viral polymerase (i.e. PA + PB1 + PB2) available to initiate this process. One possible initial source is that some free viral polymerases are already stored inside the nucleocapsid. Later in the infection cycle, these viral polymerases are replaced with fresh copies.

In forming the viral cRNA, the NP coating is either somehow removed or used with the polymerase heterotrimer to produce the cRNA transcript, where the cRNA is assembled into a RNP (cRNA with the polymerase and scaffold of NPs) that is the template used to replicate the negative-strand sequence (ccRNA; the new copy of vRNA). Both cRNA and ccRNA bind the viral polymerase heterotrimer (PA-PB1-PB2). Hence, cRNA and ccRNA have the NP complex coating onto which the newly synthesized polymerase heterotrimer can bind. As both cRNA and its replication product ccRNA bind the viral polymerase, nuclear export requires a way to distinguish between them and export only the replicated vRNPs. The apparent factor determining this is the binding of M1 to the vRNPs, and not the cRNPs [31, 76]. The export of the newly created vRNPs is further mediated by NEP that directly interacts with NEPs CRM1 (chromosomal maintenance 1, also known as exportin 1) [19, 51, 79], where NEP may enlarge the nuclear pores [80]. Somewhere during the transit (in the cytosol) or somewhere near the point where the individual vRNPs arrive at the cellular membrane, the vRNPs manage to group into the eight segments that form the genome in the 7 + 1 pattern [46]. The selective aggregation of the vRNPs may be because of mutual interactions, proteins factors or RNA-RNA interactions [81], where the latter appears to have some support (see next section). However, the precise secondary structure and the interactions between neighboring vRNPs are unknown.

At the end of the cycle, the vRNP cluster with M1 sets into the budding site, the HA-trimers and NA-tetramers position themselves at the apical positions. Evenly distributed on the surface of the budding virion, M1 and M2 also become positioned. Based on the structure in [4, 53], the NA-tetramer appears on only one side, perhaps, to help cleave the budding virion. How this assembly mechanism might lead to organization of HA and NA into complete, budding virions and how it leads to the selective incorporation of viral genome segments will require further investigation [82]. Other proteins such as NEP and NS1 would be buried within the nucleocapsid. The envelop layer encloses the viral material enabling the virion to leave the infected cell [28, 52, 54, 71].

The role of RNA secondary structure

The vRNA is wrapped up in a scaffold composed of the NP proteins [42, 43]; therefore, only a few examples of known RNA base-pair formation have been reported for the packaged vRNPs. One observed interaction occurs at the terminal ends of each vRNA segment. In all eight vRNA gene segments, the terminal 13 and 12 nt fragments at the 5′- and 3′-most ends of the vRNA segments are conserved and known to form the vRNA promoter. This region also exhibits short hairpins of RNA secondary structure [45].

The second observed interactions were inferred from the fact that the majority of virions contain exactly one complete genome packed in the daisy chain structures [9, 46, 81, 83]. This daisy chain connectivity has long suggested that RNA/RNA interactions occur between the vRNPs that help bind them together. Recent observations have shown that there appear to be some distinct RNA/RNA interactions between different vRNPs [83, 84]. In particular, one study [84] revealed that a vRNA hairpin from Segment 2 (the PB1 Segment) and a vRNA hairpin from Segment 8 (the NP Segment) appear to mutually combine via the complementary base-pair sequences of the loops. (Two RNA hairpins bound together by complementary sequences found in their loops are called a ‘kissing loop’ [85, 86].) The evidence for the kissing loop is supported by sequence alignment methods and in vitro experiments involving disrupting the complementary sequences and measuring the relevant fractions of products via native agarose gel electrophoresis [84]. Hence, the linking of the segments is most likely because of RNA/RNA interactions between adjacent segments. However, these interaction networks of kissing loops also appear to be strain-dependent [84]. This suggests that for a particular subtype of IAV to jump between species, not only does the successful subtype has to acquire specific mutations in HA, access effective nuclear import/export factors and produce proteins that successfully block host innate immunity but, additionally, the subtype requires a proper connectivity between vRNPs. This might partly explain why reassortment of various subtypes is not a common occurrence and how the strains manage to remain unique.

Kobayashi et al. [87] attempted to identify the actual RNA secondary structure (hairpins) in Segment 7 (the M Segment) and test their stability at different temperatures. (RNA secondary structure differs from protein secondary structure in that it shows the base pairing between diverse parts of the RNA sequence or between different RNA sequences; therefore, more topology information can be inferred from RNA secondary structure data.) In this analysis, the functionality of conserved stem-loop structures was directly tested using reverse genetic experiments to introduce synonymous mutations designed to disrupt secondary structures predicted at three locations. These disruptions tended to attenuate the infectivity of the recombinant virus. The disruption of the stem-loop structure at nucleotide positions 219–240, 950–964 and 967-994 of Segment 7 (M Segment) showed a significant temperature dependence, which is characteristic of decrease in RNA stability in these regions. This suggests that the conserved secondary structures predicted in Segment 7 are involved in the production of infectious viral particles during IAV replication.

The abovementioned studies strongly suggest that base pairing in the vRNP may play a significant role in the mutual binding of the vRNP structures to one another. The data from Kobayashi et al. is stable. In Figure 3, an actual scan of the free energy and variation in optimal and suboptimal structures of Segment 7 (Segment M) is shown for four different viral strains in the following order: A/Baltimore/0026/2016(H1N1), A/England/195/2009(H1N1), A/udorn/1972(H3N2) and A/ruddy turnstone/DE/509531/2007(H5N1). Sequences of this segment were obtained from the influenza database at https://www.fludb.org [88]. The spectra represent scanning for RNA secondary structure over a window size of 200 nt and, based on the predicted suboptimal structures obtained in the scan, detecting regions where a particular secondary structure motif of RNA repeatedly appears in many (or most) of the optimal and suboptimal structures. The very dark gray/blue boxes represent regions of stable RNA secondary structure that are largely independent of the scanning window size, and segments of highly variable secondary structure are shown by their lack of any color. (The very short gray/magenta boxes indicate the center of the hairpin in the stable (gray/blue) regions of the spectra.) Covering the regions at loci 219–240, 950–964 and 967–994 of Segment 7 reported by Kobayashi et al. [87], Figure 3 shows three long semi-opaque gray/green boxes (also marked by the dark arrows, bottom) overlapping these respective fragments. The scan was performed using vswindow; an expanded application of vs_subopt [89] (http://www.rna.it-chiba.ac.jp/∼vsfold/vs_subopt/) with default settings (window size 200 nt, Kuhn length 4 nt), and the spectrum is generated using Genepoem [90] (http://www.rna.it-chiba.ac.jp/∼vswindow/cgi-bin/index.cgi). Modifying the settings—such as changing the scanning window length or the Kuhn length (a measure of the stiffness of the RNA [89])—generates similar results. As the RNP (RNA wrapped in NPs) is largely exposed [42], the stable RNA secondary structure is also likely to appear in the RNP complex.

However, it is important to note that the calculations by Kobayashi et al. were for the cRNA sequence. Therefore, in Figure 4, the results for the actual vRNA sequence are shown. Regions at positions 950–964 and 967–994 appear to be weaker (represented in Figure by less intense gray/dark-blue), indicating that there is more variation in possible suboptimal structures for these regions of the vRNA. However, the region at position 219–240 appears prominently in both Figures 3 and 4. Perhaps, these variations in RNA stability help to distinguish the cRNA from vRNA in the steps toward nuclear export.

The somewhat regular pattern of gray/dark-blue bands in strong secondary structure appears on all of the segments of cRNA and vRNA. However, Figures 3 and 4 also suggest that these gray/blue bands are somewhat strain-dependent—as also stressed in reports such as Gavazzi et al. [84] and Kobayashi et al. [87].

The results from Kobayashi et al. only look at Segment 7; they do not specify where these hairpins might bind to other segments of the vRNP. However, the kissing loop reported in Gavazzi et al. [84] can also be examined in this way. Figure 5 shows a scan of Segments 2 and 8 for the same corresponding vRNA sequence reported in [84], strain A/finch/England/2051/1991(H5N2). The same default settings are used as above for vswindow/Genepoem. Likewise, in Figure 5, the two connected semi-opaque gray/green boxes is used to indicate the corresponding sequence ranges reported in Gavazzi et al. [84], and the dashed lines indicate the connection between Segment 2 and Segment 8. Though both Segment 2 and Segment 8 are clearly visible in the spectrum—indicating that the structures are relatively stable—the intensity of the gray/blue in these regions is less than some other cases in the same spectrum. This lower intensity gray/blue indicates that there is a greater degree of variation in the RNA secondary structure in these locations. In line with the observation that the linking together of vRNPs is the result of kissing hairpins, this may suggest that there are many such contacts regularly dispersed along the vRNA segments that act like a bar code that defines how the pieces should be fit together with one another and identify the vRNPs specific to a particular strain. Though this may diminish some of the long-held fears of facile reassortment of the common flu with such strains as H5N1, it also could explain how individual strains evolve and compete with one another too.

Funding

WKD, ML, DP were supported by grants from the Polish National Science Centre (grant numbers 2014/15/B/ST6/05082 and 2013/09/B/NZ2/00121); ML, DP were funded by Foundation for Polish Science (TEAM to D.P.). DP was supported by grant 1U54DK107967-01 ‘Nucleome Positioning System for Spatiotemporal Genome Organization and Regulation’ within 4DNucleome NIH program.

Wayne Dawson is a contributing researcher at the Bio-information Lab in Yayoi campus at the University of Tokyo and also associated with Center of New Technologies (CeNT) in Warsaw, Poland. His research is primarily directed to RNA structure prediction and structural biology and recently focused on chromatin structure.

Michał Łaźniewski is a researcher at the Laboratory of Functional and Structural Genomics at the University of Warsaw in the Center of New Technologies (CeNT) and Department of Physical Chemistry in the Faculty of Pharmacy at the Medical University of Warsaw. His research is in bioinformatics and structural biology with a focus on medical applications using bioinformatics.

Dariusz Plewczyński is a professor and head of the Laboratory of Functional and Structural Genomics at the University of Warsaw in the Center of New Technologies (CeNT) in Warsaw, Poland and also the affiliated with the Faculty of Mathematics and Information Science at Warsaw University of Technology in Warsaw. His main expertise covers computational genomics, biostatistics and bioinformatics. He is actively applying computational modeling to various interdisciplinary problems in Human Genomics.

References