Novel Genes, Ancient Genes, and Gene Co-Option Contributed to the Genetic Basis of the Radula, a Molluscan Innovation

Exploring proteins within the coccolith matrix

Coccolithophores comprise a major component of the oceanic carbon cycle. These unicellular algae produce ornate structures made of calcium carbonate, termed coccoliths, representing ~ 50% of calcite production in the open ocean. The exact molecular mechanisms which direct and control coccolith formation are unknown. In this study, we report on the presence and functional features of proteins within the coccoliths produced by a range of model coccolithophore species including: the globally abundant and well-studied Gephyrocapsa huxleyi (formerly Emiliania huxleyi) and related Gephyrocapsa oceanica, as well as the larger and more heavily calcified Coccolithus braarudii. Protein features were compared between species and against biomineralisation proteins previously identified in other marine calcifying organisms. Notably, several protein features were consistently seen across the examined coccolithophore species, including the cell signalling 14-3-3 domain, chromosome segregation SMC ATPase domain, as well as proteins involved in protein processing and protease inhibition. The copper-binding cupredoxin domain was observed in both Gephyrocapsa species, as well as other marine calcifiers, suggestive of a requirement of Cu in biomineralisation. Building consensus with existing work, we highlight the pentapeptide repeat as a feature which is associated with the coccolith matrix, being identified in all three examined species, and propose that this structural motif may play a role in controlling coccolith growth. This preliminary study provides insight towards the functional diversity of calcification machinery in coccolithophores and presents a number of candidates for future research towards understanding the biochemical controls which direct coccolithogenesis.

The role of heterochronic gene expression and regulatory architecture in early developmental divergence

New developmental programs can evolve through adaptive changes to gene expression. The annelid Streblospio benedicti has a developmental dimorphism, which provides a unique intraspecific framework for understanding the earliest genetic changes that take place during developmental divergence. Using comparative RNAseq through ontogeny, we find that only a small proportion of genes are differentially expressed at any time, despite major differences in larval development and life history. These genes shift expression profiles across morphs by either turning off any expression in one morph or changing the timing or amount of gene expression. We directly connect the contributions of these mechanisms to differences in developmental processes. We examine F1 offspring - using reciprocal crosses - to determine maternal mRNA inheritance and the regulatory architecture of gene expression. These results highlight the importance of both novel gene expression and heterochronic shifts in developmental evolution, as well as the trans-acting regulatory factors in initiating divergence.

The paradise fish, an advanced animal model for behavioral genetics and evolutionary developmental biology

Paradise fish (Macropodus opercularis) is an air-breathing freshwater fish species with a signature labyrinth organ capable of extracting oxygen from the air that helps these fish to survive in hypoxic environments. The appearance of this evolutionary innovation in anabantoids resulted in a rewired circulatory system, but also in the emergence of species-specific behaviors, such as territorial display, courtship and parental care in the case of the paradise fish. Early zoologists were intrigued by the structure and function of the labyrinth apparatus and a series of detailed descriptive histological studies at the beginning of the 20th century revealed the ontogenesis and function of this specialized system. A few decades later, these fish became the subject of numerous ethological studies, and detailed ethograms of their behavior were constructed. These latter studies also demonstrated a strong genetic component underlying their behavior, but due to lack of adequate molecular tools, the fine genetic dissection of the behavior was not possible at the time. The technological breakthroughs that transformed developmental biology and behavioral genetics in the past decades, however, give us now a unique opportunity to revisit these old questions. Building on the classic descriptive studies, the new methodologies will allow us to follow the development of the labyrinth apparatus at a cellular resolution, reveal the genes involved in this process and also the genetic architecture behind the complex behaviors that we can observe in this species.

Mitogenomes do not substantially improve phylogenetic resolution in a young non-model adaptive radiation of freshwater gastropods

BACKGROUND

Species flocks in ancient lakes, and particularly those arising from adaptive radiation, make up the bulk of overall taxonomic and morphological diversity in these insular ecosystems. For these mostly young species assemblages, classical mitochondrial barcoding markers have so far been key to disentangle interspecific relationships. However, with the rise and further development of next-generation sequencing (NGS) methods and mapping tools, genome-wide data have become an increasingly important source of information even for non-model groups.

RESULTS

