Conflicting phylogenetic signals at the base of the metazoan tree

Shifts in evolutionary lability underlie independent gains and losses of root-nodule symbiosis in a single clade of plants

Root nodule symbiosis (RNS) is a complex trait that enables plants to access atmospheric nitrogen converted into usable forms through a mutualistic relationship with soil bacteria. Pinpointing the evolutionary origins of RNS is critical for understanding its genetic basis, but building this evolutionary context is complicated by data limitations and the intermittent presence of RNS in a single clade of ca. 30,000 species of flowering plants, i.e., the nitrogen-fixing clade (NFC). We developed the most extensive de novo phylogeny for the NFC and an RNS trait database to reconstruct the evolution of RNS. Our analysis identifies evolutionary rate heterogeneity associated with a two-step process: An ancestral precursor state transitioned to a more labile state from which RNS was rapidly gained at multiple points in the NFC. We illustrate how a two-step process could explain multiple independent gains and losses of RNS, contrary to recent hypotheses suggesting one gain and numerous losses, and suggest a broader phylogenetic and genetic scope may be required for genome-phenome mapping.

Phylogenomic analyses of Camellia support reticulate evolution among major clades

Camellia (Theaceae) is a morphologically highly diverse genus of flowering plants and includes many famous species with high economic value, and the phylogeny of this genus is not fully resolved. We used 95 transcriptomes from 87 Camellia species and identified 1481 low-copy genes to conduct a detailed analysis of the phylogeny of this genus according to various data-screening criteria. The results show that, very different from the two existing classification systems of Camellia, 87 species are grouped into 8 main clades and two independent species, and that all 8 clades except Clade 8 were strongly supported by almost all the coalescent or concatenated trees using different gene subsets. However, the relationships among these clades were weakly supported and different from analyses using different gene subsets; furthermore, they do not agree with the phylogeny from chloroplast genomes of Camellia. Additional analyses support reticulate evolution (probably resulting from introgression or hybridization) among some major Camellia lineages, providing explanation for extensive gene tree conflicts. Furthermore, we inferred that together with the formation of East Asian subtropical evergreen broad-leaved forests, Camellia underwent a radiative divergence of major clades at 23 ∼ 19 Ma in the late Miocene then had a subsequent species burst at 10 ∼ 5 Ma. Principal component and cluster analyses provides new insights into morphological changes underlying the evolution of Camellia and a reference to further clarify subgenus and sections of this genus. The comprehensive study here including a nuclear phylogeny and other analyses reveal the rapid evolutionary history of Camellia.

The performance of outgroup-free rooting under evolutionary radiations

Tree rooting implies a temporal dimension to phylogenies. Only after defining the position of the root node is that the ancestral-descendant relationship between branches can be fully deduced. Rooting has been usually carried out by employing evolutionarily close outgroup lineages, which is a drawback when these lineages are unavailable or unknown. Alternatively, outgroup-free rooting methods were proposed, which rely on the constancy of evolutionary rates to varying degrees. In this work we analyzed the performance of two of these methods, the midpoint rooting (MPR) and the minimal ancestor deviation (MAD), in rooting topologies evolved under challenging scenarios of fast evolutionary radiations derived from empirical data, characterized by short internal branches near the crown node. Considering all branch length combinations investigated, both methods exhibited average success rates below 50%, although MAD slightly outperformed MPR. Moreover, tree balance significantly impacted the relative performance of the methods. We found that, in four-taxa unrooted trees, the outcome of whether both methodologies will correctly root the tree can be roughly predicted by two simple dimensionless metrics: the coefficient of variation of the external branch lengths, and the ratio between the internal branch length to the total sum of branch lengths, which were employed to devise a general linear model that allowed calculating the probability of correct placing the root node for any four-taxa tree. We predicted that the performance of both outgroup-free rooting methods on loci representing the placental mammal radiation ranged between 50% and 75%.

Multilocus Sequence Typing Reveals Extensive Genetic Diversity of the Emerging Fungal Pathogen Scedosporium aurantiacum

Scedosporium spp. are the second most prevalent filamentous fungi after Aspergillus spp. recovered from cystic fibrosis (CF) patients in various regions of the world. Although invasive infection is uncommon prior to lung transplantation, fungal colonization may be a risk factor for invasive disease with attendant high mortality post-transplantation. Abundant in the environment, Scedosporium aurantiacum has emerged as an important fungal pathogen in a range of clinical settings. To investigate the population genetic structure of S. aurantiacum, a MultiLocus Sequence Typing (MLST) scheme was developed, screening 24 genetic loci for polymorphisms on a tester strain set. The six most polymorphic loci were selected to form the S. aurantiacum MLST scheme: actin (ACT), calmodulin (CAL), elongation factor-1α (EF1α), RNA polymerase subunit II (RPB2), manganese superoxide dismutase (SOD2), and β-tubulin (TUB). Among 188 global clinical, veterinary, and environmental strains, 5 to 18 variable sites per locus were revealed, resulting in 8 to 23 alleles per locus. MLST analysis observed a markedly high genetic diversity, reflected by 159 unique sequence types. Network analysis revealed a separation between Australian and non-Australian strains. Phylogenetic analysis showed two major clusters, indicating correlation with geographic origin. Linkage disequilibrium analysis revealed evidence of recombination. There was no clustering according to the source of the strains: clinical, veterinary, or environmental. The high diversity, especially amongst the Australian strains, suggests that S. aurantiacum may have originated within the Australian continent and was subsequently dispersed to other regions, as shown by the close phylogenetic relationships between some of the Australian sequence types and those found in other parts of the world. The MLST data are accessible at http://mlst.mycologylab.org. This is a joined publication of the ISHAM/ECMM working groups on “Scedosporium/Pseudallescheria Infections” and “Fungal Respiratory Infections in Cystic Fibrosis”.

Genome-Scale Profiling Reveals Noncoding Loci Carry Higher Proportions of Concordant Data

Many evolutionary relationships remain controversial despite whole-genome sequencing data. These controversies arise, in part, due to challenges associated with accurately modeling the complex phylogenetic signal coming from genomic regions experiencing distinct evolutionary forces. Here, we examine how different regions of the genome support or contradict well-established relationships among three mammal groups using millions of orthologous parsimony-informative biallelic sites (PIBS) distributed across primate, rodent, and Pecora genomes. We compared PIBS concordance percentages among locus types (e.g. coding sequences (CDS), introns, intergenic regions), and contrasted PIBS utility over evolutionary timescales. Sites derived from noncoding sequences provided more data and proportionally more concordant sites compared with those from CDS in all clades. CDS PIBS were also predominant drivers of tree incongruence in two cases of topological conflict. PIBS derived from most locus types provided surprisingly consistent support for splitting events spread across the timescales we examined, although we find evidence that CDS and intronic PIBS may, respectively and to a limited degree, inform disproportionately about older and younger splits. In this era of accessible wholegenome sequence data, these results:1) suggest benefits to more intentionally focusing on noncoding loci as robust data for tree inference and 2) reinforce the importance of accurate modeling, especially when using CDS data.

