Phylogenomic Insights into Animal Evolution

Doublecortin-like kinase is required for cnidocyte development in Nematostella vectensis

The complex morphology of neurons requires precise control of their microtubule cytoskeleton. This is achieved by microtubule-associated proteins (MAPs) that regulate the assembly and stability of microtubules, and transport of molecules and vesicles along them. While many of these MAPs function in all cells, some are specifically or predominantly involved in regulating microtubules in neurons. Here we use the sea anemone Nematostella vectensis as a model organism to provide new insights into the early evolution of neural microtubule regulation. As a cnidarian, Nematostella belongs to an outgroup to all bilaterians and thus occupies an informative phylogenetic position for reconstructing the evolution of nervous system development. We identified an ortholog of the microtubule-binding protein doublecortin-like kinase (NvDclk1) as a gene that is predominantly expressed in neurons and cnidocytes (stinging cells), two classes of cells belonging to the neural lineage in cnidarians. A transgenic NvDclk1 reporter line revealed an elaborate network of neurite-like processes emerging from cnidocytes in the tentacles and the body column. A transgene expressing NvDclk1 under the control of the NvDclk1 promoter suggests that NvDclk1 localizes to microtubules and therefore likely functions as a microtubule-binding protein. Further, we generated a mutant for NvDclk1 using CRISPR/Cas9 and show that the mutants fail to generate mature cnidocytes. Our results support the hypothesis that the elaboration of programs for microtubule regulation occurred early in the evolution of nervous systems.

Chemical communication and its role in sexual selection across Animalia

Sexual selection has been studied as a major evolutionary driver of animal diversity for roughly 50 years. Much evidence indicates that competition for mates favors elaborate signaling traits. However, this evidence comes primarily from a few taxa, leaving sexual selection as a salient evolutionary force across Animalia largely untested. Here, we reviewed the evidence for sexual selection on communication across all animal phyla, classes, and orders with emphasis on chemoreception, the only sense shared across lifeforms. An exhaustive literature review documented evidence for sexual selection on chemosensory traits in 10 of 34 animal phyla and indications of sexual selection on chemosensory traits in an additional 13 phyla. Potential targets of sexual selection include structures and processes involved in production, delivery, and detection of chemical signals. Our review suggests sexual selection plays a widespread role in the evolution of communication and highlights the need for research that better reflects animal diversity.

Single-cell transcriptomics refuels the exploration of spiralian biology

Spiralians represent the least studied superclade of bilaterian animals, despite exhibiting the widest diversity of organisms. Although spiralians include iconic organisms, such as octopus, earthworms and clams, a lot remains to be discovered regarding their phylogeny and biology. Here, we review recent attempts to apply single-cell transcriptomics, a new pioneering technology enabling the classification of cell types and the characterisation of their gene expression profiles, to several spiralian taxa. We discuss the methodological challenges and requirements for applying this approach to marine organisms and explore the insights that can be brought by such studies, both from a biomedical and evolutionary perspective. For instance, we show that single-cell sequencing might help solve the riddle of the homology of larval forms across spiralians, but also to better characterise and compare the processes of regeneration across taxa. We highlight the capacity of single-cell to investigate the origin of evolutionary novelties, as the mollusc shell or the cephalopod visual system, but also to interrogate the conservation of the molecular fingerprint of cell types at long evolutionary distances. We hope that single-cell sequencing will open a new window in understanding the biology of spiralians, and help renew the interest for these overlooked but captivating organisms.

DNA Sequences Are as Useful as Protein Sequences for Inferring Deep Phylogenies

Inference of deep phylogenies has almost exclusively used protein rather than DNA sequences based on the perception that protein sequences are less prone to homoplasy and saturation or to issues of compositional heterogeneity than DNA sequences. Here, we analyze a model of codon evolution under an idealized genetic code and demonstrate that those perceptions may be misconceptions. We conduct a simulation study to assess the utility of protein versus DNA sequences for inferring deep phylogenies, with protein-coding data generated under models of heterogeneous substitution processes across sites in the sequence and among lineages on the tree, and then analyzed using nucleotide, amino acid, and codon models. Analysis of DNA sequences under nucleotide-substitution models (possibly with the third codon positions excluded) recovered the correct tree at least as often as analysis of the corresponding protein sequences under modern amino acid models. We also applied the different data-analysis strategies to an empirical dataset to infer the metazoan phylogeny. Our results from both simulated and real data suggest that DNA sequences may be as useful as proteins for inferring deep phylogenies and should not be excluded from such analyses. Analysis of DNA data under nucleotide models has a major computational advantage over protein-data analysis, potentially making it feasible to use advanced models that account for among-site and among-lineage heterogeneity in the nucleotide-substitution process in inference of deep phylogenies.

Amine-recognizing domain in diverse receptors from bacteria and archaea evolved from the universal amino acid sensor

Bacteria possess various receptors that sense different signals and transmit information to enable an optimal adaptation to the environment. A major limitation in microbiology is the lack of information on the signal molecules that activate receptors. Signals recognized by sensor domains are poorly reflected in overall sequence identity, and therefore, the identification of signals from the amino acid sequence of the sensor alone presents a challenge. Biogenic amines are of great physiological importance for microorganisms and humans. They serve as substrates for aerobic and anaerobic growth and play a role of neurotransmitters and osmoprotectants. Here, we report the identification of a sequence motif that is specific for amine-sensing sensor domains that belong to the Cache superfamily of the most abundant extracellular sensors in prokaryotes. We identified approximately 13,000 sensor histidine kinases, chemoreceptors, receptors involved in second messenger homeostasis and Ser/Thr phosphatases from 8,000 bacterial and archaeal species that contain the amine-recognizing motif. The screening of compound libraries and microcalorimetric titrations of selected sensor domains confirmed their ability to specifically bind biogenic amines. Mutants in the amine-binding motif or domains that contain a single mismatch in the binding motif had either no or a largely reduced affinity for amines. We demonstrate that the amine-recognizing domain originated from the universal amino acid-sensing Cache domain, thus providing insight into receptor evolution. Our approach enables precise "wet"-lab experiments to define the function of regulatory systems and therefore holds a strong promise to enable the identification of signals stimulating numerous receptors.

First putative occurrence in the fossil record of choanoflagellates, the sister group of Metazoa

Choanoflagellates are microeukaryotes that inhabit freshwater and marine environments and have long been regarded as the closest living relatives of Metazoa. Knowledge on the evolution of choanoflagellates is key for the understanding of the ancestry of animals, and although molecular clock evidence suggests the appearance of choanoflagellates by late Neoproterozoic, no specimens of choanoflagellates are known to occur in the fossil record. Here the first putative occurrence of choanoflagellates in sediments from the Cretaceous (Cenomanian-Turonian) is described by means of several cutting-edge petrographic techniques, and a discussion of its paleoenvironmental significance is performed. Furthermore, their placement in the organic matter classification systems is argued, with a placement in the Zoomorph Subgroup (Palynomorph Group) of the dispersed organic matter classification system being proposed. Regarding the ICCP System 1994, incorporation of choanoflagellates is, at a first glance, straightforward within the liptinite group, but the definition of a new maceral may be necessary to accommodate the genetic origin of these organisms. While modern choanoflagellates may bring light to the cellular foundations of animal origins, this discovery may provide an older term of comparison to their extant specimens and provide guidelines for possible identification of these organic components in other locations and ages throughout the geological record.

A novel, conserved and possibly functional motif "WHWGHTW" in mitochondrial transcription across Bilateria

