Insights into bilaterian evolution from three spiralian genomes

Exploring the innate immune system of Urechis unicinctus: Insights from full-length transcriptome analysis

The Echiura worm Urechis unicinctus refers to a common benthic invertebrate found in the intertidal zone of Huanghai as well as Bohai Bay. U. unicinctus is known to contain various physiologically active substances, making it highly valuable in terms of its edibility, medicinal properties, and economic potential. Nonetheless, the limited study on the immune system of U. unicinctus poses difficulties for its aquaculture and artificial reproduction. Marine invertebrates, including shellfish and U. unicinctus, are thought to primarily depend on their innate immune system for disease protection, owing to the severalinnate immune molecules they possess. Herein, we employed PacBio single-molecule real-time (SMRT) sequencing technology to perform the full-length transcriptome analysis of U. unicinctus individuals under five different conditions (room temperature (RT), low temperature (LT), high temperature (HT), without water (DRY), ultraviolet irradiation (UV)). Concequently, we identified 59,371 unigenes that had a 2,779 bp average length, 2,613 long non-coding RNAs (lncRNAs), 59,190 coding sequences (CDSs), 35,166 simple sequence repeats (SSRs), and 1,733 transcription factors (TFs), successfully annotating 90.58 % (53,778) of the unigenes. Subsequently, key factors associated with immune-related processes, such as non-self-recognition, cellular immune defenses, and humoral immune defenses, were searched. Our study also identified pattern recognition receptors (PRRs) that included 17 peptidoglycan recognition proteins (PGRPs), 13 Gram-negative binding proteins (GNBPs), 18 scavenger receptors (SRs), 74 toll-like receptors (TLRs), and 89 C-type lectins (CLTs). Altogether, the high-quality transcriptome obtained data will offer valuable insights for further investigations into U. unicinctus innate immune response, laying the foundation for subsequent molecular biology studies and aquaculture.

The energy metabolism of the freshwater leech Whitmania pigra in response to fasting

Non-blood-feeding leeches, Whitmania pigra, have evolved unique digestive structures and physiological mechanisms to cope with fasting. However, the metabolic changes and molecular mechanisms induced by fasting remain unclear. Therefore, this study recorded the weights of leeches during the fasting process. The weight changes were divided into two stages: a rapid decline period (1-9 weeks) and a fluctuating decline period (9-24 weeks). Leeches fasted for 4 (H4), 11 (H11), and 24 (H24) weeks were selected for transcriptome sequencing. Compared to the control group (H0), 436, 1157, and 337 differentially expressed genes (DEGs) were identified, which were mainly related to glycolysis/gluconeogenesis, amino acid metabolism, and the lipid metabolism pathway. The 6-phosphofructokinase (Pfk), pyruvate kinase (PK), and phosphoenolpyruvate carboxykinase (Pck) transcription levels revealed glycolysis/gluconeogenesis activation during the early stage of fasting and peaked at 11 weeks. Decreased expression of the rate-limiting enzyme acetyl-CoA carboxylase (ACC) in fatty acid synthesis during fasting may impede fatty acid synthesis. These results indicated that the nutrient storage and energy-supplying pathways in W. pigra were modified to improve fasting resistance. The findings of this study provided guidance for exploring the mechanism underlying fasting metabolism and laid a foundation for artificial breeding to improve the resistance of leeches.

Annelid Comparative Genomics and the Evolution of Massive Lineage-Specific Genome Rearrangement in Bilaterians

The organization of genomes into chromosomes is critical for processes such as genetic recombination, environmental adaptation, and speciation. All animals with bilateral symmetry inherited a genome structure from their last common ancestor that has been highly conserved in some taxa but seemingly unconstrained in others. However, the evolutionary forces driving these differences and the processes by which they emerge have remained largely uncharacterized. Here, we analyze genome organization across the phylum Annelida using 23 chromosome-level annelid genomes. We find that while many annelid lineages have maintained the conserved bilaterian genome structure, the Clitellata, a group containing leeches and earthworms, possesses completely scrambled genomes. We develop a rearrangement index to quantify the extent of genome structure evolution and show that, compared to the last common ancestor of bilaterians, leeches and earthworms have among the most highly rearranged genomes of any currently sampled species. We further show that bilaterian genomes can be classified into two distinct categories-high and low rearrangement-largely influenced by the presence or absence, respectively, of chromosome fission events. Our findings demonstrate that animal genome structure can be highly variable within a phylum and reveal that genome rearrangement can occur both in a gradual, stepwise fashion, or rapid, all-encompassing changes over short evolutionary timescales.

Diversification of von Willebrand Factor A and Chitin-Binding Domains in Pif/BMSPs Among Mollusks

Pif is a shell matrix protein (SMP) identified in the nacreous layer of Pinctada fucata (Pfu) comprised two proteins, Pif97 and Pif 80. Pif97 contains a von Willebrand factor A (VWA) and chitin-binding domains, whereas Pif80 can bind calcium carbonate crystals. The VWA domain is conserved in the SMPs of various mollusk species; however, their phylogenetic relationship remains obscure. Furthermore, although the VWA domain participates in protein-protein interactions, its role in shell formation has not been established. Accordingly, in the current study, we investigate the phylogenetic relationship between PfuPif and other VWA domain-containing proteins in major mollusk species. The shell-related proteins containing VWA domains formed a large clade (the Pif/BMSP family) and were classified into eight subfamilies with unique sequential features, expression patterns, and taxa diversity. Furthermore, a pull-down assay using recombinant proteins containing the VWA domain of PfuPif 97 revealed that the VWA domain interacts with five nacreous layer-related SMPs of P. fucata, including Pif 80 and nacrein. Collectively, these results suggest that the VWA domain is important in the formation of organic complexes and participates in shell mineralisation.

Is there potential for estradiol receptor signaling in lophotrochozoans?

Estrogen receptors (ERs) are thought to be the ancestor of all steroid receptors and are present in most lophotrochozoans studied to date, including molluscs, annelids, and rotifers. A number of studies have investigated the functional role of estrogen receptors in invertebrate species, although most are in molluscs, where the receptor is constitutively active. In vitro experiments provided evidence for ligand-activated estrogen receptors in annelids, raising important questions about the role of estrogen signalling in lophotrochozoan lineages. Here, we review the concordant and discordant evidence of estradiol receptor signalling in lophotrochozoans, with a focus on annelids and rotifers. We explore the de novo synthesis of estrogens, the evolution and expression of estrogen receptors, and physiological responses to activation of estrogen receptors in the lophotrochozoan phyla Annelida and Rotifera. Key data are missing to determine if de novo biosynthesis of estradiol in non-molluscan lophotrochozoans is likely. For example, an ortholog for the CYP11 gene is present, but confirmation of substrate conversion and measured tissue products is lacking. Orthologs CYP17 and CYP19 are lacking, yet intermediates or products (e.g. estradiol) in tissues have been measured. Estrogen receptors are present in multiple species, and for a limited number, in vitro data show agonist binding of estradiol and/or transcriptional activation. The expression patterns of the lophotrochozoan ERs suggest developmental, reproductive, and digestive roles but are highly species dependent. E2 exposures suggest that lophotrochozoan ERs may play a role in reproduction, but no strong dose-response relationship has been established. Therefore, we expect most lophotrochozoan species, outside of perhaps platyhelminths, to have an ER but their physiological role remains elusive. Mining genomes for orthologs gene families responsible for steroidogenesis, coupled with in vitro and in vivo studies of the steroid pathway are needed to better assess whether lophotrochozoans are capable of estradiol biosynthesis. One major challenge is that much of the data are divided across a diversity of species. We propose that the polychaetes Capitella teleta or Platyneris dumerilii, and rotifer Brachionus manjavacas may be strong species choices for studies of estrogen receptor signalling, because of available genomic data, established laboratory culture techniques, and gene knockout potential.

