Spliced aligner for 454 reads?

sulicon · 05-18-2011, 01:43 PM

Hi all,
I am trying to assemble 454 reads by aligning them to reference genome.

I have tried the official gsMapper before, which just gave me exon sequences, rather than transcripts sequences. Now, I am using BWA-SW to align the reads, and Cufflinks/Scripture to reconstruct the transcripts. Still, I can only get the individual exon sequences, instead of expected transcript sequences.

It seems the problem here is that BWA-SW is not a spliced aligner so that the splicing junction reads will be lost during the alignment. A spliced aligner, Tophat, is used in the tutorial of both Cufflinks and Scripture. However, this aligner is based on Bowtie, an aligner designed for short reads only.

Could anyone give me some suggestion of the spliced aligner suitable for 454 reads? I think BLAT would be an option, but still want to test some other methods developed recently.

Thanks,
Shuli

lh3 · 05-18-2011, 02:41 PM

For a read bridging a splice junction, bwasw should give two or more hits unless one of them is too short. Perhaps Cufflinks is expecting some tophat/bowtie specific information to group local hits to a transcript. I do not know.

Nonetheless, I agree bwa-sw would not work well because it is a local aligner. For RNA-seq/ESTs, a dedicated splicing-aware glocal aligner is more appropriate. In addition to blat, you may also try gmap.

Jose Blanca · 05-22-2011, 11:48 PM

gmap can align 454 ESTs against a genome taking into account the introns.

jochensupper · 05-23-2011, 12:50 AM

Hi,

our mapper (Genomatix) has a local spliced alignment mode that allows to align complete transcripts to the genome. Attached is a screenshot of assembled (velvet) RNA-Seq reads mapped to the reference genome.

Depending on the organism and your objective you could also consider mapping your reads against a transcriptome library (with no worries about splicing your reads). Then, however, you wouldn't be able to discover novel transcript variants.

sphil · 05-23-2011, 01:29 AM

Blat is also able to handle 454-splice reads but be aware of long runtime....

sulicon · 05-24-2011, 02:14 PM

Thanks Li Heng and Jose. I have tried GMAP and checked the output. It looks good. According the a document I have found, GMAP outperforms BLAT for gene structure identification in both speed and accuracy. And the latest version of GMAP can generate output in SAM format directly, which facilitates the subsequent analysis.

Finally, I got some transcripts constructed by running Cufflinks on GMAP output. However, it seems that Cufflinks has modified the original alignment and generated some artificial exons/splicing junctions I've never seen in the GMAP output...

brdido · 11-01-2011, 09:39 AM

Quote:

Originally Posted by sulicon

However, it seems that Cufflinks has modified the original alignment and generated some artificial exons/splicing junctions I've never seen in the GMAP output...

Maybe this is why (from http://cufflinks.cbcb.umd.edu/manual.html):

-g/--GTF-guide <reference_annotation.(gtf/gff)> : Tells Cufflinks to use the supplied reference annotation (GFF) to guide RABT assembly. Reference transcripts will be tiled with faux-reads to provide additional information in assembly. Output will include all reference transcripts as well as any novel genes and isoforms that are assembled.

I'll try to use gmap+cufflinks for 454 data! Thanks.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
BLAST for spliced isoforms?	tboothby	Bioinformatics	8	09-07-2011 09:20 AM
Spliced alignment with BWA	telos	SOLiD	7	10-06-2010 07:32 AM
SOLiD WTP alignment file: representation of spliced reads	Simon Anders	Bioinformatics	0	08-19-2010 10:29 AM
Best tool to map 454 reads onto sanger reads?	dan	Bioinformatics	3	07-27-2009 09:51 AM
Alignment of ABI solid reads and 454 reads	baohua100	Bioinformatics	2	02-23-2009 05:58 PM

05-18-2011, 01:43 PM	#1
sulicon Member Location: Los Angeles Join Date: Aug 2010 Posts: 41	Spliced aligner for 454 reads? Hi all, I am trying to assemble 454 reads by aligning them to reference genome. I have tried the official gsMapper before, which just gave me exon sequences, rather than transcripts sequences. Now, I am using BWA-SW to align the reads, and Cufflinks/Scripture to reconstruct the transcripts. Still, I can only get the individual exon sequences, instead of expected transcript sequences. It seems the problem here is that BWA-SW is not a spliced aligner so that the splicing junction reads will be lost during the alignment. A spliced aligner, Tophat, is used in the tutorial of both Cufflinks and Scripture. However, this aligner is based on Bowtie, an aligner designed for short reads only. Could anyone give me some suggestion of the spliced aligner suitable for 454 reads? I think BLAT would be an option, but still want to test some other methods developed recently. Thanks, Shuli

05-18-2011, 02:41 PM	#2
lh3 Senior Member Location: Boston Join Date: Feb 2008 Posts: 693	For a read bridging a splice junction, bwasw should give two or more hits unless one of them is too short. Perhaps Cufflinks is expecting some tophat/bowtie specific information to group local hits to a transcript. I do not know. Nonetheless, I agree bwa-sw would not work well because it is a local aligner. For RNA-seq/ESTs, a dedicated splicing-aware glocal aligner is more appropriate. In addition to blat, you may also try gmap.

05-22-2011, 11:48 PM	#3
Jose Blanca Member Location: Valencia, Spain Join Date: Aug 2009 Posts: 70	gmap can align 454 ESTs against a genome taking into account the introns.

05-23-2011, 01:29 AM	#5
sphil Senior Member Location: Stuttgart, Germany Join Date: Apr 2010 Posts: 192	Blat is also able to handle 454-splice reads but be aware of long runtime....

05-24-2011, 02:14 PM	#6
sulicon Member Location: Los Angeles Join Date: Aug 2010 Posts: 41	Thanks Li Heng and Jose. I have tried GMAP and checked the output. It looks good. According the a document I have found, GMAP outperforms BLAT for gene structure identification in both speed and accuracy. And the latest version of GMAP can generate output in SAM format directly, which facilitates the subsequent analysis. Finally, I got some transcripts constructed by running Cufflinks on GMAP output. However, it seems that Cufflinks has modified the original alignment and generated some artificial exons/splicing junctions I've never seen in the GMAP output...