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.

The serine repeat antigen (SERA) gene family phylogeny in Plasmodium: the impact of GC content and reconciliation of gene and species trees.

Plasmodium falciparum is the parasite responsible for the most acute form of malaria in humans. Recently, the serine repeat antigen (SERA) in P. falciparum has attracted attention as a potential vaccine and drug target, and it has been shown to be a member of a large gene family. To clarify the relationships among the numerous P. falciparum SERAs and to identify orthologs to SERA5 and SERA6 in Plasmodium species affecting rodents, gene trees were inferred from nucleotide and amino acid sequence data for 33 putative SERA homologs in seven different species. (A distance method for nucleotide sequences that is specifically designed to accommodate differing GC content yielded results that were largely compatible with the amino acid tree. Standard-distance and maximum-likelihood methods for nucleotide sequences, on the other hand, yielded gene trees that differed in important respects.) To infer the pattern of duplication, speciation, and gene loss events in the SERA gene family history, the resulting gene trees were then "reconciled" with two competing Plasmodium species tree topologies that have been identified by previous phylogenetic studies. Parsimony of reconciliation was used as a criterion for selecting a gene tree/species tree pair and provided (1) support for one of the two species trees and for the core topology of the amino acid-derived gene tree, (2) a basis for critiquing fine detail in a poorly resolved region of the gene tree, (3) a set of predicted "missing genes" in some species, (4) clarification of the relationship among the P. falciparum SERA, and (5) some information about SERA5 and SERA6 orthologs in the rodent malaria parasites. Parsimony of reconciliation and a second criterion--implied mutational pattern at two key active sites in the SERA proteins-were also seen to be useful supplements to standard "bootstrap" analysis for inferred topologies.

Gene and species trees of a Neotropical group of treefrogs: Genetic diversification in the Brazilian Atlantic Forest and the origin of a polyploid species

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Quest for Orthologs Entails Quest for Tree of Life: In Search of the Gene Stream.