Annelid methylomes reveal ancestral developmental and aging-associated epigenetic erosion across Bilateria

BACKGROUND

DNA methylation in the form of 5-methylcytosine (5mC) is the most abundant base modification in animals. However, 5mC levels vary widely across taxa. While vertebrate genomes are hypermethylated, in most invertebrates, 5mC concentrates on constantly and highly transcribed genes (gene body methylation; GbM) and, in some species, on transposable elements (TEs), a pattern known as "mosaic". Yet, the role and developmental dynamics of 5mC and how these explain interspecies differences in DNA methylation patterns remain poorly understood, especially in Spiralia, a large clade of invertebrates comprising nearly half of the animal phyla.

RESULTS

Here, we generate base-resolution methylomes for three species with distinct genomic features and phylogenetic positions in Annelida, a major spiralian phylum. All possible 5mC patterns occur in annelids, from typical invertebrate intermediate levels in a mosaic distribution to hypermethylation and methylation loss. GbM is common to annelids with 5mC, and methylation differences across species are explained by taxon-specific transcriptional dynamics or the presence of intronic TEs. Notably, the link between GbM and transcription decays during development, alongside a gradual and global, age-dependent demethylation in adult stages. Additionally, reducing 5mC levels with cytidine analogs during early development impairs normal embryogenesis and reactivates TEs in the annelid Owenia fusiformis.

CONCLUSIONS

Our study indicates that global epigenetic erosion during development and aging is an ancestral feature of bilateral animals. However, the tight link between transcription and gene body methylation is likely more important in early embryonic stages, and 5mC-mediated TE silencing probably emerged convergently across animal lineages.

Emerging trends in the study of spiralian larvae

Many animals undergo indirect development, where their embryogenesis produces an intermediate life stage, or larva, that is often free-living and later metamorphoses into an adult. As their adult counterparts, larvae can have unique and diverse morphologies and occupy various ecological niches. Given their broad phylogenetic distribution, larvae have been central to hypotheses about animal evolution. However, the evolution of these intermediate forms and the developmental mechanisms diversifying animal life cycles are still debated. This review focuses on Spiralia, a large and diverse clade of bilaterally symmetrical animals with a fascinating array of larval forms, most notably the archetypical trochophore larva. We explore how classic research and modern advances have improved our understanding of spiralian larvae, their development, and evolution. Specifically, we examine three morphological features of spiralian larvae: the anterior neural system, the ciliary bands, and the posterior hyposphere. The combination of molecular and developmental evidence with modern high-throughput techniques, such as comparative genomics, single-cell transcriptomics, and epigenomics, is a promising strategy that will lead to new testable hypotheses about the mechanisms behind the evolution of larvae and life cycles in Spiralia and animals in general. We predict that the increasing number of available genomes for Spiralia and the optimization of genome-wide and single-cell approaches will unlock the study of many emerging spiralian taxa, transforming our views of the evolution of this animal group and their larvae.

A genome resource for the marine annelid Platynereis dumerilii

The marine annelid Platynereis dumerilii is a model organism used in many research areas including evolution and development, neurobiology, ecology and regeneration. Here we present the genomes of P. dumerilii and of the closely related P. massiliensis and P. megalops, to facilitate comparative genomic approaches and help explore Platynereis biology. We used long-read sequencing technology and chromosomal-conformation capture along with extensive transcriptomic resources to obtain and annotate a draft genome assembly of ~1.47 Gbp for P. dumerilii, of which more than half represent repeat elements. We predict around 29,000 protein-coding genes, with relatively large intron sizes, over 38,000 non-coding genes, and 580 miRNA loci. We further explore the high genetic variation (~3% heterozygosity) within the Platynereis species complex. Gene ontology reveals the most variable loci to be associated with pigmentation, development and immunity. The current work sets the stage for further development of Platynereis genomic resources.

Discovery and characterization of a transient chaetal gland during the development of Capitella teleta (Sedentaria: Annelida)

Chaetae are among the most extensively studied structures in polychaetes, serving as a defining morphological trait for annelids. Capitella teleta stands out as one of the few established annelid models for developmental and morphological studies, thus receiving significant scholarly attention. In this study, we unveil a previously unnoticed glandular structure associated with chaetae within the larvae of C. teleta. Our investigations demonstrate the absence of comparable structures in the chaetal follicles of adults and juveniles (older than 1 week), as well as during active chaetogenesis, underscoring the transient nature of these glands. This indicates that larval chaetal follicles transform into a gland that later disappears. Utilizing histology and transmission electron microscopy, we characterized these glands. Our findings underscore the diversity of chaetal ultrastructure in annelids and show that, even in well-studied species, novel morphological details can be found. We emphasize the importance of examining various life-history stages to capture such transient morphological features. This work lays a crucial morphological foundation and deepens our understanding of chaetae and chaetogenesis in C. teleta, paving the way for more accurate interpretations of future experimental studies on chaetogenesis in this species.

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

Revisiting the evolution of Family B1 GPCRs and ligands: Insights from mollusca

Family B1 G protein-coupled receptors (GPCRs) are one of the most well studied neuropeptide receptor families since they play a central role in many biological processes including endocrine, gastrointestinal, cardiovascular and reproduction in animals. The genes for these receptors emerged from a common ancestral gene in bilaterian genomes and evolved via gene/genome duplications and deletions in vertebrate and invertebrate genomes. Their existence and function have mostly been characterized in vertebrates and few studies exist in invertebrate species. Recently, an increased interest in molluscs, means a series of genomes have become available, and since they are less modified than insect and nematode genomes, they are ideal to explore the origin and evolution of neuropeptide gene families. This review provides an overview of Family B1 GPCRs and their peptide ligands and incorporates new data obtained from Mollusca genomes and taking a comparative approach challenges existing models on their origin and evolution.

Acceleration of genome rearrangement in clitellate annelids

Comparisons of multiple metazoan genomes have revealed the existence of ancestral linkage groups (ALGs), genomic scaffolds sharing sets of orthologous genes that have been inherited from ancestral animals for hundreds of millions of years (Simakov et al. 2022; Schultz et al. 2023) These ALGs have persisted across major animal taxa including Cnidaria, Deuterostomia, Ecdysozoa and Spiralia. Notwithstanding this general trend of chromosome-scale conservation, ALGs have been obliterated by extensive genome rearrangements in certain groups, most notably including Clitellata (oligochaetes and leeches), a group of easily overlooked invertebrates that is of tremendous ecological, agricultural and economic importance (Charles 2019; Barrett 2016). To further investigate these rearrangements, we have undertaken a comparison of 12 clitellate genomes (including four newly sequenced species) and 11 outgroup representatives. We show that these rearrangements began at the base of the Clitellata (rather than progressing gradually throughout polychaete annelids), that the inter-chromosomal rearrangements continue in several clitellate lineages and that these events have substantially shaped the evolution of the otherwise highly conserved Hox cluster.

The Satellite DNA PcH-Sat, Isolated and Characterized in the Limpet Patella caerulea (Mollusca, Gastropoda), Suggests the Origin from a Nin-SINE Transposable Element

