Generic delimitation and phylogenetic uncertainty: An example from a group that has undergone an explosive radiation

Phylogenetic trees can provide a stable basis for a higher-level classification of organisms that reflects evolutionary relationships. However, some lineages have a complex evolutionary history that involves explosive radiation or hybridisation. Such histories have become increasingly apparent with the use of DNA sequence data for phylogeny estimation and explain, in part, past difficulties in producing stable morphology-based classifications for some groups. We illustrate this situation by using the example of tribe Mirbelieae (Fabaceae), whose generic classification has been fraught for decades. In particular, we discuss a recent proposal to combine 19 of the 25 Mirbelieae genera into a single genus, Pultenaea sens. lat., and how we might find stable and consistent ways to squeeze something as complex as life into little boxes for our own convenience. © CSIRO.

Conceptual framework and pilot study to benchmark phylogenomic databases based on reference gene trees.

Phylogenomic databases provide orthology predictions for species with fully sequenced genomes. Although the goal seems well-defined, the content of these databases differs greatly. Seven ortholog databases (Ensembl Compara, eggNOG, HOGENOM, InParanoid, OMA, OrthoDB, Panther) were compared on the basis of reference trees. For three well-conserved protein families, we observed a generally high specificity of orthology assignments for these databases. We show that differences in the completeness of predicted gene relationships and in the phylogenetic information are, for the great majority, not due to the methods used, but to differences in the underlying database concepts. According to our metrics, none of the databases provides a fully correct and comprehensive protein classification. Our results provide a framework for meaningful and systematic comparisons of phylogenomic databases. In the future, a sustainable set of 'Gold standard' phylogenetic trees could provide a robust method for phylogenomic databases to assess their current quality status, measure changes following new database releases and diagnose improvements subsequent to an upgrade of the analysis procedure.

A phylogenetic study of the phosphoenolpyruvate carboxylase multigene family in Poaceae: understanding the molecular changes linked to C4 photosynthesis evolution

Phylogenetic reconstructions have supported several independent appearances of C₄ photosynthesis within grasses (Poaceae). These recurrent appearances appear to contradict the large number of biochemical and morphological changes required to change from C₃ to C₄, a paradox that leads to questions about the genetic changes underlying C₄ evolution. In this study, we analysed sequences encoding phosphoenolpyruvate carboxylases (PEPCs) in grasses in order to gain insights into the origin of the ppc-C₄ gene, which encodes a key enzyme in the C₄ pathway. We screened databanks for PEPC genes or cDNAs in grasses. A coding sequence of 1130 base pairs was used to build phylogenetic trees that supported the existence of four distinct PEPC gene lineages. Ppc-C₄ present in all C₄ grasses was also found in two C₃ species. The ppc-C₄ clade was congruent with the species tree, suggesting orthologous evolution. This result would imply that ppc-C₄ appeared without any duplication event. Nevertheless, caution is needed since the sampling of our study is still far from comprehensive. Further investigation with an increased sampling is recommended to elucidate the evolutionary changes underlying ppc-C₄ gene evolution in grasses.

Phylogenetic relationships and divergence times among dormice (Rodentia, Gliridae) based on three nuclear genes

We examined phylogenetic relationships among six species representing three subfamilies, Glirinae, Graphiurinae and Leithiinae with sequences from three nuclear protein-coding genes (apolipoprotein B, APOB; interphotoreceptor retinoid-binding protein, IRBP; recombination-activating gene 1, RAG1). Phylogenetic trees reconstructed from maximum-parsimony (MP), maximum-likelihood (ML) and Bayesian-inference (BI) analyses showed the monophyly of Glirinae (Glis and Glirulus) and Leithiinae (Dryomys, Eliomys and Muscardinus) with strong support, although the branch length maintaining this relationship was very short, implying rapid diversification among the three subfamilies. Divergence time estimates were calculated from ML (local clock model) and Bayesian-dating method using a calibration point of 25 Myr (million years) ago for the divergence between Glis and Glirulus, and 55 Myr ago for the split between lineages of Gliridae and Sciuridae on the basis of fossil records. The results showed that each lineage of Graphiuros, Glis, Glirulus and Muscardinus dates from the Late Oligocene to the Early Miocene period, which is mostly in agreement with fossil records. Taking into account that warm climate harbouring a glirid-favoured forest dominated from Europe to Asia during this period, it is considered that this warm environment triggered the prosperity of the glirid species through the rapid diversification. Glirulus japonicas is suggested to be a relict of this ancient diversification during the warm period.

