SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Bowtie option --best --strata for unique mapper in the genome StephaniePi83 Bioinformatics 1 03-27-2012 01:29 PM
Long Read Mapper mitochy Bioinformatics 9 03-22-2012 05:33 AM
454 Mapper nitinkumar 454 Pyrosequencing 3 05-31-2011 12:57 PM
Mapper with annotation support dbrami Bioinformatics 1 07-21-2010 02:14 PM
Mapper to extned contigs dbrami Bioinformatics 1 07-02-2010 05:10 AM

Reply
 
Thread Tools
Old 08-07-2012, 12:58 AM   #1
StephaniePi83
Member
 
Location: France

Join Date: Sep 2011
Posts: 52
Default Unique mapper on genome - multi mapper on transcriptome

Hi everybody,
I'm analyzing iCLIP data (technique to study protein-RNA interactions on a genome-wide scale). I have two populations bound by the protein of interest - small RNA and mRNA.
I take all reads that map uniquely to the genome and i wanted to know all transcript that are bound by our protein of interest.
But i have a problem, i use bowtie to map unique mapper reads on the genome tothe transcriptome, but unique mapper on the genome became multi mapper on transcriptome, i'm confused, how to deal with it ? I have some very small reads ( the smallest are 13 nt).
How can i assign reads to transcript without mistake, i mean how to deal with reads aligning to multiple genes ?
do i have to associate each read with all transcripts it aligns to?
Thanks in advance for your help.

Last edited by StephaniePi83; 08-07-2012 at 01:49 AM.
StephaniePi83 is offline   Reply With Quote
Old 08-29-2012, 10:32 PM   #2
paulsydney
Junior Member
 
Location: Sydney

Join Date: Nov 2011
Posts: 6
Default

Hi Stephanie,

I wonder if you have published your iCLIP study recently. I am keen to learn about how to analyse iCLIP data.
paulsydney is offline   Reply With Quote
Old 08-30-2012, 03:48 AM   #3
StephaniePi83
Member
 
Location: France

Join Date: Sep 2011
Posts: 52
Default

Hi paulsydney,
we haven't published yet.
Some paper help me :
http://www.ncbi.nlm.nih.gov.gate1.in...ubmed/21048981
http://www.ncbi.nlm.nih.gov.gate1.in...ubmed/20601959
http://www.ncbi.nlm.nih.gov.gate1.in...ubmed/21358640
StephaniePi83 is offline   Reply With Quote
Old 08-30-2012, 06:50 PM   #4
paulsydney
Junior Member
 
Location: Sydney

Join Date: Nov 2011
Posts: 6
Default

Thanks Stephanie. I am still trying to decode many parts of the analysis section of that paper. Did you end up using any other bioinformatics tools apart from bowtie? I am not to certain how to extract annotation data from my map.
paulsydney is offline   Reply With Quote
Old 08-31-2012, 12:10 AM   #5
StephaniePi83
Member
 
Location: France

Join Date: Sep 2011
Posts: 52
Default

After removing adapter from my reads, i mapped all read on the genome allowing up to 2 mismatches with bowtie. I translate bowtie output file to bed file; then i use Bedtools (intersectBed function) to annotate all reads.
StephaniePi83 is offline   Reply With Quote
Old 09-04-2012, 12:44 AM   #6
paulsydney
Junior Member
 
Location: Sydney

Join Date: Nov 2011
Posts: 6
Default

Hi Stephanie,

What tools did you use to randomise iCLIP position when doing the pentamer z-score analysis?
paulsydney is offline   Reply With Quote
Old 09-04-2012, 12:55 AM   #7
StephaniePi83
Member
 
Location: France

Join Date: Sep 2011
Posts: 52
Default

Hi Paul,
I didn't do this analysis neither the statistical analysis. My supervisor send our data to a person who did a pipeline to analyse iclip data (he did all bioinformatic analysis in the paper i send you).

Last edited by StephaniePi83; 09-04-2012 at 01:44 AM.
StephaniePi83 is offline   Reply With Quote
Old 09-04-2012, 01:27 AM   #8
kopi-o
Senior Member
 
Location: Stockholm, Sweden

Join Date: Feb 2008
Posts: 319
Default

Unless I'm misunderstanding something, there isn't anything weird about unique genome mappers that are also multiple transcriptome mappers. Consider a read that spans a splice junction and can be uniquely mapped to the genome. If the gene containing the splice junction has multiple isoforms containing that splice junction, the read will map to multiple isoforms in the transcriptome reference. Does that answer the question?
kopi-o is offline   Reply With Quote
Old 09-04-2012, 01:30 AM   #9
kopi-o
Senior Member
 
Location: Stockholm, Sweden

Join Date: Feb 2008
Posts: 319
Default

Unless I'm misunderstanding something, there isn't anything weird about unique genome mappers that are also multiple transcriptome mappers. Consider a read that spans a splice junction and can be uniquely mapped to the genome. If the gene containing the splice junction has multiple isoforms containing that splice junction, the read will map to multiple isoforms in the transcriptome reference. Does that answer the question?
kopi-o is offline   Reply With Quote
Old 09-04-2012, 01:32 AM   #10
paulsydney
Junior Member
 
Location: Sydney

Join Date: Nov 2011
Posts: 6
Default

Ah ic... That was a great idea. I should suggest that to my supervisor too actually.
paulsydney is offline   Reply With Quote
Old 09-04-2012, 01:49 AM   #11
StephaniePi83
Member
 
Location: France

Join Date: Sep 2011
Posts: 52
Default

Hi kopi-o,
This is what i was thinking ! it map to all transcript of the same gene on transcriptome, and some reads will map to exon-exon juction that are not present on the genome.
StephaniePi83 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 05:38 AM.


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