Satellite DNA (sat-DNA) was previously described as junk and selfish DNA in the cellular economy, without a clear functional role. However, during the last two decades, evidence has been accumulated about the roles of sat-DNA in different cellular functions and its probable involvement in tumorigenesis and adaptation to environmental changes. In molluscs, studies on sat-DNAs have been performed mainly on bivalve species, especially those of economic interest. Conversely, in Gastropoda (which includes about 80% of the currently described molluscs species), studies on sat-DNA have been largely neglected. In this study, we isolated and characterized a sat-DNA, here named PcH-sat, in the limpet Patella caerulea using the restriction enzyme method, particularly HaeIII. Monomeric units of PcH-sat are 179 bp long, AT-rich (58.7%), and with an identity among monomers ranging from 91.6 to 99.8%. Southern blot showed that PcH-sat is conserved in P. depressa and P. ulyssiponensis, while a smeared signal of hybridization was present in the other three investigated limpets (P. ferruginea, P. rustica and P. vulgata). Dot blot showed that PcH-sat represents about 10% of the genome of P. caerulea, 5% of that of P. depressa, and 0.3% of that of P. ulyssiponensis. FISH showed that PcH-sat was mainly localized on pericentromeric regions of chromosome pairs 2 and 4-7 of P. caerulea (2n = 18). A database search showed that PcH-sat contains a large segment (of 118 bp) showing high identity with a homologous trait of the Nin-SINE transposable element (TE) of the patellogastropod Lottia gigantea, supporting the hypothesis that TEs are involved in the rising and tandemization processes of sat-DNAs.

Chromosome-scale Genome Assembly of the Rough Periwinkle Littorina saxatilis

The intertidal gastropod Littorina saxatilis is a model system to study speciation and local adaptation. The repeated occurrence of distinct ecotypes showing different levels of genetic divergence makes L. saxatilis particularly suited to study different stages of the speciation continuum in the same lineage. A major finding is the presence of several large chromosomal inversions associated with the divergence of ecotypes and, specifically, the species offers a system to study the role of inversions in this divergence. The genome of L. saxatilis is 1.35 Gb and composed of 17 chromosomes. The first reference genome of the species was assembled using Illumina data, was highly fragmented (N50 of 44 kb), and was quite incomplete, with a BUSCO completeness of 80.1% on the Metazoan dataset. A linkage map of one full-sibling family enabled the placement of 587 Mbp of the genome into 17 linkage groups corresponding to the haploid number of chromosomes, but the fragmented nature of this reference genome limited the understanding of the interplay between divergent selection and gene flow during ecotype formation. Here, we present a newly generated reference genome that is highly contiguous, with a N50 of 67 Mb and 90.4% of the total assembly length placed in 17 super-scaffolds. It is also highly complete with a BUSCO completeness of 94.1% of the Metazoa dataset. This new reference will allow for investigations into the genomic regions implicated in ecotype formation as well as better characterization of the inversions and their role in speciation.

Chromosome-level genome assembly of Oncomelania hupensis: the intermediate snail host of Schistosoma japonicum

BACKGROUND

Schistosoma japonicum is a parasitic flatworm that causes human schistosomiasis, which is a significant cause of morbidity in China, the Philippines and Indonesia. Oncomelania hupensis (Gastropoda: Pomatiopsidae) is the unique intermediate host of S. japonicum. A complete genome sequence of O. hupensis will enable the fundamental understanding of snail biology as well as its co-evolution with the S. japonicum parasite. Assembling a high-quality reference genome of O. hupehensis will provide data for further research on the snail biology and controlling the spread of S. japonicum.

METHODS

The draft genome was de novo assembly using the long-read sequencing technology (PacBio Sequel II) and corrected with Illumina sequencing data. Then, using Hi-C sequencing data, the genome was assembled at the chromosomal level. CAFE was used to do analysis of contraction and expansion of the gene family and CodeML module in PAML was used for positive selection analysis in protein coding sequences.

RESULTS

A total length of 1.46 Gb high-quality O. hupensis genome with 17 unique full-length chromosomes (2n = 34) of the individual including a contig N50 of 1.35 Mb and a scaffold N50 of 75.08 Mb. Additionally, 95.03% of these contig sequences were anchored in 17 chromosomes. After scanning the assembled genome, a total of 30,604 protein-coding genes were predicted. Among them, 86.67% were functionally annotated. Further phylogenetic analysis revealed that O. hupensis was separated from a common ancestor of Pomacea canaliculata and Bellamya purificata approximately 170 million years ago. Comparing the genome of O. hupensis with its most recent common ancestor, it showed 266 significantly expanded and 58 significantly contracted gene families (P < 0.05). Functional enrichment of the expanded gene families indicated that they were mainly involved with intracellular, DNA-mediated transposition, DNA integration and transposase activity.

CONCLUSIONS

Integrated use of multiple sequencing technologies, we have successfully constructed the genome at the chromosomal-level of O. hupensis. These data will not only provide the compressive genomic information, but also benefit future work on population genetics of this snail as well as evolutional studies between S. japonicum and the snail host.

The development of the adult nervous system in the annelid Owenia fusiformis

BACKGROUND

The evolutionary origins of animal nervous systems remain contentious because we still have a limited understanding of neural development in most major animal clades. Annelids - a species-rich group with centralised nervous systems - have played central roles in hypotheses about the origins of animal nervous systems. However, most studies have focused on adults of deeply nested species in the annelid tree. Recently, Owenia fusiformis has emerged as an informative species to reconstruct ancestral traits in Annelida, given its phylogenetic position within the sister clade to all remaining annelids.

METHODS

Combining immunohistochemistry of the conserved neuropeptides FVamide-lir, RYamide-lir, RGWamide-lir and MIP-lir with gene expression, we comprehensively characterise neural development from larva to adulthood in Owenia fusiformis.

RESULTS

The early larval nervous system comprises a neuropeptide-rich apical organ connected through peripheral nerves to a prototroch ring and the chaetal sac. There are seven sensory neurons in the prototroch. A bilobed brain forms below the apical organ and connects to the ventral nerve cord of the developing juvenile. During metamorphosis, the brain compresses, becoming ring-shaped, and the trunk nervous system develops several longitudinal cords and segmented lateral nerves.

CONCLUSIONS

Our findings reveal the formation and reorganisation of the nervous system during the life cycle of O. fusiformis, an early-branching annelid. Despite its apparent neuroanatomical simplicity, this species has a diverse peptidergic nervous system, exhibiting morphological similarities with other annelids, particularly at the larval stages. Our work supports the importance of neuropeptides in animal nervous systems and highlights how neuropeptides are differentially used throughout development.

De novo transcriptome sequencing and annotation of the Antarctic polychaete Microspio moorei (Spionidae) with its characterization of the heat stress-related proteins (HSP, SOD & CAT)

We present a de novo transcriptome assembly for the non-model Antarctic polychaete worm Microspio moorei (Spionidae) collected during Antarctic field expedition in Fildes Bay, King George Island, Antarctic Peninsula, in 2017. Here, we report the first transcriptome reference array for Microspio spp. The gene sequences of the spionid worm were annotated from a wide range of functions (i.e., biological, and metabolic processes, catalytic processes, and catalytic activity). HSP70, HSP90 SOD and CAT families were compared to reported annelid transcriptomes and proteomes. The phylogenetic analysis using COI, 16S, and 18S markers effectively clusters the species within the family. However, it also casts uncertainty on the monophyletic nature of the Microspio genera, indicating the necessity for additional data and potentially requiring a reevaluation of its grouping. Within these protein families, 3D model software was used to create one representative of their protein structures. Structural predictions were compared with related reported annelids living at different temperatures and a human X-ray reference. We found structural differences (RMSE >1.8) between the human HSP proteins but no significant differences between the polychaete-predicted proteins (RMSE <1.2). These results encourage further research of heat stress-related proteins, the development of genetic markers for climate change-induced temperature stress, and the study of the underlying mechanisms of the heat response. Moreover, these results motivate the extension of these findings to congeneric species.

Coupling DNA barcodes and exon-capture to resolve the phylogeny of Turridae (Gastropoda, Conoidea)