The animal mitogenomes which undergone a reductive evolution has an obvious loss of coding capacity compared to their known closest relatives, but it has not yet been fully investigated why and how the intergenic regions do not encode protein and have no known functions, are stably maintained, replicated, and transmitted by the genome. These relatively small intergenic regions may not be under neutral evolution and they may have functional and/or regulatory roles that have yet to be identified. Here, the distribution pattern, sequence content and location of a novel sequence motif of 'WWWGHTW' were bioinformatically investigated and characterised by constructing a sampling mitogenome dataset of 1889 species from 14 phyla representing the clade of Bilateria. This motif is reverse complementary of the previously described DmTTF binding sequence and found in the nd4L- (X) -trnT gene cluster. This cluster commonly exhibits a strand displacement region and an intergenic region among the bilaterian superphylums, particularly in Ecdysozoa. This motif may be accepted as a substrate providing binding sites for the specific interaction with transcription factors because of (i) its reverse complementarity of previously described DmTTF binding sequence, and (ii) the possession of G and T nucleotides in the fourth and sixth positions, (iii) the bias on T and G nucleotides instead of C and A in the degenerated positions. This suggestion is also supported by the presence of a strand displacement region in the nd4L- (X) -trnT gene cluster, particularly in Ecdysozoa consisting of the most rearranged mitogenomes among the bilaterian superphylums.

Evolutionary Divergence of Phosphorylation to Regulate Interactive Protein Networks in Lower and Higher Species

The phosphorylation of proteins affects their functions in extensively documented circumstances. However, the role of phosphorylation in many interactive networks of proteins remains very elusive due to the experimental limits of exploring the transient interaction in a large complex of assembled proteins induced by stimulation. Previous studies have suggested that phosphorylation is a recent evolutionary process that differently regulates ortholog proteins in numerous lineages of living organisms to create new functions. Despite the fact that numerous phospho-proteins have been compared between species, little is known about the organization of the full phospho-proteome, the role of phosphorylation to orchestrate large interactive networks of proteins, and the intertwined phospho-landscape in these networks. In this report, we aimed to investigate the acquired role of phosphate addition in the phenomenon of protein networking in different orders of living organisms. Our data highlighted the acquired status of phosphorylation in organizing large, connected assemblages in Homo sapiens. The protein networking guided by phosphorylation turned out to be prominent in humans, chaotic in yeast, and weak in flies. Furthermore, the molecular functions of GO annotation enrichment regulated by phosphorylation were found to be drastically different between flies, yeast, and humans, suggesting an evolutionary drift specific to each species.

The origin and distribution of the main oxygen sensing mechanism across metazoans

Oxygen sensing mechanisms are essential for metazoans, their origin and evolution in the context of oxygen in Earth history are of interest. To trace the evolution of a main oxygen sensing mechanism among metazoans, the hypoxia induced factor, HIF, we investigated the phylogenetic distribution and phylogeny of 11 of its components across 566 eukaryote genomes. The HIF based oxygen sensing machinery in eukaryotes can be traced as far back as 800 million years (Ma) ago, likely to the last metazoan common ancestor (LMCA), and arose at a time when the atmospheric oxygen content corresponded roughly to the Pasteur point, or roughly 1% of present atmospheric level (PAL). By the time of the Cambrian explosion (541–485 Ma) as oxygen levels started to approach those of the modern atmosphere, the HIF system with its key components HIF1α, HIF1β, PHD1, PHD4, FIH and VHL was well established across metazoan lineages. HIF1α is more widely distributed and therefore may have evolved earlier than HIF2α and HIF3α, and HIF1β and is more widely distributed than HIF2β in invertebrates. PHD1, PHD4, FIH, and VHL appear in all 13 metazoan phyla. The O₂ consuming enzymes of the pathway, PHDs and FIH, have a lower substrate affinity, K_m, for O₂ than terminal oxidases in the mitochondrial respiratory chain, in line with their function as an environmental signal to switch to anaerobic energy metabolic pathways. The ancient HIF system has been conserved and widespread during the period when metazoans evolved and diversified together with O₂ during Earth history.

Law and order in plants - the origin and functional relevance of phyllotaxis

The regular arrangement of organs (phyllotaxis) in vegetative shoots and flowers is one of the most stunning features of plants. Spiral patterns characterized by Fibonacci numbers have attracted the particular interest of natural scientists and mathematicians. Numerous reviews have dealt with the molecular genetic mechanisms underlying phyllotaxis, and modeling studies have sought to recreate phyllotaxis according to mathematical, biochemical, or physical laws. However, what is the functional significance of regular plant architecture, and how did it evolve? We discuss the developmental constraints and selective forces that may have favored the selection of phyllotaxis, and we argue that a central driver of regular phyllotaxis may have been limitations in the allocation of founder cells and metabolic resources to the different tissues in the shoot apex.

Distribution and evolution of the serine/aspartate racemase family in invertebrates. II. Frequent and widespread parallel evolution of aspartate racemase

Our previous studies showed that invertebrate animal serine racemase (SerR) and aspartate racemase (AspR) evolved from a common ancestral gene and are widely distributed. However, the overall molecular evolutionary background of these genes has remained unclear. In the present study we have cloned, expressed and characterized five SerR and three AspR genes from six invertebrate species. The coexistence of SerR and AspR paralogs has been observed in some species, and the presence of both SerR and AspR is here confirmed in the flatworm Macrostomum lignano, the feather star Anneissia japonica, the ark shell Anadara broughtonii and the sea hare Aplysia californica. Comparison of the gene structures revealed the evolution of SerR and AspR. The ancestral species of metazoans probably had a single SerR gene, and the first gene duplication in the common ancestor species of the eumetazoans occurred after the divergence of porifera and eumetazoans, yielding two SerR genes. Most eumetazoans lost one of the two SerR genes, while the echinoderm Anneissia japonica retained both genes. Furthermore, it is clear that invertebrate AspR genes arose through parallel evolution by duplication of the SerR gene followed by substitution of amino acid residues necessary for substrate recognition in multiple lineages.

Renewed perspectives on the sedentary-pelagic last common bilaterian ancestor

Various evaluations of the last common bilaterian ancestor (lcba) currently suggest that it resembled either a microscopic, non-segmented motile adult; or, on the contrary, a complex segmented adult motile urbilaterian. These fundamental inconsistencies remain largely unexplained. A majority of multidisciplinary data regarding sedentary adult ancestral bilaterian organization is overlooked. The sedentary-pelagic model is supported now by a number of novel developmental, paleontological and molecular phylogenetic data: (1) data in support of sedentary sponges, in the adult stage, as sister to all other Metazoa; (2) a similarity of molecular developmental pathways in both adults and larvae across sedentary sponges, cnidarians, and bilaterians; (3) a cnidarian-bilaterian relationship, including a unique sharing of a bona fide Hox-gene cluster, of which the evolutionary appearance does not connect directly to a bilaterian motile organization; (4) the presence of sedentary and tube-dwelling representatives of the main bilaterian clades in the early Cambrian; (5) an absence of definite taxonomic attribution of Ediacaran taxa reconstructed as motile to any true bilaterian phyla; (6) a similarity of tube morphology (and the clear presence of a protoconch-like apical structure of the Ediacaran sedentary Cloudinidae) among shells of the early Cambrian, and later true bilaterians, such as semi-sedentary hyoliths and motile molluscs; (7) recent data that provide growing evidence for a complex urbilaterian, despite a continuous molecular phylogenetic controversy. The present review compares the main existing models and reconciles the sedentary model of an urbilaterian and the model of a larva-like lcba with a unified sedentary(adult)-pelagic(larva) model of the lcba.

Different phylogenomic methods support monophyly of enigmatic 'Mesozoa' (Dicyemida + Orthonectida, Lophotrochozoa)

