Inferring phylogenies of evolving sequences without multiple sequence alignment

Alignment-free methods, in which shared properties of sub-sequences (e.g. identity or match length) are extracted and used to compute a distance matrix, have recently been explored for phylogenetic inference. However, the scalability and robustness of these methods to key evolutionary processes remain to be investigated. Here, using simulated sequence sets of various sizes in both nucleotides and amino acids, we systematically assess the accuracy of phylogenetic inference using an alignment-free approach, based on D2 statistics, under different evolutionary scenarios. We find that compared to a multiple sequence alignment approach, D2 methods are more robust against among-site rate heterogeneity, compositional biases, genetic rearrangements and insertions/deletions, but are more sensitive to recent sequence divergence and sequence truncation. Across diverse empirical datasets, the alignment-free methods perform well for sequences sharing low divergence, at greater computation speed. Our findings provide strong evidence for the scalability and the potential use of alignment-free methods in large-scale phylogenomics.

Evidence of host specificity and congruence between phylogenies of bitterling and freshwater mussels

Evidence of host specificity and congruence between phylogenies of bitterling and freshwater mussels. Zoological Studies 45(3): 428-434. Bitterling (Cyprinidae: Acheilognathinae) are freshwater fishes with a unique spawning relationship with freshwater mussels on whose gills they lay their eggs. During the breeding season of bitterling fishes, we collected 843 mussels belonging to 16 species from Lake Qinglan, central China and examined their gill chambers for the presence of bitterling larvae. Three species of bitterling larvae were identified; Acheilognathus tonkinensis, Ach. cf. meridianus, and Ach. barbatulus, in 3 species of mussel: Unio douglasiae, Lamprotula caveata, and L. tortuosa, suggesting host specialization. Using our own and other published data, we compared the respective phylogenies of bitterling and mussels, but failed to show clear congruence. However, broad specializations are evident, with Acheilognathus and Tanakia showing preferences for mussels with a relatively simple gill structure (Ableminae), and Rhodeus spp. showing preferences for mussels of the Anodontinae and Unioninae, which have more-complex gill structures.

Four independent gene phylogenies confirm the ancient relationships of Madagascar endemic species in the Papilio demoleus group (Lepidoptera, Papilionidae).

Madagascar is home to numerous endemic species and lineages, but the processes that have contributed to its endangered diversity are still poorly understood. Evidence is accumulating to demonstrate the importance of Tertiary dispersal across varying distances of oceanic barriers, supplementing vicariance relationships dating back to the Cretaceous, but these hypotheses remain tentative in the absence of well-supported phylogenies. In the Papilio demoleus group of swallowtail butterflies, three of the five recognized species are restricted to Madagascar, whereas the remaining two species range across the Afrotropical zone and southern Asia plus Australia. We reconstructed phylogenetic relationships for all species in the P. demoleus group, as well as 11 outgroup Papilio species, using 60 morphological characters and about 4 kb of nucleotide sequences from two mitochondrial (cytochrome oxidase I and II) and two nuclear (wg and EF-1a) genes. Of the three endemic Malagasy species, the two that are formally listed as endangered or at risk represented the most basal divergences in the group, while the more common third endemic was clearly related to African P. demodocus. The fifth species, P. demoleus, showed little differentiation across southern Asia, but showed divergence from its subspecies sthenelus in Australia. Dispersal-vicariance analysis using cladograms derived from morphology and three independent genes indicated a Malagasy diversification of lime swallowtails in the middle Miocene. Thus, diversification processes on the island of Madagascar may have contributed to the origin of common butterflies that now occur throughout much of the Old World tropical and subtemperate regions. An alternative hypothesis, that Madagascar is a refuge for ancient lineages resulting from successive colonizations from Africa, is less parsimonious and does not explain the relatively low continental diversity of the group.

Molecular Phylogenies Support Homoplasy of Multiple Morphological Characters Used in the Taxonomy of Heteroscleromorpha (Porifera: Demospongiae)

