Evolutionary transcriptomics of metazoan biphasic life cycle supports a single intercalation origin of metazoan larvae

Molluscan systematics: historical perspectives and the way ahead

Mollusca, the second-most diverse animal phylum, is estimated to have over 100,000 living species with great genetic and phenotypic diversity, a rich fossil record, and a considerable evolutionary significance. Early work on molluscan systematics was grounded in morphological and anatomical studies. With the transition from oligo gene Sanger sequencing to cutting-edge genomic sequencing technologies, molecular data has been increasingly utilised, providing abundant information for reconstructing the molluscan phylogenetic tree. However, relationships among and within most major lineages of Mollusca have long been contentious, often due to limited genetic markers, insufficient taxon sampling and phylogenetic conflict. Fortunately, remarkable progress in molluscan systematics has been made in recent years, which has shed light on how major molluscan groups have evolved. In this review of molluscan systematics, we first synthesise the current understanding of the molluscan Tree of Life at higher taxonomic levels. We then discuss how micromolluscs, which have adult individuals with a body size smaller than 5 mm, offer unique insights into Mollusca's vast diversity and deep phylogeny. Despite recent advancements, our knowledge of molluscan systematics and phylogeny still needs refinement. Further advancements in molluscan systematics will arise from integrating comprehensive data sets, including genome-scale data, exceptional fossils, and digital morphological data (including internal structures). Enhanced access to these data sets, combined with increased collaboration among morphologists, palaeontologists, evolutionary developmental biologists, and molecular phylogeneticists, will significantly advance this field.

Neurons with larval synaptic targets pioneer the later nervous system in the annelid Malacoceros fuliginosus

Comparative studies on the development of nervous systems have a significant impact on understanding animal nervous system evolution. Nevertheless, an important question is to what degree neuronal structures, which play an important role in early stages, become part of the adult nervous system or are relevant for its formation. This is likely in many direct developers, but it is not the case in forms with catastrophic metamorphosis. It is not clear in many forms with gradual metamorphosis. This introduces uncertainty in tracing the evolution of nervous systems and of larval forms. One of the prominent larval characteristics of numerous planktonic marine organisms is the epidermal ciliation used for swimming and steering, which disappears during metamorphosis. Therefore, the neuronal elements controlling the ciliary beating are often assumed to vanish with the cilia and regarded as purely larval specializations. With volume EM, we followed the neuronal targets of the very first pioneer neurons at the apical and posterior ends of the larva of the annelid Malacoceros fuliginosus. We observed that all of these pioneers appear to have a dual function. Although they are laying the paths for the later adult nervous system, they also make synaptic contact with the main ciliated ring of the larva. We propose that the formation of the later adult nervous system and the innervation of the larval locomotory organ are indeed closely linked to each other. This has implications for understanding the early nervous system development of marine larvae and for existing hypotheses on nervous system evolution.

MolluscDB 2.0: a comprehensive functional and evolutionary genomics database for over 1400 molluscan species