Here, we provide, for the first time, a comprehensive mitogenome dataset of freshwater gastropods endemic to Sulawesi and thus of an ancient lake invertebrate species flock in general. We applied low-coverage whole-genome sequencing for a total of 78 individuals including 27 out of the 28 Tylomelania morphospecies from the Malili lake system as well as selected representatives from Lake Poso and adjacent catchments. Our aim was to assess whether mitogenomes considerably contribute to the phylogenetic resolution within this young species flock. Interestingly, we identified a high number of variable and parsimony-informative sites across the other 'non-traditional' mitochondrial loci. However, although the overall support was very high, the topology obtained was largely congruent with previously published single-locus phylogenies. Several clades remained unresolved and a large number of species was recovered polyphyletic, indicative of both rapid diversification and mitochondrial introgression.

CONCLUSIONS

This once again illustrates that, despite the higher number of characters available, mitogenomes behave like a single locus and thus can only make a limited contribution to resolving species boundaries, particularly when introgression events are involved.

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.

Co-option of an Astacin Metalloprotease Is Associated with an Evolutionarily Novel Feeding Morphology in a Predatory Nematode

Animals consume a wide variety of food sources to adapt to different environments. However, the genetic mechanisms underlying the acquisition of evolutionarily novel feeding morphology remain largely unknown. While the nematode Caenorhabditis elegans feeds on bacteria, the satellite species Pristionchus pacificus exhibits predatory feeding behavior toward other nematodes, which is an evolutionarily novel feeding habit. Here, we found that the astacin metalloprotease Ppa-NAS-6 is required for the predatory killing by P. pacificus. Ppa-nas-6 mutants were defective in predation-associated characteristics, specifically the tooth morphogenesis and tooth movement during predation. Comparison of expression patterns and rescue experiments of nas-6 in P. pacificus and C. elegans suggested that alteration of the spatial expression patterns of NAS-6 may be vital for acquiring predation-related traits. Reporter analysis of the Ppa-nas-6 promoter in C. elegans revealed that the alteration in expression patterns was caused by evolutionary changes in cis- and trans-regulatory elements. This study suggests that the co-option of a metalloprotease is involved in an evolutionarily novel feeding morphology.

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

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 has 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. 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. 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.

The draft genome of the microscopic Nemertoderma westbladi sheds light on the evolution of Acoelomorpha genomes

Background: Xenacoelomorpha is a marine clade of microscopic worms that is an important model system for understanding the evolution of key bilaterian novelties, such as the excretory system. Nevertheless, Xenacoelomorpha genomics has been restricted to a few species that either can be cultured in the lab or are centimetres long. Thus far, no genomes are available for Nemertodermatida, one of the group's main clades and whose origin has been dated more than 400 million years ago. Methods: DNA was extracted from a single specimen and sequenced with HiFi following the PacBio Ultra-Low DNA Input protocol. After genome assembly, decontamination, and annotation, the genome quality was benchmarked using two acoel genomes and one Illumina genome as reference. The gene content of three cnidarians, three acoelomorphs, four deuterostomes, and eight protostomes was clustered in orthogroups to make inferences of gene content evolution. Finally, we focused on the genes related to the ultrafiltration excretory system to compare patterns of presence/absence and gene architecture among these clades. Results: We present the first nemertodermatid genome sequenced from a single specimen of Nemertoderma westbladi. Although genome contiguity remains challenging (N50: 60 kb), it is very complete (BUSCO: 80.2%, Metazoa; 88.6%, Eukaryota) and the quality of the annotation allows fine-detail analyses of genome evolution. Acoelomorph genomes seem to be relatively conserved in terms of the percentage of repeats, number of genes, number of exons per gene and intron size. In addition, a high fraction of genes present in both protostomes and deuterostomes are absent in Acoelomorpha. Interestingly, we show that all genes related to the excretory system are present in Xenacoelomorpha except Osr, a key element in the development of these organs and whose acquisition seems to be interconnected with the origin of the specialised excretory system. Conclusion: Overall, these analyses highlight the potential of the Ultra-Low Input DNA protocol and HiFi to generate high-quality genomes from single animals, even for relatively large genomes, making it a feasible option for sequencing challenging taxa, which will be an exciting resource for comparative genomics analyses.

Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

Multiple gene co-options underlie the rapid evolution of sexually deceptive flowers in Gorteria diffusa

Gene co-option, the redeployment of an existing gene in an unrelated developmental context, is an important mechanism underlying the evolution of morphological novelty. In most cases described to date, novel traits emerged by co-option of a single gene or genetic network. Here, we show that the integration of multiple co-opted genetic elements facilitated the rapid evolution of complex petal spots that mimic female bee-fly pollinators in the sexually deceptive South African daisy Gorteria diffusa. First, co-option of iron homeostasis genes altered petal spot pigmentation, producing a color similar to that of female pollinators. Second, co-option of the root hair gene GdEXPA7 enabled the formation of enlarged papillate petal epidermal cells, eliciting copulation responses from male flies. Third, co-option of the miR156-GdSPL1 transcription factor module altered petal spot placement, resulting in better mimicry of female flies resting on the flower. The three genetic elements were likely co-opted sequentially, and strength of sexual deception in different G. diffusa floral forms strongly correlates with the presence of the three corresponding morphological alterations. Our findings suggest that gene co-options can combine in a modular fashion, enabling rapid evolution of novel complex traits.

The ontogeny of elements: distinct ontogenetic patterns in the radular tooth mineralization of gastropods

The molluscan phylum is characterized by the radula, used for the gathering and processing of food. This structure can consist of a chitinous membrane with embedded rows of teeth, which show structural, chemical, and biomechanical adaptations to the preferred ingesta. With regard to the chemical composition of teeth, some taxa (Polyplacophora and Patellogastropoda) were extensively studied, and high proportions of incorporated iron, calcium, and silicon were previously reported. However, outside these two groups, there is an immense lack of knowledge about the elemental composition of radular teeth. The here presented work aims at shedding some light on the radular composition by performing energy-dispersive X-ray spectroscopy (EDX) on six non-patelliform gastropod species (Anentome helena, Cornu aspersum, Lavigeria nassa, Littorina littorea, Reymondia horei, and Vittina turrita), with the focus on the ontogeny of the elemental composition. Proportions of elements, which are not part of chitin and other purely organic molecules, were documented for overall 1027 individual teeth of all ontogenetic radular stages, i.e., for the building zone, the maturation zone, and the working zone. We detected that the proportions of these elements increased from the building to the maturation zone. However, from the maturation to the working zone, two general trends are visible: either the proportions of the elements increased or decreased. The latter trend could potentially be explained by the acidic pH of the gastropod saliva, which awaits further investigations.

Gene transcriptional profiles in gonads of Bacillus taxa (Phasmida) with different cytological mechanisms of automictic parthenogenesis

The evolution of automixis - i.e., meiotic parthenogenesis - requires several features, including ploidy restoration after meiosis and maintenance of fertility. Characterizing the relative contribution of novel versus pre-existing genes and the similarities in their expression and sequence evolution is fundamental to understand the evolution of reproductive novelties. Here we identify gonads-biased genes in two Bacillus automictic stick-insects and compare their expression profile and sequence evolution with a bisexual congeneric species. The two parthenogens restore ploidy through different cytological mechanisms: in Bacillus atticus, nuclei derived from the first meiotic division fuse to restore a diploid egg nucleus, while in Bacillus rossius, diploidization occurs in some cells of the haploid blastula through anaphase restitution. Parthenogens' gonads transcriptional program is found to be largely assembled from genes that were already present before the establishment of automixis. The three species transcriptional profiles largely reflect their phyletic relationships, yet we identify a shared core of genes with gonad-biased patterns of expression in parthenogens which are either male gonads-biased in the sexual species or are not differentially expressed there. At the sequence level, just a handful of gonads-biased genes were inferred to have undergone instances of positive selection exclusively in the parthenogen species. This work is the first to explore the molecular underpinnings of automixis in a comparative framework: it delineates how reproductive novelties can be sustained by genes whose origin precedes the establishment of the novelty itself and shows that different meiotic mechanisms of reproduction can be associated with a shared molecular ground plan.

Biomimetic generation of the strongest known biomaterial found in limpet tooth

The biomaterial with the highest known tensile strength is a unique composite of chitin and goethite (α-FeO(OH)) present in teeth from the Common Limpet (Patella vulgata). A biomimetic based on limpet tooth, with corresponding high-performance mechanical properties is highly desirable. Here we report on the replication of limpet tooth developmental processes ex vivo, where isolated limpet tissue and cells in culture generate new biomimetic structures. Transcriptomic analysis of each developmental stage of the radula, the organ from which limpet teeth originate, identifies sequential changes in expression of genes related to chitin and iron processing. We quantify iron and chitin metabolic processes in the radula and grow isolated radula cells in vitro. Bioinspired material can be developed with electrospun chitin mineralised by conditioned media from cultured radula cells. Our results inform molecular processes behind the generation of limpet tooth and establish a platform for development of a novel biomimetic with comparable properties.

The highest tensile strength biomaterial known exists in limpet teeth and replicating this material is of interest. Here, the authors report on the ex vivo growth of teeth and use of isolated limpet tissue and cells providing foundations for the development of this high-tensile biomaterial.

Evolutionary Divergence and Radula Diversification in Two Ecomorphs from an Adaptive Radiation of Freshwater Snails