Confronting Sources of Systematic Error to Resolve Historically Contentious Relationships: A Case Study Using Gadiform Fishes (Teleostei, Paracanthopterygii, Gadiformes)

Reliable estimation of phylogeny is central to avoid inaccuracy in downstream macroevolutionary inferences. However, limitations exist in the implementation of concatenated and summary coalescent approaches, and Bayesian and full coalescent inference methods may not yet be feasible for computation of phylogeny using complicated models and large data sets. Here, we explored methodological (e.g., optimality criteria, character sampling, model selection) and biological (e.g., heterotachy, branch length heterogeneity) sources of systematic error that can result in biased or incorrect parameter estimates when reconstructing phylogeny by using the gadiform fishes as a model clade. Gadiformes include some of the most economically important fishes in the world (e.g., Cods, Hakes, and Rattails). Despite many attempts, a robust higher-level phylogenetic framework was lacking due to limited character and taxonomic sampling, particularly from several species-poor families that have been recalcitrant to phylogenetic placement. We compiled the first phylogenomic data set, including 14,208 loci (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$>$\end{document}2.8 M bp) from 58 species representing all recognized gadiform families, to infer a time-calibrated phylogeny for the group. Data were generated with a gene-capture approach targeting coding DNA sequences from single-copy protein-coding genes. Species-tree and concatenated maximum-likelihood (ML) analyses resolved all family-level relationships within Gadiformes. While there were a few differences between topologies produced by the DNA and the amino acid data sets, most of the historically unresolved relationships among gadiform lineages were consistently well resolved with high support in our analyses regardless of the methodological and biological approaches used. However, at deeper levels, we observed inconsistency in branch support estimates between bootstrap and gene and site coefficient factors (gCF, sCF). Despite numerous short internodes, all relationships received unequivocal bootstrap support while gCF and sCF had very little support, reflecting hidden conflict across loci. Most of the gene-tree and species-tree discordance in our study is a result of short divergence times, and consequent lack of informative characters at deep levels, rather than incomplete lineage sorting. We use this phylogeny to establish a new higher-level classification of Gadiformes as a way of clarifying the evolutionary diversification of the order. We recognize 17 families in five suborders: Bregmacerotoidei, Gadoidei, Ranicipitoidei, Merluccioidei, and Macrouroidei (including two subclades). A time-calibrated analysis using 15 fossil taxa suggests that Gadiformes evolved \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim $\end{document}79.5 Ma in the late Cretaceous, but that most extant lineages diverged after the Cretaceous–Paleogene (K-Pg) mass extinction (66 Ma). Our results reiterate the importance of examining phylogenomic analyses for evidence of systematic error that can emerge as a result of unsuitable modeling of biological factors and/or methodological issues, even when data sets are large and yield high support for phylogenetic relationships. [Branch length heterogeneity; Codfishes; commercial fish species; Cretaceous-Paleogene (K-Pg); heterotachy; systematic error; target enrichment.]

Embracing Uncertainty in Reconstructing Early Animal Evolution

The origin of animals, one of the major transitions in evolution, remains mysterious. Many key aspects of animal evolution can be reconstructed by comparing living species within a robust phylogenetic framework. However, uncertainty remains regarding the evolutionary relationships between two ancient animal lineages - sponges and ctenophores - and the remaining animal phyla. Comparative morphology and some phylogenomic analyses support the view that sponges represent the sister lineage to the rest of the animals, while other phylogenomic analyses support ctenophores, a phylum of carnivorous, gelatinous marine organisms, as the sister lineage. Here, we explore why different studies yield different answers and discuss the implications of the two alternative hypotheses for understanding the origin of animals. Reconstruction of ancient evolutionary radiations is devilishly difficult and will likely require broader sampling of sponge and ctenophore genomes, improved analytical strategies and critical analyses of the phylogenetic distribution and molecular mechanisms underlying apparently conserved traits. Rather than staking out positions in favor of the ctenophores-sister or the sponges-sister hypothesis, we submit that research programs aimed at understanding the biology of the first animals should instead embrace the uncertainty surrounding early animal evolution in their experimental designs.

Reconstructing Yeasts Phylogenies and Ancestors from Whole Genome Data

Phylogenetic studies aim to discover evolutionary relationships and histories. These studies are based on similarities of morphological characters and molecular sequences. Currently, widely accepted phylogenetic approaches are based on multiple sequence alignments, which analyze shared gene datasets and concatenate/coalesce these results to a final phylogeny with maximum support. However, these approaches still have limitations, and often have conflicting results with each other. Reconstructing ancestral genomes helps us understand mechanisms and corresponding consequences of evolution. Most existing genome level phylogeny and ancestor reconstruction methods can only process simplified real genome datasets or simulated datasets with identical genome content, unique genome markers, and limited types of evolutionary events. Here, we provide an alternative way to resolve phylogenetic problems based on analyses of real genome data. We use phylogenetic signals from all types of genome level evolutionary events, and overcome the conflicting issues existing in traditional phylogenetic approaches. Further, we build an automated computational pipeline to reconstruct phylogenies and ancestral genomes for two high-resolution real yeast genome datasets. Comparison results with recent studies and publications show that we reconstruct very accurate and robust phylogenies and ancestors. Finally, we identify and analyze the conserved syntenic blocks among reconstructed ancestral genomes and present yeast species.

Frequent gene flow blurred taxonomic boundaries of sections in Lilium L. (Liliaceae)

Gene flow between species may last a long time in plants. Reticulation inevitably causes difficulties in phylogenetic reconstruction. In this study, we looked into the genetic divergence and phylogeny of 20 Lilium species based on multilocus analyses of 8 genes of chloroplast DNA (cpDNA), the internally transcribed nuclear ribosomal DNA (nrITS) spacer and 20 loci extracted from the expressed sequence tag (EST) libraries of L. longiflorum Thunb. and L. formosanum Wallace. The phylogeny based on the combined data of the maternally inherited cpDNA and nrITS was largely consistent with the taxonomy of Lilium sections. This phylogeny was deemed the hypothetical species tree and uncovered three groups, i.e., Cluster A consisting of 4 taxa from the sections Pseudolirium and Liriotypus, Cluster B consisting of the 4 taxa from the sections Leucolirion, Archelirion and Daurolirion, and Cluster C comprising 10 taxa mostly from the sections Martagon and Sinomartagon. In contrast, systematic inconsistency occurred across the EST loci, with up to 19 genes (95%) displaying tree topologies deviating from the hypothetical species tree. The phylogenetic incongruence was likely attributable to the frequent genetic exchanges between species/sections, as indicated by the high levels of genetic recombination and the IMa analyses with the EST loci. Nevertheless, multilocus analysis could provide complementary information among the loci on the species split and the extent of gene flow between the species. In conclusion, this study not only detected frequent gene flow among Lilium sections that resulted in phylogenetic incongruence but also reconstructed a hypothetical species tree that gave insights into the nature of the complex relationships among Lilium species.