Phylogenetic analyses of Cladophora vagabunda (L.) C. Hoek (Cladophorales, Chlorophyta) from Brazil based on SSU rDNA sequences

The complete SSU rDNA was sequenced for 10 individuals of Cladophora vagabunda collected along the coast of Brazil. For C. rupestris (L.) Kütz. a partial SSU rDNA sequence (1634 bp) was obtained. Phylogenetic trees indicate that Cladophora is paraphyletic, but the section Glomeratae sensu lato including C. vagabunda from Brazil, Japan and France, C. albida (Nees) Kütz., C. sericea (Hudson) Kütz., and C. glomerata (L.) Kütz. is monophyletic. Within this group C. vagabunda is paraphyletic. The sequence identity for the SSU rDNA varied from 98.9% to 100% for the Brazilian C. vagabunda, and from 98.3% to 99.7% comparing the Brazilian individuals to the ones from France and Japan. Sequence identity of the Brazilian C. vagabunda to C. albida and C. sericea vary from 98.0% to 98.6%. The SSU rDNA phylogeny support partially the morphological characteristics presented by Brazilian populations of C. vagabunda. On the other hand, C. rupestris from Brazil does not group with C. rupestris from France, both sequences presenting only 96.9% of identity. The inclusion of sequences of individuals from Brazil reinforces the need of taxonomical revision for the genus Cladophora and for the complex C. vagabunda.

Modelling heterotachy in phylogenetic inference by reversible-jump Markov chain Monte Carlo

The rate at which a given site in a gene sequence alignment evolves over time may vary. This phenomenon-known as heterotachy-can bias or distort phylogenetic trees inferred from models of sequence evolution that assume rates of evolution are constant. Here, we describe a phylogenetic mixture model designed to accommodate heterotachy. The method sums the likelihood of the data at each site over more than one set of branch lengths on the same tree topology. A branch-length set that is best for one site may differ from the branch-length set that is best for some other site, thereby allowing different sites to have different rates of change throughout the tree. Because rate variation may not be present in all branches, we use a reversible-jump Markov chain Monte Carlo algorithm to identify those branches in which reliable amounts of heterotachy occur. We implement the method in combination with our 'pattern-heterogeneity' mixture model, applying it to simulated data and five published datasets. We find that complex evolutionary signals of heterotachy are routinely present over and above variation in the rate or pattern of evolution across sites, that the reversible-jump method requires far fewer parameters than conventional mixture models to describe it, and serves to identify the regions of the tree in which heterotachy is most pronounced. The reversible-jump procedure also removes the need for a posteriori tests of 'significance' such as the Akaike or Bayesian information criterion tests, or Bayes factors. Heterotachy has important consequences for the correct reconstruction of phylogenies as well as for tests of hypotheses that rely on accurate branch-length information. These include molecular clocks, analyses of tempo and mode of evolution, comparative studies and ancestral state reconstruction. The model is available from the authors' website, and can be used for the analysis of both nucleotide and morphological data.

Predicting functional gene links from phylogenetic-statistical analyses of whole genomes

