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,549 records · Page 86Linked to original sources
PubMed2025

RNA sequencing offers new diagnostic opportunities in neurodevelopmental disorders: A systematic review.

PURPOSE: Transcriptomics by way of RNA sequencing (RNAseq) has emerged as a means to increase the diagnostic yield in genetic conditions. In this systematic review, we focus on the contribution of transcriptomics to improve the diagnostic yield in neurodevelopmental disorders. METHODS: We performed a systematic literature search in PubMed until January 2024, including articles describing diagnostic RNAseq on at least 1 individual with a primary neurodevelopmental phenotype. We extracted data on cohort size, phenotype, sample tissue, previously used diagnostic methods, added diagnostic yield of RNAseq, the use of control samples, and technical aspects of the RNA sequencing methodology. RESULTS: A total of 17 articles were eligible for inclusion in the systematic review. We found an average added diagnostic yield of 15.5% through RNA sequencing for individuals with neurodevelopmental disorders. There is heterogeneity in the tissue type, reported quality measures, and the computational pipeline. CONCLUSION: The significantly increased diagnostic yield demonstrates the value of this novel tool in the diagnostic setting of neurodevelopmental disorders. Our results offer an overview of common methodologies for RNAseq and allow us to formulate recommendations for genetic labs and clinicians when implementing RNAseq as a diagnostic tool. Lastly, we provide recommendations for future publications to increase transparency and reproducibility.

Humans↗
PubMed1997

Kinetic analysis of DNA and RNA strand transfer reactions catalyzed by vaccinia topoisomerase.

