Question 20 of 20 Advanced Troubleshooting
A gene tree built from your alignment gives an unexpected topology. Which alignment problems would you check?
Short answer
Check the input set for paralogs, contaminants and fragments, the alignment for misaligned blocks from early progressive errors, and the columns for gap-rich or saturated regions. Then compare trees built with another aligner or trimming setting before believing the surprise.
Why
An unexpected tree can be real biology, but alignment problems are more common causes. First, the input. If a paralog sits among orthologs, the tree mixes a duplication with speciation and species appear in the wrong place, the error the OrthoFinder authors describe for best-hit methods.
Fragments and contaminant sequences also draw odd branches. Second, the alignment itself. The Clustal Omega paper notes that progressive alignment is greedy, so an early mistake is never corrected, which can misplace whole blocks in divergent members.
Inspecting known motifs often reveals this. Third, the columns. Poorly aligned, gap-rich regions add noise, and the trimAl authors note that removing poorly aligned regions improves later analyses.
Trimming choices can change topology, so test more than one. A useful robustness check is to rebuild the tree from alignments made with a different aligner and different trimming. Splits that persist across choices deserve belief.
Splits that appear with only one combination are artifacts until shown otherwise, and low support values should be reported as unresolved rather than interpreted.
A strong answer covers
- Paralogs and contaminants distort topology
- Progressive alignment errors persist and misplace blocks
- Gap-rich, poorly aligned columns add noise
- Test robustness with other aligners and trimming
- Report low-support splits as unresolved
Common mistakes
- Publishing a surprising topology from a single alignment and setting
- Ignoring possible paralogs in the input set
- Interpreting branches with low support as evidence
The interviewer might ask next
How would you tell real horizontal gene transfer from an alignment artifact?