I am trying to understand how data that was uploaded to SRA (https://www.ncbi.nlm.nih.gov/sra?LinkName=biosample_sra&from_uid=4510743) can be analyzed with the assumption that the FASTQ file should contain replicate information (3 replicates within the sample). I downloaded the file and ran it with Nextflow, and it was mapped to the human genome and determined to contain unstranded reads of 50bp, both directions.

Which steps can I take to infer if the FASTQ file indeed is not demultiplexed?

Many thanks, Geo

The header of the FASTQ looks like this:

@SRR3467208.1 1/1
@SRR3467208.2 2/1
@SRR3467208.3 3/1
@SRR3467208.4 4/1

I would suggest looking at sample metadata:


Name: TruSeq RNA

Instrument: Illumina HiSeq 2000

Strategy: RNA-Seq


Selection: PolyA

Layout: PAIRED

You can see that it's paired. You can look at the spot descriptor at the page you link to see more information about the design.

There is no multiplexing information available there. I would instead look at the overall project and its metadata to find the different replicates. It looks like there are 2 samples x 3 replicates each, deposited independently. So it is likely to already be demultiplexed, and you can just download them individually.

"Unstranded" isn't really a meaningful term for this, if I understand correctly. Forward and reverse read pairs don't correspond to forward and reverse on the chromosome, but to the way that sequencing-by-synthesis happens on the machine. So the mapping orientation is not really relevant.

I infer that you are working from a single FASTQ file. If you are working from a single FASTQ file for paired-end out of SRA, then the FASTQ file is not yet resolved into F+R reads. I would suggest rerunning the FASTQ conversion with --split-3 or similar option.

  • $\begingroup$ Thanks! The term 'unstranded' was determined by RSeQC infer_experiment.py. The reads were 47% forward, and 47% reversely mapped. So the resulting counts from the single FASTQ file worked even w/o the FASTQ conversion. The experiments from the entire project are all 6 independent (6 x 1 sample), not 2 x 3 samples. I worry that the reported replicates are actually missing in the uploads and was hoping to confirm this from the FASTQ files somehow. $\endgroup$ – Geo Vogler Mar 24 at 21:02
  • $\begingroup$ @GeoVogler That is helpful background for trying to answer the question. Where do you see the 6 x 1 info? I was not able to track down the actual study that originated these samples, linking that would be helpful. All I am seeing is the 3 samples x 2 treatments on SRA. $\endgroup$ – Maximilian Press Mar 25 at 3:17
  • $\begingroup$ The three samples are from three individuals (mother/father/proband), not a single one (see Biosamples). The layout of SRA with Project/BioSamples/Reads could really be improved. If there is no indication of replicates I will need to contact the lab directly, I just thought maybe there is an easy way to gauge the presence of Index sequences inside a FASTQ file. $\endgroup$ – Geo Vogler Mar 25 at 3:27
  • $\begingroup$ @GeoVogler the index reads were probably split out in initial processing, but they are often preserved in FASTQ headers by e.g. bcl2fastq. SRA unfortunately has a tendency to strip out useful information from the FASTQ header, which does sometimes include the index sequence for each read. The authors should have the original files in theory, yes. $\endgroup$ – Maximilian Press Mar 25 at 13:21

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.