Comparison of Bayesian and maximum likelihood bootstrap measures of phylogenetic reliability

Owing to the exponential growth of genome databases, phylogenetic trees are now widely used to test a variety of evolutionary hypotheses. Nevertheless, computation time burden limits the application of methods such as maximum likelihood nonparametric bootstrap to assess reliability of evolutionary trees. As an alternative, the much faster Bayesian inference of phylogeny, which expresses branch support as posterior probabilities, has been introduced. However, marked discrepancies exist between nonparametric bootstrap proportions and Bayesian posterior probabilities, leading to difficulties in the interpretation of sometimes strongly conflicting results. As an attempt to reconcile these two indices of node reliability, we apply the nonparametric bootstrap resampling procedure to the Bayesian approach. The correlation between posterior probabilities, bootstrap maximum likelihood percentages, and bootstrapped posterior probabilities was studied for eight highly diverse empirical data sets and were also investigated using experimental simulation. Our results show that the relation between posterior probabilities and bootstrapped maximum likelihood percentages is highly variable but that very strong correlations always exist when Bayesian node support is estimated on bootstrapped character matrices. Moreover, simulations corroborate empirical observations in suggesting that, being more conservative, the bootstrap approach might be less prone to strongly supporting a false phylogenetic hypothesis. Thus, apparent conflicts in topology recovered by the Bayesian approach were reduced after bootstrapping. Both posterior probabilities and bootstrap supports are of great interest to phylogeny as potential upper and lower bounds of node reliability, but they are surely not interchangeable and cannot be directly compared.     

A contribution to sedentary polychaete phylogeny using 18S rRNA sequence data

The phylogenetic position of Annelida as well as its ingroup relationships are a matter of ongoing debate. A molecular phylogenetic study of sedentary polychaete relationships was conducted based on 70 sequences of 18S rRNA, including unpublished sequences of 18 polychaete species. The data set was analysed with maximum parsimony and maximum likelihood methods. Clade robustness was estimated by parsimony-bootstrapping and jackknifing, decay index, and clade support, as well as a posteriori probability tests using Bayesian inference. Irrespective of the applied method, some traditional sedentary polychaete taxa, such as Cirratulidae, Opheliidae, Orbiniidae, Siboglinidae and Spionidae, were recovered by our phylogenetic reconstruction. A close relationship between Orbiniidae and Questa received a particularly strong support. Echiura appears to be a polychaete ingroup taxon which is closely related to Dasybranchus (Capitellidae). As in previous molecular analyses, no support was found for the monophyly of Annelida nor for that of Polychaeta. However, we suggest that an increase in taxon sampling may yield additional resolution in the reconstruction of polychaete ingroup phylogeny, although the difficulties in reconstructing the basal phylogenetic relationships within Annelida may be due to their rapid radiation.     

Inference of bacterial microevolution using multilocus sequence data

We describe a model-based method for using multilocus sequence data to infer the clonal relationships of bacteria and the chromosomal position of homologous recombination events that disrupt a clonal pattern of inheritance. The key assumption of our model is that recombination events introduce a constant rate of substitutions to a contiguous region of sequence. The method is applicable both to multilocus sequence typing (MLST) data from a few loci and to alignments of multiple bacterial genomes. It can be used to decide whether a subset of isolates share common ancestry, to estimate the age of the common ancestor, and hence to address a variety of epidemiological and ecological questions that hinge on the pattern of bacterial spread. It should also be useful in associating particular genetic events with the changes in phenotype that they cause. We show that the model outperforms existing methods of subdividing recombinogenic bacteria using MLST data and provide examples from Salmonella and Bacillus. The software used in this article, ClonalFrame, is available from http://bacteria.stats.ox.ac.uk/.     

A multigene phylogeny of the Gelidiales including nuclear large-subunit rRNA sequence data

Partial sequences (1032 bp) of the nuclear-encoded large ribosomal RNA gene (LSU) were determined for 16 gelidialean species, and analyzed separately and in combination with plastid rbcL and nuclear SSU gene sequences. The number of informative characters and levels of sequence divergence among taxa are intermediate in LSU sequences as compared to that for rbcL and SSU. Analyses of the separate LSU, and a combined LSU, SSU, and rbcL data sets have identified early-diverging lineages within the Gelidiales including Gelidiella, Pterocladia, Pterocladiella, and a lineage including Gelidium and species classified in other genera. The relationships among most gelidialean taxa are well-resolved and well-supported by analyses of the combined data; however, the relationships of Ptilophora and Capreolia remain unclear. It is speculated that these two lineages have diverged from a common ancestor over an evolutionarily short period of time. This revised version was published online in June 2006 with corrections to the Cover Date.     