Taxon sampling in most phylogenomic studies is often based on known taxa and/or morphospecies, thus ignoring undescribed diversity and/or cryptic lineages. The family Turridae is a group of venomous snails within the hyperdiverse superfamily Conoidea that includes many undescribed and cryptic species. Therefore 'traditional' taxon sampling could constitute a strong risk of undersampling or oversampling Turridae lineages. To minimize potential biases, we establish a robust sampling strategy, from species delimitation to phylogenomics. More than 3,000 cox-1 "barcode" sequences were used to propose 201 primary species hypotheses, nearly half of them corresponding to species potentially new to science, including several cryptic species. A 110-taxa exon-capture tree, including species representatives of the diversity uncovered with the cox-1 dataset, was build using up to 4,178 loci. Our results show the polyphyly of the genus Gemmula, that is split into up to 10 separate lineages, of which half would not have been detected if the sampling strategy was based only on described species. Our results strongly suggest that the use of blind, exploratory and intensive barcode sampling is necessary to avoid sampling biases in phylogenomic studies.

A Chromosome-Level Genome Assembly of the Non-Hematophagous Leech Whitmania pigra (Whitman 1884): Identification and Expression Analysis of Antithrombotic Genes

Despite being a non-hematophagous leech, Whitmania pigra is widely used in traditional Chinese medicine for the treatment of antithrombotic diseases. In this study, we provide a high quality genome of W. pigra and based on which, we performed a systematic identification of the potential antithrombotic genes and their corresponding proteins. We identified twenty antithrombotic gene families including thirteen coagulation inhibitors, three platelet aggregation inhibitors, three fibrinolysis enhancers, and one tissue penetration enhancer. Unexpectedly, a total of 79 antithrombotic genes were identified, more than a typical blood-feeding Hirudinaria manillensis, which had only 72 antithrombotic genes. In addition, combining with the RNA-seq data of W. pigra and H. manillensis, we calculated the expression levels of antithrombotic genes of the two species. Five and four gene families had significantly higher and lower expression levels in W. pigra than in H. manillensis, respectively. These results showed that the number and expression level of antithrombotic genes of a non-hematophagous leech are not always less than those of a hematophagous leech. Our study provides the most comprehensive collection of antithrombotic biomacromolecules from a non-hematophagous leech to date and will significantly enhance the investigation and utilization of leech derivatives in thrombosis therapy research and pharmaceutical applications.

A chromosome-level genome for the nudibranch gastropod Berghia stephanieae helps parse clade-specific gene expression in novel and conserved phenotypes

BACKGROUND

How novel phenotypes originate from conserved genes, processes, and tissues remains a major question in biology. Research that sets out to answer this question often focuses on the conserved genes and processes involved, an approach that explicitly excludes the impact of genetic elements that may be classified as clade-specific, even though many of these genes are known to be important for many novel, or clade-restricted, phenotypes. This is especially true for understudied phyla such as mollusks, where limited genomic and functional biology resources for members of this phylum have long hindered assessments of genetic homology and function. To address this gap, we constructed a chromosome-level genome for the gastropod Berghia stephanieae (Valdés, 2005) to investigate the expression of clade-specific genes across both novel and conserved tissue types in this species.

RESULTS

The final assembled and filtered Berghia genome is comparable to other high-quality mollusk genomes in terms of size (1.05 Gb) and number of predicted genes (24,960 genes) and is highly contiguous. The proportion of upregulated, clade-specific genes varied across tissues, but with no clear trend between the proportion of clade-specific genes and the novelty of the tissue. However, more complex tissue like the brain had the highest total number of upregulated, clade-specific genes, though the ratio of upregulated clade-specific genes to the total number of upregulated genes was low.

CONCLUSIONS

Our results, when combined with previous research on the impact of novel genes on phenotypic evolution, highlight the fact that the complexity of the novel tissue or behavior, the type of novelty, and the developmental timing of evolutionary modifications will all influence how novel and conserved genes interact to generate diversity.

Chromosome-level genome assemblies of two littorinid marine snails indicate genetic basis of intertidal adaptation and ancient karyotype evolved from bilaterian ancestors

Living in the intertidal environment, littorinid snails are excellent models for understanding genetic mechanisms underlying adaptation to harsh fluctuating environments. Furthermore, the karyotypes of littorinid snails, with the same chromosome number as the presumed bilaterian ancestor, make them valuable for investigating karyotype evolution from the bilaterian ancestor to mollusks. Here, we generated high-quality, chromosome-scale genome assemblies for 2 littorinid marine snails, Littorina brevicula (927.94 Mb) and Littoraria sinensis (882.51 Mb), with contig N50 of 3.43 Mb and 2.31 Mb, respectively. Comparative genomic analyses identified 92 expanded gene families and 85 positively selected genes as potential candidates possibly associated with intertidal adaptation in the littorinid lineage, which were functionally enriched in stimulus responses, innate immunity, and apoptosis process regulation and might be involved in cellular homeostasis maintenance in stressful intertidal environments. Genome macrosynteny analyses indicated that 4 fissions and 4 fusions led to the evolution from the 17 presumed bilaterian ancestral chromosomes to the 17 littorinid chromosomes, implying that the littorinid snails have a highly conserved karyotype with the bilaterian ancestor. Based on the most parsimonious reconstruction of the common ancestral karyotype of scallops and littorinid snails, 3 chromosomal fissions and 1 chromosomal fusion from the bilaterian ancient linkage groups were shared by the bivalve scallop and gastropoda littorinid snails, indicating that the chromosome-scale ancient gene linkages were generally preserved in the mollusk genomes for over 500 million years. The highly conserved karyotype makes the littorinid snail genomes valuable resources for understanding early bilaterian evolution and biology.

A second unveiling: Haplotig masking of the eastern oyster genome improves population-level inference

Genome assembly can be challenging for species that are characterized by high amounts of polymorphism, heterozygosity, and large effective population sizes. High levels of heterozygosity can result in genome mis-assemblies and a larger than expected genome size due to the haplotig versions of a single locus being assembled as separate loci. Here, we describe the first chromosome-level genome for the eastern oyster, Crassostrea virginica. Publicly released and annotated in 2017, the assembly has a scaffold N50 of 54 mb and is over 97.3% complete based on BUSCO analysis. The genome assembly for the eastern oyster is a critical resource for foundational research into molluscan adaptation to a changing environment and for selective breeding for the aquaculture industry. Subsequent resequencing data suggested the presence of haplotigs in the original assembly, and we developed a post hoc method to break up chimeric contigs and mask haplotigs in published heterozygous genomes and evaluated improvements to the accuracy of downstream analysis. Masking haplotigs had a large impact on SNP discovery and estimates of nucleotide diversity and had more subtle and nuanced effects on estimates of heterozygosity, population structure analysis, and outlier detection. We show that haplotig masking can be a powerful tool for improving genomic inference, and we present an open, reproducible resource for the masking of haplotigs in any published genome.

Hox, homology, and parsimony: An organismal perspective

Hox genes are important regulators in animal development. They often show a mosaic of conserved (e.g., longitudinal axis patterning) and lineage-specific novel functions (e.g., development of skeletal, sensory, or locomotory systems). Despite extensive research over the past decades, it remains controversial at which node in the animal tree of life the Hox cluster evolved. Its presence already in the last common metazoan ancestor has been proposed, although the genomes of both putative earliest extant metazoan offshoots, the ctenophores and the poriferans, are devoid of Hox sequences. The lack of Hox genes in the supposedly "simple"-built poriferans and their low number in cnidarians and the basally branching bilaterians, the xenacoelomorphs, seems to support the classical notion that the number of Hox genes is correlated with the degree of animal complexity. However, the 4-fold increase of the Hox cluster in xiphosurans, a basally branching chelicerate clade, as well as the situation in some teleost fishes that show a multitude of Hox genes compared to, e.g., human, demonstrates, that there is no per se direct correlation between organismal complexity and Hox number. Traditional approaches have tried to base homology on the morphological level on shared expression profiles of individual genes, but recent data have shown that, in particular with respect to Hox and other regulatory genes, complex gene-gene interactions rather than expression signatures of individual genes alone are responsible for shaping morphological traits during ontogeny. Accordingly, for sound homology assessments and reconstructions of character evolution on organ system level, additional independent datasets (e.g., morphological, developmental) need to be included in any such analyses. If supported by solid data, proposed structural homology should be regarded as valid and not be rejected solely on the grounds of non-parsimonious distribution of the character over a given phylogenetic topology.

