A score system for quality evaluation of RNA sequence tags: an improvement for gene expression profiling

Background: High-throughput molecular approaches for gene expression profiling, such as Serial Analysis of Gene Expression (SAGE), Massively Parallel Signature Sequencing (MPSS) or Sequencing-by-Synthesis (SBS) represent powerful techniques that provide global transcription profiles of different cell types through sequencing of short fragments of transcripts, denominated sequence tags. These techniques have improved our understanding about the relationships between these expression profiles and cellular phenotypes. Despite this, more reliable datasets are still necessary. In this work, we present a web-based tool named S3T: Score System for Sequence Tags, to index sequenced tags in accordance with their reliability. This is made through a series of evaluations based on a defined rule set. S3T allows the identification/selection of tags, considered more reliable for further gene expression analysis. Results: This methodology was applied to a public SAGE dataset. In order to compare data before and after filtering, a hierarchical clustering analysis was performed in samples from the same type of tissue, in distinct biological conditions, using these two datasets. Our results provide evidences suggesting that it is possible to find more congruous clusters after using S3T scoring system. Conclusion: These results substantiate the proposed application to generate more reliable data. This is a significant contribution for determination of global gene expression profiles. The library analysis with S3T is freely available at http://gdm.fmrp.usp.br/s3t/.S3T source code and datasets can also be downloaded from the aforementioned website.

Translation of the flavivirus Kunjin NS3 gene in cis but not its RNA sequence or secondary structure is essential for efficient RNA packaging

Our previous studies using trans-complementation analysis of Kunjin virus (KUN) full-length cDNA clones harboring in-frame deletions in the NS3 gene demonstrated the inability of these defective complemented RNAs to be packaged into virus particles (W. J. Liu, P. L. Sedlak, N. Kondratieva, and A. A. Khromykh, J. Virol. 76:10766-10775). In this study we aimed to establish whether this requirement for NS3 in RNA packaging is determined by the secondary RNA structure of the NS3 gene or by the essential role of the translated NS3 gene product. Multiple silent mutations of three computer-predicted stable RNA structures in the NS3 coding region of KUN replicon RNA aimed at disrupting RNA secondary structure without affecting amino acid sequence did not affect RNA replication and packaging into virus-like particles in the packaging cell line, thus demonstrating that the predicted conserved RNA structures in the NS3 gene do not play a role in RNA replication and/or packaging. In contrast, double frameshift mutations in the NS3 coding region of full-length KUN RNA, producing scrambled NS3 protein but retaining secondary RNA structure, resulted in the loss of ability of these defective RNAs to be packaged into virus particles in complementation experiments in KUN replicon-expressing cells. Furthermore, the more robust complementation-packaging system based on established stable cell lines producing large amounts of complemented replicating NS3-deficient replicon RNAs and infection with KUN virus to provide structural proteins also failed to detect any secreted virus-like particles containing packaged NS3-deficient replicon RNAs. These results have now firmly established the requirement of KUN NS3 protein translated in cis for genome packaging into virus particles.

Characterization of a Full-Length TTP Family Member Association with RNA Sequence Elements