Ten genes and two topologies: an exploration of higher relationships in skipper butterflies (Hesperiidae)

Despite multiple attempts to infer the higher-level phylogenetic relationships of skipper butterflies (Family Hesperiidae), uncertainties in the deep clade relationships persist. The most recent phylogenetic analysis included fewer than 30% of known genera and data from three gene markers. Here we reconstruct the higher-level relationships with a rich sampling of ten nuclear and mitochondrial markers (7,726 bp) from 270 genera and find two distinct but equally plausible topologies among subfamilies at the base of the tree. In one set of analyses, the nuclear markers suggest two contrasting topologies, one of which is supported by the mitochondrial dataset. However, another set of analyses suggests mito-nuclear conflict as the reason for topological incongruence. Neither topology is strongly supported, and we conclude that there is insufficient phylogenetic evidence in the molecular dataset to resolve these relationships. Nevertheless, taking morphological characters into consideration, we suggest that one of the topologies is more likely.

SICLE: a high-throughput tool for extracting evolutionary relationships from phylogenetic trees

We present the phylogeny analysis software SICLE (Sister Clade Extractor), an easy-to-use, high-throughput tool to describe the nearest neighbors to a node of interest in a phylogenetic tree as well as the support value for the relationship. The application is a command line utility that can be embedded into a phylogenetic analysis pipeline or can be used as a subroutine within another C++ program. As a test case, we applied this new tool to the published phylome of Salinibacter ruber, a species of halophilic Bacteriodetes, identifying 13 unique sister relationships to S. ruber across the 4,589 gene phylogenies. S. ruber grouped with bacteria, most often other Bacteriodetes, in the majority of phylogenies, but 91 phylogenies showed a branch-supported sister association between S. ruber and Archaea, an evolutionarily intriguing relationship indicative of horizontal gene transfer. This test case demonstrates how SICLE makes it possible to summarize the phylogenetic information produced by automated phylogenetic pipelines to rapidly identify and quantify the possible evolutionary relationships that merit further investigation. SICLE is available for free for noncommercial use at http://eebweb.arizona.edu/sicle/.

Microsymbiont diversity and phylogeny of native bradyrhizobia associated with soybean (Glycine max L. Merr.) nodulation in South African soils

The genetic diversity and identification of slow- and fast-growing soybean root nodule bacterial isolates from different agro-climatic regions in Mpumalanga, Limpopo and Gauteng Provinces of South Africa were evaluated. The 16S-rDNA-RFLP analysis of 100 rhizobial isolates and eight reference type strains placed the isolates into six major clusters, and revealed their site-dependent genomic diversity. Sequence analysis of single and concatenated housekeeping genes (atpD, glnII and gyrB), as well as the symbiotic gene nifH captured a considerably higher level of genetic diversity and indicated the dominance of Bradyrhizobium diazoefficiens and Bradyrhizobium japonicum in Mpumalanga, Limpopo and Gauteng Provinces. Gene sequence similarities of isolates with type strains of Bradyrhizobium ranged from 97.3 to 100% for the 16S rDNA, and 83.4 to 100% for the housekeeping genes. The glnII gene phylogeny showed discordance with the other genes, suggesting lateral gene transfer or recombination events. Concatenated gene sequence analysis showed that most of the isolates did not align with known type strains and might represent new species from South Africa. This underscores the high genetic variability associated with soybean Bradyrhizobium in South African soils, and the presence of an important reservoir of novel soybean-nodulating bradyrhizobia in the country. In this study, the grouping of isolates was influenced by site origin, with Group I isolates originating from Limpopo Province and Groups II and III from Mpumlanga Province in the 16S rDNA-RFLP analysis.

Molecular signatures that are distinctive characteristics of the vertebrates and chordates and supporting a grouping of vertebrates with the tunicates

Members of the phylum Chordata and the subphylum Vertebrata are presently distinguished solely on the basis of morphological characteristics. The relationship of the vertebrates to the two non-vertebrate chordate subphyla is also a subject of debate. Analyses of protein sequences have identified multiple conserved signature indels (CSIs) that are specific for Chordata or for Vertebrata. Five CSIs in 4 important proteins are specific for the Vertebrata, whereas two other CSIs are uniquely found in all sequenced chordate species including Ciona intestinalis and Oikapleura dioica (Tunicates) as well as Branchiostoma floridae (Cephalochordates). The shared presence of these molecular signatures by all vertebrates/chordate species, but in no other animal taxa, strongly indicates that the genetic changes represented by the identified CSIs diagnose monophyletic groups. Two other discovered CSIs are uniquely shared by different vertebrate species and by either one (Ciona intestinalis) or both tunicate (Ciona and Oikapleura) species, but they are not found in Branchiostoma or other animal species. Specific presence of these CSIs in different vertebrates and either one or both tunicate species provides strong independent evidence that the vertebrate species are more closely related to the urochordates (tunicates) than to the cephalochordates.

Bacterial influences on animal origins

Animals evolved in seas teeming with bacteria, yet the influences of bacteria on animal origins are poorly understood. Comparisons among modern animals and their closest living relatives, the choanoflagellates, suggest that the first animals used flagellated collar cells to capture bacterial prey. The cell biology of prey capture, such as cell adhesion between predator and prey, involves mechanisms that may have been co-opted to mediate intercellular interactions during the evolution of animal multicellularity. Moreover, a history of bacterivory may have influenced the evolution of animal genomes by driving the evolution of genetic pathways for immunity and facilitating lateral gene transfer. Understanding the interactions between bacteria and the progenitors of animals may help to explain the myriad ways in which bacteria shape the biology of modern animals, including ourselves.

Horizontal gene transfer is a significant driver of gene innovation in dinoflagellates

The dinoflagellates are an evolutionarily and ecologically important group of microbial eukaryotes. Previous work suggests that horizontal gene transfer (HGT) is an important source of gene innovation in these organisms. However, dinoflagellate genomes are notoriously large and complex, making genomic investigation of this phenomenon impractical with currently available sequencing technology. Fortunately, de novo transcriptome sequencing and assembly provides an alternative approach for investigating HGT. We sequenced the transcriptome of the dinoflagellate Alexandrium tamarense Group IV to investigate how HGT has contributed to gene innovation in this group. Our comprehensive A. tamarense Group IV gene set was compared with those of 16 other eukaryotic genomes. Ancestral gene content reconstruction of ortholog groups shows that A. tamarense Group IV has the largest number of gene families gained (314-1,563 depending on inference method) relative to all other organisms in the analysis (0-782). Phylogenomic analysis indicates that genes horizontally acquired from bacteria are a significant proportion of this gene influx, as are genes transferred from other eukaryotes either through HGT or endosymbiosis. The dinoflagellates also display curious cases of gene loss associated with mitochondrial metabolism including the entire Complex I of oxidative phosphorylation. Some of these missing genes have been functionally replaced by bacterial and eukaryotic xenologs. The transcriptome of A. tamarense Group IV lends strong support to a growing body of evidence that dinoflagellate genomes are extraordinarily impacted by HGT.