An important element of the developing field of proteomics is to understand protein-protein interactions and other functional links amongst genes. Across-species correlation methods for detecting functional links work on the premise that functionally linked proteins will tend to show a common pattern of presence and absence across a range of genomes. We describe a maximum likelihood statistical model for predicting functional gene linkages. The method detects independent instances of the correlated gain or loss of pairs of proteins on phylogenetic trees, reducing the high rates of false positives observed in conventional across-species methods that do not explicitly incorporate a phylogeny. We show, in a dataset of 10,551 protein pairs, that the phylogenetic method improves by up to 35% on across-species analyses at identifying known functionally linked proteins. The method shows that protein pairs with at least two to three correlated events of gain or loss are almost certainly functionally linked. Contingent evolution, in which one gene's presence or absence depends upon the presence of another, can also be detected phylogenetically, and may identify genes whose functional significance depends upon its interaction with other genes. Incorporating phylogenetic information improves the prediction of functional linkages. The improvement derives from having a lower rate of false positives and from detecting trends that across-species analyses miss. Phylogenetic methods can easily be incorporated into the screening of large-scale bioinformatics datasets to identify sets of protein links and to characterise gene networks.

A phylogenetic mixture model for detecting pattern-heterogeneity in gene sequence or character-state data

We describe a general likelihood-based 'mixture model' for inferring phylogenetic trees from gene-sequence or other character-state data. The model accommodates cases in which different sites in the alignment evolve in qualitatively distinct ways, but does not require prior knowledge of these patterns or partitioning of the data. We call this qualitative variability in the pattern of evolution across sites "pattern-heterogeneity" to distinguish it from both a homogenous process of evolution and from one characterized principally by differences in rates of evolution. We present studies to show that the model correctly retrieves the signals of pattern-heterogeneity from simulated gene-sequence data, and we apply the method to protein-coding genes and to a ribosomal 12S data set. The mixture model outperforms conventional partitioning in both these data sets. We implement the mixture model such that it can simultaneously detect rate- and pattern-heterogeneity. The model simplifies to a homogeneous model or a rate- variability model as special cases, and therefore always performs at least as well as these two approaches, and often considerably improves upon them. We make the model available within a Bayesian Markov-chain Monte Carlo framework for phylogenetic inference, as an easy-to-use computer program.

Phylogenetic Validation of the Genera Angomonas and Strigomonas of Trypanosomatids Harboring Bacterial Endosymbionts with the Description of New Species of Trypanosomatids and of Proteobacterial Symbionts

We comparatively examined the nutritional, molecular and optical and electron microscopical characteristics of reference species and new isolates of trypanosomatids harboring bacterial endosymbionts. Sequencing of the V7V8 region of the small subunit of the ribosomal RNA (SSU rRNA) gene distinguished six major genotypes among the 13 isolates examined. The entire sequences of the SSU rRNA and glycosomal glyceraldehyde phosphate dehydrogenase (gGAPDH) genes were obtained for phylogenetic analyses. In the resulting phylogenetic trees, the symbiont-harboring species clustered as a major clade comprising two subclades that corresponded to the proposed genera Angomonas and Strigomonas. The genus Angomonas comprised 10 flagellates including former Crithidia deanei and C. desouzai plus a new species. The genus Strigomonas included former Crithidia oncopelti and Blastocrithidia cuiicis plus a new species. Sequences from the internal transcribed spacer of ribosomal DNA (ITS rDNA) and size polymorphism of kinetoplast DNA (kDNA) minicircles revealed considerable genetic heterogeneity within the genera Angomonas and Strigomonas. Phylogenetic analyses based on 16S rDNA and ITS rDNA sequences demonstrated that all of the endosymbionts belonged to the Betaproteobacteria and revealed three new species. The congruence of the phylogenetic trees of trypanosomatids and their symbionts support a co-divergent host-symbiont evolutionary history. (C) 2011 Elsevier GmbH. All rights reserved.

Genes of cathepsin L-like proteases in Trypanosoma rangeli isolates: Markers for diagnosis, genotyping and phylogenetic relationships

