The evolution of animal genomes

Phenogenetics of cortical granule dynamics during zebrafish oocyte-to-embryo transition

Fertilization is a critical process in sexual reproduction that involves the fusion of a capacitated sperm with a mature oocyte to form a zygote. Polyspermy, the fertilization of an oocyte by multiple sperm, leads to polyploidy and embryo lethality. Mammalian and non-mammalian oocytes have evolved mechanisms to prevent polyspermy, including fast and slow blocks. The fast block comprises membrane depolarization post-sperm fusion, temporarily preventing additional sperm fusion. The slow block, triggered by cortical granule (CG) exocytosis, involves the release of proteins that modify the zona pellucida to form a permanent barrier, avoiding the fertilization by additional sperm. The evidence shows that immature oocytes often fail to prevent polyspermy due to ineffective CG exocytosis, attributed to impaired intracellular calcium increases, lower content of this ion, and incomplete CG migration. The study of how genetic variations lead to observable phenotypes (phenogenetics) during the oocyte-to-embryo transition, have identified several maternal-effect genes in zebrafish involved in CG behavior. These genes regulate various stages of CG biology, including biosynthesis, maturation, and exocytosis. Mutations in these genes disrupt these processes, highlighting the maternal genetic control over CG properties. Zebrafish has emerged as a pivotal model for understanding the evolving genetic regulation and molecular mechanisms underlying CG biology, providing valuable insights into fertility and early embryonic development.

Morphological and transcriptional effects of crude oil and dispersant exposure on the marine sponge Cinachyrella alloclada

Marine sponges play important roles in benthic ecosystems. More than providing shelter and food to other species, they help maintain water quality by regulating nitrogen and ammonium levels in the water, and bioaccumulate heavy metals. This system, however, is particularly sensitive to sudden environmental changes including catastrophic pollution event such as oil spills. Hundreds of oil platforms are currently actively extracting oil and gas in the Gulf of Mexico. To test the vulnerability of the benthic ecosystems to oil spills, we utilized the Caribbean reef sponge, Cinachyrella alloclada, as a novel experimental indicator. We have exposed organisms to crude oil and oil dispersant for up to 24 h and measured resultant gene expression changes. Our findings indicate that 1-hour exposure to water accommodated fractions (WAF) was enough to elicit massive shifts in gene expression in sponges and host bacterial communities (8052 differentially expressed transcripts) with the up-regulation of stress related pathways, cancer related pathways, and cell integrity pathways. Genes that were upregulated included heat shock proteins, apoptosis, oncogenes (Rab/Ras, Src, CMYC), and several E3 ubiquitin ligases. 24-hour exposure of chemically enhanced WAF (CE-WAF) had the greatest impact to benthic communities, resulting in mostly downregulation of gene expression (4248 differentially expressed transcripts). Gene deregulation from 1-hour treatments follow this decreasing trend of toxicity: WAF > CE-WAF > Dispersant, while the 24-hour treatment showed a shift to CE-WAF > Dispersant > WAF in our experiments. Thus, this study supports the development of Cinachyrella alloclada as a research model organism and bioindicator species for Florida reefs and underscores the importance of developing more efficient and safer ways to remove oil in the event of a spill catastrophe.

The compact genome of the sponge Oopsacas minuta (Hexactinellida) is lacking key metazoan core genes

BACKGROUND

Explaining the emergence of the hallmarks of bilaterians is a central focus of evolutionary developmental biology-evodevo-and evolutionary genomics. For this purpose, we must both expand and also refine our knowledge of non-bilaterian genomes, especially by studying early branching animals, in particular those in the metazoan phylum Porifera.

RESULTS

We present a comprehensive analysis of the first whole genome of a glass sponge, Oopsacas minuta, a member of the Hexactinellida. Studying this class of sponge is evolutionary relevant because it differs from the three other Porifera classes in terms of development, tissue organization, ecology, and physiology. Although O. minuta does not exhibit drastic body simplifications, its genome is among the smallest of animal genomes sequenced so far, and surprisingly lacks several metazoan core genes (including Wnt and several key transcription factors). Our study also provides the complete genome of a symbiotic Archaea dominating the associated microbial community: a new Thaumarchaeota species.

