Improving animal phylogenies with genomic data

Evolution of calcitonin/calcitonin gene-related peptide family in chordates: Identification of CT/CGRP family peptides in cartilaginous fish genome

The calcitonin (CT)/CT gene-related peptide (CGRP) family is a peptide gene family that is widely found in bilaterians. CT, CGRP, adrenomedullin (AM), amylin (AMY), and CT receptor-stimulating peptide (CRSP) are members of the CT/CGRP family. In mammals, CT is involved in calcium homeostasis, while CGRP and AM primarily function in vasodilation. AMY and CRSP are associated with anorectic effects. Diversification of the molecular features and physiological functions of the CT/CGRP family in vertebrate lineages have been extensively reported. However, the origin and diversification mechanisms of the vertebrate CT/CGRP family of peptides remain unclear. In this review, the molecular characteristics of CT/CGRP family peptides and their receptors, along with their major physiological functions in mammals and teleosts, are introduced. Furthermore, novel candidates of the CT/CGRP family in cartilaginous fish are presented based on genomic information. The CT/CGRP family peptides and receptors in urochordates and cephalochordates, which are closely related to vertebrates, are also described. Finally, a putative evolutionary scenario of the CT/CGRP family peptides and receptors in chordates is discussed.

G protein-coupled receptors (GPCRs) in rotifers and cladocerans: Potential applications in ecotoxicology, ecophysiology, comparative endocrinology, and pharmacology

The G protein-coupled receptor (GPCR) superfamily plays a fundamental role in both sensory functions and the regulation of homeostasis, and is highly conserved across the eukaryote taxa. Its functional diversity is related to a conserved seven-transmembrane core and invariant set of intracellular signaling mechanisms. The interplay between these properties is key to the evolutionary success of GPCR. As this superfamily originated from a common ancestor, GPCR genes have evolved via lineage-specific duplications through the process of adaptation. Here we summarized information on GPCR gene families in rotifers and cladocerans based on their evolutionary position in aquatic invertebrates and their potential application in ecotoxicology, ecophysiology, comparative endocrinology, and pharmacology. Phylogenetic analyses were conducted to examine the evolutionary significance of GPCR gene families and to provide structural insight on their role in aquatic invertebrates. In particular, most GPCR gene families have undergone sporadic evolutionary processes, but some GPCRs are highly conserved across species despite the dynamics of GPCR evolution. Overall, this review provides a better understanding of GPCR evolution in aquatic invertebrates and expand our knowledge of the potential application of these receptors in various fields.

Snake Venomics: Fundamentals, Recent Updates, and a Look to the Next Decade

Venomic research, powered by techniques adapted from proteomics, transcriptomics, and genomics, seeks to unravel the diversity and complexity of venom through which knowledge can be applied in the treatment of envenoming, biodiscovery, and conservation. Snake venom proteomics is most extensively studied, but the methods varied widely, creating a massive amount of information which complicates data comparison and interpretation. Advancement in mass spectrometry technology, accompanied by growing databases and sophisticated bioinformatic tools, has overcome earlier limitations of protein identification. The progress, however, remains challenged by limited accessibility to samples, non-standardized quantitative methods, and biased interpretation of -omic data. Next-generation sequencing (NGS) technologies enable high-throughput venom-gland transcriptomics and genomics, complementing venom proteomics by providing deeper insights into the structural diversity, differential expression, regulation and functional interaction of the toxin genes. Venomic tissue sampling is, however, difficult due to strict regulations on wildlife use and transfer of biological materials in some countries. Limited resources for techniques and funding are among other pertinent issues that impede the progress of venomics, particularly in less developed regions and for neglected species. Genuine collaboration between international researchers, due recognition of regional experts by global organizations (e.g., WHO), and improved distribution of research support, should be embraced.

OUP accepted manuscript

Homology of highly divergent genes often cannot be determined from sequence similarity alone. For example, we recently identified in the aphid Hormaphis cornu a family of rapidly evolving bicycle genes, which encode novel proteins implicated as plant gall effectors, and sequence similarity search methods yielded few putative bicycle homologs in other species. Coding sequence-independent features of genes, such as intron-exon boundaries, often evolve more slowly than coding sequences, however, and can provide complementary evidence for homology. We found that a linear logistic regression classifier using only structural features of bicycle genes identified many putative bicycle homologs in other species. Independent evidence from sequence features and intron locations supported homology assignments. To test the potential roles of bicycle genes in other aphids, we sequenced the genome of a second gall-forming aphid, Tetraneura nigriabdominalis and found that many bicycle genes are strongly expressed in the salivary glands of the gall forming foundress. In addition, bicycle genes are strongly overexpressed in the salivary glands of a non-gall forming aphid, Acyrthosiphon pisum, and in the non-gall forming generations of H. cornu. These observations suggest that Bicycle proteins may be used by multiple aphid species to manipulate plants in diverse ways. Incorporation of gene structural features into sequence search algorithms may aid identification of deeply divergent homologs, especially of rapidly evolving genes involved in host-parasite interactions.

The community-curated Pristionchus pacificus genome facilitates automated gene annotation improvement in related nematodes

Background

The nematode Pristionchus pacificus is an established model organism for comparative studies with Caenorhabditis elegans. Over the past years, it developed into an independent animal model organism for elucidating the genetic basis of phenotypic plasticity. Community-based curations were employed recently to improve the quality of gene annotations of P. pacificus and to more easily facilitate reverse genetic studies using candidate genes from C. elegans.

Results

