Direct repeats in the flavivirus 3' untranslated region; a strategy for survival in the environment?

Previously, direct repeats (DRs) of 20-70 nucleotides were identified in the 3' untranslated regions (3'UTR) of flavivirus sequences. To address their functional significance, we have manually generated a pan-flavivirus 3'UTR alignment and correlated it with the corresponding predicted RNA secondary structures. This approach revealed that intra-group-conserved DRs evolved from six long repeated sequences (LRSs) which, as approximately 200-nucleotide domains were preserved only in the genomes of the slowly evolving tick-borne flaviviruses. We propose that short DRs represent the evolutionary remnants of LRSs rather than distinct molecular duplications. The relevance of DRs to virus replication enhancer function, and thus survival, is discussed.

Molecular archaeology of Flaviviridae untranslated regions: duplicated RNA structures in the replication enhancer of flaviviruses and pestiviruses emerged via convergent evolution.

RNA secondary structures in the 3'untranslated regions (3'UTR) of the viruses of the family Flaviviridae, previously identified as essential (promoters) or beneficial (enhancers) for replication, have been analysed. Duplicated enhancer elements are revealed as a global feature in the evolution of the 3'UTR of distantly related viruses within the genera Flavivirus and Pestivirus. For the flaviviruses, duplicated structures occur in the 3'UTR of all four distantly related ecological virus subgroups (tick-borne, mosquito-borne, no known vector and insect-specific flaviviruses (ISFV). RNA structural differences distinguish tick-borne flaviviruses with discrete pathogenetic characteristics. For Aedes- and Culex-associated ISFV, secondary RNA structures with different conformations display numerous short ssRNA direct repeats, exposed as loops and bulges. Long quadruplicate regions comprise almost the entire 3'UTR of Culex-associated ISFV. Extended duplicated sequence and associated RNA structures were also discovered in the 3'UTR of pestiviruses. In both the Flavivirus and Pestivirus genera, duplicated RNA structures were localized to the enhancer regions of the 3'UTR suggesting an adaptive role predominantly in wild-type viruses. We propose sequence reiteration might act as a scaffold for dimerization of proteins involved in assembly of viral replicase complexes. Numerous nucleotide repeats exposed as loops/bulges might also interfere with host immune responses acting as a molecular sponge to sequester key host proteins or microRNAs.

Plant U13 orthologues and orphan snoRNAs identified by RNomics of RNA from Arabidopsis nucleoli

Small nucleolar RNAs (snoRNAs) and small Cajal body-specific RNAs (scaRNAs) are non-coding RNAs whose main function in eukaryotes is to guide the modification of nucleotides in ribosomal and spliceosomal small nuclear RNAs, respectively. Full-length sequences of Arabidopsis snoRNAs and scaRNAs have been obtained from cDNA libraries of capped and uncapped small RNAs using RNA from isolated nucleoli from Arabidopsis cell cultures. We have identified 31 novel snoRNA genes (9 box C/D and 22 box H/ACA) and 15 new variants of previously described snoRNAs. Three related capped snoRNAs with a distinct gene organization and structure were identified as orthologues of animal U13snoRNAs. In addition, eight of the novel genes had no complementarity to rRNAs or snRNAs and are therefore putative orphan snoRNAs potentially reflecting wider functions for these RNAs. The nucleolar localization of a number of the snoRNAs and the localization to nuclear bodies of two putative scaRNAs was confirmed by in situ hybridization. The majority of the novel snoRNA genes were found in new gene clusters or as part of previously described clusters. These results expand the repertoire of Arabidopsis snoRNAs to 188 snoRNA genes with 294 gene variants.

Origin and evolution of 3'UTR of flaviviruses: long direct repeats as a basis for the formation of secondary structures and their significance for virus transmission

The 3' untranslated regions (3'UTRs) of flaviviruses are reviewed and analyzed in relation to short sequences conserved as direct repeats (DRs). Previously, alignments of the 3'UTRs have been constructed for three of the four recognized flavivirus groups, namely mosquito-borne, tick-borne, and nonclassified flaviviruses (MBFV, TBFV, and NCFV, respectively). This revealed (1) six long repeat sequences (LRSs) in the 3'UTR and open-reading frame (ORF) of the TBFV, (2) duplication of the 3'UTR of the NCFV by intramolecular recombination, and (3) the possibility of a common origin for all DRs within the MBFV. We have now extended this analysis and review it in the context of all previous published analyses. This has been achieved by constructing a robust alignment between all flaviviruses using the published DRs and secondary RNA structures as "anchors" to reveal additional homologies along the 3'UTR. This approach identified nucleotide regions within the MBFV, NKV (no-known vector viruses), and NCFV 3'UTRs that are homologous to different LRSs in the TBFV 3'UTR and ORF. The analysis revealed that some of the DRs and secondary RNA structures described individually within each flavivirus group share common evolutionary origins. The 3'UTR of flaviviruses, and possibly the ORF, therefore probably evolved through multiple duplication of an RNA domain, homologous to the LRS previously identified only in the TBFV. The short DRs in all virus groups appear to represent the evolutionary remnants of these domains rather than resulting from new duplications. The relevance of these flavivirus DRs to evolution, diversity, 3'UTR enhancer function, and virus transmission is reviewed.