CONCLUSIONS

The genome of the glass sponge O. minuta differs from all other available sponge genomes by its compactness and smaller number of encoded proteins. The unexpected loss of numerous genes previously considered ancestral and pivotal for metazoan morphogenetic processes most likely reflects the peculiar syncytial tissue organization in this group. Our work further documents the importance of convergence during animal evolution, with multiple convergent evolution of septate-like junctions, electrical-signaling and multiciliated cells in metazoans.

Unsupervised Deep Learning Can Identify Protein Functional Groups from Unaligned Sequences

Interpreting protein function from sequence data is a fundamental goal of bioinformatics. However, our current understanding of protein diversity is bottlenecked by the fact that most proteins have only been functionally validated in model organisms, limiting our understanding of how function varies with gene sequence diversity. Thus, accuracy of inferences in clades without model representatives is questionable. Unsupervised learning may help to ameliorate this bias by identifying highly complex patterns and structure from large datasets without external labels. Here we present DeepSeqProt, an unsupervised deep learning program for exploring large protein sequence datasets. DeepSeqProt is a clustering tool capable of distinguishing between broad classes of proteins while learning local and global structure of functional space. DeepSeqProt is capable of learning salient biological features from unaligned, unannotated sequences. DeepSeqProt is more likely to capture complete protein families and statistically significant shared ontologies within proteomes than other clustering methods. We hope this framework will prove of use to researchers and provide a preliminary step in further developing unsupervised deep learning in molecular biology.

Induced Immune Reaction in the Acorn Worm, Saccoglossus kowalevskii, Informs the Evolution of Antiviral Immunity

Evolutionary perspectives on the deployment of immune factors following infection have been shaped by studies on a limited number of biomedical model systems with a heavy emphasis on vertebrate species. Although their contributions to contemporary immunology cannot be understated, a broader phylogenetic perspective is needed to understand the evolution of immune systems across Metazoa. In our study, we leverage differential gene expression analyses to identify genes implicated in the antiviral immune response of the acorn worm hemichordate, Saccoglossus kowalevskii, and place them in the context of immunity evolution within deuterostomes-the animal clade composed of chordates, hemichordates, and echinoderms. Following acute exposure to the synthetic viral double-stranded RNA analog, poly(I:C), we show that S. kowalevskii responds by regulating the transcription of genes associated with canonical innate immunity signaling pathways (e.g., nuclear factor κB and interferon regulatory factor signaling) and metabolic processes (e.g., lipid metabolism), as well as many genes without clear evidence of orthology with those of model species. Aggregated across all experimental time point contrasts, we identify 423 genes that are differentially expressed in response to poly(I:C). We also identify 147 genes with altered temporal patterns of expression in response to immune challenge. By characterizing the molecular toolkit involved in hemichordate antiviral immunity, our findings provide vital evolutionary context for understanding the origins of immune systems within Deuterostomia.

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.

Phylogenomics and the first higher taxonomy of Placozoa, an ancient and enigmatic animal phylum

Placozoa is an ancient phylum of extraordinarily unusual animals: miniscule, ameboid creatures that lack most fundamental animal features. Despite high genetic diversity, only recently have the second and third species been named. While prior genomic studies suffer from incomplete placozoan taxon sampling, we more than double the count with protein sequences from seven key genomes and produce the first nuclear phylogenomic reconstruction of all major placozoan lineages. This leads us to the first complete Linnaean taxonomic classification of Placozoa, over a century after its discovery: This may be the only time in the 21st century when an entire higher taxonomy for a whole animal phylum is formalized. Our classification establishes 2 new classes, 4 new orders, 3 new families, 1 new genus, and 1 new species, namely classes Polyplacotomia and Uniplacotomia; orders Polyplacotomea, Trichoplacea, Cladhexea, and Hoilungea; families Polyplacotomidae, Cladtertiidae, and Hoilungidae; and genus Cladtertia with species Cladtertia collaboinventa, nov. Our likelihood and gene content tree topologies refine the relationships determined in previous studies. Adding morphological data into our phylogenomic matrices suggests sponges (Porifera) as the sister to other animals, indicating that modest data addition shifts this node away from comb jellies (Ctenophora). Furthermore, by adding the first genomic protein data of the exceptionally distinct and branching Polyplacotoma mediterranea, we solidify its position as sister to all other placozoans; a divergence we estimate to be over 400 million years old. Yet even this deep split sits on a long branch to other animals, suggesting a bottleneck event followed by diversification. Ancestral state reconstructions indicate large shifts in gene content within Placozoa, with Hoilungia hongkongensis and its closest relatives having the most unique genetics.

