SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
hg19 genome reference for short read mapping yh253 Bioinformatics 4 12-29-2013 10:11 PM
Assisted de novo genome assembly? Create new consensus mapping reads to reference? zmartine Bioinformatics 8 02-10-2012 01:31 AM
Mapping reference genome to ensembl identifier bnfoguy Bioinformatics 0 06-13-2011 07:04 PM
mapping 454 reads to a reference genome query Bioinformatics 33 02-09-2011 07:36 AM
Reference genome for MAQ - split reference genome by chromosome or not? inesdesantiago Bioinformatics 4 02-18-2009 09:44 AM

Reply
 
Thread Tools
Old 11-04-2011, 05:34 AM   #1
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 838
Default Transcript mapping without a reference genome

How should I map transcript data in the absence of a good reference genome sequence?

I'm playing around with NGS data for an organism that has a fairly bad reference genome sequence (Schmidtea mediterranea, FWIW). In order to work out how bad, I split the reference sequence up to 100bp sequences, filtered out low-complexity sequences, then remapped back to the original sequence using bowtie2, and found quite a few places that mapped to >20 different contigs (~2000 by my rough guess using uniq/sort). This means that I'm not particularly comfortable using this genome for mapping our RNAseq data.

We now have heaps of RNAseq data and would like to look at differential expression and splice variants. We have some reference transcriptome sequences to map to, but due to some variation in reported sequences, can't be sure if two contigs from the transcriptome are different gene copies, or different isoforms of the same gene. So I have a few questions about the analysis:

Would it be better to create a new reference transcriptome from our data (about 450M reads, 35/50 paired-end on a SOLiD4, expected transcriptome size is about 20-30Mb), or use the previously published transcriptomes?

I like what cufflinks can do in identifying isoforms, but aren't sure how it responds to transcriptome mappings. Can it combine multiple contigs and notice that they are actually the same gene?

Is there any point in using tophat on a transcriptome, given that there shouldn't be any large breaks in the reference transcript sequences?

Other comments / questions would be appreciated, because this will be the first "real" analysis that I have done.
gringer is offline   Reply With Quote
Old 11-04-2011, 09:07 AM   #2
HESmith
Senior Member
 
Location: Bethesda MD

Join Date: Oct 2009
Posts: 509
Default

I would recommend Trinity, which is designed specifically for your application. The reference is http://www.nature.com/nbt/journal/v2.../nbt.1883.html.
HESmith is offline   Reply With Quote
Old 11-04-2011, 09:49 AM   #3
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 838
Default

Trinity would be nice, but it's not working so well on my little computer with 24GB memory. This seems a little strange to me because there is at least one genome assembler (Ray) that has a significantly lower memory footprint. I'll keep trying though (probably involving some code modification).

I have easy access to small clusters with distributed memory (I think something like 120GB over 5 nodes), but the Inchworm process is not distributed, so that's probably not going to help.
gringer is offline   Reply With Quote
Old 09-03-2020, 12:07 AM   #4
aalyaphillips
Junior Member
 
Location: USA

Join Date: Sep 2020
Posts: 1
Default

Moreover, the Third Generation Company from Japan started the production of game consoles recently named Famicom and made NES (Nintendo Entertainment System) game consoles. Beginning then on, support was begun to get the re-perceived and welcomed energetically until at last an early game was grown Super Mario Brothers. 7xRight now, amazing game comfort seemed to contend Sega and Nintendo. Jupiter Ed App After afierce rivalry among Sega and Nintendo, another contender, Sony, accompanied the Playstation. Sony supplanted cartridge with a CD and delivered most extreme outcomes and achievement. After Sony Play station got mainstream, the prevalence of Sega and Nintendo were slowly down. Beginning from age 6, sony playstation was formed into Sony Playstation 2, which utilizes DVD rather than CD.
aalyaphillips is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 01:00 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO