Increased taxon sampling reveals thousands of hidden orthologs in flatworms

Allatotropin (AT) related peptides L-ATRP and D2-ATRP diastereomers activate an endogenous receptor and suppress heart rate in the Pacific abalone Haliotis discus hannai

Allatotropin (AT) has been identified in many insects and plays important roles in the regulation of their intestinal contraction, heart rate, ion transport, and digestive enzyme secretion. However, information on AT-related bioinformatics in other animal phyla is scarce. In this study, we cloned a full-length cDNA encoding the AT-related peptide receptor (ATRPR) of the abalone Haliotis discus hannai (Hdh) and further characterized Hdh-ATRPR with its potential ligands, Hdh-ATRPs. In luciferase reporter and Ca2+ mobilization assays, Hdh-ATRPs, including a D-type Phe at the second amino acid position, Hdh-D2-ATRP, activated Hdh-ATRPR in a dose-dependent manner, whereas all-L-type Hdh-ATRP was a more potent ligand than Hdh-D2-ATRP. Furthermore, Hdh-ATRPs induced ERK1/2 phosphorylation in Hdh-ATRPR-expressing HEK293 cells, which was dose-dependently abolished by the PKC inhibitor Gö6983. The heart rate decreased significantly within 10 min when Hdh-D2-ATRP was injected into the adduct muscle sinus of abalone (0.2 or 1.0 µg/g body weight), while the abalone injected with a high concentration of Hdh-D2-ATRP (1.5 μg/g body weight) were sublethal within 5 h. Thus, Hdh-ATRP signaling is primarily linked to the Gαq/PKC and is possibly associated with heart rate regulation in abalone.

A comparative analysis of planarian genomes reveals regulatory conservation in the face of rapid structural divergence

The planarian Schmidtea mediterranea is being studied as a model species for regeneration, but the assembly of planarian genomes remains challenging. Here, we report a high-quality haplotype-phased, chromosome-scale genome assembly of the sexual S2 strain of S. mediterranea and high-quality chromosome-scale assemblies of its three close relatives, S. polychroa, S. nova, and S. lugubris. Using hybrid gene annotations and optimized ATAC-seq and ChIP-seq protocols for regulatory element annotation, we provide valuable genome resources for the planarian research community and a first comparative perspective on planarian genome evolution. Our analyses reveal substantial divergence in protein-coding sequences and regulatory regions but considerable conservation within promoter and enhancer annotations. We also find frequent retrotransposon-associated chromosomal inversions and interchromosomal translocations within the genus Schmidtea and, remarkably, independent and nearly complete losses of ancestral metazoan synteny in Schmidtea and two other flatworm groups. Overall, our results suggest that platyhelminth genomes can evolve without syntenic constraints.

Hookworm genomics: dusk or dawn?

Hookworms are parasites, closely related to the model nematode Caenorhabditis elegans, that are a major economic and health burden worldwide. Primarily three hookworm species (Necator americanus, Ancylostoma duodenale, and Ancylostoma ceylanicum) infect humans. Another 100 hookworm species from 19 genera infect primates, ruminants, and carnivores. Genetic data exist for only seven of these species. Genome sequences are available from only four of these species in two genera, leaving 96 others (particularly those parasitizing wildlife) without any genomic data. The most recent hookworm genomes were published 5 years ago, leaving the field in a dusk. However, assembling genomes from single hookworms may bring a new dawn. Here we summarize advances, challenges, and opportunities for studying these neglected but important parasitic nematodes.

The promise and pitfalls of synteny in phylogenomics

Reconstructing the tree of life remains a central goal in biology. Early methods, which relied on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the tree of life, 2 recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in 2 or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. This Essay identifies and discusses specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the tree of life.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Spiralian genomics and the evolution of animal genome architecture

Recent developments in sequencing technologies have greatly improved our knowledge of phylogenetic relationships and genomic architectures throughout the tree of life. Spiralia, a diverse clade within Protostomia, is essential for understanding the evolutionary history of parasitism, gene conversion, nervous systems and animal body plans. In this review, we focus on the current hypotheses of spiralian phylogeny and investigate the impact of long-read sequencing on the quality of genome assemblies. We examine chromosome-level assemblies to highlight key genomic features that have driven spiralian evolution, including karyotype, synteny and the Hox gene organization. In addition, we show how chromosome rearrangement has influenced spiralian genomic structures. Although spiralian genomes have undergone substantial changes, they exhibit both conserved and lineage-specific features. We recommend increasing sequencing efforts and expanding functional genomics research to deepen insights into spiralian biology.

The NIH Comparative Genomics Resource: addressing the promises and challenges of comparative genomics on human health

