Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “RNA sequencing analysis”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,351 records · Page 75Linked to original sources

STR2: a structure to string approach for locating G-box riboswitch shapes in pre-selected genes.

Traditional sequence-based search methods such as BLAST and FASTA can be used to identify sequence similarities. Recently, there is a growing interest in performing RNA shape similarity searches inside selected genes to locate RNA structure motifs that are known to possess functionally important roles. For example, in the newly discovered RNA genetic control elements called "riboswitches", the box domain is known to be highly conserved among various bacterial species in both its nucleotide composition and shape. However, in non-bacterial species, shape conservation is likely to become more important than sequence conservation when searching for riboswitch patterns. For this purpose, we present an approach tailored for detecting RNA shape similarities. We extend the Structure to String (ST R2) method that was initially proposed to locate shape similarities in proteins to identify predicted secondary structures of RNAs. The ST R2 for RNAs is a translation of a secondary structure to a string of characters, after which known sequence-based search algorithms with an efficient implementation are being used. We validate that the ST R2 succeeds to locate G-box riboswitches in prokaryotes, as expected. Subsequently we show running examples when attempting to detect G-box riboswitch candidates in eukaryotes.

Algorithms↗

Exact determination of UV-induced crosslinks in 16S ribosomal RNA in 30S ribosomal subunits.

Escherichia coli 30S ribosomal subunits were UV-irradiated to induce intramolecular crosslinks in the 16S rRNA. Intact 16S rRNA was purified and subjected to gel electrophoresis, under denaturing conditions, to separate molecules on the basis of the crosslinked loop size. Molecules separated this way were enriched for specific crosslinks and could be analyzed by the reverse transcription arrest assay to determine exact crosslinking sites. Thirteen crosslinking sites have been identified at single nucleotide resolution. Of these, eight are within or adjacent to secondary structure elements: one of these (C582 x G760) involves an interaction between nucleotides within an interior loop, one (C1402 x X1501) involves an interaction between nucleosides in adjacent base pairs, and the others involve interactions between nucleotides that are within junction regions (A441 x G494, U562 x U884, C934 x U1345, and U991 x U1212) or are interactions between nucleotides (C54 x A353 and U1052 x C1200) that somehow cross known base pairs. Five other crosslinks connect sites distant in the secondary structure and provide global constraints for the arrangement of RNA regions within RNA domains I and II (U244 x G894, G894 x A1468, C967 x C1400) and within domain III (U1126 x C1281 and A1093 x G1182). These crosslinks, known at single-nucleotide resolution, are useful in the prediction of local RNA regions, as well as the global structure.

Base Sequence↗

RT-PCR in foot-and-mouth disease diagnosis.

A RT-PCR assay for the specific detection of RNA sequences from foot-and-mouth disease virus (FMDV) has been developed. The procedure permits also the detection of sequences that correlate with established FMDV serotypes. A computer program that allows selection of genotype-specific primers for RT-PCR amplification was used for the identification of FMDV specific sequences for PCR amplification on RNA replicase (3D) gene regions. Specific, rapid and highly sensitive detection was achieved for a wide collection of RNA samples from FMDV types C, A and O, either purified from tissue culture or extracted from lesions of infected animals. Similarly, serotype-specific primers were designed to amplify the carboxy-terminal end of the VP1 gene of FMDV types either C, A or O. The results of PCR amplification of different FMDV RNAs using type-specific primers are in agreement with the serological typing of the corresponding viruses. A combination of this approach with a simplified sample processing, carried out following direct adsorption of viral suspensions to microtiter plates, provides a rapid, reliable method of viral diagnosis.

Animals↗

A kissing-loop interaction in a hammerhead viroid RNA critical for its in vitro folding and in vivo viability.

Chrysanthemum chlorotic mottle viroid (CChMVd) RNA (398-401 nucleotides) can form hammerhead ribozymes that play a functional role in its replication through a rolling-circle mechanism. In contrast to most other viroids, which adopt rod-like or quasi-rod-like secondary structures of minimal free energy, the computer-predicted conformations of CChMVd and Peach latent mosaic viroid (PLMVd) RNAs are branched. Moreover, the covariations found in a number of natural CChMVd variants support that the same or a closely related conformation exists in vivo. Here we report that the CChMVd natural variability also supports that the branched conformation is additionally stabilized by a kissing-loop interaction resembling another one proposed in PLMVd from in vitro assays. Moreover, site-directed mutagenesis combined with bioassays and progeny analysis showed that: (1) single CChMVd mutants affecting the kissing loops had low or no infectivity at all, whereas infectivity was recovered in double mutants restoring the interaction; (2) mutations affecting the structure of the regions adjacent to the kissing loops reverted to wild type or led to rearranged stems, also supporting their interaction; and (3) the interchange between 4 nucleotides of each of the two kissing loops generated a viable CChMVd variant with eight mutations. PAGE analysis under denaturing and nondenaturing conditions revealed that the kissing-loop interaction determines proper in vitro folding of CChMVd RNA. Preservation of a similar kissing-loop interaction in two hammerhead viroids with an overall low sequence similarity suggests that it facilitates in vivo the adoption and stabilization of a compact folding critical for viroid viability.

