Expression of Patella vulgata orthologs of engrailed and dpp-BMP2/4 in adjacent domains during molluscan shell development suggests a conserved compartment boundary mechanism

The adult shell matrix protein repertoire of the marine snail Crepidula is dominated by conserved genes that are also expressed in larvae

Mollusca is a morphologically diverse phylum, exhibiting an immense variety of calcium carbonate structures. Proteomic studies of adult shells often report high levels of rapidly-evolving, 'novel' shell matrix proteins (SMPs), which are hypothesized to drive shell diversification. However, relatively little is known about the phylogenetic distribution of SMPs, or about the function of individual SMPs in shell construction. To understand how SMPs contribute to shell diversification a thorough characterization of SMPs is required. Here, we build tools and a foundational understanding of SMPs in the marine gastropod species Crepidula fornicata and Crepidula atrasolea because they are genetically-enabled mollusc model organisms. First, we established a staging system of shell development in C. atrasolea for the first time. Next, we leveraged previous findings in C. fornicata combined with phylogenomic analyses of 95 metazoan species to determine the evolutionary lineage of its adult SMP repertoire. We found that 55% of C. fornicata's SMPs belong to molluscan orthogroups, with 27% restricted to Gastropoda, and only 5% restricted at the species level. The low percentage of species-restricted SMPs underscores the importance of broad-taxon sampling and orthology inference approaches when determining homology of SMPs. From our transcriptome analysis, we found that the majority of C. fornicata SMPs that were found conserved in C. atrasolea were expressed in both larval and adult stages. We then selected a subset of SMPs of varying evolutionary ages for spatial-temporal analysis using in situ hybridization chain reaction (HCR) during larval shell development in C. atrasolea. Out of the 18 SMPs analyzed, 12 were detected in the larval shell field. These results suggest overlapping larval vs. adult SMP repertoires. Using multiplexed HCR, we observed five SMP expression patterns and three distinct cell populations within the shell field. These patterns support the idea that modular expression of SMPs could facilitate divergence of shell morphological characteristics. Collectively, these data establish an evolutionary and developmental framework in Crepidula that enables future comparisons of molluscan biomineralization to reveal mechanisms of shell diversification.

The gastropod Lottia peitaihoensis as a model to study the body patterning of trochophore larvae

The body patterning of trochophore larvae is important for understanding spiralian evolution and the origin of the bilateral body plan. However, considerable variations are observed among spiralian lineages, which have adopted varied strategies to develop trochophore larvae or even omit a trochophore stage. Some spiralians, such as patellogastropod mollusks, are suggested to exhibit ancestral traits by producing equal-cleaving fertilized eggs and possessing "typical" trochophore larvae. In recent years, we developed a potential model system using the patellogastropod Lottia peitaihoensis (= Lottia goshimai). Here, we introduce how the species were selected and establish sources and techniques, including gene knockdown, ectopic gene expression, and genome editing. Investigations on this species reveal essential aspects of trochophore body patterning, including organizer signaling, molecular and cellular processes connecting the various developmental functions of the organizer, the specification and behaviors of the endomesoderm and ectomesoderm, and the characteristic dorsoventral decoupling of Hox expression. These findings enrich the knowledge of trochophore body patterning and have important implications regarding the evolution of spiralians as well as bilateral body plans.

Expression of bone morphogenetic protein 10 and its role in biomineralization in Hyriopsis cumingii

Shells and pearls are the products of biomineralization of shellfish after ingesting external mineral ions. Bone morphogenetic proteins (BMPs) play a role in a variety of biological function, and the genes that encode them, are considered important shell-forming genes in mollusks and are associated with shell and pearl formation, embryonic development, and other functions, but bone morphogenetic protein 10 (BMP10) is poorly understood in Hyriopsis cumingii. In this study, we cloned Hc-BMP10 and obtained a 2477 bp full-length sequence encoding 460 amino acids with a conserved TGF-β structural domain. During the embryonic developmental stages, the cleavage stage had the highest expression of Hc-BMP10, followed by juvenile clams; the expression in the mantle gradually decreased with increasing mussel age. A strong signal was detected on epidermal cells on the mantle edge by in situ hybridization. In both the shell notching and inserting operations of the pearl fragment assay, we found that the expression of Hc-BMP10 increased after the above treatments. RNA interference assays showed that the silencing of Hc-BMP10 resulted in a change in the morphology of the prismatic layer and nacreous layer, with the prismatic layer less closely aligned and the disordered aragonite flakes in the nacreous layer. These findings indicate that Hc-BMP10 is involved in the growth and development of H. cumingii, as well as the formation of shells and pearls. Therefore, this study provides some reference for selecting superior species for growth and pearl breeding of H. cumingii at a molecular level and further investigation of the molecular mechanism for biomineralization of Hc-BMP10.

Cell type and gene regulatory network approaches in the evolution of spiralian biomineralisation

Biomineralisation is the process by which living organisms produce hard structures such as shells and bone. There are multiple independent origins of biomineralised skeletons across the tree of life. This review gives a glimpse into the diversity of spiralian biominerals and what they can teach us about the evolution of novelty. It discusses different levels of biological organisation that may be informative to understand the evolution of biomineralisation and considers the relationship between skeletal and non-skeletal biominerals. More specifically, this review explores if cell type and gene regulatory network approaches could enhance our understanding of the evolutionary origins of biomineralisation.

Molecular and phenotypic effects of early exposure to an environmentally relevant pesticide mixture in the Pacific oyster, Crassostrea gigas

Early life stages are crucial for organism development, especially for those displaying external fertilization, whose gametes and early stages face environmental stressors such as xenobiotics. The pacific oyster, Crassostrea gigas, is considered a model species in ecotoxicology because of its ecological characteristics (benthic, sessile, filter feeding). So far studies have investigated the impact of xenobiotics at embryotoxic, genotoxic and physiological endpoints, sometimes at the multigenerational scale, highlighting the role of epigenetic mechanisms in transmitting alterations induced by exposure to single xenobiotics. However, to date, little is known about the impact of environmentally-mimicking contaminants cocktails. Thus, we examined the impact of an early exposure to environmentally relevant mixture on the Pacific oyster life history. We studied transcriptomic, epigenetic and physiological alterations induced in oysters exposed to 18 pesticides and metals at environmental concentration (nominal sum concentration: 2.85 μg.L-1, measured sum concentration: 3.74 ± 0.013 μg.L-1) during embryo-larval stage (0-48 h post fertilization, hpf). No significant differences in embryo-larval abnormalities at 24 hpf were observed during larval and spat rearing; the swimming behaviour of exposed individuals was disturbed, while they were longer and heavier at specific time points, and exhibited a lower epinephrine-induced metamorphosis rate as well as a higher survival rate in the field. In addition, RNA-seq analyses of gastrula embryos revealed the differential expression of development-related genes (e.g. Hox orthologues and cell cycle regulators) between control and exposed oysters. Whole-genome DNA methylation analyses demonstrated a significant modification of DNA methylation in exposed larvae marked by a demethylation trend. Those findings suggest that early exposure to an environmentally relevant pesticide mixture induces multi-scale latent effects possibly affecting life history traits in the Pacific oyster.