Comparative genomics is the comparison of genetic information within and across organisms to understand the evolution, structure, and function of genes, proteins, and non-coding regions (Sivashankari and Shanmughavel, Bioinformation 1:376-8, 2007). Advances in sequencing technology and assembly algorithms have resulted in the ability to sequence large genomes and provided a wealth of data that are being used in comparative genomic analyses. Comparative analysis can be leveraged to systematically explore and evaluate the biological relationships and evolution between species, aid in understanding the structure and function of genes, and gain a better understanding of disease and potential drug targets. As our knowledge of genetics expands, comparative genomics can help identify emerging model organisms among a broader span of the tree of life, positively impacting human health. This impact includes, but is not limited to, zoonotic disease research, therapeutics development, microbiome research, xenotransplantation, oncology, and toxicology. Despite advancements in comparative genomics, new challenges have arisen around the quantity, quality assurance, annotation, and interoperability of genomic data and metadata. New tools and approaches are required to meet these challenges and fulfill the needs of researchers. This paper focuses on how the National Institutes of Health (NIH) Comparative Genomics Resource (CGR) can address both the opportunities for comparative genomics to further impact human health and confront an increasingly complex set of challenges facing researchers.

Distinct genomic routes underlie transitions to specialised symbiotic lifestyles in deep-sea annelid worms

Bacterial symbioses allow annelids to colonise extreme ecological niches, such as hydrothermal vents and whale falls. Yet, the genetic principles sustaining these symbioses remain unclear. Here, we show that different genomic adaptations underpin the symbioses of phylogenetically related annelids with distinct nutritional strategies. Genome compaction and extensive gene losses distinguish the heterotrophic symbiosis of the bone-eating worm Osedax frankpressi from the chemoautotrophic symbiosis of deep-sea Vestimentifera. Osedax's endosymbionts complement many of the host's metabolic deficiencies, including the loss of pathways to recycle nitrogen and synthesise some amino acids. Osedax's endosymbionts possess the glyoxylate cycle, which could allow more efficient catabolism of bone-derived nutrients and the production of carbohydrates from fatty acids. Unlike in most Vestimentifera, innate immunity genes are reduced in O. frankpressi, which, however, has an expansion of matrix metalloproteases to digest collagen. Our study supports that distinct nutritional interactions influence host genome evolution differently in highly specialised symbioses.

Genome-Scale Analysis Reveals Extensive Diversification of Voltage-Gated K+ Channels in Stem Cnidarians

Ion channels are highly diverse in the cnidarian model organism Nematostella vectensis (Anthozoa), but little is known about the evolutionary origins of this channel diversity and its conservation across Cnidaria. Here, we examined the evolution of voltage-gated K+ channels in Cnidaria by comparing genomes and transcriptomes of diverse cnidarian species from Anthozoa and Medusozoa. We found an average of over 40 voltage-gated K+ channel genes per species, and a phylogenetic reconstruction of the Kv, KCNQ, and Ether-a-go-go (EAG) gene families identified 28 voltage-gated K+ channels present in the last common ancestor of Anthozoa and Medusozoa (23 Kv, 1 KCNQ, and 4 EAG). Thus, much of the diversification of these channels took place in the stem cnidarian lineage prior to the emergence of modern cnidarian classes. In contrast, the stem bilaterian lineage, from which humans evolved, contained no more than nine voltage-gated K+ channels. These results hint at a complexity to electrical signaling in all cnidarians that contrasts with the perceived anatomical simplicity of their neuromuscular systems. These data provide a foundation from which the function of these cnidarian channels can be investigated, which will undoubtedly provide important insights into cnidarian physiology.

Marine animal evolutionary developmental biology-Advances through technology development

Evolutionary developmental biology, the interdisciplinary effort of illuminating the conserved similarities and differences during animal development across all phylogenetic clades, has gained renewed interest in the past decades. As technology (immunohistochemistry, next-generation sequencing, advanced imaging, and computational resources) has advanced, so has our ability of resolving fundamental hypotheses and overcoming the genotype-phenotype gap. This rapid progress, however, has also exposed gaps in the collective knowledge around the choice and representation of model organisms. It has become clear that evo-devo requires a comparative, large-scale approach including marine invertebrates to resolve some of the most urgent questions about the phylogenetic positioning and character traits of the last common ancestors. Many invertebrates at the base of the tree of life inhabit marine environments and have been used for some years due to their accessibility, husbandry, and morphology. Here, we briefly review the major concepts of evolutionary developmental biology and discuss the suitability of established model organisms to address current research questions, before focussing on the importance, application, and state-of-the-art of marine evo-devo. We highlight novel technical advances that progress evo-devo as a whole.

Building Phylogenies from Transcriptomic Data

Transcriptomic data (obtained from RNA sequencing) has become a very powerful source of information to reconstruct the evolutionary relationships among organisms. Although phylogenetic inference using transcriptomes retains the same core steps as when working with few molecular markers (viz., nucleic acid extraction and sequencing, sequence treatment, and tree inference), all of them show significant differences. First, the needed quantity and quality of the extracted RNA has to be very high. Although this may not represent a challenge when working with certain organisms, it may well be a headache with others, especially for those with small body sizes. Second, the tremendous increase in the quantity of sequences obtained requires a high computational power for both treating the sequences and inferring the subsequent phylogenies. This means that transcriptomic data can no longer be analyzed using personal computers nor local programs with a graphical interface. This, in turn, implies the requirement of an increased set of bioinformatic skills from the researchers. Finally, the genomic peculiarities of each group of organisms, such as the level of heterozygosity or the percentage of base composition, also need to be considered when inferring phylogenies using transcriptomic data.