(1) Background: Adaptive diversification of complex traits plays a pivotal role in the evolution of organismal diversity. In the freshwater snail genus Tylomelania, adaptive radiations were likely promoted by trophic specialization via diversification of their key foraging organ, the radula. (2) Methods: To investigate the molecular basis of radula diversification and its contribution to lineage divergence, we used tissue-specific transcriptomes of two sympatric Tylomelania sarasinorum ecomorphs. (3) Results: We show that ecomorphs are genetically divergent lineages with habitat-correlated abundances. Sequence divergence and the proportion of highly differentially expressed genes are significantly higher between radula transcriptomes compared to the mantle and foot. However, the same is not true when all differentially expressed genes or only non-synonymous SNPs are considered. Finally, putative homologs of some candidate genes for radula diversification (hh, arx, gbb) were also found to contribute to trophic specialization in cichlids and Darwin’s finches. (4) Conclusions: Our results are in line with diversifying selection on the radula driving Tylomelania ecomorph divergence and indicate that some molecular pathways may be especially prone to adaptive diversification, even across phylogenetically distant animal groups.

Shell features and anatomy of the springsnail genus Radomaniola (Caenogastropoda: Hydrobiidae) show a different pace and mode of evolution over five million years

Shell features are key factors for systematics and evolutionary biology studies of freshwater gastropods. Based mostly on shell morphology, 19 species of the springsnail genus Radomaniola (family Hydrobiidae) have been described. Although the scarce differentiation of various shell dimensions suggests morphological stasis, the evolutionary dynamics of shell and other anatomical structures have not yet been explored fully. By inferring a phylogeny and a time-calibrated species tree from mitochondrial (COI and 16S) and nuclear (28S) sequences of 15 recognized species, we initially examine the species diversity of the Radomaniola dataset and then, through phylogenetic comparative methods, assess the evolutionary correlation, pace and mode of 40 continuous shell and anatomical characters. By synthesizing the results of four species delimitation methods and through morphological examinations, we recognize 21 groups, for which taxonomy is discussed here. Seven new species are described. We reveal a high degree of correlated evolution between characters of the shell, which are constrained by a single morphological optimum, consistent with the morphological stasis model. Anatomical traits diverged rapidly in an unconstrained manner or toward multiple optima. These findings indicate that in order to understand the morphological evolution of springsnails, it is essential to examine different organs in detail.

Gene family evolution underlies cell-type diversification in the hypothalamus of teleosts

Hundreds of cell types form the vertebrate brain but it is largely unknown how similar cellular repertoires are between or within species or how cell-type diversity evolves. To examine cell-type diversity across and within species, we performed single-cell RNA sequencing of ~130,000 hypothalamic cells from zebrafish (Danio rerio) and surface and cave morphs of Mexican tetra (Astyanax mexicanus). We found that over 75% of cell types were shared between zebrafish and Mexican tetra, which diverged from a common ancestor over 150 million years ago. Shared cell types displayed shifts in paralogue expression that were generated by subfunctionalization after genome duplication. Expression of terminal effector genes, such as neuropeptides, was more conserved than the expression of their associated transcriptional regulators. Species-specific cell types were enriched for the expression of species-specific genes and characterized by the neofunctionalization of expression patterns of members of recently expanded or contracted gene families. Comparisons between surface and cave morphs revealed differences in immune repertoires and transcriptional changes in neuropeptidergic cell types associated with genomic differences. The single-cell atlases presented here are a powerful resource to explore hypothalamic cell types and reveal how gene family evolution and shifts in paralogue expression contribute to cellular diversity.

Inflammation and convergent placenta gene co-option contributed to a novel reproductive tissue

The evolution of pregnancy exposes parental tissues to new, potentially stressful conditions, which can trigger inflammation.¹ Inflammation is costly^2,3 and can induce embryo rejection, which constrains the evolution of pregnancy.¹ In contrast, inflammation can also promote morphological innovation at the maternal-embryonic interface as exemplified by co-option of pro-inflammatory signaling for eutherian embryo implantation.^1,4,5 Given its dual function, inflammation could be a key process explaining how innovations such as pregnancy and placentation evolved many times convergently. Pelvic brooding ricefishes evolved a novel “plug” tissue,^6,7 which forms inside the female gonoduct after spawning, anchors egg-attaching filaments, and enables pelvic brooders to carry eggs externally until hatching.^6,8 Compared to pregnancy, i.e., internal bearing of embryos, external bearing should alleviate constraints on inflammation in the reproductive tract. We thus hypothesized that an ancestral inflammation triggered by the retention of attaching filaments gave rise to pathways orchestrating plug formation. In line with our hypothesis, histological sections of the developing plug revealed signs of gonoduct injuries by egg-attaching filaments in the pelvic brooding ricefish Oryzias eversi. Tissue-specific transcriptomes showed that inflammatory signaling dominates the plug transcriptome and inflammation-induced genes controlling vital processes for plug development such as tissue growth and angiogenesis were overexpressed in the plug. Finally, mammalian placenta genes were enriched in the plug transcriptome, indicating convergent gene co-option for building, attaching, and sustaining a transient tissue in the female reproductive tract. This study highlights the role of gene co-option and suggests that recruiting inflammatory signaling into physiological processes provides a fast-track to evolutionary innovation.

