Computationally predicted gene regulatory networks in molluscan biomineralization identify extracellular matrix production and ion transportation pathways

Abstract Motivation The molecular processes regulating molluscan shell production remain relatively uncharacterized, despite the clear evolutionary and societal importance of biomineralization. Results Here we built the first computationally predicted gene regulatory network (GRN) for molluscan biomineralization using Antarctic clam (Laternula elliptica) mantle gene expression data produced over an age-categorized shell damage-repair time-course. We used previously published in vivo in situ hybridization expression data to ground truth gene interactions predicted by the GRN and show that candidate biomineralization genes from different shell layers, and hence microstructures, were connected in unique modules. We characterized two biomineralization modules of the GRN and hypothesize that one module is responsible for translating the extracellular proteins required for growing, repairing or remodelling the nacreous shell layer, whereas the second module orchestrates the transport of both ions and proteins to the shell secretion site, which are required during normal shell growth, and repair. Our findings demonstrate that unbiased computational methods are particularly valuable for studying fundamental biological processes and gene interactions in non-model species where rich sources of gene expression data exist, but annotation rates are poor and the ability to carry out true functional tests are still lacking. Availability and implementation The raw RNA-Seq data is freely available for download from NCBI SRA (Accession: PRJNA398984), the assembled and annotated transcriptome can be viewed and downloaded from molluscDB (ensembl.molluscdb.org) and in addition, the assembled transcripts, reconstructed GRN, modules and detailed annotations are all available as Supplementary Files. Supplementary information Supplementary data are available at Bioinformatics online.

[1]  Rongqing Zhang,et al.  Novel Basic Protein, PfN23, Functions as Key Macromolecule during Nacre Formation* , 2012, The Journal of Biological Chemistry.

[2]  Melody S Clark,et al.  Transcriptomic response to shell damage in the Antarctic clam, Laternula elliptica: time scales and spatial localisation. , 2015, Marine genomics.

[3]  J. Henry,et al.  CRISPR/Cas9‐mediated genome modification in the mollusc, Crepidula fornicata , 2015, Genesis.

[4]  Subburaman Mohan,et al.  The Roles and Mechanisms of Actions of Vitamin C in Bone: New Developments , 2015, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[5]  Benjamin Marie,et al.  Different secretory repertoires control the biomineralization processes of prism and nacre deposition of the pearl oyster shell , 2012, Proceedings of the National Academy of Sciences.

[6]  Rachel L. Flores,et al.  The skeletal proteome of the sea star Patiria miniata and evolution of biomineralization in echinoderms , 2017, BMC Evolutionary Biology.

[7]  T Morita,et al.  A carbonic anhydrase from the nacreous layer in oyster pearls. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[8]  Frédéric Marin,et al.  The Role of Enzymes in Biomineralization Processes , 2010 .

[9]  Philip C. J. Donoghue,et al.  Evolving form and function: fossils and development , 2005 .

[10]  Helmut Sigel,et al.  Biomineralization : from nature to application , 2010 .

[11]  Pauline M. Ross,et al.  Impacts of ocean acidification on marine shelled molluscs , 2013, Marine Biology.

[12]  Lloyd S. Peck,et al.  Cellular stress responses to chronic heat shock and shell damage in temperate Mya truncata , 2018, Cell Stress and Chaperones.

[13]  Benjamin Marie,et al.  Insights from the Shell Proteome: Biomineralization to Adaptation , 2016, Molecular biology and evolution.

[14]  Chris Wiggins,et al.  ARACNE: An Algorithm for the Reconstruction of Gene Regulatory Networks in a Mammalian Cellular Context , 2004, BMC Bioinformatics.

[15]  Qi Wang,et al.  Functional characterization and molecular mechanism exploration of three granulin epithelin precursor splice variants in biomineralization of the pearl oyster Pinctada fucata , 2015, Molecular Genetics and Genomics.

[16]  Melody S Clark,et al.  Characterization of the mantle transcriptome in bivalves: Pecten maximus, Mytilus edulis and Crassostrea gigas. , 2016, Marine genomics.

[17]  Benjamin Marie,et al.  Molecular modularity and asymmetry of the molluscan mantle revealed by a gene expression atlas , 2018, GigaScience.

[18]  Brian Biehs,et al.  From molecules to mastication: the development and evolution of teeth , 2013, Wiley interdisciplinary reviews. Developmental biology.

[19]  Reiko Kuroda,et al.  The development of CRISPR for a mollusc establishes the formin Lsdia1 as the long-sought gene for snail dextral/sinistral coiling , 2019, Development.

