detection positive selection using multiple copy genes in orthologs
Entering edit mode
12 months ago
Mehmet ▴ 670

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?

genome sequence gene • 303 views
Entering edit mode

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.

Entering edit mode

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, on), the names don't match with the names in the tree. This is another problem.

Entering edit mode
11 months ago
h.mon 32k

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.


Login before adding your answer.

Traffic: 2785 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6