As of May 31, 2023, we have updated our Code of Conduct.

Questions tagged [rna]

Use this tag to refer questions that are related to rna sequence.

Filter by
Sorted by
Tagged with
12 votes
2 answers
491 views

What to use to edit RNA alignments?

I have many alignments from Rfam Database, and I would like to edit them. I saw that many tools are used for Protein sequence alignments, but there is something specific to edit RNA alignments ? e....
Peter's user avatar
  • 343
11 votes
3 answers
1k views

How to extract RNA sequence and secondary structure restrains from a PDB file

I'm trying to find a programmatic way to automatically extract the following information from a PDB file: RNA sequence Secondary structure restraints in bracket format, e.g. ...
Peter's user avatar
  • 343
8 votes
1 answer
143 views

Coronavirus RNA structures?

Is there anything known about the RNA structures of coronaviruses? More specifically - do they have any interesting known structures in the translatable region, like RRE of HIV or the double loops in ...
Roger Vadim's user avatar
6 votes
2 answers
4k views

5' and 3' bias in Rna-seq data

I'm working with rna-seq samples. I see 5' bias and also 3' bias in the per-base sequence content plot. From this link I see that the bias at the start of the sequences appears to be the result of ...
stack_learner's user avatar
6 votes
2 answers
298 views

Estimate the length of poly-A tails from randomly-primed RNAseq data

So a poly-A tail is a long chain of adenine nucleotides that is added to a messenger RNA (mRNA) molecule during RNA processing to increase the stability of the molecule. For my project, I would like ...
ShanZhengYang's user avatar
5 votes
1 answer
86 views

Why does the SARS-Cov2 genome has letter t [duplicate]

ATTAAAGGTT TATACCTTCC CAGGTAACAA ACCAACCAAC TTTCGAT... is part of the 5'UTR of genome of an RNA virus SARS-Cov-2. RNA contains letters ...
Pavel Pelan's user avatar
5 votes
2 answers
1k views

Reads mapped to exonic, intronic and intergenic regions

After the alignment step I checked the rnaseq metrics of all the samples. Among 40 samples three samples show high percentage of reads mapped to intronic regions. What could be the reason? ...
stack_learner's user avatar
5 votes
1 answer
95 views

How to compute the chance of failing to detect a gene given the detection limit of a protocol

In Shapiro et al., when discussing about loss of molecules as source of error in single-cell sequencing, it is written that: Another source of error is losses, which can be severe. The detection ...
gc5's user avatar
  • 1,773
5 votes
1 answer
110 views

How can I dock a protein to a nucleic acid?

I have a protein of interest and I would like to now how it interacts with RNA. I have structures of both molecules. What tool can I use?
John Deo's user avatar
5 votes
2 answers
92 views

Integrating bulk RNA-Seq data with different sequencing depths and from different sources

I am attempting to integrate different bulk RNA-Seq datasets. While this is not ideal, I'm trying to reduce the technical variability in these datasets by using data generated by similar protocols (...
CodingSquirrel's user avatar
5 votes
0 answers
47 views

Parse RNA variant effect annotations ("r." format)

I've got annotations for splicing variants in a format like this (this is one variant): Variant: NM_004092.3:c.88+5G>A Effect: Retention; r.87_88ins1_88+10:p.(Ala31Glufs*23) I want to extract ...
Ron's user avatar
  • 51
4 votes
1 answer
184 views

Why is ribosomal RNA difficult to remove even with Poly(A) selection?

In this answer (actually in a comment), it is stated that: As you've noticed from your own analysis, the ribosomal genes have quite variable expression across cells. They're expressed everywhere, ...
gc5's user avatar
  • 1,773
4 votes
1 answer
44 views

How to calculate the Maximum Ladder Distance (MLD) of ssRNA from Dot-Bracket notation?

The ladder distance (LD) between residue i and residue j is the minimum number of base pairs (created due to secondary structure of the ssRNA) one has to cross when going from i to j. The Maximum ...
user avatar
4 votes
3 answers
430 views

Expression of a gene in different groups

I would like to check the expression of a gene in different groups like Disease vs Normal samples. I want to make a plot out of that to check whether it is significant or not. From this paper lncRNA ...
beginner's user avatar
  • 621
4 votes
1 answer
379 views

Plot to show the expression of genes between tumor and normal

I have RNA-seq raw counts data for 50 samples. 20 Normal and 30 tumor. After differential analysis I got 30 DEGs. I want to make a violin plot showing the expression of each gene. I transformed counts ...
beginner's user avatar
  • 621
4 votes
1 answer
77 views

RNA folding at specific temperature with ViennaRNA in python