[20]  Felipe Aguilera,et al.  Evolution of the tyrosinase gene family in bivalve molluscs: independent expansion of the mantle gene repertoire. , 2014, Acta biomaterialia.

[21]  Melody S Clark,et al.  Insights into shell deposition in the Antarctic bivalve Laternula elliptica: gene discovery in the mantle transcriptome using 454 pyrosequencing , 2010, BMC Genomics.

[22]  Benjamin Marie,et al.  An Antarctic molluscan biomineralisation tool-kit , 2016, Scientific Reports.

[23]  Alfonso Valencia,et al.  A hierarchical unsupervised growing neural network for clustering gene expression patterns , 2001, Bioinform..

[24]  Ullrich Steiner,et al.  Biomimetic layer-by-layer assembly of artificial nacre , 2012, Nature Communications.

[25]  Philip Rosenstiel,et al.  A shell regeneration assay to identify biomineralization candidate genes in mytilid mussels. , 2016, Marine genomics.

[26]  K. Boesze-Battaglia,et al.  The role of alkaline phosphatase in mineralization , 2007 .

[27]  Lingming Zhang,et al.  Toll-like receptors in invertebrate innate immunity , 2005 .

[28]  Antje H L Fischer,et al.  Deployment of regulatory genes during gastrulation and germ layer specification in a model spiralian mollusc Crepidula , 2015, Developmental dynamics : an official publication of the American Association of Anatomists.

[29]  Lloyd S. Peck,et al.  Characterisation of the mantle transcriptome and biomineralisation genes in the blunt-gaper clam, Mya truncata. , 2016, Marine genomics.

[30]  T. Rapoport Protein translocation across the eukaryotic endoplasmic reticulum and bacterial plasma membranes , 2007, Nature.

[31]  D. Jackson,et al.  An ancient process in a modern mollusc: early development of the shell in Lymnaea stagnalis , 2013, BMC Developmental Biology.

[32]  Benjamin Marie,et al.  The evolution of metazoan α-carbonic anhydrases and their roles in calcium carbonate biomineralization , 2014, Frontiers in Zoology.

[33]  R. Zimmermann,et al.  Let’s talk about Secs: Sec61, Sec62 and Sec63 in signal transduction, oncology and personalized medicine , 2017, Signal Transduction and Targeted Therapy.

[34]  Michael Banf,et al.  Computational inference of gene regulatory networks: Approaches, limitations and opportunities. , 2017, Biochimica et biophysica acta. Gene regulatory mechanisms.

[35]  Benjamin Marie,et al.  Proteomic analysis of the organic matrix of the abalone Haliotis asinina calcified shell , 2010, Proteome Science.

[36]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[37]  Masayoshi Yamaguchi,et al.  The role of regucalcin in bone homeostasis: involvement as a novel cytokine. , 2014, Integrative biology : quantitative biosciences from nano to macro.

[38]  Takashi Kato,et al.  An Acidic Matrix Protein, Pif, Is a Key Macromolecule for Nacre Formation , 2009, Science.

[39]  Benjamin Marie,et al.  The shell‐forming proteome of Lottia gigantea reveals both deep conservations and lineage‐specific novelties , 2013, The FEBS journal.

[40]  Gary W. Harding,et al.  Skeletal overgrowth and deafness in mice lacking fibroblast growth factor receptor 3 , 1996, Nature Genetics.

[41]  A. Wanninger,et al.  Comparative transcriptomics enlarges the toolkit of known developmental genes in mollusks , 2016, BMC Genomics.

[42]  Amiel A. Dror,et al.  Atrophic thyroid follicles and inner ear defects reminiscent of cochlear hypothyroidism in Slc26a4-related deafness , 2014, Mammalian Genome.

[43]  K. Sundell,et al.  Calcium transfer across the outer mantle epithelium in the Pacific oyster, Crassostrea gigas , 2018, Proceedings of the Royal Society B.

[44]  Ross Crawford,et al.  Notch expressed by osteocytes plays a critical role in mineralisation , 2018, Journal of Molecular Medicine.

[45]  Bo Zhang,et al.  Identification and Differential Expression of Biomineralization Genes in the Mantle of Pearl Oyster Pinctada fucata , 2017, Marine Biotechnology.

[46]  Mei Chen,et al.  Extracellular and Non-Chaperone Function of Heat Shock Protein-90α Is Required for Skin Wound Healing. , 2017, The Journal of investigative dermatology.

[47]  Nadav Ahituv,et al.  Identification of novel Fgf enhancers and their role in dental evolution , 2016, Evolution & development.