Asteraceae↗

Dimeric transfer RNA precursors in S. pombe.

Sequence analysis of a Schizosaccharomyces pombe DNA fragment revealed two tRNA coding regions separated by a seven nucleotide spacer. the 5'-proximal tRNA gene encodes a tRNAUCGSer sequence, which is interrupted by a 16 nucleotide intron at the 3' side of the base adjacent to the anticodon. The second tRNA gene encodes an initiator tRNAMet sequence. This DNA fragment, cloned into pBR322, was used as template for in vitro transcription in a nuclear extract of Xenopus oocytes. The tRNA genes were transcribed into one RNA precursor which contained both tRNA sequences. The primary transcription product initiates with pppG, contains a 9 nucleotide leader sequence, a 16 nucleotide intron, a 7 nucleotide spacer between the two tRNA molecules and an 8--9 nucleotide trailer sequence. RNA initiation was only observed upstream of the 5'-proximal tRNASer. We used RNA analysis to establish a sequence of the enzymatic steps of tRNA maturation in the nuclear extract. The first step in processing the dimeric precursor is an endonuclease cleavage which generates the mature 5' end of the tRNAMet. Further steps include the removal of the flanking sequences and addition of the CCAOH 3' terminus. The last step is the splicing of the tRNASer precursor to remove the intervening sequence.

Ascomycota↗

Molecular localisation of a G-protein mRNA using differential display and in situ hybridization.

Defining molecular repertoires within virally infected tissues of the nervous system may provide insight into the pathogenesis of, and immunity to, neurotropic viruses. Here we report the application of such a method, namely mRNA differential display (DD), to the identification of mRNAs that are expressed at different levels in herpes simplex virus (HSV) infected nervous tissue from immunocompetent and CD8(+) lymphocyte depleted mice. Small amounts of input RNA can be used by DD, making the method ideal for experiments based on murine sensory ganglia (DRG), which on average yield less than 0.5 microg of total RNA. In the current work, DD facilitated the identification of a mRNA whose abundance in HSV-infected ganglia, based on Northern blot analysis, was reduced in mice depleted of CD8(+) cells. The cloned product of this mRNA was of particular interest to our research as sequence data strongly suggested that it represented the murine homologue of the alpha chain of a G protein termed Golf. This G protein had not previously been reported from dorsal root ganglial tissue. RT-PCR confirmed the presence of Golf in DRG and in situ hybridization studies localised this molecule to primary sensory neurons. These data indicate that DD is sufficiently robust to be applied to the study of virus pathogenesis within the nervous system.

Animals↗

Enhancing pan-cancer spatial transcriptomics at single-cell resolution with stPainter.

Subcellular spatial transcriptomics can resolve tissue architecture at cellular scale, but sparse gene panels and limited detection sensitivity constrain downstream analysis. Existing enhancement methods often require tissue-matched single-cell RNA sequencing (scRNA-seq) references and dataset-specific retraining. Here we show that stPainter, a conditional generative model pretrained on a pan-cancer scRNA-seq atlas, can enhance spatial transcriptomics data without matched references or retraining. Using a latent diffusion architecture guided by Stochastic Differential Equations (SDE), stPainter reconstructs expanded expression profiles from sparse measurements and produces latent representations for clustering and cell-state analysis. When we apply stPainter upon 6 spatial transcriptomics datasets of different cancer types, we demonstrate that our model empowers downstream biological analyses, including fine-grained subpopulation clustering and pathway enrichment. Comparison with spatially resolved proteomics (CODEX) provided independent support for regional agreement between imputed cellular compositions and protein-level tissue organization. These results establish stPainter as a scalable approach for analyzing tumor microenvironments without auxiliary sequencing data.

Spatial Transcriptomics↗

Cloning and expression of a cDNA copy of the viral K28 killer toxin gene in yeast.