Origin and evolution of flavivirus 5'UTRs and panhandles: trans-terminal duplications?

Flavivirus replication is mediated by interactions between complementary ssRNA sequences of the 5'- and 3'-termini that form dsRNA cyclisation stems or panhandles, varying in length, sequence and specific location in the mosquito-borne, tick-borne, non-vectored and non-classified flaviviruses. In this manuscript we manually aligned the flavivirus 5'UTRs and adjacent capsid genes and revealed significantly more homology than has hitherto been identified. Analysis of the alignments revealed that the panhandles represent evolutionary remnants of a long cyclisation domain that probably emerged through duplication of one of the UTR termini.

A simple and rapid method for preserving RNA of aquatic invertebrates for ecotoxicogenomics

Here we describe a novel, inexpensive and simple method for preserving RNA that reduces handling stress in aquatic invertebrates following ecotoxicogenomic experimentation. The application of the method is based on transcriptomic experiments conducted on Daphnia magna, but may easily be applied on a range of other aquatic organisms of a particular size with e.g. amphipod Gammarus pulex representing an upper size limit. We explain in detail how to apply this new method, named the "Cylindrical Sieve (CS) system", and highlight its advantages and disadvantages.

Mass spectroscopic characterization of the coronavirus infectious bronchitis virus nucleoprotein and elucidation of the role of phosphorylation in RNA binding by using surface plasmon resonance

Phosphorylation of the coronavirus nucleoprotein (N protein) has been predicted to play a role in RNA binding. To investigate this hypothesis, we examined the kinetics of RNA binding between nonphosphorylated and phosphorylated infectious bronchitis virus N protein with nonviral and viral RNA by surface plasmon resonance (Biacore). Mass spectroscopic analysis of N protein identified phosphorylation sites that were proximal to RNA binding domains. Kinetic analysis, by surface plasmon resonance, indicated that nonphospborylated N protein bound with the same affinity to viral RNA as phosphorylated N protein. However, phosphorylated N protein bound to viral RNA with a higher binding affinity than nonviral RNA, suggesting that phosphorylation of N protein determined the recognition of virus RNA. The data also indicated that a known N protein binding site (involved in transcriptional regulation) consisting of a conserved core sequence present near the 5' end of the genome (in the leader sequence) functioned by promoting high association rates of N protein binding. Further analysis of the leader sequence indicated that the core element was not the only binding site for N protein and that other regions functioned to promote high-affinity binding.

Factors required for the uridylylation of the foot-and-mouth disease virus 3B1, 3B2, and 3B3 peptides by the RNA-dependent RNA polymerase (3D(pol)) in vitro

The 5' terminus of picornavirus genomic RNA is covalently linked to the virus-encoded peptide 313 (VTg). Foot-and-mouth disease virus (FMDV) is unique in encoding and using 3 distinct forms of this peptide. These peptides each act as primers for RNA synthesis by the virus-encoded RNA polymerase 3D(pol). To act as the primer for positive-strand RNA synthesis, the 3B peptides have to be uridylylated to form VPgpU(pU). For certain picornaviruses, it has been shown that this reaction is achieved by the 3D(pol) in the presence of the 3CD precursor plus an internal RNA sequence termed a cis-acting replication element (cre). The FMDV ere has been identified previously to be within the 5' untranslated region, whereas all other picornavirus cre structures are within the viral coding region. The requirements for the in vitro uridylylation of each of the FMDV 313 peptides has now been determined, and the role of the FMDV ere (also known as the 3B-uridylylation site, or bus) in this reaction has been analyzed. The poly(A) tail does not act as a significant template for FMDV 3B uridylylation.

