Question: detection positive selection using multiple copy genes in orthologs
0
gravatar for Mehmet
7 weeks ago by
Mehmet540
Japan
Mehmet540 wrote:

Dear all,

I have done positive selection analysis using one-to-one copy ortholog groups via codeml tool. But I have seen some papers in which multiple copy gene ortholog families were also used to test positive selection. When I tried to do same analysis (multiple copy approach), ETE tool used only one sequence (one gene) from each species although each species has multiple copies in the same ortholog family.

Is there any way to deal with gene duplication event because ETE uses only the first sequence of each species?

Codeml in PAML says Seq occurs more than once in the alignment.

Any of you has done this kind of analysis (using multiple copy orthologs) before?

sequence gene genome • 115 views
ADD COMMENTlink modified 7 weeks ago by h.mon30k • written 7 weeks ago by Mehmet540

Are the sequence names duplicated? Or are the names unique? Are there repeated sequences? ETE and codeml doesn't care about species, but they probably don't like duplicated names and / or duplicated sequences.

ADD REPLYlink written 7 weeks ago by h.mon30k

Sequence names are duplicated because the sequence names must match with the names in the tree. For instance;

in the sequence file; names are :

SpeciesA SpeciesA SpeciesA

SpeciesB SpeciesB

This is the name order in the tree as well. (SpeciesA, SpeciesB so on). When I use sequence name as different (SpeciesA_seq1, SpeciesA_seq2...so on), the names don't match with the names in the tree. This is another problem.

ADD REPLYlink modified 7 weeks ago • written 7 weeks ago by Mehmet540
0
gravatar for h.mon
7 weeks ago by
h.mon30k
Brazil
h.mon30k wrote:

Each gene needs a unique identifier, and alignment and tree gene names must match. So your genes should be named SpeciesA_seq1, SpeciesA_seq2, SpeciesB_seq1, and so on, at both the nucleotide alignment file and tree file.

ADD COMMENTlink written 7 weeks ago by h.mon30k
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: 1518 users visited in the last hour