Shell field morphogenesis in the polyplacophoran mollusk Acanthochitona rubrolineata

BACKGROUND

The polyplacophoran mollusks (chitons) possess serially arranged shell plates. This feature is unique among mollusks and believed to be essential to explore the evolution of mollusks as well as their shells. Previous studies revealed several cell populations in the dorsal epithelium (shell field) of polyplacophoran larvae and their roles in the formation of shell plates. Nevertheless, they provide limited molecular information, and shell field morphogenesis remains largely uninvestigated.

RESULTS

In the present study, we investigated shell field development in the chiton Acanthochitona rubrolineata based on morphological characteristics and molecular patterns. A total of four types of tissue could be recognized from the shell field of A. rubrolineata. The shell field comprised not only the centrally located, alternatively arranged plate fields and ridges, but also the tissues surrounding them, which were the precursors of the girdle and we termed as the girdle field. The girdle field exhibited a concentric organization composed of two circularly arranged tissues, and spicules were only developed in the outer circle. Dynamic engrailed expression and F-actin (filamentous actin) distributions revealed relatively complicated morphogenesis of the shell field. The repeated units (plate fields and ridges) were gradually established in the shell field, seemingly different from the manners used in the segmentation of Drosophila or vertebrates. The seven repeated ridges also experienced different modes of ontogenesis from each other. In the girdle field, the presumptive spicule-formation cells exhibited different patterns of F-actin aggregations as they differentiate.

CONCLUSIONS

These results reveal the details concerning the structure of polyplacophoran shell field as well as its morphogenesis. They would contribute to exploring the mechanisms of polyplacophoran shell development and molluscan shell evolution.

Spatial-temporal expression analysis of lineage-restricted shell matrix proteins reveals shell field regionalization and distinct cell populations in the slipper snail Crepidula atrasolea

Molluscs are one of the most morphologically diverse clades of metazoans, exhibiting an immense diversification of calcium carbonate structures, such as the shell. Biomineralization of the calcified shell is dependent on shell matrix proteins (SMPs). While SMP diversity is hypothesized to drive molluscan shell diversity, we are just starting to unravel SMP evolutionary history and biology. Here we leveraged two complementary model mollusc systems, Crepidula fornicata and Crepidula atrasolea , to determine the lineage-specificity of 185 Crepidula SMPs. We found that 95% of the adult C. fornicata shell proteome belongs to conserved metazoan and molluscan orthogroups, with molluscan-restricted orthogroups containing half of all SMPs in the shell proteome. The low number of C. fornicata -restricted SMPs contradicts the generally-held notion that an animal’s biomineralization toolkit is dominated by mostly novel genes. Next, we selected a subset of lineage-restricted SMPs for spatial-temporal analysis using in situ hybridization chain reaction (HCR) during larval stages in C. atrasolea . We found that 12 out of 18 SMPs analyzed are expressed in the shell field. Notably, these genes are present in 5 expression patterns, which define at least three distinct cell populations within the shell field. These results represent the most comprehensive analysis of gastropod SMP evolutionary age and shell field expression patterns to date. Collectively, these data lay the foundation for future work to interrogate the molecular mechanisms and cell fate decisions underlying molluscan mantle specification and diversification.

Evolutionary conservation and divergence of the transcriptional regulation of bivalve shell secretion across life-history stages

Adult molluscs produce shells with diverse morphologies and ornamentations, different colour patterns and microstructures. The larval shell, however, is a phenotypically more conserved structure. How do developmental and evolutionary processes generate varying diversity at different life-history stages within a species? Using live imaging, histology, scanning electron microscopy and transcriptomic profiling, we have described shell development in a heteroconchian bivalve, the Antarctic clam, Laternula elliptica, and compared it to adult shell secretion processes in the same species. Adult downstream shell genes, such as those encoding extracellular matrix proteins and biomineralization enzymes, were largely not expressed during shell development. Instead, a development-specific downstream gene repertoire was expressed. Upstream regulatory genes such as transcription factors and signalling molecules were largely conserved between developmental and adult shell secretion. Comparing heteroconchian data with recently reported pteriomorphian larval shell development data suggests that, despite being phenotypically more conserved, the downstream effectors constituting the larval shell 'tool-kit' may be as diverse as that of adults. Overall, our new data suggest that a larval shell formed using development-specific downstream effector genes is a conserved and ancestral feature of the bivalve lineage, and possibly more broadly across the molluscs.

Identification and functional analysis of two drosophila mothers against decapentaplegic protein(Smad)genes and their involvement in immune responses in the red swamp crayfish Procambarus clarkii

Drosophila mothers against decapentaplegic proteins (Smads), the crucial signal transducers in activating downstream gene transcription through transforming growth factor beta (TGF-β) receptors, are the pleiotropic factors with important role in mediating cell proliferation, homeostasis, differentiation, apoptosis and immune response. However, whether Smads are involved in immune response in crustaceans remains unexplored. In the present study, the Smad3 and Smad4 were firstly identified and functionally characterized from the Red Swamp Crayfish Procambarus clarkii. The full-length cDNAs of pcSmad3 and pcSmad4 were 1, 670 bp and 3, 060 bp with 1, 326 bp and 1, 875 bp open reading frame (ORF), respectively. Real-time PCR analysis of the expression profiles demonstrated that pcSmad3 and pcSmad4 were predominantly expressed at in stomach, heart, and hemocytes. Notably, the expression levels of pcSmad3 and pcSmad4 both Aeromonas hydrophila and WSSV challenges were significantly altered, suggesting the involvement of pcSmad3 and pcSmad4 in innate immune responses. Knockdown of pcSmad3 and pcSmad4 in vivo dramatically activated the transcriptions of NF-κB signaling genes and anti-lipopolysaccharide factor genes. The overexpression of pcSmad3 and pcSmad4 could significantly activate NF-κB signaling in HEK293T cells. Meanwhile, the clearance of bacteria was significantly reduced with knockdown of pcSmad3 and pcSmad4 in vivo. Results indicated that pcSmad3 and pcSmad4 played an immune-regulatory role in crayfish's innate immunity, which might pave the for a better understanding of the TGF-β superfamily members in crustacean.

Brachiopod and mollusc biomineralisation is a conserved process that was lost in the phoronid-bryozoan stem lineage

BACKGROUND

