Question: ERCC count matrix smart-seq2
1
gravatar for VHahaut
10 months ago by
VHahaut1.1k
Belgium
VHahaut1.1k wrote:

Hi!

I received some data produced by single-cell sequencing (smart-seq2). After aligning onto the corresponding genome (+ERCC) I got a matrix of count which look quite sparsed:

ID  start   end strand  length  count1  count2  count3
ERCC-00004  1   523 +   523 1902    145 2328
ERCC-00009  1   984 +   984 318 29  428
ERCC-00012  1   994 +   994 0   0   0
ERCC-00013  1   808 +   808 0   0   0
ERCC-00014  1   1957    +   1957    0   0   0

With typically higher concentration ERCCs being detected but not lower ones.

I know for sure that those samples were not sequenced deep enough but I wanted to know if for ERCCs in single-cell sequencing I should expect to see every ERCC detected or is it common to see some drop-off ?

Thanks in advance!

rna-seq • 574 views
ADD COMMENTlink modified 10 months ago by Charles Plessy2.5k • written 10 months ago by VHahaut1.1k
4
gravatar for Charles Plessy
10 months ago by
Charles Plessy2.5k
Japan
Charles Plessy2.5k wrote:

The ERCC RNA spike-in mixes manufactured by Thermo Fisher cover six orders of magnitude of concentrations, therefore it is expected by design that a large number of them will not be detected. You can have a look at Svensson et al., 2007 for an example on how to analyse the spike data.

Side comment: while people have been for a long time aligning their reads to the sequence of the plasmid inserts used to produce the spikes, the NIST now distributes the sequence of the transcribed products, which are more accurate at the 5′ end. It may not make a big difference for a RNA-seq analysis, but I recommend to use these reference sequences (disclaimer: I contributed to prepare the file, see Patch ERCC spike sequences to get their real 5-prime ends. for details).

ADD COMMENTlink modified 10 months ago • written 10 months ago by Charles Plessy2.5k

Isn't ERCC-00012 one of the ERCC controls with amongst the highest concentrations, though? I would have expected that to show some counts.

ADD REPLYlink modified 10 months ago • written 10 months ago by Kevin Blighe28k
1

Thanks for your comments! According to thermoSc. website ERCC-00012 is one of the lowest:

ERCC-00012  C   0.11444092 attomoles/ul
ADD REPLYlink written 10 months ago by VHahaut1.1k

Then it makes sense what Charles wrote. Great, thank for confirming.

ADD REPLYlink written 10 months ago by Kevin Blighe28k
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: 1628 users visited in the last hour