Here, I demonstrate that the reannotation of phylogenomic data from nine related nematode species using the community-curated P. pacificus gene set as homology data substantially improves the quality of gene annotations. Benchmarking of universal single copy orthologs (BUSCO) estimates a median completeness of 84% which corresponds to a 9% increase over previous annotations. Nevertheless, the ability to infer gene models based on homology already drops beyond the genus level reflecting the rapid evolution of nematode lineages. This also indicates that the highly curated C. elegans genome is not optimally suited for annotating non-Caenorhabditis genomes based on homology. Furthermore, comparative genomic analysis of apparently missing BUSCO genes indicates a failure of ortholog detection by the BUSCO pipeline due to the insufficient sample size and phylogenetic breadth of the underlying OrthoDB data set. As a consequence, the quality of multiple divergent nematode genomes might be underestimated.

Conclusions

This study highlights the need for optimizing gene annotation protocols and it demonstrates the benefit of a high quality genome for phylogenomic data of related species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07529-x.

Rethinking the Role of the Nervous System: Lessons From the Hydra Holobiont

Here we evaluate our current understanding of the function of the nervous system in Hydra, a non-bilaterian animal which is among the first metazoans that contain neurons. We highlight growing evidence that the nervous system, with its rich repertoire of neuropeptides, is involved in controlling resident beneficial microbes. We also review observations that indicate that microbes affect the animal's behavior by directly interfering with neuronal receptors. These findings provide new insight into the original role of the nervous system, and suggest that it emerged to orchestrate multiple functions including host-microbiome interactions. The excitement of future research in the Hydra model now relies on uncovering the common rules and principles that govern the interaction between neurons and microbes and the extent to which such laws might apply to other and more complex organisms.

Orthonectids Are Highly Degenerate Annelid Worms

The animal groups of Orthonectida and Dicyemida are tiny, extremely simple, vermiform endoparasites of various marine animals and have been linked in the Mesozoa (Figure 1). The Orthonectida (Figures 1A and 1B) have a few hundred cells, including a nervous system of just ten cells [2], and the Dicyemida (Figure 1C) are even simpler, with ∼40 cells [3]. They are classic "Problematica" [4]-the name Mesozoa suggests an evolutionary position intermediate between Protozoa and Metazoa (animals) [5] and implies that their simplicity is a primitive state, but molecular data have shown they are members of Lophotrochozoa within Bilateria [6-9], which means that they derive from a more complex ancestor. Their precise affinities remain uncertain, however, and it is disputed whether they even constitute a clade. Ascertaining their affinities is complicated by the very fast evolution observed in their genes, potentially leading to the common systematic error of long-branch attraction (LBA) [10]. Here, we use mitochondrial and nuclear gene sequence data and show that both dicyemids and orthonectids are members of the Lophotrochozoa. Carefully addressing the effects of unequal rates of evolution, we show that the Mesozoa is polyphyletic. While the precise position of dicyemids remains unresolved within Lophotrochozoa, we identify orthonectids as members of the phylum Annelida. This result reveals one of the most extreme cases of body-plan simplification in the animal kingdom; our finding makes sense of an annelid-like cuticle in orthonectids [2] and suggests that the circular muscle cells repeated along their body [11] may be segmental in origin.

A comprehensive and quantitative exploration of thousands of viral genomes

The complete assembly of viral genomes from metagenomic datasets (short genomic sequences gathered from environmental samples) has proven to be challenging, so there are significant blind spots when we view viral genomes through the lens of metagenomics. One approach to overcoming this problem is to leverage the thousands of complete viral genomes that are publicly available. Here we describe our efforts to assemble a comprehensive resource that provides a quantitative snapshot of viral genomic trends – such as gene density, noncoding percentage, and abundances of functional gene categories – across thousands of viral genomes. We have also developed a coarse-grained method for visualizing viral genome organization for hundreds of genomes at once, and have explored the extent of the overlap between bacterial and bacteriophage gene pools. Existing viral classification systems were developed prior to the sequencing era, so we present our analysis in a way that allows us to assess the utility of the different classification systems for capturing genomic trends.

Zoology: War of the Worms

The phylogenetic affinities of Xenacoelomorpha - the phylum comprising Xenoturbella bocki and acoelomorph worms - are debated. Two recent studies conclude they represent the earliest branching bilaterally symmetrical animals, but additional tests may be needed to confirm this notion.

The Calcitonin/Calcitonin Gene-Related Peptide Family in Invertebrate Deuterostomes

Calcitonin (CT)/CT gene-related peptide (CGRP) family peptides (CT/CGRP family peptides) including CT, CGRP, adrenomedullin, amylin, and CT receptor-stimulating peptide have been identified from various vertebrates and perform a variety of important physiological functions. These peptides bind to two types of receptors including CT receptor (CTR) and CTR-like receptor (CLR). Receptor recognition of CT/CGRP family peptides is determined by the heterodimer between CTR/CLR and receptor activity-modifying protein (RAMP). Comparative studies of the CT/CGRP family have been exclusively performed in vertebrates from teleost fishes to mammals and strongly manifest that the CGRP family system containing peptides, their receptors, and RAMPs was derived from a common ancestor. In addition, CT/CGRP family peptides and their receptors are also identified and inferred from various invertebrate species. However, the evolutionary process of the CT/CGRP family from invertebrates to vertebrates remains enigmatic. In this review, I principally summarize the CT/CGRP family peptides and their receptors in invertebrate deuterostomes, highlighting the study of invertebrate chordates including ascidians and amphioxi. The CT/CGRP family peptide that shows similar molecular structure and function with that of vertebrate CT has been identified from ascidian, Ciona intestinalis. Amphioxus, Branchiostoma floridae also possessed three CT/CGRP family peptides, one CTR/CLR receptor, and three RAMP-like proteins. The molecular function of the receptor complex formed by amphioxus CTR/CLR and a RAMP-like protein was clarified. Moreover, CT/CGRP family peptides have been identified in the superphylum Ambulacraria, which is close to Chordata. Finally, this review provides potential hypotheses of the evolution of CGRP family peptides and their receptors from invertebrates to vertebrates.

New Insights Into the Roles of Retinoic Acid Signaling in Nervous System Development and the Establishment of Neurotransmitter Systems