The killer toxin K28, secreted by certain killer strains of the yeast Saccharomyces cerevisiae is genetically encoded by a 1.9 kb double-stranded RNA, M-dsRNA (M28), that is present within the cell as a cytoplasmically inherited virus-like particle (VLP). For stable maintenance and replication, M28-VLPs depend on a second dsRNA virus (LA), which has been shown to encode the major capsid protein (cap) and a capsid-polymerase fusion protein (cap-pol) that provides the toxin-coding M-satellites with their transcription and replicase functions. K28 toxin-coding M28-VLPs were isolated, purified and used in vitro for the synthesis of the single-stranded M28 transcript, which was shown to be of plus strand polarity and to bind to oligo(dT)-cellulose, indicating that M28(+)ssRNA contains an internal A-rich tract. Strand separation of the 1.9 kb M28-dsRNA and direct RNA sequencing of its 3' ends was performed in order to obtain specific DNA oligonucleotides that could be used as primers for cDNA synthesis. The nucleotide sequence of the toxin-coding M28-cDNA identified a single open reading frame (ORF) coding for a polypeptide of 345 amino acids, which contained two potential Kex2p/Kex1p processing sites and three potential sites for protein N-glycosylation. The toxin-coding cDNA was cloned and expressed in sensitive non-killer strains under the control of the yeast PGK promoter. Upon transformation, this construct conferred the complete K28 phenotype, demonstrating that both toxin and immunity determinants are contained within the cloned cDNA. In vitro translational analysis of the M28(+)ssRNA in vitro transcript identified the primary gene product of M28 as a K28 preprotoxin of 38 kDa (M-p38).

Amino Acid Sequence↗

Development of a microarray assay that measures hybridization stoichiometry in moles.

Microarray data is most useful when it can be compared with other genetic detection technologies. In this report, we designed a microarray assay format that transforms raw data into a defined scientific unit (i.e., moles) by measuring the amount of array feature present and the cDNA sequence hybridized. This study profiles a mouse reference universal RNA sample on a microarray consisting of PCR products. In measuring array features, a labeled DNA sequence was designed that hybridizes to a conserved sequence that is present in every array feature. To measure the amount of cDNA sample hybridized, the RNA sample was processed to ensure consistent dye to DNA ratio for every labeled target cDNA molecule, using labeled branched dendrimers rather than by incorporation. A dye printing assay was then performed in order to correlate molecules of cyanine dye to signal intensity. We demonstrate that by using this microarray assay design, raw data can be transformed into defined scientific units, which will facilitate interpretation of other experiments, such as data deposited at the Gene Expression Omnibus and ArrayExpress.

Animals↗

Murine gammaherpesvirus 68 encodes tRNA-like sequences which are expressed during latency.

Murine gammaherpesvirus 68 (MHV-68) is a virus of wild rodents and is a convenient small animal model for studies of gammaherpesvirus pathogenesis. We have sequenced 6162 bp at the left end of the MHV-68 genome and identified two unique open reading frames (ORFs) (ORF2 and ORF3) and an ORF (ORF1) which displays similarity to poxvirus members of the serpin family. Interspersed with the ORFs is a family of eight novel tRNA-like sequences sharing tRNA-like predicted secondary structures and RNA polymerase III promoter elements. These sequences are expressed to high levels during lytic infection and are processed into mature tRNAs with post-transcriptionally added 3' CCA termini, indicating their recognition as tRNAs by cellular machinery. Acidic Northern analysis of four tRNAs tested has demonstrated that they are not aminoacylated by aminoacyl-tRNA synthetases present in the infected cell. Thus, it is currently unclear what biological function these uncharged viral tRNA-like sequences may fulfil. In situ hybridization analysis has shown that in addition to being expressed within productively infected tissues during acute stages of infection, the tRNA-like sequences are abundantly expressed within splenic germinal centres of latently infected mice. Therefore, the MHV-68 viral tRNAs represent a marker for latent infection and constitute the first report of tRNA-like sequences encoded by a virus of eukaryotes.

Acylation↗

Robust analysis of 5'-transcript ends (5'-RATE): a novel technique for transcriptome analysis and genome annotation.