•

Pelvic brooding induces tissue-specific changes in gene expression

•

Inflammatory signaling characterizes transcriptome of the egg-anchoring plug

•

Similar to embryo implantation, the plug likely evolved from an inflammatory response

•

Mammalian placenta genes were independently co-opted into the plug

Molecular and morphological analysis of the developing nemertean brain indicates convergent evolution of complex brains in Spiralia

BACKGROUND

The brain anatomy in the clade Spiralia can vary from simple, commissural brains (e.g., gastrotrichs, rotifers) to rather complex, partitioned structures (e.g., in cephalopods and annelids). How often and in which lineages complex brains evolved still remains unclear. Nemerteans are a clade of worm-like spiralians, which possess a complex central nervous system (CNS) with a prominent brain, and elaborated chemosensory and neuroglandular cerebral organs, which have been previously suggested as homologs to the annelid mushroom bodies. To understand the developmental and evolutionary origins of the complex brain in nemerteans and spiralians in general, we investigated details of the neuroanatomy and gene expression in the brain and cerebral organs of the juveniles of nemertean Lineus ruber.

RESULTS

In the juveniles, the CNS is already composed of all major elements present in the adults, including the brain, paired longitudinal lateral nerve cords, and an unpaired dorsal nerve cord, which suggests that further neural development is mostly related with increase in the size but not in complexity. The ultrastructure of the juvenile cerebral organ revealed that it is composed of several distinct cell types present also in the adults. The 12 transcription factors commonly used as brain cell type markers in bilaterians show region-specific expression in the nemertean brain and divide the entire organ into several molecularly distinct areas, partially overlapping with the morphological compartments. Additionally, several of the mushroom body-specific genes are expressed in the developing cerebral organs.

CONCLUSIONS

The dissimilar expression of molecular brain markers between L. ruber and the annelid Platynereis dumerilii indicates that the complex brains present in those two species evolved convergently by independent expansions of non-homologous regions of a simpler brain present in their last common ancestor. Although the same genes are expressed in mushroom bodies and cerebral organs, their spatial expression within organs shows apparent differences between annelids and nemerteans, indicating convergent recruitment of the same genes into patterning of non-homologous organs or hint toward a more complicated evolutionary process, in which conserved and novel cell types contribute to the non-homologous structures.

Genomic signatures supporting the symbiosis and formation of chitinous tube in the deep-sea tubeworm Paraescarpia echinospica

Vestimentiferan tubeworms are iconic animals that present as large habitat-forming chitinized tube bushes in deep-sea chemosynthetic ecosystems. They are gutless and depend entirely on their endosymbiotic sulfide-oxidizing chemoautotrophic bacteria for nutrition. Information on the genomes of several siboglinid endosymbionts has improved our understanding of their nutritional supplies. However, the interactions between tubeworms and their endosymbionts remain largely unclear due to a paucity of host genomes. Here, we report the chromosome-level genome of the vestimentiferan tubeworm Paraescarpia echinospica. We found that the genome has been remodeled to facilitate symbiosis through the expansion of gene families related to substrate transfer and innate immunity, suppression of apoptosis, regulation of lysosomal digestion, and protection against oxidative stress. Furthermore, the genome encodes a programmed cell death pathway that potentially controls the endosymbiont population. Our integrated genomic, transcriptomic, and proteomic analyses uncovered matrix proteins required for the formation of the chitinous tube and revealed gene family expansion and co-option as evolutionary mechanisms driving the acquisition of this unique supporting structure for deep-sea tubeworms. Overall, our study provides novel insights into the host’s support system that has enabled tubeworms to establish symbiosis, thrive in deep-sea hot vents and cold seeps, and produce the unique chitinous tubes in the deep sea.

The genome of Nautilus pompilius illuminates eye evolution and biomineralization