Secreted chiefly from the underlying mesoderm, the morphogen retinoic acid (RA) is well known to contribute to the specification, patterning, and differentiation of neural progenitors in the developing vertebrate nervous system. Furthermore, RA influences the subtype identity and neurotransmitter phenotype of subsets of maturing neurons, although relatively little is known about how these functions are mediated. This review provides a comprehensive overview of the roles played by RA signaling during the formation of the central and peripheral nervous systems of vertebrates and highlights its effects on the differentiation of several neurotransmitter systems. In addition, the evolutionary history of the RA signaling system is discussed, revealing both conserved properties and alternate modes of RA action. It is proposed that comparative approaches should be employed systematically to expand our knowledge of the context-dependent cellular mechanisms controlled by the multifunctional signaling molecule RA.

Snake Genome Sequencing: Results and Future Prospects

Snake genome sequencing is in its infancy—very much behind the progress made in sequencing the genomes of humans, model organisms and pathogens relevant to biomedical research, and agricultural species. We provide here an overview of some of the snake genome projects in progress, and discuss the biological findings, with special emphasis on toxinology, from the small number of draft snake genomes already published. We discuss the future of snake genomics, pointing out that new sequencing technologies will help overcome the problem of repetitive sequences in assembling snake genomes. Genome sequences are also likely to be valuable in examining the clustering of toxin genes on the chromosomes, in designing recombinant antivenoms and in studying the epigenetic regulation of toxin gene expression.

Refuting the hypothesis that the acquisition of germ plasm accelerates animal evolution

Primordial germ cells (PGCs) give rise to the germ line in animals. PGCs are specified during embryogenesis either by an ancestral mechanism of cell-cell signalling (induction) or by a derived mechanism of maternally provided germ plasm (preformation). Recently, a hypothesis was set forth purporting that germ plasm liberates selective constraint and accelerates an organism's protein sequence evolution, especially for genes from early developmental stages, thereby leading to animal species radiations; empirical validation has been claimed in vertebrates. Here we present findings from global rates of protein evolution in vertebrates and invertebrates refuting this hypothesis. Contrary to assertions of the hypothesis, we find no effect of preformation on protein sequence evolution, the evolutionary rates of early-stage developmental genes, or on species diversification. We conclude that the hypothesis is mechanistically implausible, and our multi-faceted analysis shows no empirical support for any of its predictions.

Phylogenomic Insights into Animal Evolution

Animals make up only a small fraction of the eukaryotic tree of life, yet, from our vantage point as members of the animal kingdom, the evolution of the bewildering diversity of animal forms is endlessly fascinating. In the century following the publication of Darwin's Origin of Species, hypotheses regarding the evolution of the major branches of the animal kingdom - their relationships to each other and the evolution of their body plans - was based on a consideration of the morphological and developmental characteristics of the different animal groups. This morphology-based approach had many successes but important aspects of the evolutionary tree remained disputed. In the past three decades, molecular data, most obviously primary sequences of DNA and proteins, have provided an estimate of animal phylogeny largely independent of the morphological evolution we would ultimately like to understand. The molecular tree that has evolved over the past three decades has drastically altered our view of animal phylogeny and many aspects of the tree are no longer contentious. The focus of molecular studies on relationships between animal groups means, however, that the discipline has become somewhat divorced from the underlying biology and from the morphological characteristics whose evolution we aim to understand. Here, we consider what we currently know of animal phylogeny; what aspects we are still uncertain about and what our improved understanding of animal phylogeny can tell us about the evolution of the great diversity of animal life.

A New Orthology Assessment Method for Phylogenomic Data: Unrooted Phylogenetic Orthology

Current sequencing technologies are making available unprecedented amounts of genetic data for a large variety of species including nonmodel organisms. Although many phylogenomic surveys spend considerable time finding orthologs from the wealth of sequence data, these results do not transcend the original study and after being processed for specific phylogenetic purposes these orthologs do not become stable orthology hypotheses. We describe a procedure to detect and document the phylogenetic distribution of orthologs allowing researchers to use this information to guide selection of loci best suited to test specific evolutionary questions. At the core of this pipeline is a new phylogenetic orthology method that is neither affected by the position of the root nor requires explicit assignment of outgroups. We discuss the properties of this new orthology assessment method and exemplify its utility for phylogenomics using a small insects dataset. In addition, we exemplify the pipeline to identify and document stable orthologs for the group of orb-weaving spiders (Araneoidea) using RNAseq data. The scripts used in this study, along with sample files and additional documentation, are available at https://github.com/ballesterus/UPhO.

Histories of molecules: Reconciling the past

Molecular data and methods have become centrally important to evolutionary analysis, largely because they have enabled global phylogenetic reconstructions of the relationships between organisms in the tree of life. Often, however, molecular stories conflict dramatically with morphology-based histories of lineages. The evolutionary origin of animal groups provides one such case. In other instances, different molecular analyses have so far proved irreconcilable. The ancient and major divergence of eukaryotes from prokaryotic ancestors is an example of this sort of problem. Efforts to overcome these conflicts highlight the role models play in phylogenetic reconstruction. One crucial model is the molecular clock; another is that of 'simple-to-complex' modification. I will examine animal and eukaryote evolution against a backdrop of increasing methodological sophistication in molecular phylogeny, and conclude with some reflections on the nature of historical science in the molecular era of phylogeny.

Phylogenetic Relationships of Five Asian Schilbid Genera Including Clupisoma (Siluriformes: Schilbeidae)

