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.

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.

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.

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.

Analysis of wheat SAGE tags reveals evidence for widespread antisense transcription

BACKGROUND: Serial Analysis of Gene Expression (SAGE) is a powerful tool for genome-wide transcription studies. Unlike microarrays, it has the ability to detect novel forms of RNA such as alternatively spliced and antisense transcripts, without the need for prior knowledge of their existence. One limitation of using SAGE on an organism with a complex genome and lacking detailed sequence information, such as the hexaploid bread wheat Triticum aestivum, is accurate annotation of the tags generated. Without accurate annotation it is impossible to fully understand the dynamic processes involved in such complex polyploid organisms. Hence we have developed and utilised novel procedures to characterise, in detail, SAGE tags generated from the whole grain transcriptome of hexaploid wheat. RESULTS: Examination of 71,930 Long SAGE tags generated from six libraries derived from two wheat genotypes grown under two different conditions suggested that SAGE is a reliable and reproducible technique for use in studying the hexaploid wheat transcriptome. However, our results also showed that in poorly annotated and/or poorly sequenced genomes, such as hexaploid wheat, considerably more information can be extracted from SAGE data by carrying out a systematic analysis of both perfect and "fuzzy" (partially matched) tags. This detailed analysis of the SAGE data shows first that while there is evidence of alternative polyadenylation this appears to occur exclusively within the 3' untranslated regions. Secondly, we found no strong evidence for widespread alternative splicing in the developing wheat grain transcriptome. However, analysis of our SAGE data shows that antisense transcripts are probably widespread within the transcriptome and appear to be derived from numerous locations within the genome. Examination of antisense transcripts showing sequence similarity to the Puroindoline a and Puroindoline b genes suggests that such antisense transcripts might have a role in the regulation of gene expression. CONCLUSION: Our results indicate that the detailed analysis of transcriptome data, such as SAGE tags, is essential to understand fully the factors that regulate gene expression and that such analysis of the wheat grain transcriptome reveals that antisense transcripts maybe widespread and hence probably play a significant role in the regulation of gene expression during grain development.

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.

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.

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.

Regulation of the cardiomyocyte transcriptome vs translatome by endothelin-1 and insulin: translational regulation of 5' terminal oligopyrimidine tract (TOP) mRNAs by insulin

Background: Changes in cellular phenotype result from underlying changes in mRNA transcription and translation. Endothelin-1 stimulates cardiomyocyte hypertrophy with associated changes in mRNA/protein expression and an increase in the rate of protein synthesis. Insulin also increases the rate of translation but does not promote overt cardiomyocyte hypertrophy. One mechanism of translational regulation is through 5' terminal oligopyrimidine tracts (TOPs) that, in response to growth stimuli, promote mRNA recruitment to polysomes for increased translation. TOP mRNAs include those encoding ribosomal proteins, but the full panoply remains to be established. Here, we used microarrays to compare the effects of endothelin-1 and insulin on the global transcriptome of neonatal rat cardiomyocytes, and on mRNA recruitment to polysomes (i.e. the translatome). Results: Globally, endothelin-1 and insulin (1 h) promoted >1.5-fold significant (false discovery rate < 0.05) changes in expression of 341 and 38 RNAs, respectively. For these transcripts with this level of change there was little evidence of translational regulation. However, 1336 and 712 RNAs had >1.25-fold significant changes in expression in total and/or polysomal RNA induced by endothelin-1 or insulin, respectively, of which ~35% of endothelin-1-responsive and ~56% of insulin-responsive transcripts were translationally regulated. Of mRNAs for established proteins recruited to polysomes in response to insulin, 49 were known TOP mRNAs with a further 15 probable/possible TOP mRNAs, but 49 had no identifiable TOP sequences or other consistent features in the 5' untranslated region. Conclusions: Endothelin-1, rather than insulin, substantially affects global transcript expression to promote cardiomyocyte hypertrophy. Effects on RNA recruitment to polysomes are subtle, with differential effects of endothelin-1 and insulin on specific transcripts. Furthermore, although insulin promotes recruitment of TOP mRNAs to cardiomyocyte polysomes, not all recruited mRNAs are TOP mRNAs.

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.