Origin and evolution of glutamyl-prolyl tRNA synthetase WHEP domains reveal evolutionary relationships within Holozoa

Repeated domains in proteins that have undergone duplication or loss, and sequence divergence, are especially informative about phylogenetic relationships. We have exploited divergent repeats of the highly structured, 50-amino acid WHEP domains that join the catalytic subunits of bifunctional glutamyl-prolyl tRNA synthetase (EPRS) as a sequence-informed repeat (SIR) to trace the origin and evolution of EPRS in holozoa. EPRS is the only fused tRNA synthetase, with two distinct aminoacylation activities, and a non-canonical translation regulatory function mediated by the WHEP domains in the linker. Investigating the duplications, deletions and divergence of WHEP domains, we traced the bifunctional EPRS to choanozoans and identified the fusion event leading to its origin at the divergence of ichthyosporea and emergence of filozoa nearly a billion years ago. Distribution of WHEP domains from a single species in two or more distinct clades suggested common descent, allowing the identification of linking organisms. The discrete assortment of choanoflagellate WHEP domains with choanozoan domains as well as with those in metazoans supported the phylogenetic position of choanoflagellates as the closest sister group to metazoans. Analysis of clustering and assortment of WHEP domains provided unexpected insights into phylogenetic relationships amongst holozoan taxa. Furthermore, observed gaps in the transition between WHEP domain groupings in distant taxa allowed the prediction of undiscovered or extinct evolutionary intermediates. Analysis based on SIR domains can provide a phylogenetic counterpart to palaentological approaches of discovering “missing links” in the tree of life.

Inferring explicit weighted consensus networks to represent alternative evolutionary histories

BACKGROUND

The advent of molecular biology techniques and constant increase in availability of genetic material have triggered the development of many phylogenetic tree inference methods. However, several reticulate evolution processes, such as horizontal gene transfer and hybridization, have been shown to blur the species evolutionary history by causing discordance among phylogenies inferred from different genes.

METHODS

To tackle this problem, we hereby describe a new method for inferring and representing alternative (reticulate) evolutionary histories of species as an explicit weighted consensus network which can be constructed from a collection of gene trees with or without prior knowledge of the species phylogeny.

RESULTS

We provide a way of building a weighted phylogenetic network for each of the following reticulation mechanisms: diploid hybridization, intragenic recombination and complete or partial horizontal gene transfer. We successfully tested our method on some synthetic and real datasets to infer the above-mentioned evolutionary events which may have influenced the evolution of many species.

CONCLUSIONS

Our weighted consensus network inference method allows one to infer, visualize and validate statistically major conflicting signals induced by the mechanisms of reticulate evolution. The results provided by the new method can be used to represent the inferred conflicting signals by means of explicit and easy-to-interpret phylogenetic networks.

Phylogenetic relationships and genome divergence among the AA- genome species of the genus Oryza as revealed by 53 nuclear genes and 16 intergenic regions

Rapid radiations have long been regarded as the most challenging issue for elucidating poorly resolved phylogenies in evolutionary biology. The eight diploid AA- genome species in the genus Oryza represent a typical example of a closely spaced series of recent speciation events in plants. However, questions regarding when and how they diversified have long been an issue of extensive interest but remain a mystery. Here, a data set comprising >60 kb of 53 singleton fragments and 16 intergenic regions is used to perform phylogenomic analyses of all eight AA- genome species plus four diploid Oryza species with BB-, CC-, EE- and GG- genomes. We fully reconstruct phylogenetic relationships of AA- genome species with confidence. Oryza meridionalis, native to Australia, is found to be the earliest divergent lineage around 2.93 mya, whereas O. punctata, a BB- genome species, serves as the best outgroup to distinguish their phylogenetic relationships. They separated from O. punctata approximately 9.11 mya during the Miocene epoch, and subsequently radiated to generate the entire AA- genome lineage diversity. The success in resolving the phylogeny of AA- genome species highlights the potential of phylogenomics to determine their divergence and evolutionary histories.

Tubulin posttranslational modifications induced by cadmium in the sponge Clathrina clathrus

As sessile filter feeders, sponges are exposed to environmental stress due to pollutants of both anthropogenic and natural origins and are able to accumulate harmful substances. Thus, sponges are considered a good tool for the biomonitoring of coastal areas. In this study, we used biochemical and immunocytochemical analyses to provide new data on the cadmium-related changes in sponge cells. In particular, we analyzed the effects of different concentrations of cadmium on the microtubule network in the calcisponge Clathrina clathrus. Quantitative densitometry of the immunoblots showed that, while the levels of α- and β-tubulin remained relatively constant in C. clathrus when exposed to 1 and 5 μM CdCl2, there were progressive shifts in the levels of some tubulin isoforms. Exposure for 24h to sublethal concentrations of cadmium reduced the level of tyrosinated α-tubulin and enhanced the levels of acetylated and detyrosinated α-tubulin relative to the levels in controls. Confocal microscopy analysis of immunolabeled tissue sections showed that the inhibitory effect of cadmium was associated with a decrease in the labeling of the cells with a monoclonal antibody that recognizes tyrosinated α-tubulin. By contrast, the reactivity with a monoclonal antibody that recognizes acetylated α-tubulin and with a polyclonal antibody specific for detyrosinated α-tubulin was enhanced at the same time points. Because the acetylation and detyrosination of α-tubulin occur on stable microtubules, the marked enhancement of α-tubulin acetylation and detyrosination in Cd(2+)-treated cells indicates that divalent Cd ions stabilize microtubules. The possibility that Cd(2+) may increase the stability of cytoplasmic microtubules was tested by exposing Cd(2+)-treated cells to a cold temperature (0°C). As shown, the microtubule bundles induced by Cd(2+), which were labeled by the monoclonal antibodies against acetylated and detyrosinated α-tubulin, were resistant to cold.

Insights into bilaterian evolution from three spiralian genomes