The phylogenetic relationships of Asian schilbid catfishes of the genera Clupisoma, Ailia, Horabagrus, Laides and Pseudeutropius are poorly understood, especially those of Clupisoma. Herein, we reconstruct the phylogeny of 38 species of catfishes belonging to 28 genera and 14 families using the concatenated mitochondrial genes COI, cytb, and 16S rRNA, as well as the nuclear genes RAG1 and RAG2. The resulting phylogenetic trees consistently place Clupisoma as the sister taxon of Laides, and the five representative Asian schilbid genera form two monophyletic groups with the relationships (Ailia (Laides, Clupisoma)) and (Horabagrus, Pseudeutropius). The so-called "Big Asia" lineage relates distantly to African schilbids. Independent analyses of the mitochondrial and nuclear DNA data yield differing trees for the two Asian schilbid groups. Analyses of the mitochondrial gene data support a sister-group relationship for (Ailia (Laides, Clupisoma)) and the Sisoroidea and a sister-taxon association of (Horabagrus, Pseudeutropius) and the Bagridae. In contrast, analyses of the combined nuclear data indicate (Ailia (Laides, Clupisoma)) to be the sister group to (Horabagrus, Pseudeutropius). Our results indicate that the Horabagridae, recognized by some authors as consisting of Horabagrus, Pseudeutropius and Clupisoma does not include the latter genus. We formally erect a new family, Ailiidae fam. nov. for a monophyletic Asian group comprised of the genera Ailia, Laides and Clupisoma.

Evolution of Prdm Genes in Animals: Insights from Comparative Genomics

Prdm genes encode transcription factors with a subtype of SET domain known as the PRDF1-RIZ (PR) homology domain and a variable number of zinc finger motifs. These genes are involved in a wide variety of functions during animal development. As most Prdm genes have been studied in vertebrates, especially in mice, little is known about the evolution of this gene family. We searched for Prdm genes in the fully sequenced genomes of 93 different species representative of all the main metazoan lineages. A total of 976 Prdm genes were identified in these species. The number of Prdm genes per species ranges from 2 to 19. To better understand how the Prdm gene family has evolved in metazoans, we performed phylogenetic analyses using this large set of identified Prdm genes. These analyses allowed us to define 14 different subfamilies of Prdm genes and to establish, through ancestral state reconstruction, that 11 of them are ancestral to bilaterian animals. Three additional subfamilies were acquired during early vertebrate evolution (Prdm5, Prdm11, and Prdm17). Several gene duplication and gene loss events were identified and mapped onto the metazoan phylogenetic tree. By studying a large number of nonmetazoan genomes, we confirmed that Prdm genes likely constitute a metazoan-specific gene family. Our data also suggest that Prdm genes originated before the diversification of animals through the association of a single ancestral SET domain encoding gene with one or several zinc finger encoding genes.

The evolution of animal genomes

Genome sequences are now available for hundreds of species sampled across the animal phylogeny, bringing key features of animal genome evolution into sharper focus. The field of animal evolutionary genomics has focused on identifying and classifying the diversity genomic features, reconstructing the history of evolutionary changes in animal genomes, and testing hypotheses about the evolutionary relationships of animals. The grand challenges moving forward are to connect evolutionary changes in genomes with particular evolutionary changes in phenotypes, and to determine which changes are driven by selection. This will require far greater genome sampling both across and within species, extensive phenotype data, a well resolved animal phylogeny, and advances in comparative methods.

Harnessing cross-species alignment to discover SNPs and generate a draft genome sequence of a bighorn sheep (Ovis canadensis)

Background

Whole genome sequences (WGS) have proliferated as sequencing technology continues to improve and costs decline. While many WGS of model or domestic organisms have been produced, a growing number of non-model species are also being sequenced. In the absence of a reference, construction of a genome sequence necessitates de novo assembly which may be beyond the ability of many labs due to the large volumes of raw sequence data and extensive bioinformatics required. In contrast, the presence of a reference WGS allows for alignment which is more tractable than assembly. Recent work has highlighted that the reference need not come from the same species, potentially enabling a wide array of species WGS to be constructed using cross-species alignment. Here we report on the creation a draft WGS from a single bighorn sheep (Ovis canadensis) using alignment to the closely related domestic sheep (Ovis aries).

Results

Two sequencing libraries on SOLiD platforms yielded over 865 million reads, and combined alignment to the domestic sheep reference resulted in a nearly complete sequence (95% coverage of the reference) at an average of 12x read depth (104 SD). From this we discovered over 15 million variants and annotated them relative to the domestic sheep reference. We then conducted an enrichment analysis of those SNPs showing fixed differences between the reference and sequenced individual and found significant differences in a number of gene ontology (GO) terms, including those associated with reproduction, muscle properties, and bone deposition.

Conclusion

Our results demonstrate that cross-species alignment enables the creation of novel WGS for non-model organisms. The bighorn sheep WGS will provide a resource for future resequencing studies or comparative genomics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1618-x) contains supplementary material, which is available to authorized users.

Phylogenomic analysis of echinoderm class relationships supports Asterozoa

While some aspects of the phylogeny of the five living echinoderm classes are clear, the position of the ophiuroids (brittlestars) relative to asteroids (starfish), echinoids (sea urchins) and holothurians (sea cucumbers) is controversial. Ophiuroids have a pluteus-type larva in common with echinoids giving some support to an ophiuroid/echinoid/holothurian clade named Cryptosyringida. Most molecular phylogenetic studies, however, support an ophiuroid/asteroid clade (Asterozoa) implying either convergent evolution of the pluteus or reversals to an auricularia-type larva in asteroids and holothurians. A recent study of 10 genes from four of the five echinoderm classes used 'phylogenetic signal dissection' to separate alignment positions into subsets of (i) suboptimal, heterogeneously evolving sites (invariant plus rapidly changing) and (ii) the remaining optimal, homogeneously evolving sites. Along with most previous molecular phylogenetic studies, their set of heterogeneous sites, expected to be more prone to systematic error, support Asterozoa. The homogeneous sites, in contrast, support an ophiuroid/echinoid grouping, consistent with the cryptosyringid clade, leading them to posit homology of the ophiopluteus and echinopluteus. Our new dataset comprises 219 genes from all echinoderm classes; analyses using probabilistic Bayesian phylogenetic methods strongly support Asterozoa. The most reliable, slowly evolving quartile of genes also gives highest support for Asterozoa; this support diminishes in second and third quartiles and the fastest changing quartile places the ophiuroids close to the root. Using phylogenetic signal dissection, we find heterogenous sites support an unlikely grouping of Ophiuroidea + Holothuria while homogeneous sites again strongly support Asterozoa. Our large and taxonomically complete dataset finds no support for the cryptosyringid hypothesis; in showing strong support for the Asterozoa, our preferred topology leaves the question of homology of pluteus larvae open.