Post-transcriptional regulation of cytoplasmic mRNAs is an efficient mechanism of regulating the amounts of active protein within a eukaryotic cell. RNA sequence elements located in the untranslated regions of mRNAs can influence transcript degradation or translation through associations with RNA-binding proteins. Tristetraprolin (TTP) is the best known member of a family of CCCH zinc finger proteins that targets adenosine-uridine rich element (ARE) binding sites in the 3’ untranslated regions (UTRs) of mRNAs, promoting transcript deadenylation through the recruitment of deadenylases. More specifically, TTP has been shown to bind AREs located in the 3’-UTRs of transcripts with known roles in the inflammatory response. The mRNA-binding region of the protein is the highly conserved CCCH tandem zinc finger (TZF) domain. The synthetic TTP TZF domain has been shown to bind with high affinity to the 13-mer sequence of UUUUAUUUAUUUU. However, the binding affinities of full-length TTP family members to the same sequence and its variants are unknown. Furthermore, the distance needed between two overlapping or neighboring UUAUUUAUU 9-mers for tandem binding events of a full-length TTP family member to a target transcript has not been explored. To address these questions, we recombinantly expressed and purified the full-length C. albicans TTP family member Zfs1. Using full-length Zfs1, tagged at the N-terminus with maltose binding protein (MBP), we determined the binding affinities of the protein to the optimal TTP binding sequence, UUAUUUAUU. Fluorescence anisotropy experiments determined that the binding affinities of MBP-Zfs1 to non-canonical AREs were influenced by ionic buffer strength, suggesting that transcript selectivity may be affected by intracellular conditions. Furthermore, electrophoretic mobility shift assays (EMSAs) revealed that separation of two core AUUUA sequences by two uridines is sufficient for tandem binding of MBP-Zfs1. Finally, we found evidence for tandem binding of MBP-Zfs1 to a 27-base RNA oligonucleotide containing only a single ARE-binding site, and showed that this was concentration and RNA length dependent; this phenomenon had not been seen previously. These data suggest that the association of the TTP TZF domain and the TZF domains of other species, to ARE-binding sites is highly conserved. Domains outside of the TZF domain may mediate transcript selectivity in changing cellular conditions, and promote protein-RNA interactions not associated with the ARE-binding TZF domain.

In summary, the evidence presented here suggests that Zfs1-mediated decay of mRNA targets may require additional interactions, in addition to ARE-TZF domain associations, to promote transcript destabilization and degradation. These studies further our understanding of post-transcriptional steps in gene regulation.

RNA pentaloop structures as effective targets of regulators belonging to the RsmA/CsrA protein family.

In the Gac/Rsm signal transduction pathway of Pseudomonas fluorescens CHA0, the dimeric RNA-binding proteins RsmA and RsmE, which belong to the vast bacterial RsmA/CsrA family, effectively repress translation of target mRNAs containing a typical recognition sequence near the translation start site. Three small RNAs (RsmX, RsmY, RsmZ) with clustered recognition sequences can sequester RsmA and RsmE and thereby relieve translational repression. According to a previously established structural model, the RsmE protein makes optimal contacts with an RNA sequence 5'- (A)/(U)CANGGANG(U)/(A)-3', in which the central ribonucleotides form a hexaloop. Here, we questioned the relevance of the hexaloop structure in target RNAs. We found that two predicted pentaloop structures, AGGGA (in pltA mRNA encoding a pyoluteorin biosynthetic enzyme) and AAGGA (in mutated pltA mRNA), allowed effective interaction with the RsmE protein in vivo. By contrast, ACGGA and AUGGA were poor targets. Isothermal titration calorimetry measurements confirmed the strong binding of RsmE to the AGGGA pentaloop structure in an RNA oligomer. Modeling studies highlighted the crucial role of the second ribonucleotide in the loop structure. In conclusion, a refined structural model of RsmE-RNA interaction accommodates certain pentaloop RNAs among the preferred hexaloop RNAs.

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.

Cytogenetic mapping of H1 histone and ribosomal RNA genes in hybrids between catfish species Pseudoplatystoma corruscans and Pseudoplatystoma reticulatum

A physical chromosome mapping of the H1 histone and 5S and 18S ribosomal RNA (rRNA) genes was performed in interspecific hybrids of Pseudoplatystoma corruscans and P. reticulatum. The results showed that 5S rRNA clusters were located in the terminal region of 2 chromosomes. H1 histone and 18S ribosomal genes were co-localized in the terminal portion of 2 chromosomes (distinct from the chromosomes bearing 5S clusters). These results represent the first report of association between H1 histone and 18S genes in fish genomes. The chromosome clustering of ribosomal and histone genes was already reported for different organisms and suggests a possible selective pressure for the maintenance of this association. © 2012 S. Karger AG, Basel.

Induktion und Analyse spezifischer Immunantworten nach intranodaler Immunisierung mit RNA im murinen Modell