Current genomic perspectives on animal diversity neglect two prominent phyla, the molluscs and annelids, that together account for nearly one-third of known marine species and are important both ecologically and as experimental systems in classical embryology^1–3. Here we describe the draft genomes of the owl limpet (Lottia gigantea), a marine polychaete (Capitella teleta) and a freshwater leech (Helobdella robusta), and compare them with other animal genomes to investigate the origin and diversification of bilaterians from a genomic perspective. We find that the genome organization, gene structure and functional content of these species are more similar to those of some invertebrate deuterostome genomes (for example, amphioxus and sea urchin) than those of other protostomes that have been sequenced to date (flies, nematodes and flatworms). The conservation of these genomic features enables us to expand the inventory of genes present in the last common bilaterian ancestor, establish the tripartite diversification of bilaterians using multiple genomic characteristics and identify ancient conserved long- and short-range genetic linkages across metazoans. Superimposed on this broadly conserved pan-bilaterian background we find examples of lineage-specific genome evolution, including varying rates of rearrangement, intron gain and loss, expansions and contractions of gene families, and the evolution of clade-specific genes that produce the unique content of each genome.

Rhizome of life, catastrophes, sequence exchanges, gene creations, and giant viruses: how microbial genomics challenges Darwin

Darwin's theory about the evolution of species has been the object of considerable dispute. In this review, we have described seven key principles in Darwin's book The Origin of Species and tried to present how genomics challenge each of these concepts and improve our knowledge about evolution. Darwin believed that species evolution consists on a positive directional selection ensuring the “survival of the fittest.” The most developed state of the species is characterized by increasing complexity. Darwin proposed the theory of “descent with modification” according to which all species evolve from a single common ancestor through a gradual process of small modification of their vertical inheritance. Finally, the process of evolution can be depicted in the form of a tree. However, microbial genomics showed that evolution is better described as the “biological changes over time.” The mode of change is not unidirectional and does not necessarily favors advantageous mutations to increase fitness it is rather subject to random selection as a result of catastrophic stochastic processes. Complexity is not necessarily the completion of development: several complex organisms have gone extinct and many microbes including bacteria with intracellular lifestyle have streamlined highly effective genomes. Genomes evolve through large events of gene deletions, duplications, insertions, and genomes rearrangements rather than a gradual adaptative process. Genomes are dynamic and chimeric entities with gene repertoires that result from vertical and horizontal acquisitions as well as de novo gene creation. The chimeric character of microbial genomes excludes the possibility of finding a single common ancestor for all the genes recorded currently. Genomes are collections of genes with different evolutionary histories that cannot be represented by a single tree of life (TOL). A forest, a network or a rhizome of life may be more accurate to represent evolutionary relationships among species.

The impact of taxon sampling on phylogenetic inference: a review of two decades of controversy

Over the past two decades, there has been a long-standing debate about the impact of taxon sampling on phylogenetic inference. Studies have been based on both real and simulated data sets, within actual and theoretical contexts, and using different inference methods, to study the impact of taxon sampling. In some cases, conflicting conclusions have been drawn for the same data set. The main questions explored in studies to date have been about the effects of using sparse data, adding new taxa, including more characters from genome sequences and using different (or concatenated) locus regions. These questions can be reduced to more fundamental ones about the assessment of data quality and the design guidelines of taxon sampling in phylogenetic inference experiments. This review summarizes progress to date in understanding the impact of taxon sampling on the accuracy of phylogenetic analysis.

Deep phylogeny and evolution of sponges (phylum Porifera)

Sponges (phylum Porifera) are a diverse taxon of benthic aquatic animals of great ecological, commercial, and biopharmaceutical importance. They are arguably the earliest-branching metazoan taxon, and therefore, they have great significance in the reconstruction of early metazoan evolution. Yet, the phylogeny and systematics of sponges are to some extent still unresolved, and there is an on-going debate about the exact branching pattern of their main clades and their relationships to the other non-bilaterian animals. Here, we review the current state of the deep phylogeny of sponges. Several studies have suggested that sponges are paraphyletic. However, based on recent phylogenomic analyses, we suggest that the phylum Porifera could well be monophyletic, in accordance with cladistic analyses based on morphology. This finding has many implications for the evolutionary interpretation of early animal traits and sponge development. We further review the contribution that mitochondrial genes and genomes have made to sponge phylogenetics and explore the current state of the molecular phylogenies of the four main sponge lineages (Classes), that is, Demospongiae, Hexactinellida, Calcarea, and Homoscleromorpha, in detail. While classical systematic systems are largely congruent with molecular phylogenies in the class Hexactinellida and in certain parts of Demospongiae and Homoscleromorpha, the high degree of incongruence in the class Calcarea still represents a challenge. We highlight future areas of research to fill existing gaps in our knowledge. By reviewing sponge development in an evolutionary and phylogenetic context, we support previous suggestions that sponge larvae share traits and complexity with eumetazoans and that the simple sedentary adult lifestyle of sponges probably reflects some degree of secondary simplification. In summary, while deep sponge phylogenetics has made many advances in the past years, considerable efforts are still required to achieve a comprehensive understanding of the relationships among and within the main sponge lineages to fully appreciate the evolution of this extraordinary metazoan phylum.

Statistics and truth in phylogenomics

Phylogenomics refers to the inference of historical relationships among species using genome-scale sequence data and to the use of phylogenetic analysis to infer protein function in multigene families. With rapidly decreasing sequencing costs, phylogenomics is becoming synonymous with evolutionary analysis of genome-scale and taxonomically densely sampled data sets. In phylogenetic inference applications, this translates into very large data sets that yield evolutionary and functional inferences with extremely small variances and high statistical confidence (P value). However, reports of highly significant P values are increasing even for contrasting phylogenetic hypotheses depending on the evolutionary model and inference method used, making it difficult to establish true relationships. We argue that the assessment of the robustness of results to biological factors, that may systematically mislead (bias) the outcomes of statistical estimation, will be a key to avoiding incorrect phylogenomic inferences. In fact, there is a need for increased emphasis on the magnitude of differences (effect sizes) in addition to the P values of the statistical test of the null hypothesis. On the other hand, the amount of sequence data available will likely always remain inadequate for some phylogenomic applications, for example, those involving episodic positive selection at individual codon positions and in specific lineages. Again, a focus on effect size and biological relevance, rather than the P value, may be warranted. Here, we present a theoretical overview and discuss practical aspects of the interplay between effect sizes, bias, and P values as it relates to the statistical inference of evolutionary truth in phylogenomics.

Optimal data partitioning and a test case for ray-finned fishes (Actinopterygii) based on ten nuclear loci