Genomic repeat abundances contain phylogenetic signal

A large proportion of genomic information, particularly repetitive elements, is usually ignored when researchers are using next-generation sequencing. Here we demonstrate the usefulness of this repetitive fraction in phylogenetic analyses, utilizing comparative graph-based clustering of next-generation sequence reads, which results in abundance estimates of different classes of genomic repeats. Phylogenetic trees are then inferred based on the genome-wide abundance of different repeat types treated as continuously varying characters; such repeats are scattered across chromosomes and in angiosperms can constitute a majority of nuclear genomic DNA. In six diverse examples, five angiosperms and one insect, this method provides generally well-supported relationships at interspecific and intergeneric levels that agree with results from more standard phylogenetic analyses of commonly used markers. We propose that this methodology may prove especially useful in groups where there is little genetic differentiation in standard phylogenetic markers. At the same time as providing data for phylogenetic inference, this method additionally yields a wealth of data for comparative studies of genome evolution.

A transcriptome approach to ecdysozoan phylogeny

The monophyly of Ecdysozoa, which comprise molting phyla, has received strong support from several lines of evidence. However, the internal relationships of Ecdysozoa are still contended. We generated expressed sequence tags from a priapulid (penis worm), a kinorhynch (mud dragon), a tardigrade (water bear) and five chelicerate taxa by 454 transcriptome sequencing. A multigene alignment was assembled from 63 taxa, which comprised after matrix optimization 24,249 amino acid positions with high data density (2.6% gaps, 19.1% missing data). Phylogenetic analyses employing various models support the monophyly of Ecdysozoa. A clade combining Priapulida and Kinorhyncha (i.e. Scalidophora) was recovered as the earliest branch among Ecdysozoa. We conclude that Cycloneuralia, a taxon erected to combine Priapulida, Kinorhyncha and Nematoda (and others), are paraphyletic. Rather Arthropoda (including Onychophora) are allied with Nematoda and Tardigrada. Within Arthropoda, we found strong support for most clades, including monophyletic Mandibulata and Pancrustacea. The phylogeny within the Euchelicerata remained largely unresolved. There is conflicting evidence on the position of tardigrades: While Bayesian and maximum likelihood analyses of only slowly evolving genes recovered Tardigrada as a sister group to Arthropoda, analyses of the full data set, and of subsets containing genes evolving at fast and intermediate rates identified a clade of Tardigrada and Nematoda. Notably, the latter topology is also supported by the analyses of indel patterns.

Naturally occurring tumours in the basal metazoan Hydra

The molecular nature of tumours is well studied in vertebrates, although their evolutionary origin remains unknown. In particular, there is no evidence for naturally occurring tumours in pre-bilaterian animals, such as sponges and cnidarians. This is somewhat surprising given that recent computational studies have predicted that most metazoans might be prone to develop tumours. Here we provide first evidence for naturally occurring tumours in two species of Hydra. Histological, cellular and molecular data reveal that these tumours are transplantable and might originate by differentiation arrest of female gametes. Growth of tumour cells is independent from the cellular environment. Tumour-bearing polyps have significantly reduced fitness. In addition, Hydra tumours show a greatly altered transcriptome that mimics expression shifts in vertebrate cancers. Therefore, this study shows that spontaneous tumours have deep evolutionary roots and that early branching animals may be informative in revealing the fundamental mechanisms of tumorigenesis.

Genomic support for a moa-tinamou clade and adaptive morphological convergence in flightless ratites

One of the most startling discoveries in avian molecular phylogenetics is that the volant tinamous are embedded in the flightless ratites, but this topology remains controversial because recent morphological phylogenies place tinamous as the closest relative of a monophyletic ratite clade. Here, we integrate new phylogenomic sequences from 1,448 nuclear DNA loci totaling almost 1 million bp from the extinct little bush moa, Chilean tinamou, and emu with available sequences from ostrich, elegant crested tinamou, four neognaths, and the green anole. Phylogenetic analysis using standard homogeneous models and heterogeneous models robust to common topological artifacts recovered compelling support for ratite paraphyly with the little bush moa closest to tinamous within ratites. Ratite paraphyly was further corroborated by eight independent CR1 retroposon insertions. Analysis of morphological characters reinterpreted on a 27-gene paleognath topology indicates that many characters are convergent in the ratites, probably as the result of adaptation to a cursorial life style.

Selecting optimal partitioning schemes for phylogenomic datasets

BACKGROUND

Partitioning involves estimating independent models of molecular evolution for different subsets of sites in a sequence alignment, and has been shown to improve phylogenetic inference. Current methods for estimating best-fit partitioning schemes, however, are only computationally feasible with datasets of fewer than 100 loci. This is a problem because datasets with thousands of loci are increasingly common in phylogenetics.

METHODS

We develop two novel methods for estimating best-fit partitioning schemes on large phylogenomic datasets: strict and relaxed hierarchical clustering. These methods use information from the underlying data to cluster together similar subsets of sites in an alignment, and build on clustering approaches that have been proposed elsewhere.

RESULTS

We compare the performance of our methods to each other, and to existing methods for selecting partitioning schemes. We demonstrate that while strict hierarchical clustering has the best computational efficiency on very large datasets, relaxed hierarchical clustering provides scalable efficiency and returns dramatically better partitioning schemes as assessed by common criteria such as AICc and BIC scores.

