Selection of Conserved Blocks from Multiple Alignments for Their Use in Phylogenetic Analysis

Abstract

The use of some multiple-sequence alignments in phylogenetic analysis, particularly those that are not very well conserved, requires the elimination of poorly aligned positions and divergent regions, since they may not be homologous or may have been saturated by multiple substitutions. A computerized method that eliminates such positions and at the same time tries to minimize the loss of informative sites is presented here. The method is based on the selection of blocks of positions that fulfill a simple set of requirements with respect to the number of contiguous conserved positions, lack of gaps, and high conservation of flanking positions, making the final alignment more suitable for phylogenetic analysis. To illustrate the efficiency of this method, alignments of 10 mitochondrial proteins from several completely sequenced mitochondrial genomes belonging to diverse eukaryotes were used as examples. The percentages of removed positions were higher in the most divergent alignments. After removing divergent segments, the amino acid composition of the different sequences was more uniform, and pairwise distances became much smaller. Phylogenetic trees show that topologies can be different after removing conserved blocks, particularly when there are several poorly resolved nodes. Strong support was found for the grouping of animals and fungi but not for the position of more basal eukaryotes. The use of a computerized method such as the one presented here reduces to a certain extent the necessity of manually editing multiple alignments, makes the automation of phylogenetic analysis of large data sets feasible, and facilitates the reproduction of the final alignment by other researchers.

Comparative transcriptomics across 14 Drosophila species reveals signatures of longevity

“…Briefly, the ortholog sequences were aligned with Clustal Omega v1.2.0 (Sievers et al., 2011) and concatenated gap‐free with Gblocks v0.91 (Castresana, 2000). The tree was constructed using neighbor‐joining method (Saitou & Nei, 1987) in Mega 6.06 (Tamura, Stecher, Peterson, Filipski & Kumar, 2013) and calibrated using the estimates of divergence time in the literature.…”

The first complete mitogenome of the South China deep‐sea giant isopod Bathynomus sp. (Crustacea: Isopoda: Cirolanidae) allows insights into the early mitogenomic evolution of isopods

“…Both the nucleotides and amino acids of the 13 protein‐coding genes were subjected to concatenated alignments using MEGA 5 (Tamura et al., 2011). Poorly aligned positions and divergent regions were removed using Gblocks version 0.91b (Castresana, 2000) set at the default block parameters. The final nucleotide and amino acid datasets consisted of 10,071 nt and 3,313 aa, respectively.…”

Comparative mitochondrial genomic analyses of three chemosynthetic vesicomyid clams from deep‐sea habitats

“…Summarizing, we included in our dataset 32 Veneroida, four Myoida, three Lucinoida, one Pholadomyoida, and the outgroup two Ostreoida. Each PCG, with the exception of atp8 , was separately aligned with codon‐based multiple alignments implemented in Mega7.0, and the gap and divergent regions were removed using Gblocks (version 0.91b) (Castresana, 2000) under default (stringent) settings. The atp8 gene was not included in the analyses due to the missing of this gene in several bivalves.…”