SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
how to compute all restriction enzyme sites in the human genome? lunacab General 10 07-09-2017 04:36 PM
PubMed: Establishing the baseline level of repetitive element expression in the human Newsbot! Literature Watch 0 02-09-2012 03:00 AM
What is a Probe Set? What is a Probe pair? Pixel08 General 1 12-31-2011 02:15 PM
Human transcriptome array for high-throughput clinical studies NGSfan Literature Watch 0 02-24-2011 07:27 AM
targets of exon enrichment probe johnshembb General 2 02-18-2011 01:23 AM

Reply
 
Thread Tools
Old 10-20-2010, 09:31 AM   #1
tujchl
Member
 
Location: BEIJING, CHINA

Join Date: Sep 2009
Posts: 74
Default How to compute exression ratio of Human Exon 1.0 ST array on probe level

Hi guys:
Recently I search some chip data (Human Exon 1.0 ST array) from GEO, I found the annotation file on GEO is based on either probeset level (usually 4 probes) or even transcripts level. My question is how can I get expression value on probe level. I search on Bioconductor and want to find some packages. Unfortunately I don`t find any package suitable for me.
could you please recommend some packages or solutions for me?

thank you in advance
tujchl is offline   Reply With Quote
Old 10-20-2010, 05:32 PM   #2
malachig
Senior Member
 
Location: WashU

Join Date: Aug 2010
Posts: 117
Default Getting probe-level data out of Affymetrix exon arrays

If you have the original .CEL values, one option it to get RAW probe values out of these using 'Affymetrix Power Tools'. This software is a command line tool that works in windows (I know crazy...). It does work though. In particular the 'apt-cel-extract' function should allow you to get the individual probe level values. Probe level annotations are available here: GeneChip Human Exon 1.0 ST Array
malachig is offline   Reply With Quote
Old 10-20-2010, 05:42 PM   #3
malachig
Senior Member
 
Location: WashU

Join Date: Aug 2010
Posts: 117
Default

If you just want to re-process the data to the exon or gene level (including using your own custom annotation) then there are several options including: RMA Express, Expression Console, affy/BioConductor, Affymetrix Power Tools and others.
malachig is offline   Reply With Quote
Old 10-20-2010, 05:55 PM   #4
kenosaki
Member
 
Location: Tokyo, Japan

Join Date: May 2010
Posts: 27
Default

What data set did you download from? Tell us GEO accession # and downloaded dataset. It may help us to tell you how to get expression value.
kenosaki is offline   Reply With Quote
Old 10-20-2010, 06:29 PM   #5
tujchl
Member
 
Location: BEIJING, CHINA

Join Date: Sep 2009
Posts: 74
Default

Thank you malachig! Thank you kenosaki!
the GEO accession NO. is GSE12941. Its annotation file is GPL5175, you can see the results directly download from GEO is transcripts based. I want to see the ratio based on probe level.
tujchl is offline   Reply With Quote
Old 10-20-2010, 08:23 PM   #6
kenosaki
Member
 
Location: Tokyo, Japan

Join Date: May 2010
Posts: 27
Default

Ok, you can get CEL file. I usually download the CEL file when I want to probe level signals.
I like Expression Console to analyze it.

You need registration (free) to get the Expression Console. Once you installed it, you download library file and annotation file. Analysis will be run either Core, Extend, or Full. Choose Core, because Core probes are linked to the validated genes. Then RMA or whatever you like anlaysis will be done on probe level.
Export probe level result with certain annotations.
And also run gene level analysis and export gene level result.
Be carefull, "transcript ID of probe level result" is equal to "probe ID of gene level result".
Probe level results and Gene level results are now be able to be combined.

Well, this may not be the best short-cut for your purpose.
Take care.
kenosaki is offline   Reply With Quote
Old 10-20-2010, 10:22 PM   #7
tujchl
Member
 
Location: BEIJING, CHINA

Join Date: Sep 2009
Posts: 74
Default

Thank you kenosaki, you help me a lot. I have dowloaded Expression Console and Now installing library, I will try it, thank you again!
One more question~~
Exon array could be analyzed on probeset, exon and transcripts level respectively. To my knowledge each probeset , exon and transcripts DO have their own ID which can be serached on NetAffy.
My question is whether the ID is unique? (I mean whether ID could overlappe among probeset , exon and transcripts ID) and don`t change when version update (I mean one id could not point to different probe seq in different chip version).
In the GSE12941, Its probeset ID such as 3553690, could be transcript cluster ID NOT probeset ID when I search on NetAffy. could you please tell me why?