Die Wirksamkeit einer Vakzine ist von vielen Parametern abhängig. Dazu gehören unter anderen: das ausgewählte Antigen, die Formulation in der das Antigen benutzt wird sowie die Applikationsroute. Antigen-kodierende Ribonukleinsäuren (RNA) gilt heutzutage als eine sichere und effiziente Alternative zu traditionellen Impfstoff-Formulierungen, wie Peptiden, rekombinanten Proteinen, viralen Systemen oder DNA basierten Impfstoffen. Bezüglich des Applikationsortes repräsentiert der Lymphknoten ein optimales Milieu für die Interaktion zwischen antigenpräsentierenden Zellen und T-Zellen. Vor diesem Hintergrund war die Zielsetzung dieser Arbeit, ein auf direktem in vivo Transfer von Antigen-kodierender in vitro transkribierter RNA (IVT-RNA) basierendes Impfverfahren zu entwickeln, zu charakterisieren und auf seine anti-tumorale Wirksamkeit zu testen. In der vorliegenden Arbeit konnte gezeigt werden, dass dendritische Zellen (DCs) in vitro hocheffizient mit IVT-RNA transfiziert werden können und eine hohe stimulatorische Kapazität besitzen. Durch Sequenzmodifikation der IVT-RNA konnten wir die Transkriptstabilität und Translationseffizienz erhöhen was zu einer Steigerung der stimulatorischen Kapazität in vivo führte. Darüber hinaus untersuchten wir die Auswirkung der Insertion eines Signalpeptides 5’ sowie einer C-terminalen transmembran- und zytosolischen-Domäne eines MHC-Klasse-I-Moleküls am 3’ der Antigen-kodierenden Sequenz auf die Effizienz der MHC-Klasse-I und -II Präsentation. Wir konnten in vitro und in vivo nachweisen, dass diese Modifikation zu einer gesteigerten, simultanen Stimulation von antigenspezifischen CD4+ und CD8+ T-Zellen führt. Auf der Basis der optimierten Vektorkassetten etablierten wir die intranodale (i.n.) Transfektion von antigenpräsentierenden Zellen in der Maus. Dazu nutzten wir verschiedene Reportersysteme (eGFP-RNA, fluoreszensmarkierte RNA) und konnten zeigen, dass die intranodale Applikation von IVT-RNA zu selektiven Transfektion und Maturation lymphknotenresidenter DCs führt. Zur Untersuchung der immunologischen Effekte wurden in erster Linie auf Influenza-Hemagglutinin-A und Ovalbumin basierende Modellantigensysteme verwendet. Beide Antigene wurden als Antigen-MHC-Fusionskonstrukte genutzt. Als Responderzellen wurden TCR-transgene Lymphozyten verwendet, die MHC-Klasse-I oder -Klasse-II restringierte Epitope des Influenza-Hemagglutinin-A bzw. des Ovalbumin-Proteins erkennen. Wir konnten in vivo zeigen, dass die intranodale Immunisierung mit IVT-RNA zu einer effizienten Stimulation und Expansion von antigenspezifischen CD4+ und CD8+ T-Zellen in einer dosisabhängigen Weise führt. Funktionell konnte gezeigt werden, dass diese T-Zellen Zytokine sezernieren und zur Zytolyse befähigt sind. Wir waren in der Lage durch repetitive i.n. RNA Immunisierung ein ‚Priming’ CD8+ T-Zellen in naiven Mäusen sowohl gegen virale als auch gegen Tumor assoziierte Antigene zu erreichen. Die geprimten T-Zellen waren befähigt eine zytolytische Aktivität gegen mit spezifischem Peptid beladene Targetzellen zu generieren. Darüber hinaus waren wir in der Lage Gedächtnisszellen expandieren zu können. Abschließend konnten wir in Tumormodellen sowohl in prophylaktischen als auch in therapeutischen Experimenten zeigen dass die i.n. RNA Vakzination die Potenz zur Induktion einer anti-tumoralen Immunität besitzt.

Differences between RNA and DNA due to RNA editing in temporal lobe epilepsy