A missense mutation in ISPD contributes to maintain muscle fiber stability

Background

Results

Conclusion

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.

Knockin' on Egg's Door: Maternal Control of Egg Activation That Influences Cortical Granule Exocytosis in Animal Species

Fertilization by multiple sperm leads to lethal chromosomal number abnormalities, failed embryo development, and miscarriage. In some vertebrate and invertebrate eggs, the so-called cortical reaction contributes to their activation and prevents polyspermy during fertilization. This process involves biogenesis, redistribution, and subsequent accumulation of cortical granules (CGs) at the female gamete cortex during oogenesis. CGs are oocyte- and egg-specific secretory vesicles whose content is discharged during fertilization to block polyspermy. Here, we summarize the molecular mechanisms controlling critical aspects of CG biology prior to and after the gametes interaction. This allows to block polyspermy and provide protection to the developing embryo. We also examine how CGs form and are spatially redistributed during oogenesis. During egg activation, CG exocytosis (CGE) and content release are triggered by increases in intracellular calcium and relies on the function of maternally-loaded proteins. We also discuss how mutations in these factors impact CG dynamics, providing unprecedented models to investigate the genetic program executing fertilization. We further explore the phylogenetic distribution of maternal proteins and signaling pathways contributing to CGE and egg activation. We conclude that many important biological questions and genotype–phenotype relationships during fertilization remain unresolved, and therefore, novel molecular players of CG biology need to be discovered. Future functional and image-based studies are expected to elucidate the identity of genetic candidates and components of the molecular machinery involved in the egg activation. This, will open new therapeutic avenues for treating infertility in humans.

Dissecting the chromosome-level genome of the Asian Clam (Corbicula fluminea)

The Asian Clam (Corbicula fluminea) is a valuable commercial and medicinal bivalve, which is widely distributed in East and Southeast Asia. As a natural nutrient source, the clam is rich in protein, amino acids, and microelements. The genome of C. fluminea has not yet been characterized; therefore, genome-assisted breeding and improvements cannot yet be implemented. In this work, we present a de novo chromosome-scale genome assembly of C. fluminea using PacBio and Hi-C sequencing technologies. The assembled genome comprised 4728 contigs, with a contig N50 of 521.06 Kb, and 1,215 scaffolds with a scaffold N50 of 70.62 Mb. More than 1.51 Gb (99.17%) of genomic sequences were anchored to 18 chromosomes, of which 1.40 Gb (92.81%) of genomic sequences were ordered and oriented. The genome contains 38,841 coding genes, 32,591 (83.91%) of which were annotated in at least one functional database. Compared with related species, C. fluminea had 851 expanded gene families and 191 contracted gene families. The phylogenetic tree showed that C. fluminea diverged from Ruditapes philippinarum, ~ 228.89 million years ago (Mya), and the genomes of C. fluminea and R. philippinarum shared 244 syntenic blocks. Additionally, we identified 2 MITF members and 99 NLRP members in C. fluminea genome. The high-quality and chromosomal Asian Clam genome will be a valuable resource for a range of development and breeding studies of C. fluminea in future research.

A stony coral cell atlas illuminates the molecular and cellular basis of coral symbiosis, calcification, and immunity

Stony corals are colonial cnidarians that sustain the most biodiverse marine ecosystems on Earth: coral reefs. Despite their ecological importance, little is known about the cell types and molecular pathways that underpin the biology of reef-building corals. Using single-cell RNA sequencing, we define over 40 cell types across the life cycle of Stylophora pistillata. We discover specialized immune cells, and we uncover the developmental gene expression dynamics of calcium-carbonate skeleton formation. By simultaneously measuring the transcriptomes of coral cells and the algae within them, we characterize the metabolic programs involved in symbiosis in both partners. We also trace the evolution of these coral cell specializations by phylogenetic integration of multiple cnidarian cell type atlases. Overall, this study reveals the molecular and cellular basis of stony coral biology.