The state of Medusozoa genomics: current evidence and future challenges

Medusozoa is a widely distributed ancient lineage that harbors one-third of Cnidaria diversity divided into 4 classes. This clade is characterized by the succession of stages and modes of reproduction during metagenic lifecycles, and includes some of the most plastic body plans and life cycles among animals. The characterization of traditional genomic features, such as chromosome numbers and genome sizes, was rather overlooked in Medusozoa and many evolutionary questions still remain unanswered. Modern genomic DNA sequencing in this group started in 2010 with the publication of the Hydra vulgaris genome and has experienced an exponential increase in the past 3 years. Therefore, an update of the state of Medusozoa genomics is warranted. We reviewed different sources of evidence, including cytogenetic records and high-throughput sequencing projects. We focused on 4 main topics that would be relevant for the broad Cnidaria research community: (i) taxonomic coverage of genomic information; (ii) continuity, quality, and completeness of high-throughput sequencing datasets; (iii) overview of the Medusozoa specific research questions approached with genomics; and (iv) the accessibility of data and metadata. We highlight a lack of standardization in genomic projects and their reports, and reinforce a series of recommendations to enhance future collaborative research.

Harmonizing model organism data in the Alliance of Genome Resources

The Alliance of Genome Resources (the Alliance) is a combined effort of 7 knowledgebase projects: Saccharomyces Genome Database, WormBase, FlyBase, Mouse Genome Database, the Zebrafish Information Network, Rat Genome Database, and the Gene Ontology Resource. The Alliance seeks to provide several benefits: better service to the various communities served by these projects; a harmonized view of data for all biomedical researchers, bioinformaticians, clinicians, and students; and a more sustainable infrastructure. The Alliance has harmonized cross-organism data to provide useful comparative views of gene function, gene expression, and human disease relevance. The basis of the comparative views is shared calls of orthology relationships and the use of common ontologies. The key types of data are alleles and variants, gene function based on gene ontology annotations, phenotypes, association to human disease, gene expression, protein-protein and genetic interactions, and participation in pathways. The information is presented on uniform gene pages that allow facile summarization of information about each gene in each of the 7 organisms covered (budding yeast, roundworm Caenorhabditis elegans, fruit fly, house mouse, zebrafish, brown rat, and human). The harmonized knowledge is freely available on the alliancegenome.org portal, as downloadable files, and by APIs. We expect other existing and emerging knowledge bases to join in the effort to provide the union of useful data and features that each knowledge base currently provides.

Regeneration of Planarian Auricles and Reestablishment of Chemotactic Ability

Detection of chemical stimuli is crucial for living systems and also contributes to quality of life in humans. Since loss of olfaction becomes more prevalent with aging, longer life expectancies have fueled interest in understanding the molecular mechanisms behind the development and maintenance of chemical sensing. Planarian flatworms possess an unsurpassed ability for stem cell-driven regeneration that allows them to restore any damaged or removed part of their bodies. This includes anteriorly-positioned lateral flaps known as auricles, which have long been thought to play a central role in chemotaxis. The contribution of auricles to the detection of positive chemical stimuli was tested in this study using Girardia dorotocephala, a North American planarian species known for its morphologically prominent auricles. Behavioral experiments staged under laboratory conditions revealed that removal of auricles by amputation leads to a significant decrease in the ability of planarians to find food. However, full chemotactic capacity is observed as early as 2 days post-amputation, which is days prior from restoration of auricle morphology, but correlative with accumulation of ciliated cells in the position of auricle regeneration. Planarians subjected to x-ray irradiation prior to auricle amputation were unable to restore auricle morphology, but were still able to restore chemotactic capacity. These results indicate that although regeneration of auricle morphology requires stem cells, some restoration of chemotactic ability can still be achieved in the absence of normal auricle morphology, corroborating with the initial observation that chemotactic success is reestablished 2-days post-amputation in our assays. Transcriptome profiles of excised auricles were obtained to facilitate molecular characterization of these structures, as well as the identification of genes that contribute to chemotaxis and auricle development. A significant overlap was found between genes with preferential expression in auricles of G. dorotocephala and genes with reduced expression upon SoxB1 knockdown in Schmidtea mediterranea, suggesting that SoxB1 has a conserved role in regulating auricle development and function. Models that distinguish between possible contributions to chemotactic behavior obtained from cellular composition, as compared to anatomical morphology of the auricles, are discussed.

Tracing Eukaryotic Ribosome Biogenesis Factors Into the Archaeal Domain Sheds Light on the Evolution of Functional Complexity