Brachiopods and molluscs are lophotrochozoans with hard external shells which are often believed to have evolved convergently. While palaeontological data indicate that both groups are descended from biomineralising Cambrian ancestors, the closest relatives of brachiopods, phoronids and bryozoans, are mineralised to a much lower extent and are comparatively poorly represented in the Palaeozoic fossil record. Although brachiopod and mollusc shells are structurally analogous, genomic and proteomic evidence indicates that their formation involves a complement of conserved, orthologous genes. Here, we study a set of genes comprised of 3 homeodomain transcription factors, one signalling molecule and 6 structural proteins which are implicated in mollusc and brachiopod shell formation, search for their orthologs in transcriptomes or genomes of brachiopods, phoronids and bryozoans, and present expression patterns of 8 of the genes in postmetamorphic juveniles of the rhynchonelliform brachiopod T. transversa.

RESULTS

Transcriptome and genome searches for the 10 target genes in the brachiopods Terebratalia transversa, Lingula anatina, Novocrania anomala, the bryozoans Bugula neritina and Membranipora membranacea, and the phoronids Phoronis australis and Phoronopsis harmeri resulted in the recovery of orthologs of the majority of the genes in all taxa. While the full complement of genes was present in all brachiopods with a single exception in L. anatina, a bloc of four genes could consistently not be retrieved from bryozoans and phoronids. The genes engrailed, distal-less, ferritin, perlucin, sp1 and sp2 were shown to be expressed in the biomineralising mantle margin of T. transversa juveniles.

CONCLUSIONS

The gene expression patterns we recovered indicate that while mineralised shells in brachiopods and molluscs are structurally analogous, their formation builds on a homologous process that involves a conserved complement of orthologous genes. Losses of some of the genes related to biomineralisation in bryozoans and phoronids indicate that loss of the capacity to form mineralised structures occurred already in the phoronid-bryozoan stem group and supports the idea that mineralised skeletons evolved secondarily in some of the bryozoan subclades.

Comparative Single-Cell Transcriptomics Reveals Novel Genes Involved in Bivalve Embryonic Shell Formation and Questions Ontogenetic Homology of Molluscan Shell Types

Mollusks are known for their highly diverse repertoire of body plans that often includes external armor in form of mineralized hardparts. Representatives of the Conchifera, one of the two major lineages that comprises taxa which originated from a uni-shelled ancestor (Monoplacophora, Gastropoda, Cephalopoda, Scaphopoda, Bivalvia), are particularly relevant regarding the evolution of mollusk shells. Previous studies have found that the shell matrix of the adult shell (teleoconch) is rapidly evolving and that the gene set involved in shell formation is highly taxon-specific. However, detailed annotation of genes expressed in tissues involved in the formation of the embryonic shell (protoconch I) or the larval shell (protoconch II) are currently lacking. Here, we analyzed the genetic toolbox involved in embryonic and larval shell formation in the quagga mussel Dreissena rostriformis using single cell RNA sequencing. We found significant differences in genes expressed during embryonic and larval shell secretion, calling into question ontogenetic homology of these transitory bivalve shell types. Further ortholog comparisons throughout Metazoa indicates that a common genetic biomineralization toolbox, that was secondarily co-opted into molluscan shell formation, was already present in the last common metazoan ancestor. Genes included are engrailed, carbonic anhydrase, and tyrosinase homologs. However, we found that 25% of the genes expressed in the embryonic shell field of D. rostriformis lack an ortholog match with any other metazoan. This indicates that not only adult but also embryonic mollusk shells may be fast-evolving structures. We raise the question as to what degree, and on which taxonomic level, the gene complement involved in conchiferan protoconch formation may be lineage-specific or conserved across taxa.

The role of Smad1/5 in mantle immunity of the pearl oyster Pinctada fucata martensii

The Smad protein family is an important medium for transducing BMP-Smads signals, and which have been proved that their important role in regulating shell biomineralization in Pinctada fucata martensii in our previous study. The members of TGF-β superfamily were involved in innate immunity in vertebrates and invertebrates, and Smad regulatory networks construct a balanced immune system. However, little is known about the role of Smad1/5 in immunity in P. f. martensii. The present study shows that the tissue distribution and the expression profiles of Smad1/5 at developmental stages suggested its wide distribution and crucial role in development at embryonic stages other than larval stage; the increased expression of bone morphogenetic proteins 2 (BMP2), Smad4, Smad1/5 and MSX mRNAs at mantle tissue after LPS and Poly (I:C) challenged implied the potential immune role of Smad1/5 and BMP2-Smad signals to defense against bacterial and virus infections; the reduced expression of immune gene nuclear factor kappa-B (NF-κB), matrix metalloproteinase (MMP), interleukin 17 (IL-17), CuZn-superoxide dismutase (CuZn-SOD), tissue inhibitors of metalloproteinase (TIMP) and lipopolysaccharide-induced TNF-α factor (LITAF) mRNA following knockdown of Smad1/5 indicated that Smad1/5 can regulate their expression via BMP2-Smads pathway in the immunity process; the up-regulated expression of Smad1/5 and BMP2-Smad signals genes, and immune genes during wound healing indicated that Smad1/5 and BMP2-Smad signals genes may be involved in wound healing collaborated with immune genes via a different and complex Smads signaling pathway. These results indicated Smad1/5 could regulate innate immunity via BMP2-Smads signal pathway, and which provided new insights into the relationship between BMP2-Smads signal pathway and mantle immunity.

Mantle Modularity Underlies the Plasticity of the Molluscan Shell: Supporting Data From Cepaea nemoralis

Molluscs have evolved the capacity to fabricate a wide variety of shells over their 540+ million-year history. While modern sequencing and proteomic technologies continue to expand the catalog of molluscan shell-forming proteins, a complete functional understanding of how any mollusc constructs its shell remains an ambitious goal. This lack of understanding also constrains our understanding of how evolution has generated a plethora of molluscan shell morphologies. Taking advantage of a previous expression atlas for shell-forming genes in Lymnaea stagnalis, I have characterized the spatial expression patterns of seven shell-forming genes in the terrestrial gastropod Cepaea nemoralis, with the aim of comparing and contrasting their expression patterns between the two species. Four of these genes were selected from a previous proteomic screen of the C. nemoralis shell, two were targeted by bioinformatics criteria designed to identify likely shell-forming gene products, and the final one was a clear homolog of a peroxidase sequence in the L. stagnalis dataset. While the spatial expression patterns of all seven C. nemoralis genes could be recognized as falling into distinct zones within the mantle tissue similar to those established in L. stagnalis, some zones have apparently been modified. These similarities and differences hint at a modularity to the molluscan mantle that may provide a mechanistic explanation as to how evolution has efficiently generated a diversity of molluscan shells.

The pond snail Lymnaea stagnalis

The freshwater snail Lymnaea stagnalis has a long research history, but only relatively recently has it emerged as an attractive model organism to study molecular mechanisms in the areas of developmental biology and translational medicine such as learning/memory and neurodegenerative diseases. The species has the advantage of being a hermaphrodite and can both cross- and self-mate, which greatly facilitates genetic approaches. The establishment of body-handedness, or chiromorphogenesis, is a major topic of study, since chirality is evident in the shell coiling. Chirality is maternally inherited, and only recently a gene-editing approach identified the actin-related gene Lsdia1 as the key handedness determinant. This short article reviews the natural habitat, life cycle, major research questions and interests, and experimental approaches.

