dN/dS ratio (omega)
2
0
Entering edit mode
6.5 years ago
dav • 0

So, I have a group of genes (proteins) - that arent necesarely related and I am interested in dn/ds ratio, or rather the conservation. They are all from the same species (A.thaliana) and I have data on around 600 individuals (mostly from 1001genomes).

What software should I stick to?

I was thinking of PAML (Phylogenetic Analysis by Maximum Likelihood), more specifically codeml. I would give it aligned sequences, so 600 sequences of the same protein and a tree that goes with it, but I saw people say it is best suited for more distant species.

Is my approach wrong? What should I stick to for this?

dn-ds-ratio PAML VCF selection sequence • 4.9k views
ADD COMMENT
2
Entering edit mode

if they are not related (i.e. not homologous) you CANNOT use phylogenetic based methods

ADD REPLY
4
Entering edit mode
6.5 years ago
natasha.sernova ★ 3.9k

You are on the correct track, but PAML description is very hard.

This is the manual.

I would recommend to read all about codeml program (pp. from 28th to 38th.). You will understand what these parameters represent and all that is necessary to specify to run PAML.

I suggest using runmode=-2. What this mode means is described on p. 34.

Good luck!

Natasha

ADD COMMENT
0
Entering edit mode

Ah, so runmode=-2 gives me the exact lineage ratio I need!

Thank you a lot!

ADD REPLY
0
Entering edit mode

But wait , if I run it pairwise , (runmode -2) , then it compares all sequences to the first one?

Or will it compare 3rd to 2nd , 4th to 3rd etc?

Should I leave runmode = 0 to get it globaly?

ADD REPLY
0
Entering edit mode
6.5 years ago
natasha.sernova ★ 3.9k

Dear dav,

I think the first url contains the answer to your question about Dn/Ds for a global alignment.

This file is about runmode-option.

In any case, PAML is a very complicated program as well as its manual.

Below there are a manual how to read PAML-description and other useful urls

Good luck!

Natasha

ADD COMMENT

Login before adding your answer.

Traffic: 1364 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6