Evolutionary Cell Type Mapping with Single-Cell Genomics

A fundamental characteristic of animal multicellularity is the spatial coexistence of functionally specialized cell types that are all encoded by a single genome sequence. Cell type transcriptional programs are deployed and maintained by regulatory mechanisms that control the asymmetric, differential access to genomic information in each cell. This genome regulation ultimately results in specific cellular phenotypes. However, the emergence, diversity, and evolutionary dynamics of animal cell types remain almost completely unexplored beyond a few species. Single-cell genomics is emerging as a powerful tool to build comprehensive catalogs of cell types and their associated gene regulatory programs in non-traditional model species. We review the current state of sampling efforts across the animal tree of life and challenges ahead for the comparative study of cell type programs. We also discuss how the phylogenetic integration of cell atlases can lead to the development of models of cell type evolution and a phylogenetic taxonomy of cells.

Symmetry Transformations in Metazoan Evolution and Development

In this review, we consider transformations of axial symmetry in metazoan evolution and development, the genetic basis, and phenotypic expressions of different axial body plans. In addition to the main symmetry types in metazoan body plans, such as rotation (radial symmetry), reflection (mirror and glide reflection symmetry), and translation (metamerism), many biological objects show scale (fractal) symmetry as well as some symmetry-type combinations. Some genetic mechanisms of axial pattern establishment, creating a coordinate system of a metazoan body plan, bilaterian segmentation, and left–right symmetry/asymmetry, are analysed. Data on the crucial contribution of coupled functions of the Wnt, BMP, Notch, and Hedgehog signaling pathways (all pathways are designated according to the abbreviated or full names of genes or their protein products; for details, see below) and the axial Hox-code in the formation and maintenance of metazoan body plans are necessary for an understanding of the evolutionary diversification and phenotypic expression of various types of axial symmetry. The lost body plans of some extinct Ediacaran and early Cambrian metazoans are also considered in comparison with axial body plans and posterior growth in living animals.

Draft genome of Bugula neritina, a colonial animal packing powerful symbionts and potential medicines

Many animal phyla have no representatives within the catalog of whole metazoan genome sequences. This dataset fills in one gap in the genome knowledge of animal phyla with a draft genome of Bugula neritina (phylum Bryozoa). Interest in this species spans ecology and biomedical sciences because B. neritina is the natural source of bioactive compounds called bryostatins. Here we present a draft assembly of the B. neritina genome obtained from PacBio and Illumina HiSeq data, as well as genes and proteins predicted de novo and verified using transcriptome data, along with the functional annotation. These sequences will permit a better understanding of host-symbiont interactions at the genomic level, and also contribute additional phylogenomic markers to evaluate Lophophorate or Lophotrochozoa phylogenetic relationships. The effort also fits well with plans to ultimately sequence all orders of the Metazoa.

Measurement(s)	DNA • genome • sequence_assembly • sequence feature annotation
Technology Type(s)	DNA sequencing • sequence assembly process • sequence annotation
Sample Characteristic - Organism	Bugula neritina

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12988355

The biology of Lonp1: More than a mitochondrial protease

Initially discovered as a protease responsible for degradation of misfolded or damaged proteins, the mitochondrial Lon protease (Lonp1) turned out to be a multifaceted enzyme, that displays at least three different functions (proteolysis, chaperone activity, binding of mtDNA) and that finely regulates several cellular processes, within and without mitochondria. Indeed, LONP1 in humans is ubiquitously expressed, and is involved in regulation of response to oxidative stress and, heat shock, in the maintenance of mtDNA, in the regulation of mitophagy. Furthermore, its proteolytic activity can regulate several biochemical pathways occurring totally or partially within mitochondria, such as TCA cycle, oxidative phosphorylation, steroid and heme biosynthesis and glutamine production. Because of these multiple activities, Lon protease is highly conserved throughout evolution, and mutations occurring in its gene determines severe diseases in humans, including a rare syndrome characterized by Cerebral, Ocular, Dental, Auricular and Skeletal anomalies (CODAS). Finally, alterations of LONP1 regulation in humans can favor tumor progression and aggressiveness, further highlighting the crucial role of this enzyme in mitochondrial and cellular homeostasis.

