EPoS: a modular software framework for phylogenetic analysis

Estimating Phylogenies of Species (EPoS) is a modular software framework for phylogenetic analysis, visualization and data management. It provides a plugin-based system that integrates a storage facility, a rich user interface and the ability to easily incorporate new methods, functions and visualizations. EPoS ships with persistent data management, a set of well-known phylogenetic algorithms and a multitude of tree visualization methods and layouts. Implemented algorithms cover distance-based tree construction, consensus trees and various graph-based supertree methods. The rendering system can be customized for, say, different edge and node styles.     

Multilocus sequence analysis reveals multiple symbiovars within Mesorhizobium species

The genus Mesorhizobium includes species nodulating several legumes, such as chickpea, which has a high agronomic importance. Chickpea rhizobia were originally described as either Mesorhizobium ciceri or M. mediterraneum. However, rhizobia able to nodulate chickpea have been shown to belong to several different species within the genus Mesorhizobium. The present study used a multilocus sequence analysis approach to infer a high resolution phylogeny of the genus Mesorhizobium and to confirm the existence of a new chickpea nodulating genospecies. The phylogenetic structure of the Mesorhizobium clade was evaluated by sequence analysis of the 16S rRNA gene, ITS region and the five core genes atpD, dnaJ, glnA, gyrB, and recA. Phylogenies obtained with the different genes are in overall good agreement and a well-supported, almost fully resolved, phylogenetic tree was obtained using the combined data. Our phylogenetic analyses of core genes sequences and their comparison with the symbiosis gene nodC, corroborate the existence of one new chickpea Mesorhizobium genospecies and one new symbiovar, M. opportunistum sv. ciceri. Furthermore, our results show that symbiovar ciceri spreads over six species of mesorhizobia. To our knowledge this study shows the most complete Mesorhizobium multilocus phylogeny to date and contributes to the understanding of how a symbiovar may be present in different species.     

Euclidean nature of phylogenetic distance matrices

Phylogenies are fundamental to comparative biology as they help to identify independent events on which statistical tests rely. Two groups of phylogenetic comparative methods (PCMs) can be distinguished: those that take phylogenies into account by introducing explicit models of evolution and those that only consider phylogenies as a statistical constraint and aim at partitioning trait values into a phylogenetic component (phylogenetic inertia) and one or multiple specific components related to adaptive evolution. The way phylogenetic information is incorporated into the PCMs depends on the method used. For the first group of methods, phylogenies are converted into variance-covariance matrices of traits following a given model of evolution such as Brownian motion (BM). For the second group of methods, phylogenies are converted into distance matrices that are subsequently transformed into Euclidean distances to perform principal coordinate analyses. Here, we show that simply taking the elementwise square root of a distance matrix extracted from a phylogenetic tree ensures having a Euclidean distance matrix. This is true for any type of distances between species (patristic or nodal) and also for trees harboring multifurcating nodes. Moreover, we illustrate that this simple transformation using the square root imposes less geometric distortion than more complex transformations classically used in the literature such as the Cailliez method. Given the Euclidean nature of the elementwise square root of phylogenetic distance matrices, the positive semidefinitiveness of the phylogenetic variance-covariance matrix of a trait following a BM model, or related models of trait evolution, can be established. In that way, we build a bridge between the two groups of statistical methods widely used in comparative analysis. These results should be of great interest for ecologists and evolutionary biologists performing statistical analyses incorporating phylogenies.     

The evolution of myiasis in blowflies (Calliphoridae)