The Increased Expression of an Engrailed to Sustain Shell Formation in Response to Ocean Acidification

Engrailed is a transcription factor required in numerous species for important developmental steps such as neurogenesis, segment formation, preblastoderm organization, and compartment formation. Recent study has proved that engrailed is also a key gene related to shell formation in marine bivalves. In the present study, the expression pattern of an engrailed gene (Cgengrailed-1) in Pacific oyster Crassostrea gigas under CO2-driven acidification was investigated to understand its possible role in the regulation of shell formation and adaptation to ocean acidification (OA). The open reading frame (ORF) of Cgengrailed-1 was obtained, which was of 690 bp encoding a polypeptide of 229 amino acids with a HOX domain. Phylogenetic analysis indicated that the deduced amino acid sequence of Cgengrailed-1 shared high homology with other engraileds from Drosophila melanogaster, Mizuhopecten yessoensi, and Crassostrea virginica. The mRNA transcripts of Cgengrailed-1 were constitutively expressed in various tissues with the highest expression levels detected in labial palp and mantle, which were 86.83-fold (p < 0.05) and 75.87-fold (p < 0.05) higher than that in hepatopancreas. The mRNA expression of Cgengrailed-1 in mantle decreased dramatically after moderate (pH 7.8) and severe (pH 7.4) acidification treatment (0.75- and 0.15-fold of that in control group, p < 0.05). The results of immunofluorescence assay demonstrated that the expression level of Cgengrailed-1 in the middle fold of mantle increased significantly upon moderate and severe acidification treatment. Moreover, after the oyster larvae received acidification treatment at trochophore stage, the mRNA expression levels of Cgengrailed-1 increased significantly in D-shape larvae stages, which was 3.11- (pH 7.8) and 4.39-fold (pH 7.4) of that in control group (p < 0.05). The whole-mount immunofluorescence assay showed that Cgengrailed-1 was mainly expressed on the margin of shell gland, and the periostracum in trochophore, early D-shape larvae and D-shape larvae in both control and acidification treatment groups, and the intensity of positive signals in early D-shape larvae and D-shape larvae increased dramatically under acidification treatment. These results collectively suggested that the expression of Cgengrailed-1 could be triggered by CO2-driven acidification treatment, which might contribute to induce the initial shell formation in oyster larvae and the formation of periostracum in adult oyster to adapt to the acidifying marine environment.

Ocean acidification inhibits initial shell formation of oyster larvae by suppressing the biosynthesis of serotonin and dopamine

Ocean acidification has severely affected the initial shell formation of marine bivalves during their larval stages. In the present study, it was found that dopamine (DA) content in early D-shape larvae was significantly higher than that in trochophore and D-shape larvae, while the serotonin (5-HT) content in early D-shape larvae and D-shape larvae was obviously higher than that in trochophore. Incubation of trochophore with 5-HT or DA could accelerate the formation of calcified shell, and the treatments with selective antagonists of receptors for 5-HT and DA (Cg5-HTR-1 and CgD1DR-1) obviously inhibited the formation of calcified shells. When oyster larvae were subjected to an experimental acidified treatment (pH 7.4), the biosynthesis of 5-HT and DA was inhibited, while the mRNA expression levels of the components in TGF-β pathway were significantly up-regulated in D-shape larvae. Moreover, the phosphorylation of TIR and the translocation of smad4 were hindered upon acidification treatments, and the expression patterns of chitinase and tyrosinase were completely reverted. These results collectively suggested that monoamine neurotransmitters 5-HT and DA could modulate the initial shell formation in oyster larvae through TGF-β smad pathway by regulating the expression of tyrosinase and chitinase to guarantee the chitin synthesis for shell formation. CO₂-induced seawater acidification could suppress the biosynthesis of 5-HT and DA, as well as the activation of TGF-β smad pathway, which would subvert the expression patterns of chitinase and tyrosinase and cause the failure of initial shell formation in oyster early D-shape larvae.

Early shell field morphogenesis of a patellogastropod mollusk predominantly relies on cell movement and F-actin dynamics

Background

The morphogenesis of the shell field is an essential step of molluscan shell formation, which exhibits both conserved features and interlineage variations. As one major gastropod lineage, the patellogastropods show different characters in its shell field morphogenesis compared to other gastropods (e.g., the pulmonate gastropod Lymnaea stagnalis), likely related to its epibolic gastrulation. The investigation on the shell field morphogenesis of patellogastropods would be useful to reveal the lineage-specific characters in the process and explore the deep conservation among different molluscan lineages.

Results

We investigated the early shell field morphogenesis in the patellogastropod Lottia goshimai using multiple techniques. Electron microscopy revealed distinct morphological characters for the central and peripheral cells of the characteristic rosette-like shell field. Gene expression analysis and F-actin staining suggested that the shell field morphogenesis in this species predominantly relied on cell movement and F-actin dynamics, while BrdU assay revealed that cell proliferation contributed little to the process. We found constant contacts between ectodermal and meso/endodermal tissues during the early stages of shell field morphogenesis, which did not support the induction of shell field by endodermal tissues in general, but a potential stage-specific induction was indicated.

Conclusions

Our results emphasize the roles of cell movement and F-actin dynamics during the morphogenesis of the shell field in Lo. goshimai, and suggest potential regulators such as diffusible factors and F-actin modulators. These findings reflect the differences in shell field morphogenesis of different gastropods, and add to the knowledge of molluscan larval shell formation.

Establishment of the novel bivalve body plan through modification of early developmental events in mollusks

Mollusks have a wide variety of body plans, which develop through conserved early embryogenesis, namely spiral embryonic development and trochophore larvae. Although the comparative study of mollusks has attracted the interest of evolutionary developmental biology researchers, less attention has been paid to bivalves. In this review, we focused on the evolutionary process from single-shell ancestors to bivalves, which possess bilaterally separated shells. Our study tracing the lineage of shell field cells in bivalves did not support the old hypothesis that shell plate morphology is due to modification of the spiral cleavage pattern. Rather, we suggest that modification of the shell field induction process is the key to understanding the evolution of shell morphology. The novel body plan of bivalves cannot be established solely via separating shell plates, but rather requires the formation of additional organs, such as adductor muscles. The evolutionary biology of bivalves offers a unique view on how multiple organs evolve in a coordinated manner to establish a novel body plan.

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.

PfSMAD1/5 Can Interact with PfSMAD4 to Inhibit PfMSX to Regulate Shell Biomineralization in Pinctada fucata martensii