Nautilus is the sole surviving externally shelled cephalopod from the Palaeozoic. It is unique within cephalopod genealogy and critical to understanding the evolutionary novelties of cephalopods. Here, we present a complete Nautilus pompilius genome as a fundamental genomic reference on cephalopod innovations, such as the pinhole eye and biomineralization. Nautilus shows a compact, minimalist genome with few encoding genes and slow evolutionary rates in both non-coding and coding regions among known cephalopods. Importantly, multiple genomic innovations including gene losses, independent contraction and expansion of specific gene families and their associated regulatory networks likely moulded the evolution of the nautilus pinhole eye. The conserved molluscan biomineralization toolkit and lineage-specific repetitive low-complexity domains are essential to the construction of the nautilus shell. The nautilus genome constitutes a valuable resource for reconstructing the evolutionary scenarios and genomic innovations that shape the extant cephalopods.

Comparative Proteomics of Octocoral and Scleractinian Skeletomes and the Evolution of Coral Calcification

Corals are the ecosystem engineers of coral reefs, one of the most biodiverse marine ecosystems. The ability of corals to form reefs depends on the precipitation of calcium carbonate (CaCO₃) under biological control. However, several mechanisms underlying coral biomineralization remain elusive, for example, whether corals employ different molecular machineries to deposit different CaCO₃ polymorphs (i.e., aragonite or calcite). Here, we used tandem mass spectrometry (MS/MS) to compare the proteins occluded in the skeleton of three octocoral and one scleractinian species: Tubipora musica and Sinularia cf. cruciata (calcite sclerites), the blue coral Heliopora coerulea (aragonitic skeleton), and the scleractinian aragonitic Montipora digitata. Reciprocal Blast analysis revealed extremely low overlap between aragonitic and calcitic species, while a core set of proteins is shared between octocorals producing calcite sclerites. However, the carbonic anhydrase CruCA4 is present in the skeletons of both polymorphs. Phylogenetic analysis highlighted several possible instances of protein co-option in octocorals. These include acidic proteins and scleritin, which appear to have been secondarily recruited for calcification and likely derive from proteins playing different functions. Similarities between octocorals and scleractinians included presence of a galaxin-related protein, carbonic anhydrases, and one hephaestin-like protein. Although the first two appear to have been independently recruited, the third appear to share a common origin. This work represents the first attempt to identify and compare proteins associated with coral skeleton polymorph diversity, providing several new research targets and enabling both future functional and evolutionary studies aimed at elucidating the origin and evolution of coral biomineralization.

The hard clam genome reveals massive expansion and diversification of inhibitors of apoptosis in Bivalvia

BACKGROUND

Inhibitors of apoptosis (IAPs) are critical regulators of programmed cell death that are essential for development, oncogenesis, and immune and stress responses. However, available knowledge regarding IAP is largely biased toward humans and model species, while the distribution, function, and evolutionary novelties of this gene family remain poorly understood in many taxa, including Mollusca, the second most speciose phylum of Metazoa.

RESULTS

Here, we present a chromosome-level genome assembly of an economically significant bivalve, the hard clam Mercenaria mercenaria, which reveals an unexpected and dramatic expansion of the IAP gene family to 159 members, the largest IAP gene repertoire observed in any metazoan. Comparative genome analysis reveals that this massive expansion is characteristic of bivalves more generally. Reconstruction of the evolutionary history of molluscan IAP genes indicates that most originated in early metazoans and greatly expanded in Bivalvia through both lineage-specific tandem duplication and retroposition, with 37.1% of hard clam IAPs located on a single chromosome. The expanded IAPs have been subjected to frequent domain shuffling, which has in turn shaped their architectural diversity. Further, we observed that extant IAPs exhibit dynamic and orchestrated expression patterns among tissues and in response to different environmental stressors.

CONCLUSIONS

Our results suggest that sophisticated regulation of apoptosis enabled by the massive expansion and diversification of IAPs has been crucial for the evolutionary success of hard clam and other molluscan lineages, allowing them to cope with local environmental stresses. This study broadens our understanding of IAP proteins and expression diversity and provides novel resources for studying molluscan biology and IAP function and evolution.

The Scaly-foot Snail genome and implications for the origins of biomineralised armour