I am trying to get dot-bracket-notations of single stranded RNA via the ViennaRNA python package (https://pypi.org/project/ViennaRNA/) at different temperatures. I have read in the docs (https://www....
gojih's user avatar
  • 41
3 votes
4 answers
716 views

Basic questions about GSEA

I'm currently learning about Gene Set Enrichment Analysis (GSEA) in the hopes of using it in my analysis of differentially expressed genes, and I just had a few questions about the program, ...
J0HN_TIT0R's user avatar
3 votes
2 answers
6k views

How to convert featureCounts to FPKM?

I have seen many posts regarding counts to RPKM and TPM. I haven't seen any post for counts to FPKM. I have RNA-Seq data which is paired-end reads. Extracted the counts using featureCounts for all ...
beginner's user avatar
  • 621
3 votes
2 answers
942 views

Detecting differentially expressed genes with foldchange >= 2 and FDR < 0.05

I'm using edgeR for differential analysis. Using glmTreat function I'm detecting differentially expressed genes between Tumor and Normal. I have set the arguments like below: ...
stack_learner's user avatar
3 votes
1 answer
324 views

How to select only RNA with Hetero atoms from pdb file with python?

I'm trying to separate RNA from protein in a complex protein/RNA PDB file and I want all RNA info with the hetero atoms in between the bases BUT without H20 etc. In short I want RNA part of pdb file ...
Raph's user avatar
  • 61
3 votes
3 answers
202 views

How to visualize genome track of gene in specific cell-lines?

I'm trying to make a plot showing genome tracks of specific genes in specific cell-lines of RNA-seq and Chip-seq data. It should look something like this I have recently seen this encode, but in ...
stack_learner's user avatar
3 votes
1 answer
2k views

Strange per sequence GC content results

I would be grateful if someone could take a quick look at these FASTQC results. This is rna-seq paired-end data. From the FASTQC manual, an unusual distribution seems to be suggestive of contamination ...
beginner's user avatar
  • 621
3 votes
1 answer
321 views

Why is bulk RNA sequencing reflecting AVERAGE expression but not TOTAL expression of all cells?

When I am reading papers that compares bulk RNA sequencing and single-cell RNA sequencing, we often see papers describe bulk RNA seq measures the average cell expression. For example, in this paper ...
benson23's user avatar
  • 151
2 votes
3 answers
243 views

Question about Co-expression analysis and finding targets for lncRNAs

I have a dataset with 88 tumor and 113 normal samples. Among the 50k genes after filtering there are a total of 28k genes (both mRNAs and lncRNAs) I wanted to do co-expression analysis between ...
stack_learner's user avatar
2 votes
2 answers
86 views

Which software is used to run Molecular Dynamics simulation of a small RNA hairpin?

I am currently working in a project that needs MD simulation of a small RNA hairpin. I am currently trying it on GROMACS but it seems to be protein based, hence am facing several issues. Is there any ...
user avatar
2 votes
3 answers
320 views

Selection of differential expressed genes

I'm working with RNA-seq data. I have 40 tumor samples and 5 Normal samples. Differential analysis with Deseq2 based on Fold change > 1.2 and alpha < 0.05 gave ...
beginner's user avatar
  • 621
2 votes
2 answers
38 views

Is there a computational tool or possibility to identify mRNA isoforms from the count matrix of a bulk RNA sequencing dataset?

I have the counts matrix of an RNA sequencing dataset of fibroblasts and I wish to identify isoforms of a particular gene of interest in it. Can anyone please hint me on a bioinformatics method to ...
Charles's user avatar
  • 516
2 votes
1 answer
30 views

Real time transcript profiles

Do any methods exist (or are in the process of development) for investigating transcript data without lysing the cells, i.e, destroying the sample?
libby's user avatar
  • 121
2 votes
1 answer
44 views

Getting VCF file that contain common SNPs from 6 VCF file using isec

I have 6 VCF files, where I would like to obtain the SNPs that are common (by position) in all the 6 files. I have tried this command ...
Mohamed Samir's user avatar
2 votes
1 answer
30 views

rRNA genes missing from metagenomic bins; is there a way to recruit them?

I have a couple of metagenomic bins which are okay in quality, but often missing rRNAs (16S, 23S...). I assume this is due to high population variability, combined with the high conservation of rRNAs, ...
Laura's user avatar
  • 787
2 votes
3 answers
75 views

beginner RNA-seq Replicate papers

have a good R and statistical analysis background (also with machine learning). in addition, I'm a fresh biotechnology grad. I would like to try to replicate some Rna-seq analysis with R papers (with ...
Layan 's user avatar
  • 31
2 votes
1 answer
64 views

Understanding ViennaRNA RNAdistance scoring table

I'm trying to compare the output of 2 different algorithms of RNA structure prediction (my implementation of Nussinov vs RNA-mfold algorithm) using the RNAdistance algorithm that is part of ViennaRNA ...
Eli Zatlawy's user avatar
2 votes
1 answer
57 views

I have used TER to break the long bonds of a chain in my PDB

I’m now not sure what I need to alter in my PDB to get it to work in leap. I know breaking the bonds turns the formerly connected residues into terminal residues. It keeps saying that 3 of my atoms no ...
user avatar
2 votes
1 answer
57 views

Protein misfolding

I am looking for literature on protein and RNA mis-folding. I am on the computational/biophysics side, and what interests me are the practical applications: diseases caused by misfolding, drug ...
Roger Vadim's user avatar
2 votes
2 answers
237 views

Tool for rna/lna melting temperature prediction

Is there any tool available to predict the melting temperature for rna/lna oligos with its perfect complement ? I know this amazing website from qiagen/geneglobe that estimates the tm value. But ...
Paul's user avatar
  • 297
2 votes
1 answer
62 views

How can I interpret gene expression data from Bioconductor packages?

I am currently looking at microarray data from a bioconductor microarray dataset. Specifically, I have data (a snippet) which looks like the following: ...
Larson's user avatar
  • 21
2 votes
1 answer
338 views

How to get rid of the temp files while using "featureCounts" for extracting readcounts from bam files?

I used this command for extracting read counts for almost 60 samples using featureCounts. ...
beginner's user avatar
  • 621
2 votes
1 answer
57 views

Is there a graphical/interactive 16S rRNA clustering method?

I've been doing phylogenetics with large (hundreds) 16S rRNA sequences lately. Usually I'm focusing on one order, and using a combination of trees and sequence similarity to assess stuff like 'is this ...
Laura's user avatar
  • 787
2 votes
1 answer
130 views

How does FDRtool work?

I have a question about using FDRtool. In the below code (on RNA seq data whose p values were acquired using Deseq2), the FDRtool was first used and thereafter p.adjust using the benjamini hochberg ...
Emma's user avatar
  • 31
2 votes
1 answer
43 views

Simple Affy Differential Expression

I think I have a fairly easy-to-solve problem. I have an expression table generated from an Affy set. It has one row of gene names (already symbols) and then 180 sample rows. 181 rows total. So I ...
julianstanley's user avatar
2 votes
0 answers
101 views

ExecutiveAlign error: atomic alignment failed (mismatched identifiers? )

So I have a molecule (a dsRNA) , that has two conformations. I want to align it in pymol with molecules but I get the following error: ExecutiveAlign error: atomic alignment failed (mismatched ...
majumderS's user avatar
2 votes
0 answers
77 views

Adapter trimming

I am trying to do adapter trimming, alignment and sorting for a range of large scale paired end fastq files. The code I am using is given below: ...
Aranyak Goswami's user avatar
2 votes
0 answers
45 views

Explaining the algorithm of RNA folding: what each symbol & value represents?

I need someone to explain to me from the nuts and bolts how algorithms/ maths is used to work out RNA folding. Explain it to me like I am an alien or child. I am looking at this paper - https://eprint....
confused_novice's user avatar
2 votes
0 answers
103 views

How is that ensemble Free energy is lower than the MFE in RNAfold from ViennaRNA

In the tutorial of RNAfold it states, RNAfold reads RNA sequences from stdin, calculates their minimum free energy (MFE) structure, prints the MFE structure in ...
aerijman's user avatar
  • 585
1 vote
2 answers
5k views

How to filter ribosomal RNA from scRNA-seq data

I want to filter out ribosomal RNA from scRNA-seq data (downloaded from here). Is there a list of known ribosomal RNA? The only solution I found is SortMeRNA, however it works with raw sequencing ...
gc5's user avatar
  • 1,773
1 vote
2 answers
627 views

Why the t-test for a specific gene shows different value compared to differential analysis?

I have RNA-Seq data for LUNG cancer. 370 tumor and 50 Normal. For differential analysis initially I did some filtering and kept approx. 19k genes for further analysis. I used edgeR. With a FC > or &...
beginner's user avatar
  • 621
1 vote
1 answer
54 views

Gene symbol list for all protein coding genes in mice

How can I get a csv or a list of the gene symbol names for all the protein coding genes in mice? I have RNA sequencing data and I'm not interested in the non-coding stuff. I'm worried it could mess ...
Angus Campbell's user avatar
1 vote
2 answers
96 views

Basic RNA Differential Expression in R

I have two matrices, one for individuals before treatment and one for the same individuals after treatment. Both matrices are raw read counts of RNA expression. ...
aurelius_37770's user avatar
1 vote
2 answers
177 views

Are codons in RNA layered? Are we misinterpreting RNA codons?

I am analyzing nucleotide base-pair patterns in RNA and DNA, and had a thought about RNA and DNA (Let me first state though, I am not a biologist; I am an algorithmatician, employing a sort of ...
Jonathan Charlton's user avatar
1 vote
3 answers
332 views

Deep learning RNA sequences

Currently I'm working on a project, which combines deep learning with RNA sequences. I'll try to predict pseudotorsion angles [1] from raw rna sequence. The ideas is to train a neural network with raw ...
Patrick's user avatar
  • 11