Complicated cloning procedures and the high cost of sequencing have inhibited the wide application of serial analysis of gene expression and massively parallel signature sequencing for genome-wide transcriptome profiling of complex genomes. Here we describe a new method called robust analysis of 5'-transcript ends (5'-RATE) for rapid and cost-effective isolation of long 5' transcript ends (approximately 80 bp). It consists of three major steps including 5'-oligocapping of mRNA, NlaIII tag and ditag generation, and pyrosequencing of NlaIII tags. Complicated steps, such as purification and cloning of concatemers, colony picking and plasmid DNA purification, are eliminated and the conventional Sanger sequencing method is replaced with the newly developed pyrosequencing method. Sequence analysis of a maize 5'-RATE library revealed complex alternative transcription start sites and a 5' poly(A) tail in maize transcripts. Our results demonstrate that 5'-RATE is a simple, fast and cost-effective method for transcriptome analysis and genome annotation of complex genomes.

5' Untranslated Regions↗

Gene discovery in the embryonic chick retina.

PURPOSE: The chick embryo is a powerful model system for the study of retinal development. However, analysis of gene expression in the chick retina has lagged behind biological studies. The purpose of this study was to identity and characterize genes expressed in the chick embryo retina as candidate molecules involved in the development and function of photoreceptors and other retinal cell types. METHODS: RNA from embryonic day (ED) 18 White Leghorn chick embryo retinae was used to generate an oligo dT-primed cDNA library. Bacterial colonies representing five thousand individual clones were arrayed onto nylon membranes using a microarray robot. Replicate membranes were hybridized with cDNA probes synthesized from ED 18 retina, brain and liver. Clones that appeared preferentially expressed in retina were identified by homology searches, and their spatial and temporal expression patterns were analyzed by in situ hybridization. RESULTS: Two hundred and seventy-two clones were identified. Approximately forty percent of the clones represented potential novel genes, including ESTs, hypothetical proteins and clones with no assigned identities. Furthermore, many genes were identified that are the putative chick orthologues of genes cloned from other species. We determined the expression pattern of several clones for which sequence homologies suggested possible roles in transcriptional regulation, apoptosis or intercellular signaling. Their corresponding mRNAs were expressed in the embryonic retina in topographically specific, developmentally regulated patterns. CONCLUSIONS: We identified and characterized genes in the chick embryo retina using a combination of microarray analysis and in situ hybridization. Analysis of the expression patterns suggests involvement of several of these genes in key events during embryogenesis.

Animals↗

The diversity of rabies virus in Russia demonstrated by anti-nucleocapsid monoclonal antibody application and limited gene sequencing.

A comparative analysis of anti-nucleocapsid monoclonal antibody (N-mAb) reaction profiles and DNA sequences was performed on 97 selected rabies virus (RABV) isolates from Russia and neighbouring states. A panel of 73 N-mAb from the Wistar Institute (Philadelphia, PA, USA), Veterinary Laboratories Agency (Weybridge, UK) and P-41 (Tübingen, Germany) was used. The sequence of a 400 bp fragment of the N gene was generated for all available isolates and a phylogenetic analysis of the fragment was carried out. The results of genetic and antigenic typing were in concordance, except for a few deviations. Three RABV variants with distinct antigenic patterns corresponded to well characterised genetic groups. The phylogenetic analysis also allowed the discrimination of RABVs with similar antigenic profiles originating from distant geographical locations in Russia and other states of the former Soviet Union. In some instances, more than one specific antigenic variant was detected within a phylogenetic group. A number of antigenic and genetic variants were associated with a specific host species (such as red fox, raccoon dog, wolf or polar fox) in different parts of Russia. A single, antigenically divergent, P-41-positive, rodent isolate from Siberia was described.

Animals↗

An application for assessing quality of RNA hybridized to Affymetrix GeneChips.

UNLABELLED: This paper describes a stand-alone application for estimating the 3' to 5' ratio by fitting a mixed effects model to the interior pixel intensities of perfect match probes for selected control probe sets from an Affymetrix *.DAT file. The effectiveness of this method was demonstrated previously by an application of the method to two microarray datasets for which external verification of RNA quality was known. This application provides a more objective assessment of sample quality in that both a point estimate and 95% confidence interval about the 3' to 5' ratio are provided. AVAILABILITY: The software and installation instructions are freely available for download at http://www.people.vcu.edu/~kjarcher/Research/Data.htm

Algorithms↗

Distribution of SR protein exonic splicing enhancer motifs in human protein-coding genes.

