Singletons in MSA
0
0
Entering edit mode
5.0 years ago

I am reconstructing gene tree using Maximum Likelihood. Multiple Sequences Alignment of Bacterial strains consists of singletons. Should I remove them from alignment then reconstruct phylogeny or should I include them in alignment and then reconstruct tree?

gene sequence alignment • 2.2k views
ADD COMMENT
0
Entering edit mode

single tones

Do you mean "singletons"?

ADD REPLY
0
Entering edit mode

Sorry. Of course its singletons

ADD REPLY
0
Entering edit mode

I'm not sure what you mean by singletons in this context? Can you provide some example data or elaborate your question?

ADD REPLY
0
Entering edit mode

Singletons are sites which contains at least 2 types of nucleotides and occurs multiple times in an alignment.

ADD REPLY
0
Entering edit mode

In the context of multiple sequence alignment and phylogeny reconstruction, singleton are sites at which only one of the species has a distinct nucleotide, whereas all others are the same. For parsimony methods, they are not informative, so sometimes are removed - though it doesn't change the estimated phylogeny, it just saves computation cost. This is not the case for maximum-likelihood (and probably also Bayesian) methods.

ADD REPLY
0
Entering edit mode

Ah, think I’m understanding the question now.

I don’t profess to be an expert but to my knowledge, most ML techniques/tools ignore uninformative sites if necessary, so I don’t see a compelling reason to remove them - you’re just throwing away potentially real evolutionary signal as far as I can work out.

ADD REPLY
0
Entering edit mode

Multiple Sequences Alignment of Bacterial strains consists of singletons

Can you clarify what you mean by that? What exactly are you aligning? Genes, Whole genomes?

ADD REPLY
0
Entering edit mode

I am aligning genes present in different strains of bacteria. After alignment some sites were shown as singletons. I am using Maximum likelihood method for reconstructing tree. Dose these sites affects the phylogeny?

ADD REPLY

Login before adding your answer.

Traffic: 2969 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