The BMP2 signal transduced by SMAD1/5 plays an important role in osteoblast differentiation and bone formation. Shell formation of Pinctada fucata martensii is a typical biomineralization process that is similar to that of teeth/bone formation. However, whether the Pinctada fucata BMP2 (PfBMP2) signal transduced by PfSMAD1/5 occurs in P. f. martensii, how the PfBMP2 signal is transduced by PfSMAD1/5, and how PfSMAD1/5 regulates the biomineralization process in this species and other shellfish are poorly understood. Therefore, injection experiments of recombinant PfBMP2 and inhibitor dorsomorphin revealed that PfSMAD1/5 can transduce PfBMP2 signals. Subcellular localization and bimolecular fluorescence complementation assays indicated that PfSMAD1/5 phosphorylated by PfBMPR1b interacts with PfSMAD4 in the cytoplasm to form a complex, which translocates to the nucleus to transduce PfBMP2 signals. Co-immunoprecipitation and luciferase assays revealed that PfSMAD1/5 may interact with PfMSX to dislodge it from its binding element, resulting in initiation of mantle gene transcription. The in vivo functional assay showed that knockdown of PfMSAD1/5 decreased expression of shell matrix genes and disordered the nacreous layer, and the correlation assay of shell regeneration showed the concomitant expression pattern of PfSMAD1/5 and shell matrix genes. Together, these data showed that PfSMAD1/5 can transduce PfBMP2 signals to regulate shell biomineralization in P. f. martensii, which illustrated conservation of the BMP2-SMAD signal pathway among invertebrates. Particularly, the results suggest that there is only one PfMSX gene, which functions like the Hox gene in vertebrates, that interacts with PfSMAD1/5 in a protein-protein action form and plays the role of transcription repressor.

Identification of three cell populations from the shell gland of a bivalve mollusc

The molluscan larval shell formation is a complicated process. There is evidence that the mantle of the primary larva (trochophore) contains functionally different cell populations with distinct gene expression profiles. However, it remains unclear how these cells are specified. In the present study, we identified three cell populations from the shell gland in earlier stages (gastrula) from the bivalve mollusc Crassostrea gigas. These cell populations were determined by analyzing the co-expression relationships among six potential shell formation (pSF) genes using two-color hybridization. The three cell populations, which we designated as SGCPs (shell gland cell populations), formed a concentric-circle pattern from outside to inside of the shell gland. SGCP I was located in the outer edge of the shell gland and the cells expressed pax2/5/8, gata2/3, and bmp2/4. SGCP II was located more internally and the cells expressed two engrailed genes. The last population, SGCP III, was located in the central region of the shell gland and the cells expressed lox4. Determination of the gene expression profiles of SGCPs would help trace their origins and fates and elucidate how these cell populations are specified. Moreover, potential roles of the SGCPs, e.g., development of sensory cells and shell biogenesis, are suggested. Our results reveal the internal organization of the embryonic shell gland at the molecular level and add to the knowledge of larval shell formation.

Growth and morphogenesis of the gastropod shell

Gastropod shell morphologies are famously diverse but generally share a common geometry, the logarithmic coil. Variations on this morphology have been modeled mathematically and computationally but the developmental biology of shell morphogenesis remains poorly understood. Here we characterize the organization and growth patterns of the shell-secreting epithelium of the larval shell of the basket whelk Tritia (also known as Ilyanassa). Despite the larval shell's relative simplicity, we find a surprisingly complex organization of the shell margin in terms of rows and zones of cells. We examined cell division patterns with EdU incorporation assays and found two growth zones within the shell margin. In the more anterior aperture growth zone, we find that inferred division angles are biased to lie parallel to the shell edge, and these divisions occur more on the margin's left side. In the more posterior mantle epithelium growth zone, inferred divisions are significantly biased to the right, relative to the anterior-posterior axis. These growth zones, and the left-right asymmetries in cleavage patterns they display, can explain the major modes of shell morphogenesis at the level of cellular behavior. In a gastropod with a different coiling geometry, Planorbella sp., we find similar shell margin organization and growth zones as Tritia, but different left-right asymmetries than we observed in the helically coiled shell of Tritia These results indicate that differential growth patterns in the mantle edge epithelium contribute to shell shape in gastropod shells and identify cellular mechanisms that may vary to generate shell diversity in evolution.

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.

The evolution of molluscs

Molluscs are extremely diverse invertebrate animals with a rich fossil record, highly divergent life cycles, and considerable economical and ecological importance. Key representatives include worm-like aplacophorans, armoured groups (e.g. polyplacophorans, gastropods, bivalves) and the highly complex cephalopods. Molluscan origins and evolution of their different phenotypes have largely remained unresolved, but significant progress has been made over recent years. Phylogenomic studies revealed a dichotomy of the phylum, resulting in Aculifera (shell-less aplacophorans and multi-shelled polyplacophorans) and Conchifera (all other, primarily uni-shelled groups). This challenged traditional hypotheses that proposed that molluscs gradually evolved complex phenotypes from simple, worm-like animals, a view that is corroborated by developmental studies that showed that aplacophorans are secondarily simplified. Gene expression data indicate that key regulators involved in anterior-posterior patterning (the homeobox-containing Hox genes) lost this function and were co-opted into the evolution of taxon-specific novelties in conchiferans. While the bone morphogenetic protein (BMP)/decapentaplegic (Dpp) signalling pathway, that mediates dorso-ventral axis formation, and molecular components that establish chirality appear to be more conserved between molluscs and other metazoans, variations from the common scheme occur within molluscan sublineages. The deviation of various molluscs from developmental pathways that otherwise appear widely conserved among metazoans provides novel hypotheses on molluscan evolution that can be tested with genome editing tools such as the CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/clustered regularly interspaced short palindromic repeats-associated protein9) system.

GSK3β controls the timing and pattern of the fifth spiral cleavage at the 2-4 cell stage in Lymnaea stagnalis

Establishment of the body plan of multicellular organisms by the primary body axis determination and cell-fate specification is a key issue in biology. We have examined the mRNA localization of three Wnt pathway components gsk3β, β-catenin, and disheveled and investigated the effects of four selective inhibitors of these proteins on the early developmental stages of the spiral cleavage embryo of the fresh water snail Lymnaea (L.) stagnalis. mRNAs for gsk3β and β-catenin were distributed uniformly throughout the embryo during development whereas disheveled mRNA showed specific localization with intra- and inter-blastomere differences in concentration along the A-V axis during spiral cleavages. Remarkably, through inhibitor studies, we identified a short sensitive period from the 2- to 4-cell stage in which GSK3β inhibition by the highly specific 1-azakenpaullone (AZ) and by LiCl induced a subsequent dramatic developmental delay and alteration of the cleavage patterns of blastomeres at the fifth cleavage (16- to 24-cell stage) resulting in exogastrulation and other abnormalities in later stages. Inhibition of β-Catenin or Disheveled had no effect. Our inhibitor experiments establish a novel role for GSK3β in the developmental timing and orientated cell division of the snail embryo. Further work will be needed to identify the downstream targets of the kinase.