Dicyemids and orthonectids were traditionally classified in a group called Mesozoa, but their placement in a single clade has been contested and their position(s) within Metazoa is uncertain. Here, we assembled a comprehensive matrix of Lophotrochozoa (Metazoa) and investigated the position of Dicyemida (= Rhombozoa) and Orthonectida, employing multiple phylogenomic approaches. We sequenced seven new transcriptomes and one draft genome from dicyemids (Dicyema, Dicyemennea) and two transcriptomes from orthonectids (Rhopalura). Using these and published data, we assembled and analysed contamination-filtered datasets with up to 987 genes. Our results recover Mesozoa monophyletic and as a close relative of Platyhelminthes or Gnathifera. Because of the tendency of the long-branch mesozoans to group with other long-branch taxa in our analyses, we explored the impact of approaches purported to help alleviate long-branch attraction (e.g. taxon removal, coalescent inference, gene targeting). None of these were able to break the association of Orthonectida with Dicyemida in the maximum-likelihood trees. Contrastingly, the Bayesian analysis and site-specific frequency model in maximum-likelihood did not recover a monophyletic Mesozoa (but only when using a specific 50 gene matrix). The classic hypothesis on monophyletic Mesozoa is possibly reborn and should be further tested.

A phylogenomic backbone for gastropod molluscs

Gastropods have survived several mass extinctions during their evolutionary history resulting in extraordinary diversity in morphology, ecology, and developmental modes, which complicate the reconstruction of a robust phylogeny. Currently, gastropods are divided into six subclasses: Caenogastropoda, Heterobranchia, Neomphaliones, Neritimorpha, Patellogastropoda, and Vetigastropoda. Phylogenetic relationships among these taxa historically lack consensus, despite numerous efforts using morphological and molecular information. We generated sequence data for transcriptomes derived from twelve taxa belonging to clades with little or no prior representation in previous studies in order to infer the deeper cladogenetic events within Gastropoda and, for the first time, infer the position of the deep-sea Neomphaliones using a phylogenomic approach. We explored the impact of missing data, homoplasy, and compositional heterogeneity on the inferred phylogenetic hypotheses. We recovered a highly supported backbone for gastropod relationships that is congruent with morphological and mitogenomic evidence, in which Patellogastropoda, true limpets, are the sister lineage to all other gastropods (Orthogastropoda) which are divided into two main clades (i) Vetigastropoda s.l. (including Pleurotomariida + Neomphaliones) and (ii) Neritimorpha + (Caenogastropoda + Heterobranchia). As such, our results support the recognition of five subclasses (or infraclasses) in Gastropoda: Patellogastropoda, Vetigastropoda, Neritimorpha, Caenogastropoda and Heterobranchia.

CView: A network based tool for enhanced alignment visualization

To date basic visualization of sequence alignments have largely focused on displaying per-site columns of nucleotide, or amino acid, residues along with associated frequency summarizations. The persistence of this tendency to the recent tools designed for viewing mapped read data indicates that such a perspective not only provides a reliable visualization of per-site alterations, but also offers implicit reassurance to the end-user in relation to data accessibility. However, the initial insight gained is limited, something that is especially true when viewing alignments consisting of many sequences representing differing factors such as location, date and subtype. A basic alignment viewer can have potential to increase initial insight through visual enhancement, whilst not delving into the realms of complex sequence analysis. We present CView, a visualizer that expands on the per-site representation of residues through the incorporation of a dynamic network that is based on the summarization of diversity present across different regions of the alignment. Within the network, nodes are based on the clustering of sequence fragments that span windows placed consecutively along the alignment. Edges are placed between nodes of neighbouring windows where they share sequence identification(s), i.e. different regions of the same sequence(s). Thus, if a node is selected on the network, then the relationship that sequences passing through that node have to other regions of diversity within the alignment can be observed through path tracing. In addition to augmenting visual insight, CView provides export features including variant summarization, per-site residue and kmer frequencies, consensus sequence, alignment dissection as well as clustering; each useful across a range of research areas. The software has been designed to be user friendly, intuitive and interactive. It is open source and an executable jar, source code, quick start, usage tutorial and test data are available (under the GNU General Public License) from https://sourceforge.net/projects/cview/.

Multigenerational laboratory culture of pelagic ctenophores and CRISPR-Cas9 genome editing in the lobate Mnemiopsis leidyi

Despite long-standing experimental interest in ctenophores due to their unique biology, ecological influence and evolutionary status, previous work has largely been constrained by the periodic seasonal availability of wild-caught animals and difficulty in reliably closing the life cycle. To address this problem, we have developed straightforward protocols that can be easily implemented to establish long-term multigenerational cultures for biological experimentation in the laboratory. In this protocol, we describe the continuous culture of the Atlantic lobate ctenophore Mnemiopsis leidyi. A rapid 3-week egg-to-egg generation time makes Mnemiopsis suitable for a wide range of experimental genetic, cellular, embryological, physiological, developmental, ecological and evolutionary studies. We provide recommendations for general husbandry to close the life cycle of Mnemiopsis in the laboratory, including feeding requirements, light-induced spawning, collection of embryos and rearing of juveniles to adults. These protocols have been successfully applied to maintain long-term multigenerational cultures of several species of pelagic ctenophores, and can be utilized by laboratories lacking easy access to the ocean. We also provide protocols for targeted genome editing via microinjection with CRISPR-Cas9 that can be completed within ~2 weeks, including single-guide RNA synthesis, early embryo microinjection, phenotype assessment and sequence validation of genome edits. These protocols provide a foundation for using Mnemiopsis as a model organism for functional genomic analyses in ctenophores.

Evolutionary plasticity in the requirement for force exerted by ligand endocytosis to activate C. elegans Notch proteins

The conserved transmembrane receptor Notch has diverse and profound roles in controlling cell fate during animal development. In the absence of ligand, a negative regulatory region (NRR) in the Notch ectodomain adopts an autoinhibited confirmation, masking an ADAM protease cleavage site;1,2 ligand binding induces cleavage of the NRR, leading to Notch ectodomain shedding as the first step of signal transduction.3,4 In Drosophila and vertebrates, recruitment of transmembrane Delta/Serrate/LAG-2 (DSL) ligands by the endocytic adaptor Epsin, and their subsequent internalization by Clathrin-mediated endocytosis, exerts a "pulling force" on Notch that is essential to expose the cleavage site in the NRR.4-6 Here, we show that Epsin-mediated endocytosis of transmembrane ligands is not essential to activate the two C. elegans Notch proteins, LIN-12 and GLP-1. Using an in vivo force sensing assay in Drosophila,6 we present evidence (1) that the LIN-12 and GLP-1 NRRs are tuned to lower force thresholds than the NRR of Drosophila Notch, and (2) that this difference depends on the absence of a "leucine plug" that occludes the cleavage site in the Drosophila and vertebrate Notch NRRs.1,2 Our results thus establish an unexpected evolutionary plasticity in the force-dependent mechanism of Notch activation and implicate a specific structural element, the leucine plug, as a determinant.

Insm1-expressing neurons and secretory cells develop from a common pool of progenitors in the sea anemone Nematostella vectensis

Neurons are highly specialized cells present in nearly all animals, but their evolutionary origin and relationship to other cell types are not well understood. We use here the sea anemone Nematostella vectensis as a model system for early-branching animals to gain fresh insights into the evolutionary history of neurons. We generated a transgenic reporter line to show that the transcription factor NvInsm1 is expressed in postmitotic cells that give rise to various types of neurons and secretory cells. Expression analyses, double transgenics, and gene knockdown experiments show that the NvInsm1-expressing neurons and secretory cells derive from a common pool of NvSoxB(2)-positive progenitor cells. These findings, together with the requirement for Insm1 for the development of neurons and endocrine cells in vertebrates, support a close evolutionary relationship of neurons and secretory cells.

Metazoans and Intrinsic Apoptosis: An Evolutionary Analysis of the Bcl-2 Family