Ribosome assembly is an essential and carefully choreographed cellular process. In eukaryotes, several 100 proteins, distributed across the nucleolus, nucleus, and cytoplasm, co-ordinate the step-wise assembly of four ribosomal RNAs (rRNAs) and approximately 80 ribosomal proteins (RPs) into the mature ribosomal subunits. Due to the inherent complexity of the assembly process, functional studies identifying ribosome biogenesis factors and, more importantly, their precise functions and interplay are confined to a few and very well-established model organisms. Although best characterized in yeast (Saccharomyces cerevisiae), emerging links to disease and the discovery of additional layers of regulation have recently encouraged deeper analysis of the pathway in human cells. In archaea, ribosome biogenesis is less well-understood. However, their simpler sub-cellular structure should allow a less elaborated assembly procedure, potentially providing insights into the functional essentials of ribosome biogenesis that evolved long before the diversification of archaea and eukaryotes. Here, we use a comprehensive phylogenetic profiling setup, integrating targeted ortholog searches with automated scoring of protein domain architecture similarities and an assessment of when search sensitivity becomes limiting, to trace 301 curated eukaryotic ribosome biogenesis factors across 982 taxa spanning the tree of life and including 727 archaea. We show that both factor loss and lineage-specific modifications of factor function modulate ribosome biogenesis, and we highlight that limited sensitivity of the ortholog search can confound evolutionary conclusions. Projecting into the archaeal domain, we find that only few factors are consistently present across the analyzed taxa, and lineage-specific loss is common. While members of the Asgard group are not special with respect to their inventory of ribosome biogenesis factors (RBFs), they unite the highest number of orthologs to eukaryotic RBFs in one taxon. Using large ribosomal subunit maturation as an example, we demonstrate that archaea pursue a simplified version of the corresponding steps in eukaryotes. Much of the complexity of this process evolved on the eukaryotic lineage by the duplication of ribosomal proteins and their subsequent functional diversification into ribosome biogenesis factors. This highlights that studying ribosome biogenesis in archaea provides fundamental information also for understanding the process in eukaryotes.

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.

Systematic errors in orthology inference and their effects on evolutionary analyses

The availability of complete sets of genes from many organisms makes it possible to identify genes unique to (or lost from) certain clades. This information is used to reconstruct phylogenetic trees; identify genes involved in the evolution of clade specific novelties; and for phylostratigraphy—identifying ages of genes in a given species. These investigations rely on accurately predicted orthologs. Here we use simulation to produce sets of orthologs that experience no gains or losses. We show that errors in identifying orthologs increase with higher rates of evolution. We use the predicted sets of orthologs, with errors, to reconstruct phylogenetic trees; to count gains and losses; and for phylostratigraphy. Our simulated data, containing information only from errors in orthology prediction, closely recapitulate findings from empirical data. We suggest published downstream analyses must be informed to a large extent by errors in orthology prediction that mimic expected patterns of gene evolution.

•

Presence of shared orthologs across species is used for evolutionary analyses

•

We simulated realistic sets of orthologs with no gains or losses

•

Errors predicting shared orthologs correlate with phylogenetic relationships

•

Presence/absence datasets based on errors recapitulate findings from empirical data

Weighted gene co-expression network analysis identifies modules and functionally enriched pathways in the lactation process

The exponential growth in knowledge has resulted in a better understanding of the lactation process in a wide variety of animals. However, the underlying genetic mechanisms are not yet clearly known. In order to identify the mechanisms involved in the lactation process, various mehods, including meta-analysis, weighted gene co-express network analysis (WGCNA), hub genes identification, gene ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment at before peak (BP), peak (P), and after peak (AP) stages of the lactation processes have been employed. A total of 104, 85, and 26 differentially expressed genes were identified based on PB vs. P, BP vs. AP, and P vs. AP comparisons, respectively. GO and KEGG pathway enrichment analysis revealed that DEGs were significantly enriched in the "ubiquitin-dependent ERAD" and the "chaperone cofactor-dependent protein refolding" in BP vs. P and P vs. P, respectively. WGCNA identified five significant functional modules related to the lactation process. Moreover, GJA1, AP2A2, and NPAS3 were defined as hub genes in the identified modules, highlighting the importance of their regulatory impacts on the lactation process. The findings of this study provide new insights into the complex regulatory networks of the lactation process at three distinct stages, while suggesting several candidate genes that may be useful for future animal breeding programs. Furthermore, this study supports the notion that in combination with a meta-analysis, the WGCNA represents an opportunity to achieve a higher resolution analysis that can better predict the most important functional genes that might provide a more robust bio-signature for phenotypic traits, thus providing more suitable biomarker candidates for future studies.

Prospective enzymes for omega-3 PUFA biosynthesis found in endoparasitic classes within the phylum Platyhelminthes