Vaccinia topoisomerase binds duplex DNA and forms a covalent DNA-(3'-phosphotyrosyl) protein adduct at the sequence 5'-CCCTT downward arrow. The enzyme reacts readily with a 36-mer CCCTT strand (DNA-p-RNA) composed of DNA 5' and RNA 3' of the scissile bond. However, a 36-mer composed of RNA 5' and DNA 3' of the scissile phosphate (RNA-p-DNA) is a poor substrate for covalent adduct formation. Vaccinia topoisomerase efficiently transfers covalently held CCCTT-containing DNA to 5'-OH-terminated RNA acceptors; the topoisomerase can therefore be used to tag the 5' end of RNA in vitro. Religation of the covalently bound CCCTT-containing DNA strand to a 5'-OH-terminated DNA acceptor is efficient and rapid (krel > 0.5 s-1), provided that the acceptor DNA is capable of base pairing to the noncleaved DNA strand of the topoisomerase-DNA donor complex. The rate of strand transfer to DNA is not detectably affected by base mismatches at the 5' nucleotide of the acceptor strand. Nucleotide deletions and insertions at the 5' end of the acceptor slow the rate of religation; the observed hierarchy of reaction rates is as follows: +1 insertion > -1 deletion > +2 insertion >> -2 deletion. These findings underscore the importance of a properly positioned 5'-OH terminus in transesterification reaction chemistry, but they also raise the possibility that topoisomerase may generate mutations by sealing DNA molecules with mispaired or unpaired ends.

Catalysis↗
PubMed2004

Identification of a novel MAGE D2 antisense RNA transcript in human tissues.

Using cDNA microarray analysis, we identified a cDNA clone, DD74, from primary human bronchial epithelial cells, which exhibits increased expression in vitro after treatment with all-trans retinoic acid. This clone corresponded to MAGE D2 mRNA, a gene previously identified to be upregulated in several cancer tissues. Surprisingly, in situ hybridization of lung tissue demonstrated positive hybridization signals with sense, but not antisense, MAGE D2-specific cRNA probes. Examination of several cell lines by Northern blot hybridization confirmed significant expression of two RNA bands. With strand-specific riboprobes, we identified a 2.0kb RNA transcript with the antisense probe as expected and identified a 4.1kb transcript by the sense probe. Further sequence analysis of the 4.1kb transcript revealed at least a 509 nucleotide sequence exactly complementary to the 2.0kb MAGE D2 mRNA sequence. This MAGE D2i sequence contains unique structural features not shared with those of previously described antisense transcripts. Identification of this transcript potentially has important implications for future studies examining MAGE D2 expression patterns in cancer and normal tissues.

Adaptor Proteins, Signal Transducing↗
PubMed2004

Alignment of RNA base pairing probability matrices.

MOTIVATION: Many classes of functional RNA molecules are characterized by highly conserved secondary structures but little detectable sequence similarity. Reliable multiple alignments can therefore be constructed only when the shared structural features are taken into account. Since multiple alignments are used as input for many subsequent methods of data analysis, structure-based alignments are an indispensable necessity in RNA bioinformatics. RESULTS: We present here a method to compute pairwise and progressive multiple alignments from the direct comparison of base pairing probability matrices. Instead of attempting to solve the folding and the alignment problem simultaneously as in the classical Sankoff's algorithm, we use McCaskill's approach to compute base pairing probability matrices which effectively incorporate the information on the energetics of each sequences. A novel, simplified variant of Sankoff's algorithms can then be employed to extract the maximum-weight common secondary structure and an associated alignment. AVAILABILITY: The programs pmcomp and pmmulti described in this contribution are implemented in Perl and can be downloaded together with the example datasets from http://www.tbi.univie.ac.at/RNA/PMcomp/. A web server is available at http://rna.tbi.univie.ac.at/cgi-bin/pmcgi.pl

Algorithms↗
PubMed1999

Sequence determination and phylogenetic analysis of RNA-dependent RNA polymerase (RdRp) of the porcine enterovirus 1 (PEV-1) Talfan strain.

The nucleotide sequence of the region including the RNA-dependent RNA polymerase (RdRp) of the porcine enterovirus 1 (PEV-1) Talfan strain was determined. Amino acid identities of Talfan RdRp with those of other picornaviruses were significantly lower than those seen among the viruses of the same genus of picornaviruses. The phylogenetic analysis of RdRp indicated that picornaviruses were divided into 8 clusters, in which Talfan was genetically distinct from other picornaviruses in contrast to the current taxonomy. In this study, we propose that PEV-1 Talfan be regarded as a new genus for the family Picorna- viridae.

Amino Acid Sequence↗
PubMed2026

An end-to-end computational framework for "Record-seq" transcriptional recording data.

MOTIVATION: Record-seq captures cumulative transcriptional activity over time in engineered Escherichia coli by integrating cellular RNA-derived spacer sequences into clustered regularly interspaced short palindromic repeats (CRISPR) arrays, which are read out by sequencing. Unlike the approximately uniform transcript sampling of RNA-seq, Record-seq records biological signal as spacers sampled by the CRISPR spacer acquisition machinery. Consequently, standard RNA-seq analysis strategies are not directly applicable, limiting sensitivity and interpretability. Our previous pipeline addressed these challenges only partially, retained inherited RNA-seq assumptions, and had limited algorithmic efficiency. RESULTS: Here, we present an end-to-end computational framework for Record-seq data. To address the primary computational bottleneck of spacer sequence extraction, we implemented a wavefront alignment approach for efficient quasi-local pattern matching, achieving an approximately 30-fold speedup. We introduce transcription unit-based feature counting as an alternative to gene-body quantification to better represent prokaryotic transcription and increase statistical power by capturing signal from untranslated regions, which are spacer acquisition hotspots. For downstream analyses, we incorporate multiple normalization strategies and a nonparametric differential expression testing framework designed for sparse datasets. Further, we analyze spacer acquisition patterns and train sequence-based neural models that predict acquisition propensity from genomic sequence and annotations, providing a framework for assessing whether acquisition rules generalize as Record-seq is extended to new microbial hosts. AVAILABILITY AND IMPLEMENTATION: The primary analysis workflow, the recoRdseq package, acquisition modeling repository, and relevant data are all linked at https://github.com/plattlab/Record-seq-Framework. Acquisition models and training data are on Zenodo at https://doi.org/10.5281/zenodo.18891434.

Escherichia coli↗
PubMed1997

Complete large subunit ribosomal RNA sequences from the heterokont algae Ochromonas danica, Nannochloropsis salina, and Tribonema aequale, and phylogenetic analysis.

The large subunit ribosomal RNA sequences from the heterokont algae Ochromonas danica, Nannochloropsis salina, and Tribonema aequale were determined. These sequences were combined with small subunit ribosomal RNA sequences in order to carry out a phylogenetic analysis based on neighbor-joining, maximum parsimony, and maximum likelihood methods. Our results indicate that heterokont fungi and heterokont algae each are monophyletic, and confirm that they together form a monophyletic group called "stramenopiles." Within the heterokont algae, the eustigmatophyte Nannochloropsis salina either clusters with the chrysophyte Ochromonas danica or forms a sister group to a cluster comprising the phaeophyte Scytosiphon lomentaria and the xanthophyte Tribonema aequale. The alveolates were identified as the closest relatives of the stramenopiles, but the exact order of divergence between the eukaryotic crown taxa could not be established with confidence.

Animals↗
PubMed2025

Permutation tests to assess sex differences in omics data.

It is common to sex-stratify analyses of omics data and to report effects as 'sex-specific' when they are significant in only one sex. However, when analysing hundreds or thousands of molecules, this approach will yield many spurious 'sex-specific' effects if not supported by significant interactions. I illustrate this problem using an RNA sequencing dataset showing almost no significant sex by treatment interactions, but where sex-stratified analyses yield hundreds of 'sex-specific' effects of treatment. These 'sex-specific' effects could be spurious or could be real but not show interactions due to low statistical power. To distinguish these possibilities, I describe permutation tests, which provide an intuitive way to determine if a pattern of observations differs from what would be expected due to chance. For this dataset, assigning sex at random often generates more 'sex-specific' effects than the real data, demonstrating that there is little evidence of sex differences. Next, I simulate an RNA sequencing dataset that includes genes modelled to have sex-specific effects of a condition. As expected, analysis of this simulated dataset yields both significant interactions and sex-specific effects in sex-stratified analyses. While stratified analyses detect a higher number of sex-specific effects than the analysis of interactions, they erroneously identify genes not modelled to show sex-specific effects more often than interactions. A permutation test confirms that the number of sex-specific effects observed in the simulated dataset is greater than expected due to chance. Permutation tests can be applied to omics studies of sex differences, simultaneously providing (i) a clear and simple demonstration of the problems of sex-stratified analyses, and (ii) additional evidence of sex-specific effects where these are present. R code is provided for permutations, simulations, and plots to visualize potential sex-specific effects, which can be adapted to other types of data.

Female↗
PubMed2004

Estrogen receptor-alpha messenger RNA variants that lack exon 5 or exon 7 are coexpressed with wild-type form in human endometrium during all phases of the menstrual cycle.

OBJECTIVE: We have assessed the expression levels of messenger RNA for estrogen receptor-alpha and splice variants lacking exon 5 or exon 7 that presumably exert dominant positive (splice variants lacking exon 5) and negative (splice variants lacking exon 7) effects, respectively, on estrogen responses in the human endometrium. STUDY DESIGN: This was a prospective study that was conducted at an academic community-based hospital. The patients, aged 18 to 40 years, underwent hysterectomy for benign gynecologic causes. Eighty-one endometrial specimens (46 proliferative, 35 secretory) were analyzed with the use of reverse transcription-polymerase chain reaction assay for the messenger RNA levels of estrogen receptor-alpha, and splice variants lacking exon 5 and exon 7. RESULTS: Wild-type estrogen receptor-alpha and splice variants splice variants lacking exon 5 and lacking exon 7 messenger RNAs were detected in all endometrial specimens throughout the menstrual cycle. In addition, a double-splice estrogen receptor-alpha messenger RNA variant (splice variants lacking exon 5 and exon 7) was detected at constant low levels of expression. Semiquantitative analysis showed higher levels of estrogen receptor-alpha messenger RNA in the early and mid proliferative endometrial phases than in late proliferative and secretory endometrium ( P <.05). The splice variant lacking exon 7 messenger RNA expression level was about 10-fold higher than the splice variant lacking exon 5 messenger RNA relative to wild-type estrogen receptor-alpha messenger RNA ( P <.001). The expression of splice variants lacking exon 5 compared with wild-type estrogen receptor-alpha messenger RNA is relatively constant throughout endometrial development. In contrast, an examination of the ratio of the levels of splice variants lacking exon 7 to wild-type estrogen receptor-alpha messenger RNA indicated a small, but significantly higher, splice variant lacking exon 7 level in the mid secretory phase (postovulatory days 5-8) than the mid proliferative and early secretory phases ( P <.05). CONCLUSION: We found no evidence of differential coexpression of the positive dominant estrogen receptor variant, splice variants lacking exon 5, with wild-type estrogen receptor-alpha. We did find that the dominant negative splice variant lacking exon 7 was slightly increased relative to wild-type estrogen receptor-alpha in the postovulatory phase. Future investigation is required to suggest the biologic significance of the observed increased relative expression of the splice variants lacking exon 7 in secretory endometrium and to determine the function of splice variants lacking exon 5 and splice variants lacking exon 7.

Alternative Splicing↗
PubMed2005

A novel strategy to identify the location of necessary and sufficient cis-acting regulatory mRNA elements in trypanosomes.

Expression of nearly all protein coding genes in trypanosomes is regulated post-transcriptionally, predominantly at the level of mRNA half-life. The identification of cis-acting elements involved in mRNA stability has been hindered by a lack of ability to screen for loss-of-regulation mutants. The method described in this article allows the region containing the necessary and sufficient elements within a mRNA to be identified and uses antibiotic resistance genes as both selectable markers and reporters. In the case of unstable mRNAs, the strategy can be extended by performing a screen for spontaneous loss-of-function mutants in regulatory parts of a mRNA. The method was validated by using the GPI-PLC mRNA, which is unstable in procyclic form trypanosomes and showed that the 3'UTR of the GPI-PLC mRNA contains all elements required for developmentally regulated instability. Loss-of-instability mutants all contained deletions within the 2300-nucleotide-long 3'UTR, and their analysis showed that a deletion including the last 800 nt of the gene stabilized the mRNA. The method is nonpresumptive, allows far more rapid screening for cis-elements than existing procedures, and has the advantage of identifying functional mutants. It is applicable to all eukaryotes using polycistronic transcription.

3' Untranslated Regions↗
PubMed2004

Computational prediction of RNA editing sites.

MOTIVATION: Some organisms edit their messenger RNA resulting in differences between the genomic sequence for a gene and the corresponding messenger RNA sequence. This difference complicates experimental and computational attempts to find and study genes in organisms with RNA editing even if the full genomic sequence is known. Nevertheless, knowledge of these editing sites is crucial for understanding the editing machinery of these organisms. RESULTS: We present a computational technique that predicts the position of editing sites in the genomic sequence. It uses a statistical approach drawing on the protein sequences of related genes and general features of editing sites of the organism. We apply the method to the mitochondrion of the slime mold Physarum polycephalum. It correctly predicts over 90% of the amino acids and over 70% of the editing sites.

Algorithms↗
PubMed2005

Sequence characteristics of functional siRNAs.

RNA interference in mammalian cells is actively used to conduct genetic screens, to identify and to validate targets, and to elucidate regulators and modifiers of cellular pathways. To ensure the specificity and efficacy of the active 21mer siRNA molecules, it is pertinent to develop a strategy for their rational design. Here we show that most functional siRNAs have characteristic sequence features. We tested 601 siRNAs targeting one exogenous and three endogenous genes. The efficacy of the siRNAs was determined at the protein level. Using a decision tree algorithm in combination with information analysis, our analyses revealed four sets of rules with a mean knockdown efficacy ranging from 60% to 73%. (To distinguish between percentages used to describe the quality of an siRNA and the percentages used to describe parts of data sets we underlined the former throughout this paper.) The best rule comprises an A/U at positions 10 and 19, a G/C at position 1, and more than three A/Us between positions 13 and 19, in the sense strand of the siRNA sequence. Using these rules, there is a 99.9% chance of designing an effective siRNA in a set of three with more than 50% knockdown efficiency in a biological readout.

Base Sequence↗
PubMed2003

[Optimization of T7-based RNA amplification system for cDNA microarray].

cDNA microarrays are powerful parallel tools for gene expression profiling analysis, which help us to understand the molecular mechanism of diseases and to identify potential targets for therapeutic intervention. However, their broader application are hampered by the large amount of RNA required: up to 200 microg of total RNA or 5 microg of mRNA for one chip, making analysis of small samples difficult. In this work, combined with a template switching effect, the T7 RNA linear amplification procedure was optimized, providing multiple copies of anti-sense RNA with reduced sample inputs: no more than 3 microg total RNA for one chip. Using the same RNA sample in all cases, the anti-sense RNA labeling method was compared with standard total RNA and mRNA methods by two sets of self-comparison experiments. Furthermore, the three methods in profiling analysis were compared with the same pair RNA samples. All the results indicated that the fidelity, reproducibility, and reliability showed no significant difference with conventional total RNA or mRNA microarrays.

DNA-Directed RNA Polymerases↗
PubMed2005

Analysis of hypoxia-related gene expression in sarcomas and effect of hypoxia on RNA interference of vascular endothelial cell growth factor A.

Vascular endothelial cell growth factor A (VEGF-A) and hypoxia play important roles in tumor angiogenesis. VEGF-A gene expression is up-regulated in tumors under hypoxic conditions, yet it is unclear how such up-regulation will affect the efficacy of RNA interference strategies targeting VEGF-A. Four potential short interfering RNA (siRNA) sequences for the VEGF-A gene were cloned into expression plasmids and transfected into HT1080 human fibrosarcoma cells. Stable transfection of these plasmids decreased VEGF-A mRNA levels and protein secretion by up to 99%. Our analysis of >100 hypoxia-related genes using oligonucleotide microarrays of 38 human sarcoma samples and 14 normal tissues identified distinctly different patterns of expression between sarcomas and normal tissues as assessed by hierarchical clustering analysis. Numerous hypoxia-related genes were significantly up-regulated in sarcomas including hypoxia-inducible factor 1alpha (HIF-1alpha). Exposure of wild-type HT1080 cells to 1% hypoxia resulted in HIF-1alpha up-regulation and a 74% increase in VEGF-A secretion as compared with secretion under normoxic conditions. Surprisingly, stable cell lines expressing VEGF-A siRNAs silenced VEGF-A expression equally well in hypoxia and normoxia. S.c. injection of cells with VEGF-A siRNAs into athymic nude mice led to slower-growing tumors, decreased blood vessel density, and greater apoptosis when compared with controls. Immunofluorescence analysis of tumor sections revealed areas of HIF-1alpha nuclear expression, suggesting areas of hypoxia, in both control tumors and VEGF-suppressed tumors. We conclude that hypoxia plays an important role in human sarcomas but hypoxic up-regulation of VEGF-A expression does not attenuate the efficacy of VEGF-A RNA interference.

Animals↗
PubMed2003

Sequence confirmation of the EWS-WT1 fusion gene transcript in the peritoneal effusion of a patient with desmoplastic small round cell tumor.

Desmoplastic small round cell tumor (DSRCT) is a rare undifferentiated neoplasm. The prognosis is poor, even if therapy is instituted promptly, and thus it is important to differentiate it from other histologically and cytologically similar-looking malignancies of the young adult. We present a case of DSRCT in a 17-yr-old male with disseminated peritoneal disease and peritoneal effusion. The cytology sample showed a malignant small round cell tumor, the classical cytological features of DSRCT, and immunohistochemistry performed in the prepared cell block exhibited an antibody expression profile in keeping with DSRCT. Further material from the effusion was prepared for RNA extraction, following which a reverse-transcriptase polymerase chain reaction (RT-PCR) and sequencing of the t(11;22)(p13;q11 or q12) were carried out. The result showed the presence of the reciprocal translocation and thus confirmed the diagnosis of DSRCT. This case shows how molecular techniques (including sequencing) can be applied to cytology in clarifying and confirming certain difficult diagnosis of undifferentiated neoplasms, DSRCT in this particular case.

Adolescent↗
PubMed1998

Rhodopsins from three frog and toad species: sequences and functional comparisons.

The frequency of thermal 'dark events' in the membrane current of rhodopsin rods of the bullfrog, Rana catesbeiana, is considerably lower than observed in rods of two toad species, even though all three rhodopsins have approximately the same absorbance characteristics. In order to map amino acid substitutions possibly associated with thermal stability in the genus Rana, the cDNA's coding for the rhodopsins of Bufo bufo, B. marinus and R. temporaria were sequenced and the conceptually translated protein sequences aligned to the previously sequenced rhodopsins of R. catesbeiana, R. pipiens and Xenopus laevis. Across the six anuran species studied, there are sixteen non-conserved substitutions and six changes that include gain or loss of a hydroxyl group. Serine or threonine at position 220 is unique to the three Rana species, phenylalanine at position 270 is unique to all three Ranas and to X. laevis, and phenylalanine at position 274 is unique to both species of the genus Bufo. This investigation produces a list of substitutions that are candidates for future studies of thermal stability. In addition, a number of amino acids are identified that apparently do not influence absorbance characteristics, at least not cumulatively.

Amino Acid Sequence↗
PubMed2006

Construction of shRNA targeted to the rat angiotensin II type 1 receptors and its RNAi in cytoplasma.

The expression vector of shRNA targeted to the rat angiotensin II receptor gene was constructed and the efficacy of siRNAs to modulate the expression of target gene in the in vitro cultured mammalian cells was investigated for antihypertensive therapy in spontaneous hypertensive rat (SHR) at post-transcriptional level. The sense and antisense RNA oligonucleotides strands targeting angiotensin II receptor mRNA were synthesized individually according to the sequence of the rat angiotensin II receptor. For preparation of duplexes, sense- and antisense-stranded oligonucleotides were mixed and annealed, and the annealed duplexes were cloned into the pGenesil-1 vector. . The rat glioma cells were transfected with constructed pGenesil-1-shRNA plasmid and scrambled plasmid. The cultured cells were collected at different phases. RT-PCR and Western blot were performed. The AT1 mRNA and protein levels behaved ultimately same. Compared to control after 48 h, AT1 mRNA levels were decreased to 35.5% +/- 3.0%, and the levels reached their lowest point after 72 h (20.7% +/- 4% of control). At 24 and 48 h, AT1 protein was reduced to 46.9% +/- 4. 2% and 36.98% +/- 3.7% respectively compared to control and a maximum reduction was observed after 72 h of incubation (28.1% +/- 4% compared to controls). Plasmid-based shRNA expression systems targeted against the rat angiotensin II receptor gene were generated successfully. The shRNAs with a 22-nt stem and a short loop were cleaved into small interfering dsRNA (siRNA) by the Dicer. The in vitro transcribed siRNA enables the effective silencing of gene expression to the target mRNA and leads to effective inhibition of translation of proteins and will be lay the foundation of application of gene silencing technology to hypertensive rats.

Angiotensin II↗