Sequencing Disparity in the Genomic Era

Advances in sequencing technology have resulted in the expectation that genomic studies will become more representative of organismal diversity. To test this expectation, we explored species representation of nonhuman eukaryotes in the Sequence Read Archive. Though species richness has been increasing steadily, species evenness is decreasing over time. Moreover, the top 1% most studied organisms increasingly represent a larger proportion of total experiments, demonstrating growing bias in favor of a small minority of species. To better understand molecular processes and patterns, genomic studies should reverse current trends by adopting more comparative approaches.

The Genome and mRNA Transcriptome of the Cosmopolitan Calanoid Copepod Acartia tonsa Dana Improve the Understanding of Copepod Genome Size Evolution

Members of the crustacean subclass Copepoda are likely the most abundant metazoans worldwide. Pelagic marine species are critical in converting planktonic microalgae to animal biomass, supporting oceanic food webs. Despite their abundance and ecological importance, only six copepod genomes are publicly available, owing to a number of factors including large genome size, repetitiveness, GC-content, and small animal size. Here, we report the seventh representative copepod genome and the first genome and the first transcriptome from the calanoid copepod species Acartia tonsa Dana, which is among the most numerous mesozooplankton in boreal coastal and estuarine waters. The ecology, physiology, and behavior of A. tonsa have been studied extensively. The genetic resources contributed in this work will allow researchers to link experimental results to molecular mechanisms. From PCR-free whole genome sequence and mRNA Illumina data, we assemble the largest copepod genome to date. We estimate that A. tonsa has a total genome size of 2.5 Gb including repetitive elements we could not resolve. The nonrepetitive fraction of the genome assembly is estimated to be 566 Mb. Our DNA sequencing-based analyses suggest there is a 14-fold difference in genome size between the six members of Copepoda with available genomic information. This finding complements nucleus staining genome size estimations, where 100-fold difference has been reported within 70 species. We briefly analyze the repeat structure in the existing copepod whole genome sequence data sets. The information presented here confirms the evolution of genome size in Copepoda and expands the scope for evolutionary inferences in Copepoda by providing several levels of genetic information from a key planktonic crustacean species.

Catecholic Compounds in Ctenophore Colloblast and Nerve Net Proteins Suggest a Structural Role for DOPA-Like Molecules in an Early-Diverging Animal Lineage

Ctenophores, or comb jellies, are among the earliest-diverging extant animal lineages. Several recent phylogenomic studies suggest that they may even be the sister group to all other animals. This unexpected finding remains difficult to contextualize, particularly given ctenophores' unique and sometimes poorly understood physiology. Colloblasts, a ctenophore-specific cell type found on the surface of these animals' tentacles, are emblematic of this difficulty. The exterior of the colloblast is dotted with granules that burst and release an adhesive on contact with prey, ensnaring it for consumption. To date, little is known about the fast-acting underwater adhesive that these cells secrete or its biochemistry. We present evidence that proteins in the colloblasts of the ctenophore Pleurobrachia bachei incorporate catecholic compounds similar to the amino acid l-3,4-dihydroxyphenylalanine. These compounds are associated with adhesive-containing granules on the surface of colloblasts, suggesting that they may play a role in prey capture, akin to dihydroxyphenylalanine-based adhesives in mussel byssus. We also present unexpected evidence of similar catecholic compounds in association with the subepithelial nerve net. There, catecholic compounds are present in spatial patterns similar to those of l-3,4-dihydroxyphenylalanine and its derivatives in cnidarian nerves, where they are associated with membranes and possess unknown functionality. This "structural" use of catecholic molecules in ctenophores represents the earliest-diverging animal lineage in which this trait has been observed, though it remains unclear whether structural catechols are deeply rooted in animals or whether they have arisen multiple times.

Multiple Facets of Marine Invertebrate Conservation Genomics