Sponge classification has long been based mainly on morphocladistic analyses but is now being greatly challenged by more than 12 years of accumulated analyses of molecular data analyses. The current study used phylogenetic hypotheses based on sequence data from 18S rRNA, 28S rRNA, and the CO1 barcoding fragment, combined with morphology to justify the resurrection of the order Axinellida Lévi, 1953. Axinellida occupies a key position in different morphologically derived topologies. The abandonment of Axinellida and the establishment of Halichondrida Vosmaer, 1887 sensu lato to contain Halichondriidae Gray, 1867, Axinellidae Carter, 1875, Bubaridae Topsent, 1894, Heteroxyidae Dendy, 1905, and a new family Dictyonellidae van Soest et al., 1990 was based on the conclusion that an axially condensed skeleton evolved independently in separate lineages in preference to the less parsimonious assumption that asters (star-shaped spicules), acanthostyles (club-shaped spicules with spines), and sigmata (C-shaped spicules) each evolved more than once. Our new molecular trees are congruent and contrast with the earlier, morphologically based, trees. The results show that axially condensed skeletons, asters, acanthostyles, and sigmata are all homoplasious characters. The unrecognized homoplasious nature of these characters explains much of the incongruence between molecular-based and morphology-based phylogenies. We use the molecular trees presented here as a basis for re-interpreting the morphological characters within Heteroscleromorpha. The implications for the classification of Heteroscleromorpha are discussed and a new order Biemnida ord. nov. is erected.

Consensus and Confusion in Molluscan Trees: Evaluating Morphological and Molecular Phylogenies

Mollusks are the most morphologically disparate living animal phylum, they have diversified into all habitats, and have a deep fossil record. Monophyly and identity of their eight living classes is undisputed, but relationships between these groups and patterns of their early radiation have remained elusive. Arguments about traditional morphological phylogeny focus on a small number of topological concepts but often without regard to proximity of the individual classes. In contrast, molecular studies have proposed a number of radically different, inherently contradictory, and controversial sister relationships. Here, we assembled a dataset of 42 unique published trees describing molluscan interrelationships. We used these data to ask several questions about the state of resolution of molluscan phylogeny compared to a null model of the variation possible in random trees constructed from a monophyletic assemblage of eight terminals. Although 27 different unique trees have been proposed from morphological inference, the majority of these are not statistically different from each other. Within the available molecular topologies, only four studies to date have included the deep-sea class Monoplacophora; but 36.4% of all trees are not significantly different. We also present supertrees derived from 2 data partitions and 3 methods, including all available molecular molluscan phylogenies, which will form the basis for future hypothesis testing. The supertrees presented here were not constructed to provide yet another hypothesis of molluscan relationships, but rather to algorithmically evaluate the relationships present in the disparate published topologies. Based on the totality of available evidence, certain patterns of relatedness among constituent taxa become clear. The internodal distance is consistently short between a few taxon pairs, particularly supporting the relatedness of Monoplacophora and the chitons, Polyplacophora. Other taxon pairs are rarely or never found in close proximity, such as the vermiform Caudofoveata and Bivalvia. Our results have specific utility for guiding constructive research planning in order to better test relationships in Mollusca as well as other problematic groups. Taxa with consistently proximate relationships should be the focus of a combined approach in a concerted assessment of potential genetic and anatomical homology, while unequivocally distant taxa will make the most constructive choices for exemplar selection in higher-level phylogenomic analyses.

The evolution of mobile DNAs: when will transposons create phylogenies that look as if there is a master gene?

Some families of mammalian interspersed repetitive DNA, such as the Alu SINE sequence, appear to have evolved by the serial replacement of one active sequence with another, consistent with there being a single source of transposition: the "master gene." Alternative models, in which multiple source sequences are simultaneously active, have been called "transposon models." Transposon models differ in the proportion of elements that are active and in whether inactivation occurs at the moment of transposition or later. Here we examine the predictions of various types of transposon model regarding the patterns of sequence variation expected at an equilibrium between transposition, inactivation, and deletion. Under the master gene model, all bifurcations in the true tree of elements occur in a single lineage. We show that this property will also hold approximately for transposon models in which most elements are inactive and where at least some of the inactivation events occur after transposition. Such tree shapes are therefore not conclusive evidence for a single source of transposition.

Bayesian estimation of ancestral character states on phylogenies