We have sequenced genes encoding cathepsin L-like (CatL-like) cysteine proteases from isolates of Trypanosoma rangeli from humans, wild mammals and Rhodnius species of Central and South America. Phylogenetic trees of sequences encoding mature CatL-like enzymes of T rangeli and homologous genes from other trypanosomes, Leishmania spp. and bodonids positioned sequences of T rangeli (rangelipain) closest to T cruzi (cruzipain). Phylogenetic tree of kinetoplastids based on sequences of CatL-like was totally congruent with those derived from SSU rRNA and gGAPDH genes. Analysis of sequences from the CatL-like catalytic domains of 17 isolates representative of the overall phylogenetic diversity and geographical range of T rangeli supported all the lineages (A-D) previously defined using ribosomal and spliced leader genes. Comparison of the proteolytic activities of T rangeli isolates revealed heterogeneous banding profiles of cysteine proteases in gelatin gels, with differences even among isolates of the same lineage. CatL-like sequences proved to be excellent targets for diagnosis and genotyping of T rangeli by PCR. Data from CatL-like encoding genes agreed with results from previous studies of kDNA markers, and ribosomal and spliced leader genes, thereby corroborating clonal evolution, independent transmission cycles and the divergence of T rangeli lineages associated with sympatric species of Rhodnius. (c) 2009 Elsevier B.V. All rights reserved.

Phylogenetic Analyses Based on Small Subunit rRNA and Glycosomal Glyceraldehyde-3-Phosphate Dehydrogenase Genes and Ultrastructural Characterization of Two Snake Trypanosomes: Trypanosoma serpentis n. sp from Pseudoboa nigra and Trypanosoma cascavelli from Crotalus durissus terrificus

We sequenced the small subunit (SSU) rRNA and glycosomal glyceraldehyde-3-phosphate dehydrogenase (gGAPDH) genes of two trypanosomes isolated from the Brazilian snakes Pseudoboa nigra and Crotalus durissus terrificus. Trypanosomes were cultured and their morphometrical and ultrastructural features were characterized by light microscopy and scanning and transmission electron microscopy. Phylogenetic trees inferred using independent or combined SSU rRNA and gGAPDH data sets always clustered the snake trypanosomes together in a clade closest to lizard trypanosomes, forming a strongly supported monophyletic assemblage (i.e. lizard-snake clade). The positioning in the phylogenetic trees and the barcoding based on the variable V7-V8 region of the SSU rRNA, which showed high sequence divergences, allowed us to classify the isolates from distinct snake species as separate species. The isolate from P. nigra is described as a new species, Trypanosoma serpentis n. sp., whereas the isolate from C. d. terrificus is redescribed here as Trypanosoma cascavelli.

A Phylogenetic Lineage of Closely Related Trypanosomes (Trypanosomatidae, Kinetoplastida) of Anurans and Sand Flies (Psychodidae, Diptera) Sharing the Same Ecotopes in Brazilian Amazonia

Analysis of the phylogenetic relationships among trypanosomes from vertebrates and invertebrates disclosed a new lineage of trypanosomes circulating among anurans and sand flies that share the same ecotopes in Brazilian Amazonia. This assemblage of closely related trypanosomes was determined by comparing whole SSU rDNA sequences of anuran trypanosomes from the Brazilian biomes of Amazonia, the Pantanal, and the Atlantic Forest and from Europe, North America, and Africa, and from trypanosomes of sand flies from Amazonia. Phylogenetic trees based on maximum likelihood and parsimony corroborated the positioning of all new anuran trypanosomes in the aquatic clade but did not support the monophyly of anuran trypanosomes. However, all analyses always supported four major clades (An01-04) of anuran trypanosomes. Clade An04 is composed of trypanosomes from exotic anurans. Isolates in clades An01 and An02 were from Brazilian frogs and toads captured in the three biomes studied, Amazonia, the Pantanal and the Atlantic Forest. Clade An01 contains mostly isolates from Hylidae whereas clade An02 comprises mostly isolates from Bufonidae; and clade An03 contains trypanosomes from sand flies and anurans of Bufonidae, Leptodactylidae, and Leiuperidae exclusively from Amazonia. To our knowledge, this is the first study describing morphological and growth features, and molecular phylogenetic affiliation of trypanosomes from anurans and phlebotomines, incriminating these flies as invertebrate hosts and probably also as important vectors of Amazonian terrestrial anuran trypanosomes.