Quest for Orthologs (QfO) is a community effort with the goal to improve and benchmark orthology predictions. As quality assessment assumes prior knowledge on species phylogenies, we investigated the congruency between existing species trees by comparing the relationships of 147 QfO reference organisms from six Tree of Life (ToL)/species tree projects: The National Center for Biotechnology Information (NCBI) taxonomy, Opentree of Life, the sequenced species/species ToL, the 16S ribosomal RNA (rRNA) database, and trees published by Ciccarelli et al. (Ciccarelli FD, et al. 2006. Toward automatic reconstruction of a highly resolved tree of life. Science 311:1283-1287) and by Huerta-Cepas et al. (Huerta-Cepas J, Marcet-Houben M, Gabaldon T. 2014. A nested phylogenetic reconstruction approach provides scalable resolution in the eukaryotic Tree Of Life. PeerJ PrePrints 2:223) Our study reveals that each species tree suggests a different phylogeny: 87 of the 146 (60%) possible splits of a dichotomous and rooted tree are congruent, while all other splits are incongruent in at least one of the species trees. Topological differences are observed not only at deep speciation events, but also within younger clades, such as Hominidae, Rodentia, Laurasiatheria, or rosids. The evolutionary relationships of 27 archaea and bacteria are highly inconsistent. By assessing 458,108 gene trees from 65 genomes, we show that consistent species topologies are more often supported by gene phylogenies than contradicting ones. The largest concordant species tree includes 77 of the QfO reference organisms at the most. Results are summarized in the form of a consensus ToL (http://swisstree.vital-it.ch/species_tree) that can serve different benchmarking purposes.

Phylogenetic relationships among species of Anopheles (Nyssorhynchus) (Diptera, Culicidae) based on nuclear and mitochondrial gene sequences

Phylogenetic relationships among 21 species of mosquitoes in subgenus Nyssorhynchus were inferred from the nuclear white and mitochondrial NADH dehydrogenase subunit 6 (ND6) genes. Bayestan phylogenetic methods found that none of the three Sections within Nyssorhynchus (Albimanus, Argyritarsis, Myzorhynchella) were supported in all analyses, although Myzorhynchella was found to be monophyletic at the combined genes Within the Albimanus Section the monophyly of the Stroder Subgroup was strongly supported and within the Myzorhynchella Section Anopheles anrunesi and An lutzu formed a strongly supported monophyletic group The epidemiologically significant Albitarsis Complex showed evidence of paraphyly (relative to An lanet-Myzorhynchella) and discordance across gene trees, and the previously synonomized species of An. dunhami and An goeldii were recovered as sister species Finally, there was evidence of complexes in several species, including An antunesi, An deaneorum, and An. strodei (c) 2010 Elsevier B.V. All rights reserved

Tracking molecular evolution of photosynthesis by characterization of a major photosynthesis gene cluster from Heliobacillus mobilis

A DNA sequence has been obtained for a 35.6-kb genomic segment from Heliobacillus mobilis that contains a major cluster of photosynthesis genes. A total of 30 ORFs were identified, 20 of which encode enzymes for bacteriochlorophyll and carotenoid biosynthesis, reaction-center (RC) apoprotein, and cytochromes for cyclic electron transport. Donor side electron-transfer components to the RC include a putative RC-associated cytochrome c553 and a unique four-large-subunit cytochrome bc complex consisting of Rieske Fe-S protein (encoded by petC), cytochrome b6 (petB), subunit IV (petD), and a diheme cytochrome c (petX). Phylogenetic analysis of various photosynthesis gene products indicates a consistent grouping of oxygenic lineages that are distinct and descendent from anoxygenic lineages. In addition, H. mobilis was placed as the closest relative to cyanobacteria, which form a monophyletic origin to chloroplast-based photosynthetic lineages. The consensus of the photosynthesis gene trees also indicates that purple bacteria are the earliest emerging photosynthetic lineage. Our analysis also indicates that an ancient gene-duplication event giving rise to the paralogous bchI and bchD genes predates the divergence of all photosynthetic groups. In addition, our analysis of gene duplication of the photosystem I and photosystem II core polypeptides supports a “heterologous fusion model” for the origin and evolution of oxygenic photosynthesis.

Root of the universal tree of life based on ancient aminoacyl-tRNA synthetase gene duplications.

Universal trees based on sequences of single gene homologs cannot be rooted. Iwabe et al. [Iwabe, N., Kuma, K.-I., Hasegawa, M., Osawa, S. & Miyata, T. (1989) Proc. Natl. Acad. Sci. USA 86, 9355-9359] circumvented this problem by using ancient gene duplications that predated the last common ancestor of all living things. Their separate, reciprocally rooted gene trees for elongation factors and ATPase subunits showed Bacteria (eubacteria) as branching first from the universal tree with Archaea (archaebacteria) and Eucarya (eukaryotes) as sister groups. Given its topical importance to evolutionary biology and concerns about the appropriateness of the ATPase data set, an evaluation of the universal tree root using other ancient gene duplications is essential. In this study, we derive a rooting for the universal tree using aminoacyl-tRNA synthetase genes, an extensive multigene family whose divergence likely preceded that of prokaryotes and eukaryotes. An approximately 1600-bp conserved region was sequenced from the isoleucyl-tRNA synthetases of several species representing deep evolutionary branches of eukaryotes (Nosema locustae), Bacteria (Aquifex pyrophilus and Thermotoga maritima) and Archaea (Pyrococcus furiosus and Sulfolobus acidocaldarius). In addition, a new valyl-tRNA synthetase was characterized from the protist Trichomonas vaginalis. Different phylogenetic methods were used to generate trees of isoleucyl-tRNA synthetases rooted by valyl- and leucyl-tRNA synthetases. All isoleucyl-tRNA synthetase trees showed Archaea and Eucarya as sister groups, providing strong confirmation for the universal tree rooting reported by Iwabe et al. As well, there was strong support for the monophyly (sensu Hennig) of Archaea. The valyl-tRNA synthetase gene from Tr. vaginalis clustered with other eukaryotic ValRS genes, which may have been transferred from the mitochondrial genome to the nuclear genome, suggesting that this amitochondrial trichomonad once harbored an endosymbiotic bacterium.

Phylogeography of the common vampire bat (Desmodus rotundus): Marked population structure, Neotropical Pleistocene vicariance and incongruence between nuclear and mtDNA markers

Background: The common vampire bat Desmodus rotundus is an excellent model organism for studying ecological vicariance in the Neotropics due to its broad geographic range and its preference for forested areas as roosting sites. With the objective of testing for Pleistocene ecological vicariance, we sequenced a mitocondrial DNA (mtDNA) marker and two nuclear markers (RAG2 and DRB) to try to understand how Pleistocene glaciations affected the distribution of intraspecific lineages in this bat. Results: Five reciprocally monophyletic clades were evident in the mitochondrial gene tree, and in most cases with high bootstrap support: Central America (CA), Amazon and Cerrado (AMC), Pantanal (PAN), Northern Atlantic Forest (NAF) and Southern Atlantic Forest (SAF). The Atlantic forest clades formed a monophyletic clade with high bootstrap support, creating an east/west division for this species in South America. On the one hand, all coalescent and non-coalescent estimates point to a Pleistocene time of divergence between the clades. On the other hand, the nuclear markers showed extensive sharing of haplotypes between distant localities, a result compatible with male-biased gene flow. In order to test if the disparity between the mitochondrial and nuclear markers was due to the difference in mutation rate and effective size, we performed a coalescent simulation to examine the feasibility that, given the time of separation between the observed lineages, even with a gene flow rate close to zero, there would not be reciprocal monophyly for a neutral nuclear marker. We used the observed values of theta and an estimated mutation rate for the nuclear marker gene to perform 1000 iterations of the simulation. The results of this simulation were inconclusive: the number of iterations with and without reciprocal monophyly of one or more clades are similar. Conclusions: We therefore conclude that the pattern exhibited by the common vampire bat, with marked geographical structure for a mitochondrial marker and no phylogeographic structure for nuclear markers is compatible with a historical scenario of complete isolation of refuge-like populations during the Pleistocene. The results on demographic history on this species is compatible with the Carnaval-Moritz model of Pleistocene vicariance, with demographic expansions in the southern Atlantic forest.

Diversification and Species Boundaries of Rhinebothrium (Cestoda; Rhinebothriidea) in South American Freshwater Stingrays (Batoidea; Potamotrygonidae)

Background: Neotropical freshwater stingrays (Batoidea: Potamotrygonidae) host a diverse parasite fauna, including cestodes. Both cestodes and their stingray hosts are marine-derived, but the taxonomy of this host/parasite system is poorly understood. Methodology: Morphological and molecular (Cytochrome oxidase I) data were used to investigate diversity in freshwater lineages of the cestode genus Rhinebothrium Linton, 1890. Results were based on a phylogenetic hypothesis for 74 COI sequences and morphological analysis of over 400 specimens. Cestodes studied were obtained from 888 individual potamotrygonids, representing 14 recognized and 18 potentially undescribed species from most river systems of South America. Results: Morphological species boundaries were based mainly on microthrix characters observed with scanning electron microscopy, and were supported by COI data. Four species were recognized, including two redescribed (Rhinebothrium copianullum and R. paratrygoni), and two newly described (R. brooksi n. sp. and R. fulbrighti n. sp.). Rhinebothrium paranaensis Menoret & Ivanov, 2009 is considered a junior synonym of R. paratrygoni because the morphological features of the two species overlap substantially. The diagnosis of Rhinebothrium Linton, 1890 is emended to accommodate the presence of marginal longitudinal septa observed in R. copianullum and R. brooksi n. sp. Patterns of host specificity and distribution ranged from use of few host species in few river basins, to use of as many as eight host species in multiple river basins. Significance: The level of intra-specific morphological variation observed in features such as total length and number of proglottids is unparalleled among other elasmobranch cestodes. This is attributed to the large representation of host and biogeographical samples. It is unclear whether the intra-specific morphological variation observed is unique to this freshwater system. Nonetheless, caution is urged when using morphological discontinuities to delimit elasmobranch cestode species because the amount of variation encountered is highly dependent on sample size and/or biogeographical representation.

Molecular phylogenetics of Chaetodon and the Chaetodontidae (Teleostei : Perciformes) with reference to morphology

Butterflyfish are colourful, pan-tropical coastal fish that are important and distinctive members of coral reef communities. A successful systematic scheme and a robust phylogeny is considered essential in understanding further their biogeography and ecology, although recent cladistic treatments of butterflyfish phylogeny, based on soft tissue and bone morphology and coded at the generic and subgeneric levels, differ in character coding and subsequently tree topology. This study provides an independent test of the morphologically based hypotheses, using molecular systematic data from two partial mitochondrial gene fragments, cytochrome b (cytb) and small subunit rRNA (rrnS), for 52 ingroup chaetodontids and seven pomacanthids used to root the molecular trees. Individual gene trees were largely compatible and a combined molecular phylogeny, inferred from Bayesian analysis, was used to test alternative hypotheses suggested by morphological analyses. The tree was also used to map the latest morphological matrix in order to evaluate potential synapomorphies for various nodes defining butterflyfish interrelationships. A clade comprised of Chelmon and Coradion was sister group to other chaetodontids. Heniochus and Hemitaurichthys were each resolved as monophyletic groups, and as sister taxa Of the taxa sampled, Prognothodes was resolved as the sister genus to Chaeotodon. Of the ten Chaetodon subgenera sampled, all were monophyletic but their interrelationships differed significantly from that inferred from morphological characters. Lepidochaetodon was the most basal subgenus followed by Exornator and the remaining subgenera. Molecular data support the sister group relationship between Corallochaetodon and Citharoedus suggested by morphology, but major differences occur among the remaining more derived taxa. Chaetodon trifascialis and C. oligacanthus were resolved as sister taxa adding weight to the inclusion of the latter in C. Megaprotodon. Of those pairs of taxa known to hybridize and sampled with molecular data, all were closely related phylogenetically, except those hybrids known to occur in the Rabdophorus subgenus. Two base changes separated C. pelewensis from C. paucifasciatus which have been regarded previously as a single species. Cytb provided greater resolution than rrnS and will likely provide additional resolution with greater taxon sampling.

Platyhelminth phylogenetics - a key to understanding parasitism?

The comparative method, the inference of biological processes from phylogenetic patterns, is founded on the reliability of the phylogenetic tree. In attempting to apply the comparative method to the understanding of the evolution of parasitism in the phylum Platyhelminthes, we have highlighted several points we consider to be of value along with many problems. We discuss four of these topics. Firstly, we view the group at a phylum level, in particular discussing the importance of establishing the sister taxon to the obligate parasite group, the Neodermata, for addressing such questions as the monophyly, parasitism or the endo or ectoparasitic nature of the early parasites. The variety of non-congruent phylogenetic trees presented so far, utilising either or both morphological and molecular data, gives rise to the suggestion that any evolutionary scenarios presented at this stage be treated as interesting hypotheses rather than well-supported theories. Our second point of discussion is the conflict between morphological and molecular estimates of monogenean evolution. The Monogenea presents several well-established morphological autapomorphies, such that morphology consistently estimates the group as monophyletic, whereas molecular sequence analyses indicate paraphyly, with different genes giving different topologies. We discuss the problem of reconciling gene and species trees. Thirdly, we use recent phylogenetic results on the tapeworms to interpret the evolution of strobilation, proglottization, segmentation and scolex structure. In relation to the latter, the results presented indicate that the higher cestodes are diphyletic, with one branch difossate and the other tetrafossate. Finally, we use a SSU rDNA phylogenetic tree of the Trematoda as a basis for the discussion of an aspect of the digenean life-cycle, namely the nature of the first intermediate host. Frequent episodes of host-switching, between gastropod and bivalve hosts or even into annelids, are indicated.

Tapping the wealth of microbial data in high-throughput metabolic model reconstruction

Genome-scale metabolic models are valuable tools in the metabolic engineering process, based on the ability of these models to integrate diverse sources of data to produce global predictions of organism behavior. At the most basic level, these models require only a genome sequence to construct, and once built, they may be used to predict essential genes, culture conditions, pathway utilization, and the modifications required to enhance a desired organism behavior. In this chapter, we address two key challenges associated with the reconstruction of metabolic models: (a) leveraging existing knowledge of microbiology, biochemistry, and available omics data to produce the best possible model; and (b) applying available tools and data to automate the reconstruction process. We consider these challenges as we progress through the model reconstruction process, beginning with genome assembly, and culminating in the integration of constraints to capture the impact of transcriptional regulation. We divide the reconstruction process into ten distinct steps: (1) genome assembly from sequenced reads; (2) automated structural and functional annotation; (3) phylogenetic tree-based curation of genome annotations; (4) assembly and standardization of biochemistry database; (5) genome-scale metabolic reconstruction; (6) generation of core metabolic model; (7) generation of biomass composition reaction; (8) completion of draft metabolic model; (9) curation of metabolic model; and (10) integration of regulatory constraints. Each of these ten steps is documented in detail.

Future trypanosomatid phylogenies: refined homologies, supertrees and networks

There has been good progress in inferring the evolutionary relationships within trypanosomes from DNA data as until relatively recently, many relationships have remained rather speculative. Ongoing molecular studies have provided data that have adequately shown Trypanosoma to be monophyletic and, rather surprisingly, that there are sharply contrasting levels of genetic variation within and between the major trypanosomatid groups. There are still, however, areas of research that could benefit from further development and resolution that broadly fall upon three questions. Are the current statements of evolutionary homology within ribosomal small sub-unit genes in need of refinement? Can the published phylograms be expanded upon to form `supertrees' depicting further relationships? Does a bifurcating tree structure impose an untenable dogma upon trypanosomatid phylogeny where hybridisation or reticulate evolutionary steps have played a part? This article briefly addresses these three questions and, in so doing, hopes to stimulate further interest in the molecular evolution of the group.

Evolution of ascariasis in humans and pigs: a multi-disciplinary approach

The nematode parasite Ascaris lumbricoides infects the digestive tracts of over 1.4 billion people worldwide, and its sister species, Ascaris suum, has infected a countless number of domesticated and feral pigs. It is generally thought that the putative ancestor to these worms infected either humans or pigs, but with the advent of domestication, they had ample opportunity to jump to a new host and subsequently specialize and evolve into a new species. While nuclear DNA markers decisively separate the two populations, mitochondrial sequences reveal that three major haplotypes are found in A. suum and in A. lumbricoides, indicating either occasional hybridization, causing introgression of gene trees, or retention of polymorphism dating back to the original ancestral species. This article provides an illustration of the combined contribution of parasitology, archaeoparasitology, genetics and paleogenetics to the history of ascariasis. We specifically investigate the molecular history of ascariasis in humans by sequencing DNA from the eggs of Ascaris found among ancient archeological remains. The findings of this paleogenetic survey will explain whether the three mitochondrial haplotypes result from recent hybridization and introgression, due to intensive human-pig interaction, or whether their co-occurrence predates pig husbandry, perhaps dating back to the common ancestor. We hope to show how human-pig interaction has shaped the recent evolutionary history of this disease, perhaps revealing the identity of the ancestral host.

Selectome update: quality control and computational improvements to a database of positive selection.

Selectome (http://selectome.unil.ch/) is a database of positive selection, based on a branch-site likelihood test. This model estimates the number of nonsynonymous substitutions (dN) and synonymous substitutions (dS) to evaluate the variation in selective pressure (dN/dS ratio) over branches and over sites. Since the original release of Selectome, we have benchmarked and implemented a thorough quality control procedure on multiple sequence alignments, aiming to provide minimum false-positive results. We have also improved the computational efficiency of the branch-site test implementation, allowing larger data sets and more frequent updates. Release 6 of Selectome includes all gene trees from Ensembl for Primates and Glires, as well as a large set of vertebrate gene trees. A total of 6810 gene trees have some evidence of positive selection. Finally, the web interface has been improved to be more responsive and to facilitate searches and browsing.