Proteinaceous secretion of bioadhesive produced during crawling and settlement of Crassostrea gigas larvae

Bioadhesion of marine organisms has been intensively studied over the last decade because of their ability to attach in various wet environmental conditions and the potential this offers for biotechnology applications. Many marine mollusc species are characterized by a two-phase life history: pelagic larvae settle prior to metamorphosis to a benthic stage. The oyster Crassostrea gigas has been extensively studied for its economic and ecological importance. However, the bioadhesive produced by ready to settle larvae of this species has been little studied. The pediveliger stage of oysters is characterized by the genesis of a specific organ essential for adhesion, the foot. Our scanning electron microscopy and histology analysis revealed that in C. gigas the adhesive is produced by several foot glands. This adhesive is composed of numerous fibres of differing structure, suggesting differences in chemical composition and function. Fourier transformed infrared spectroscopy indicated a mainly proteinaceous composition. Proteomic analysis of footprints was able to identify 42 proteins, among which, one uncharacterized protein was selected on the basis of its pediveliger transcriptome specificity and then located by mRNA in situ hybridization, revealing its potential role during substrate exploration before oyster larva settlement.

Regulation of Extracellular Matrix Synthesis by Shell Extracts from the Marine Bivalve Pecten maximus in Human Articular Chondrocytes- Application for Cartilage Engineering

The shells of the bivalve mollusks are organo-mineral structures predominantly composed of calcium carbonate, but also of a minor organic matrix, a mixture of proteins, glycoproteins, and polysaccharides. These proteins are involved in mineral deposition and, more generally, in the spatial organization of the shell crystallites in well-defined microstructures. In this work, we extracted different organic shell extracts (acid-soluble matrix, acid-insoluble matrix, water-soluble matrix, guanidine HCl/EDTA-extracted matrix, referred as ASM, AIM, WSM, and EDTAM, respectively) from the shell of the scallop Pecten maximus and studied their biological activities on human articular chondrocytes (HACs). We found that these extracts differentially modulate the biological activities of HACs, depending on the type of extraction and the concentration used. Furthermore, we showed that, unlike ASM and AIM, WSM promotes maintenance of the chondrocyte phenotype in monolayer culture. WSM increased the expression of chondrocyte-specific markers (aggrecan and type II collagen), without enhancing that of the main chondrocyte dedifferentiation marker (type I collagen). We also demonstrated that WSM could favor redifferentiation of chondrocyte in collagen sponge scaffold in hypoxia. Thus, this study suggests that the organic matrix of Pecten maximus, particularly WSM, may contain interesting molecules with chondrogenic effects. Our research emphasizes the potential use of WSM of Pecten maximus for cell therapy of cartilage.

Possible co-option of engrailed during brachiopod and mollusc shell development

In molluscs, two homeobox genes, engrailed (en) and distal-less (dlx), are transcription factors that are expressed in correlation with shell development. They are expressed in the regions between shell-forming and non-shell-forming cells, likely defining the boundaries of shell-forming fields. Here we investigate the expression of two transcription factors in the brachiopod Lingula anatina We find that en is expressed in larval mantle lobes, whereas dlx is expressed in larval tentacles. We also demonstrate that the embryonic shell marker mantle peroxidase (mpox) is specifically expressed in mantle lobes. Our results suggest that en and mpox are possibly involved in brachiopod embryonic shell development. We discuss the evolutionary developmental origin of lophotrochozoan biomineralization through independent gene co-option.

A SoxC gene related to larval shell development and co-expression analysis of different shell formation genes in early larvae of oyster

Among the potential larval shell formation genes in mollusks, most are expressed in cells surrounding the shell field during the early phase of shell formation. The only exception (cgi-tyr1) is expressed in the whole larval mantle and thus represents a novel type of expression pattern. This study reports another gene with such an expression pattern. The gene encoded a SoxC homolog of the Pacific oyster Crassostrea gigas and was named cgi-soxc. Whole-mount in situ hybridization revealed that the gene was highly expressed in the whole larval mantle of early larvae. Based on its spatiotemporal expression, cgi-soxc is hypothesized to be involved in periostracum biogenesis, biomineralization, and regulation of cell proliferation. Furthermore, we investigated the interrelationship between cgi-soxc expression and two additional potential shell formation genes, cgi-tyr1 and cgi-gata2/3. The results confirmed co-expression of the three genes in the larval mantle of early D-veliger. Nevertheless, cgi-gata2/3 was only expressed in the mantle edge, and the other two genes were expressed in all mantle cells. Based on the spatial expression patterns of the three genes, two cell groups were identified from the larval mantle (tyr1 ⁺/soxc ⁺/gata2/3 ⁺ cells and tyr1 ⁺/soxc ⁺/gata2/3 ^- cells) and are important to study the differentiation and function of this tissue. The results of this study enrich our knowledge on the structure and function of larval mantle and provide important information to understand the molecular mechanisms of larval shell formation.

Expression patterns indicate that BMP2/4 and Chordin, not BMP5-8 and Gremlin, mediate dorsal-ventral patterning in the mollusk Crassostrea gigas

Though several bilaterian animals use a conserved BMP2/4-Chordin antagonism to pattern the dorsal-ventral (DV) axis, the only lophotrochozoan species in which early DV patterning has been studied to date, the leech Helobdella robusta, appears to employ BMP5-8 and Gremlin. These findings call into question the conservation of a common DV patterning mechanism among bilaterian animals. To explore whether the unusual DV patterning mechanism in H. robusta is also used in other lophotrochozoan species, we investigated the expression of orthologous genes in the early embryo of a bivalve mollusk, Crassostrea gigas. Searching of the genome and phylogenetic analysis revealed that C. gigas possesses single orthologs of BMP2/4, Chordin, and BMP5-8 and no Gremlin homolog. Whole mount in situ hybridization revealed mRNA localization of BMP2/4 and Chordin on the opposite sides of embryos, suggesting the potential involvement of a BMP2/4-Chordin antagonism in DV patterning in this species. Furthermore, universal BMP5-8 expression and the absence of a Gremlin homolog in the C. gigas genome called into question any major contribution by BMP5-8 and Gremlin to early DV patterning in this species. Additionally, we identified seven genes showing asymmetric expression along the DV axis, providing further insight into DV patterning in C. gigas. We present the first report of a Chordin gene in a lophotrochozoan species and of the opposite expression of BMP2/4 (dorsal) and Chordin (ventral) along the D/V axis of a lophotrochozoan embryo. The findings of this study further the knowledge of axis formation in lophotrochozoan species and provide insight into the evolution of the animal DV patterning mechanism.

Diaphanous gene mutation affects spiral cleavage and chirality in snails

