Questions tagged [transcriptome]

The set of all transcripts of a biological entity, e.g. the transcriptome of a single cell, of a tissue, of an organ, of an individual, of a species. Can also be used to refer to a subset of the transcriptome, e.g. the polyadenylated transcriptome, given that often experiments ignore the 99% of transcripts that are ribosomal RNAs.

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

How to visualise /add Arabidopsis T-DNA Lines As Tracks On IGV

We have done some RNA-Seq on Arabidopsis T-DNA lines and mapped the reads on to Arabidopsis genome. Now we wanted to confirm the T-DNA deletion of those lines by loading them on IGV. As a control we ...
0
votes
0answers
13 views

What Galaxy tools can annotate a de novo assembly RNA seq plant transcripts?

What Galaxy tools can annotate a de novo assembly RNA seq plant transcripts? How to do it using Python API?
2
votes
1answer
241 views

How to find novel transcripts using GFFcompare?

I am trying to find novel transcripts from an RNA-seq database. Based on the advice I got, it seemed that using Stringtie for transcript assembly is a good way to go, and it supports novel transcript ...
0
votes
2answers
65 views

Proteomics data Vs Transcriptomics data?

I want to use either of Proteomics or Transcriptomics data for integrating it into my kinetic model. Before proceeding, I want to know what are the advantages of using either of them so that I could ...
0
votes
1answer
145 views

What are the right parameters to trim a small RNA transcriptome with trimmomatic?

I'm having some problems with finding the right parameters to trim my small RNA Illumina reads (51 nt long) with Trimmomatic. Before trimming, one of the samples (21M reads) looks like this: So for ...
0
votes
0answers
8 views

How to integrate the transcriptomics data with kinetic metabolic models?

I have created a convenience kinetics model. Now, I want to integrate the transcriptomic data with my convenience kinetics* model for altering/weighing the kinetic parameter values. I read some ...
2
votes
1answer
258 views

How to highlight specific cells in Seurat 2.4

I used Seurat 2.4 on our scRNA dataset to obtain the following tSNE plot.I was able to successfully extract cell IDs from the different clusters, and generate gene expression profiles. The analysis, ...
1
vote
1answer
69 views

Transcription factor targets for Human and Mouse genes

I want to get transcription factor targets for my list of genes (For Human and Mouse) I have found ENCODE Transcription Factor Targets database (The link is given below). But I am looking for some ...
3
votes
1answer
74 views

Merging transcriptomes coming from different experiments

I'm planning to build a transcriptome by pooling all existing transcriptomes in SRA for a non-model species (which has no reference genome) to study differentially expressed genes and the like. It ...
0
votes
1answer
108 views

Assembly by stringtie

I run this cmd ./stringtie G1_sorted.bam -B -o G1.gtf -G Triticum_aestivum.IWGSC.42.gtf -p 4 -C G1.refs.gtf -A G1.abund.tab Error is: ...
0
votes
1answer
55 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 ...
2
votes
1answer
82 views

How I know which gene is a good predictor in this neural network or not?

I have 25 highly differentially expressed genes among and patients to chemotherapy. I have made a neural network of these genes. Accuracy of model is 0.73 but I don't know from 25 genes how I could ...
2
votes
2answers
78 views

scRNA-seq differential transcript usage

Many of the modern gene-quantification tools (Salmon/Kallisto) output transcript-level (as opposed to gene-level) data. All of the scRNA-seq analysis I have seen just uses the gene-level values. I ...
2
votes
1answer
27 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?
3
votes
0answers
34 views

Simulating 3' end tag-based scRNA-seq reads

Are there any tools that simulate 3' end tag-based single-cell RNA-seq reads? For example, Drop-seq, 10X Chromium, CEL-seq, etc.. I see there are tools that simulate scRNA-seq gene count data (e.g. ...
1
vote
0answers
16 views

Finding gene name from human genome using SP1 transcrition factor binding site from Postion Weight Matrix

I have a list binding site motif of SP1 transcription factor which collects from PWMScan tools. Actually this tool scan whole human genome and give the entire binding motif list. From that list of ...
4
votes
1answer
103 views

At what stage of a transcriptome assembly is it better to perform read contaminant filter?

I'm trying to assemble a bivalve transcriptome. Since bivalves are filter feeders, their transcriptomes tend to be highly contaminated by bacteria, algae and whatnot. Since I pooled several ...
3
votes
0answers
56 views

RNA-Seq type and and optimal fusion detection

There are several popular types of RNA-Seq library prep which are frequently used: total RNA (with and without ribosomal depletion), mRNA-Seq/poly-A, and targeted mRNA-Seq/RNA exome. What would the ...
2
votes
0answers
165 views

Reproducing GTEx transcriptome analysis

I am willing to reproduce part of the analysis from "The human transcriptome across tissues and individuals" (Melé et all, 2015). I downloaded GTEx v6 FPKM data in txt format from GTEx portal. I want ...
4
votes
0answers
44 views

What is the current state-of-the-art in assembling hybrid transcriptomes?

We are considering attempting de novo assembly of a species transcriptomes (i.e. without a reference genome) using the combined NGS outputs of Iso-seq and Illumina. One example I saw (Li et al 2017),...
2
votes
0answers
37 views

Merge / Reconciliate several de novo transcriptome assemblies with different kmers

I am trying to build a De Novo transcriptome reference assembly for an eukaryotic organism. I've created several assemblies with rnaSpades using different kmer sizes (19 to 69 with step 10). Now I ...
2
votes
2answers
774 views

gene-level versus transcript-level analysis

Traditionally, RNA-seq data was quantified on gene level. Newer methods quantify on transcript/isoform level. For example, Kallisto only outputs transcript-level abundances. From the DESeq2 vignette: ...
1
vote
0answers
95 views

