Questions tagged [rna]

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

Filter by
Sorted by
Tagged with
0
votes
0answers
19 views

What's the difference between Seurat Merge and Integration [closed]

Note: this question has also been asked on Biostars I have been given some scRNAseq 10X data, including 3 groups: 4 normal, 8 disease, and 8 treatment samples. My goal is to find out how many cell ...
1
vote
1answer
9 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 ...
2
votes
0answers
26 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 ...
0
votes
1answer
51 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 ...
2
votes
1answer
93 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 ...
0
votes
1answer
40 views

How to count RNA sequencereads using custom made python scripts?

I am trying to do RNA seq analysis and my goal is to filter gene counts less than 5 using custom made python scripts. The code chunk goes as follows ...
0
votes
1answer
147 views

How to handle control samples in CLIP-Seq

I have a CLIP-Seq dataset I'm processing, which includes control samples and no inputs. This is the second CLIP analysis I've performed to help out users of our genome core facility and the first one ...
3
votes
1answer
55 views

Practical use of RNA structure

I have spent some time studying RNA structures. While there is a range of interesting phenomena and functions that certainly deserve understanding from the scientific point of view, I have never ...
1
vote
2answers
88 views

Bacterial DNA at the tail of transcriptome reads. What does that mean?

I am assembling a transcriptome obtained from the Internet. The transcriptome was extracted from a human cancer tissue that had been previously grafted into a mouse. I have detected that many ...
7
votes
1answer
137 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 ...
0
votes
1answer
24 views

Is is possible to predict ncRNAs from sequence and homology alone?

I'm working with a set of homologous genes (let's call it gene A) from several bacterial species. I know (from previously published research) that in gene B (a close paralogue of gene A), there is a ...
2
votes
1answer
28 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 ...
1
vote
2answers
158 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 ...
0
votes
1answer
22 views

smallest length of reads in 3'quantseq?

I am studying about RNA seq, and especially about 3'Quantseq(tagseq, 3primeseq). I wonder if there is a cutoff for the reads length. By this I mean that given that 3'quantseq targets the end of the ...
0
votes
0answers
10 views

RNA loop prediction: how to interpret differences between RNAfold and Mfold?

I'm trying to find out the structure of a conserved ITS loop. As my initial prediction (using only bases homologous to those present in the loop in a closely-related species) was fairly low-scoring, I ...
1
vote
0answers
39 views

Can someone explain how this algorithm describes RNA folding to me in terms of 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 a child, an alien-child. I am looking at this ...
0
votes
0answers
18 views

Detecting Transcript Isoforms with HISAT2/Stringtie RNAseq Output

I have the following Stringtie RNA-seq output files for several samples in the image below. For additional context, I have utilized as part of my RNA-seq workflow. HISAT2, samtools, and Stringtie ...
0
votes
1answer
75 views

about the cpm normalization after using normalization factor

Is it okay to use CPM normalization (with/without log transform) after using TMM normalization? Why do we need both? ...
0
votes
1answer
51 views

RNA seq .counts.txt to bigwig conversion

I've downloaded neuronal RNA seq data from GEO. The files are in .counts.txt format. Could I convert them to bigwig? If so, how?
3
votes
2answers
297 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: ...
0
votes
0answers
84 views

about the scaling after normalization and transformation in RNA-seq data

When we use count data in RNA-seq analysis, we usually use normalization and sometimes vst, rlog transformation (DESeq2)or log2(CMP+4) transformation (edgeR) to perform K-means clustering. Can we use ...
0
votes
1answer
48 views

about the normalization of RNAseq data for calculating distance for unsupervised learning

I have been working on clustering using RNAseq data. To compute distance, what kinds of normalization is optimal? Can we use normalization using relative log normalization? I understand this should ...
6
votes
2answers
2k 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 ...
0
votes
1answer
41 views

How to understand and analyse RNA-seq data (for a beginner)?

I am trying to understand expression of a certain protein across Pseudomonas species. I downloaded an SRA file from NCBI and converted it to a fastaq file. I am not able to understand how to interpret ...
2
votes
1answer
42 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 ...
2
votes
0answers
68 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 ...
2
votes
2answers
35 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 ...
0
votes
0answers
24 views

What are the meanings of these transcript ids?

I have three types of transcripts for Rosa chinensis "Old Blush" homozygous genome v2.0 from GDR genome browser. ...
0
votes
1answer
696 views

WGCNA: Problem with selecting soft threshold

I have 25 tumor samples with counts data. Initially, I filtered out low expressed genes and then converted counts to varianceStabilizingTransformation using ...
1
vote
1answer
64 views

Writing a python module to run Infernal on a list of fasta files

I am very new to using python3 and am hoping to use it to reduce the time I need to spend running Infernal to look for a variety of RNAs in a large collection of bacterial genomes. I am hoping to ...
1
vote
1answer
34 views

how to calculate coverage of each base of mapped reads?

I have a RNA aligned file without any header info. I need to calculate the coverage of each mapped base for per chromosome. Is it possible to do it with perl ? ...
0
votes
0answers
15 views

Arm-level SCNA does not change mRNA expression level?

It is almost common sense that local amplification of gene increase the its mRNA expression . I recently heard that large scale SCNA (such as arm-level amplification) does not increase the mRNA ...
5
votes
2answers
665 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? ...
3
votes
1answer
125 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 ...
1
vote
3answers
272 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 ...
0
votes
1answer
70 views

I need some tips and suggestions for further analysis of NGS expression data (log2cpm)

I am a PhD student who inherited some log2cpm data of expression data from bulk kidney tissue from a UUO(unilateral urethral obstruction) experiment that tests a new drug. The sample material consists ...
3
votes
2answers
4k 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 ...
2
votes
1answer
60 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: ...
1
vote
2answers
60 views

Retrieve RNA sequencing data for human p53 colon cancer cell lines

I want to download the next generation RNA sequencing data in different cell lines. I want the data (alongwith normal) for p53 wild type, p53 knock-out and p53 null cell lines (human colon cancer). ...
2
votes
3answers
220 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 ...
0
votes
3answers
76 views

Building RNAz on MacOS

Having trouble building RNAz 2.1 on MacOS. Following instructions at https://www.tbi.univie.ac.at/software/RNAz/ ...
2
votes
2answers
51 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 ...
11
votes
2answers
470 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....
2
votes
1answer
29 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?
1
vote
2answers
286 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 &...
1
vote
0answers
85 views

How to correctly parallelise RSeQC scripts with GNU parallel?

I have a .bam, as ouput from STAR aligner, from which I need to extract some info using RSeQC while using all the computational resources available to increase ...
4
votes
1answer
256 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 ...
1
vote
1answer
2k views

How to calculate logCPM across all samples?

Using edgeR for differential analysis between Tumor and Normal gave me differential expressed genes with logFC, logCPM, PValue and FDR. From the details of glmTreat function I see that logCPM is ...
2
votes
3answers
229 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 ...
1
vote
0answers
105 views

How to create a SAF file for virus sequence?

I am trying to create a SAF file for the virus sequence. I can download .gff3 and .gb format for the virus gene sequence but they are not supported by featureCounts. I need to create my own SAF. I ...