L-R (left and right) symmetry breaking during embryogenesis and the establishment of asymmetric body plan are key issues in developmental biology, but the onset including the handedness-determining gene locus still remains unknown. Using pure dextral (DD) and sinistral (dd) strains of the pond snail Lymnaea stagnalis as well as its F₂ through to F₁₀ backcrossed lines, the single handedness-determining-gene locus was mapped by genetic linkage analysis, BAC cloning and chromosome walking. We have identified the actin-related diaphanous gene Lsdia1 as the strongest candidate. Although the cDNA and derived amino acid sequences of the tandemly duplicated Lsdia1 and Lsdia2 genes are very similar, we could discriminate the two genes/proteins in our molecular biology experiments. The Lsdia1 gene of the sinistral strain carries a frameshift mutation that abrogates full-length LsDia1 protein expression. In the dextral strain, it is already translated prior to oviposition. Expression of Lsdia1 (only in the dextral strain) and Lsdia2 (in both chirality) decreases after the 1-cell stage, with no asymmetric localization throughout. The evolutionary relationships among body handedness, SD/SI (spiral deformation/spindle inclination) at the third cleavage, and expression of diaphanous proteins are discussed in comparison with three other pond snails (L. peregra, Physa acuta and Indoplanorbis exustus).

Assessment of housekeeping genes as internal references in quantitative expression analysis during early development of oyster

The early development of mollusks exhibits important characteristics from the developmental and evolutionary perspective. With the increasing number of genome-wide studies, accurate analyses of quantitative gene expression during development are impeded by the lack of validated reference genes. To improve the situation, in this study, we analyzed the expression stability of seven candidate housekeeping genes during early development of the Pacific oyster Crassostrea gigas: actin, glyceraldehyde-3-phosphate dehydrogenase (gapdh), α subunit of elongation factor 1 (elf1α), adp-ribosylation factor 1 (arf1), heterogeneous nuclear ribonucleoprotein q, ubiquitin-conjugating enzyme e2d2 and ribosomal protein s18. We focused on 11 stages from oocyte to D-veliger, which include crucial developmental processes such as axis determination, gastrulation and shell formation. Gene expression stabilities were assessed with the three commonly used programs geNorm, NormFinder and BestKeeper. Although the results obtained with the three programs varied to some extent, in general, arf1, elf1α and gapdh were highly ranked and actin was poorly ranked. This analysis also indicated that multiple genes should be used for normalization, and we concluded that arf1-elf1α-gapdh should be used as internal references. The findings of this study will help researchers to obtain accurate results in future quantitative gene expression analysis of development in bivalve mollusks.

Sea shell diversity and rapidly evolving secretomes: insights into the evolution of biomineralization

An external skeleton is an essential part of the body plan of many animals and is thought to be one of the key factors that enabled the great expansion in animal diversity and disparity during the Cambrian explosion. Molluscs are considered ideal to study the evolution of biomineralization because of their diversity of highly complex, robust and patterned shells. The molluscan shell forms externally at the interface of animal and environment, and involves controlled deposition of calcium carbonate within a framework of macromolecules that are secreted from the dorsal mantle epithelium. Despite its deep conservation within Mollusca, the mantle is capable of producing an incredible diversity of shell patterns, and macro- and micro-architectures. Here we review recent developments within the field of molluscan biomineralization, focusing on the genes expressed in the mantle that encode secreted proteins. The so-called mantle secretome appears to regulate shell deposition and patterning and in some cases becomes part of the shell matrix. Recent transcriptomic and proteomic studies have revealed marked differences in the mantle secretomes of even closely-related molluscs; these typically exceed expected differences based on characteristics of the external shell. All mantle secretomes surveyed to date include novel genes encoding lineage-restricted proteins and unique combinations of co-opted ancient genes. A surprisingly large proportion of both ancient and novel secreted proteins containing simple repetitive motifs or domains that are often modular in construction. These repetitive low complexity domains (RLCDs) appear to further promote the evolvability of the mantle secretome, resulting in domain shuffling, expansion and loss. RLCD families further evolve via slippage and other mechanisms associated with repetitive sequences. As analogous types of secreted proteins are expressed in biomineralizing tissues in other animals, insights into the evolution of the genes underlying molluscan shell formation may be applied more broadly to understanding the evolution of metazoan biomineralization.

PfSMAD4 plays a role in biomineralization and can transduce bone morphogenetic protein-2 signals in the pearl oyster Pinctada fucata

Background

Mollusca is the second largest phylum in nature. The shell of molluscs is a remarkable example of a natural composite biomaterial. Biomineralization and how it affects mollusks is a popular research topic. The BMP-2 signaling pathway plays a canonical role in biomineralization. SMAD4 is an intracellular transmitter in the BMP signaling pathway in mammals, and some genomic data show SMAD4’s involvment in BMP signaling in invertbrates, but whether SMAD4 plays a conservative role in pearl oyster, Pinctada fucata, still need to be tested.

Results

In this study, we identified a SMAD4 gene (hereafter designated PfSMAD4) in pearl oyster Pinctada fucata. Bioinformatics analysis of PfSMAD4 showed high identity with its orthologs. PfSMAD4 was located in the cytoplasm in immunofluorescence assays and analyses of PfSMAD4 mRNA in tissues and developmental stages showed high expression in ovaries and D-shaped larvae. An RNA interference experiment, performed by PfSMAD4 double-stranded RNA (dsRNA) injection, demonstrated inhibition not only of nacre growth but also organic sheet formation with a decrease in PfSMAD4 expression. A knockdown experiment using PfBMP2 dsRNA showed decreased PfBMP2 and PfSMAD4 mRNA and irregular crystallization of the nacreous layer using scanning electron microscopy. In co-transfection experiments, PfBMP2-transactivated reporter constructs contained PfSMAD4 promoter sequences.

Conclusions

Our results suggest that PfSMAD4 plays a role in biomineralization and can transduce BMP signals in P. fucata. Our data provides important clues about the molecular mechanisms that regulate biomineralization in pearl oyster.

Molecular characterization of the BMP7 gene and its potential role in shell formation in Pinctada martensii

Bone morphogenetic protein 7 (BMP7), also called osteogenetic protein-1, can induce bone formation. In this study, the obtained full-length cDNA of BMP7 from Pinctada martensii (Pm-BMP7) was 2972 bp, including a 5'-untranslated region (UTR) of 294 bp, an open reading fragment of 1290 bp encoding a 429 amino acid polypeptide and a 3'-UTR of 1388 bp. The deduced protein sequence of Pm-BMP7 contained a signal peptide, a pro-domain and a mature peptide. The mature peptide consisted of 135 amino acids and included a transforming growth factor β family domain with six shared cysteine residues. The protein sequence of Pm-BMP7 showed 66% identity with that from Crassostrea gigas. Two unigenes encoding Pm-BMPRI (Pm-BMP receptor I) and Pm-BMPRII were obtained from the transcriptome database of P. martensii. Tissue expression analysis demonstrated Pm-BMP7 and Pm-BMPRI were highly expressed in the mantle (shell formation related-tissue), while Pm-BMPRII was highly expressed in the foot. After inhibiting Pm-BMP7 expression using RNA interference (RNAi) technology, Pm-BMP7 mRNA was significantly down-regulated (p < 0.05) in the mantle pallium (nacre formation related-tissue) and the mantle edge (prismatic layer formation related-tissue). The microstructure, observed using a scanning electron microscope, indicated a disordered growth status in the nacre and obvious holes in the prismatic layer in the dsRNA-Pm-BMP7 injected-group. These results suggest that Pm-BMP7 plays a crucial role in the nacre and prismatic layer formation process of the shell.