Blowflies (Calliphoridae) are characterised by the ability of their larvae to develop in animal flesh. Where the host is a living vertebrate, such parasitism by dipterous larvae is known as myiasis. However, the evolutionary origins of the myiasis habit in the Calliphoridae, a family which includes the blowflies and screwworm flies, remain unclear. Species associated with an ectoparasitic lifestyle can be divided generally into three groups based on their larval feeding habits: saprophagy, facultative ectoparasitism, and obligate parasitism, and it has been proposed that this functional division may reflect the progressive evolution of parasitism in the Calliphoridae. In order to evaluate this hypothesis, phylogenetic analysis of 32 blowfly species displaying a range of forms of ectoparasitism from key subfamilies, i.e. Calliphorinae, Luciliinae, Chrysomyinae, Auchmeromyiinae and Polleniinae, was undertaken using likelihood and parsimony methods. Phylogenies were constructed from the nuclear 28S large subunit ribosomal RNA gene (28S rRNA), sequenced from each of the 32 calliphorid species, together with suitable outgroup taxa, and mitochondrial cytochrome oxidase subunit I and II (COI+II) sequences, derived primarily from published data. Phylogenies derived from each of the two markers (28S rRNA, COI+II) were largely (though not completely) congruent, as determined by incongruence-length difference and Kishino-Hasegawa tests. However, the phylogenetic relationships of blowfly subfamilies based on molecular data did not concur with the pattern of relationships defined by previous morphological analysis; significantly, molecular analysis supported the monophyly of blowflies (Calliphoridae), distinct from the bot and warble flies (Oestridae). Comparative analysis of the myiasis habit based primarily on the 28S rRNA phylogeny indicated that obligate parasitism, and the ability to initiate myiasis in higher vertebrates, has multiple independent origins across myiasis-causing flies (Calliphoridae and Oestridae) and in at least three subfamilies of blowfly (Calliphoridae). Finally, the general association of various blowfly genera and subfamily clades with particular continental and geographical regions suggests that these groups probably came into existence in the Late Cretaceous period, following the break-up of Gondwana.     

Phylogenetic significance of morphological characters in the taxonomy of Pestalotiopsis species

There has been considerable disagreement regarding the relationships among Pestalotiopsis species and their delimitations. A molecular phylogenetic analysis was conducted on 32 species of Pestalotiopsis in order to evaluate the utility of morphological characters currently used in their taxonomy. Phylogenetic relationships were inferred from nucleotide sequences in the ITS regions and 5.8S gene of the rDNA under four optimality criteria: maximum parsimony, weighted parsimony, maximum likelihood, and neighbor joining. Phylogenies estimated from all analyses yielded trees of essentially similar topology and revealed 3 major groups that correspond with morphology-based classification systems. Molecular data indicated that the genus contains two distinct lineages based on pigmentation of median cells and four distinct groupings based on morphology of apical appendages. The analyses did not support reliability of other phenotypic characters of this genus, such as spore dimensions. Characters with particular phylogenetic significance are discussed in relation to the taxonomy of Pestalotiopsis.     

The Use (and Misuse) of Phylogenetic Trees in Comparative Behavioral Analyses

Phylogenetic comparative methods play a critical role in our understanding of the adaptive origin of primate behaviors. To incorporate evolutionary history directly into comparative behavioral research, behavioral ecologists rely on strong, well-resolved phylogenetic trees. Phylogenies provide the framework on which behaviors can be compared and homologies can be distinguished from similarities due to convergent or parallel evolution. Phylogenetic reconstructions are also of critical importance when inferring the ancestral state of behavioral patterns and when suggesting the evolutionary changes that behavior has undergone. Improvements in genome sequencing technologies have increased the amount of data available to researchers. Recently, several primate phylogenetic studies have used multiple loci to produce robust phylogenetic trees that include hundreds of primate species. These trees are now commonly used in comparative analyses and there is a perception that we have a complete picture of the primate tree. But how confident can we be in those phylogenies? And how reliable are comparative analyses based on such trees? Herein, we argue that even recent molecular phylogenies should be treated cautiously because they rely on many assumptions and have many shortcomings. Most phylogenetic studies do not model gene tree diversity and can produce misleading results, such as strong support for an incorrect species tree, especially in the case of rapid and recent radiations. We discuss implications that incorrect phylogenies can have for reconstructing the evolution of primate behaviors and we urge primatologists to be aware of the current limitations of phylogenetic reconstructions when applying phylogenetic comparative methods.     