A method for molecular phylogeny construction by direct use of nucleotide sequence data

Summary A method for molecular phylogeny construction is newly developed. The method, called the stepwise ancestral sequence method, estimates molecular phylogenetic trees and ancestral sequences simultaneously on the basis of parsimony and sequence homology. For simplicity the emphasis is placed more on parsiomony than on sequence homology in the present study, though both are certainly important. Because parsimony alone will sometimes generate plural candidate trees, the method retains not one but five candidates from which one can then single out the final tree taking other criteria into account.The properties and performance of the method are then examined by simulating an evolving gene along a model phylogenetic tree. The estimated trees are found to lie in a narrow range of the parsimony criteria used in the present study. Thus, other criteria such as biological evidence and likelihood are necessary to single out the correct tree among them, with biological evidence taking precedence over any other criterion. The computer simulation also reveals that the method satisfactorily estimates both tree topology and ancestral sequences, at least for the evolutionary model used in the present study.     

Molecular phylogeny of Pottiaceae (Musci) based on chloroplast rps4 sequence data

Comparative sequencing of the chloroplast rps4 gene was used to reconstruct the phylogenetic relationships within the family Pottiaceae (Musci). The results confirm that Ephemerum spinulosum, Splachnobryum obtusum, Goniomitrium acuminatum and Cinclidotus fontinaloides are clearly positioned within the Pottiaceae and that Hypodontium dregei is not a member. At subfamily level, the data support the subfamily Pottioideae as being a monophyletic clade. The Trichostomoideae are probably paraphyletic. Neither the subfamily Chionolomideae, represented in this study by Pseudosymblepharis schimperiana, nor the subfamily Erythrophyllopsoideae, represented by both known species, Erythrophyllastrum andinum and Erythrophyllopsis fuscula, are supported by the sequence data. The Timmielloideae should be excluded from the Pottiaceae. The Merceyoideae, represented in this study by Scopelophila cataractae, might form a sister clade to all other Pottiaceae, but their position is not fully resolved. At the genus level, Barbula is clearly polyphyletic since Barbula bolleana and Barbula indica appear in a clade clearly separated from Barbula unguiculata. Pottiopsis caespitosa and Leptobarbula berica are placed within the Trichostomoideae. Likewise, the genera Gymnostomum and Anoectangium are excluded from the Pottioideae and placed within the Trichostomoideae. Leptophascum leptophyllum is closely related with Syntrichia; Aloina is not closely related to Tortula or Crossidium. Evidence of a clade within the Pottioideae, formed of Leptodontium and Triquetrella, is provided.     

Reconstruction of molecular phylogeny of extant hominoids from DNA sequence data

Evolutionary distance matrices of the extant hominoids are computed from DNA sequence data, and hominoid DNA phylogenies are reconstructed by applying the neighbor-joining method to these distance matrices. The chimpanzee is clustered with the human in most of the phylogenetic trees thus obtained. The proportion of the distance between human and chimpanzee to that between human/chimpanzee and orangutan is estimated. Both mitochondrial DNA and nuclear DNA show a similar value (0.44), which is close to values derived from DNA-DNA hybridization data.     

Prokaryote phylogeny based on ribosomal proteins and aminoacyl tRNA synthetases by using the compositional distance approach

In order to show that the newly developed K-string composition distance method, based on counting oligopeptide frequencies, for inferring phylogenetic relations of prokaryotes works equally well without requiring the whole proteome data, we used all ribosomal proteins and the set of aminoacyl tRNA synthetases for each species. The latter group has been known to yield inconsistent trees if used individually. Our trees are obtained without making any sequence alignment. Altogether 16 Archaea, 105 Bacteria and 2 Eucarya are represented on the tree. Most of the lower branchings agree well with the latest, 2003, Outline of the second edition of the Bergey’s Manual of Systematic Bacteriology and the trees also suggest some relationships among higher taxa.     

Revaluation of deuterostome phylogeny and evolutionary relationships among chordate subphyla using mitogenome data

The traditional knowledge in textbooks indicated that cephalochordates were the closest relatives to vertebrates among all extant organisms. However, this opinion was challenged by several recent phylogenetic studies using hundreds of nuclear genes. The researchers suggested that urochordates, but not cephalochordates, should be the closest living relatives to vertebrates. In the present study, by using data generated from hundreds of mtDNA sequences, we revalue the deuterostome phylogeny in terms of whole mitochondrial genomes （mitogenomes）. Our results firmly demonstrate that each of extant deuterostome phyla and chordate subphyla is monophyletic. But the results present several alternative phylogenetic trees depending on different sequence datasets used in the analysis. Although no clear phylogenetic relationships are obtained, those trees indicate that the ancient common ancestor diversified rapidly soon after their appearance in the early Cambrian and generated all major deuterostome lineages during a short historical period, which is consistent with ＂Cambrian explosion＂ revealed by paleontologists. It was the 520-million-year＇s evolution that obscured the phylogenetic relationships of extant deuterostomes. Thus, we conclude that an integrative analysis approach rather than simply using more DNA sequences should be employed to address the distant evolutionary relationship.     