CONCLUSIONS

These two methods provide the best current approaches to inferring partitioning schemes for very large datasets. We provide free open-source implementations of the methods in the PartitionFinder software. We hope that the use of these methods will help to improve the inferences made from large phylogenomic datasets.

Selecting optimal partitioning schemes for phylogenomic datasets

BACKGROUND

Partitioning involves estimating independent models of molecular evolution for different subsets of sites in a sequence alignment, and has been shown to improve phylogenetic inference. Current methods for estimating best-fit partitioning schemes, however, are only computationally feasible with datasets of fewer than 100 loci. This is a problem because datasets with thousands of loci are increasingly common in phylogenetics.

METHODS

We develop two novel methods for estimating best-fit partitioning schemes on large phylogenomic datasets: strict and relaxed hierarchical clustering. These methods use information from the underlying data to cluster together similar subsets of sites in an alignment, and build on clustering approaches that have been proposed elsewhere.

RESULTS

We compare the performance of our methods to each other, and to existing methods for selecting partitioning schemes. We demonstrate that while strict hierarchical clustering has the best computational efficiency on very large datasets, relaxed hierarchical clustering provides scalable efficiency and returns dramatically better partitioning schemes as assessed by common criteria such as AICc and BIC scores.

CONCLUSIONS

These two methods provide the best current approaches to inferring partitioning schemes for very large datasets. We provide free open-source implementations of the methods in the PartitionFinder software. We hope that the use of these methods will help to improve the inferences made from large phylogenomic datasets.

Selecting optimal partitioning schemes for phylogenomic datasets

BACKGROUND

Partitioning involves estimating independent models of molecular evolution for different subsets of sites in a sequence alignment, and has been shown to improve phylogenetic inference. Current methods for estimating best-fit partitioning schemes, however, are only computationally feasible with datasets of fewer than 100 loci. This is a problem because datasets with thousands of loci are increasingly common in phylogenetics.

METHODS

We develop two novel methods for estimating best-fit partitioning schemes on large phylogenomic datasets: strict and relaxed hierarchical clustering. These methods use information from the underlying data to cluster together similar subsets of sites in an alignment, and build on clustering approaches that have been proposed elsewhere.

RESULTS

We compare the performance of our methods to each other, and to existing methods for selecting partitioning schemes. We demonstrate that while strict hierarchical clustering has the best computational efficiency on very large datasets, relaxed hierarchical clustering provides scalable efficiency and returns dramatically better partitioning schemes as assessed by common criteria such as AICc and BIC scores.

CONCLUSIONS

These two methods provide the best current approaches to inferring partitioning schemes for very large datasets. We provide free open-source implementations of the methods in the PartitionFinder software. We hope that the use of these methods will help to improve the inferences made from large phylogenomic datasets.

Ivermectin sensitivity is an ancient trait affecting all ecdysozoa but shows phylogenetic clustering among sepsid flies

Avermectins are potent and popular veterinary pharmaceuticals used globally to fight parasites of livestock and humans. By disturbing ion channel transport through the membrane, avermectins are effective against endo- and ectoparasitic round and horsehair worms (Nematoida), insects, or ticks (Arthropoda), but not against Plathelminthes, including flatworms (Trematoda) and tapeworms (Cestoda), or segmented worms (Annelida). Unfortunately, excreted avermectins have strong nontarget effects on beneficial arthropods such as the insect community decomposing livestock dung, ultimately impeding this important ecosystem function to the extent that regulators mandate standardized eco-toxicological tests of dung organisms worldwide. We show that the ancient phylogenetic pattern and qualitative mechanism of avermectin sensitivity is conserved and compatible with most recent phylogenomic hypotheses grouping the Nematoida with the Arthropoda as Ecdysozoa (molting animals). At the species level, we demonstrate phylogenetic clustering in ivermectin sensitivities of 23 species of sepsid dung flies (Diptera: Sepsidae). This clustered 500-fold quantitative variation in sensitivity may indicate recent lineage-specific responses to selection, but more likely reflects pre-existing genetic variation with pleiotropic effects on eco-toxicological responses to pollutants. Regardless, our results question the common practice in eco-toxicology of choosing single test species to infer detrimental effects on entire species communities, which should ideally assess a representative taxonomic sample.

When integration fails: Prokaryote phylogeny and the tree of life

Much is being written these days about integration, its desirability and even its necessity when complex research problems are to be addressed. Seldom, however, do we hear much about the failure of such efforts. Because integration is an ongoing activity rather than a final achievement, and because today's literature about integration consists mostly of manifesto statements rather than precise descriptions, an examination of unsuccessful integration could be illuminating to understand better how it works. This paper will examine the case of prokaryote phylogeny and its apparent failure to achieve integration within broader tree-of-life accounts of evolutionary history (often called 'universal phylogeny'). Despite the fact that integrated databases exist of molecules pertinent to the phylogenetic reconstruction of all lineages of life, and even though the same methods can be used to construct phylogenies wherever the organisms fall on the tree of life, prokaryote phylogeny remains at best only partly integrated within tree-of-life efforts. I will examine why integration does not occur, compare it with integrative practices in animal and other eukaryote phylogeny, and reflect on whether there might be different expectations of what integration should achieve. Finally, I will draw some general conclusions about integration and its function as a 'meta-heuristic' in the normative commitments guiding scientific practice.

Homology versus convergence in resolving transphyletic correspondences of brain organization