Are phylogenies resolved at the genus level appropriate for studies on phylogenetic structure of species assemblages?

Phylogenies are essential to studies investigating the effect of evolutionary history on assembly of species in ecological communities and geographical and ecological patterns of phylogenetic structure of species assemblages. Because phylogenies well resolved at the species level are lacking for many major groups of organisms such as vascular plants, researchers often generate a species-level phylogenies using a phylogeny well resolved at the genus level as a backbone and attaching species to their respective genera in the phylogeny as polytomies or by using a megaphylogeny well resolved at the genus level as a backbone and adding additional species to the megaphylogeny as polytomies of their respective genera. However, whether the result of a study using species-level phylogenies generated in these ways is robust, compared to that based on phylogenies fully resolved at the species level, has not been assessed. Here, we use 1093 angiosperm tree assemblages (each in a 110 × 110 km quadrat) in North America as a model system to address this question, by examining six commonly used metrics of phylogenetic structure (phylogenetic diversity and phylogenetic relatedness) and six climate variables commonly used in ecology. Our results showed that (1) the scores of phylogenetic metrics derived from species-level phylogenies resolved at the genus level with species being attached to their respective genera as polytomies are very strongly or perfectly correlated to those derived from a phylogeny fully resolved at the species level (the mean of correlation coefficients is 0.973), and (2) the relationships between the scores of phylogenetic metrics and climate variables are consistent between the two sets of analyses based on the two types of phylogeny. Our study suggests that using species-level phylogenies resolved at the genus level with species being attached to their genera as polytomies is appropriate in studies exploring patterns of phylogenetic structure of species in ecological communities across geographical and ecological gradients.     

Phylogenetic--evolutionary approaches to bioinformatics

Phylogenies of organisms are essential to investigating a range of evolutionary questions of interest to researchers in the field of bioinformatics. Phylogenies not only help to define how to study many evolutionary questions, they must also be taken into account when conducting statistical analyses. Here it is shown how phylogenies can be used to investigate variability along the sites of a gene, reconstruct ancestral states of ancient genes and proteins, identify and characterise events of parallel and convergent evolution, find events of gene duplication, analyse predictions from molecular clocks, seek evidence for correlated changes among different parts of the same gene or genome, and test theories of molecular evolution. A table of statistical and phylogenetic methods is presented.     

Highly incomplete taxa can rescue phylogenetic analyses from the negative impacts of limited taxon sampling

Background

Phylogenies are essential to many areas of biology, but phylogenetic methods may give incorrect estimates under some conditions. A potentially common scenario of this type is when few taxa are sampled and terminal branches for the sampled taxa are relatively long. However, the best solution in such cases (i.e., sampling more taxa versus more characters) has been highly controversial. A widespread assumption in this debate is that added taxa must be complete (no missing data) in order to save analyses from the negative impacts of limited taxon sampling. Here, we evaluate whether incomplete taxa can also rescue analyses under these conditions (empirically testing predictions from an earlier simulation study).

Methodology/Principal Findings

We utilize DNA sequence data from 16 vertebrate species with well-established phylogenetic relationships. In each replicate, we randomly sample 4 species, estimate their phylogeny (using Bayesian, likelihood, and parsimony methods), and then evaluate whether adding in the remaining 12 species (which have 50, 75, or 90% of their data replaced with missing data cells) can improve phylogenetic accuracy relative to analyzing the 4 complete taxa alone. We find that in those cases where sampling few taxa yields an incorrect estimate, adding taxa with 50% or 75% missing data can frequently (>75% of relevant replicates) rescue Bayesian and likelihood analyses, recovering accurate phylogenies for the original 4 taxa. Even taxa with 90% missing data can sometimes be beneficial.

Conclusions