The Scaly-foot Snail, Chrysomallon squamiferum, presents a combination of biomineralised features, reminiscent of enigmatic early fossil taxa with complex shells and sclerites such as sachtids, but in a recently-diverged living species which even has iron-infused hard parts. Thus the Scaly-foot Snail is an ideal model to study the genomic mechanisms underlying the evolutionary diversification of biomineralised armour. Here, we present a high-quality whole-genome assembly and tissue-specific transcriptomic data, and show that scale and shell formation in the Scaly-foot Snail employ independent subsets of 25 highly-expressed transcription factors. Comparisons with other lophotrochozoan genomes imply that this biomineralisation toolkit is ancient, though expression patterns differ across major lineages. We suggest that the ability of lophotrochozoan lineages to generate a wide range of hard parts, exemplified by the remarkable morphological disparity in Mollusca, draws on a capacity for dynamic modification of the expression and positioning of toolkit elements across the genome.

Genome-wide analysis of short interspersed nuclear elements provides insight into gene and genome evolution in citrus

Short interspersed nuclear elements (SINEs) are non-autonomous retrotransposons that are highly abundant, but not well annotated, in plant genomes. In this study, we identified 41,573 copies of SINEs in seven citrus genomes, including 11,275 full-length copies. The citrus SINEs were distributed among 12 families, with an average full-length rate of 0.27, and were dispersed throughout the chromosomes, preferentially in AT-rich areas. Approximately 18.4% of citrus SINEs were found in close proximity (≤1 kb upstream) to genes, indicating a significant enrichment of SINEs in promoter regions. Citrus SINEs promote gene and genome evolution by offering exons as well as splice sites and start and stop codons, creating novel genes and forming tandem and dispersed repeat structures. Comparative analysis of unique homologous SINE-containing loci (HSCLs) revealed chromosome rearrangements in sweet orange, pummelo, and mandarin, suggesting that unique HSCLs might be valuable for understanding chromosomal abnormalities. This study of SINEs provides us with new perspectives and new avenues by which to understand the evolution of citrus genes and genomes.

Signatures of Divergence, Invasiveness, and Terrestrialization Revealed by Four Apple Snail Genomes

The family Ampullariidae includes both aquatic and amphibious apple snails. They are an emerging model for evolutionary studies due to the high diversity, ancient history, and wide geographical distribution. Insight into drivers of ampullariid evolution is hampered, however, by the lack of genomic resources. Here, we report the genomes of four ampullariids spanning the Old World (Lanistes nyassanus) and New World (Pomacea canaliculata, P. maculata, and Marisa cornuarietis) clades. The ampullariid genomes have conserved ancient bilaterial karyotype features and a novel Hox gene cluster rearrangement, making them valuable in comparative genomic studies. They have expanded gene families related to environmental sensing and cellulose digestion, which may have facilitated some ampullarids to become notorious invasive pests. In the amphibious Pomacea, novel acquisition of an egg neurotoxin and a protein for making the calcareous eggshell may have been key adaptations enabling their transition from underwater to terrestrial egg deposition.

Adaptive divergence in shell morphology in an ongoing gastropod radiation from Lake Malawi

Background

Ecological speciation is a prominent mechanism of diversification but in many evolutionary radiations, particularly in invertebrates, it remains unclear whether supposedly critical ecological traits drove or facilitated diversification. As a result, we lack accurate knowledge on the drivers of diversification for most evolutionary radiations along the tree of life. Freshwater mollusks present an enigmatic example: Putatively adaptive radiations are being described in various families, typically from long-lived lakes, whereas other taxa represent celebrated model systems in the study of ecophenotypic plasticity. Here we examine determinants of shell-shape variation in three nominal species of an ongoing ampullariid radiation in the Malawi Basin (Lanistes nyassanus, L. solidus and Lanistes sp. (ovum-like)) with a common garden experiment and semi-landmark morphometrics.

Results

We found significant differences in survival and fecundity among these species in contrasting habitats. Morphological differences observed in the wild persisted in our experiments for L. nyassanus versus L. solidus and L. sp. (ovum-like), but differences between L. solidus and L. sp. (ovum-like) disappeared and re-emerged in the F₁ and F₂ generations, respectively. These results indicate that plasticity occurred, but that it is not solely responsible for the observed differences. Our experiments provide the first unambiguous evidence for genetic divergence in shell morphology in an ongoing freshwater gastropod radiation in association with marked fitness differences among species under controlled habitat conditions.

Conclusions

Our results indicate that differences in shell morphology among Lanistes species occupying different habitats have an adaptive value. These results also facilitate an accurate reinterpretation of morphological variation in fossil Lanistes radiations, and thus macroevolutionary dynamics. Finally, our work testifies that the shells of freshwater gastropods may retain signatures of adaptation at low taxonomic levels, beyond representing an evolutionary novelty responsible for much of the diversity and disparity in mollusks altogether.

Insights into the Evolution of Shells and Love Darts of Land Snails Revealed from Their Matrix Proteins

