Search
Question: Normalization for comparative transcriptomics between species
0
gravatar for jcfnielsen
5 months ago by
jcfnielsen0
jcfnielsen0 wrote:

I would like to compare expression levels of orthologous genes in different species for the purpose of conducting a co-expression analysis such as WGCNA.

I am using DESeq2 with the the following RNA-seq data:

6 different species (the species are from the same genus and phylogentically related)
2 different conditions
3 replicates

What would be the ideal experimental design to generate comparable expression levels between the orthologous genes? And How do I test whether the normalized expression levels are indeed comparable between species.

What I did so far:
- Used orthoMCL to identify single copy orthologous genes
- Generated a matrix with raw counts of all single copy orthologous genes for each of the 36 samples
- Normalized the counts using DESeq2, using the following design: ~ condition + species

I'm, thinking that this might be an issue since the model doesn't know that there is quite some variation in genome size and number of genes between the species, and it doesn't know the sequencing depth since only a subset of the RNA-seq data is being used.

Any suggestions on how to improve this?

ADD COMMENTlink modified 5 months ago by Michael Love15k • written 5 months ago by jcfnielsen0
0
gravatar for Michael Love
5 months ago by
Michael Love15k
United States
Michael Love15k wrote:
I don't have any great ideas for comparing RNAseq across species, it seems bioinformatically / conceptually difficult to me.
ADD COMMENTlink written 5 months ago by Michael Love15k

This kind of comparative analysis between species has been done before e.g. between ant species and between rat and human. However, nothing is stated in the methods section in these two paper about what they did to ensure proper normalization to account for e.g. varying genome size across species. Maybe they didn't do anything and maybe it doesn't have a great effect on the final result.

ADD REPLYlink written 5 months ago by jcfnielsen0
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.2.0
Traffic: 145 users visited in the last hour