We show that adding taxa that are highly incomplete can improve phylogenetic accuracy in cases where analyses are misled by limited taxon sampling. These surprising empirical results confirm those from simulations, and show that the benefits of adding taxa may be obtained with unexpectedly small amounts of data. These findings have important implications for the debate on sampling taxa versus characters, and for studies attempting to resolve difficult phylogenetic problems.     

Phylogenetic selection of Narcissus species for drug discovery

Plants are of immense importance in providing healthcare worldwide. With over 250 000 species of angiosperms alone, the potential for finding new medicinal plants and lead compounds for drug development is enormous. Some way of selecting plants for drug discovery programs is necessary. Phylogenies have great explanatory power and also enable a predictive perspective not offered by previous classifications of plants. Phylogenetic selection of target species is a new approach to drug discovery and the present study is the first attempt to correlate acetylcholine esterase (AChE) inhibitory activity and alkaloid distribution with a molecular phylogenetic hypothesis of Narcissus. The distribution of alkaloids with AChE inhibitory activity is significantly constrained by the phylogeny. Simultaneous evaluation of all available information of alkaloids and AChE inhibitory activity in a phylogenetic framework allowed us to discuss various strategies for selection of target species for further studies of AChE inhibitory activity.     

Missing something? Codon aversion as a new character system in phylogenetics

Although many studies have documented codon usage bias in different species, the importance of codon usage in a phylogenetic framework remains largely unknown. We demonstrate that a phylogenetic signal is present in the codon usage and non‐usage biases of 17 717 orthologues evaluated across 72 tetrapod species using a simple parsimony analysis of a binary matrix of codon characters. Phylogenies estimated using stop codons were more congruent with previous hypotheses than phylogenies based on any other single codon or a combination of codons. Although each codon is present in every species, specific genes have different codon preferences and may or may not use every possible codon. This observation allowed us to map the pattern of codon usage and non‐usage across the topology. These results suggest that codon usage is phylogenetically conserved across shallow and deep levels within tetrapods.     

A new dynamic null model for phylogenetic community structure

Phylogenies are increasingly applied to identify the mechanisms structuring ecological communities but progress has been hindered by a reliance on statistical null models that ignore the historical process of community assembly. Here, we address this, and develop a dynamic null model of assembly by allopatric speciation, colonisation and local extinction. Incorporating these processes fundamentally alters the structure of communities expected due to chance, with speciation leading to phylogenetic overdispersion compared to a classical statistical null model assuming equal probabilities of community membership. Applying this method to bird and primate communities in South America we show that patterns of phylogenetic overdispersion – often attributed to negative biotic interactions – are instead consistent with a species neutral model of allopatric speciation, colonisation and local extinction. Our findings provide a new null expectation for phylogenetic community patterns and highlight the importance of explicitly accounting for the dynamic history of assembly when testing the mechanisms governing community structure.     

Phylogenetic relationships and biogeography of xantusiid lizards, inferred from mitochondrial DNA sequences.

Portions of two mitochondrial genes (12S ribosomal RNA and cytochrome b) were sequenced in seven species to examine phylogenetic relationships within the lizard family Xantusiidae. Phylogenies derived from these sequences (709 total bp) are concordant and indicate that the Cuban species Cricosaura typica is the sister group to all other xantusiids. The Middle American genus Lepidophyma is the closest relative of Xantusia, and X. riversiana (California Islands) the closest relative of X. vigilis (mainland). These findings are not in agreement either with the results of a recent morphological analysis that united Cricosaura and Lepidophyma as closest relatives or with past studies that have recognized X. riversiana as a separate genus. Levels of sequence divergence, as well as the age and affinities of some mainland fossil taxa, suggest that the origin of Cricosaura was associated with the tectonic evolution of the Greater Antilles in the late Cretaceous. These results further demonstrate that significant resolution of phylogenies can be obtained with relatively short DNA sequences and that these mitochondrial genes are concordant in their estimation of phylogeny.     

Early contributions of molecular phylogenetics to understanding the evolution of Rotifera