Data partitioning, the combined phylogenetic analysis of homogeneous blocks of data, is a common strategy used to accommodate heterogeneities in complex multilocus data sets. Variation in evolutionary rates and substitution patterns among sites are typically addressed by partitioning data by gene, codon position, or both. Excessive partitioning of the data, however, could lead to overparameterization; therefore, it seems critical to define the minimum numbers of partitions necessary to improve the overall fit of the model. We propose a new method, based on cluster analysis, to find an optimal partitioning strategy for multilocus protein-coding data sets. A heuristic exploration of alternative partitioning schemes, based on Bayesian and maximum likelihood (ML) criteria, is shown here to produce an optimal number of partitions. We tested this method using sequence data of 10 nuclear genes collected from 52 ray-finned fish (Actinopterygii) and four tetrapods. The concatenated sequences included 7995 nucleotide sites maximally split into 30 partitions defined a priori based on gene and codon position. Our results show that a model based on only 10 partitions defined by cluster analysis performed better than partitioning by both gene and codon position. Alternative data partitioning schemes also are shown to affect the topologies resulting from phylogenetic analysis, especially when Bayesian methods are used, suggesting that overpartitioning may be of major concern. The phylogenetic relationships among the major clades of ray-finned fish were assessed using the best data-partitioning schemes under ML and Bayesian methods. Some significant results include the monophyly of "Holostei" (Amia and Lepisosteus), the sister-group relationships between (1) esociforms and salmoniforms and (2) osmeriforms and stomiiforms, the polyphyly of Perciformes, and a close relationship of cichlids and atherinomorphs.

Molecular inference of phylogenetic relationships among Decapodiformes (Mollusca: Cephalopoda) with special focus on the squid order Oegopsida

Squids, cuttlefish and bobtail squids comprise the molluscan superorder Decapodiformes (Mollusca: Cephalopoda). Although these animals exemplify the morphological and ecological diversity seen in Cephalopoda, no previous study has focused resolving decapodiform relationships, particularly within Oegopsida, a large order comprised of pelagic squid. To further clarify the phylogenetic history of Decapodiformes, and Oegopsida in particular, molecular data for five genes (18S rRNA, 28S rRNA, Histone H3, 16S rRNA, COI) was collected for 90 taxa representing all major lineages and families and evaluated using parsimony, likelihood, and Bayesian analysis. Although ordinal relationships were sensitive to analytical method, several conclusions can be inferred: the pelagic order Myopsida is closely related to the benthic sepioids, whose relationships were ambiguous, and Bathyteuthoidea is distinct from Oegopsida. Within Oegopsida several clades are consistently recovered, some with previous morphological support (e.g. chiroteuthid, lepidoteuthid, histioteuthid families) while others suggest novel relationships (e.g. Architeuthidae+Neoteuthidae). This study, with its broad coverage of taxa, provides the first in-depth analysis of Decapodiformes with special focus on the morphologically and biogeographically diverse Oegopsida, confirms several sister-taxon relationships, and provides new hypotheses of cephalopod evolution in the open ocean.

Leveraging skewed transcript abundance by RNA-Seq to increase the genomic depth of the tree of life

Assembling the tree of life is a major goal of biology, but progress has been hindered by the difficulty and expense of obtaining the orthologous DNA required for accurate and fully resolved phylogenies. Next-generation DNA sequencing technologies promise to accelerate progress, but sequencing the genomes of hundreds of thousands of eukaryotic species remains impractical. Eukaryotic transcriptomes, which are smaller than genomes and biased toward highly expressed genes that tend to be conserved, could potentially provide a rich set of phylogenetic characters. We sampled the transcriptomes of 10 mosquito species by assembling 36-bp sequence reads into phylogenomic data matrices containing hundreds of thousands of orthologous nucleotides from hundreds of genes. Analysis of these data matrices yielded robust phylogenetic inferences, even with data matrices constructed from surprisingly few sequence reads. This approach is more efficient, data-rich, and economical than traditional PCR-based and EST-based methods and provides a scalable strategy for generating phylogenomic data matrices to infer the branches and twigs of the tree of life.

A phylogenomic approach to resolve the basal pterygote divergence

One of the most fascinating Bauplan transitions in the animal kingdom was the invention of insect wings, a change that also contributed to the success and enormous diversity of this animal group. However, the origin of insect flight and the relationships of basal winged insect orders are still controversial. Three hypotheses have been proposed to explain the phylogeny of winged insects: 1) the traditional Palaeoptera hypothesis (Ephemeroptera + Odonata, Neoptera), 2) the Metapterygota hypothesis (Ephemeroptera, Odonata + Neoptera), and 3) the Chiastomyaria hypothesis (Odonata, Ephemeroptera + Neoptera). Neither phylogenetic analyses of single genes nor even multiple marker systems (e.g., molecular markers + morphological characters) have yet been able to conclusively resolve basal pterygote divergences. A possible explanation for the lack of resolution is that the divergences took place in the mid-Devonian within a short period of time and attempts to solve this problem have been confounded by the major challenge of finding molecular markers to accurately track these short ancient internodes. Although phylogenomic data are available for Neoptera and some wingless (apterygote) orders, they are lacking for the crucial Odonata and Ephemeroptera orders. We adopt a multigene approach including data from two new expressed sequence tag projects-from the orders Ephemeroptera (Baetis sp.) and Odonata (Ischnura elegans)-to evaluate the potential of phylogenomic analyses in clarifying this unresolved issue. We analyzed two data sets that differed in represented taxa, genes, and overall sequence lengths: maxspe (15 taxa, 125 genes, and 31,643 amino acid positions) and maxgen (8 taxa, 150 genes, and 42,541 amino acid positions). Maximum likelihood and Bayesian inference analyses both place the Odonata at the base of the winged insects. Furthermore, statistical hypotheses testing rejected both the Palaeoptera and the Metapterygota hypotheses. The comprehensive molecular data set developed here provides conclusive support for odonates as the most basal winged insect order (Chiastomyaria hypothesis). Data quality assessment indicates that proteins involved in cellular processes and signaling harbor the most informative phylogenetic signal.

P-glycoprotein (multi-xenobiotic resistance) and heat shock protein gene expression in the reef coral Montastraea franksi in response to environmental toxicants

The deleterious impacts of marine pollutants on reef corals and their symbiotic algae are an important element of global coral reef decline. In the current study we examined the impacts of toxicants on the reef coral Montastraea franksi by analysing the expression of three stress-related genes belonging to the coral host, using Taqman real-time quantitative reverse transcription-PCR. Gene expression profiles of P-glycoprotein (or multi-xenobiotic resistance protein) (P-gp); heat shock protein 70 (Hsp70) and heat shock protein 90 (Hsp90) were examined following 4 and 8h exposures to the heavy metal copper (3, 10, 30 and 100 microgL(-1)) or the third generation oil dispersant Corexit9527 (1, 5, 10 and 50 ppm). Additionally, the expression of P-gp was examined following exposure to 0.5 and 5 microM concentrations of the chemotherapeutic drug vinblastine, a classic substrate of P-gp. The expression of P-gp increased significantly in corals treated with vinblastine and also increased following exposure to copper and Corexit9527. Hsp70, and to a lesser extent Hsp90 expression increased following exposure to copper and Corexit9527 indicating a general cellular stress response. Densities of symbiotic algae in the tissues of the corals did not change significantly during the experiments, nor was any loss or paling of coral tissues observed. These findings provide important insight into how corals defend themselves against pollution and complement ongoing initiatives developing molecular biomarkers of stress in reef-building corals.