Conservation genomics aims to preserve the viability of populations and the biodiversity of living organisms. Invertebrate organisms represent 95% of animal biodiversity; however, few genomic resources currently exist for the group. The subset of marine invertebrates includes the most ancient metazoan lineages and possesses codes for unique gene products and possible keys to adaptation. The benefits of supporting invertebrate conservation genomics research (e.g., likely discovery of novel genes, protein regulatory mechanisms, genomic innovations, and transposable elements) outweigh the various hurdles (rare, small, or polymorphic starting materials). Here we review best conservation genomics practices in the laboratory and in silico when applied to marine invertebrates and also showcase unique features in several case studies of acroporid corals, crown-of-thorns starfish, apple snails, and abalone. Marine conservation genomics should also address how diversity can lead to unique marine innovations, the impact of deleterious variation, and how genomic monitoring and profiling could positively affect broader conservation goals (e.g., value of baseline data for in situ/ex situ genomic stocks).

Getting Nervous: An Evolutionary Overhaul for Communication

The evolution of a nervous system as a control system of the body's functions is a key innovation of animals. Its fundamental units are neurons, highly specialized cells dedicated to fast cell-cell communication. Neurons pass signals to other neurons, muscle cells, or gland cells at specialized junctions, the synapses, where transmitters are released from vesicles in a Ca²⁺-dependent fashion to activate receptors in the membrane of the target cell. Reconstructing the origins of neuronal communication out of a more simple process remains a central challenge in biology. Recent genomic comparisons have revealed that all animals, including the nerveless poriferans and placozoans, share a basic set of genes for neuronal communication. This suggests that the first animal, the Urmetazoan, was already endowed with neurosecretory cells that probably started to connect into neuronal networks soon afterward. Here, we discuss scenarios for this pivotal transition in animal evolution.

Advancing Genomics through the Global Invertebrate Genomics Alliance (GIGA)

The Global Invertebrate Genomics Alliance (GIGA), a collaborative network of diverse scientists, marked its second anniversary with a workshop in Munich, Germany, where international attendees focused on discussing current progress, milestones and bioinformatics resources. The community determined the recruitment and training talented researchers as one of the most pressing future needs and identified opportunities for network funding. GIGA also promotes future research efforts to prioritize taxonomic diversity and create new synergies. Here, we announce the generation of a central and simple data repository portal with a wide coverage of available sequence data, via the compagen platform, in parallel with more focused and specialized organism databases to globally advance invertebrate genomics. Therefore this article serves the objectives of GIGA by disseminating current progress and future prospects in the science of invertebrate genomics with the aim of promotion and facilitation of interdisciplinary and international research.

An annotated draft genome for Radix auricularia (Gastropoda, Mollusca)

Molluscs are the second most species-rich phylum in the animal kingdom, yet only 11 genomes of this group have been published so far. Here, we present the draft genome sequence of the pulmonate freshwater snail Radix auricularia . Six whole genome shotgun libraries with different layouts were sequenced. The resulting assembly comprises 4,823 scaffolds with a cumulative length of 910 Mb and an overall read coverage of 72×. The assembly contains 94.6% of a metazoan core gene collection, indicating an almost complete coverage of the coding fraction. The discrepancy of ∼690 Mb compared with the estimated genome size of R. auricularia (1.6 Gb) results from a high repeat content of 70% mainly comprising DNA transposons. The annotation of 17,338 protein coding genes was supported by the use of publicly available transcriptome data. This draft will serve as starting point for further genomic and population genetic research in this scientifically important phylum.

Animal Mitochondrial DNA as We Do Not Know It: mt-Genome Organization and Evolution in Nonbilaterian Lineages

Animal mitochondrial DNA (mtDNA) is commonly described as a small, circular molecule that is conserved in size, gene content, and organization. Data collected in the last decade have challenged this view by revealing considerable diversity in animal mitochondrial genome organization. Much of this diversity has been found in nonbilaterian animals (phyla Cnidaria, Ctenophora, Placozoa, and Porifera), which, from a phylogenetic perspective, form the main branches of the animal tree along with Bilateria. Within these groups, mt-genomes are characterized by varying numbers of both linear and circular chromosomes, extra genes (e.g. atp9, polB, tatC), large variation in the number of encoded mitochondrial transfer RNAs (tRNAs) (0-25), at least seven different genetic codes, presence/absence of introns, tRNA and mRNA editing, fragmented ribosomal RNA genes, translational frameshifting, highly variable substitution rates, and a large range of genome sizes. This newly discovered diversity allows a better understanding of the evolutionary plasticity and conservation of animal mtDNA and provides insights into the molecular and evolutionary mechanisms shaping mitochondrial genomes.