Due to the largely absent fossil record, phylogenetic comparisons of brain structures rely on the analysis of nervous systems in extant taxa, many of which appear to have distinctive and dissimilar neural arrangements. The use of a multitude of comparative criteria, including developmental genetics, phylogenomics and neural circuit architecture, has recently resolved a highly conserved structural and functional ground pattern organization in the arthropod central complex and vertebrate basal ganglia. The minuteness of resemblance is exemplified by orthologous action selection circuits that are formed by homologous gene networks and which can lead to similar pathologies and behavioral disorders. It has been argued, however, that these similarities of brain centers can only be due to convergent evolution. What is still missing is a plausible scenario to explain how convergence could result in such a multitude of similarities and minuteness of resemblances, including gene expression, functional attributes and pathologies. In contrast, homology by common descent is the more parsimonious explanation. Moreover, the divergent elaboration of arthropod central complex and vertebrate basal ganglia does not obscure their shared ground pattern organization and thus genealogical correspondence.

miRNAs: small genes with big potential in metazoan phylogenetics

microRNAs (miRNAs) are a key component of gene regulatory networks and have been implicated in the regulation of virtually every biological process found in multicellular eukaryotes. What makes them interesting from a phylogenetic perspective is the high conservation of primary sequence between taxa, their accrual in metazoan genomes through evolutionary time, and the rarity of secondary loss in most metazoan taxa. Despite these properties, the use of miRNAs as phylogenetic markers has not yet been discussed within a clear conceptual framework. Here we highlight five properties of miRNAs that underlie their utility in phylogenetics: 1) The processes of miRNA biogenesis enable the identification of novel miRNAs without prior knowledge of sequence; 2) The continuous addition of miRNA families to metazoan genomes through evolutionary time; 3) The low level of secondary gene loss in most metazoan taxa; 4) The low substitution rate in the mature miRNA sequence; and 5) The small probability of convergent evolution of two miRNAs. Phylogenetic analyses using both Bayesian and parsimony methods on a eumetazoan miRNA data set highlight the potential of miRNAs to become an invaluable new tool, especially when used as an additional line of evidence, to resolve previously intractable nodes within the tree of life.

How are comparative genomics and the study of microRNAs changing our views on arthropod endocrinology and adaptations to the environment?

As the last few decades of work has shown, precise regulation of biosynthesis and release of arthropod hormones is essential to cope with environmental stresses and challenges. In crustaceans and insects, the sesquiterpenoids methyl farnesoate (MF), farnesoic acid (FA) and juvenile hormone (JH) regulate many developmental, physiological, and reproductive processes. In this review, we discuss how comparative genomics has and will impact our views on arthropod endocrinology. We will also highlight the current knowledge of regulation of genes involved in arthropod hormone biosynthesis by microRNAs, and describe the potential insights into arthropod endocrinology, evolution, and adaptation that are likely to come from the study of microRNAs.

Evolution of signal multiplexing by 14-3-3-binding 2R-ohnologue protein families in the vertebrates

14-3-3 proteins regulate cellular responses to stimuli by docking onto pairs of phosphorylated residues on target proteins. The present study shows that the human 14-3-3-binding phosphoproteome is highly enriched in 2R-ohnologues, which are proteins in families of two to four members that were generated by two rounds of whole genome duplication at the origin of the vertebrates. We identify 2R-ohnologue families whose members share a ‘lynchpin’, defined as a 14-3-3-binding phosphosite that is conserved across members of a given family, and aligns with a Ser/Thr residue in pro-orthologues from the invertebrate chordates. For example, the human receptor expression enhancing protein (REEP) 1–4 family has the commonest type of lynchpin motif in current datasets, with a phosphorylatable serine in the –2 position relative to the 14-3-3-binding phosphosite. In contrast, the second 14-3-3-binding sites of REEPs 1–4 differ and are phosphorylated by different kinases, and hence the REEPs display different affinities for 14-3-3 dimers. We suggest a conceptual model for intracellular regulation involving protein families whose evolution into signal multiplexing systems was facilitated by 14-3-3 dimer binding to lynchpins, which gave freedom for other regulatory sites to evolve. While increased signalling complexity was needed for vertebrate life, these systems also generate vulnerability to genetic disorders.

Field et Al. Redux