The B-cell lymphoma-2 (Bcl-2) family is a group of genes regulating intrinsic apoptosis, a process controlling events such as development, homeostasis and the innate and adaptive immune responses in metazoans. In higher organisms, Bcl-2 proteins coordinate intrinsic apoptosis through their regulation of the integrity of the mitochondrial outer membrane; this function appears to have originated in the basal metazoans. Bcl-2 genes predate the cnidarian-bilaterian split and have been identified in porifera, placozoans and cnidarians but not ctenophores and some nematodes. The Bcl-2 family is composed of two groups of proteins, one with an α-helical Bcl-2 fold that has been identified in porifera, placozoans, cnidarians, and almost all higher bilaterians. The second group of proteins, the BH3-only group, has little sequence conservation and less well-defined structures and is found in cnidarians and most bilaterians, but not porifera or placozoans. Here we examine the evolutionary relationships between Bcl-2 proteins. We show that the structures of the Bcl-2-fold proteins are highly conserved over evolutionary time. Some metazoans such as the urochordate Oikopleura dioica have lost all Bcl-2 family members. This gene loss indicates that Bcl-2 regulated apoptosis is not an absolute requirement in metazoans, a finding mirrored in recent gene deletion studies in mice. Sequence analysis suggests that at least some Bcl-2 proteins lack the ability to bind BH3-only antagonists and therefore potentially have other non-apoptotic functions. By examining the foundations of the Bcl-2 regulated apoptosis, functional relationships may be clarified that allow us to understand the role of specific Bcl-2 proteins in evolution and disease.

The crustacean model Parhyale hawaiensis

Arthropods are the most abundant and diverse animals on earth. Among them, pancrustaceans are an ancient and morphologically diverse group, comprising a wide range of aquatic and semi-aquatic crustaceans as well as the insects, which emerged from crustacean ancestors to colonize most terrestrial habitats. Within insects, Drosophila stands out as one of the most powerful animal models, making major contributions to our understanding of development, physiology and behavior. Given these attributes, crustaceans provide a fertile ground for exploring biological diversity through comparative studies. However, beyond insects, few crustaceans are developed sufficiently as experimental models to enable such studies. The marine amphipod Parhyale hawaiensis is currently the best established crustacean system, offering year-round accessibility to developmental stages, transgenic tools, genomic resources, and established genetics and imaging approaches. The Parhyale research community is small but diverse, investigating the evolution of development, regeneration, aspects of sensory biology, chronobiology, bioprocessing and ecotoxicology.

Functional Conservation and Genetic Divergence of Chordate Glycinergic Neurotransmission: Insights from Amphioxus Glycine Transporters

Glycine is an important neurotransmitter in vertebrates, performing both excitatory and inhibitory actions. Synaptic levels of glycine are tightly controlled by the action of two glycine transporters, GlyT1 and GlyT2, located on the surface of glial cells and neurons, respectively. Only limited information is available on glycinergic neurotransmission in invertebrates, and the evolution of glycinergic neurotransmission is poorly understood. Here, by combining phylogenetic and gene expression analyses, we characterized the glycine transporter complement of amphioxus, an important invertebrate model for studying the evolution of chordates. We show that amphioxus possess three glycine transporter genes. Two of these (GlyT2.1 and GlyT2.2) are closely related to GlyT2 of vertebrates, whereas the third (GlyT) is a member of an ancestral clade of deuterostome glycine transporters. GlyT2.2 expression is predominantly non-neural, whereas GlyT and GlyT2.1 are widely expressed in the amphioxus nervous system and are differentially expressed, respectively, in neurons and glia. Vertebrate glycinergic neurons express GlyT2 and glia GlyT1, suggesting that the evolution of the chordate glycinergic system was accompanied by a paralog-specific inversion of gene expression. Despite this genetic divergence between amphioxus and vertebrates, we found strong evidence for conservation in the role glycinergic neurotransmission plays during larval swimming, the implication being that the neural networks controlling the rhythmic movement of chordate bodies may be homologous.

Whole-Genome Phylogenetic Reconstruction as a Powerful Tool to Reveal Homoplasy and Ancient Rapid Radiation in Waterflea Evolution

Although phylogeny estimation is notoriously difficult in radiations that occurred several hundred million years ago, phylogenomic approaches offer new ways to examine relationships among ancient lineages and evaluate hypotheses that are key to evolutionary biology. Here, we reconstruct the deep-rooted relationships of one of the oldest living arthropod clades, the branchiopod crustaceans, using a kaleidoscopic approach. We use concatenation and coalescent tree-building methods to analyze a large multigene data set at the nucleotide and amino acid level and examine gene tree versus species tree discordance. We unequivocally resolve long-debated relationships among extant orders of the Cladocera, the waterfleas, an ecologically relevant zooplankton group in global aquatic and marine ecosystems that is famous for its model systems in ecology and evolution. To build the data set, we assembled eight de novo genomes of key taxa including representatives of all extant cladoceran orders and suborders. Our phylogenetic analysis focused on a BUSCO-based set of 823 conserved single-copy orthologs shared among 23 representative taxa spanning all living branchiopod orders, including 11 cladoceran families. Our analysis supports the monophyly of the Cladocera and reveals remarkable homoplasy in their body plans. We found large phylogenetic distances between lineages with similar ecological specializations, indicating independent evolution in major body plans, such as in the pelagic predatory orders Haplopoda and Onychopoda (the “Gymnomera”). In addition, we assessed rapid cladogenesis by estimating relative timings of divergence in major lineages using reliable fossil-calibrated priors on eight nodes in the branchiopod tree, suggesting a Paleozoic origin around 325 Ma for the cladoceran ancestor and an ancient rapid radiation around 252 Ma at the Perm/Triassic boundary. These findings raise new questions about the roles of homoplasy and rapid radiation in the diversification of the cladocerans and help examine trait evolution from a genomic perspective in a functionally well understood, ancient arthropod group. [Cladocera; Daphnia; evolution; homoplasy; molecular clock; phylogenomics; systematics; waterfleas.]

Ancient role of sulfakinin/cholecystokinin-type signalling in inhibitory regulation of feeding processes revealed in an echinoderm

Sulfakinin (SK)/cholecystokinin (CCK)-type neuropeptides regulate feeding and digestion in protostomes (e.g. insects) and chordates. Here, we characterised SK/CCK-type signalling for the first time in a non-chordate deuterostome - the starfish Asterias rubens (phylum Echinodermata). In this species, two neuropeptides (ArSK/CCK1, ArSK/CCK2) derived from the precursor protein ArSK/CCKP act as ligands for an SK/CCK-type receptor (ArSK/CCKR) and these peptides/proteins are expressed in the nervous system, digestive system, tube feet, and body wall. Furthermore, ArSK/CCK1 and ArSK/CCK2 cause dose-dependent contraction of cardiac stomach, tube foot, and apical muscle preparations in vitro, and injection of these neuropeptides in vivo triggers cardiac stomach retraction and inhibition of the onset of feeding in A. rubens. Thus, an evolutionarily ancient role of SK/CCK-type neuropeptides as inhibitory regulators of feeding-related processes in the Bilateria has been conserved in the unusual and unique context of the extra-oral feeding behaviour and pentaradial body plan of an echinoderm.

DNA methylation atlas and machinery in the developing and regenerating annelid Platynereis dumerilii

BACKGROUND

Methylation of cytosines in DNA (5mC methylation) is a major epigenetic modification that modulates gene expression and constitutes the basis for mechanisms regulating multiple aspects of embryonic development and cell reprogramming in vertebrates. In mammals, 5mC methylation of promoter regions is linked to transcriptional repression. Transcription regulation by 5mC methylation notably involves the nucleosome remodeling and deacetylase complex (NuRD complex) which bridges DNA methylation and histone modifications. However, less is known about regulatory mechanisms involving 5mC methylation and their function in non-vertebrate animals. In this paper, we study 5mC methylation in the marine annelid worm Platynereis dumerilii, an emerging evolutionary and developmental biology model capable of regenerating the posterior part of its body post-amputation.