Estimating the phylogeny in mollusc Littorina saxatilis (Olivi) from enzyme data: methodological considerations

The evolutionary history of 19 populations of Littorina saxatilis (Olivi) was estimated by four different approaches. Three of these operate upon a population by population matrix of genetic distances: average linkage clustering, and two versions of the Fitch-Margoliash method. The fourth method was a maximum likelihood estimate based on differences in allele frequencies between populations. The study aims to assess how well each method estimates the phylogeny by including seven populations of the closely related species L. arcana Hannaford Ellis. The rationale behind this is that a good estimation technique should be able to separate these two monophyletic taxa.The results show that, by our criteria, the maximum likelihood method yields the best estimate and the unconstrained Fitch-Margoliash technique gives reasonable estimates. Both average-linkage clustering and the Fitch-Margoliash method with evolutionary clock perform less well. We argue that this is expected since both these techniques are based on probably unrealistic assumptions such as the overall rate of evolutionary divergence being homogeneous over phyletic lines.     

Molecular phylogeny and evogenomics of heterocystous cyanobacteria using rbcl gene sequence data

Taxonomic affiliations and molecular diversity of 41 heterocystous cyanobacteria representing 12 genera have been assessed on an evolutionary landscape using rbcl gene sequence data-based phylogenomics and evogenomics approaches. Phylogenetic affiliations have clearly demonstrated the polyphyly of the true branching cyanobacteria, along with a frequent intermixing amongst the heterocystous cyanobacteria. The monophyletic origin of the heterocystous cyanobacteria was also quite evident from maximum parsimony and neighbor joining analyses. Incongruency with the traditional scheme of cyanobacterial taxonomy was frequently observed, thus advocating towards some re-amendments in the cyanobacterial classificatory schemes. Evogenomics analyses of gene sequence data gave a clear indication about the greater evolutionary pace of the unbranched cyanobacteria as compared to the branched forms. It was evident that the order Nostocales would be controlling the future pace of evolution of heterocystous cyanobacteria. The cyanobacteria Nostoc was found to have the greatest genetic heterogeneity amongst the studied genera, along with some evidence towards events of lateral gene transfer amongst the heterocystous cyanobacteria in case of the rbcl gene. Thus, heterocystous cyanobacteria were found to be a fast evolving group, with estimates of gene conversion tracts pointing towards the unbranched heterocystous cyanobacteria being at the base of evolutionary diversifications of the complete heterocystous lineage.     

Mitochondrial sequence data reveals the phylogeny of the Asian Helice group of crabs (Decapoda: Brachyura: Varunidae)

The taxonomy of the Helice s. l. crabs has been recently revised, subdividing the Asian species in three genera, Helice , Helicana and Pseudohelice . The last one is widespread throughout the Indo-West Pacific Ocean and a single species, Pseudohelice subquadrata , is found in Asia. The other two genera appear to be endemic to Asia. However, phylogenetic affinities among the Asian representatives of Helice s. l. are poorly understood, and different morphology-based hypotheses of relationships remain untested. In this study, we used sequence data from a mitochondrial gene-block spanning three different genes [16S rDNA- trnL ₁ - IGS -nad1 ] to perform phylogenetic analyses. The results provide independent molecular evidence supporting the splitting of the Asian species into Helice and Helicana . Helice is suggested to be paraphyletic with respect to P. subquadrata . In addition, the heterogeneity of the genus Pseudohelice is shown when including the East Asian and East African counterparts. We propose that the two genera, Helice and Pseudohelice , are not monophyletic in their current composition.     

The phylogeny of Chinese indigenous pig breeds inferred from microsatellite markers

A genetic study of 32 local Chinese, three foreign pig breeds [Duroc (DU), Landrace and Yorkshire], and two types of wild boar (Hainan and Dongbei wild boar) based on 34 microsatellite loci was carried out to clarify the phylogeny of Chinese indigenous pig breeds. The allele frequencies, effective numbers of alleles, and the average heterozygosity within populations were calculated. The results showed that the genetic variability of the Lingao pig was the largest, while the Jiaxing pig was the lowest. The greatest distance between domestic pigs was found between Shanggao and DU pig and the shortest was found between Wuzhishan and Lingao pig, respectively. A neighbour-joining tree constructed from Modified Cavalli-Sforza genetic distances divided Chinese pigs into two clusters; four subclusters were also identified. Our results only partly agree with the traditional types of classification and also provide a new relationship among Chinese local pig breeds. Our data also confirmed that Chinese pig breeds have a different origin from European/American breeds and can be utilized in programmes that aim to maintain Chinese indigenous pig breeds.     