MarFERReT, an open-source, version-controlled reference library of marine microbial eukaryote functional genes

Metatranscriptomics generates large volumes of sequence data about transcribed genes in natural environments. Taxonomic annotation of these datasets depends on availability of curated reference sequences. For marine microbial eukaryotes, current reference libraries are limited by gaps in sequenced organism diversity and barriers to updating libraries with new sequence data, resulting in taxonomic annotation of about half of eukaryotic environmental transcripts. Here, we introduce Marine Functional EukaRyotic Reference Taxa (MarFERReT), a marine microbial eukaryotic sequence library designed for use with taxonomic annotation of eukaryotic metatranscriptomes. We gathered 902 publicly accessible marine eukaryote genomes and transcriptomes and assessed their sequence quality and cross-contamination issues, selecting 800 validated entries for inclusion in MarFERReT. Version 1.1 of MarFERReT contains reference sequences from 800 marine eukaryotic genomes and transcriptomes, covering 453 species- and strain-level taxa, totaling nearly 28 million protein sequences with associated NCBI and PR2 Taxonomy identifiers and Pfam functional annotations. The MarFERReT project repository hosts containerized build scripts, documentation on installation and use case examples, and information on new versions of MarFERReT.

Topological structures and syntenic conservation in sea anemone genomes

There is currently little information about the evolution of gene clusters, genome architectures and karyotypes in early branching animals. Slowly evolving anthozoan cnidarians can be particularly informative about the evolution of these genome features. Here we report chromosome-level genome assemblies of two related anthozoans, the sea anemones Nematostella vectensis and Scolanthus callimorphus. We find a robust set of 15 chromosomes with a clear one-to-one correspondence between the two species. Both genomes show chromosomal conservation, allowing us to reconstruct ancestral cnidarian and metazoan chromosomal blocks, consisting of at least 19 and 16 ancestral linkage groups, respectively. We show that, in contrast to Bilateria, the Hox and NK clusters of investigated cnidarians are largely disintegrated, despite the presence of staggered hox/gbx expression in Nematostella. This loss of microsynteny conservation may be facilitated by shorter distances between cis-regulatory sequences and their cognate transcriptional start sites. We find no clear evidence for topologically associated domains, suggesting fundamental differences in long-range gene regulation compared to vertebrates. These data suggest that large sets of ancestral metazoan genes have been retained in ancestral linkage groups of some extant lineages; yet, higher order gene regulation with associated 3D architecture may have evolved only after the cnidarian-bilaterian split.

Comparative Analysis of Bivalve and Sea Urchin Genetics and Development: Investigating the Dichotomy in Bilateria

This comprehensive review presents a comparative analysis of early embryogenesis in Protostomia and Deuterostomia, the first of which exhibit a mosaic pattern of development, where cells are fated deterministically, while Deuterostomia display a regulatory pattern of development, where the fate of cells is indeterminate. Despite these fundamental differences, there are common transcriptional mechanisms that underline their evolutionary linkages, particularly in the field of functional genomics. By elucidating both conserved and unique regulatory strategies, this review provides essential insights into the comparative embryology and developmental dynamics of these groups. The objective of this review is to clarify the shared and distinctive characteristics of transcriptional regulatory mechanisms. This will contribute to the extensive areas of functional genomics, evolutionary biology and developmental biology, and possibly lay the foundation for future research and discussion on this seminal topic.

Genomic and transcriptomic survey of bryozoan Hox and ParaHox genes with emphasis on phylactolaemate bryozoans

BACKGROUND

Bryozoans are mostly sessile aquatic colonial invertebrates belonging to the clade Lophotrochozoa, which unites many protostome bilaterian phyla such as molluscs, annelids and brachiopods. While Hox and ParaHox genes have been extensively studied in various lophotrochozoan lineages, investigations on Hox and ParaHox gene complements in bryozoans are scarce.

RESULTS

Herein, we present the most comprehensive survey of Hox and ParaHox gene complements in bryozoans using four genomes and 35 transcriptomes representing all bryozoan clades: Cheilostomata, Ctenostomata, Cyclostomata and Phylactolaemata. Using similarity searches, phylogenetic analyses and detailed manual curation, we have identified five Hox genes in bryozoans (pb, Dfd, Lox5, Lox4 and Post2) and one ParaHox gene (Cdx). Interestingly, we observed lineage-specific duplication of certain Hox and ParaHox genes (Dfd, Lox5 and Cdx) in some bryozoan lineages.

CONCLUSIONS

The bryozoan Hox cluster does not retain the ancestral lophotrochozoan condition but appears relatively simple (includes only five genes) and broken into two genomic regions, characterized by the loss and duplication of serval genes. Importantly, bryozoans share the lack of two Hox genes (Post1 and Scr) with their proposed sister-taxon, Phoronida, which suggests that those genes were missing in the most common ancestor of bryozoans and phoronids.

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.

Emerging questions on the mechanisms and dynamics of 3D genome evolution in spiralians

Information on how 3D genome topology emerged in animal evolution, how stable it is during development, its role in the evolution of phenotypic novelties and how exactly it affects gene expression is highly debated. So far, data to address these questions are lacking with the exception of a few key model species. Several gene regulatory mechanisms have been proposed, including scenarios where genome topology has little to no impact on gene expression, and vice versa. The ancient and diverse clade of spiralians may provide a crucial testing ground for such mechanisms. Sprialians have followed distinct evolutionary trajectories, with some clades experiencing genome expansions and/or large-scale genome rearrangements, and others undergoing genome contraction, substantially impacting their size and organisation. These changes have been associated with many phenotypic innovations in this clade. In this review, we describe how emerging genome topology data, along with functional tools, allow for testing these scenarios and discuss their predicted outcomes.

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.

Chromosome-level genome assembly of the deep-sea snail Phymorhynchus buccinoides provides insights into the adaptation to the cold seep habitat

BACKGROUND

The deep-sea snail Phymorhynchus buccinoides belongs to the genus Phymorhynchus (Neogastropoda: Raphitomidae), and it is a dominant specie in the cold seep habitat. As the environment of the cold seep is characterized by darkness, hypoxia and high concentrations of toxic substances such as hydrogen sulfide (H2S), exploration of the diverse fauna living around cold seeps will help to uncover the adaptive mechanisms to this unique habitat. In the present study, a chromosome-level genome of P. buccinoides was constructed and a series of genomic and transcriptomic analyses were conducted to explore its molecular adaptation mechanisms to the cold seep environments.

RESULTS

The assembled genome size of the P. buccinoides was approximately 2.1 Gb, which is larger than most of the reported snail genomes, possibly due to the high proportion of repetitive elements. About 92.0% of the assembled base pairs of contigs were anchored to 34 pseudo-chromosomes with a scaffold N50 size of 60.0 Mb. Compared with relative specie in the shallow water, the glutamate regulative and related genes were expanded in P. buccinoides, which contributes to the acclimation to hypoxia and coldness. Besides, the relatively high mRNA expression levels of the olfactory/chemosensory genes in osphradium indicate that P. buccinoides might have evolved a highly developed and sensitive olfactory organ for its orientation and predation. Moreover, the genome and transcriptome analyses demonstrate that P. buccinoides has evolved a sulfite-tolerance mechanism by performing H2S detoxification. Many genes involved in H2S detoxification were highly expressed in ctenidium and hepatopancreas, suggesting that these tissues might be critical for H2S detoxification and sulfite tolerance.

CONCLUSIONS