An RNA molecule that specifically inhibits G-protein-coupled receptor kinase 2 in vitro

G-protein-coupled receptors are desensitized by a two-step process. In a first step, G-protein-coupled receptor kinases (GRKs) phosphorylate agonist-activated receptors that subsequently bind to a second class of proteins, the arrestins. GRKs can be classified into three subfamilies, which have been implicated in various diseases. The physiological role(s) of GRKs have been difficult to study as selective inhibitors are not available. We have used SELEX (systematic evolution of ligands by exponential enrichment) to develop RNA aptamers that potently and selectively inhibit GRK2. This process has yielded an aptamer, C13, which bound to GRK2 with a high affinity and inhibited GRK2-catalyzed rhodopsin phosphorylation with an IC50 of 4.1 nM. Phosphorylation of rhodopsin catalyzed by GRK5 was also inhibited, albeit with 20-fold lower potency (IC50 of 79 nM). Furthermore, C13 reveals significant specificity, since almost no inhibitory activity was detectable testing it against a panel of 14 other kinases. The aptamer is two orders of magnitude more potent than the best GRK2 inhibitors described previously and shows high selectivity for the GRK family of protein kinases.

Induction of aquaporin 1 but not aquaporin 4 messenger RNA in rat primary brain microvessel endothelial cells in culture

Aquaporins (AQPs) are a family of proteins that mediate water transport across cells, but the extent to which they are involved in water transport across endothelial cells of the blood-brain barrier is not clear. Expression of AQP1 and AQP4 in rat brain microvessel endothelial cells was investigated in order to determine whether these isoforms were present and, in particular, to examine the hypothesis that brain endothelial expression of AQPs is dynamic and regulated by astrocytic influences. Reverse-transcriptase-polymerase chain reaction (RT-PCR) and immunocytochemistry showed that AQP1 mRNA and protein are present at very low levels in primary rat brain microvessel endothelial cells, and are up-regulated in passaged cells. Upon passage, endothelial cell expression of mdr1a mRNA is decreased, indicating loss of blood-brain barrier phenotype. In passage 4 endothelial cells, AQP1 mRNA levels are reduced by coculture above rat astrocytes, demonstrating that astrocytic influences are important in maintaining the low levels of AQP1 characteristic of the blood-brain barrier endothelium. Reverse-transcriptase-PCR revealed very low levels of AQP1 mRNA present in the RBE4 rat brain microvessel endothelial cell line, with no expression detected in primary cultures of rat astrocytes or in the C6 rat glioma cell line. In contrast, AQP4 mRNA is strongly expressed in astrocytes, but no expression is found in primary or passaged brain microvessel endothelial cells, or in RBE4 or C6 cells. Our results support the concept that expression of AQP1, which is seen in many non-brain endothelia, is suppressed in the specialized endothelium of the blood-brain barrier.

Ribose 2′-O-methylation provides a molecular signature for the distinction of self and non-self mRNA dependent on the RNA sensor Mda5

The 5'-cap-structures of higher eukaryote mRNAs are ribose 2'-O-methylated. Likewise, a number of viruses replicating in the cytoplasm of eukayotes have evolved 2'-O-methyltransferases to modify autonomously their mRNAs. However, a defined biological role of mRNA 2'-O-methylation remains elusive. Here we show that viral mRNA 2'-O-methylation is critically involved in subversion of type-I-interferon (IFN-I) induction. We demonstrate that human and murine coronavirus 2'-O-methyltransferase mutants induce increased IFN-I expression, and are highly IFN-I sensitive. Importantly, IFN-I induction by 2'-O-methyltransferase-deficient viruses is dependent on the cytoplasmic RNA sensor melanoma differentiation-associated gene 5 (MDA5). This link between MDA5-mediated sensing of viral RNA and mRNA 2'-O-methylation suggests that RNA modifications, such as 2'-O-methylation, provide a molecular signature for the discrimination of self and non-self mRNA.

SARS coronavirus unique domain: three-domain molecular architecture in solution and RNA binding