Biologists frequently attempt to infer the character states at ancestral nodes of a phylogeny from the distribution of traits observed in contemporary organisms. Because phylogenies are normally inferences from data, it is desirable to account for the uncertainty in estimates of the tree and its branch lengths when making inferences about ancestral states or other comparative parameters. Here we present a general Bayesian approach for testing comparative hypotheses across statistically justified samples of phylogenies, focusing on the specific issue of reconstructing ancestral states. The method uses Markov chain Monte Carlo techniques for sampling phylogenetic trees and for investigating the parameters of a statistical model of trait evolution. We describe how to combine information about the uncertainty of the phylogeny with uncertainty in the estimate of the ancestral state. Our approach does not constrain the sample of trees only to those that contain the ancestral node or nodes of interest, and we show how to reconstruct ancestral states of uncertain nodes using a most-recent-common-ancestor approach. We illustrate the methods with data on ribonuclease evolution in the Artiodactyla. Software implementing the methods ( BayesMultiState) is available from the authors.

Phylogenies reveal new interpretation of speciation and the Red Queen

The Red Queen metaphor has species accumulating small changes to keep up with a continually changing environment, with speciation occurring at a constant rate. This constant-rate claim is now tested against four competing models, using 101 phylogenies of animal, plant and fungal taxa. The results provide a new interpretation of the Red Queen; a view linking speciation to rare stochastic events that cause reproductive isolation.

Phylogenies reveal predictive power of traditional medicine in bioprospecting

There is controversy about whether traditional medicine can guide drug discovery, and investment in ethnobotanically led research has fluctuated. One view is that traditionally used plants are not necessarily efficacious and there are no robust methods for distinguishing the ones that are most likely to be bioactive when selecting species for further testing. Here, we reconstruct a genus-level molecular phylogeny representing the 20,000 species found in the floras of three disparate biodiversity hotspots: Nepal, New Zealand and the Cape of South Africa. Borrowing phylogenetic methods from community ecology, we reveal significant clustering of the 1,500 traditionally used species, and provide a direct measure of the relatedness of the three medicinal floras. We demonstrate shared phylogenetic patterns across the floras: related plants from these regions are used to treat medical conditions in the same therapeutic areas. This strongly suggests independent discovery of plant efficacy, an interpretation corroborated by the presence of a significantly greater proportion of known bioactive species in these plant groups than in a random sample. Phylogenetic cross-cultural comparison can focus screening efforts on a subset of traditionally used plants that are richer in bioactive compounds, and could revitalise the use of traditional knowledge in bioprospecting.

Positive phenotypic selection inferred from phylogenies

Rates of phenotypic evolution vary widely in nature and these rates may often reflect the intensity of natural selection. Here we outline an approach for detecting exceptional shifts in the rate of phenotypic evolution across phylogenies. We introduce a simple new branch-specific metric ∆V/∆B that divides observed phenotypic change along a branch into two components: (1) that attributable to the background rate (∆B), and (2) that attributable to departures from the background rate (∆V). Where the amount of expected change derived from variation in the rate of morphological evolution doubles that explained by to the background rate (∆V/∆B > 2), we identify this as positive phenotypic selection. We apply our approach to six datasets, finding multiple instances of positive selection in each. Our results support the growing appreciation that the traditional gradual view of phenotypic evolution is rarely upheld, with a more episodic view taking its place. This moves focus away from viewing phenotypic evolution as a simple homogeneous process and facilitates reconciliation with macroevolutionary interpretations from a genetic perspective, paving the way to novel insights into the link between genotype and phenotype. The ability to detect positive selection when genetic data are unavailable or unobtainable represents an attractive prospect for extant species, but when applied to fossil data it can reveal patterns of natural selection in deep time that would otherwise be impossible.

Molecular phylogenies, chromosomes and dispersion in Brazilian akodontines (Rodentia, Sigmodontinae)

Uma nova filogenia molecular para roedores akodontinos do Brasil é proposta. A árvore filogenética foi enriquecida com a área de ocorrência e com informações sobre o cariótipo das amostras. Baseado nisto, e com a metodologia descrita, foram propostas hipóteses sobre as características do cariótipo e sobre a área de ocorrência dos ancestrais de cada clado. Assim, foi possível discutir hipóteses sobre evolução cromossômica do grupo, e sobre eventos de dispersão a partir da área de diferenciação original dos akodontinos nos Andes. A evolução cromossômica começou com números diplóides altos (2n=52) e mostrou uma tendência a redução (até 2n=14 em clados mais recentes). Ramos independentes da árvore mostraram redução do 2n e num caso aumentou o numero diplóide. Foram propostos pelo menos quatro eventos de dispersão dos Andes até o Brasil Sul-Oriental. Os resultados indicam a direção de novos estudos em cariologia comparada.