Exonic splicing enhancers (ESEs) are pre-mRNA cis-acting elements required for splice-site recognition. We previously developed a web-based program called ESEfinder that scores any sequence for the presence of ESE motifs recognized by the human SR proteins SF2/ASF, SRp40, SRp55 and SC35 (http://rulai.cshl.edu/tools/ESE/). Using ESEfinder, we have undertaken a large-scale analysis of ESE motif distribution in human protein-coding genes. Significantly higher frequencies of ESE motifs were observed in constitutive internal protein-coding exons, compared with both their flanking intronic regions and with pseudo exons. Statistical analysis of ESE motif frequency distributions revealed a complex relationship between splice-site strength and increased or decreased frequencies of particular SR protein motifs. Comparison of constitutively and alternatively spliced exons demonstrated slightly weaker splice-site scores, as well as significantly fewer ESE motifs, in the alternatively spliced group. Our results underline the importance of ESE-mediated SR protein function in the process of exon definition, in the context of both constitutive splicing and regulated alternative splicing.

Alternative Splicing↗

Comparative analysis of orthologous eukaryotic mRNAs: potential hidden functional signals.

Sequencing of multiple, nearly complete eukaryotic genomes creates opportunities for detecting previously unnoticed, subtle functional signals in non-coding regions. A genome-wide comparative analysis of orthologous sets of mammalian and yeast mRNAs revealed distinct patterns of evolutionary conservation at the boundaries of the untranslated regions (UTRs) and the coding region (CDS). Elevated sequence conservation was detected in approximately 30 nt regions around the start codon. There seems to be a complementary relationship between sequence conservation in the approximately 30 nt regions of the 5'-UTR immediately upstream of the start codon and that in the synonymous positions of the 5'-terminal 30 nt of the CDS: in mammalian mRNAs, the 5'-UTR shows a greater conservation than the CDS, whereas the opposite trend holds for yeast mRNAs. Unexpectedly, a approximately 30 nt region downstream of the stop codon shows a substantially lower level of sequence conservation than the downstream portions of the 3'-UTRs. However, the sequence in this poorly conserved 30 nt portion of the 3'-UTR is non-random in that it has a higher GC content than the rest of the UTR. It is hypothesized that the elevated sequence conservation in the region immediately upstream of the start codon is related to the requirement for initiation factor binding during pre-initiation ribosomal scanning. In contrast, the poorly conserved region downstream of the stop codon could be involved in the post- termination scanning and dissociation of the ribosomes from the mRNA, which requires only the mRNA-ribosome interaction. Additionally, it was found that the choice of the stop codon in mammals, but not in yeasts, and the context in the immediate vicinity of the stop codons in both mammals and yeasts are subject to strong selection. Thus, genome-wide analysis of orthologous gene sets allows detection of previously unrecognized patterns of sequence conservation, which are likely to reflect hidden functional signals, such as ribosomal filters that could regulate translation by modulating the interaction between the mRNA and ribosomes.

3' Untranslated Regions↗

The RNA world of the nucleolus: two major families of small RNAs defined by different box elements with related functions.

We have discovered that all known yeast and vertebrate small nucleolar RNAs (snoRNAs), except for the MRP/7-2 RNA, fall into two major classes. One class is defined by conserved boxes C and D and the other by a novel element: a consensus ACA triplet positioned 3 nt before the 3' end of the RNA. A role for the ACA box is snoRNA stability has been established by mutational analysis of a yeast ACA snoRNA (snR 11). Full function of the box depends on the integrity of an adjacent upstream stem. All members of the yeast ACA family are associated with the GAR1 protein. Binding of this or another common small nucleolar ribonucleoprotein particle protein is predicted to be a critical entry point to snoRNA posttranscriptional life, including precise formation of the snoRNA 3' end.

Animals↗

Verification of bluetongue virus S9 segment nucleotide sequences.

During the course of our bluetongue virus (BTV) nucleic acid sequence investigations, conflicts among United States (US) prototype BTV S9 genome segment sequences deposited in GenBank were noted. In order to rectify these inter-laboratory discrepancies, the S9 segments of Arthropod-borne Animal Diseases Research Laboratory (ABADRL)-stored US prototype BTV 2, BTV 10, BTV 11, BTV 13, and BTV 17 isolates were resequenced. Our S9 sequences, determined by direct sequencing of full-length reverse transcriptase-polymerase chain reaction (RT-PCR) generated amplicons, shared 99% or greater nucleotide identity with one or more respective S9 sequences previously reported. Possible sources of remaining unsupported US prototype BTV S9 sequences were evaluated by amplifying and sequencing the S9 segments of BTV 2 Ona A strain, South African (SA) prototype BTV 1, BTV 2, and BTV 4 strains, and the North American (NA) prototype epizootic hemorrhagic disease virus (EHDV) serotype 2 (Alberta) strain. Comparative analysis using these S9 sequences, as well as sequences of US BTV 2 field isolates, identified potential contributors to inter-laboratory sequence disagreements.

Base Sequence↗