RESULTS

Using in silico and experimental approaches, we show that P. dumerilii displays a high level of DNA methylation comparable to that of mammalian somatic cells. 5mC methylation in P. dumerilii is dynamic along the life cycle of the animal and markedly decreases at the transition between larval to post-larval stages. We identify a full repertoire of mainly single-copy genes encoding the machinery associated with 5mC methylation or members of the NuRD complex in P. dumerilii and show that this repertoire is close to the one inferred for the last common ancestor of bilaterians. These genes are dynamically expressed during P. dumerilii development and regeneration. Treatment with the DNA hypomethylating agent Decitabine impairs P. dumerilii larval development and regeneration and has long-term effects on post-regenerative growth.

CONCLUSIONS

Our data reveal high levels of 5mC methylation in the annelid P. dumerilii, highlighting that this feature is not specific to vertebrates in the bilaterian clade. Analysis of DNA methylation levels and machinery gene expression during development and regeneration, as well as the use of a chemical inhibitor of DNA methylation, suggest an involvement of 5mC methylation in P. dumerilii development and regeneration. We also present data indicating that P. dumerilii constitutes a promising model to study biological roles and mechanisms of DNA methylation in non-vertebrate bilaterians and to provide new knowledge about evolution of the functions of this key epigenetic modification in bilaterian animals.

Nemertean, brachiopod and phoronid neuropeptidomics reveals ancestral spiralian signalling systems

Neuropeptides are diverse signaling molecules in animals commonly acting through G-protein coupled receptors (GPCRs). Neuropeptides and their receptors underwent extensive diversification in bilaterians and the relationships of many peptide–receptor systems have been clarified. However, we lack a detailed picture of neuropeptide evolution in lophotrochozoans as in-depth studies only exist for mollusks and annelids. Here, we analyze peptidergic systems in Nemertea, Brachiopoda, and Phoronida. We screened transcriptomes from 13 nemertean, 6 brachiopod, and 4 phoronid species for proneuropeptides and neuropeptide GPCRs. With mass spectrometry from the nemertean Lineus longissimus, we validated several predicted peptides and identified novel ones. Molecular phylogeny combined with peptide-sequence and gene-structure comparisons allowed us to comprehensively map spiralian neuropeptide evolution. We found most mollusk and annelid peptidergic systems also in nemerteans, brachiopods, and phoronids. We uncovered previously hidden relationships including the orthologies of spiralian CCWamides to arthropod agatoxin-like peptides and of mollusk APGWamides to RGWamides from annelids, with ortholog systems in nemerteans, brachiopods, and phoronids. We found that pleurin neuropeptides previously only found in mollusks are also present in nemerteans and brachiopods. We also identified cases of gene family duplications and losses. These include a protostome-specific expansion of RFamide/Wamide signaling, a spiralian expansion of GnRH-related peptides, and duplications of vasopressin/oxytocin before the divergence of brachiopods, phoronids, and nemerteans. This analysis expands our knowledge of peptidergic signaling in spiralians and other protostomes. Our annotated data set of nearly 1,300 proneuropeptide sequences and 600 GPCRs presents a useful resource for further studies of neuropeptide signaling.

18S rRNA variability maps reveal three highly divergent, conserved motifs within Rotifera

BACKGROUND

18S rRNA is a major component of the small subunit of the eukaryotic ribosome and an important phylogenetic marker for many groups, often to the point of being the only marker available for some. A core structure across eukaryotes exists for this molecule that can help to inform about its evolution in different groups. Using an alignment of 18S rDNA for Rotifera as traditionally recognized (=Bdelloidea, Monogononta, and Seisonacea, but not Acanthocephala), I fitted sequences for three exemplar species (Adineta vaga, Brachionus plicatilis, and Seison nebaliae, respectively) to the core structure and used these maps to reveal patterns of evolution for the remainder of this diverse group of microscopic animals.

RESULTS

The obtained variability maps of the 18S rRNA molecule revealed a pattern of high diversity among the three major rotifer clades coupled with strong conservation within each of bdelloids and monogononts. A majority of individual sites (ca. 60%) were constant even across rotifers as a whole with variable sites showing only intermediate rates of evolution. Although the three structural maps each showed good agreement with the inferred core structure for eukaryotic 18S rRNA and so were highly similar to one another at the secondary and tertiary levels, the overall pattern is of three highly distinct, but conserved motifs within the group at the primary sequence level. A novel finding was that of a variably expressed deletion at the 3' end of the V3 hypervariable region among some bdelloid species that occasionally extended into and included the pseudoknot structure following this region as well as the central "square" of the 18S rRNA molecule. Compared to other groups, levels of variation and rates of evolution for 18S rRNA in Rotifera roughly matched those for Gastropoda and Acanthocephala, despite increasing evidence for the latter being a clade within Rotifera.

CONCLUSIONS

The lack of comparative data for comparable groups makes interpretation of the results (i.e., very low variation within each of the three major rotifer clades, but high variation between them) and their potential novelty difficult. However, these findings in combination with the high morphological diversity within rotifers potentially help to explain why no clear consensus has been reached to date with regard to the phylogenetic relationships among the major groups.

Six-state Amino Acid Recoding is not an Effective Strategy to Offset Compositional Heterogeneity and Saturation in Phylogenetic Analyses

Six-state amino acid recoding strategies are commonly applied to combat the effects of compositional heterogeneity and substitution saturation in phylogenetic analyses. While these methods have been endorsed from a theoretical perspective, their performance has never been extensively tested. Here, we test the effectiveness of six-state recoding approaches by comparing the performance of analyses on recoded and non-recoded data sets that have been simulated under gradients of compositional heterogeneity or saturation. In our simulation analyses, non-recoding approaches consistently outperform six-state recoding approaches. Our results suggest that six-state recoding strategies are not effective in the face of high saturation. Furthermore, while recoding strategies do buffer the effects of compositional heterogeneity, the loss of information that accompanies six-state recoding outweighs its benefits. In addition, we evaluate recoding schemes with 9, 12, 15, and 18 states and show that these consistently outperform six-state recoding. Our analyses of other recoding schemes suggest that under conditions of very high compositional heterogeneity, it may be advantageous to apply recoding using more than six states, but we caution that applying any recoding should include sufficient justification. Our results have important implications for the more than 90 published papers that have incorporated six-state recoding, many of which have significant bearing on relationships across the tree of life. [Compositional heterogeneity; Dayhoff 6-state recoding; S&R 6-state recoding; six-state amino acid recoding; substitution saturation.]

Ancient DNA and paleogenetics: risks and potentiality

Paleopathology, the science that studies the diseases of the past, has always been addressed to the future in the use of new diagnostic methods. One of its relatively recent branches is paleogenetics, which is the study of genetic material from the past. Nuclear and mitochondrial DNA recovered from archaeological and paleontological specimens is called ancient DNA (aDNA), which can be extracted from a large variety of biological materials, of different origin, state of preservation and age, such as bones, teeth, coprolites, mummified tissues and hairs. There are many applications for ancient DNA research in the field of archaeology and paleopathology: population demography, genealogy, disease studies, archaeological reconstruction of plant vegetation, calibration of the molecular clock, phylogenetic relationship between different mammals and interpretation of the paleoclimate. However, the study of ancient genetic material is extremely difficult due to its poor quality and quantity, as well possible contamination with modern DNA. New advanced methods will allow extracting DNA from a greater variety of materials, and improvements in sequencing techniques will unveil data that are currently concealed.

The aim of this paper is to provide initial insights into paleogenetics and ancient DNA study and to illustrate the limits, risks and potentiality of the research on the genetic material of ancient specimens, whose results have a strong impact on the present and future medicine.