In summary, our report of this chromosome-level deep-sea snail genome provides a comprehensive genomic basis for the understanding of the adaptation strategy of P. buccinoides to the extreme environment at the deep-sea cold seeps.

Identification and functional analysis of myeloid differentiation factor 88 (MyD88) in early development of Haliotis diversicolor

Myeloid differentiation factor 88 (MyD88) is a universal adaptor protein and plays an important role in the signal transduction of Toll like receptors (TLR) family. In this study, the MyD88 gene from the Haliotis diversicolor (hdMyD88) was identified. The full-length cDNA of hdMyD88 has a 1927 base pairs (bp), with an open reading frame of 1314 bp encoding 437 amino acids including a death domain (DD) at the N-terminus and TIR domain at the C-terminus which are typical features of MyD88 family proteins. Three conserved boxes are also found in the hdMyD88, which are similar to MyD88 in vertebrates. The expression levels of hdMyD88 mRNA at different early embryonic developmental stages of abalone were measured by qPCR revealed that their constitutive expression at all developmental stages analyzed with the considerably highest values at 8 cell stage and the lowest level at the trochosphere stage. Additionally, the mRNA expression of hdMyD88 decreased significantly (P < 0.05) after MyD88-dsRNA soak in the stage of trochosphere and veliger than EGFP-dsRNA group and blank control group. Whole embryo in situ hybridization showed that the positive signals of hdMyD88 were in visceral mass of trochophore larvae and veliger larvae. These results indicate hdMyD88 may could respond to pathogenic infection and may play an important role in early innate immunity in the process of abalone larval development.

Collagen IV of basement membranes: I. Origin and diversification of COL4 genes enabling animal evolution

Collagen IV is a primordial component of basement membranes, a specialized form of extracellular matrix that enabled multi-cellular epithelial tissues. In mammals, collagen IV assembles from a family of six α-chains (α1 to α6), encoded by six genes (COL4A1 to COL4A6), into three distinct scaffolds: the α121, the α345 and a mixed scaffold containing both α121 and α565. The six mammalian COL4A genes occur in pairs that occur in a head-to-head arrangement on three distinct chromosomes. In Alport syndrome, variants in the COL4A3, 4 or 5 genes cause either loss or defective assembly of the collagen IV α345 scaffold which results in a dysfunctional glomerular basement membrane, proteinuria and progression to renal failure in millions of people worldwide. Here, we determine the evolutionary emergence and diversification of the COL4A genes using comparative genomics and biochemical analyses. Using syntenic relationships to genes closely linked to the COL4A genes, we determine that the COL4A3 and COL4A4 gene pair appeared in cyclostomes (hagfish and lampreys) while the COL4A5 and COL4A6 gene pair emerged in gnathostomes, jawed vertebrates. The more basal chordate species, lancelets and tunicates, do not have discrete kidneys and have a single COL4A gene pair, though often with single isolated COL4 genes similar to those found in C elegans . Remarkably, while the six COL4A genes are conserved in vertebrates, amphibians have lost the COL4A3 and COL4A4 genes. Our findings of the evolutionary emergence of these genes, together with the amphibian double-knockout, opens an experimental window to gain insights into functionality of the Col IV α345 scaffold.

Scaphopoda is the sister taxon to Bivalvia: Evidence of ancient incomplete lineage sorting

The almost simultaneous emergence of major animal phyla during the early Cambrian shaped modern animal biodiversity. Reconstructing evolutionary relationships among such closely spaced branches in the animal tree of life has proven to be a major challenge, hindering understanding of early animal evolution and the fossil record. This is particularly true in the species-rich and highly varied Mollusca where dramatic inconsistency among paleontological, morphological, and molecular evidence has led to a long-standing debate about the group's phylogeny and the nature of dozens of enigmatic fossil taxa. A critical step needed to overcome this issue is to supplement available genomic data, which is plentiful for well-studied lineages, with genomes from rare but key lineages, such as Scaphopoda. Here, by presenting chromosome-level genomes from both extant scaphopod orders and leveraging complete genomes spanning Mollusca, we provide strong support for Scaphopoda as the sister taxon of Bivalvia, revitalizing the morphology-based Diasoma hypothesis originally proposed 50 years ago. Our molecular clock analysis confidently dates the split between Bivalvia and Scaphopoda at ~520 Ma, prompting a reinterpretation of controversial laterally compressed Early Cambrian fossils, including Anabarella, Watsonella, and Mellopegma, as stem diasomes. Moreover, we show that incongruence in the phylogenetic placement of Scaphopoda in previous phylogenomic studies was due to ancient incomplete lineage sorting (ILS) that occurred during the rapid radiation of Conchifera. Our findings highlight the need to consider ILS as a potential source of error in deep phylogeny reconstruction, especially in the context of the unique nature of the Cambrian Explosion.

Functionally diverse front-end desaturases are widespread in the phylum Annelida

Aquatic single-cell organisms have long been believed to be unique primary producers of omega-3 long-chain (≥C20) polyunsaturated fatty acids (ω3 LC-PUFA). Multiple invertebrates including annelids have been discovered to possess methyl-end desaturases enabling key steps in the de novo synthesis of ω3 LC-PUFA, and thus potentially contributing to their production in the ocean. Along methyl-end desaturases, the repertoire and function of further LC-PUFA biosynthesising enzymes is largely missing in Annelida. In this study we examined the front-end desaturase gene repertoire across the phylum Annelida, from Polychaeta and Clitellata, major classes of annelids comprising most annelid diversity. We further characterised the functions of the encoded enzymes in selected representative species by using a heterologous expression system based in yeast, demonstrating that functions of Annelida front-end desaturases have highly diversified during their expansion in both terrestrial and aquatic ecosystems. We concluded that annelids possess at least two front-end desaturases with Δ5 and Δ6Δ8 desaturase regioselectivities, enabling all the desaturation reactions required to convert the C18 precursors into the physiologically relevant LC-PUFA such as eicosapentaenoic and arachidonic acids, but not docosahexaenoic acid. Such a gene complement is conserved across the different taxonomic groups within Annelida.

Functional Characterization of Gonadotropin-Releasing Hormone and Corazonin Signaling Systems in Pacific Abalone: Toward Reclassification of Invertebrate Neuropeptides

INTRODUCTION

The proposed evolutionary origins and corresponding nomenclature of bilaterian gonadotropin-releasing hormone (GnRH)-related neuropeptides have changed tremendously with the aid of receptor deorphanization. However, the reclassification of the GnRH and corazonin (CRZ) signaling systems in Lophotrochozoa remains unclear.

METHODS

We characterized GnRH and CRZ receptors in the mollusk Pacific abalone, Haliotis discus hannai (Hdh), by phylogenetic and gene expression analyses, bioluminescence-based reporter, Western blotting, substitution of peptide amino acids, in vivo neuropeptide injection, and RNA interference assays.

RESULTS

Two Hdh CRZ-like receptors (Hdh-CRZR-A and Hdh-CRZR-B) and three Hdh GnRH-like receptors (Hdh-GnRHR1-A, Hdh-GnRHR1-B, and Hdh-GnRHR2) were identified. In phylogenetic analysis, Hdh-CRZR-A and -B grouped within the CRZ-type receptors, whereas Hdh-GnRHR1-A/-B and Hdh-GnRHR2 clustered within the GnRH/adipokinetic hormone (AKH)/CRZ-related peptide-type receptors. Hdh-CRZR-A/-B and Hdh-GnRHR1-A were activated by Hdh-CRZ (pQNYHFSNGWHA-NH2) and Hdh-GnRH (pQISFSPNWGT-NH2), respectively. Hdh-CRZR-A/-B dually coupled with the Gαq and Gαs signaling pathways, whereas Hdh-GnRHR1-A was linked only with Gαq signaling. Analysis of substituted peptides, [I2S3]Hdh-CRZ and [N2Y3H4]Hdh-GnRH, and in silico docking models revealed that the N-terminal amino acids of the peptides are critical for the selectivity of Hdh-CRZR and Hdh-GnRHR. Two precursor transcripts for Hdh-CRZ and Hdh-GnRH peptides and their receptors were mainly expressed in the neural ganglia, and their levels increased in starved abalones. Injection of Hdh-CRZ peptide into abalones decreased food consumption, whereas Hdh-CRZR knockdown increased food consumption. Moreover, Hdh-CRZ induced germinal vesicle breakdown in mature oocytes.