The origins of multicellularity and the early history of the genetic toolkit for animal development

Multicellularity appeared early and repeatedly in life's history; its instantiations presumably required the confluence of environmental, ecological, and genetic factors. Comparisons of several independently evolved pairs of multicellular and unicellular relatives indicate that transitions to multicellularity are typically associated with increases in the numbers of genes involved in cell differentiation, cell-cell communication, and adhesion. Further examination of the DNA record suggests that these increases in gene complexity are the product of evolutionary innovation, tinkering, and expansion of genetic material. Arguably, the most decisive multicellular transition was the emergence of animals. Decades of developmental work have demarcated the genetic toolkit for animal multicellularity, a select set of a few hundred genes from a few dozen gene families involved in adhesion, communication, and differentiation. Examination of the DNA records of the earliest-branching animal phyla and their closest protist relatives has begun to shed light on the origins and assembly of this toolkit. Emerging data favor a model of gradual assembly, with components originating and diversifying at different time points prior to or shortly after the origin of animals.

The phylogenetic informativeness of nucleotide and amino acid sequences for reconstructing the vertebrate tree

To aid in future efforts to accurately reconstruct the vertebrate tree, a quantitative measure of phylogenetic informativeness was applied to nucleotide and amino acid sequences for a set of 11 genes. We identified orthologues and assembled published fossil-calibrated divergence times between taxa that had been sequenced for each gene. Rates of molecular evolution for each site were estimated to characterize the molecular evolutionary pattern of genes and to calculate the phylogenetic informativeness. The fast-evolving gene albumin yielded the highest informativeness over the period from 60 million years ago to 500 million years ago. In contrast, calmodulin yielded the lowest informativeness, presumably because functional constraint minimized substitutions in the amino acid sequence. The gene c-myc showed an intermediate level of informativeness. The nucleotide sequence of cytochrome b showed extremely high utility for recent epochs, but low utility for times before 100 million years ago. We ranked nine other genes for their utility during the epochs of the divergence of the muroid rodents, early placental mammals, early vertebrates, and early metazoa, yielding results consistent with, but more precise than, previous studies. Interestingly, DNA sequence always exceeded amino acid sequence in informativeness over all time scales, yet support values were at best moderately higher. For epochs not subject to strong phylogenetic conflict due to convergence, we advocate gleaning the additional power of the threefold increase in number of characters that is present for DNA sequences over resorting to the less noisy but less informative amino acid sequences.

From gene-scale to genome-scale phylogenetics: the data flood in, but the challenges remain

An important goal of phylogenetics is to be able to consistently and accurately reconstruct the historical patterns of cladogenesis among major organismic groups. Gene-scale phylogenetics is insufficient to attain this goal owing to the presence of poor resolution and incongruence in single--and few--gene phylogenies. The increasing availability of genome-scale amounts of data promises to overcome the insufficiency of gene-scale phylogenetics and uncover the genealogical tapestry uniting all living organisms with unprecedented accuracy. Here, we argue that a vast increase in data size alone--although necessary--may not be sufficient to achieve the desired accuracy for three reasons: (i) the existence of short stems in the tree of life, (ii) the saturation of phylogenetic signal in molecular sequences, and (iii) the effect of systematic error on phylogenetic inference. Devising strategies to ameliorate the effect of such challenges on sequence evolution will be critical to the success of current efforts to reconstruct the tree of life.

Phylogenetic placement of the enigmatic parasite, Polypodium hydriforme, within the Phylum Cnidaria

BACKGROUND

Polypodium hydriforme is a parasite with an unusual life cycle and peculiar morphology, both of which have made its systematic position uncertain. Polypodium has traditionally been considered a cnidarian because it possesses nematocysts, the stinging structures characteristic of this phylum. However, recent molecular phylogenetic studies using 18S rDNA sequence data have challenged this interpretation, and have shown that Polypodium is a close relative to myxozoans and together they share a closer affinity to bilaterians than cnidarians. Due to the variable rates of 18S rDNA sequences, these results have been suggested to be an artifact of long-branch attraction (LBA). A recent study, using multiple protein coding markers, shows that the myxozoan Buddenbrockia, is nested within cnidarians. Polypodium was not included in this study. To further investigate the phylogenetic placement of Polypodium, we have performed phylogenetic analyses of metazoans with 18S and partial 28S rDNA sequences in a large dataset that includes Polypodium and a comprehensive sampling of cnidarian taxa.

RESULTS

Analyses of a combined dataset of 18S and partial 28S sequences, and partial 28S alone, support the placement of Polypodium within Cnidaria. Removal of the long-branched myxozoans from the 18S dataset also results in Polypodium being nested within Cnidaria. These results suggest that previous reports showing that Polypodium and Myxozoa form a sister group to Bilateria were an artifact of long-branch attraction.

CONCLUSION

By including 28S rDNA sequences and a comprehensive sampling of cnidarian taxa, we demonstrate that previously conflicting hypotheses concerning the phylogenetic placement of Polypodium can be reconciled. Specifically, the data presented provide evidence that Polypodium is indeed a cnidarian and is either the sister taxon to Hydrozoa, or part of the hydrozoan clade, Leptothecata. The former hypothesis is consistent with the traditional view that Polypodium should be placed in its own cnidarian class, Polypodiozoa.

The early evolution of the phosphagen kinases--insights from choanoflagellate and poriferan arginine kinases

Arginine kinase (AK) is a member of a large family of phosphoryl transfer enzymes called phosphagen (guanidino) kinases. AKs are present in certain protozoans, sponges, cnidarians, and both lophotrochozoan and ecdysozoan protostomes. Another phosphagen kinase, creatine kinase (CK), is found in sponges, cnidarians, and both deuterostome and protostome groups but does not appear to be present in protozoans. To probe the early evolution of phosphagen kinases, we have amplified the cDNAs for AKs from three choanoflagellates and from the hexactinellid sponge Aphrocallistes beatrix and the demosponges Suberites fuscus and Microciona prolifera. Phylogenetic analysis using maximum likelihood of these choanoflagellate and sponge AKs with other AK sequences revealed that the AK from the choanoflagellate Monosiga brevicollis clusters with the AK from the glass sponge Aphrocallistes and is part of a larger cluster containing AKs from the demosponges Suberites and Microciona as well as basal and protostome invertebrates. In contrast, AKs from Codonosiga gracilis and Monosiga ovata form a distinct cluster apart from all other AK sequences. tBLASTn searches of the recently released M. brevicollis genome database showed that this species has three unique AK genes-one virtually identical to the M. brevicollis cDNA and the other two showing great similarity to C. gracilis and M. ovata AKs. Three distinct AK genes are likely present in choanoflagellates. Two of these AKs display extensive similarity to both CKs and an AK from sponges. Previous work has shown CK evolved from an AK-like ancestor prior to the divergence of sponges. The present results provide evidence suggesting that the initial gene duplication event(s) leading to the CK lineage may have occurred before the divergence of the choanoflagellate and animal lineages.