Variability in genes across platforms

I have a matrix of gene expression for single cell (cells in column and genes in row) in 9 time points with about to 200 cells for each time point. I have also a matrix of gene expression in bulk RNA-...
3
votes
1answer
28 views

Batch detection of CRISP proteins in fasta file

Probably a naive question. I am inexperienced. I am interested in identifying potential CRISP (Cysteine-rich secretory proteins) in a certain tissue transcriptome (ca. 20k sequences in fasta). I ...
5
votes
2answers
172 views

Reduce number of transcripts in a highly variable de novo transcriptome assembly

I have a de novo assembly using both multiple SRA and locally sequenced transcriptomes. I started with 270M PE reads from 9 tissues. Here are the assembly stats generated with ...
2
votes
1answer
50 views

Correct way to assemble reference transcriptome - what's --samples_file option in Trinity?

I have 101 samples from 9 tissues of a nonmodel species, coming either from SRA or sequenced by my lab. I want to generate a reference transcriptome with Trinity for further analysis. In order to do ...
7
votes
2answers
139 views

What is “Possible polymerase run-on fragment”?

I used gffcompare recently to compare a newly assembled transcriptome and a published transcriptome. The class code "p" transcripts are "possible polymerase run on". What does that mean? I tried very ...
2
votes
1answer
372 views

What process and input data is required for a cellranger reference transcriptome?

I'm analysing single-cell RNA-Seq data using the 10X Genomics cellranger platform. While they provide reference data for Mouse and Humans, other species require a ...
3
votes
1answer
179 views

Hierarchical models with limma?

I have a dataset with (microarray) gene expression data that was sampled from the same individuals at multiple timepoints. Our exposure is a continuous variable, and because this was an observational ...
-1
votes
1answer
190 views

Command not found error on HISAT2

I have been attempting to map reads to a reference genome using HISAT2 using the Pertea, et al 2016 Nature Protocols paper. I am able to successfully generate an index (using the hisat2-build command),...
2
votes
0answers
77 views

How do I compare transcriptomes of the same species across different developmental stages?

I want to identify the most highly expressed genes across several developmental stages of a given insect, in hopes of identifying a protein that can be used to design an insecticide. I have a few ...
5
votes
2answers
444 views

Difference between de novo transcriptome assembly methods

I have been looking around (including read the original papers) to understand what is essentially the difference between StringTie in non-reference based mode (de novo) and Trinity de novo assembly. I ...
4
votes
3answers
491 views

How to compare transcriptomic profiles of two cell types (single cell RNA-seq)?

I found this interesting Single RNA-seq data set in GEO, but I am not sure how to analyze it appropriately. They have deposited transcriptomic profiles of human and mouse pancreatic islets (...
5
votes
2answers
132 views

RapMap: reference transcriptome for simulated reads

I am very new to bioinformatics and trying to repeat the benchmark in the RapMap paper with an experimental tool working in a similar but different fashion. In the paper (taken from their github) ...
7
votes
1answer
483 views

Coverage calculation: long reads (RNA-seq)

Say your aim is to calculate the coverage of an RNA-seq experiment generated with long-read sequencing (so, uneven read length). Up to now, I relied on the Lander/Waterman equation: $$C = L*N / G$$...
7
votes
0answers
74 views

Run cuffcompare in strand-agnostic mode

Is there a way to run Cufflinks' cuffcompare in a strand-agnostic mode? I would like to do this because I have some RNA-seq datasets derived from an unstranded run, that should be compared to a ...
5
votes
1answer
108 views

Documentation for Transfrag class codes (cuffcompare)

Is there any extensive documentation or description for each class of Transfrag class codes as reported by the cuffcompare tool in the Cufflinks package? Official doc might not be the best. E.g. ...
5
votes
3answers
292 views

tools to reconcile experimental transcripts with reference annotation

Looking for tools to reconcile alignment file of experimental transcripts mapped to genome (SAM/BAM) with the reference transcriptome annotation (GTF) from Ensembl (organism: D. melanogaster). The ...
5
votes
1answer
166 views

PASA pipeline: compare experimental transcripts to the reference annotation

I would like to ask if anyone has experience in running a subset of the PASA pipeline, in particular for the reconciliation of some experimental 'transcripts' with the reference annotation. In more ...
4
votes
3answers
324 views

Analysis of differential transcript usage (DTU)

Recent breakthroughs in bioinformatics tools for quantification (e.g. Cufflinks/Kallisto/Salmon etc.) and tools which can identify differential transcript usage (DTU) (e.g. DRIMSeq, Cufflinks etc.) ...
15
votes
4answers
4k views

How to compute RPKM in R?

I have the following data of fragment counts for each gene in 16 samples: ...
10
votes
1answer
201 views

Linear models of complex diseases

A popular framework to analyze differences between groups, either experiments or diseases, in transcriptomics is using linear models (limma is a popular choice). For instance we have a disease D ...
13
votes
2answers
2k views

Normalization methods with RNA-Seq ERCC spike in?

ERCC spike-in is a set of synthetic controls developed for RNA-Seq. I'm interested in using it to normalize my RNA-Seq samples. In particular, I'd like to use the spike-ins to remove technical bias ...
14
votes
2answers
907 views

Alignment based vs reference-free (transcriptome analysis)?

I want to focus on transcriptome analysis. We know it's possible to analyze RNA-Seq experiment based on alignment or k-mers. Possible alignment workflow: Align sequence reads with TopHat2 Quantify ...
0
votes
2answers
7k views

What is the difference between a transcriptome and a genome? [closed]

I have a computer engineering background, not biology. I started working on a bioinformatics project recently, which involves de-novo assembly. I came to know the terms ...