This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Comparing gene expression profiles between species

I am looking for input on ways to compare gene expression profiles between species. The scenario I am thinking of is:

  • I have a set of phenotypes in species A, and associated gene expression profiles.
  • Now I obtain a gene expression profile for species B.
  • Which phenotype in A looks most like the gene expression profile measured for species B?

For example, let's say we knock out a mouse gene and get a gene expression profile which is perturbed relative to the wild type. Which human disease looks most like this knock-out phenotype, based on gene expression?

I suppose this could be done on the gene level (by mapping genes to orthologs), on the pathway level (by assessing which pathways have been perturbed in each case), or on other levels of abstraction. Any success stories out there?

phenotype species gene-expression
• 6,820 views • link
updated by Ram Ram • written by Mikael Huss Mikael Huss

May be I can suggest you something, tell me what you have, microarray or RNA-seq?

• 0 views • link
updated by Ram Ram • written by Manvendra Singh Manvendra Singh

RNA-seq in this case, but I don't exclude suggestions or previous work based on microarrays!

• 0 views • link
written by Mikael Huss Mikael Huss

I have written my answer, but I think that your question is hot and should be thoroughly discussed on this forum, I am looking forward to have more suggestions/improvements/comments/criticisms

• 0 views • link
updated by Ram Ram • written by Manvendra Singh Manvendra Singh

1 answer

This is what I would for RNA-seq

  1. Map the reads on both genomes (human and mouse)
  2. Take only those reads for further analysis which mapped on both genomes (biasness of insertions and deletions of sequences between the genomes are removed and moreover, you get orthologous regions from your reads)
  3. Now count the reads over gene features and remove those genes which has low counts in all samples. (you would lose lot of them)
  4. Assign mean of counts over different transcripts to their respective gene, transform it on log scale.
  5. Now you have rownames as your genes colnames as your samples, now merge both species data into one dataframe
  6. Normalize them by their quantiles or surrogate variances.
  7. Calculate relative expression of each gene across the sample (assign the relative value to the rowmeans to each gene of each sample)
  8. Calculate spearman's correlation between the samples, and see which of them are forming clusters.

For microarray data I already have done and it was published in this paper, just read the methods section where cross-species and cross-platform expression analysis is mentioned, its also there in this publication as well.

HTH

• 0 views • link
updated by Ram Ram • written by Manvendra Singh Manvendra Singh

Thanks, that is useful. I will need to try this and think about it and perhaps continue the discussion here.

• 0 views • link
updated by Ram Ram • written by Mikael Huss Mikael Huss

Yes, I am very happy to be involved in discussion :)

• 0 views • link
updated by Ram Ram • written by Manvendra Singh Manvendra Singh

What if you have differential expression profiles (treated vs untreated) for each of say, 9 species? Could this method extrapolate to this instance?

• 1 views • link
updated by Ram Ram • written by rleach rleach

@Manvendra..what do you mean by 4) assign mean of counts over different transcripts to their respective gene, transform it on log scale?

• 0 views • link
updated by Ram Ram • written by ifudontmind_plzz ifudontmind_plzz

It's just mean of all transcripts to assign one value to gene

PS: Its quite old post, now this job is done by featureCount where we calculate counts on gene_id but not on transcript_id given in gtf files

• 1 views • link
updated by Ram Ram • written by Manvendra Singh Manvendra Singh

Log in to answer this question.

More posts like this