Does DESeq2 need "scale_counts" for GTEx and TCGA data from recount2 (TCGAquery_recount2)?
1
0
Entering edit mode
@xiaofeiwang18266-13498
Last seen 15 months ago
Singapore

Dear there,

From the TCGAbiolinks paper (https://journals.plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1006701), I see it need scaling Recount2 data because it was sequenced using Rail-RNA. Here is the code they used as below for downstream differential expression analysis (DEA) using "TCGAanalyze_DEA". They also used "TCGAanalyze_Normalization" and "TCGAanalyze_Filtering" to do the normalization and filtering. Since DESeq2 takes the raw read counts as input, I am wondering do I need any of these 3 (scaling, normalization, and filtering) if I used "TCGAquery_recount2" to download GTEx and TCGA data? Thanks a lot!

Xiaofei

#####Preparing/scaling Recount2 data because it was sequenced using Rail-RNA

eset.gtex<-assays(scale_counts(ucs.recount.gtex$GTEX_uterus, round = TRUE))$counts
eset.tcga<-assays(scale_counts(ucs.recount.tcga$TCGA_uterus, round = TRUE))$counts
TCGAbiolinks RECOUNT2 DESeq2 • 1.6k views
ADD COMMENT
1
Entering edit mode
@mikelove
Last seen 1 day ago
United States

recount has a function scale_counts to create a count matrix appropriate for running DESeq2, edgeR, limma-voom or other count based packages. From their vignette, in the code just before running DESeq2:

## Scale counts by taking into account the total coverage per sample
rse <- scale_counts(rse_gene)

If that function has been run, then you are set.

ADD COMMENT

Login before adding your answer.

Traffic: 739 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6