Question: Differential expression without technical replicates
3
gravatar for jakalssj3
3.8 years ago by
jakalssj380
Muenster
jakalssj380 wrote:

Hello everyone! I've got a question about a sense of conducting a differential expression analysis without having technical replicates. What I have now, is RNA-seq data from four developmental stages of a parasite (1 biological replicate per sample, no technical replicates). The genome size is ~40 Mb and for each sample there is between 100 to 170 million paired-end reads, so sequencing depth is really big. I have used both DESeq and Cufflinks to estimate expression on a gene level. In the case of having no replicates, should I just focus on fold change between the samples? I'm not sure if the calculated p-values in both DESeq and Cufflinks are meaningful in case of no replicates.. Thanks for any suggestions!

ADD COMMENTlink modified 3.8 years ago by Gary450 • written 3.8 years ago by jakalssj380
4
gravatar for Devon Ryan
3.8 years ago by
Devon Ryan88k
Freiburg, Germany
Devon Ryan88k wrote:

Technical replicates are largely useless in RNAseq, biological replicates, in contrast, are vital. If you have only a single sample per condition (it's somewhat unclear if your "1 biological replicate" means 1 or 2 samples per condition/timepoint) then the yes, you're largely stuck looking at fold-changes. You might try GFold, though I don't know if it allows time-course designs. Regardless of the tool, be highly suspect of the results, they'll only be vaguely useful.

ADD COMMENTlink written 3.8 years ago by Devon Ryan88k

Sorry, what I have is all together 4 RNA-seq datasets - 1 dataset per developmental stage of the parasite. Additionally, couple of genes were tested with RT-PCR by my colleague, so I can always use it as sort of benchmark to see if the observed fold changes from RNA-seq are being reflected in RT-PCR. Also, there is another very closely related species, for which gene expression was estimated in two (ouf of four) developmental stages, but with several biological and technical replicates this time. I'm just assuming that these results could also serve for checking if the observed expression patterns from my analyses are correct.

ADD REPLYlink written 3.8 years ago by jakalssj380

That seems reasonable.

ADD REPLYlink written 3.8 years ago by Devon Ryan88k
1
gravatar for Gary
3.8 years ago by
Gary450
Taiwan/Taichung/China Medical University Hospital
Gary450 wrote:

For a time-series RNA-Seq data without replicates, STEM (Short Time-series Expression Miner) is a pretty well option. STEM cannot identify significantly expressed genes. It groups genes based on their expression patterns along the time-course, and tells us whether a group of genes is significantly enriched for particular GO terms. Of course, this software can provide more promising results with biological replicates.

STEM was developed by Jason Ernst for microarray analysis in the beginning. This tool can analyze our RNA-Seq data very well, too. Now, Jason is an assistant professor in UCLA, and friendly to respond questions about STEM based on my limited experience.

 

Short Time-series Expression Miner (STEM)

http://www.cs.cmu.edu/~jernst/stem/

Jason Ernst Lab

http://www.biolchem.ucla.edu/labs/ernst/

ADD COMMENTlink written 3.8 years ago by Gary450
0
gravatar for SmallChess
3.8 years ago by
SmallChess480
Australia
SmallChess480 wrote:

Without any replicate, you can't really estimate the dispersion factor. However, you mention in your comments that your colleague has done some testing before, maybe you can estimate a dispersion from it? Note that the main point to have replicates is to estimate dispersion, if you have some kind of data to estimate it, you might get away without any replicate. Please take a look at the EdgeR user manual.

If you don't have the information, I'm afraid you should just look at log fold-change.

ADD COMMENTlink written 3.8 years ago by SmallChess480
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1318 users visited in the last hour