Molecular systematics of fifteen old world cichlid species using allozyme data

The evolutionary relationships among 15 cichlid species, from six different genera, were investigated using allozyme data from 40 protein coding loci in an attempt to elucidate relationships between species which could previously not be established on morphological grounds. Dendrograms were constructed by using both genetic distances and binary coded presence or absence of alleles. The results obtained using cladistical methods and consensus trees were more consistent with hypothesized morphological-based phylogenies, compared to results obtained using phenetic methods. It is suggested that the predatory and substratum-spawning Hemichromis elongatus (Guichenot) is separated from the herbivorous (and substratum-spawning) Tilapia species and that the latter group of fish are separated from the mouth-brooding and omnivorous Oreochromis species. It is also evident that the above-mentioned genera are consistently separated from the mouthbrooding genera Chetia, Pharyngochromis and Serranochromis, which are predators.     

A Poisson model of sequence comparison and its application to coronavirus phylogeny

In this paper, we propose two metrics to compare DNA and protein sequences based on a Poisson model of word occurrences. Instead of comparing the frequencies of all fixed-length words in two sequences, we consider (1) the probability of ‘generating’ one sequence under the Poisson model estimated from the other; (2) their different expression levels of words. Phylogenetic trees of 25 viruses including SARS-CoVs are constructed to illustrate our approach.     

Bayesian estimation of movement and survival probabilities from capture-recapture data

The Arnason–Schwarz model is usually used for estimatingsurvival and movement probabilities of animal populations fromcapture-recapture data. The missing data structure of this capture-recapturemodel is exhibited and summarised via a directed graph representation.Taking advantage of this structure we implement a Gibbs samplingalgorithm from which Bayesian estimates and credible intervalsfor survival and movement probabilities are derived. Convergenceof the algorithm is proved using a duality principle. We illustrateour approach through a real example.     

The phylogeny of Patrinieae sensu Graebner (Valerianaceae) revisited: additional evidence from ndhF sequence data

 DNA sequences of both 5′ and 3′ regions of the plastid ndhF gene were generated in order to study the position of Patrinia and Nardostachys, to check the potential paraphyletic nature of Patrinieae, and to evaluate the possible link between the tribe and Linnaeaceae. Parsimony analysis showed very strong support for Patrinia as sister to all members of Valerianaceae (including Nardostachys) and indicated the paraphyletic nature of the tribe Patrinieae. Additionally, trees were constructed from available rbcL data separately and supplemented with ndhF sequences. Topologies of these combined cladograms are in agreement with the ndhF phylogeny, suggesting that the traditionally circumscribed Patrinieae can no longer be recognized but must be considered as part of a basal grade in Valerianaceae. Parsimony analysis based on a morphological data set supported a monophyletic Patrinieae; combination with the molecular data showed a paraphyletic Patrinieae. Furthermore, the possible link between Patrinieae and Linnaeaceae is evaluated. Received July 12, 2001 Accepted February 25, 2002     

Population genetic analysis of Pseudomonas aeruginosa using multilocus sequence typing

To study the population genetic structure of Pseudomonas aeruginosa, we developed a multilocus sequence typing scheme. The sequences of internal fragments of seven housekeeping genes were obtained for 34 P. aeruginosa isolates from patients hospitalized in five different European cities. Twenty-six different allelic profiles were identified. The mean allelic diversity was 0.854 (range: 0.606-0.978), which was about six times greater than the results obtained with the multilocus enzyme electrophoresis method. Linkage disequilibrium was measured with the index of association. An index of 1.95+/-0.24 was calculated when all the strains were considered. This index was 1.76+/-0.27 when only one strain per sequence type was considered. Both results were different from 0, indicating linkage among loci, which means that the population structure of our set of P. aeruginosa isolates is clonal. The clonal structure of the population was also suggested by the congruence of the topology of the different trees obtained from the seven housekeeping genes. These results are in contrast to previous studies, finding a non clonal population structure. Since a small number of isolates was analyzed in this study, there might be a bias of selection which includes the possibility that they belong to widely disseminated epidemic clones. Another possibility is that recombination did not occurred homogeneously throughout the genome of P. aeruginosa, so that part of it has a clonal structure, while the remaining part of the genome is more frequently subject to recombination.