Quantification and geometric analysis of coiling patterns in gastropod shells based on 3D and 2D image data

The morphology of gastropod shells has been a focus of analyses in ecology and evolution. It has recently emerged as an important issue in developmental biology, thanks to recent advancements in molecular biological techniques. The growing tube model is a theoretical morphological model for describing various coiling patterns of molluscan shells, and it is a useful theoretical tool to relate local tissue growth with global shell morphology. However, the growing tube model has rarely been adopted in empirical research owing to the difficulty in estimating the parameters of the model from morphological data. In this article, I solve this problem by developing methods of parameter estimation when (1) 3D Computed Tomography (CT) data are available and (2) only 2D image data (such as photographs) are available. When 3D CT data are available, the parameters can be estimated by fitting an analytical solution of the growing tube model to the data. When only 2D image data are available, we first fit Raup׳s model to the 2D image data and then convert the parameters of Raup׳s model to those of the growing tube model. To illustrate the use of these methods, I apply them to data generated by a computer simulation of the model. Both methods work well, except when shells grow without coiling. I also demonstrate the effectiveness of the methods by applying the model to actual 3D CT data and 2D image data of land snails. I conclude that the method proposed in this article can reconstruct the coiling pattern from observed data.

Molecular evolution of calcification genes in morphologically similar but phylogenetically unrelated scleractinian corals

Molecular phylogenies of scleractinian corals often fail to agree with traditional phylogenies derived from morphological characters. These discrepancies are generally attributed to non-homologous or morphologically plastic characters used in taxonomic descriptions. Consequently, morphological convergence of coral skeletons among phylogenetically unrelated groups is considered to be the major evolutionary process confounding molecular and morphological hypotheses. A strategy that may help identify cases of convergence and/or diversification in coral morphology is to compare phylogenies of existing "neutral" genetic markers used to estimate genealogic phylogenetic history with phylogenies generated from non-neutral genes involved in calcification (biomineralization). We tested the hypothesis that differences among calcification gene phylogenies with respect to the "neutral" trees may represent convergent or divergent functional strategies among calcification gene proteins that may correlate to aspects of coral skeletal morphology. Partial sequences of two nuclear genes previously determined to be involved in the calcification process in corals, "Cnidaria-III" membrane-bound/secreted α-carbonic anhydrase (CIII-MBSα-CA) and bone morphogenic protein (BMP) 2/4, were PCR-amplified, cloned and sequenced from 31 scleractinian coral species in 26 genera and 9 families. For comparison, "neutral" gene phylogenies were generated from sequences from two protein-coding "non-calcification" genes, one nuclear (β-tubulin) and one mitochondrial (cytochrome b), from the same individuals. Cloned CIII-MBSα-CA sequences were found to be non-neutral, and phylogenetic analyses revealed CIII-MBSα-CAs to exhibit a complex evolutionary history with clones distributed between at least 2 putative gene copies. However, for several coral taxa only one gene copy was recovered. With CIII-MBSα-CA, several recovered clades grouped taxa that differed from the "non-calcification" loci. In some cases, these taxa shared aspects of their skeletal morphology (i.e., convergence or diversification relative to the "non-calcification" loci), but in other cases they did not. For example, the "non-calcification" loci recovered Atlantic and Pacific mussids as separate evolutionary lineages, whereas with CIII-MBSα-CA, clones of two species of Atlantic mussids (Isophyllia sinuosa and Mycetophyllia sp.) and two species of Pacific mussids (Acanthastrea echinata and Lobophyllia hemprichii) were united in a distinct clade (except for one individual of Mycetophyllia). However, this clade also contained other taxa which were not unambiguously correlated with morphological features. BMP2/4 also contained clones that likely represent different gene copies. However, many of the sequences showed no significant deviation from neutrality, and reconstructed phylogenies were similar to the "non-calcification" tree topologies with a few exceptions. Although individual calcification genes are unlikely to precisely explain the diverse morphological features exhibited by scleractinian corals, this study demonstrates an approach for identifying cases where morphological taxonomy may have been misled by convergent and/or divergent molecular evolutionary processes in corals. Studies such as this may help illuminate our understanding of the likely complex evolution of genes involved in the calcification process, and enhance our knowledge of the natural history and biodiversity within this central ecological group.

Transcriptome analysis of the scleractinian coral Stylophora pistillata

The principal architects of coral reefs are the scleractinian corals; these species are divided in two major clades referred to as “robust” and “complex” corals. Although the molecular diversity of the “complex” clade has received considerable attention, with several expressed sequence tag (EST) libraries and a complete genome sequence having been constructed, the “robust” corals have received far less attention, despite the fact that robust corals have been prominent focal points for ecological and physiological studies. Filling this gap affords important opportunities to extend these studies and to improve our understanding of the differences between the two major clades. Here, we present an EST library from Stylophora pistillata (Esper 1797) and systematically analyze the assembled transcripts compared to putative homologs from the complete proteomes of six well-characterized metazoans: Nematostella vectensis, Hydra magnipapillata, Caenorhabditis elegans, Drosophila melanogaster, Strongylocentrotus purpuratus, Ciona intestinalis and Homo sapiens. Furthermore, comparative analyses of the Stylophora pistillata ESTs were performed against several Cnidaria from the Scleractinia, Actiniaria and Hydrozoa, as well as against other stony corals separately. Functional characterization of S. pistillata transcripts into KOG/COG categories and further description of Wnt and bone morphogenetic protein (BMP) signaling pathways showed that the assembled EST library provides sufficient data and coverage. These features of this new library suggest considerable opportunities for extending our understanding of the molecular and physiological behavior of “robust” corals.

A genome-wide survey of genes encoding transcription factors in the Japanese pearl oyster, Pinctada fucata: I. homeobox genes

Homeobox genes are involved in various aspects of the development of multicellular animals, including anterior-posterior patterning of the body plan. We performed a genomic survey of homeobox genes in the Japanese pearl oyster, Pinctada fucata, and annotated 92 homeobox-containing genes and five homeobox-less Pax genes. This species possesses 10 or 11 Hox genes. We annotated another homeobox genes that cover 77 out of the 111 gene families identified in the amphioxus genome. Investigation of these repertoires of homeobox genes will shed new light on the comparatively less well-understood lophotrochozoan development.