Over the past decade, many skeletal matrix proteins that are possibly related to calcification have been reported in various calcifying animals. Molluscs are among the most diverse calcifying animals and some gastropods have adapted to terrestrial ecological niches. Although many shell matrix proteins (SMPs) have already been reported in molluscs, most reports have focused on marine molluscs, and the SMPs of terrestrial snails remain unclear. In addition, some terrestrial stylommatophoran snails have evolved an additional unique calcified character, called a "love dart," used for mating behavior. We identified 54 SMPs in the terrestrial snail Euhadra quaesita, and found that they contain specific domains that are widely conserved in molluscan SMPs. However, our results also suggest that some of them possibly have evolved independently by domain shuffling, domain recruitment, or gene co-option. We then identified four dart matrix proteins, and found that two of them are the same proteins as those identified as SMPs. Our results suggest that some dart matrix proteins possibly have evolved by independent gene co-option from SMPs during dart evolution events. These results provide a new perspective on the evolution of SMPs and "love darts" in land snails.

Hox gene expression in postmetamorphic juveniles of the brachiopod Terebratalia transversa

BACKGROUND

Hox genes encode a family of homeodomain containing transcription factors that are clustered together on chromosomes of many Bilateria. Some bilaterian lineages express these genes during embryogenesis in spatial and/or temporal order according to their arrangement in the cluster, a phenomenon referred to as collinearity. Expression of Hox genes is well studied during embryonic and larval development of numerous species; however, relatively few studies focus on the comparison of pre- and postmetamorphic expression of Hox genes in animals with biphasic life cycle. Recently, the expression of Hox genes was described for embryos and larvae of Terebratalia transversa, a rhynchonelliformean brachiopod, which possesses distinct metamorphosis from planktonic larvae to sessile juveniles. During premetamorphic development, T. transversa does not exhibit spatial collinearity and several of its Hox genes are recruited for the morphogenesis of novel structures. In our study, we determined the expression of Hox genes in postmetamorphic juveniles of T. transversa in order to examine metamorphosis-related changes of expression patterns and to test whether Hox genes are expressed in the spatially collinear way in the postmetamorphic juveniles.

RESULTS

Hox genes are expressed in a spatially non-collinear manner in juveniles, generally showing similar patterns as ones observed in competent larvae: genes labial and post1 are expressed in chaetae-related structures, sex combs reduced in the shell-forming epithelium, whereas lox5 and lox4 in dorso-posterior epidermis. After metamorphosis, expression of genes proboscipedia, hox3, deformed and antennapedia becomes restricted to, respectively, shell musculature, prospective hinge rudiments and pedicle musculature and epidermis.

CONCLUSIONS

All developmental stages of T. transversa, including postmetamorphic juveniles, exhibit a spatial non-collinear Hox genes expression with only minor changes observed between pre- and postmetamorphic stages. Our results are concordant with morphological observation that metamorphosis in rhynchonelliformean brachiopods, despite being rapid, is rather gradual. The most drastic changes in Hox gene expression patterns observed during metamorphosis could be explained by the inversion of the mantle lobe, which relocates some of the more posterior larval structures into the anterior edge of the juveniles. Co-option of Hox genes for the morphogenesis of novel structures is even more pronounced in postmetamorphic brachiopods when compared to larvae.

CRISPR/Cas9 deletions in a conserved exon of Distal-less generates gains and losses in a recently acquired morphological novelty in flies

Distal-less has been repeatedly co-opted for the development of many novel traits. Here, we document its curious role in the development of a novel abdominal appendage (“sternite brushes”) in sepsid flies. CRISPR/Cas9 deletions in the homeodomain result in losses of sternite brushes, demonstrating that Distal-less is necessary for their development. However, deletions in the upstream coding exon (Exon 2) produce losses or gains of brushes. A dissection of Exon 2 reveals that the likely mechanism for gains involves a deletion in an exon-splicing enhancer site that leads to exon skipping. Such contradictory phenotypes are also observed in butterflies, suggesting that mutations in the conserved upstream regions have the potential to generate phenotypic variability in insects that diverged 300 million years ago. Our results demonstrate the importance of Distal-less for the development of a novel abdominal appendage in insects and highlight how site-specific mutations in the same exon can produce contradictory phenotypes.

•

Distal-less is necessary for the development of a novel abdominal appendage

•

CRISPR/Cas9 editing produced both losses and gains of novel abdominal appendages

•

Gains of appendages result from mutations in exonic splicing enhancer (ESEs) sites

•

ESE mutations likely led to exon skipping and an altered Distal-less protein