To investigate whether alterations in RNA editing (an enzymatic base-specific change to the RNA sequence during primary transcript formation from DNA) of neurotransmitter receptor genes and of transmembrane ion channel genes play a role in human temporal lobe epilepsy (TLE), this exploratory study analyzed 14 known cerebral editing sites in RNA extracted from the brain tissue of 41 patients who underwent surgery for mesial TLE, 23 with hippocampal sclerosis (MTLE+HS). Because intraoperatively sampled RNA cannot be obtained from healthy controls and the best feasible control is identically sampled RNA from patients with a clinically shorter history of epilepsy, the primary aim of the study was to assess the correlation between epilepsy duration and RNA editing in the homogenous group of MTLE+HS. At the functionally relevant I/V site of the voltage-gated potassium channel Kv1.1, an inverse correlation of RNA editing was found with epilepsy duration (r=-0.52, p=0.01) but not with patient age at surgery, suggesting a specific association with either the epileptic process itself or its antiepileptic medication history. No significant correlations were found between RNA editing and clinical parameters at other sites within glutamate receptor or serotonin 2C receptor gene transcripts. An "all-or-none" (≥95% or ≤5%) editing pattern at most or all sites was discovered in 2 patients. As a secondary part of the study, RNA editing was also analyzed as in the previous literature where up to now, few single editing sites were compared with differently obtained RNA from inhomogenous patient groups and autopsies, and by measuring editing changes in our mouse model. The present screening study is first to identify an editing site correlating with a clinical parameter, and to also provide an estimate of the possible effect size at other sites, which is a prerequisite for power analysis needed in planning future studies.

Isolation of an active gene and of two pseudogenes for mouse U7 small nuclear RNA

Three U7 RNA-related sequences were isolated from mouse genomic DNA libraries. Only one of the sequences completely matches the published mouse U7 RNA sequence, whereas the other two apparently represent pseudogenes. The matching sequence represents a functional gene, as it is expressed after microinjection into Xenopus laevis oocytes. Sequence variations of the conserved cis-acting 5' and 3' elements of U RNA genes may partly explain the low abundance of U7 RNA.

Structures of four human pseudogenes for U7 small nuclear RNA.

Four U7 RNA-related sequences were isolated from a human genomic DNA library. None of the sequences completely match the published human U7 RNA sequence and all of them contain features typical of reverse-transcribed pseudogenes.

Assembly of a catalytic unit for RNA microhelix aminoacylation using nonspecific RNA binding domains

An assembly of a catalytic unit for aminoacylation of an RNA microhelix is demonstrated here. This assembly may recapitulate a step in the historical development of tRNA synthetases. The class-defining domain of a tRNA synthetase is closely related to the primordial enzyme that catalyzed synthesis of aminoacyl adenylate. RNA binding elements are imagined to have been added so that early RNA substrates could be docked proximal to the activated amino acid. RNA microhelices that recapitulate the acceptor stem of modern tRNAs are potential examples of early substrates. In this work, we examined a fragment of Escherichia coli alanyl-tRNA synthetase, which catalyzes aminoacyl adenylate formation but is virtually inactive for catalysis of RNA microhelix aminoacylation. Fusion to the fragment of either of two unrelated nonspecific RNA binding domains activated microhelix aminoacylation. Although the fusion proteins lacked the RNA sequence specificity of the natural enzyme, their activity was within 1–2 kcal⋅mol−1 of a truncated alanyl-tRNA synthetase that has aminoacylation activity sufficient to sustain cell growth. These results show that, starting with an activity for adenylate synthesis, barriers are relatively low for building catalytic units for aminoacylation of RNA helices.

A general procedure for locating and analyzing protein-binding sequence motifs in nucleic acids

In the last decade, two tools, one drawn from information theory and the other from artificial neural networks, have proven particularly useful in many different areas of sequence analysis. The work presented herein indicates that these two approaches can be joined in a general fashion to produce a very powerful search engine that is capable of locating members of a given nucleic acid sequence family in either local or global sequence searches. This program can, in turn, be queried for its definition of the motif under investigation, ranking each base in context for its contribution to membership in the motif family. In principle, the method used can be applied to any binding motif, including both DNA and RNA sequence families, given sufficient family size.