The past decade has seen the application of DNA sequence data to phylogenetic investigations of Rotifera, both expanding and challenging our understanding of the evolution of the phylum. Evidence that Acanthocephala, long regarded as a separate but closely related phylum, is a highly derived class of Rotifera demonstrates the potential of molecular analyses to suggest relationships not obvious from morphological analysis. Phylogenies based on the sequence of the gene for the small ribosomal RNA suggest that rotifers and acanthocephalans are associated with Platyhelminthes and Gastrotricha, perhaps in a clade with Gnathostomula and Cycliophora; at present, this group lacks a clear morphological synapomorphy. A more complete resolution of the molecular phylogeny of Rotifera will require surveying multiple genes and several species from each clade under investigation.     

Molecular phylogenetics of the moth genus Omiodes Guenée (Crambidae: Spilomelinae), and the origins of the Hawaiian lineage

The moth genus Omiodes (Crambidae) comprises about 80 species and has a circumtropical distribution, with the type species, O. humeralis, occurring in Central America. In Hawaii, there are 23 native species currently placed in Omiodes, but this classification has been disputed, and they were previously placed in various other genera. We used molecular phylogenetic analyses to assess the monophyly of Omiodes as a whole, and specifically of the Hawaiian species, as well as their geographic origins and possible ancestral host plants. Mitochondrial (COI) and nuclear (wingless, EF1α, CAD, and RPS5) DNA was sequenced for Omiodes from Hawaii, South America, and Australasia, along with many other putative outgroup spilomeline genera. Phylogenies were estimated using maximum likelihood and Bayesian inference, and various taxon and character datasets. With the exception of two paleotropical species (O. basalticalis and O. odontosticta, whose placement was unresolved) all Hawaiian, paleotropical and neotropical Omiodes, including the type species, fell within a well-supported, monophyletic clade. Although the center of diversity for Omiodes is in the Neotropics, its center of origin was ambiguous, due to poor resolution of the basal splits between paleotropical and neotropical Omiodes. Very low genetic divergence within the Hawaiian Omiodes suggests a relatively recent colonization of the Hawaiian Islands. Phylogenies constructed using all codon positions were poorly resolved at intergeneric levels, and did not reveal a sister taxon for Omiodes, but phylogenies constructed using only first and second codon positions suggested a close relationship with Cnaphalocrocis. The monophyly of several other spilomeline genera is also discussed.     

Nuclear and mtDNA phylogenies of the Trimeresurus complex: implications for the gene versus species tree debate

Phylogenies based on mitochondrial DNA (mtDNA) may represent gene trees that may not be congruent with the equivalent species tree. One solution to this problem is to include additional, independent loci from the nuclear genome. Sequence data from the seventh intron of the beta-fibrinogen gene were generated for 25 specimens of vipers, including 8 nominal species of the Trimeresurus complex of Asian pit vipers. Phylogenetic trees were generated using maximum-parsimony and maximum-likelihood methods. The taxonomic level at which the intron provided significant phylogenetic information was examined and the trees were compared to those produced from previously obtained mtDNA cytochrome b sequences. A variety of different approaches (separate analyses, conditional data combination, and consensus) were used in an attempt to provide a sound organismal phylogeny based on both nuclear and mtDNA data sets. We discuss the implications for the gene tree-species tree debate and its particular relevance to medically important organisms.     

Phylogeny and cryptic diversity in geckos (Phyllopezus; Phyllodactylidae; Gekkota) from South America's open biomes