The free-living infectious stages of macroparasites, specifically, the cercariae of trematodes (flatworms), are likely to be significant (albeit underappreciated) vectors of nutritionally important polyunsaturated fatty acids (PUFA) to consumers within aquatic food webs, and other macroparasites could serve similar roles. In the context of de novo omega-3 (n-3) PUFA biosynthesis, it was thought that most animals lack the fatty acid (FA) desaturase enzymes that convert stearic acid (18:0) into ɑ-linolenic acid (ALA; 18:3n-3), the main FA precursor for n-3 long-chain PUFA. Recently, novel sequences of these enzymes were recovered from 80 species from six invertebrate phyla, with experimental confirmation of gene function in five phyla. Given this wide distribution, and the unusual attributes of flatworm genomes, we conducted an additional search for genes for de novo n-3 PUFA in the phylum Platyhelminthes. Searches with experimentally confirmed sequences from Rotifera recovered nine relevant FA desaturase sequences from eight species in four genera in the two exclusively endoparasite classes (Trematoda and Cestoda). These results could indicate adaptations of these particular parasite species, or may reflect the uneven taxonomic coverage of sequence databases. Although additional genomic data and, particularly, experimental study of gene functionality are important future validation steps, our results indicate endoparasitic platyhelminths may have enzymes for de novo n-3 PUFA biosynthesis, thereby contributing to global PUFA production, but also representing a potential target for clinical antihelmintic applications.

PlanMine 3.0-improvements to a mineable resource of flatworm biology and biodiversity