Finally, Thank you again

Last edited by tujchl; 10-20-2010 at 10:28 PM.
tujchl is offline   Reply With Quote
Old 10-20-2010, 10:59 PM   #8
tujchl
Member
 
Location: BEIJING, CHINA

Join Date: Sep 2009
Posts: 74
Default

One more question, kenosaki
when I load CEL file and Expression Console tell me to download "HuEx-1_0-st-v1.cdf" file (for Expression arrays) or "HuEx-1_0-st-v1.exon_analysis_configuration" (for Exon arrays). I just can find v2 on Affy, could you please help me again?
tujchl is offline   Reply With Quote
Old 10-20-2010, 11:37 PM   #9
tujchl
Member
 
Location: BEIJING, CHINA

Join Date: Sep 2009
Posts: 74
Default

malachig, I use the package RMAExpress, and do you know how to export expression value on probe level? Thank you
tujchl is offline   Reply With Quote
Old 10-21-2010, 01:57 AM   #10
kenosaki
Member
 
Location: Tokyo, Japan

Join Date: May 2010
Posts: 27
Default

Question 1,
A: Probeset ID is unique, as long as of the same chip. I don't think there is a change on different version. Every probe was designed from a transcript cluster (RefSeq RNA, EST,etc), so this cluster ID should be decided at the time of array design, it means they don't change at all. You better to double check to Affy.

For Question 2,
Yes I see the same error.
You have "HuEx-1_0-st-v2.exon_analysis_configuration" file in your library folder, and you need "HuEx-1_0-st-v1.exon_analysis_configuration" insted.
This file is not officially supported, so you need to ask Affy to get this file.
This is small 6kb file. I can send you by e-mail if you want.

Last edited by kenosaki; 10-21-2010 at 02:00 AM.
kenosaki is offline   Reply With Quote
Old 10-21-2010, 06:52 PM   #11
tujchl
Member
 
Location: BEIJING, CHINA

Join Date: Sep 2009
Posts: 74
Default

It is very kinds of you to offer me the annotation file-_-, My e-mail is delete[/email].
a little question about Expression Console:
After I load the CEL file, there are only two categories of available analysis methods: Gene Level and Exon Level (I choose exon level core), whithout probe level. After analysis, I can only export the results based on probesets....
could you please tell me more about operation?

Thank you

Last edited by tujchl; 10-21-2010 at 10:54 PM.
tujchl is offline   Reply With Quote
Old 10-21-2010, 09:58 PM   #12
kenosaki
Member
 
Location: Tokyo, Japan

Join Date: May 2010
Posts: 27
Default

I sent the conf file to you by e-mail. Locate this file into the libray folder of expression console
I suggest that you delete your e-mail address from this thread now.

By the way, exon level means probe level, as you expected.
Operation? Depends on what you want to do. Could you look at Help or tutorial that Affy provided? It helps the operation pretty much.
I know they have other technical documents about exon array.
And I believe they support you.

If you can't find any help sitll, feel freet to ask me
kenosaki is offline   Reply With Quote
Old 10-21-2010, 10:54 PM   #13
tujchl
Member
 
Location: BEIJING, CHINA

Join Date: Sep 2009
Posts: 74
Default

OK, I have the the e-mail, Thank you very very much
tujchl is offline   Reply With Quote
Old 12-22-2010, 12:30 AM   #14
honey
Senior Member
 
Location: Pittsburgh

Join Date: Feb 2010
Posts: 151
Default

I was wondering if from the exon array can I get exon level expression analysis or transcript level. I am interested in more of exon level expression.
Has anyone successfully completed such analysis with exon arrays.

Thanks
honey is offline   Reply With Quote
Old 10-30-2012, 08:30 AM   #15
gargvin2003
Junior Member
 
Location: south carolina

Join Date: Oct 2012
Posts: 1
Default

Hi kenosaki,
Could you please send me the "HuEx-1_0-st-v1.exon_analysis_configuration" file. I am having the same problem with opening the "CEL" file in expression console.

gargvin2003@aol.in
gargvin2003 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 03:39 AM.


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