Bayesian coestimation of phylogeny and sequence alignment

Background  

Two central problems in computational biology are the determination of the alignment and phylogeny of a set of biological sequences. The traditional approach to this problem is to first build a multiple alignment of these sequences, followed by a phylogenetic reconstruction step based on this multiple alignment. However, alignment and phylogenetic inference are fundamentally interdependent, and ignoring this fact leads to biased and overconfident estimations. Whether the main interest be in sequence alignment or phylogeny, a major goal of computational biology is the co-estimation of both.     

Incorporating indel information into phylogeny estimation for rapidly emerging pathogens

Background  

Phylogenies of rapidly evolving pathogens can be difficult to resolve because of the small number of substitutions that accumulate in the short times since divergence. To improve resolution of such phylogenies we propose using insertion and deletion (indel) information in addition to substitution information. We accomplish this through joint estimation of alignment and phylogeny in a Bayesian framework, drawing inference using Markov chain Monte Carlo. Joint estimation of alignment and phylogeny sidesteps biases that stem from conditioning on a single alignment by taking into account the ensemble of near-optimal alignments.     

A practical approach to phylogenomics: the phylogeny of ray-finned fish (Actinopterygii) as a case study

Background  

Molecular systematics occupies one of the central stages in biology in the genomic era, ushered in by unprecedented progress in DNA technology. The inference of organismal phylogeny is now based on many independent genetic loci, a widely accepted approach to assemble the tree of life. Surprisingly, this approach is hindered by lack of appropriate nuclear gene markers for many taxonomic groups especially at high taxonomic level, partially due to the lack of tools for efficiently developing new phylogenetic makers. We report here a genome-comparison strategy to identifying nuclear gene markers for phylogenetic inference and apply it to the ray-finned fishes – the largest vertebrate clade in need of phylogenetic resolution.     

A DNA and morphology based phylogenetic framework of the ant genus <Emphasis Type="Italic">Lasius</Emphasis> with hypotheses for the evolution of social parasitism and fungiculture

Background  

Ants of the genus Lasius are ecologically important and an important system for evolutionary research. Progress in evolutionary research has been hindered by the lack of a well-founded phylogeny of the subgenera, with three previous attempts disagreeing. Here we employed two mitochondrial genes (cytochrome c oxidase subunit I, 16S ribosomal RNA), comprising 1,265 bp, together with 64 morphological characters, to recover the phylogeny of Lasius by Bayesian and Maximum Parsimony inference after exploration of potential causes of phylogenetic distortion. We use the resulting framework to infer evolutionary pathways for social parasitism and fungiculture.     

ConStruct: Improved construction of RNA consensus structures

Background  

Aligning homologous non-coding RNAs (ncRNAs) correctly in terms of sequence and structure is an unresolved problem, due to both mathematical complexity and imperfect scoring functions. High quality alignments, however, are a prerequisite for most consensus structure prediction approaches, homology searches, and tools for phylogeny inference. Automatically created ncRNA alignments often need manual corrections, yet this manual refinement is tedious and error-prone.     

The accuracy of several multiple sequence alignment programs for proteins

Background  

There have been many algorithms and software programs implemented for the inference of multiple sequence alignments of protein and DNA sequences. The "true" alignment is usually unknown due to the incomplete knowledge of the evolutionary history of the sequences, making it difficult to gauge the relative accuracy of the programs.     

Bayesian random local clocks,or one rate to rule them all

Background  

Relaxed molecular clock models allow divergence time dating and "relaxed phylogenetic" inference, in which a time tree is estimated in the face of unequal rates across lineages. We present a new method for relaxing the assumption of a strict molecular clock using Markov chain Monte Carlo to implement Bayesian modeling averaging over random local molecular clocks. The new method approaches the problem of rate variation among lineages by proposing a series of local molecular clocks, each extending over a subregion of the full phylogeny. Each branch in a phylogeny (subtending a clade) is a possible location for a change of rate from one local clock to a new one. Thus, including both the global molecular clock and the unconstrained model results, there are a total of 2^2n-2 possible rate models available for averaging with 1, 2, ..., 2n - 2 different rate categories.     

