Question: How to align .fastq files to an Addgene vector sequence?
0
gravatar for Kristin Muench
15 months ago by
United States
Kristin Muench410 wrote:

Hello,

I would like to find out if any reads in my .fastq file were transcribed from a vector sequence that is not in the human reference genome.

The sequence is available here: https://www.addgene.org/27080/sequences/#depositor-full

I thought I might try aligning the .fastq files to the sequence using TopHat, as if the sequence were a human reference genome, and seeing if any alignments pop up.

However, I'm not sure how to go about doing this.

Should I make the sequence above into a .gtf file somehow?

How do I make the corresponding annotation (.gff) file?

Is there an easier way to go about doing this? E.g., isolating every sequence in the .fastq file and using grep to search for it in the vector sequence?

Thank you!

EDIT: I am aligning Illumina RNA-Seq .fastq files. Also, I'd appreciate any resources folks have regarding how to modify or add chromosomes to a reference file!

rna-seq alignment • 562 views
ADD COMMENTlink modified 15 months ago • written 15 months ago by Kristin Muench410
1

My standard answer when someone mentions Tophat:
You should know that the old 'Tuxedo' pipeline of Tophat(2) and Cufflinks is no longer the "advisable" tool for RNA-seq analysis. The software is deprecated/ in low maintenance and should be replaced by HISAT2, StringTie and ballgown. See this paper: Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. (If you can't get access to that publication, let me know and I'll -cough- help you.) There are also other alternatives, including alignment with STAR and bbmap, or pseudo-alignment using salmon.

ADD REPLYlink written 15 months ago by WouterDeCoster39k

Eek, thank you for the heads-up! I'll take a look at those other methods.

ADD REPLYlink written 15 months ago by Kristin Muench410
2
gravatar for WouterDeCoster
15 months ago by
Belgium
WouterDeCoster39k wrote:

I believe the most correct would be to add an additional chromosome, your vector, to the human genome and index that for alignment using e.g. STAR or HISAT2 (assuming you have Illumina RNA-seq data, which you did not specify).

Aligning only to the vector (without the human genome) or vice versa might lead to false positive/negative alignments.

ADD COMMENTlink written 15 months ago by WouterDeCoster39k

Great point! And yes, I am using Illumina RNA-Seq data.

ADD REPLYlink written 15 months ago by Kristin Muench410
1
gravatar for genomax
15 months ago by
genomax67k
United States
genomax67k wrote:

You could use bbsplit.sh from BBMap suite to bin the reads to quickly find out which ones are from human genome and which are from vector. You can decide how to handle those reads that multi-map (both within and across the genomes).

ADD COMMENTlink written 15 months ago by genomax67k

Oh, interesting! Thank you!

ADD REPLYlink written 15 months ago by Kristin Muench410
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: 1563 users visited in the last hour