NvPOU4/Brain3 Functions as a Terminal Selector Gene in the Nervous System of the Cnidarian Nematostella vectensis

Terminal selectors are transcription factors that control the morphological, physiological, and molecular features that characterize distinct cell types. Here, we show that, in the sea anemone Nematostella vectensis, NvPOU4 is expressed in post-mitotic cells that give rise to a diverse set of neural cell types, including cnidocytes and NvElav1-expressing neurons. Morphological analyses of NvPOU4 mutants crossed to transgenic reporter lines show that the loss of NvPOU4 does not affect the initial specification of neural cells. Transcriptomes derived from the mutants and from different neural cell populations reveal that NvPOU4 is required for the execution of the terminal differentiation program of these neural cells. These findings suggest that POU4 genes have ancient functions as terminal selectors for morphologically and functionally disparate types of neurons and they provide experimental support for the relevance of terminal selectors for understanding the evolution of cell types.

Evidence for sponges as sister to all other animals from partitioned phylogenomics with mixture models and recoding

Resolving the relationships between the major lineages in the animal tree of life is necessary to understand the origin and evolution of key animal traits. Sponges, characterized by their simple body plan, were traditionally considered the sister group of all other animal lineages, implying a gradual increase in animal complexity from unicellularity to complex multicellularity. However, the availability of genomic data has sparked tremendous controversy as some phylogenomic studies support comb jellies taking this position, requiring secondary loss or independent origins of complex traits. Here we show that incorporating site-heterogeneous mixture models and recoding into partitioned phylogenomics alleviates systematic errors that hamper commonly-applied phylogenetic models. Testing on real datasets, we show a great improvement in model-fit that attenuates branching artefacts induced by systematic error. We reanalyse key datasets and show that partitioned phylogenomics does not support comb jellies as sister to other animals at either the supermatrix or partition-specific level.

Lack of support for Deuterostomia prompts reinterpretation of the first Bilateria

The bilaterally symmetric animals (Bilateria) are considered to comprise two monophyletic groups, Protostomia (Ecdysozoa and the Lophotrochozoa) and Deuterostomia (Chordata and the Xenambulacraria). Recent molecular phylogenetic studies have not consistently supported deuterostome monophyly. Here, we compare support for Protostomia and Deuterostomia using multiple, independent phylogenomic datasets. As expected, Protostomia is always strongly supported, especially by longer and higher-quality genes. Support for Deuterostomia, however, is always equivocal and barely higher than support for paraphyletic alternatives. Conditions that cause tree reconstruction errors-inadequate models, short internal branches, faster evolving genes, and unequal branch lengths-coincide with support for monophyletic deuterostomes. Simulation experiments show that support for Deuterostomia could be explained by systematic error. The branch between bilaterian and deuterostome common ancestors is, at best, very short, supporting the idea that the bilaterian ancestor may have been deuterostome-like. Our findings have important implications for the understanding of early animal evolution.

The origin of animals: an ancestral reconstruction of the unicellular-to-multicellular transition

How animals evolved from a single-celled ancestor, transitioning from a unicellular lifestyle to a coordinated multicellular entity, remains a fascinating question. Key events in this transition involved the emergence of processes related to cell adhesion, cell–cell communication and gene regulation. To understand how these capacities evolved, we need to reconstruct the features of both the last common multicellular ancestor of animals and the last unicellular ancestor of animals. In this review, we summarize recent advances in the characterization of these ancestors, inferred by comparative genomic analyses between the earliest branching animals and those radiating later, and between animals and their closest unicellular relatives. We also provide an updated hypothesis regarding the transition to animal multicellularity, which was likely gradual and involved the use of gene regulatory mechanisms in the emergence of early developmental and morphogenetic plans. Finally, we discuss some new avenues of research that will complement these studies in the coming years.

From Species to Regional and Local Specialization of Intestinal Macrophages

Initially intended for nutrient uptake, phagocytosis represents a central mechanism of debris removal and host defense against invading pathogens through the entire animal kingdom. In vertebrates and also many invertebrates, macrophages (MFs) and MF-like cells (e.g., coelomocytes and hemocytes) are professional phagocytic cells that seed tissues to maintain homeostasis through pathogen killing, efferocytosis and tissue shaping, repair, and remodeling. Some MF functions are common to all species and tissues, whereas others are specific to their homing tissue. Indeed, shaped by their microenvironment, MFs become adapted to perform particular functions, highlighting their great plasticity and giving rise to high population diversity. Interestingly, the gut displays several anatomic and functional compartments with large pools of strikingly diversified MF populations. This review focuses on recent advances on intestinal MFs in several species, which have allowed to infer their specificity and functions.

Evolution of new cell types at the lateral neural border

During the course of evolution, animals have become increasingly complex by the addition of novel cell types and regulatory mechanisms. A prime example is represented by the lateral neural border, known as the neural plate border in vertebrates, a region of the developing ectoderm where presumptive neural and non-neural tissue meet. This region has been intensively studied as the source of two important embryonic cell types unique to vertebrates-the neural crest and the ectodermal placodes-which contribute to diverse differentiated cell types including the peripheral nervous system, pigment cells, bone, and cartilage. How did these multipotent progenitors originate in animal evolution? What triggered the elaboration of the border during the course of chordate evolution? How is the lateral neural border patterned in various bilaterians and what is its fate? Here, we review and compare the development and fate of the lateral neural border in vertebrates and invertebrates and we speculate about its evolutionary origin. Taken together, the data suggest that the lateral neural border existed in bilaterian ancestors prior to the origin of vertebrates and became a developmental source of exquisite evolutionary change that frequently enabled the acquisition of new cell types.

Analysis of Fox genes in Schmidtea mediterranea reveals new families and a conserved role of Smed-foxO in controlling cell death

The forkhead box (Fox) genes encode transcription factors that control several key aspects of development. Present in the ancestor of all eukaryotes, Fox genes underwent several duplications followed by loss and diversification events that gave rise to the current 25 families. However, few Fox members have been identified from the Lophotrochozoa clade, and specifically from planarians, which are a unique model for understanding development, due to the striking plasticity of the adult. The aim of this study was to identify and perform evolutionary and functional studies of the Fox genes of lophotrochozoan species and, specifically, of the planarian Schmidtea mediterranea. Generating a pipeline for identifying Forkhead domains and using phylogenetics allowed us the phylogenetic reconstruction of Fox genes. We corrected the annotation for misannotated genes and uncovered a new family, the QD, present in all metazoans. According to the new phylogeny, the 27 Fox genes found in Schmidtea mediterranea were classified into 12 families. In Platyhelminthes, family losses were accompanied by extensive gene diversification and the appearance of specific families, the A(P) and N(P). Among the newly identified planarian Fox genes, we found a single copy of foxO, which shows an evolutionary conserved role in controlling cell death.

Sea Urchin as a Universal Model for Studies of Gene Networks

The purple sea urchin Strongylocentrotus purpuratus has been used for over 150 years as a model organism in developmental biology. Using this model species, scientists have been able to describe, in detail, the mechanisms of cell cycle control and cell adhesion, fertilization, calcium signaling, cell differentiation, and death. Massive parallel sequencing of the sea urchin genome enabled the deciphering of the main components of gene regulatory networks during the activation of embryonic signaling pathways. This knowledge helped to extrapolate aberrations in somatic cells that may lead to diseases, including cancer in humans. Furthermore, since many, if not all, developmental signaling pathways were shown to be controlled by non-coding RNAs (ncRNAs), the sea urchin organism represents an attractive experimental model. In this review, we discuss the main discoveries in the genetics, genomics, and transcriptomics of sea urchins during embryogenesis with the main focus on the role of ncRNAs. This information may be useful for comparative studies between different organisms, and may help identify new regulatory networks controlled by ncRNAs.