Nonstructural protein 3 of the severe acute respiratory syndrome (SARS) coronavirus includes a "SARS-unique domain" (SUD) consisting of three globular domains separated by short linker peptide segments. This work reports NMR structure determinations of the C-terminal domain (SUD-C) and a two-domain construct (SUD-MC) containing the middle domain (SUD-M) and the C-terminal domain, and NMR data on the conformational states of the N-terminal domain (SUD-N) and the SUD-NM two-domain construct. Both SUD-N and SUD-NM are monomeric and globular in solution; in SUD-NM, there is high mobility in the two-residue interdomain linking sequence, with no preferred relative orientation of the two domains. SUD-C adopts a frataxin like fold and has structural similarity to DNA-binding domains of DNA-modifying enzymes. The structures of both SUD-M (previously determined) and SUD-C (from the present study) are maintained in SUD-MC, where the two domains are flexibly linked. Gel-shift experiments showed that both SUD-C and SUD-MC bind to single-stranded RNA and recognize purine bases more strongly than pyrimidine bases, whereby SUD-MC binds to a more restricted set of purine-containing RNA sequences than SUD-M. NMR chemical shift perturbation experiments with observations of (15)N-labeled proteins further resulted in delineation of RNA binding sites (i.e., in SUD-M, a positively charged surface area with a pronounced cavity, and in SUD-C, several residues of an anti-parallel beta-sheet). Overall, the present data provide evidence for molecular mechanisms involving the concerted actions of SUD-M and SUD-C, which result in specific RNA binding that might be unique to the SUD and, thus, to the SARS coronavirus.

Roles of the translationally controlled tumour protein (TCTP) and the double-stranded RNA-dependent protein kinase, PKR, in cellular stress responses

Translationally controlled tumour protein (TCTP) is a highly conserved protein present in all eukaryotic organisms. Various cellular functions and molecular interactions have been ascribed to this protein, many related to its growth-promoting and antiapoptotic properties. TCTP levels are highly regulated in response to various cellular stimuli and stresses. We have shown recently that the double-stranded RNA-dependent protein kinase, PKR, is involved in translational regulation of TCTP. Here we extend these studies by demonstrating that TCTP is downregulated in response to various proapoptotic treatments, in particular agents that induce Ca++ stress, in a PKR-dependent manner. This regulation requires phosphorylation of protein synthesis factor eIF2α. Since TCTP has been characterized as an antiapoptotic and Ca++-binding protein, we asked whether it is involved in protecting cells from Ca++-stress-induced apoptosis. Overexpression of TCTP partially protects cells against thapsigargin-induced apoptosis, as measured using caspase-3 activation assays, a nuclear fragmentation assay, using fluorescence-activated cell sorting analysis, and time-lapse video microscopy. TCTP also protects cells against the proapoptotic effects of tunicamycin and etoposide, but not against those of arsenite. Our results imply that cellular TCTP levels influence sensitivity to apoptosis and that PKR may exert its proapoptotic effects at least in part through downregulation of TCTP via eIF2α phosphorylation.

Replication enhancer elements within the open reading frame of tick-borne encephalitis virus and their evolution within the Flavivirus genus

We provide experimental evidence of a replication enhancer element (REE) within the capsid gene of tick-borne encephalitis virus (TBEV, genus Flavivirus). Thermodynamic and phylogenetic analyses predicted that the REE folds as a long stable stem–loop (designated SL6), conserved among all tick-borne flaviviruses (TBFV). Homologous sequences and potential base pairing were found in the corresponding regions of mosquito-borne flaviviruses, but not in more genetically distant flaviviruses. To investigate the role of SL6, nucleotide substitutions were introduced which changed a conserved hexanucleotide motif, the conformation of the terminal loop and the base-paired dsRNA stacking. Substitutions were made within a TBEV reverse genetic system and recovered mutants were compared for plaque morphology, single-step replication kinetics and cytopathic effect. The greatest phenotypic changes were observed in mutants with a destabilized stem. Point mutations in the conserved hexanucleotide motif of the terminal loop caused moderate virus attenuation. However, all mutants eventually reached the titre of wild-type virus late post-infection. Thus, although not essential for growth in tissue culture, the SL6 REE acts to up-regulate virus replication. We hypothesize that this modulatory role may be important for TBEV survival in nature, where the virus circulates by non-viraemic transmission between infected and non-infected ticks, during co-feeding on local rodents.