Controlling for Phylogenetic Relatedness and Evolutionary Rates Improves the Discovery of Associations Between Species' Phenotypic and Genomic Differences

The growing number of sequenced genomes allows us now to address a key question in genetics and evolutionary biology: which genomic changes underlie particular phenotypic changes between species? Previously, we developed a computational framework called Forward Genomics that associates phenotypic to genomic differences by focusing on phenotypes that are independently lost in different lineages. However, our previous implementation had three main limitations. Here, we present two new Forward Genomics methods that overcome these limitations by (1) directly controlling for phylogenetic relatedness, (2) controlling for differences in evolutionary rates, and (3) computing a statistical significance. We demonstrate on large-scale simulated data and on real data that both new methods substantially improve the sensitivity to detect associations between phenotypic and genomic differences. We applied these new methods to detect genomic differences involved in the loss of vision in the blind mole rat and the cape golden mole, two independent subterranean mammals. Forward Genomics identified several genes that are enriched in functions related to eye development and the perception of light, as well as genes involved in the circadian rhythm. These new Forward Genomics methods represent a significant advance in our ability to discover the genomic basis underlying phenotypic differences between species. Source code: https://github.com/hillerlab/ForwardGenomics/.

Serum amyloid A in marine bivalves: An acute phase and innate immunity protein

Serum amyloid A (SAA) is among the most potent acute phase proteins (APP) in vertebrates. After injury, its early expression can dramatically increase to promote the recruitment of immuno-competent cells, expression of pro-inflammatory proteins and the activation of the innate immune defences. Although APP have been studied in many vertebrates, only recently their search was extended to invertebrates and the finding of SAA-like molecules has opened new questions on the immune-regulatory functions of these soluble proteins in the animal kingdom. Taking advantage of the considerable amount of genomic and transcriptomic data currently available, we retrieved 51 SAA-like proteins in several protostome taxa comprising 21 marine bivalve species and basal metazoans. In addition to vertebrate-like SAAs, we identified a second protein type with peculiar features. In the bivalves Crassostrea gigas and Mytilus galloprovincialis, both digital expression analysis and qPCR data indicated an induction of the classical SAA after bacterial challenge.

Predicting RAD-seq Marker Numbers across the Eukaryotic Tree of Life

High-throughput sequencing of reduced representation libraries obtained through digestion with restriction enzymes--generically known as restriction site associated DNA sequencing (RAD-seq)--is a common strategy to generate genome-wide genotypic and sequence data from eukaryotes. A critical design element of any RAD-seq study is knowledge of the approximate number of genetic markers that can be obtained for a taxon using different restriction enzymes, as this number determines the scope of a project, and ultimately defines its success. This number can only be directly determined if a reference genome sequence is available, or it can be estimated if the genome size and restriction recognition sequence probabilities are known. However, both scenarios are uncommon for nonmodel species. Here, we performed systematic in silico surveys of recognition sequences, for diverse and commonly used type II restriction enzymes across the eukaryotic tree of life. Our observations reveal that recognition sequence frequencies for a given restriction enzyme are strikingly variable among broad eukaryotic taxonomic groups, being largely determined by phylogenetic relatedness. We demonstrate that genome sizes can be predicted from cleavage frequency data obtained with restriction enzymes targeting "neutral" elements. Models based on genomic compositions are also effective tools to accurately calculate probabilities of recognition sequences across taxa, and can be applied to species for which reduced representation data are available (including transcriptomes and neutral RAD-seq data sets). The analytical pipeline developed in this study, PredRAD (https://github.com/phrh/PredRAD), and the resulting databases constitute valuable resources that will help guide the design of any study using RAD-seq or related methods.