CONCLUSION

Characterization of Hdh-CRZRs and Hdh-GnRHRs and their cognate peptides provides new insight into the evolutionary route of GnRH-related signaling systems in bilaterians.

Updated knowledge and a proposed nomenclature for nuclear receptors with two DNA binding domains (2DBD-NRs)

Nuclear receptors (NRs) are important transcriptional modulators in metazoans. Typical NRs possess a conserved DNA binding domain (DBD) and a ligand binding domain (LBD). Since we discovered a type of novel NRs each of them has two DBDs and single LBD (2DBD-NRs) more than decade ago, there has been very few studies about 2DBD-NRs. Recently, 2DBD-NRs have been only reported in Platyhelminths and Mollusca and are thought to be specific NRs to lophotrochozoan. In this study, we searched different databases and identified 2DBD-NRs in different animals from both protostomes and deuterostomes. Phylogenetic analysis shows that at least two ancient 2DBD-NR genes were present in the urbilaterian, a common ancestor of protostomes and deuterostomes. 2DBD-NRs underwent gene duplication and loss after the split of different animal phyla, most of them in a certain animal phylum are paralogues, rather than orthologues, like in other animal phyla. Amino acid sequence analysis shows that the conserved motifs in typical NRs are also present in 2DBD-NRs and they are gene specific. From our phylogenetic analysis of 2DBD-NRs and following the rule of Nomenclature System for the Nuclear Receptors, a nomenclature for 2DBD-NRs is proposed.

Chromosome-level genome assembly of the caenogastropod snail Rapana venosa

The carnivorous gastropod Rapana venosa (Valenciennes, 1846) is one of the most notorious ecological invaders worldwide. Here, we present the first high-quality chromosome-scale reference R. venosa genome obtained via PacBio sequencing, Illumina paired-end sequencing, and high-throughput chromosome conformation capture scaffolding. The assembled genome has a size of 2.30 Gb, with a scaffold N50 length of 64.63 Mb, and is anchored to 35 chromosomes. It contains 29,649 protein-coding genes, 77.22% of which were functionally annotated. Given its high heterozygosity (1.41%) and large proportion of repeat sequences (57.72%), it is one of the most complex genome assemblies. This chromosome-level genome assembly of R. venosa is an important resource for understanding molluscan evolutionary adaption and provides a genetic basis for its biological invasion control.

Genome assembly of the Korean intertidal mud-creeper Batillaria attramentaria

Batillaridae is a common gastropod family that occurs abundantly in the shallow coastal zone of the intertidal mudflats of the northwest Pacific Ocean, Australasia, and North America. In this family, Batillaria attramentaria is known for its biological invasion and colonization in estuarine and intertidal zones. It can endure and adapt the harsh intertidal conditions such as frequent temperature alteration, salinity, and air exposure. Therefore, we sequenced and assembled this Korean batillariid genome to get insight into its intertidal adaptive features. Approximately 53 Gb of DNA sequences were generated, and 863 scaffolds were assembled into a draft genome of 0.715 Gb with 97.1% BUSCO completeness value. A total of 40,596 genes were predicted. We estimated that B. attramentaria and Conus consors diverged about 230 million years ago (MYA) based on the phylogenetic analysis of closely related gastropod species. This genome study sets the footstep for genomics studies among native and introduced Batillaria populations and the Batillaridae family members.

Spatial Colinear but Broken Temporal Expression of Duplicated ParaHox Genes in Asexually Reproducing Annelids, Nais communis and Pristina longiseta

ParaHox genes are key developmental regulators involved in the patterning of the digestive tract along the anteroposterior axis and the development of the nervous system. Most studies have focused on the function of these genes in embryogenesis, while their expression patterns in postembryonic development often remain unknown. In this study, we identified for the first time all ParaHox orthologs in two naidid oligochaetes, N. communis and P. longiseta, and described their expression patterns during normal growth and fission in these animals. We showed that Gsx and Cdx are presented by two paralogs, while Xlox is a single copy gene in both species. Using whole-mount in situ hybridization, we also found that orthologs, except for the Xlox gene, have similar activity patterns with minor differences in details, while the expression patterns of paralogs can differ significantly. However, all these genes are involved in axial patterning and/or in tissue remodeling during growth and asexual reproduction in naidids. Moreover, during paratomic fission, these genes are expressed with spatial colinearity but temporal colinearity is broken. The results of this study may be evidence of the functional diversification of duplicated genes and suggest involvement of the ParaHox genes in whole-body patterning during growth and asexual reproduction in annelids.

A High-Quality Chromosome-Level Genome Assembly of a Snail Cipangopaludina cathayensis (Gastropoda: Viviparidae)

Cipangopaludina cathayensis (Gastropoda: Prosobranchia; Mesogastropoda; Viviparidae) is widely distributed in the freshwater habitats of China. It is an economically important snail with high edible and medicinal value. However, the genomic resources and the reference genome of this snail are lacking. In this study, we assembled the first chromosome-level genome of C. cathayensis. The preliminary assembly genome was 1.48 Gb in size, with a contig N50 size of 93.49 Mb. The assembled sequences were anchored to nine pseudochromosomes using Hi-C data. The final genome after Hi-C correction was 1.48 Gb, with a contig N50 of 98.49 Mb and scaffold N50 of 195.21 Mb. The anchored rate of the chromosome was 99.99%. A total of 22,702 protein-coding genes were predicted. Phylogenetic analyses indicated that C. cathayensis diverged with Bellamya purificata approximately 158.10 million years ago. There were 268 expanded and 505 contracted gene families in C. cathayensis when compared with its most recent common ancestor. Five putative genes under positive selection in C. cathayensis were identified (false discovery rate <0.05). These genome data provide a valuable resource for evolutionary studies of the family Viviparidae, and for the genetic improvement of C. cathayensis.

Pyroptotic gasdermin exists in Mollusca and is vital to eliminating bacterial infection

Gasdermin (GSDM) is a family of proteins that execute pyroptosis in vertebrate. In invertebrate, pyroptotic GSDM was documented only in coral. Recent studies identified abundant GSDM structural homologs in Mollusca, but their functions are unclear. Herein, we report a functional GSDM from Pacific abalone Haliotis discus (HdGSDME). HdGSDME is specifically activated by abalone caspase 3 (HdCASP3) cleavage at two distinct sites, generating two active isoforms with pyroptotic and cytotoxic activities. HdGSDME possesses evolutionarily conserved residues that proved to be essential to the N-terminal pore-formation and C-terminal auto-inhibition capacities. Bacterial challenge activates the HdCASP3-HdGSDME pathway and induces pyroptosis and extracellular traps in abalone. Blockage of the HdCASP3-HdGSDME axis promotes bacterial invasion and host mortality. Collectively, this study reveals the existence of functionally conserved and yet distinct-featured GSDM in Mollusca and provides insights into the function and evolution of invertebrate GSDM.

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.

Unexpected Distribution of Chitin and Chitin Synthase across Soft-Bodied Cnidarians