Uncovering DNA-PKcs ancient phylogeny, unique sequence motifs and insights for human disease

DNA-dependent protein kinase catalytic subunit (DNA-PKcs) is a key member of the phosphatidylinositol-3 kinase-like (PIKK) family of protein kinases with critical roles in DNA-double strand break repair, transcription, metastasis, mitosis, RNA processing, and innate and adaptive immunity. The absence of DNA-PKcs from many model organisms has led to the assumption that DNA-PKcs is a vertebrate-specific PIKK. Here, we find that DNA-PKcs is widely distributed in invertebrates, fungi, plants, and protists, and that threonines 2609, 2638, and 2647 of the ABCDE cluster of phosphorylation sites are highly conserved amongst most Eukaryotes. Furthermore, we identify highly conserved amino acid sequence motifs and domains that are characteristic of DNA-PKcs relative to other PIKKs. These include residues in the Forehead domain and a novel motif we have termed YRPD, located in an α helix C-terminal to the ABCDE phosphorylation site loop. Combining sequence with biochemistry plus structural data on human DNA-PKcs unveils conserved sequence and conformational features with functional insights and implications. The defined generally progressive DNA-PKcs sequence diversification uncovers conserved functionality supported by Evolutionary Trace analysis, suggesting that for many organisms both functional sites and evolutionary pressures remain identical due to fundamental cell biology. The mining of cancer genomic data and germline mutations causing human inherited disease reveal that robust DNA-PKcs activity in tumors is detrimental to patient survival, whereas germline mutations compromising function are linked to severe immunodeficiency and neuronal degeneration. We anticipate that these collective results will enable ongoing DNA-PKcs functional analyses with biological and medical implications.

Mouse WIF1 Is Only Modified with O-Fucose in Its EGF-like Domain III Despite Two Evolutionarily Conserved Consensus Sites

The Wnt Inhibitory Factor 1 (Wif1), known to inhibit Wnt signaling pathways, is composed of a WIF domain and five EGF-like domains (EGF-LDs) involved in protein interactions. Despite the presence of a potential O-fucosylation site in its EGF-LDs III and V, the O-fucose sites occupancy has never been demonstrated for WIF1. In this study, a phylogenetic analysis on the distribution, conservation and evolution of Wif1 proteins was performed, as well as biochemical approaches focusing on O-fucosylation sites occupancy of recombinant mouse WIF1. In the monophyletic group of gnathostomes, we showed that the consensus sequence for O-fucose modification by Pofut1 is highly conserved in Wif1 EGF-LD III while it was more divergent in EGF-LD V. Using click chemistry and mass spectrometry, we demonstrated that mouse WIF1 was only modified with a non-extended O-fucose on its EGF-LD III. In addition, a decreased amount of mouse WIF1 in the secretome of CHO cells was observed when the O-fucosylation site in EGF-LD III was mutated. Based on sequence comparison and automated protein modeling, we suggest that the absence of O-fucose on EGF-LD V of WIF1 in mouse and probably in most gnathostomes, could be related to EGF-LD V inability to interact with POFUT1.

Survival and selection biases in early animal evolution and a source of systematic overestimation in molecular clocks

Important evolutionary events such as the Cambrian Explosion have inspired many attempts at explanation: why do they happen when they do? What shapes them, and why do they eventually come to an end? However, much less attention has been paid to the idea of a 'null hypothesis'-that certain features of such diversifications arise simply through their statistical structure. Such statistical features also appear to influence our perception of the timing of these events. Here, we show in particular that study of unusually large clades leads to systematic overestimates of clade ages from some types of molecular clocks, and that the size of this effect may be enough to account for the puzzling mismatches seen between these molecular clocks and the fossil record. Our analysis of the fossil record of the late Ediacaran to Cambrian suggests that it is likely to be recording a true evolutionary radiation of the bilaterians at this time, and that explanations involving various sorts of cryptic origins for the bilaterians do not seem to be necessary.

The genetic factors of bilaterian evolution

The Cambrian explosion was a unique animal radiation ~540 million years ago that produced the full range of body plans across bilaterians. The genetic mechanisms underlying these events are unknown, leaving a fundamental question in evolutionary biology unanswered. Using large-scale comparative genomics and advanced orthology evaluation techniques, we identified 157 bilaterian-specific genes. They include the entire Nodal pathway, a key regulator of mesoderm development and left-right axis specification; components for nervous system development, including a suite of G-protein-coupled receptors that control physiology and behaviour, the Robo-Slit midline repulsion system, and the neurotrophin signalling system; a high number of zinc finger transcription factors; and novel factors that previously escaped attention. Contradicting the current view, our study reveals that genes with bilaterian origin are robustly associated with key features in extant bilaterians, suggesting a causal relationship.

The membrane receptors that appeared before their ligand: The different proposed scenarios

The interactions between membrane receptors and their endogenous ligands are key interactions in organisms. Recently, we have shown that a high number of genes encoding human receptors appeared at the same moment as their ligand in the animal tree of life. However, a set of receptors appeared before their present ligand. Different scenarios have been proposed to explain how a receptor can be conserved if its ligand is not yet appeared. However, these scenarios have been proposed individually and have never been studied in a global way. In this study, we investigated 30 mammalian pairs of ligand/receptor for which the first ligand appeared after its receptor in the tree of life, by using common indexes of selection, and proposed different scenarios explaining the earlier appearance of a receptor relative to its ligand. Based on 3D structural studies, our indexes allowed us to classify the evolution of these partners into different scenarios: 1) a scenario where the binding interface of the receptor is already present and under purifying selection before the appearance of the ligand; 2) a scenario where the binding interface seems to have appeared progressively, and 3) a scenario where the binding site seems to have been reshuffled since its appearance. As some scenarios were confirmed by the literature, we concluded that simple indexes can give a good highlight of the evolutive history of two partners that did not appear at the same time. Based on these scenarios, we also hypothesize that the replacement of a ligand by another is a frequent phenomenon during evolution.

Phylogenetic tree building in the genomic age

Knowing phylogenetic relationships among species is fundamental for many studies in biology. An accurate phylogenetic tree underpins our understanding of the major transitions in evolution, such as the emergence of new body plans or metabolism, and is key to inferring the origin of new genes, detecting molecular adaptation, understanding morphological character evolution and reconstructing demographic changes in recently diverged species. Although data are ever more plentiful and powerful analysis methods are available, there remain many challenges to reliable tree building. Here, we discuss the major steps of phylogenetic analysis, including identification of orthologous genes or proteins, multiple sequence alignment, and choice of substitution models and inference methodologies. Understanding the different sources of errors and the strategies to mitigate them is essential for assembling an accurate tree of life.

A gonadotropin-releasing hormone type neuropeptide with a high affinity binding site for copper(ii) and nickel(ii)

In vertebrates gonadotropin-releasing hormone I (GnRH-I) is a key regulator of reproductive development and function. The receptor-binding activity of human GnRH-I can be modified by the presence of divalent copper. Thus, copper binding to N-terminal amino acids in GnRH-I induces structural changes that influence receptor interactions and downstream intracellular signalling cascades. It is not known if copper-binding is restricted to human GnRH-I or if it is also a feature of GnRH-type peptides that have been identified in other taxa. To investigate this, we have characterised copper binding to a recently discovered GnRH-type peptide from the starfish Asterias rubens (ArGnRH). Using a range of spectroscopic and biophysical techniques we show that this peptide can bind copper(ii) and nickel(ii). Copper(ii) is bound in a square-planar, high-affinity (Kd ∼ 10-12 M) site incorporating four nitrogen donor atoms from a histidine imidazole group, two amides and the N-terminal amine group. The ArGnRH copper affinity and geometry are quite different to GnRH-I suggesting the copper sites have evolved to suit the environment the peptides are exposed to. By comparing the copper binding sites in ArGnRH and human GnRH-I and conducting a phylogenetic analysis of GnRH-type peptide sequences from a range of species, we predict that copper-binding is an evolutionarily ancient feature of GnRH-type peptides that has been retained, modified or lost in different lineages.