Fossils Impact as Hard as Living Taxa in Parsimony Analyses of Morphology

Systematists disagree whether data from fossils should be included in parsimony analyses. In a handful of well-documented cases, the addition of fossil data radically overturns a hypothesis of relationships based on extant taxa alone. Fossils can break up long branches and preserve character combinations closer in time to deep splitting events. However, fossils usually require more interpretation than extant taxa, introducing greater potential for spurious codings. Moreover, because fossils often have more "missing" codings, they are frequently accused of increasing numbers of MPTs, frustrating resolution and reducing support. Despite the controversy, remarkably little is known about the effects of fossils more generally. Here we provide the first systematic study, investigating empirically the behavior of fossil and extant taxa in 45 published morphological data sets. First-order jackknifing is used to determine the effects that each terminal has on inferred relationships, on the number of MPTs, and on CI' and RI as measures of homoplasy. Bootstrap leaf stabilities provide a proxy for the contribution of individual taxa to the branch support in the rest of the tree. There is no significant difference in the impact of fossil versus extant taxa on relationships, numbers of MPTs, and CI' or RI. However, adding individual fossil taxa is more likely to reduce the total branch support of the tree than adding extant taxa. This must be weighed against the superior taxon sampling afforded by including judiciously coded fossils, providing data from otherwise unsampled regions of the tree. We therefore recommend that investigators should include fossils, in the absence of compelling and case specific reasons for their exclusion.

Profiling phylogenetic informativeness

The resolution of four controversial topics in phylogenetic experimental design hinges upon the informativeness of characters about the historical relationships among taxa. These controversies regard the power of different classes of phylogenetic character, the relative utility of increased taxonomic versus character sampling, the differentiation between lack of phylogenetic signal and a historical rapid radiation, and the design of taxonomically broad phylogenetic studies optimized by taxonomically sparse genome-scale data. Quantification of the informativeness of characters for resolution of phylogenetic hypotheses during specified historical epochs is key to the resolution of these controversies. Here, such a measure of phylogenetic informativeness is formulated. The optimal rate of evolution of a character to resolve a dated four-taxon polytomy is derived. By scaling the asymptotic informativeness of a character evolving at a nonoptimal rate by the derived asymptotic optimum, and by normalizing so that net phylogenetic informativeness is equivalent for all rates when integrated across all of history, an informativeness profile across history is derived. Calculation of the informativeness per base pair allows estimation of the cost-effectiveness of character sampling. Calculation of the informativeness per million years allows comparison across historical radiations of the utility of a gene for the inference of rapid adaptive radiation. The theory is applied to profile the phylogenetic informativeness of the genes BRCA1, RAG1, GHR, and c-myc from a muroid rodent sequence data set. Bounded integrations of the phylogenetic profile of these genes over four epochs comprising the diversifications of the muroid rodents, the mammals, the lobe-limbed vertebrates, and the early metazoans demonstrate the differential power of these genes to resolve the branching order among ancestral lineages. This measure of phylogenetic informativeness yields a new kind of information for evaluation of phylogenetic experiments. It conveys the utility of the addition of characters a phylogenetic study and it provides a basis for deciding whether appropriate phylogenetic power has been applied to a polytomy that is proposed to be a rapid radiation. Moreover, it provides a quantitative measure of the capacity of a gene to resolve soft polytomies.

A post-synaptic scaffold at the origin of the animal kingdom

BACKGROUND

The evolution of complex sub-cellular structures such as the synapse requires the assembly of multiple proteins, each conferring added functionality to the integrated structure. Tracking the early evolution of synapses has not been possible without genomic information from the earliest branching animals. As the closest extant relatives to the Eumetazoa, Porifera (sponges) represent a pivotal group for understanding the evolution of nervous systems, because sponges lack neurons with clearly recognizable synapses, in contrast to eumetazoan animals.

METHODOLOGY/PRINCIPAL FINDINGS

We show that the genome of the demosponge Amphimedon queenslandica possesses a nearly complete set of post-synaptic protein homologs whose conserved interaction motifs suggest assembly into a complex structure. In the critical synaptic scaffold gene, dlg, residues that make hydrogen bonds and van der Waals interactions with the PDZ ligand are 100% conserved between sponge and human, as is the motif organization of the scaffolds. Expression in Amphimedon of multiple post-synaptic gene homologs in larval flask cells further supports the existence of an assembled structure. Among the few post-synaptic genes absent from Amphimedon, but present in Eumetazoa, are receptor genes including the entire ionotropic glutamate receptor family.

CONCLUSIONS/SIGNIFICANCE

Highly conserved protein interaction motifs and co-expression in sponges of multiple proteins whose homologs interact in eumetazoan synapses indicate that a complex protein scaffold was present at the origin of animals, perhaps predating nervous systems. A relatively small number of crucial innovations to this pre-existing structure may represent the founding changes that led to a post-synaptic element.

Mitochondrial Genome of Savalia savaglia (Cnidaria, Hexacorallia) and Early Metazoan Phylogeny

Mitochondrial genomes have recently become widely used in animal phylogeny, mainly to infer the relationships between vertebrates and other bilaterians. However, only 11 of 723 complete mitochondrial genomes available in the public databases are of early metazoans, including cnidarians (Anthozoa, mainly Scleractinia) and sponges. Although some cnidarians (Medusozoa) are known to possess atypical linear mitochondrial DNA, the anthozoan mitochondrial genome is circular and its organization is similar to that of other metazoans. Because the phylogenetic relationships among Anthozoa as well as their relation to other early metazoans still need to be clarified, we tested whether sequencing the complete mitochondrial genome of Savalia savaglia, an anthozoan belonging to the order Zoantharia (=Zoanthidea), could be useful to infer such relationships. Compared to other anthozoans, S. savaglia's genome is unusually long (20,766 bp) due to the presence of several noncoding intergenic regions (3691 bp). The genome contains all 13 protein coding genes commonly found in metazoans, but like other Anthozoa it lacks most of the tRNAs. Phylogenetic analyses of S. savaglia mitochondrial sequences show Zoantharia branching closely to other Hexacorallia, either as a sister group to Actiniaria or as a sister group to Actiniaria and Scleractinia. The close relationships suggested between Zoantharia and Actiniaria are reinforced by strong similarities in their gene order and the presence of similar introns in the COI and ND5 genes. Our study suggests that mitochondrial genomes can be a source of potentially valuable information on the phylogeny of Hexacorallia and may provide new insights into the evolution of early metazoans.