iGTP: A software package for large-scale gene tree parsimony analysis

Background  

The ever-increasing wealth of genomic sequence information provides an unprecedented opportunity for large-scale phylogenetic analysis. However, species phylogeny inference is obfuscated by incongruence among gene trees due to evolutionary events such as gene duplication and loss, incomplete lineage sorting (deep coalescence), and horizontal gene transfer. Gene tree parsimony (GTP) addresses this issue by seeking a species tree that requires the minimum number of evolutionary events to reconcile a given set of incongruent gene trees. Despite its promise, the use of gene tree parsimony has been limited by the fact that existing software is either not fast enough to tackle large data sets or is restricted in the range of evolutionary events it can handle.     

BEAST: Bayesian evolutionary analysis by sampling trees

Background  

The evolutionary analysis of molecular sequence variation is a statistical enterprise. This is reflected in the increased use of probabilistic models for phylogenetic inference, multiple sequence alignment, and molecular population genetics. Here we present BEAST: a fast, flexible software architecture for Bayesian analysis of molecular sequences related by an evolutionary tree. A large number of popular stochastic models of sequence evolution are provided and tree-based models suitable for both within- and between-species sequence data are implemented.     

Stability of multiple alignments and phylogenetic trees: an analysis of ABC-transporter proteins family

Background  

Sequence-based phylogeny reconstruction is a fundamental task in Bioinformatics. Practically all methods for phylogeny reconstruction are based on multiple alignments. The quality and stability of the underlying alignments is therefore crucial for phylogenetic analysis.     

Reconstructing phylogenies from noisy quartets in polynomial time with a high success probability

Background  

In recent years, quartet-based phylogeny reconstruction methods have received considerable attentions in the computational biology community. Traditionally, the accuracy of a phylogeny reconstruction method is measured by simulations on synthetic datasets with known "true" phylogenies, while little theoretical analysis has been done. In this paper, we present a new model-based approach to measuring the accuracy of a quartet-based phylogeny reconstruction method. Under this model, we propose three efficient algorithms to reconstruct the "true" phylogeny with a high success probability.     

Complex inter‐island colonization and peripatric founder speciation promote diversification of flightless Pachyrhynchus weevils in the Taiwan–Luzon volcanic belt

Aim

We investigated the spatial and temporal patterns of diversification among colourful and flightless weevils, the Pachyrhynchus orbifer complex, to test the stepping‐stone hypothesis of colonization across the Taiwan–Luzon volcanic belt.

Location

Southeast Asia.

Methods

The phylogeny of the P. orbifer complex was reconstructed from a multi‐locus data set of mitochondrial and nuclear genes using maximum likelihood in RAxML and Bayesian inference in MRBAYES. Likelihood‐based tests in CONSEL were used to evaluate alternative tree topologies. Divergence times were estimated in beast based on a range of mutation rates. Ancestral range and biogeographical history were reconstructed using Bayesian binary MCMC (BBM) methods in RASP and in BioGeoBEARS. Demographic histories were inferred using the extended Bayesian skyline plot (EBSP). Species boundaries were tested using BPP.

Results

The phylogeny of the P. orbifer complex indicated strong support for seven reciprocally monophyletic lineages grouped by current island boundaries (Camiguin, Fuga, Dalupiri, Calayan, Babuyan, Orchid and Yaeyama Islands), except for a sister Green + Itbayat lineage. Complex and stochastic colonization of P. orbifer was inferred to have involved both northward and southward directions with short‐ and long‐distance dispersal events, which are strongly inconsistent with the strict stepping‐stone hypothesis. Divergence time estimates for all extant island lineages (<1 Myr of Middle Pleistocene) are much more recent than the geological ages (22.4–1.7 Myr) and subaerial existence (c. 3 Myr) of the islands. The statistically delimited seven cryptic species imply that the diversity of Pachyrhynchus from small peripheral islands continues to be largely under‐estimated.

Main conclusions

The non‐linear, more complex spatial and temporal settings of the archipelago and stochastic dispersal were probable key factors shaping the colonization history of the P. orbifer complex. Speciation of the P. orbifer complex may have occurred only between islands, indicating that peripatric speciation through the founders of stochastic dispersals was the major evolutionary driver.     