Cnidarians are commonly recognized as sea jellies, corals, or complex colonies such as the Portuguese man-of-war. While some cnidarians possess rigid internal calcareous skeletons (e.g., corals), many are soft-bodied. Intriguingly, genes coding for the chitin-biosynthetic enzyme, chitin synthase (CHS), were recently identified in the model anemone Nematostella vectensis, a species lacking hard structures. Here we report the prevalence and diversity of CHS across Cnidaria and show that cnidarian chitin synthase genes display diverse protein domain organizations. We found that CHS is expressed in cnidarian species and/or developmental stages with no reported chitinous or rigid morphological structures. Chitin affinity histochemistry indicates that chitin is present in soft tissues of some scyphozoan and hydrozoan medusae. To further elucidate the biology of chitin in cnidarian soft tissues, we focused on CHS expression in N. vectensis. Spatial expression data show that three CHS orthologs are differentially expressed in Nematostella embryos and larvae during development, suggesting that chitin has an integral role in the biology of this species. Understanding how a non-bilaterian lineage such as Cnidaria employs chitin may provide new insight into hitherto unknown functions of polysaccharides in animals, as well as their role in the evolution of biological novelty.

Division of developmental phases of freshwater leech Whitmania pigra and key genes related to neurogenesis revealed by whole genome and transcriptome analysis

The freshwater leech Whitmania pigra (W. pigra) Whitman (Annelida phylum) is a model organism for neurodevelopmental studies. However, molecular biology research on its embryonic development is still scarce. Here, we described a series of developmental stages of the W. pigra embryos and defined five broad stages of embryogenesis: cleavage stages, blastocyst stage, gastrula stage, organogenesis and refinement, juvenile. We obtained a total of 239.64 Gb transcriptome data of eight representative developmental phases of embryos (from blastocyst stage to maturity), which was then assembled into 21,482 unigenes according to our reference genome sequenced by single-molecule real-time (SMRT) long-read sequencing. We found 3114 genes differentially expressed during the eight phases with phase-specific expression pattern. Using a comprehensive transcriptome dataset, we demonstrated that 57, 49 and 77 DEGs were respectively related to morphogenesis, signal pathways and neurogenesis. 49 DEGs related to signal pathways included 30 wnt genes, 14 notch genes, and 5 hedgehog genes. In particular, we found a cluster consisting of 7 genes related to signal pathways as well as synapses, which were essential for regulating embryonic development. Eight genes cooperatively participated in regulating neurogenesis. Our results reveal the whole picture of W. pigra development mechanism from the perspective of transcriptome and provide new clues for organogenesis and neurodevelopmental studies of Annelida species.

Ancestral protein reconstruction reveals evolutionary events governing variation in Dicer helicase function

Antiviral defense in ecdysozoan invertebrates requires Dicer with a helicase domain capable of ATP hydrolysis. But despite well-conserved ATPase motifs, human Dicer is incapable of ATP hydrolysis, consistent with a muted role in antiviral defense. To investigate this enigma, we used ancestral protein reconstruction to resurrect Dicer's helicase in animals and trace the evolutionary trajectory of ATP hydrolysis. Biochemical assays indicated ancient Dicer possessed ATPase function, that like extant invertebrate Dicers, is stimulated by dsRNA. Analyses revealed that dsRNA stimulates ATPase activity by increasing ATP affinity, reflected in Michaelis constants. Deuterostome Dicer-1 ancestor, while exhibiting lower dsRNA affinity, retained some ATPase activity; importantly, ATPase activity was undetectable in the vertebrate Dicer-1 ancestor, which had even lower dsRNA affinity. Reverting residues in the ATP hydrolysis pocket was insufficient to rescue hydrolysis, but additional substitutions distant from the pocket rescued vertebrate Dicer-1's ATPase function. Our work suggests Dicer lost ATPase function in the vertebrate ancestor due to loss of ATP affinity, involving motifs distant from the active site, important for coupling dsRNA binding to the active conformation. By competing with Dicer for viral dsRNA, RIG-I-like receptors important for interferon signaling may have allowed or actively caused loss of ATPase function.

Evolution of homology: From archetype towards a holistic concept of cell type

The concept of homology lies in the heart of comparative biological science. The distinction between homology as structure and analogy as function has shaped the evolutionary paradigm for a century and formed the axis of comparative anatomy and embryology, which accept the identity of structure as a ground measure of relatedness. The advent of single-cell genomics overturned the classical view of cell homology by establishing a backbone regulatory identity of cell types, the basic biological units bridging the molecular and phenotypic dimensions, to reveal that the cell is the most flexible unit of living matter and that many approaches of classical biology need to be revised to understand evolution and diversity at the cellular level. The emerging theory of cell types explicitly decouples cell identity from phenotype, essentially allowing for the divergence of evolutionarily related morphotypes beyond recognition, as well as it decouples ontogenetic cell lineage from cell-type phylogeny, whereby explicating that cell types can share common descent regardless of their structure, function or developmental origin. The article succinctly summarizes current progress and opinion in this field and formulates a more generalistic view of biological cell types as avatars, transient or terminal cell states deployed in a continuum of states by the developmental programme of one and the same omnipotent cell, capable of changing or combining identities with distinct evolutionary histories or inventing ad hoc identities that never existed in evolution or development. It highlights how the new logic grounded in the regulatory nature of cell identity transforms the concepts of cell homology and phenotypic stability, suggesting that cellular evolution is inherently and massively network-like, with one-to-one homologies being rather uncommon and restricted to shallower levels of the animal tree of life.

Robust 3D modeling reveals spatiosyntenic properties of animal genomes

Animal genomes are organized into chromosomes that are remarkably conserved in their gene content, forming distinct evolutionary units (synteny). Using versatile chromosomal modeling, we infer three-dimensional topology of genomes from representative clades spanning the earliest animal diversification. We apply a partitioning approach using interaction spheres to compensate for varying quality of topological data. Using comparative genomics approaches, we test whether syntenic signal at gene pair, local, and whole chromosomal scale is reflected in the reconstructed spatial organization. We identify evolutionarily conserved three-dimensional networks at all syntenic scales revealing novel evolutionarily maintained interactors associated with known conserved local gene linkages (such as hox). We thus present evidence for evolutionary constraints that are associated with three-, rather than just two-, dimensional animal genome organization, which we term spatiosynteny. As more accurate topological data become available, together with validation approaches, spatiosynteny may become relevant in understanding the functionality behind the observed conservation of animal chromosomes.

Hologenome analysis reveals independent evolution to chemosymbiosis by deep-sea bivalves

BACKGROUND

Bivalves have independently evolved a variety of symbiotic relationships with chemosynthetic bacteria. These relationships range from endo- to extracellular interactions, making them ideal for studies on symbiosis-related evolution. It is still unclear whether there are universal patterns to symbiosis across bivalves. Here, we investigate the hologenome of an extracellular symbiotic thyasirid clam that represents the early stages of symbiosis evolution.

RESULTS

We present a hologenome of Conchocele bisecta (Bivalvia: Thyasiridae) collected from deep-sea hydrothermal vents with extracellular symbionts, along with related ultrastructural evidence and expression data. Based on ultrastructural and sequencing evidence, only one dominant Thioglobaceae bacteria was densely aggregated in the large bacterial chambers of C. bisecta, and the bacterial genome shows nutritional complementarity and immune interactions with the host. Overall, gene family expansions may contribute to the symbiosis-related phenotypic variations in different bivalves. For instance, convergent expansions of gaseous substrate transport families in the endosymbiotic bivalves are absent in C. bisecta. Compared to endosymbiotic relatives, the thyasirid genome exhibits large-scale expansion in phagocytosis, which may facilitate symbiont digestion and account for extracellular symbiotic phenotypes. We also reveal that distinct immune system evolution, including expansion in lipopolysaccharide scavenging and contraction of IAP (inhibitor of apoptosis protein), may contribute to the different manners of bacterial virulence resistance in C. bisecta.

CONCLUSIONS

Thus, bivalves employ different pathways to adapt to the long-term co-existence with their bacterial symbionts, further highlighting the contribution of stochastic evolution to the independent gain of a symbiotic lifestyle in the lineage.