On 12 February 1988 (by coincidence Charles Darwin's birthday), a paper published in Science by Katherine Field, Rudy Raff, and colleagues presented the first credible molecular analysis of metazoan phylogeny based on sequences from the small subunit ribosomal RNA gene (SSU). Here I examine the main conclusions reached in this manuscript. I reconstitute their dataset and, by recompiling software available in 1988, I consider how they might have achieved a more accurate tree. I show how three common methods to avoid systematic error - more data, careful taxon sampling and superior models of evolution - overcome the errors that exist in the original paper. This approach illustrates the basis of some of the major advances of the past 25 years resulting in our current understanding of animal phylogeny.

Deep transcriptome-sequencing and proteome analysis of the hydrothermal vent annelid Alvinella pompejana identifies the CvP-bias as a robust measure of eukaryotic thermostability

BACKGROUND

Alvinella pompejana is an annelid worm that inhabits deep-sea hydrothermal vent sites in the Pacific Ocean. Living at a depth of approximately 2500 meters, these worms experience extreme environmental conditions, including high temperature and pressure as well as high levels of sulfide and heavy metals. A. pompejana is one of the most thermotolerant metazoans, making this animal a subject of great interest for studies of eukaryotic thermoadaptation.

RESULTS

In order to complement existing EST resources we performed deep sequencing of the A. pompejana transcriptome. We identified several thousand novel protein-coding transcripts, nearly doubling the sequence data for this annelid. We then performed an extensive survey of previously established prokaryotic thermoadaptation measures to search for global signals of thermoadaptation in A. pompejana in comparison with mesophilic eukaryotes. In an orthologous set of 457 proteins, we found that the best indicator of thermoadaptation was the difference in frequency of charged versus polar residues (CvP-bias), which was highest in A. pompejana. CvP-bias robustly distinguished prokaryotic thermophiles from prokaryotic mesophiles, as well as the thermophilic fungus Chaetomium thermophilum from mesophilic eukaryotes. Experimental values for thermophilic proteins supported higher CvP-bias as a measure of thermal stability when compared to their mesophilic orthologs. Proteome-wide mean CvP-bias also correlated with the body temperatures of homeothermic birds and mammals.

CONCLUSIONS

Our work extends the transcriptome resources for A. pompejana and identifies the CvP-bias as a robust and widely applicable measure of eukaryotic thermoadaptation.

Chloroplast phylogenomics indicates that Ginkgo biloba is sister to cycads

Molecular phylogenetic studies have not yet reached a consensus on the placement of Ginkgoales, which is represented by the only living species, Ginkgo biloba (common name: ginkgo). At least six discrepant placements of ginkgo have been proposed. This study aimed to use the chloroplast phylogenomic approach to examine possible factors that lead to such disagreeing placements. We found the sequence types used in the analyses as the most critical factor in the conflicting placements of ginkgo. In addition, the placement of ginkgo varied in the trees inferred from nucleotide (NU) sequences, which notably depended on breadth of taxon sampling, tree-building methods, codon positions, positions of Gnetopsida (common name: gnetophytes), and including or excluding gnetophytes in data sets. In contrast, the trees inferred from amino acid (AA) sequences congruently supported the monophyly of a ginkgo and Cycadales (common name: cycads) clade, regardless of which factors were examined. Our site-stripping analysis further revealed that the high substitution saturation of NU sequences mainly derived from the third codon positions and contributed to the variable placements of ginkgo. In summary, the factors we surveyed did not affect results inferred from analyses of AA sequences. Congruent topologies in our AA trees give more confidence in supporting the ginkgo-cycad sister-group hypothesis.

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.

Cephalopod genomics: A plan of strategies and organization

The Cephalopod Sequencing Consortium (CephSeq Consortium) was established at a NESCent Catalysis Group Meeting, “Paths to Cephalopod Genomics- Strategies, Choices, Organization,” held in Durham, North Carolina, USA on May 24-27, 2012. Twenty-eight participants representing nine countries (Austria, Australia, China, Denmark, France, Italy, Japan, Spain and the USA) met to address the pressing need for genome sequencing of cephalopod mollusks. This group, drawn from cephalopod biologists, neuroscientists, developmental and evolutionary biologists, materials scientists, bioinformaticians and researchers active in sequencing, assembling and annotating genomes, agreed on a set of cephalopod species of particular importance for initial sequencing and developed strategies and an organization (CephSeq Consortium) to promote this sequencing. The conclusions and recommendations of this meeting are described in this white paper.

Impact of missing data on phylogenies inferred from empirical phylogenomic data sets

Progress in sequencing technology allows researchers to assemble ever-larger supermatrices for phylogenomic inference. However, current phylogenomic studies often rest on patchy data sets, with some having 80% missing (or ambiguous) data or more. Though early simulations had suggested that missing data per se do not harm phylogenetic inference when using sufficiently large data sets, Lemmon et al. (Lemmon AR, Brown JM, Stanger-Hall K, Lemmon EM. 2009. The effect of ambiguous data on phylogenetic estimates obtained by maximum likelihood and Bayesian inference. Syst Biol. 58:130-145.) have recently cast doubt on this consensus in a study based on the introduction of parsimony-uninformative incomplete characters. In this work, we empirically reassess the issue of missing data in phylogenomics while exploring possible interactions with the model of sequence evolution. First, we note that parsimony-uninformative incomplete characters are actually informative in a probabilistic framework. A reanalysis of Lemmon's data set with this in mind gives a very different interpretation of their results and shows that some of their conclusions may be unfounded. Second, we investigate the effect of the progressive introduction of missing data in a complete supermatrix (126 genes × 39 species) capable of resolving animal relationships. These analyses demonstrate that missing data perturb phylogenetic inference slightly beyond the expected decrease in resolving power. In particular, they exacerbate systematic errors by reducing the number of species effectively available for the detection of multiple substitutions. Consequently, large sparse supermatrices are more sensitive to phylogenetic artifacts than smaller but less incomplete data sets, which argue for experimental designs aimed at collecting a modest number (~50) of highly covered genes. Our results further confirm that including incomplete yet short-branch taxa (i.e., slowly evolving species or close outgroups) can help to eschew artifacts, as predicted by simulations. Finally, it appears that selecting an adequate model of sequence evolution (e.g., the site-heterogeneous CAT model instead of the site-homogeneous WAG model) is more beneficial to phylogenetic accuracy than reducing the level of missing data.

The role of androdioecy and gynodioecy in mediating evolutionary transitions between dioecy and hermaphroditism in the animalia

Dioecy (gonochorism) is dominant within the Animalia, although a recent review suggests hermaphroditism is also common. Evolutionary transitions from dioecy to hermaphroditism (or vice versa) have occurred frequently in animals, but few studies suggest the advantage of such transitions. In particular, few studies assess how hermaphroditism evolves from dioecy or whether androdioecy or gynodioecy should be an "intermediate" stage, as noted in plants. Herein, these transitions are assessed by documenting the numbers of androdioecious and gynodioecious animals and inferring their ancestral reproductive mode. Both systems are rare, but androdioecy was an order of magnitude more common than gynodioecy. Transitions from dioecious ancestors were commonly to androdioecy rather than gynodioecy. Hermaphrodites evolving from sexually dimorphic dioecious ancestors appear to be constrained to those with female-biased sex allocation; such hermaphrodites replace females to coexist with males. Hermaphrodites evolving from sexually monomorphic dioecious ancestors were not similarly constrained. Species transitioning from hermaphroditic ancestors were more commonly androdioecious than gynodioecious, contrasting with similar transitions in plants. In animals, such transitions were associated with size specialization between the sexes, whereas in plants these transitions were to avoid inbreeding depression. Further research should frame these reproductive transitions in a theoretical context, similar to botanical studies.