Phylogenetic, morphological and behavioural analyses support host switching of Trypanosoma (Herpetosoma) lewisi from domestic rats to primates

We characterized four Brazilian trypanosomes isolated from domestic rats and three from captive nonhuman primates that were morphologically similar to T. lewisi, a considered non-pathogenic species restricted to rodents and transmitted by fleas, despite its potential pathogenicity for infants. These isolates were identified as T. lewisi by barcoding using V7V8 SSU rDNA sequences. In inferred phylogenetic trees, all isolates clustered tightly with reference T. lewisi and T. lewisi-like trypanosomes from Europe, Asia and Africa and despite their high sequence conservation formed a homogeneous clade separate from other species of the subgenus T. (Herpetosoma). With the aim of clearly resolving the relationships between the Brazilian isolates from domestic rats and primates, we compared sequences from more polymorphic ITS rDNA. Results corroborated that isolates from Brazilian rats and monkeys were indeed of the same species and quite close to T. lewisi isolates of humans and rats from different geographical regions. Morphology of the monkey isolates and their behaviour in culture and in experimentally infected rats were also compatible with T. lewisi. However, infection with T. lewisi is rare among monkeys. We have examined more than 200 free-ranging and 160 captive monkeys and found only three infected individuals among the monkeys held in captivity. The findings of this work suggest that proximity of monkeys and infected rats and their exposure to infected fleas may be responsible for the host switching of T. Iewisi from their natural rodent species to primates. This and previous studies reporting T. lewisi in humans suggest that this trypanosome can cause sporadic and opportunistic fleaborne infection in primates. (C) 2010 Elsevier B.V. All rights reserved.

Genetic and phylogenetic analysis of glycoprotein of rabies virus isolated from several species in Brazil

Genetic and phylogenetic analyses of the region containing the glycoprotein (G) gene, which is related to pathogenicity and antigenicity, and the G-L intergenic region were carried out in 14 Brazilian rabies virus isolates. The isolates were classified as dog-related rabies virus (DRRV) or vampire bat-related rabies virus (VRRV), by nucleoprotein (N) analysis. The nucleotide and amino acid (AA) homologies of the area containing the G protein gene and G-L intergenic region were generally lower than those of the ectodomain. In both regions, nucleotide and deduced AA homologies were lower among VRRVs than among DRRVs. There were AA differences between DRRV and VRRV at 3 antigenic sites and epitopes (IIa, WB+ and III), suggesting that DRRV and VRRV can be distinguished by differences of antigenicity. In a comparison of phylogenetic trees between the ectodomain and the area containing the G protein gene and G-L intergenic region, the branching patterns of the chiropteran and carnivoran rabies virus groups differed, whereas there were clear similarities in patterns within the DRRV and VRRV groups. Additionally, the VRRV isolates were more closely related to chiropteran strains isolated from Latin America than to Brazilian DRRV. These results indicate that Brazilian rabies virus isolates can be classified as DRRV or VRRV by analysis of the G gene and the G-L intergenic region, as well as by N gene analysis.

EthoSeq: A tool for phylogenetic analysis and data mining in behavioral sequences

This article introduces the software program called EthoSeq, which is designed to extract probabilistic behavioral sequences (tree-generated sequences, or TGSs) from observational data and to prepare a TGS-species matrix for phylogenetic analysis. The program uses Graph Theory algorithms to automatically detect behavioral patterns within the observational sessions. It includes filtering tools to adjust the search procedure to user-specified statistical needs. Preliminary analyses of data sets, such as grooming sequences in birds and foraging tactics in spiders, uncover a large number of TGSs which together yield single phylogenetic trees. An example of the use of the program is our analysis of felid grooming sequences, in which we have obtained 1,386 felid grooming TGSs for seven species, resulting in a single phylogeny. These results show that behavior is definitely useful in phylogenetic analysis. EthoSeq simplifies and automates such analyses, uncovers much of the hidden patterns of long behavioral sequences, and prepares this data for further analysis with standard phylogenetic programs. We hope it will encourage many empirical studies on the evolution of behavior.