Question: Gene co-expression analysis
0
gravatar for mjoyraj
4.9 years ago by
mjoyraj50
Taiwan
mjoyraj50 wrote:

Hi, I have a FPKM matrix of 9644 columns and 15 rows from a RNA-seq expression experiment. The columns are the gene and the rows represents expression of the gene at different time points. My basic idea is to identify transcription factor binding site (TFBS) upstream of target gene. I want to do de-novo motif discovery based on over-represented sequence search in regulatory regions of target gene along with a group of co-expressed gene. Therefore, my first idea is to search genes co-expressed with my target gene. How can I do the co-expression analysis in R. How can I visualize the results? Please help me with R scripts. Let us consider the data name is ‘Rdata.csv’

rna-seq • 3.4k views
ADD COMMENTlink modified 4.9 years ago by EagleEye6.4k • written 4.9 years ago by mjoyraj50

Is it a time-series experiment?

ADD REPLYlink modified 4.9 years ago • written 4.9 years ago by Asaf6.1k

Yes, its a time-series experiment.

ADD REPLYlink written 4.9 years ago by mjoyraj50

joyraj_07@yahoo.co.in

ADD REPLYlink written 4.9 years ago by mjoyraj50
1
gravatar for biogirl
4.9 years ago by
biogirl180
European Union
biogirl180 wrote:

If it is time-series data, I would suggest using the TCAP method (temporal clustering by affinity propagation): http://bioinformatics.oxfordjournals.org/content/26/3/355.short

ADD COMMENTlink written 4.9 years ago by biogirl180

Thank you. I will go through the paper.

ADD REPLYlink written 4.9 years ago by mjoyraj50
1
gravatar for EagleEye
4.9 years ago by
EagleEye6.4k
Sweden
EagleEye6.4k wrote:

You can do time series gene expression using your target genes along with all other genes and look for your target gene in the time series profiles. Your target gene along with other genes will be clustered in particular profile or pattern of expression, which confirms the co-expressing genes with same pattern of expression as your target.

Use this program: STEM

http://www.cs.cmu.edu/~jernst/stem/

Hope that I understood your question properly.

 

ADD COMMENTlink written 4.9 years ago by EagleEye6.4k

Thank you. I will see the link and paper. Does any option is available to do the analysis in in R? I thought of finding the correlation between gene expression and filter top 1% correlations. Let us consider the top 1% correlations as set. The set will next be segregated to subsets based on mutual correlation. The genes within a subsets are co-expressed genes and gene ontology analysis will be done to predict if they are co-regulated or not. A set of co-regulated genes are believed to be regulated by same TF and their co-regulated genes are likely to have same motifs. Since my approach is de-novo motif discovery I would like to see whether the co-expressed/co-regulated genes have over represented sequence motif in their regulatory region. Since the motif discovery is all based on finding correct set of co-expressed/co-regulated gene, I want to go with the best available approaches. Please suggest, which approach out performs others. If any analysis is available in R-platform please help me with the pipeline and script. 

ADD REPLYlink written 4.9 years ago by mjoyraj50

I have no idea whether there is an R package to do it. But my suggestion:

1) First you can select the target gene along with co-expressed genes from STEM analysis (which also provides you with significant values). If the profile or pattern is similar with significant value means the genes in the profile are well correlated. 

2) You can use DAVID for your protein coding genes to get the GO terms or if you want command line tool with more extensive results you can use GeneSCF (only for Human): Gene Set Clustering based on Functional annotation (GeneSCF)

3) Use those Co-expressed genes sequence (there are many tools available to get sequence from chromosome coordinates, like blast).

3) Get all sequence and find the motif using the tools like MEME: http://meme.nbcr.net/meme/

ADD REPLYlink written 4.9 years ago by EagleEye6.4k

Thanks for your valuable suggestion. Can you please elaborate the step 2 and 3. My genes are from bird genome, so can I use the the said packages (as you mentioned its only for human). What I got from your answer is that:

1) Find co-expressed gene by STEM

2) Find ontology between the co-expressed gene

3) Filter co-expressed gene having ontology

4) Retrieve regulatory region of co-expressed gene having ontology

5) Analyse regulatory regions for motif discovery in MEME

If I am right, and got your points correctly, please let me know.

Many thanks again..

ADD REPLYlink written 4.9 years ago by mjoyraj50
0
gravatar for mjoyraj
4.9 years ago by
mjoyraj50
Taiwan
mjoyraj50 wrote:

Thank you. 

ADD COMMENTlink written 4.9 years ago by mjoyraj50
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: 2630 users visited in the last hour