The gecko genus Phyllopezus occurs across South America's open biomes: Cerrado, Seasonally Dry Tropical Forests (SDTF, including Caatinga), and Chaco. We generated a multi-gene dataset and estimated phylogenetic relationships among described Phyllopezus taxa and related species. We included exemplars from both described Phyllopezus pollicaris subspecies, P. p. pollicaris and P. p.przewalskii. Phylogenies from the concatenated data as well as species trees constructed from individual gene trees were largely congruent. All phylogeny reconstruction methods showed Bogertia lutzae as the sister species of Phyllopezus maranjonensis, rendering Phyllopezus paraphyletic. We synonymized the monotypic genus Bogertia with Phyllopezus to maintain a taxonomy that is isomorphic with phylogenetic history. We recovered multiple, deeply divergent, cryptic lineages within P. pollicaris. These cryptic lineages possessed mtDNA distances equivalent to distances among other gekkotan sister taxa. Described P. pollicaris subspecies are not reciprocally monophyletic and current subspecific taxonomy does not accurately reflect evolutionary relationships among cryptic lineages. We highlight the conservation significance of these results in light of the ongoing habitat loss in South America's open biomes.     

Phylogenetics of the Caprinae Based on Cytochrome b Sequence

Relationships within the subfamily Caprinae have never been fully resolved. Phylogenies have been proposed based on morphological, behavioral, ecological, and some molecular comparisons. Because of the relatively poor fossil record of the Caprinae, paleontological evidence has not been extensively used in phylogenetic reconstruction for this group. Traditionally, four tribes: Saigini, Ovibovini, Rupicaprini, and Caprini, have been recognized. We investigated relationships within the Caprinae by comparing sequences of the cytochrome b gene of mitochondrial DNA from 11 species of Caprinae and 1 Bovinae species. Our analyses suggest that revisions to previous phylogenies, including the dissolution of the Ovibovini, are warranted.     

Differential rates of genic and chromosomal evolution in bats of the family Rhinolophidae

Data for nondifferentially stained chromosomes from 10 species of Rhinolophus (Chiroptera: Rhinolophidae) suggest a conserved chromosomal evolution. G-banded chromosomes for three well differentiated species (Rhinolophus hipposideros, Rhinolophus blasii, and Rhinolophus acuminatus) corroborate a low level of gross chromosomal rearrangements. Additionally, a comparison between G-banded chromosomes of Rhinolophus (Rhinolophidae) and Hipposideros (Hipposideridae) suggests extreme conservatism in chromosomal arms between these two distantly related groups. On the other hand, we report extensive genic divergence as assayed by starch gel electrophoresis among these 10 species, and between Rhinolophus and two hipposiderid genera (Hipposideros and Aselliscus). The present chromosomal data are not sufficient for phylogenetic analysis. Phylogenies based on electrophoretic data are in many aspects discordant with those based on the classical morphological criteria. Different (and as yet not clearly understood) evolutionary forces affecting chromosomal, morphologic, and electrophoretic variation may be the reason for the apparent lack of concordance in these independent data sets.     

Speciation and Species Separation in Hordeum L. (Poaceae) Resolved by Discontinuous Molecular Markers1

Abstract: Amplified fragment length polymorphisms (AFLPs) were used to evaluate the capacity of discontinuous markers to reveal genetic structure within Hordeum , a challenging higher plant genus from the standpoint of natural systematics. Phylogenies of 63 accessions encompassing nine species from four Hordeum sections were inferred from polymorphisms scored at 600 loci. Phylogenies based on sequences from the nuclear internal transcribed spacer (ITS) region were constructed for comparison, but revealed severe sampling errors inherent to single genes. Although superior by virtue of providing genome-wide estimates of genetic similarity, the adoption of AFLPs in infrageneric studies requires caution. Comigrating AFLP bands studied here could be divided on the basis of band intensity variation into two types that are □ 100 % identical and < 40 % identical in DNA sequence, respectively, in infrageneric comparisons. Thus, the careful selection of AFLP bands to be analyzed bears heavily upon their phylogenetic utility. Within the H. murinum complex, which encompasses three morphologically distinct subspecies, AFLP data from 37 accessions reveal unexpected genetic differentiation between H. murinum, glaucum populations to the east and west of Alexandria (Egypt), suggesting the presence of allopatric speciation in the wake of human settlement.