Reference counts for RNA-seq
I have TPM counts from 36 participants, but they are all diseased.
- aligned to GRCh37.75
- ensembl genes
Are there any good sources of "healthy" control TPM counts?
expression
rna
rna-seq
tpm
transcripts
• 2,131 views
•
link
updated
by
ATpoint
9.1K
• |
written
by
LayneSadler
10
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Differentially Expressed Genes between two conditions (scRNA, single GEO dataset with multiple samp…
written by prietto 1Hi guys, My aim is to identify Differentially expressed genes from GEO183837. This dataset contains Single-cell transcriptome profiling of the human endometrium from patients with …
-
Calculation of tRNA abundances from different cell types
written by nkabo 8Dear all, My research group had completed an analysis on tRNA levels before I joined to study. The aim was to compare tRNA abundances of …
-
Use TPM to compare between the samples
written by Kunal 0After going through multiple community posts and papers, I understand that TPM cannot be used for between sample analysis. But currently I am stuck with …
-
How to find trusted healthy control data for childhood disorders
written by nkabo 8Dear all, I have been working on a rare condition occurring in children. Since it is a rare disease, I have 3 samples from diseased …
-
How to evaluate gene expression
written by Lorenzo 1Hi Biostars community, Let me explain my situation: I have to determine the expression of different genes in different tissues. Each tissue has a biological …
-
How do I use t-test for genes in R?
written by bkalem4 0Hi everyone, I want to identify differentially expressed genes from a dataset I got from NCBI GEO. The data is composed of 14 .CEL samples. …
-
plotting a expression data in figure
written by rob.costa1234 32I have expression matrix of 20 genes in 5 replicates in healthy vs disease animals treated with a drug. I want to show that the …
-
Where to find open access microbiome data for case-control comparative study?
written by stats134711 0I would like to do a comparative study of some methods that require human microbiome data from two groups, preferably healthy controls and diseased cases. …
-
How to find transposable element genes in ensembl annotation
written by james.lloyd 10I have some differential genes from an RNA-seq experiment and I wanted to look at changes of transposable element genes (transposons). This data is human …
-
Statistical technique for mutation to gene expression link
written by RBee 0I have following statistical data from a set of samples (healthy and diseased tissue - for approximately 100 individuals). For each sample: - gene expression …
The only good source would have been to include matched tissue in your study. RNA-seq is a relative assay and batch effects between unrelated studies make it close to impossible to meaningfully compare them unless prepared in the same batch, same kit, same procedures, same everything. Are you aware of that? What information do you exactly seek?
That's right, you cannot do DE without samples from the same experiment. At most, you can compare the top expressed genes against the GTEX data to get an idea of whether your top expressed genes are also normally expressed in the same tissue. There are a few methods to do meta-analysis for RNASeq data, for example, this one metaseq - could the user apply one of these?
That having said, it is completely unstandardized how TPM is calculated. Correctly one would use the length information based on the factual length of the transcripts being expressed, like e.g. the salmon-tximport pipeline provides it. But some people use either the entire length of the annotated genes, the average of all transcripts, the union if exons, etc to calculate it. Hence, on top of the experimental batch effect different sources of TPM might have in silico batch effects making even top-wise comparisons difficult. Sure, if something is zero in one and skyrocketing in the other it could be true, but everything else is to be considered with much care (or not at all).
Thank you. I am constructing a diagnostic algorithm that uses gene TPMs to predict the presence of the disease. Then I want to permute the TPMs to figure out the most important genes. So I need TPMs from both cases and controls.
Without knowing details it sounds you are using "ground truth" that is confounded, so in all likelihood performance of that algorithm will suffer. Arguments have been made above, on you to follow it or not, but I personally would get a collaborator and generate case and controls in a matched study to actually have a solid ground truth.