Enhanced Bayesian modelling in BAPS software for learning genetic structures of populations

Background  

During the most recent decade many Bayesian statistical models and software for answering questions related to the genetic structure underlying population samples have appeared in the scientific literature. Most of these methods utilize molecular markers for the inferences, while some are also capable of handling DNA sequence data. In a number of earlier works, we have introduced an array of statistical methods for population genetic inference that are implemented in the software BAPS. However, the complexity of biological problems related to genetic structure analysis keeps increasing such that in many cases the current methods may provide either inappropriate or insufficient solutions.     

The phylogeny of the social wasp subfamily Polistinae: evidence from microsatellite flanking sequences,mitochondrial COI sequence,and morphological characters

Background  

Social wasps in the subfamily Polistinae (Hymenoptera: Vespidae) have been important in studies of the evolution of sociality, kin selection, and within colony conflicts of interest. These studies have generally been conducted within species, because a resolved phylogeny among species is lacking. We used nuclear DNA microsatellite flanking sequences, mitochondrial COI sequence, and morphological characters to generate a phylogeny for the Polistinae (Hymenoptera) using 69 species.     

A web-database of mammalian morphology and a reanalysis of placental phylogeny

Background  

Recent publications concerning the interordinal phylogeny of placental mammals have converged on a common signal, consisting of four major radiations with some ambiguity regarding the placental root. The DNA data with which these relationships have been reconstructed are easily accessible from public databases; access to morphological characters is much more difficult. Here, I present a graphical web-database of morphological characters focusing on placental mammals, in tandem with a combined-data phylogenetic analysis of placental mammal phylogeny.     

Molecular phylogeny and evolution of inflorescence types in Eperua

Premise

The Amazonian hyperdominant genus Eperua (Fabaceae) currently holds 20 described species and has two strongly different inflorescence and flower types, with corresponding different pollination syndrome. The evolution of these vastly different inflorescence types within this genus was unknown and the main topic in this study.

Methods

We constructed a molecular phylogeny, based on the full nuclear ribosomal DNA and partial plastome, using Bayesian inference and maximum likelihood methods, to test whether the genus is monophyletic, whether all species are monophyletic and if the shift from bat to bee pollination (or vice versa) occurred once in this genus.

Results

All but two species are well supported by the nuclear ribosomal phylogeny. The plastome phylogeny, however, shows a strong geographic signal suggesting strong local hybridization or chloroplast capture, rendering chloroplast barcodes meaningless in this genus.

Conclusions

With our data, we cannot fully resolve the backbone of the tree to clarify sister genera relationships and confirm monophyly of the genus Eperua. Within the genus, the shift from bat to bee and bee to bat pollination has occurred several times but, with the bee to bat not always leading to a pendant inflorescence.

     

Bayesian and maximum likelihood phylogenetic analyses of protein sequence data under relative branch-length differences and model violation

Background  

Bayesian phylogenetic inference holds promise as an alternative to maximum likelihood, particularly for large molecular-sequence data sets. We have investigated the performance of Bayesian inference with empirical and simulated protein-sequence data under conditions of relative branch-length differences and model violation.     

Nuclear gene phylogeography using PHASE: dealing with unresolved genotypes,lost alleles,and systematic bias in parameter estimation

Background  

A widely-used approach for screening nuclear DNA markers is to obtain sequence data and use bioinformatic algorithms to estimate which two alleles are present in heterozygous individuals. It is common practice to omit unresolved genotypes from downstream analyses, but the implications of this have not been investigated. We evaluated the haplotype reconstruction method implemented by PHASE in the context of phylogeographic applications. Empirical sequence datasets from five non-coding nuclear loci with gametic phase ascribed by molecular approaches were coupled with simulated datasets to investigate three key issues: (1) haplotype reconstruction error rates and the nature of inference errors, (2) dataset features and genotypic configurations that drive haplotype reconstruction uncertainty, and (3) impacts of omitting unresolved genotypes on levels of observed phylogenetic diversity and the accuracy of downstream phylogeographic analyses.