Comparative phylogenies of the housekeeping genes atpD, infB and rpoB and the 16S rRNA gene within the Pasteurellaceae

Phylogenies of housekeeping gene and 16S rRNA gene sequences were compared to improve the classification of the bacterial family Pasteurellaceae and knowledge of the evolutionary relationships of its members. Deduced partial protein sequences of the housekeeping genes atpD, infB and rpoB were compared in 28, 36 and 28 representative taxa of the Pasteurellaceae, respectively. The monophyly of representatives of the genus Gallibacterium was recognized by analysis of all housekeeping genes, while members of Mannheimia, Actinobacillus sensu stricto and the core group of Pasteurella sensu stricto formed monophyletic groups with two out of three housekeeping genes. Representatives of Mannheimia, Actinobacillus sensu stricto, [Haemophilus] ducreyi and [Pasteurella] trehalosi formed a monophyletic unit by analysis of all three housekeeping genes, which was in contrast to the 16S rRNA gene-derived phylogeny, where these taxa occurred at separate positions in the phylogenetic tree. Representatives of the Rodent, Avian and Aphrophilus-Haemophilus 16S rRNA gene groups were weakly supported by phylogenetic analysis of housekeeping genes. Phylogenies derived by comparison of the housekeeping genes diverged significantly from the 16S rRNA gene-derived phylogeny as evaluated by the likelihood ratio test. A low degree of congruence was also observed between the individual housekeeping gene-derived phylogenies. Estimates on speciation derived from 16S rRNA and housekeeping gene sequence comparisons resulted in quite different evolutionary scenarios for members of the Pasteurellaceae. The phylogeny based on the housekeeping genes supported observed host associations between Mannheimia, Actinobacillus sensu stricto and [Pasteurella] trehalosi and animals with paired hooves.

Using plant functional traits and phylogenies to understand patterns of plant community assembly in a seasonal tropical forest in Lao PDR

Plant functional traits reflect different evolutionary responses to environmental variation, and among extant species determine the outcomes of interactions between plants and their environment, including other plant species. Thus, combining phylogenetic and trait-based information can be a powerful approach for understanding community assembly processes across a range of spatial scales. We used this approach to investigate tree community composition at Phou Khao Khouay National Park (18°14’-18°32’N; 102°38’- 102°59’E), Laos, where several distinct forest types occur in close proximity. The aim of our study was to examine patterns of plant community assembly across the strong environmental gradients evident at our site. We hypothesized that differences in tree community composition were being driven by an underlying gradient in soil conditions. Thus, we predicted that environmental filtering would predominate at the site and that the filtering would be strongest on sandier soil with low pH, as these are the conditions least favorable to plant growth. We surveyed eleven 0.25 ha (50x50 m) plots for all trees above 10 cm dbh (1221 individual trees, including 47 families, 70 genera and 123 species) and sampled soils in each plot. For each species in the community, we measured 11 commonly studied plant functional traits covering both the leaf and wood economic spectrum traits and we reconstructed a phylogenetic tree for 115 of the species in the community using rbcL and matK sequences downloaded from Genebank (other species were not available). Finally we compared the distribution of trait values and species at two scales (among plots and 10x10m subplots) to examine trait and phylogenetic community structures. Although there was strong evidence that an underlying soil gradient was determining patterns of species composition at the site, our results did not support the hypothesis that the environmental filtering dominated community assembly processes. For the measured plant functional traits there was no consistent pattern of trait dispersion across the site, either when traits were considered individually or when combined in a multivariate analysis. However, there was a significant correlation between the degree of phylogenetic dispersion and the first principle component axis (PCA1) for the soil parameters.Moreover, the more phylogenetically clustered plots were on sandier soils with lower pH. Hence, we suggest that the community assembly processes across our sitemay reflect the influence ofmore conserved traits that we did not measure. Nevertheless, our results are equivocal and other interpretations are possible. Our study illustrates some difficulties in combining trait and phylogenetic approaches that may result from the complexities of integrating spatial and evolutionary processes that vary at different scales.