Flatworms (Platyhelminthes) are a basally branching phylum that harbours a wealth of fascinating biology, including planarians with their astonishing regenerative abilities and the parasitic tape worms and blood flukes that exert a massive impact on human health. PlanMine (http://planmine.mpi-cbg.de/) has the mission objective of providing both a mineable sequence repository for planarians and also a resource for the comparative analysis of flatworm biology. While the original PlanMine release was entirely based on transcriptomes, the current release transitions to a more genomic perspective. Building on the recent availability of a high quality genome assembly of the planarian model species Schmidtea mediterranea, we provide a gene prediction set that now assign existing transcripts to defined genomic coordinates. The addition of recent single cell and bulk RNA-seq datasets greatly expands the available gene expression information. Further, we add transcriptomes from a broad range of other flatworms and provide a phylogeny-aware interface that makes evolutionary species comparisons accessible to non-experts. At its core, PlanMine continues to utilize the powerful InterMine framework and consistent data annotations to enable meaningful inter-species comparisons. Overall, PlanMine 3.0 thus provides a host of new features that makes the fascinating biology of flatworms accessible to the wider research community.

Inferring Ancient Relationships with Genomic Data: A Commentary on Current Practices

Contemporary phylogeneticists enjoy an embarrassment of riches, not only in the volumes of data now available, but also in the diversity of bioinformatic tools for handling these data. Here, I discuss a subset of these tools I consider well-suited to the task of inferring ancient relationships with coding sequence data in particular, encompassing data generation, orthology assignment, alignment and gene tree inference, supermatrix construction, and analysis under the best-fitting models applicable to large-scale datasets. Throughout, I compare and critique methods, considering both their theoretical principles and the details of their implementation, and offering practical tips on usage where appropriate. I also entertain different motivations for analyzing what are almost always originally DNA sequence data as codons, amino acids, and higher-order recodings. Although presented in a linear order, I see value in using the diversity of tools available to us to assess the sensitivity of clades of biological interest to different gene and taxon sets and analytical modes, which can be an indication of the presence of systematic error, of which a few forms remain poorly controlled by even the best available inference methods.

The Evolutionary Traceability of a Protein

Orthologs document the evolution of genes and metabolic capacities encoded in extant and ancient genomes. However, the similarity between orthologs decays with time, and ultimately it becomes insufficient to infer common ancestry. This leaves ancient gene set reconstructions incomplete and distorted to an unknown extent. Here we introduce the “evolutionary traceability” as a measure that quantifies, for each protein, the evolutionary distance beyond which the sensitivity of the ortholog search becomes limiting. Using yeast, we show that genes that were thought to date back to the last universal common ancestor are of high traceability. Their functions mostly involve catalysis, ion transport, and ribonucleoprotein complex assembly. In turn, the fraction of yeast genes whose traceability is not sufficient to infer their presence in last universal common ancestor is enriched for regulatory functions. Computing the traceabilities of genes that have been experimentally characterized as being essential for a self-replicating cell reveals that many of the genes that lack orthologs outside bacteria have low traceability. This leaves open whether their orthologs in the eukaryotic and archaeal domains have been overlooked. Looking at the example of REC8, a protein essential for chromosome cohesion, we demonstrate how a traceability-informed adjustment of the search sensitivity identifies hitherto missed orthologs in the fast-evolving microsporidia. Taken together, the evolutionary traceability helps to differentiate between true absence and nondetection of orthologs, and thus improves our understanding about the evolutionary conservation of functional protein networks. “protTrace,” a software tool for computing evolutionary traceability, is freely available at https://github.com/BIONF/protTrace.git; last accessed February 10, 2019.

Comparative genomics of the major parasitic worms

Parasitic nematodes (roundworms) and platyhelminths (flatworms) cause debilitating chronic infections of humans and animals, decimate crop production and are a major impediment to socioeconomic development. Here we report a broad comparative study of 81 genomes of parasitic and non-parasitic worms. We have identified gene family births and hundreds of expanded gene families at key nodes in the phylogeny that are relevant to parasitism. Examples include gene families that modulate host immune responses, enable parasite migration though host tissues or allow the parasite to feed. We reveal extensive lineage-specific differences in core metabolism and protein families historically targeted for drug development. From an in silico screen, we have identified and prioritized new potential drug targets and compounds for testing. This comparative genomics resource provides a much-needed boost for the research community to understand and combat parasitic worms.

Integrating Embryonic Development and Evolutionary History to Characterize Tentacle-Specific Cell Types in a Ctenophore

The origin of novel traits can promote expansion into new niches and drive speciation. Ctenophores (comb jellies) are unified by their possession of a novel cell type: the colloblast, an adhesive cell found only in the tentacles. Although colloblast-laden tentacles are fundamental for prey capture among ctenophores, some species have tentacles lacking colloblasts and others have lost their tentacles completely. We used transcriptomes from 36 ctenophore species to identify gene losses that occurred specifically in lineages lacking colloblasts and tentacles. We cross-referenced these colloblast- and tentacle-specific candidate genes with temporal RNA-Seq during embryogenesis in Mnemiopsis leidyi and found that both sets of candidates are preferentially expressed during tentacle morphogenesis. We also demonstrate significant upregulation of candidates from both data sets in the tentacle bulb of adults. Both sets of candidates were enriched for an N-terminal signal peptide and protein domains associated with secretion; among tentacle candidates we also identified orthologs of cnidarian toxin proteins, presenting tantalizing evidence that ctenophore tentacles may secrete toxins along with their adhesive. Finally, using cell lineage tracing, we demonstrate that colloblasts and neurons share a common progenitor, suggesting the evolution of colloblasts involved co-option of a neurosecretory gene regulatory network. Together these data offer an initial glimpse into the genetic architecture underlying ctenophore cell-type diversity.

Interactive online application for the prediction, ranking and prioritisation of drug targets in Schistosoma haematobium

BACKGROUND

Human schistosomiasis is a neglected tropical disease caused by parasitic worms of the genus Schistosoma that still affects some 200 million people. The mainstay of schistosomiasis control is a single drug, praziquantel. The reliance on this drug carries a risk of resistance emerging to this anthelmintic, such that research towards alternative anti-schistosomal drugs is warranted. In this context, a number of studies have employed computational approaches to prioritise proteins for investigation as drug targets, based on extensive genomic, transcriptomic and small-molecule data now available.

METHODS

Here, we established a customisable, online application for the prioritisation of drug targets and applied it, for the first time, to the entire inferred proteome of S. haematobium. This application enables selection of weighted and ranked proteins representing potential drug targets, and integrates transcriptional data, orthology and gene essentiality information as well as drug-drug target associations and chemical properties of predicted ligands.

RESULTS

Using this application, we defined 25 potential drug targets in S. haematobium that associated with approved drugs, and 3402 targets that (although they could not be linked to any compounds) are conserved among a range of socioeconomically important flatworm species and might represent targets for new trematocides.

CONCLUSIONS

The online application developed here represents an interactive, customisable, expandable and reproducible drug target ranking and prioritisation approach that should be useful for the prediction of drug targets in schistosomes and other species of parasitic worms in the future. We have demonstrated the utility of this online application by predicting potential drug targets in S. haematobium that can now be evaluated using functional genomics tools and/or small molecules, to establish whether they are indeed essential for parasite survival, and to assist in the discovery of novel anti-schistosomal compounds.

Xenacoelomorph Neuropeptidomes Reveal a Major Expansion of Neuropeptide Systems during Early Bilaterian Evolution

Neuropeptides are neurosecretory signaling molecules in protostomes and deuterostomes (together Nephrozoa). Little, however, is known about the neuropeptide complement of the sister group of Nephrozoa, the Xenacoelomorpha, which together form the Bilateria. Because members of the xenacoelomorph clades Xenoturbella, Nemertodermatida, and Acoela differ extensively in their central nervous system anatomy, the reconstruction of the xenacoelomorph and bilaterian neuropeptide complements may provide insights into the relationship between nervous system evolution and peptidergic signaling. Here, we analyzed transcriptomes of seven acoels, four nemertodermatids, and two Xenoturbella species using motif searches, similarity searches, mass spectrometry and phylogenetic analyses to characterize neuropeptide precursors and neuropeptide receptors. Our comparison of these repertoires with previously reported nephrozoan and cnidarian sequences shows that the majority of annotated neuropeptide GPCRs in cnidarians are not orthologs of specific bilaterian neuropeptide receptors, which suggests that most of the bilaterian neuropeptide systems evolved after the cnidarian–bilaterian evolutionary split. This expansion of more than 20 peptidergic systems in the stem leading to the Bilateria predates the evolution of complex nephrozoan organs and nervous system architectures. From this ancient set of neuropeptides, acoels show frequent losses that correlate with their divergent central nervous system anatomy. We furthermore detected the emergence of novel neuropeptides in xenacoelomorphs and their expansion along the nemertodermatid and acoel lineages, the two clades that evolved nervous system condensations. Together, our study provides fundamental insights into the early evolution of the bilaterian peptidergic systems, which will guide future functional and comparative studies of bilaterian nervous systems.

Toward Reducing Phylostratigraphic Errors and Biases

Phylostratigraphy is a method for estimating gene age, usually applied to large numbers of genes in order to detect nonrandom age-distributions of gene properties that could shed light on mechanisms of gene origination and evolution. However, phylostratigraphy underestimates gene age with a nonnegligible probability. The underestimation is severer for genes with certain properties, creating spurious age distributions of these properties and those correlated with these properties. Here we explore three strategies to reduce phylostratigraphic error/bias. First, we test several alternative homology detection methods (PSIBLAST, HMMER, PHMMER, OMA, and GLAM2Scan) in phylostratigraphy, but fail to find any that noticeably outperforms the commonly used BLASTP. Second, using machine learning, we look for predictors of error-prone genes to exclude from phylostratigraphy, but cannot identify reliable predictors. Finally, we remove from phylostratigraphic analysis genes exhibiting errors in simulation, which by definition minimizes error/bias if the simulation is sufficiently realistic. Using this last approach, we show that some previously reported phylostratigraphic trends (e.g., younger proteins tend to evolve more rapidly and be shorter) disappear or even reverse, reconfirming the necessity of controlling phylostratigraphic error/bias. Taken together, our analyses demonstrate that phylostratigraphic errors/biases are refractory to several potential solutions but can be controlled at least partially by the exclusion of error-prone genes identified via realistic simulations. These results are expected to stimulate the judicious use of error-aware phylostratigraphy and reevaluation of previous phylostratigraphic findings.

Improved strategy for the curation and classification of kinases, with broad applicability to other eukaryotic protein groups

Despite the substantial amount of genomic and transcriptomic data available for a wide range of eukaryotic organisms, most genomes are still in a draft state and can have inaccurate gene predictions. To gain a sound understanding of the biology of an organism, it is crucial that inferred protein sequences are accurately identified and annotated. However, this can be challenging to achieve, particularly for organisms such as parasitic worms (helminths), as most gene prediction approaches do not account for substantial phylogenetic divergence from model organisms, such as Caenorhabditis elegans and Drosophila melanogaster, whose genomes are well-curated. In this paper, we describe a bioinformatic strategy for the curation of gene families and subsequent annotation of encoded proteins. This strategy relies on pairwise gene curation between at least two closely related species using genomic and transcriptomic data sets, and is built on recent work on kinase complements of parasitic worms. Here, we discuss salient technical aspects of this strategy and its implications for the curation of protein families more generally.

Profiling G protein-coupled receptors of Fasciola hepatica identifies orphan rhodopsins unique to phylum Platyhelminthes

G protein-coupled receptors (GPCRs) are established drug targets. Despite their considerable appeal as targets for next-generation anthelmintics, poor understanding of their diversity and function in parasitic helminths has thwarted progress towards GPCR-targeted anti-parasite drugs. This study facilitates GPCR research in the liver fluke, Fasciola hepatica, by generating the first profile of GPCRs from the F. hepatica genome. Our dataset describes 147 high confidence GPCRs, representing the largest cohort of GPCRs, and the largest set of in silico ligand-receptor predictions, yet reported in any parasitic helminth. All GPCRs fall within the established GRAFS nomenclature; comprising three glutamate, 135 rhodopsin, two adhesion, five frizzled, one smoothened, and one secretin GPCR. Stringent annotation pipelines identified 18 highly diverged rhodopsins in F. hepatica that maintained core rhodopsin signatures, but lacked significant similarity with non-flatworm sequences, providing a new sub-group of potential flukicide targets. These facilitated identification of a larger cohort of 76 related sequences from available flatworm genomes, representing new members of existing groups (PROF1/Srfb, Rho-L, Rho-R, Srfa, Srfc) of flatworm-specific rhodopsins. These receptors imply flatworm specific GPCR functions, and/or co-evolution with unique flatworm ligands, and could facilitate the development of exquisitely selective anthelmintics. Ligand binding domain sequence conservation relative to deorphanised rhodopsins enabled high confidence ligand-receptor matching of seventeen receptors activated by acetylcholine, neuropeptide F/Y, octopamine or serotonin. RNA-Seq analyses showed expression of 101 GPCRs across various developmental stages, with the majority expressed most highly in the pathogenic intra-mammalian juvenile parasites. These data identify a broad complement of GPCRs in F. hepatica, including rhodopsins likely to have key functions in neuromuscular control and sensory perception, as well as frizzled and adhesion/secretin families implicated, in other species, in growth, development and reproduction. This catalogue of liver fluke GPCRs provides a platform for new avenues into our understanding of flatworm biology and anthelmintic discovery.

Chromosome Evolution in the Free-Living Flatworms: First Evidence of Intrachromosomal Rearrangements in Karyotype Evolution of Macrostomum lignano (Platyhelminthes, Macrostomida)

The free-living flatworm Macrostomum lignano is a hidden tetraploid. Its genome was formed by a recent whole genome duplication followed by chromosome fusions. Its karyotype (2n = 8) consists of a pair of large chromosomes (MLI1), which contain regions of all other chromosomes, and three pairs of small metacentric chromosomes. Comparison of MLI1 with metacentrics was performed by painting with microdissected DNA probes and fluorescent in situ hybridization of unique DNA fragments. Regions of MLI1 homologous to small metacentrics appeared to be contiguous. Besides the loss of DNA repeat clusters (pericentromeric and telomeric repeats and the 5S rDNA cluster) from MLI1, the difference between small metacentrics MLI2 and MLI4 and regions homologous to them in MLI1 were revealed. Abnormal karyotypes found in the inbred DV1/10 subline were analyzed, and structurally rearranged chromosomes were described with the painting technique, suggesting the mechanism of their origin. The revealed chromosomal rearrangements generate additional diversity, opening the way toward massive loss of duplicated genes from a duplicated genome. Our findings suggest that the karyotype of M. lignano is in the early stage of genome diploidization after whole genome duplication, and further studies on M. lignano and closely related species can address many questions about karyotype evolution in animals.

Highly diversified expansions shaped the evolution of membrane bound proteins in metazoans

The dramatic increase in membrane proteome complexity is arguably one of the most pivotal evolutionary events that underpins the origin of multicellular animals. However, the origin of a significant number of membrane families involved in metazoan development has not been clarified. In this study, we have manually curated the membrane proteomes of 22 metazoan and 2 unicellular holozoan species. We identify 123,014 membrane proteins in these 24 eukaryotic species and classify 86% of the dataset. We determine 604 functional clusters that are present from the last holozoan common ancestor (LHCA) through many metazoan species. Intriguingly, we show that more than 70% of the metazoan membrane protein families have a premetazoan origin. The data show that enzymes are more highly represented in the LHCA and expand less than threefold throughout metazoan species; in contrast to receptors that are relatively few in the LHCA but expand nearly eight fold within metazoans. Expansions related to cell adhesion, communication, immune defence, and developmental processes are shown in conjunction with emerging biological systems, such as neuronal development, cytoskeleton organization, and the adaptive immune response. This study defines the possible LHCA membrane proteome and describes the fundamental functional clusters that underlie metazoan diversity and innovation.

Conservation and diversification of small RNA pathways within flatworms

Background

Small non-coding RNAs, including miRNAs, and gene silencing mediated by RNA interference have been described in free-living and parasitic lineages of flatworms, but only few key factors of the small RNA pathways have been exhaustively investigated in a limited number of species. The availability of flatworm draft genomes and predicted proteomes allowed us to perform an extended survey of the genes involved in small non-coding RNA pathways in this phylum.

Results

Overall, findings show that the small non-coding RNA pathways are conserved in all the analyzed flatworm linages; however notable peculiarities were identified. While Piwi genes are amplified in free-living worms they are completely absent in all parasitic species. Remarkably all flatworms share a specific Argonaute family (FL-Ago) that has been independently amplified in different lineages. Other key factors such as Dicer are also duplicated, with Dicer-2 showing structural differences between trematodes, cestodes and free-living flatworms. Similarly, a very divergent GW182 Argonaute interacting protein was identified in all flatworm linages. Contrasting to this, genes involved in the amplification of the RNAi interfering signal were detected only in the ancestral free living species Macrostomum lignano. We here described all the putative small RNA pathways present in both free living and parasitic flatworm lineages.

Conclusion

These findings highlight innovations specifically evolved in platyhelminths presumably associated with novel mechanisms of gene expression regulation mediated by small RNA pathways that differ to what has been classically described in model organisms. Understanding these phylum-specific innovations and the differences between free living and parasitic species might provide clues to adaptations to parasitism, and would be relevant for gene-silencing technology development for parasitic flatworms that infect hundreds of million people worldwide.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1061-5) contains supplementary material, which is available to authorized users.

Evolution of the EGFR pathway in Metazoa and its diversification in the planarian Schmidtea mediterranea

The EGFR pathway is an essential signaling system in animals, whose core components are the epidermal growth factors (EGF ligands) and their trans-membrane tyrosine kinase receptors (EGFRs). Despite extensive knowledge in classical model organisms, little is known of the composition and function of the EGFR pathway in most animal lineages. Here, we have performed an extensive search for the presence of EGFRs and EGF ligands in representative species of most major animal clades, with special focus on the planarian Schmidtea mediterranea. With the exception of placozoans and cnidarians, we found that the EGFR pathway is potentially present in all other analyzed animal groups, and has experienced frequent independent expansions. We further characterized the expression domains of the EGFR/EGF identified in S. mediterranea, revealing a wide variety of patterns and localization in almost all planarian tissues. Finally, functional experiments suggest an interaction between one of the previously described receptors, Smed-egfr-5, and the newly found ligand Smed-egf-6. Our findings provide the most comprehensive overview to date of the EGFR pathway, and indicate that the last common metazoan ancestor had an initial complement of one EGFR and one putative EGF ligand, which was often expanded or lost during animal evolution.