Mollusca represents the second-largest animal phylum but remains less explored genomically. The increase in high-quality genomes and diverse functional genomic data holds great promise for advancing our understanding of molluscan biology and evolution. To address the opportunities and challenges facing the molluscan research community in managing vast multi-omics resources, we developed MolluscDB 2.0 (http://mgbase.qnlm.ac), which integrates extensive functional genomic data and offers user-friendly tools for multilevel integrative and comparative analyses. MolluscDB 2.0 covers 1450 species across all eight molluscan classes and compiles ∼4200 datasets, making it the most comprehensive multi-omics resource for molluscs to date. MolluscDB 2.0 expands the layers of multi-omics data, including genomes, bulk transcriptomes, single-cell transcriptomes, proteomes, epigenomes and metagenomes. MolluscDB 2.0 also more than doubles the number of functional modules and analytical tools, updating 14 original modules and introducing 20 new, specialized modules. Overall, MolluscDB 2.0 provides highly valuable, open-access multi-omics platform for the molluscan research community, expediting scientific discoveries and deepening our understanding of molluscan biology and evolution.

Tissue-specific, temporal, and core gene-dependent expression patterns of Hsp70s reveal functional allocation in Chlamys farreri under heat stress

Heat shock proteins 70 KDa (Hsp70s) engage in a broad spectrum of cellular functions in response to various stressors. Marine bivalves face substantial threats from the rising seawater temperature attributed to global warming. In the present study, expression patterns of Hsp70s in Zhikong scallop Chlamys farreri (CfHsp70s) were determined in embryos and larvae at all developmental stages, in healthy adult tissues, and across four various tissues exposed to high temperature for acute and chronic periods through in silico analysis. Spatiotemporal expressions suggested CfHsp70s performed specific functional differentiations in scallop's development and growth. Regulatory expression patterns of CfHsp70s, characterized by predominant down-regulation in the mantle, gill and hemocytes, as well as contrasting up-regulation in the heart, suggest differential functional allocation of CfHsp70s among tissues in response to heat stress. Particularly, a core set of 14 CfHsp70s, especially the nine members of the Hsp70B2s, characterized by gene expansion, intron-less structure, shorter gene length, preference for hydrophilic amino acids, and coordinated expression profiles, was predominantly responsible for the inducible up-regulations observed across all four tissue types. Collectively, the tissue-specific, temporal and core gene-dependent expression patterns of CfHsp70s illustrate the functional allocation and molecular evolution of Hsp70 family members in Zhikong scallops.

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.

Single-cell phylotranscriptomics of developmental and cell type evolution

Single-cell phylotranscriptomics is an emerging tool to reveal the molecular and cellular mechanisms of evolution. We summarize its utility in studying the hourglass pattern of ontogenetic evolution and for understanding the evolutionary history of cell types. The developmental hourglass model suggests that the mid-embryonic stage is the most conserved period of development across species, which is supported by morphological and molecular studies. Single-cell phylotranscriptomic analysis has revealed previously underappreciated heterogeneity in transcriptome ages among lineages and cell types throughout development, and has identified the lineages and tissues that drive the whole-organism hourglass pattern. Single-cell transcriptome age analyses also provide important insights into the origin of germ layers, the different selective forces on tissues during adaptation, and the evolutionary relationships between cell types.

'Conga lines' of Ediacaran fronds: insights into the reproductive biology of early metazoans

Late Ediacaran strata from Newfoundland, Canada (~574-560 Ma) document near-census palaeocommunities of some of the earliest metazoans. Such preservation enables reproductive strategies to be inferred from the spatial distribution of populations of fossilized benthic organisms, previously revealing the existence of both propagule and stoloniferous reproductive modes among Ediacaran frondose taxa. Here, we describe 'conga lines': linear arrangements of more than three closely spaced fossil specimens. We calculate probabilistic models of point maps of 13 fossil-bearing bedding surfaces and show that four surfaces contain conga lines that are not the result of chance alignments. We then test whether these features could result from passive pelagic propagules settling in the lee of an existing frond, using computational fluid dynamics and discrete phase modelling. Under Ediacaran palaeoenvironmental conditions, preferential leeside settlement at the spatial scale of the conga lines is unlikely. We therefore conclude that these features are novel and do not reflect previously described reproductive strategies employed by Ediacaran organisms, suggesting the use of mixed reproductive strategies in the earliest animals. Such strategies enabled Ediacaran frondose taxa to act as reproductive generalists and may be an important facet of early metazoan evolution.

Identification of the principal neuropeptide MIP and its action pathway in larval settlement of the echiuran worm Urechis unicinctus

BACKGROUND

Larval settlement and metamorphosis represent critical events in the life history of marine benthic animals. Myoinhibitory peptide (MIP) plays a pivotal role in larval settlement of marine invertebrates. However, the molecular mechanisms of MIP involved in this process are not well understood.

RESULTS

In this study, we evaluated the effects of thirteen MIP mature peptides on triggering the larval settlement of Urechis unicinctus (Xenopneusta, Urechidae), and determined that MIP2 was the principal neuropeptide. Transcriptomic analysis was employed to identify differentially expressed genes (DEGs) between the MIP2-treated larvae and normal early-segmentation larvae. Both cAMP and calcium signaling pathways were enriched in the DEGs of the MIP2-treated larvae, and two neuropeptide receptor genes (Spr, Fmrfar) were up-regulated in the MIP2-treated larvae. The activation of the SPR-cAMP pathway by MIP2 was experimentally validated in HEK293T cells. Furthermore, fourteen cilia-related genes, including Tctex1d2, Cfap45, Ift43, Ift74, Ift22, Cav1 and Mns1, etc. exhibited down-regulated expression in the MIP2-treated larvae. Whole-mount in situ hybridization identified two selected ciliary genes, Tctex1d2 and Cfap45, were specially expressed in circumoral ciliary cells of the early-segmentation larvae. Knocking down Tctex1d2 mRNA levels by in vivo RNA interference significantly increased the larval settlement rate.

CONCLUSION

Our findings suggest that MIP2 inhibits the function of the cilia-related genes, such as Tctex1d2, through the SPR-cAMP-PKA pathway, thereby inducing larval settlement in U. unicinctus. The study contributes important data to the understanding of neuropeptide regulation in larval settlement.

Polymodal sensory perception drives settlement and metamorphosis of Ciona larvae

The Earth's oceans brim with an incredible diversity of microscopic lifeforms, including motile planktonic larvae, whose survival critically depends on effective dispersal in the water column and subsequent exploration of the seafloor to identify a suitable settlement site. How their nervous systems mediate sensing of diverse multimodal cues remains enigmatic. Here, we uncover that the tunicate Ciona intestinalis larvae employ ectodermal sensory cells to sense various mechanical and chemical cues. Combining whole-brain imaging and chemogenetics, we demonstrate that stimuli encoded at the periphery are sufficient to drive global brain-state changes to promote or impede both larval attachment and metamorphosis behaviors. The ability of C. intestinalis larvae to leverage polymodal sensory perception to support information coding and chemotactile behaviors may explain how marine larvae make complex decisions despite streamlined nervous systems.

The Mediterranean mussel Mytilus galloprovincialis: a novel model for developmental studies in mollusks

A model organism in developmental biology is defined by its experimental amenability and by resources created for the model system by the scientific community. For the most powerful invertebrate models, the combination of both has already yielded a thorough understanding of developmental processes. However, the number of developmental model systems is still limited, and their phylogenetic distribution heavily biased. Members of one of the largest animal lineages, the Spiralia, for example, have long been neglected. In order to remedy this shortcoming, we have produced a detailed developmental transcriptome for the bivalve mollusk Mytilus galloprovincialis, and have expanded the list of experimental protocols available for this species. Our high-quality transcriptome allowed us to identify transcriptomic signatures of developmental progression and to perform a first comparison with another bivalve mollusk: the Pacific oyster Crassostrea gigas. To allow co-labelling studies, we optimized and combined protocols for immunohistochemistry and hybridization chain reaction to create high-resolution co-expression maps of developmental genes. The resources and protocols described here represent an enormous boost for the establishment of Mytilus galloprovincialis as an alternative model system in developmental biology.

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.

Single-cell atlases of two lophotrochozoan larvae highlight their complex evolutionary histories

Pelagic larval stages are widespread across animals, yet it is unclear whether larvae were present in the last common ancestor of animals or whether they evolved multiple times due to common selective pressures. Many marine larvae are at least superficially similar; they are small, swim through the beating of bands of cilia, and sense the environment with an apical organ. To understand these similarities, we have generated single-cell atlases for marine larvae from two animal phyla and have compared their cell types. We found clear similarities among ciliary band cells and between neurons of the apical organ in the two larvae pointing to possible homology of these structures, suggesting a single origin of larvae within Spiralia. We also find several clade-specific innovations in each larva, including distinct myocytes and shell gland cells in the oyster larva. Oyster shell gland cells express many recently evolved genes that have made previous gene age estimates for the origin of trochophore larvae too young.

Morphological Observation and Transcriptome Analysis of Ciliogenesis in Urechis unicinctus (Annelida, Echiura)

During the early development of marine invertebrates, planktic larvae usually occur, and their body surfaces often form specific types of cilia that are involved in locomotion and feeding. The echiuran worm Urechis unicinctus sequentially undergoes the formation and disappearance of different types of body surface cilia during embryonic and larval development. The morphological characteristics and molecular mechanisms involved in the process remain unclear. In this study, we found that body surface cilia in U. unicinctus embryos and larvae can be distinguished into four types: body surface short cilia, apical tufts, circumoral cilia and telotrochs. Further, distribution and genesis of the body surface cilia were characterized using light microscope and electron microscope. To better understand the molecular mechanism during ciliogenesis, we revealed the embryonic and larval transcriptome profile of the key stages of ciliogenesis in U. unicinctus using RNA-Seq technology. A total of 29,158 differentially expressed genes (DEGs) were obtained from 24 cDNA libraries by RNA-Seq. KEGG pathway enrichment results showed that Notch, Wnt and Ca²⁺ signaling pathways were significantly enriched during the occurrence of apical tufts and circumoral cilia. Furthermore, all DEGs were classified according to their expression pattern, and DEGs with similar expression pattern were grouped into a module. All DEG co-expression modules were correlated with traits (body surface short cilia, apical tufts, circumoral cilia and telotrochs) by WGCNA, the results showed DEGs were divided into 13 modules by gene expression patterns and that the genes in No. 7, No. 8 and No. 10 modules were to be highly correlated with the occurrence of apical tufts, circumoral cilia and telotrochs. The top 10 hub genes in the above three modules were identified to be highly correlated with ciliogenesis, including the reported cilium-related gene Cnbd2 and unreported cilium-related candidate genes FAM181B, Capsl, Chst3, TMIE and Innexin. Notably, Innexin was included in the top10 hub genes of the two modules (No. 7 and No. 8), suggesting that Innexin may play an important role in U. unicinctus apical tufts, circumoral cilia and telotrochs genesis. This study revealed the characteristics of ciliogenesis on the body surface of U. unicinctus embryos and larvae, providing basic data for exploring the molecular mechanism of ciliogenesis on the body surface.

Identification and validation of the reference genes in the echiuran worm Urechis unicinctus based on transcriptome data

BACKGROUND

Real-time quantitative PCR (RT-qPCR) is a crucial and widely used method for gene expression analysis. Selecting suitable reference genes is extremely important for the accuracy of RT-qPCR results. Commonly used reference genes are not always stable in various organisms or under different environmental conditions. With the increasing application of high-throughput sequencing, transcriptome analysis has become an effective method for identifying novel stable reference genes.

RESULTS

In this study, we identified candidate reference genes based on transcriptome data covering embryos and larvae of early development, normal adult tissues, and the hindgut under sulfide stress using the coefficient of variation (CV) method in the echiuran Urechis unicinctus, resulting in 6834 (15.82%), 7110 (16.85%) and 13880 (35.87%) candidate reference genes, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the candidate reference genes were significantly enriched in cellular metabolic process, protein metabolic process and ribosome in early development and normal adult tissues as well as in cellular localization and endocytosis in the hindgut under sulfide stress. Subsequently, ten genes including five new candidate reference genes and five commonly used reference genes, were validated by RT-qPCR. The expression stability of the ten genes was analyzed using four methods (geNorm, NormFinder, BestKeeper, and ∆Ct). The comprehensive results indicated that the new candidate reference genes were more stable than most commonly used reference genes. The commonly used ACTB was the most unstable gene. The candidate reference genes STX12, EHMT1, and LYAG were the most stable genes in early development, normal adult tissues, and hindgut under sulfide stress, respectively. The log2(TPM) of the transcriptome data was significantly negatively correlated with the Ct values of RT-qPCR (Ct =  - 0.5405 log2(TPM) + 34.51), which made it possible to estimate the Ct value before RT-qPCR using transcriptome data.

CONCLUSION

Our study is the first to select reference genes for RT-qPCR from transcriptome data in Echiura and provides important information for future gene expression studies in U. unicinctus.

Annelid functional genomics reveal the origins of bilaterian life cycles

Indirect development with an intermediate larva exists in all major animal lineages1, which makes larvae central to most scenarios of animal evolution2-11. Yet how larvae evolved remains disputed. Here we show that temporal shifts (that is, heterochronies) in trunk formation underpin the diversification of larvae and bilaterian life cycles. We performed chromosome-scale genome sequencing in the annelid Owenia fusiformis with transcriptomic and epigenomic profiling during the life cycles of this and two other annelids. We found that trunk development is deferred to pre-metamorphic stages in the feeding larva of O. fusiformis but starts after gastrulation in the non-feeding larva with gradual metamorphosis of Capitella teleta and the direct developing embryo of Dimorphilus gyrociliatus. Accordingly, the embryos of O. fusiformis develop first into an enlarged anterior domain that forms larval tissues and the adult head12. Notably, this also occurs in the so-called 'head larvae' of other bilaterians13-17, with which the O. fusiformis larva shows extensive transcriptomic similarities. Together, our findings suggest that the temporal decoupling of head and trunk formation, as maximally observed in head larvae, facilitated larval evolution in Bilateria. This diverges from prevailing scenarios that propose either co-option9,10 or innovation11 of gene regulatory programmes to explain larva and adult origins.

From head to rootlet: comparative transcriptomic analysis of a rhizocephalan barnacle Peltogaster reticulata (Crustacea: Rhizocephala)

Background: Rhizocephalan barnacles stand out in the diverse world of metazoan parasites. The body of a rhizocephalan female is modified beyond revealing any recognizable morphological features, consisting of the interna, a system of rootlets, and the externa, a sac-like reproductive body. Moreover, rhizocephalans have an outstanding ability to control their hosts, literally turning them into "zombies". Despite all these amazing traits, there are no genomic or transcriptomic data about any Rhizocephala. Methods: We collected transcriptomes from four body parts of an adult female rhizocephalan Peltogaster reticulata: the externa, and the main, growing, and thoracic parts of the interna. We used all prepared data for the de novo assembly of the reference transcriptome. Next, a set of encoded proteins was determined, the expression levels of protein-coding genes in different parts of the parasite's body were calculated and lists of enriched bioprocesses were identified. We also in silico identified and analyzed sets of potential excretory / secretory proteins. Finally, we applied phylostratigraphy and evolutionary transcriptomics approaches to our data.  Results: The assembled reference transcriptome included transcripts of 12,620 protein-coding genes and was the first for any rhizocephalan. Based on the results obtained, the spatial heterogeneity of protein-coding gene expression in different regions of the adult female body of P. reticulata was established. The results of both transcriptomic analysis and histological studies indicated the presence of germ-like cells in the lumen of the interna. The potential molecular basis of the interaction between the nervous system of the host and the parasite's interna was also determined. Given the prolonged expression of development-associated genes, we suggest that rhizocephalans "got stuck in their metamorphosis", even at the reproductive stage. Conclusions: The results of the first comparative transcriptomic analysis for Rhizocephala not only clarified but also expanded the existing ideas about the biology of these extraordinary parasites.

EDomics: a comprehensive and comparative multi-omics database for animal evo-devo

Evolutionary developmental biology (evo-devo) has been among the most fascinating interdisciplinary fields for decades, which aims to elucidate the origin and evolution of diverse developmental processes. The rapid accumulation of omics data provides unprecedented opportunities to answer many interesting but unresolved evo-devo questions. However, the access and utilization of these resources are hindered by challenges particularly in non-model animals. Here, we establish a comparative multi-omics database for animal evo-devo (EDomics, http://edomics.qnlm.ac) containing comprehensive genomes, bulk transcriptomes, and single-cell data across 40 representative species, many of which are generally used as model organisms for animal evo-devo study. EDomics provides a systematic view of genomic/transcriptomic information from various aspects, including genome assembly statistics, gene features and families, transcription factors, transposable elements, and gene expressional profiles/networks. It also exhibits spatiotemporal gene expression profiles at a single-cell level, such as cell atlas, cell markers, and spatial-map information. Moreover, EDomics provides highly valuable, customized datasets/resources for evo-devo research, including gene family expansion/contraction, inferred core gene repertoires, macrosynteny analysis for karyotype evolution, and cell type evolution analysis. EDomics presents a comprehensive and comparative multi-omics platform for animal evo-devo community to decipher the whole history of developmental evolution across the tree of life.

Renewed perspectives on the sedentary-pelagic last common bilaterian ancestor

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

Signalling molecules inducing metamorphosis in marine organisms

Covering: findings from early 1980s until early 2022Microbial-derived cues of marine biofilms induce settlement and metamorphosis of marine organisms, a process responsible for the emergence of diverse flora and fauna in marine habitats. Although this phenomenon is known for more than 80 years, the research field has only recently gained much momentum. Here, we summarize the currently existing biochemical and microbial knowledge about microbial signalling molecules, con-specific signals, and synthetic compounds that induce or prevent recruitment, settlement, and metamorphosis in invertebrate larvae. We discuss the possible modes of action and conclude with perspectives for future research directions in the field of marine chemical ecology.

From head to rootlet: comparative transcriptomic analysis of a rhizocephalan barnacle Peltogaster reticulata (Crustacea: Rhizocephala)

Background: Rhizocephalan barnacles stand out in the diverse world of metazoan parasites. The body of a rhizocephalan female is modified beyond revealing any recognizable morphological features, consisting of the interna, a system of rootlets, and the externa, a sac-like reproductive body. Moreover, rhizocephalans have an outstanding ability to control their hosts, literally turning them into "zombies". Despite all these amazing traits, there are no genomic or transcriptomic data about any Rhizocephala. Methods: We collected transcriptomes from four body parts of an adult female rhizocephalan Peltogaster reticulata: the externa, and the main, growing, and thoracic parts of the interna. We used all prepared data for the de novo assembly of the reference transcriptome. Next, a set of encoded proteins was determined, the expression levels of protein-coding genes in different parts of the parasite's body were calculated and lists of enriched bioprocesses were identified. We also in silico identified and analyzed sets of potential excretory / secretory proteins. Finally, we applied phylostratigraphy and evolutionary transcriptomics approaches to our data.  Results: The assembled reference transcriptome included transcripts of 12,620 protein-coding genes and was the first for any rhizocephalan. Based on the results obtained, the spatial heterogeneity of protein-coding gene expression in different regions of the adult female body of P. reticulata was established. The results of both transcriptomic analysis and histological studies indicated the presence of germ-like cells in the lumen of the interna. The potential molecular basis of the interaction between the nervous system of the host and the parasite's interna was also determined. Given the prolonged expression of development-associated genes, we suggest that rhizocephalans "got stuck in their metamorphosis", even at the reproductive stage. Conclusions: The results of the first comparative transcriptomic analysis for Rhizocephala not only clarified but also expanded the existing ideas about the biology of these extraordinary parasites.

Ctenophores are direct developers that reproduce continuously beginning very early after hatching

A substantial body of literature reports that ctenophores exhibit an apparently unique life history characterized by biphasic sexual reproduction, the first phase of which is called larval reproduction or dissogeny. Whether this strategy is plastically deployed or a typical part of these species’ life history was unknown. In contrast to previous reports, we show that the ctenophore Mnemiopsis leidyi does not have separate phases of early and adult reproduction, regardless of the morphological transition to what has been considered the adult form. Rather, these ctenophores begin to reproduce at a small body size and spawn continuously from this point onward under adequate environmental conditions. They do not display a gap in productivity for metamorphosis or other physiological transition at a certain body size. Furthermore, nutritional and environmental constraints on fecundity are similar in both small and large animals. Our results provide critical parameters for understanding resource partitioning between growth and reproduction in this taxon, with implications for management of this species in its invaded range. Finally, we report an observation of similarly small-size spawning in a beroid ctenophore, which is morphologically, ecologically, and phylogenetically distinct from other ctenophores reported to spawn at small sizes. We conclude that spawning at small body size should be considered as the default, on-time developmental trajectory rather than as precocious, stress-induced, or otherwise unusual for ctenophores. The ancestral ctenophore was likely a direct developer, consistent with the hypothesis that multiphasic life cycles were introduced after the divergence of the ctenophore lineage.

Decoding the byssus fabrication by spatiotemporal secretome analysis of scallop foot

•

The first secretome about scallop byssal adhesion is profiled based on a new computational strategy.

•

Scallop byssal secretome covered almost all of the known structural elements and functional domains of aquatic adhesives.

•

The EGF-like domain containing proteins, the Tyr-rich proteins and 4C-repeats containing proteins are the main components of scallop byssus.

•

A novel “nearby secretion” model of scallop byssus secretion and adhesion is proposed.

Secretome is involved in almost all physiological, developmental, and pathological processes, but to date there is still a lack of highly-efficient research strategy to comprehensively study the secretome of invertebrates. Adhesive secretion is a ubiquitous and essential physiological process in aquatic invertebrates with complicated protein components and unresolved adhesion mechanisms, making it a good subject for secretome profiling studies. Here we proposed a computational pipeline for systematic profiling of byssal secretome based on spatiotemporal transcriptomes of scallop. A total of 186 byssus-related proteins (BRPs) were identified, which represented the first characterized secretome of scallop byssal adhesion. Scallop byssal secretome covered almost all of the known structural elements and functional domains of aquatic adhesives, which suggested this secretome-profiling strategy had both high efficiency and accuracy. We revealed the main components of scallop byssus (including EGF-like domain containing proteins, the Tyr-rich proteins and 4C-repeats containing proteins) and the related modification enzymes primarily contributing to the rapid byssus assembly and adhesion. Spatiotemporal expression and co-expression network analyses of BRPs suggested a simultaneous secretion pattern of scallop byssal proteins across the entire region of foot and revealed their diverse functions on byssus secretion. In contrast to the previously proposed “root-initiated secretion and extension-based assembly” model, our findings supported a novel “foot-wide simultaneous secretion and in situ assembly” model of scallop byssus secretion and adhesion. Systematic analysis of scallop byssal secretome provides important clues for understanding the aquatic adhesive secretion process, as well as a common framework for studying the secretome of non-model invertebrates.

Heme-Peroxidase 2 Modulated by POU2F1 and SOX5 is Involved in Pigmentation in Pacific Oyster (Crassostrea gigas)

Color polymorphism is frequently observed in molluscan shellfish, while the molecular regulation of shell pigmentation is not well understood. Peroxidase is a key enzyme involved in melanogenesis. Here, we identified a heme-peroxidase 2 gene (CgHPX2), and characterized the expression patterns and transcriptional regulation of CgHPX2 in the Pacific oyster Crassostrea gigas. Tissues expression analysis showed that CgHPX2 was a mantle-specific gene and primarily expressed in the edge mantle in black shell color oyster compared with white shell oyster. In situ hybridization showed that strong signals for CgHPX2 were detected in the both inner and outer surface of the outer fold of mantle in the black shell color oyster, whereas positive signals in white shell oyster were mainly localized in the outer surface of the outer fold of mantle. In the embryos and larvae, a high expression level of CgHPX2 was detected in the trochophore stage in both black and white shell color oysters. The temporal localization of CgHPX2 was mainly detected in the shell gland and edge mantle of trochophore and calcified shell larvae, respectively. In addition, a 2227 bp of 5' flanking region sequence of CgHPX2 was cloned, which contained a presumed core promoter region and many potential transcription factor binding sites. Further luciferase assay experiment confirmed that POU domain, class 2, transcription factor 1 (POU2F1), and SRY-box transcription factor 5 (SOX5) were involved in transcriptional regulation of CgHPX2 gene through binding to its specific promoter region. After CgPOU2F1 and CgSOX5 RNA interference, the CgHPX2 gene expression was significantly decreased. These results suggested that CgPOU2F1 and CgSOX5 might be two important transcription factors that positively regulated the expression of CgHPX2 gene, improving our understanding of the transcriptional regulation of molluscan shell pigmentation.

A chromosome-scale genome assembly and karyotype of the ctenophore Hormiphora californensis

Here, we present a karyotype, a chromosome-scale genome assembly, and a genome annotation from the ctenophore Hormiphora californensis (Ctenophora: Cydippida: Pleurobrachiidae). The assembly spans 110 Mb in 44 scaffolds and 99.47% of the bases are contained in 13 scaffolds. Chromosome micrographs and Hi-C heatmaps support a karyotype of 13 diploid chromosomes. Hi-C data reveal three large heterozygous inversions on chromosome 1, and one heterozygous inversion shares the same gene order found in the genome of the ctenophore Pleurobrachia bachei. We find evidence that H. californensis and P. bachei share thirteen homologous chromosomes, and the same karyotype of 1n = 13. The manually curated PacBio Iso-Seq-based genome annotation reveals complex gene structures, including nested genes and trans-spliced leader sequences. This chromosome-scale assembly is a useful resource for ctenophore biology and will aid future studies of metazoan evolution and phylogenetics.

Comparative and evolutionary analyses reveal conservation and divergence of the notch pathway in lophotrochozoa

Lophotrochozoan species exhibit wide morphological diversity; however, the molecular basis underlying this diversity remains unclear. Here, we explored the evolution of Notch pathway genes across 37 metazoan species via phylogenetic and molecular evolutionary studies with emphasis on the lophotrochozoans. We displayed the components of Notch pathway in metazoans and found that Delta and Hes/Hey-related genes, as well as their functional domains, are duplicated in lophotrochozoans. Comparative transcriptomics analyses allow us to pinpoint sequence divergence of multigene families in the Notch signalling pathway. We identified the duplication mechanism of a mollusc-specific gene, Delta2, and found it displayed complementary expression throughout development. Furthermore, we found the functional diversification not only in expanded genes in the Notch pathway (Delta and Hes/Hey-related genes), but also in evolutionary conservative genes (Notch, Presenilin, and Su(H)). Together, this comprehensive study demonstrates conservation and divergence within the Notch pathway, reveals evolutionary relationships among metazoans, and provides evidence for the occurrence of developmental diversity in lophotrochozoans, as well as a basis for future gene function studies.

Diversity in the Development of the Neuromuscular System of Nemertean Larvae (Nemertea, Spiralia)

In studies on the development of nervous systems and musculature, fluorescent labeling of neuroactive substances and filamentous actin (f-actin) of muscle cells and the subsequent analysis with confocal laser scanning microscopy (CLSM), has led to a broad comparative data set for the majority of the clades of the superphylum Spiralia. However, a number of clades remain understudied, which results in gaps in our knowledge that drastically hamper the formulation of broad-scale hypotheses on the evolutionary developmental biology (EvoDevo) of the structures in question. Regarding comparative data on the development of the peptidergic nervous system and the musculature of species belonging to the spiralian clade Nemertea (ribbon worms), such considerable knowledge gaps are manifest. This paper presents first findings on fluorescent labeling of the FMRFamide-like component of the nervous system and contributes additional data on the muscle development in the presently still underrepresented larvae of palaeo- and hoplonemertean species. Whereas the architecture of the FMRFamide-like nervous system is comparably uniform between the studied representatives, the formation of the musculature differs considerably, exhibiting developmental modes yet undescribed for any spiralian species. The presented results fill a significant gap in the spiralian EvoDevo data set and thus allow for further elaboration of hypotheses on the ancestral pattern of the musculature and a prominent component of the nervous system in Nemertea. However, with respect to the variety observed, it is expected that the true diversity of the developmental pathways is still to be discovered when more detailed data on other nemertean species will be available.

Systematic errors in orthology inference and their effects on evolutionary analyses

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

•

Presence of shared orthologs across species is used for evolutionary analyses

•

We simulated realistic sets of orthologs with no gains or losses

•

Errors predicting shared orthologs correlate with phylogenetic relationships

•

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

MolluscDB: an integrated functional and evolutionary genomics database for the hyper-diverse animal phylum Mollusca

Mollusca represents the second largest animal phylum but remains poorly explored from a genomic perspective. While the recent increase in genomic resources holds great promise for a deep understanding of molluscan biology and evolution, access and utilization of these resources still pose a challenge. Here, we present the first comprehensive molluscan genomics database, MolluscDB (http://mgbase.qnlm.ac), which compiles and integrates current molluscan genomic/transcriptomic resources and provides convenient tools for multi-level integrative and comparative genomic analyses. MolluscDB enables a systematic view of genomic information from various aspects, such as genome assembly statistics, genome phylogenies, fossil records, gene information, expression profiles, gene families, transcription factors, transposable elements and mitogenome organization information. Moreover, MolluscDB offers valuable customized datasets or resources, such as gene coexpression networks across various developmental stages and adult tissues/organs, core gene repertoires inferred for major molluscan lineages, and macrosynteny analysis for chromosomal evolution. MolluscDB presents an integrative and comprehensive genomics platform that will allow the molluscan community to cope with ever-growing genomic resources and will expedite new scientific discoveries for understanding molluscan biology and evolution.

How can recapitulation be reconciled with modern concepts of evolution?

To understand Haeckel's idea of recapitulation with modern evolutionary biology, one has to realize how evolutionarily conserved embryonic stages appear sequentially in developmental processes as chains of causality. Whether the idea of evolution was accepted or not, Haeckel and von Baer commonly saw an importance of a particularly conserved mid-embryonic stage in biphasic development of metazoans, the phylotype, that defines an animal phylum as the developmental source of a basic body plan. In an evolutionary context, the phylotypic stage was once understood by Haeckel to reflect the common ancestor of animal phyla, which went through hypermorphosis independently into various phyla. Recent comprehensive molecular studies, however, accumulated data to refute this idea. The conserved embryonic pattern does not reflect an ancestral adult morphology but appears to have arisen primarily as an embodiment of developmental constraints established through evolutionary processes. How the developmental burden results in a nested series of constraints will solve the recapitulative tendency of developmental programs.

Conserved Patterns in Developmental Processes and Phases, Rather than Genes, Unite the Highly Divergent Bilateria

Bilateria are the predominant clade of animals on Earth. Despite having evolved a wide variety of body plans and developmental modes, they are characterized by common morphological traits. By default, researchers have tried to link clade-specific genes to these traits, thus distinguishing bilaterians from non-bilaterians, by their gene content. Here we argue that it is rather biological processes that unite Bilateria and set them apart from their non-bilaterian sisters, with a less complex body morphology. To test this hypothesis, we compared proteomes of bilaterian and non-bilaterian species in an elaborate computational pipeline, aiming to search for a set of bilaterian-specific genes. Despite the limited confidence in their bilaterian specificity, we nevertheless detected Bilateria-specific functional and developmental patterns in the sub-set of genes conserved in distantly related Bilateria. Using a novel multi-species GO-enrichment method, we determined the functional repertoire of genes that are widely conserved among Bilateria. Analyzing expression profiles in three very distantly related model species-D. melanogaster, D. rerio and C. elegans-we find characteristic peaks at comparable stages of development and a delayed onset of expression in embryos. In particular, the expression of the conserved genes appears to peak at the phylotypic stage of different bilaterian phyla. In summary, our study illustrate how development connects distantly related Bilateria after millions of years of divergence, pointing to processes potentially separating them from non-bilaterians. We argue that evolutionary biologists should return from a purely gene-centric view of evolution and place more focus on analyzing and defining conserved developmental processes and periods.

The evo-devo of molluscs: Insights from a genomic perspective

Molluscs represent one of ancient and evolutionarily most successful groups of marine invertebrates, with a tremendous diversity of morphology, behavior, and lifestyle. Molluscs are excellent subjects for evo-devo studies; however, understanding of the evo-devo of molluscs has been largely hampered by incomplete fossil records and limited molecular data. Recent advancement of genomics and other technologies has greatly fueled the molluscan "evo-devo" field, and decoding of several molluscan genomes provides unprecedented insights into molluscan biology and evolution. Here, we review the recent progress of molluscan genome sequencing as well as novel insights gained from their genomes, by emphasizing how molluscan genomics enhances our understanding of the evo-devo of molluscs.