Cell-type diversity and regionalized gene expression in the planarian intestine

Proper function and repair of the digestive system are vital to most animals. Deciphering the mechanisms involved in these processes requires an atlas of gene expression and cell types. Here, we applied laser-capture microdissection (LCM) and RNA-seq to characterize the intestinal transcriptome of Schmidtea mediterranea, a planarian flatworm that can regenerate all organs, including the gut. We identified hundreds of genes with intestinal expression undetected by previous approaches. Systematic analyses revealed extensive conservation of digestive physiology and cell types with other animals, including humans. Furthermore, spatial LCM enabled us to uncover previously unappreciated regionalization of gene expression in the planarian intestine along the medio-lateral axis, especially among intestinal goblet cells. Finally, we identified two intestine-enriched transcription factors that specifically regulate regeneration (hedgehog signaling effector gli-1) or maintenance (RREB2) of goblet cells. Altogether, this work provides resources for further investigation of mechanisms involved in gastrointestinal function, repair and regeneration.

Phylogenetic Analyses of Sites in Different Protein Structural Environments Result in Distinct Placements of the Metazoan Root

Phylogenomics, the use of large datasets to examine phylogeny, has revolutionized the study of evolutionary relationships. However, genome-scale data have not been able to resolve all relationships in the tree of life; this could reflect, at least in part, the poor-fit of the models used to analyze heterogeneous datasets. Some of the heterogeneity may reflect the different patterns of selection on proteins based on their structures. To test that hypothesis, we developed a pipeline to divide phylogenomic protein datasets into subsets based on secondary structure and relative solvent accessibility. We then tested whether amino acids in different structural environments had distinct signals for the topology of the deepest branches in the metazoan tree. We focused on a dataset that appeared to have a mixture of signals and we found that the most striking difference in phylogenetic signal reflected relative solvent accessibility. Analyses of exposed sites (residues located on the surface of proteins) yielded a tree that placed ctenophores sister to all other animals whereas sites buried inside proteins yielded a tree with a sponge+ctenophore clade. These differences in phylogenetic signal were not ameliorated when we conducted analyses using a set of maximum-likelihood profile mixture models. These models are very similar to the Bayesian CAT model, which has been used in many analyses of deep metazoan phylogeny. In contrast, analyses conducted after recoding amino acids to limit the impact of deviations from compositional stationarity increased the congruence in the estimates of phylogeny for exposed and buried sites; after recoding amino acid trees estimated using the exposed and buried site both supported placement of ctenophores sister to all other animals. Although the central conclusion of our analyses is that sites in different structural environments yield distinct trees when analyzed using models of protein evolution, our amino acid recoding analyses also have implications for metazoan evolution. Specifically, our results add to the evidence that ctenophores are the sister group of all other animals and they further suggest that the placozoa+cnidaria clade found in some other studies deserves more attention. Taken as a whole, these results provide striking evidence that it is necessary to achieve a better understanding of the constraints due to protein structure to improve phylogenetic estimation.

More Arrows in the Ancient DNA Quiver: Use of Paleoepigenomes and Paleomicrobiomes to Investigate Animal Adaptation to Environment

Whether and how epigenetic mechanisms and the microbiome play a role in mammalian adaptation raised considerable attention and controversy, mainly because they have the potential to add new insights into the Modern Synthesis. Recent attempts to reconcile neo-Darwinism and neo-Lamarckism in a unified theory of molecular evolution give epigenetic mechanisms and microbiome a prominent role. However, supporting empirical data are still largely missing. Because experimental studies using extant animals can hardly be done over evolutionary timescales, we propose that advances in ancient DNA techniques provide a valid alternative. In this piece, we evaluate 1) the possible roles of epigenomes and microbiomes in animal adaptation, 2) advances in the retrieval of paleoepigenome and paleomicrobiome data using ancient DNA techniques, and 3) the plasticity of either and interactions between the epigenome and the microbiome, while emphasizing that it is essential to take both into account, as well as the underlying genetic factors that may confound the findings. We propose that advanced ancient DNA techniques should be applied to a wide range of past animals, so novel dynamics in animal evolution and adaption can be revealed.

Model Choice, Missing Data, and Taxon Sampling Impact Phylogenomic Inference of Deep Basidiomycota Relationships

Resolving deep divergences in the tree of life is challenging even for analyses of genome-scale phylogenetic data sets. Relationships between Basidiomycota subphyla, the rusts and allies (Pucciniomycotina), smuts and allies (Ustilaginomycotina), and mushroom-forming fungi and allies (Agaricomycotina) were found particularly recalcitrant both to traditional multigene and genome-scale phylogenetics. Here, we address basal Basidiomycota relationships using concatenated and gene tree-based analyses of various phylogenomic data sets to examine the contribution of several potential sources of bias. We evaluate the contribution of biological causes (hard polytomy, incomplete lineage sorting) versus unmodeled evolutionary processes and factors that exacerbate their effects (e.g., fast-evolving sites and long-branch taxa) to inferences of basal Basidiomycota relationships. Bayesian Markov Chain Monte Carlo and likelihood mapping analyses reject the hard polytomy with confidence. In concatenated analyses, fast-evolving sites and oversimplified models of amino acid substitution favored the grouping of smuts with mushroom-forming fungi, often leading to maximal bootstrap support in both concatenation and coalescent analyses. On the contrary, the most conserved data subsets grouped rusts and allies with mushroom-forming fungi, although this relationship proved labile, sensitive to model choice, to different data subsets and to missing data. Excluding putative long-branch taxa, genes with high proportions of missing data and/or with strong signal failed to reveal a consistent trend toward one or the other topology, suggesting that additional sources of conflict are at play. While concatenated analyses yielded strong but conflicting support, individual gene trees mostly provided poor support for any resolution of rusts, smuts, and mushroom-forming fungi, suggesting that the true Basidiomycota tree might be in a part of tree space that is difficult to access using both concatenation and gene tree-based approaches. Inference-based assessments of absolute model fit strongly reject best-fit models for the vast majority of genes, indicating a poor fit of even the most commonly used models. While this is consistent with previous assessments of site-homogenous models of amino acid evolution, this does not appear to be the sole source of confounding signal. Our analyses suggest that topologies uniting smuts with mushroom-forming fungi can arise as a result of inappropriate modeling of amino acid sites that might be prone to systematic bias. We speculate that improved models of sequence evolution could shed more light on basal splits in the Basidiomycota, which, for now, remain unresolved despite the use of whole genome data.

Fast speciations and slow genes: uncovering the root of living canids

Despite ongoing efforts relying on computationally intensive tree-building methods and large datasets, the deeper phylogenetic relationships between living canid genera remain controversial. We demonstrate that this issue arises fundamentally from the uncertainty of root placement as a consequence of the short length of the branch connecting the major canid clades, which probably resulted from a fast radiation during the early diversification of extant Canidae. Using both nuclear and mitochondrial genes, we investigate the position of the canid root and its consistency by using three rooting methods. We find that mitochondrial genomes consistently retrieve a root node separating the tribe Canini from the remaining canids, whereas nuclear data mostly recover a root that places the Urocyon foxes as the sister lineage of living canids. We demonstrate that, to resolve the canid root, the nuclear segments sequenced so far are significantly less informative than mitochondrial genomes. We also propose that short intervals between speciations obscure the place of the true root, because methods are susceptible to stochastic error in the presence of short internal branches near the root.