The SCPP gene repertoire in bony vertebrates and graded differences in mineralized tissues

Single-cell sequencing provides clues about the developmental genetic basis of evolutionary adaptations in syngnathid fishes

Seahorses, pipefishes, and seadragons are fishes from the family Syngnathidae that have evolved extraordinary traits including male pregnancy, elongated snouts, loss of teeth, and dermal bony armor. The developmental genetic and cellular changes that led to the evolution of these traits are largely unknown. Recent syngnathid genome assemblies revealed suggestive gene content differences and provided the opportunity for detailed genetic analyses. We created a single-cell RNA sequencing atlas of Gulf pipefish embryos to understand the developmental basis of four traits: derived head shape, toothlessness, dermal armor, and male pregnancy. We completed marker gene analyses, built genetic networks, and examined the spatial expression of select genes. We identified osteochondrogenic mesenchymal cells in the elongating face that express regulatory genes bmp4, sfrp1a, and prdm16. We found no evidence for tooth primordia cells, and we observed re-deployment of osteoblast genetic networks in developing dermal armor. Finally, we found that epidermal cells expressed nutrient processing and environmental sensing genes, potentially relevant for the brooding environment. The examined pipefish evolutionary innovations are composed of recognizable cell types, suggesting that derived features originate from changes within existing gene networks. Future work addressing syngnathid gene networks across multiple stages and species is essential for understanding how the novelties of these fish evolved.

Paleo-evo-devo implications of a revised conceptualization of enameloids and enamels

Understanding the origin and evolution of the mineralized skeleton is crucial for unravelling vertebrate history. However, several limitations hamper our progress. The first obstacle is the lack of uniformity and clarity in the literature for the definition of the tissues of concern, especially of enameloid(s) and enamel(s), resulting in ambiguous terminology and inconsistencies among studies. Moreover, the identification criteria currently employed to characterize hypermineralized tissues in extinct taxa, such as the presence or absence of tubules for enameloids, may lead to unsupported conclusions. We suggest that comparative developmental studies may be key to unambiguous terminology, truly diagnostic identification criteria and developmentally informed evolutionary hypotheses. We exemplify this approach by: (i) introducing a new conceptual framework for enameloid(s) and enamel(s), with clear terminologies, definitions and interactions between concepts; (ii) suggesting more rigorous ways to identify tissues, based on the observation of defining or additional properties, as well as on the comparison of developmental scenarios when possible; (iii) constructing a clear phylogenetic framework to discuss their homologies and highlighting possible transitions between these tissues; and by (iv) proposing developmental models that explain both enamel and enameloid formation, and suggest possible transitions between them.

Evolution, development, and regeneration of tooth-like epithelial appendages in sharks

Sharks and their relatives are typically covered in highly specialized epithelial appendages embedded in the skin called dermal denticles; ancient tooth-like units (odontodes) composed of dentine and enamel-like tissues. These 'skin teeth' are remarkably similar to oral teeth of vertebrates and share comparable morphological and genetic signatures. Here we review the histological and morphological data from embryonic sharks to uncover characters that unite all tooth-like elements (odontodes), including teeth and skin denticles in sharks. In addition, we review the differences between the skin and oral odontodes that reflect their varied capacity for renewal. Our observations have begun to decipher the developmental and genetic shifts that separate these seemingly similar dental units, including elements of the regenerative nature in both oral teeth and the emerging skin denticles from the small-spotted catshark (Scyliorhinus canicula) and other chondrichthyan models. Ultimately, we ask what defines a tooth at both the molecular and morphological level. These insights aim to help us understand how nature makes, replaces and evolves a vast array of odontodes.

Growth patterns of caudal fin rays are informed by both external signals from the regenerating organ and remembered identity autonomous to the local tissue

Regenerating tissues must remember or interpret their spatial position, using this information to restore original size and patterning. The external skeleton of the zebrafish caudal fin is composed of 18 rays; after any portion of the fin is amputated, position-dependent regenerative growth restores each ray to its original length. We tested for transcriptional differences during regeneration of proximal versus distal tissues and identified 489 genes that differed in proximodistal expression. Thyroid hormone directs multiple aspects of ray patterning along the proximodistal axis, and we identified 364 transcripts showing a proximodistal expression pattern that was dependent on thyroid hormone context. To test what aspects of ray positional identity are directed by extrinsic environental cues versus remembered identity autonomous to the tissue, we transplanted distal portions of rays to proximal environments and evaluated regeneration within the new location. Native regenerating proximal tissue showed robust expression of scpp7, a transcript with thyroid-regulated proximal enrichment; in contrast, regenerating rays originating from transplanted distal tissue showed reduced (distal-like) expression during outgrowth. These distal-to-proximal transplants regenerated far beyond the length of the graft itself, indicating that cues from the proximal environment promoted additional growth. Nonetheless, these transplants initiated regeneration at a much slower rate compared to controls, suggesting memory of distal identity was retained by the transplanted tissue. This early growth retardation caused rays that originated from transplants to grow noticeably shorter than neighboring native rays. While several aspects of fin ray morphology (bifurcation, segment length) were found to be determined by the environment, we found that both regeneration speed and ray length are remembered autonomously by tissues, and that persist through multiple rounds of amputation and regeneration.

Single Cell Sequencing Provides Clues about the Developmental Genetic Basis of Evolutionary Adaptations in Syngnathid Fishes

Seahorses, pipefishes, and seadragons are fishes from the family Syngnathidae that have evolved extraordinary traits including male pregnancy, elongated snouts, loss of teeth, and dermal bony armor. The developmental genetic and cellular changes that led to the evolution of these traits are largely unknown. Recent syngnathid genome assemblies revealed suggestive gene content differences and provide the opportunity for detailed genetic analyses. We created a single cell RNA sequencing atlas of Gulf pipefish embryos to understand the developmental basis of four traits: derived head shape, toothlessness, dermal armor, and male pregnancy. We completed marker gene analyses, built genetic networks, and examined spatial expression of select genes. We identified osteochondrogenic mesenchymal cells in the elongating face that express regulatory genes bmp4, sfrp1a, and prdm16. We found no evidence for tooth primordia cells, and we observed re-deployment of osteoblast genetic networks in developing dermal armor. Finally, we found that epidermal cells expressed nutrient processing and environmental sensing genes, potentially relevant for the brooding environment. The examined pipefish evolutionary innovations are composed of recognizable cell types, suggesting derived features originate from changes within existing gene networks. Future work addressing syngnathid gene networks across multiple stages and species is essential for understanding how their novelties evolved.

Xenopus tropicalis osteoblast-specific open chromatin regions reveal promoters and enhancers involved in human skeletal phenotypes and shed light on early vertebrate evolution

While understanding the genetic underpinnings of osteogenesis has far-reaching implications for skeletal diseases and evolution, a comprehensive characterization of the osteoblastic regulatory landscape in non-mammalian vertebrates is still lacking. Here, we compared the ATAC-Seq profile of Xenopus tropicalis (Xt) osteoblasts to a variety of non mineralizing control tissues, and identified osteoblast-specific nucleosome free regions (NFRs) at 527 promoters and 6747 distal regions. Sequence analyses, Gene Ontology, RNA-Seq and ChIP-Seq against four key histone marks confirmed that the distal regions correspond to bona fide osteogenic transcriptional enhancers exhibiting a shared regulatory logic with mammals. We report 425 regulatory regions conserved with human and globally associated to skeletogenic genes. Of these, 35 regions have been shown to impact human skeletal phenotypes by GWAS, including one trps1 enhancer and the runx2 promoter, two genes which are respectively involved in trichorhinophalangeal syndrome type I and cleidocranial dysplasia. Intriguingly, 60 osteoblastic NFRs also align to the genome of the elephant shark, a species lacking osteoblasts and bone tissue. To tackle this paradox, we chose to focus on dlx5 because its conserved promoter, known to integrate regulatory inputs during mammalian osteogenesis, harbours an osteoblast-specific NFR in both frog and human. Hence, we show that dlx5 is expressed in Xt and elephant shark odontoblasts, supporting a common cellular and genetic origin of bone and dentine. Taken together, our work (i) unravels the Xt osteogenic regulatory landscape, (ii) illustrates how cross-species comparisons harvest data relevant to human biology and (iii) reveals that a set of genes including bnc2, dlx5, ebf3, mir199a, nfia, runx2 and zfhx4 drove the development of a primitive form of mineralized skeletal tissue deep in the vertebrate lineage.

Growth patterns of caudal fin rays are informed by both external signals from the regenerating organ and remembered identity autonomous to the local tissue

Regenerating tissues must remember or interpret their spatial position, using this information to restore original size and patterning. The external skeleton of the zebrafish caudal fin is composed of 18 rays; after any portion of the fin is amputated, position-dependent regenerative growth restores each ray to its original length. We tested for transcriptional differences during regeneration of proximal versus distal tissues and identified 489 genes that differed in proximodistal expression. Thyroid hormone directs multiple aspects of ray patterning along the proximodistal axis, and we identified 364 transcripts showing a proximodistal expression pattern that was dependent on thyroid hormone context. To test what aspects of ray positional identity are directed by extrinsic cues versus remembered identity autonomous to the tissue itself, we transplanted distal portions of rays to proximal environments and evaluated regeneration within the new location. While neighboring proximal tissue showed robust expression of scpp7, a transcript with thyroid-regulated proximal enrichment, regenerating rays originating from transplanted distal tissue showed reduced (distal-like) expression during outgrowth. These distal-to-proximal transplants regenerated far beyond the length of the graft itself, indicating that cues from the proximal environment promoted additional growth. Nonetheless, these transplants initially regenerated at a much slower rate compared to controls, suggesting memory of distal identity was retained by the transplanted tissue. This early growth retardation caused rays that originated from transplants to become noticeably shorter than their native neighboring rays. While several aspects of fin ray morphology (bifurcation, segment length) were found to be determined by the environment, regeneration speed and ray length are remembered autonomously by tissues, persisting across multiple rounds of amputation and regeneration.

Genetic regulation of injury-induced heterotopic ossification in adult zebrafish

Heterotopic ossification is the inappropriate formation of bone in soft tissues of the body. It can manifest spontaneously in rare genetic conditions or as a response to injury, known as acquired heterotopic ossification. There are several experimental models for studying acquired heterotopic ossification from different sources of damage. However, their tenuous mechanistic relevance to the human condition, invasive and laborious nature and/or lack of amenability to chemical and genetic screens, limit their utility. To address these limitations, we developed a simple zebrafish injury model that manifests heterotopic ossification with high penetrance in response to clinically emulating injuries, as observed in human myositis ossificans traumatica. Using this model, we defined the transcriptional response to trauma, identifying differentially regulated genes. Mutant analyses revealed that an increase in the activity of the potassium channel Kcnk5b potentiates injury response, whereas loss of function of the interleukin 11 receptor paralogue (Il11ra) resulted in a drastically reduced ossification response. Based on these findings, we postulate that enhanced ionic signalling, specifically through Kcnk5b, regulates the intensity of the skeletogenic injury response, which, in part, requires immune response regulated by Il11ra.

The significant role of glycosaminoglycans in tooth development

This review delves into the roles of glycosaminoglycans (GAGs), integral components of proteoglycans, in tooth development. Proteoglycans consist of a core protein linked to GAG chains, comprised of repeating disaccharide units. GAGs are classified into several types, such as hyaluronic acid, heparan sulfate, chondroitin sulfate, dermatan sulfate, and keratan sulfate. Functioning as critical macromolecular components within the dental basement membrane, these GAGs facilitate cell adhesion and aggregation, and play key roles in regulating cell proliferation and differentiation, thereby significantly influencing tooth morphogenesis. Notably, our recent research has identified the hyaluronan-degrading enzyme Transmembrane protein 2 (Tmem2) and we have conducted functional analyses using mouse models. These studies have unveiled the essential role of Tmem2-mediated hyaluronan degradation and its involvement in hyaluronan-mediated cell adhesion during tooth formation. This review provides a comprehensive summary of the current understanding of GAG functions in tooth development, integrating insights from recent research, and discusses future directions in this field.

Focused ion beam-SEM 3D study of osteodentin in the teeth of the Atlantic wolfish Anarhichas lupus

The palette of mineralized tissues in fish is wide, and this is particularly apparent in fish dentin. While the teeth of all vertebrates except fish contain a single dentinal tissue type, called orthodentin, dentin in the teeth of fish can be one of several different tissue types. The most common dentin type in fish is orthodentin. Orthodentin is characterized by several key structural features that are fundamentally different from those of bone and from those of osteodentin. Osteodentin, the second-most common dentin type in fish (based on the tiny fraction of fish species out of ∼30,000 extant fish species in which tooth structure was so far studied), is found in most Selachians (sharks and rays) as well as in several teleost species, and is structurally different from orthodentin. Here we examine the hypothesis that osteodentin is similar to anosteocytic bone tissue in terms of its micro- and nano-structure. We use Focused Ion Beam-Scanning Electron Microscopy (FIB/SEM), as well as several other high-resolution imaging techniques, to characterize the 3D architecture of the three main components of osteodentin (denteons, inter-denteonal matrix, and the transition zone between them). We show that the matrix of osteodentin, although acellular, is extremely similar to mammalian osteonal bone matrix, both in general morphology and in the three-dimensional nano-arrangement of its mineralized collagen fibrils. We also document the presence of a complex network of nano-channels, similar to such networks recently described in bone. Finally, we document the presence of strings of hyper-mineralized small 'pearls' which surround the denteonal canals, and characterize their structure.

Emergence of saliva protein genes in the secretory calcium-binding phosphoprotein (SCPP) locus and accelerated evolution in primates

Genes within the secretory calcium-binding phosphoprotein (SCPP) family evolved in conjunction with major evolutionary milestones: the formation of a calcified skeleton in vertebrates, the emergence of tooth enamel in fish, and the introduction of lactation in mammals. The SCPP gene family also contains genes expressed primarily and abundantly in human saliva. Here, we explored the evolution of the saliva-related SCPP genes by harnessing currently available genomic and transcriptomic resources. Our findings provide insights into the expansion and diversification of SCPP genes, notably identifying previously undocumented convergent gene duplications. In primate genomes, we found additional duplication and diversification events that affected genes coding for proteins secreted in saliva. These saliva-related SCPP genes exhibit signatures of positive selection in the primate lineage while the other genes in the same locus remain conserved. We found that regulatory shifts and gene turnover events facilitated the accelerated gain of salivary expression. Collectively, our results position the SCPP gene family as a hotbed of evolutionary innovation, suggesting the potential role of dietary and pathogenic pressures in the adaptive diversification of the saliva composition in primates, including humans.

Expression of 20 SCPP genes during tooth and bone mineralization in Senegal bichir

The African bichir (Polypterus senegalus) is a living representative of Polypteriformes. P. senegalus possesses teeth composed of dentin covered by an enameloid cap and a layer of collar enamel on the tooth shaft, as in lepisosteids. A thin layer of enamel matrix can also be found covering the cap enameloid after its maturation and during the collar enamel formation. Teleosts fish do not possess enamel; teeth are protected by cap and collar enameloid, and inversely in sarcopterygians, where teeth are only covered by enamel, with the exception of the cap enameloid in teeth of larval urodeles. The presence of enameloid and enamel in the teeth of the same organism is an opportunity to solve the evolutionary history of the presence of enamel/enameloid in basal actinopterygians. In silico analyses of the jaw transcriptome of a juvenile bichir provided twenty SCPP transcripts. They included enamel, dentin, and bone-specific SCPPs known in sarcopterygians and several actinopterygian-specific SCPPs. The expression of these 20 genes was investigated by in situ hybridizations on jaw sections during tooth and dentary bone formation. A spatiotemporal expression patterns were established and compared with previous studies of SCPP gene expression during enamel/enameloid and bone formation. Similarities and differences were highlighted, and several SCPP transcripts were found specifically expressed during tooth or bone formation suggesting either conserved or new functions of these SCPPs.

Expression of secretory calcium-binding phosphoprotein (scpp) genes in medaka during the formation and replacement of pharyngeal teeth

BACKGROUND

Analyses of tooth families and tooth-forming units in medaka with regard to tooth replacement cycles and the localization of odontogenic stem cell niches in the pharyngeal dentition clearly indicate that continuous tooth replacement is maintained. The secretory calcium-binding phosphoprotein (scpp) gene cluster is involved in the formation of mineralized tissues, such as dental and bone tissues, and the genes encoding multiple SCPPs are conserved in fish, amphibians, reptiles, and mammals. In the present study, we examined the expression patterns of several scpp genes in the pharyngeal teeth of medaka to elucidate their roles during tooth formation and replacement.

METHODS

Himedaka (Japanese medaka, Oryzias latipes) of both sexes (body length: 28 to 33 mm) were used in this study. Real-time quantitative reverse transcription-polymerase chain reaction (PCR) (qPCR) data were evaluated using one-way analysis of variance for multi-group comparisons, and the significance of differences was determined by Tukey's comparison test. The expression of scpp genes was examined using in situ hybridization (ISH) with a digoxigenin-labeled, single-stranded antisense probe.

RESULTS

qPCR results showed that several scpp genes were strongly expressed in pharyngeal tissues. ISH analysis revealed specific expression of scpp1, scpp5, and sparc in tooth germ, and scpp5 was continually expressed in the odontoblasts of teeth attached to pedicles, but not in the osteoblasts of pedicles. In addition, many scpp genes were expressed in inner dental epithelium (ide), but not in odontoblasts, and scpp2 consistently showed epithelial-specific expression in the functional teeth. Taken together, these data indicate that specific expression of scpp2 and scpp5 may play a critical role in pharyngeal tooth formation in medaka.

CONCLUSION

We characterized changes in the expression patterns of scpp genes in medaka during the formation and replacement of pharyngeal teeth.

Unravelling the genomic features, phylogeny and genetic basis of tooth ontogenesis in Characiformes through analysis of four genomes

Characiformes is a diverse and evolutionarily significant order of freshwater fish encompassing over 2,300 species. Despite its diversity, our understanding of Characiformes' evolutionary relationships and adaptive mechanisms is limited due to insufficient genome sequences. In this study, we sequenced and assembled the genomes of four Characiformes species, three of which were chromosome-level assemblies. Our analyses revealed dynamic changes in gene family evolution, repeat sequences and variations in chromosomal collinearity within these genomes. With the assembled genomes, we were not only able to elucidate the evolutionary relationship of the four main orders in Otophysi but also indicated Characiformes as the paraphyletic group. Comparative genomic analysis with other available fish genomes shed light on the evolution of genes related to tooth development in Characiformes. Notably, variations in the copy number of secretory calcium-binding phosphoproteins (SCPP) genes were observed among different orders of Otophysi, indicating their potential contribution to the diversity of tooth types. Our study offers invaluable genome sequences and novel insights into Characiformes' evolution, paving the way for further genomic and evolutionary research in fish.

Transcriptomic profiling of tissue environments critical for post-embryonic patterning and morphogenesis of zebrafish skin

Pigment patterns and skin appendages are prominent features of vertebrate skin. In zebrafish, regularly patterned pigment stripes and an array of calcified scales form simultaneously in the skin during post-embryonic development. Understanding the mechanisms that regulate stripe patterning and scale morphogenesis may lead to the discovery of fundamental mechanisms that govern the development of animal form. To learn about cell types and signaling interactions that govern skin patterning and morphogenesis, we generated and analyzed single-cell transcriptomes of skin from wild-type fish as well as fish having genetic or transgenically induced defects in squamation or pigmentation. These data reveal a previously undescribed population of epidermal cells that express transcripts encoding enamel matrix proteins, suggest hormonal control of epithelial-mesenchymal signaling, clarify the signaling network that governs scale papillae development, and identify a critical role for the hypodermis in supporting pigment cell development. Additionally, these comprehensive single-cell transcriptomic data representing skin phenotypes of biomedical relevance should provide a useful resource for accelerating the discovery of mechanisms that govern skin development and homeostasis.

Amyloid fibril formation by αS1- and β-casein implies that fibril formation is a general property of casein proteins

Caseins are a diverse family of intrinsically disordered proteins present in the milks of all mammals. A property common to two cow paralogues, α_S2- and κ-casein, is their propensity in vitro to form amyloid fibrils, the highly ordered protein aggregates associated with many age-related, including neurological, diseases. In this study, we explored whether amyloid fibril-forming propensity is a general feature of casein proteins by examining the other cow caseins (α_S1 and β) as well as β-caseins from camel and goat. Small-angle X-ray scattering measurements indicated that cow α_S1- and β-casein formed large spherical aggregates at neutral pH and 20°C. Upon incubation at 65°C, α_S1- and β-casein underwent conversion to amyloid fibrils over the course of ten days, as shown by thioflavin T binding, transmission electron microscopy, and X-ray fibre diffraction. At the lower temperature of 37°C where fibril formation was more limited, camel β-casein exhibited a greater fibril-forming propensity than its cow or goat orthologues. Limited proteolysis of cow and camel β-casein fibrils and analysis by mass spectrometry indicated a common amyloidogenic sequence in the proline, glutamine-rich, C-terminal region of β-casein. These findings highlight the persistence of amyloidogenic sequences within caseins, which likely contribute to their functional, heterotypic self-assembly; in all mammalian milks, at least two caseins coalesce to form casein micelles, implying that caseins diversified partly to avoid dysfunctional amyloid fibril formation.

The versatile roles of odontogenic ameloblast-associated protein in odontogenesis, junctional epithelium regeneration and periodontal disease

Junctional epithelium (JE) is a vital epithelial component which forms an attachment to the tooth surface at the gingival sulcus by the adhesion of protein complexes from its basal layer. Disruption of the JE is associated with the development of gingivitis, periodontal disease, and alveolar bone loss. Odontogenic ameloblast-associated (ODAM) is comprised of a signal peptide and an ODAM protein with 12 putative glycosylation sites. It is expressed during odontogenesis by maturation stage ameloblasts and is incorporated into the enamel matrix during the formation of outer and surface layer enamel. ODAM, as a secreted protein which is accumulated at the interface between basal lamina and enamel, mediates the adhesion of the JE to the tooth surface; and is involved with extracellular signalling of WNT and ARHGEF5-RhoA, as well as intracellular signalling of BMP-2-BMPR-IB-ODAM. ODAM is also found to be highly expressed in salivary glands and appears to have implications for the regulation of formation, repair, and regeneration of the JE. Bioinformatics and research data have identified the anti-cancer properties of ODAM, indicating its potential both as a prognostic biomarker and therapeutic target. Understanding the biology of ODAM will help to design therapeutic strategies for periodontal and dental disorders.

Spatiotemporal characteristics of the pharyngeal teeth in interspecific distant hybrids of cyprinid fish: Phylogeny and expression of the initiation marker genes

As an important feeding organ and taxonomical characteristic, the pharyngeal teeth of cyprinid fish have very high morphological diversity and exhibit species-specific numbers and arrangements. Many genes have been verified to regulate the pharyngeal teeth development and act as the initiation marker for teeth. Six initiation marker genes for pharyngeal teeth were used as RNA probes to investigate the expression pattern, and these genes were further used to construct a phylogenetic tree for cyprinid fish including some distant hybrids. The results from in situ hybridization showed that similarities and differences existed in the expression of dlx2b, dlx4b, dlx5a, pitx2, fth1b, and scpp5 in the pharyngeal region of the hybrids (BT) by the crosses of blunt snout bream (BSB, ♀) × topmouth culter (TC, ♂). Particularly, we found a high specificity marker gene scpp5 for the early development of pharyngeal teeth. The Scpp5 expression pattern established a clear graphic representation on the spatiotemporal characteristics of the early morphogenesis of pharyngeal teeth in BT and BSB. Our results suggested that the scpp5 expression in 4V¹, 3V¹, and 5V¹ in BT occurred earlier than that in BSB, while the replacement rate of pharyngeal teeth (4V², 3V², and 5V²) was faster in BSB. Phylogenetic analysis revealed that the six marker genes were highly conserved and could be used as the molecular marker for identifying the parents of the distant hybrids in cyprinid fish. The expression patterns of the scpp5 gene was examined in various tissues, including the brain, gill, heart, liver, muscle, skin, fins, gonad, eye, and kidney, showing that the scpp5 gene was ubiquitously expressed, indicating its important role in cyprinid fish.

Evolution and developmental expression of the sodium-iodide symporter (NIS, slc5a5) gene family: Implications for perchlorate toxicology

The vertebrate sodium-iodide symporter (NIS or SLC5A5) transports iodide into the thyroid follicular cells that synthesize thyroid hormone. The SLC5A protein family includes transporters of vitamins, minerals, and nutrients. Disruption of SLC5A5 function by perchlorate, a pervasive environmental contaminant, leads to human pathologies, especially hypothyroidism. Perchlorate also disrupts the sexual development of model animals, including threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio), but the mechanism of action is unknown. To test the hypothesis that SLC5A5 paralogs are expressed in tissues necessary for the development of reproductive organs, and therefore are plausible candidates to mediate the effects of perchlorate on sexual development, we first investigated the evolutionary history of Slc5a paralogs to better understand potential functional trajectories of the gene family. We identified two clades of slc5a paralogs with respect to an outgroup of sodium/choline cotransporters (slc5a7); these clades are the NIS clade of sodium/iodide and lactate cotransporters (slc5a5, slc5a6, slc5a8, slc5a8, and slc5a12) and the SGLT clade of sodium/glucose cotransporters (slc5a1, slc5a2, slc5a3, slc5a4, slc5a10, and slc5a11). We also characterized expression patterns of slc5a genes during development. Stickleback embryos and early larvae expressed NIS clade genes in connective tissue, cartilage, teeth, and thyroid. Stickleback males and females expressed slc5a5 and its paralogs in gonads. Single-cell transcriptomics (scRNA-seq) on zebrafish sex-genotyped gonads revealed that NIS clade-expressing cells included germ cells (slc5a5, slc5a6a, and slc5a6b) and gonadal soma cells (slc5a8l). These results are consistent with the hypothesis that perchlorate exerts its effects on sexual development by interacting with slc5a5 or its paralogs in reproductive tissues. These findings show novel expression domains of slc5 genes in stickleback and zebrafish, which suggest similar functions across vertebrates including humans, and provide candidates to mediate the effects of perchlorate on sexual development.

Parallel evolution of ameloblastic scpp genes in bony and cartilaginous vertebrates

In bony vertebrates, skeletal mineralization relies on the secretory calcium-binding phosphoproteins (Scpp) family whose members are acidic extracellular proteins posttranslationally regulated by the Fam20°C kinase. As scpp genes are absent from the elephant shark genome, they are currently thought to be specific to bony fishes (osteichthyans). Here, we report a scpp gene present in elasmobranchs (sharks and rays) that evolved from local tandem duplication of sparc-L 5′ exons and show that both genes experienced recent gene conversion in sharks. The elasmobranch scpp is remarkably similar to the osteichthyan scpp members as they share syntenic and gene structure features, code for a conserved signal peptide, tyrosine-rich and aspartate/glutamate-rich regions, and harbor putative Fam20°C phosphorylation sites. In addition, the catshark scpp is coexpressed with sparc-L and fam20°C in tooth and scale ameloblasts, similarly to some osteichthyan scpp genes. Despite these strong similarities, molecular clock and phylogenetic data demonstrate that the elasmobranch scpp gene originated independently from the osteichthyan scpp gene family. Our study reveals convergent events at the sparc-L locus in the two sister clades of jawed vertebrates, leading to parallel diversification of the skeletal biomineralization toolkit. The molecular evolution of sparc-L and its coexpression with fam20°C in catshark ameloblasts provides a unifying genetic basis that suggests that all convergent scpp duplicates inherited similar features from their sparc-L precursor. This conclusion supports a single origin for the hypermineralized outer odontode layer as produced by an ancestral developmental process performed by Sparc-L, implying the homology of the enamel and enameloid tissues in all vertebrates.

Scale Development-Related Genes Identified by Transcriptome Analysis

Scales, as key structures of fish skin, play an important role in physiological function. The study of fish scale development mechanisms provides a basis for exploring the molecular-level developmental differences between scaled and non-scaled fishes. In this study, alizarin red staining was used to divide the different stages of zebrafish (Danio rerio) scale development. Four developmental stages, namely stage I (~17 dpf, scales have not started to grow), stage II (~33 dpf, the point at which scales start to grow), stage III (~41 dpf, the period in which the scales almost cover the whole body), and stage IV (~3 mpf, scales cover the whole body), were determined and used for subsequent transcriptome analysis. WGCNA (weighted correlation network analysis) and DEG (differentially expressed gene) analysis were used for screening the key genes. Based on the comparison between stage II and stage I, 54 hub-genes were identified by WGCNA analysis. Key genes including the Scpp family (Scpp7, Scpp6, Scpp5, and Scpp8), the Fgf family (Fgfr1b and Fgfr3), Tcf7, Wnt10b, Runx2b, and Il2rb were identified by DEG analysis, which indicated that these genes played important roles in the key nodes of scale development signal pathways. Combined with this analysis, the TGF-β, Wnt/β-catenin, and FGF signaling pathways were suggested to be the most important signal pathways for scales starting to grow. This study laid a foundation for exploring the scale development mechanism of other fishes. The scale development candidate genes identified in the current study will facilitate functional gene identifications in the future.

Convergent losses of SCPP genes and ganoid scales among non-teleost actinopterygians

Various secretory calcium-binding phosphoprotein (SCPP) genes are expressed in the skin and jaw during the formation of bone, teeth, and scales in osteichthyans (bony vertebrates). Among these mineralized skeletal units is the ganoid scale, found in many fossil actinopterygians (ray-finned fish) but confirmed only in Polypteriformes (bichirs, reedfish) and Lepisosteiformes (gars) among extant clades. Here, we examined SCPP genes in the genome of seven non-teleost actinopterygian species that possess or do not possess ganoid scales. As a result, 39-43 SCPP genes were identified in Polypteriformes and Lepisosteiformes, whereas 22-24 SCPP genes were found in Acipenseriformes (sturgeons, paddlefish) and Amiiformes (bowfin). Most of these genes form two clusters in the genome of Polypteriformes, Lepisosteiformes, and Amiiformes, and these two clusters are duplicated in Acipenseriformes. Despite their distant phylogenetic relationship, Polypteriformes and Lepisosteiformes retain many orthologous SCPP genes. These results imply that common ancestors of extant actinopterygians possessed a large repertoire of SCPP genes, and that many SCPP genes were lost independently in Acipenseriformes and Amiiformes. Notably, most SCPP genes originally located in one of the two SCPP gene clusters are retained in Polypteriformes and Lepisosteiformes but were secondarily lost in Acipenseriformes and Amiiformes. In Lepisosteiformes, orthologs of these lost genes show high or detectable expression levels in the skin but not in the jaw. We thus hypothesize that many SCPP genes located in this cluster are involved in the formation of ganoid scales in Polypteriformes and Lepisosteiformes, and that their orthologs and ganoid scales were convergently lost in Acipenseriformes and Amiiformes.

Regenerating zebrafish scales express a subset of evolutionary conserved genes involved in human skeletal disease

Background

Scales are mineralised exoskeletal structures that are part of the dermal skeleton. Scales have been mostly lost during evolution of terrestrial vertebrates whilst bony fish have retained a mineralised dermal skeleton in the form of fin rays and scales. Each scale is a mineralised collagen plate that is decorated with both matrix-building and resorbing cells. When removed, an ontogenetic scale is quickly replaced following differentiation of the scale pocket-lining cells that regenerate a scale. Processes promoting de novo matrix formation and mineralisation initiated during scale regeneration are poorly understood. Therefore, we performed transcriptomic analysis to determine gene networks and their pathways involved in dermal scale regeneration.

Results

We defined the transcriptomic profiles of ontogenetic and regenerating scales of zebrafish and identified 604 differentially expressed genes (DEGs). These were enriched for extracellular matrix, ossification, and cell adhesion pathways, but not in enamel or dentin formation processes indicating that scales are reminiscent to bone. Hypergeometric tests involving monogenetic skeletal disorders showed that DEGs were strongly enriched for human orthologues that are mutated in low bone mass and abnormal bone mineralisation diseases (P< 2× 10⁻³). The DEGs were also enriched for human orthologues associated with polygenetic skeletal traits, including height (P< 6× 10⁻⁴), and estimated bone mineral density (eBMD, P< 2× 10⁻⁵). Zebrafish mutants of two human orthologues that were robustly associated with height (COL11A2, P=6× 10⁻²⁴) or eBMD (SPP1, P=6× 10⁻²⁰) showed both exo- and endo- skeletal abnormalities as predicted by our genetic association analyses; col11a2^Y228X/Y228X mutants showed exoskeletal and endoskeletal features consistent with abnormal growth, whereas spp1^P160X/P160X mutants predominantly showed mineralisation defects.

Conclusion

We show that scales have a strong osteogenic expression profile comparable to other elements of the dermal skeleton, enriched in genes that favour collagen matrix growth. Despite the many differences between scale and endoskeletal developmental processes, we also show that zebrafish scales express an evolutionarily conserved sub-population of genes that are relevant to human skeletal disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01209-8.

Divergent Expression of SPARC, SPARC-L, and SCPP Genes During Jawed Vertebrate Cartilage Mineralization

While cartilage is an ancient tissue found both in protostomes and deuterostomes, its mineralization evolved more recently, within the vertebrate lineage. SPARC, SPARC-L, and the SCPP members (Secretory Calcium-binding PhosphoProtein genes which evolved from SPARC-L) are major players of dentine and bone mineralization, but their involvement in the emergence of the vertebrate mineralized cartilage remains unclear. We performed in situ hybridization on mineralizing cartilaginous skeletal elements of the frog Xenopus tropicalis (Xt) and the shark Scyliorhinus canicula (Sc) to examine the expression of SPARC (present in both species), SPARC-L (present in Sc only) and the SCPP members (present in Xt only). We show that while mineralizing cartilage expresses SPARC (but not SPARC-L) in Sc, it expresses the SCPP genes (but not SPARC) in Xt, and propose two possible evolutionary scenarios to explain these opposite expression patterns. In spite of these genetic divergences, our data draw the attention on an overlooked and evolutionarily conserved peripheral cartilage subdomain expressing SPARC or the SCPP genes and exhibiting a high propensity to mineralize.

Cells at the Edge: The Dentin-Bone Interface in Zebrafish Teeth

Bone-producing osteoblasts and dentin-producing odontoblasts are closely related cell types, a result from their shared evolutionary history in the ancient dermal skeleton. In mammals, the two cell types can be distinguished based on histological characters and the cells’ position in the pulp cavity or in the tripartite periodontal complex. Different from mammals, teleost fish feature a broad diversity in tooth attachment modes, ranging from fibrous attachment to firm ankylosis to the underlying bone. The connection between dentin and jaw bone is often mediated by a collar of mineralized tissue, a part of the dental unit that has been termed “bone of attachment”. Its nature (bone, dentin, or an intermediate tissue type) is still debated. Likewise, there is a debate about the nature of the cells secreting this tissue: osteoblasts, odontoblasts, or yet another (intermediate) type of scleroblast. Here, we use expression of the P/Q rich secretory calcium-binding phosphoprotein 5 (scpp5) to characterize the cells lining the so-called bone of attachment in the zebrafish dentition. scpp5 is expressed in late cytodifferentiation stage odontoblasts but not in the cells depositing the “bone of attachment”. nor in bona fide osteoblasts lining the supporting pharyngeal jaw bone. Together with the presence of the osteoblast marker Zns-5, and the absence of covering epithelium, this links the cells depositing the “bone of attachment” to osteoblasts rather than to odontoblasts. The presence of dentinal tubule-like cell extensions and the near absence of osteocytes, nevertheless distinguishes the “bone of attachment” from true bone. These results suggest that the “bone of attachment” in zebrafish has characters intermediate between bone and dentin, and, as a tissue, is better termed “dentinous bone”. In other teleosts, the tissue may adopt different properties. The data furthermore support the view that these two tissues are part of a continuum of mineralized tissues. Expression of scpp5 can be a valuable tool to investigate how differentiation pathways diverge between osteoblasts and odontoblasts in teleost models and help resolving the evolutionary history of tooth attachment structures in actinopterygians.

Somite Compartments in Amphioxus and Its Implications on the Evolution of the Vertebrate Skeletal Tissues

Mineralized skeletal tissues of vertebrates are an evolutionary novelty within the chordate lineage. While the progenitor cells that contribute to vertebrate skeletal tissues are known to have two embryonic origins, the mesoderm and neural crest, the evolutionary origin of their developmental process remains unclear. Using cephalochordate amphioxus as our model, we found that cells at the lateral wall of the amphioxus somite express SPARC (a crucial gene for tissue mineralization) and various collagen genes. During development, some of these cells expand medially to surround the axial structures, including the neural tube, notochord and gut, while others expand laterally and ventrally to underlie the epidermis. Eventually these cell populations are found closely associated with the collagenous matrix around the neural tube, notochord, and dorsal aorta, and also with the dense collagen sheets underneath the epidermis. Using known genetic markers for distinct vertebrate somite compartments, we showed that the lateral wall of amphioxus somite likely corresponds to the vertebrate dermomyotome and lateral plate mesoderm. Furthermore, we demonstrated a conserved role for BMP signaling pathway in somite patterning of both amphioxus and vertebrates. These results suggest that compartmentalized somites and their contribution to primitive skeletal tissues are ancient traits that date back to the chordate common ancestor. The finding of SPARC-expressing skeletal scaffold in amphioxus further supports previous hypothesis regarding SPARC gene family expansion in the elaboration of the vertebrate mineralized skeleton.

Coevolution of enamel, ganoin, enameloid, and their matrix SCPP genes in osteichthyans

We resolve debate over the evolution of vertebrate hypermineralized tissues through analyses of matrix protein-encoding secretory calcium-binding phosphoprotein (SCPP) genes and phylogenetic inference of hypermineralized tissues. Among these genes, AMBN and ENAM are found in both sarcopterygians and actinopterygians, whereas AMEL and SCPP5 are found only in sarcopterygians and actinopterygians, respectively. Actinopterygian AMBN, ENAM, and SCPP5 are expressed during the formation of hypermineralized tissues on scales and teeth: ganoin, acrodin, and collar enamel in gar, and acrodin and collar enameloid in zebrafish. Our phylogenetic analyses indicate the emergence of an ancestral enamel in stem-osteichthyans, whereas ganoin emerged in stem-actinopterygians and true enamel in stem-sarcopterygians. Thus, AMBN and ENAM originated in concert with ancestral enamel, SCPP5 evolved in association with ganoin, and AMEL evolved with true enamel. Shifts in gene expression domain and timing explain the evolution of different hypermineralized tissues. We propose that hypermineralized tissues in osteichthyans coevolved with matrix SCPP genes.

•

Ganoin emerged in actinopterygians; true enamel arose in sarcopterygians

•

Dental enamel, acrodin, and enameloid in actinopterygians are related to ganoin

•

SCPP5 evolved in association with ganoin, whereas AMEL evolved with true enamel

•

Shifts in SCPP gene expression explain the evolution of hypermineralized tissues

Grand Challenges in Comparative Tooth Biology

Teeth are a model system for integrating developmental genomics, functional morphology, and evolution. We are at the cusp of being able to address many open issues in comparative tooth biology and we outline several of these newly tractable and exciting research directions. Like never before, technological advances and methodological approaches are allowing us to investigate the developmental machinery of vertebrates and discover both conserved and excitingly novel mechanisms of diversification. Additionally, studies of the great diversity of soft tissues, replacement teeth, and non-trophic functions of teeth are providing new insights into dental diversity. Finally, we highlight several emerging model groups of organisms that are at the forefront of increasing our appreciation of the mechanisms underlying tooth diversification.

Insights into the Evolution of Neoteny from the Genome of the Asian Icefish Protosalanx chinensis

Salangids, known as Asian icefishes, represent a peculiar radiation within the bony fish order Protacanthopterygii where adult fish retain larval characteristics such as transparent and miniaturized bodies and a cartilaginous endoskeleton into adulthood. Here, we report a de novo genome of Protosalanx chinensis, the most widely distributed salangid lineage. The P. chinensis genome assembly is more contiguous and complete than a previous assembly. We estimate that P. chinensis, salmons, trouts, and pikes diverged from a common ancestor 185 million years ago. A juxtaposition with other fish genomes revealed loss of the genes encoding ectodysplasin-A receptor (EDAR), SCPP1, and four Hox proteins and likely lack of canonical fibroblast growth factor 5 (FGF5) function. We also report genomic variations of P. chinensis possibly reflecting the immune system repertoire of a species with a larval phenotype in sexually mature individuals. The new Asian icefish reference genome provides a solid foundation for future studies.

Mutant dlx3b disturbs normal tooth mineralization and bone formation in zebrafish

Background

Tricho-dento-osseous (TDO) syndrome is an autosomal dominant disorder characterized by anomalies in hair, teeth and bone (OMIM190320). Various mutations of Distal-Less 3 (DLX3) gene are found to be responsible for human TDO. The aim of this study was to investigate effects of DLX3 on tooth and bone development using a zebrafish model.

Methods

The dlx3b mutant zebrafish lines were established using the gene targeting tool transcription activator-like effector nuclease (TALEN). Micro-computed tomography was used to render the three-dimensional skeletal structures of mutant fishes. The pharyngeal bone along with connected teeth was isolated and stained by Alizarine Red S, then observed under stereomicroscope. Scanning electron microscopy (SEM) and energy dispersive spectrometer (EDS) were used to examine the tooth surface morphology and mineral composition. Quantitative real-time PCR was used to analyze gene expression.

Results

A moderate curvature of the spine toward the dorsal side was found at the early larval stages, appearing in 86 out of 100 larvae in dlx3b^-/- group as compared to 3 out of 99 in the dlx3b^+/+ group. At the adult stage, three of the thirty dlx3b^-/- homozygotes exhibited prominent abnormal curvature in the spine. SEM revealed morphological surface changes in pharyngeal teeth enameloid, accompanied by a decrease in the mineral content detected by EDS. Furthermore, specific secretory calcium-binding phosphoprotein (SCPP) genes, including odam, scpp9, spp1, scpp1, and scpp5 were significantly downregulated in dlx3b mutants.

Conclusion

The findings of this study suggest that dlx3b is critical for enamel mineralization and bone formation in zebrafish. Moreover, the discovery of the downregulation of SCPP genes in dlx3b mutants sheds new light on the molecular mechanisms underlying TDO syndrome.

Comprehensive RNA-Seq analysis of notochord-enriched genes induced during Xenopus tropicalis tail resorption

Anuran metamorphosis is perhaps the most dramatic developmental process regulated by thyroid hormone (TH). One of the unique processes that occur during metamorphosis is the complete resorption of the tail, including the notochord. Interestingly, recent gene knockout studies have shown that of the two known vertebrate TH receptors, TRα and TRβ, TRβ appears to be critical for notochord regression during tail resorption in Xenopus tropicalis. To determine the mechanisms underlying notochord regression, we carried out a comprehensive gene expression analysis in the notochord during metamorphosis by using RNA-Seq analyses of whole tail at stage 60 before any noticeable tail length reduction, whole tail at stage 63 when the tail length is reduced by about one half, and the rest of the tail at stage 63 after removing the notochord. This allowed us to identify many notochord-enriched, metamorphosis-induced genes at stage 63. Future studies on these genes should help to determine if they are regulated by TRβ and play any roles in notochord regression.

Skeletal Mineralization in Association with Type X Collagen Expression Is an Ancestral Feature for Jawed Vertebrates

In order to characterize the molecular bases of mineralizing cell evolution, we targeted type X collagen, a nonfibrillar network forming collagen encoded by the Col10a1 gene. It is involved in the process of endochondral ossification in ray-finned fishes and tetrapods (Osteichthyes), but until now unknown in cartilaginous fishes (Chondrichthyes). We show that holocephalans and elasmobranchs have respectively five and six tandemly duplicated Col10a1 gene copies that display conserved genomic synteny with osteichthyan Col10a1 genes. All Col10a1 genes in the catshark Scyliorhinus canicula are expressed in ameloblasts and/or odontoblasts of teeth and scales, during the stages of extracellular matrix protein secretion and mineralization. Only one duplicate is expressed in the endoskeletal (vertebral) mineralizing tissues. We also show that the expression of type X collagen is present in teeth of two osteichthyans, the zebrafish Danio rerio and the western clawed frog Xenopus tropicalis, indicating an ancestral jawed vertebrate involvement of type X collagen in odontode formation. Our findings push the origin of Col10a1 gene prior to the divergence of osteichthyans and chondrichthyans, and demonstrate its ancestral association with mineralization of both the odontode skeleton and the endoskeleton.

Immunolocalization of enamel matrix protein-like proteins in the tooth enameloid of spotted gar, Lepisosteus oculatus, an actinopterygian bony fish

Enameloid is a well-mineralized tissue covering the tooth surface in fish and it corresponds to the outer-most layer of dentin. It was reported that both dental epithelial cells and odontoblasts are involved in the formation of enameloid. Nevertheless, the localization and timing of secretion of ectodermal enamel matrix proteins in enameloid are unclear. In the present study, the enameloid matrix during the stages of enameloid formation in spotted gar, Lepisosteus oculatus, an actinopterygian, was examined mainly by transmission electron microscopy-based immunohistochemistry using an anti-mammalian amelogenin antibody and antiserum. Positive immunoreactivity with the antibody and antiserum was found in enameloid from the surface to the dentin-enameloid junction just before the formation of crystallites. This immunoreactivity disappeared rapidly before the full appearance of crystallites in the enameloid during the stage of mineralization. Immunolabelling was usually found along the collagen fibrils but was not seen on the electron-dense fibrous structures, which were probably derived from matrix vesicles in the previous stage. In inner dental epithelial cells, the granules in the distal cytoplasm often showed positive immunoreactivity, suggesting that the enamel matrix protein-like proteins originated from inner dental epithelial cells. Enamel matrix protein-like proteins in the enameloid matrix might be common to the enamel matrix protein-like proteins previously reported in the collar enamel of teeth and ganoine of ganoid scales, because they exhibited marked immunoreactivity with the same anti-mammalian amelogenin antibodies. It is likely that enamel matrix protein-like proteins are involved in the formation of crystallites along collagen fibrils in enameloid.

Predictable genome-wide sorting of standing genetic variation during parallel adaptation to basic versus acidic environments in stickleback fish

Genomic studies of parallel (or convergent) evolution often compare multiple populations diverged into two ecologically different habitats to search for loci repeatedly involved in adaptation. Because the shared ancestor of these populations is generally unavailable, the source of the alleles at adaptation loci, and the direction in which their frequencies were shifted during evolution, remain elusive. To shed light on these issues, we here use multiple populations of threespine stickleback fish adapted to two different types of derived freshwater habitats-basic and acidic lakes on the island of North Uist, Outer Hebrides, Scotland-and the present-day proxy of their marine ancestor. In a first step, we combine genome-wide pooled sequencing and targeted individual-level sequencing to demonstrate that ecological and phenotypic parallelism in basic-acidic divergence is reflected by genomic parallelism in dozens of genome regions. Exploiting data from the ancestor, we next show that the acidic populations, residing in ecologically more extreme derived habitats, have adapted by accumulating alleles rare in the ancestor, whereas the basic populations have retained alleles common in the ancestor. Genomic responses to selection are thus predictable from the ecological difference of each derived habitat type from the ancestral one. This asymmetric sorting of standing genetic variation at loci important to basic-acidic divergence has further resulted in more numerous selective sweeps in the acidic populations. Finally, our data suggest that the maintenance in marine fish of standing variation important to adaptive basic-acidic differentiation does not require extensive hybridization between the marine and freshwater populations. Overall, our study reveals striking genome-wide determinism in both the loci involved in parallel divergence, and in the direction in which alleles at these loci have been selected.

Evolution of dental tissue mineralization: an analysis of the jawed vertebrate SPARC and SPARC-L families

BACKGROUND

The molecular bases explaining the diversity of dental tissue mineralization across gnathostomes are still poorly understood. Odontodes, such as teeth and body denticles, are serial structures that develop through deployment of a gene regulatory network shared between all gnathostomes. Dentin, the inner odontode mineralized tissue, is produced by odontoblasts and appears well-conserved through evolution. In contrast, the odontode hypermineralized external layer (enamel or enameloid) produced by ameloblasts of epithelial origin, shows extensive structural variations. As EMP (Enamel Matrix Protein) genes are as yet only found in osteichthyans where they play a major role in the mineralization of teeth and others skeletal organs, our understanding of the molecular mechanisms leading to the mineralized odontode matrices in chondrichthyans remains virtually unknown.

RESULTS

We undertook a phylogenetic analysis of the SPARC/SPARC-L gene family, from which the EMPs are supposed to have arisen, and examined the expression patterns of its members and of major fibrillar collagens in the spotted catshark Scyliorhinus canicula, the thornback ray Raja clavata, and the clawed frog Xenopus tropicalis. Our phylogenetic analyses reveal that the single chondrichthyan SPARC-L gene is co-orthologous to the osteichthyan SPARC-L1 and SPARC-L2 paralogues. In all three species, odontoblasts co-express SPARC and collagens. In contrast, ameloblasts do not strongly express collagen genes but exhibit strikingly similar SPARC-L and EMP expression patterns at their maturation stage, in the examined chondrichthyan and osteichthyan species, respectively.

CONCLUSIONS

A well-conserved odontoblastic collagen/SPARC module across gnathostomes further confirms dentin homology. Members of the SPARC-L clade evolved faster than their SPARC paralogues, both in terms of protein sequence and gene duplication. We uncover an osteichthyan-specific duplication that produced SPARC-L1 (subsequently lost in pipidae frogs) and SPARC-L2 (independently lost in teleosts and tetrapods).Our results suggest the ameloblastic expression of the single chondrichthyan SPARC-L gene at the maturation stage reflects the ancestral gnathostome situation, and provide new evidence in favor of the homology of enamel and enameloids in all gnathostomes.

A Genomic Survey of SCPP Family Genes in Fishes Provides Novel Insights into the Evolution of Fish Scales

The family of secretory calcium-binding phosphoproteins (SCPPs) have been considered vital to skeletal tissue mineralization. However, most previous SCPP studies focused on phylogenetically distant animals but not on those closely related species. Here we provide novel insights into the coevolution of SCPP genes and fish scales in 10 species from Otophysi. According to their scale phenotypes, these fishes can be divided into three groups, i.e., scaled, sparsely scaled, and scaleless. We identified homologous SCPP genes in the genomes of these species and revealed an absence of some SCPP members in some genomes, suggesting an uneven evolutionary history of SCPP genes in fishes. In addition, most of these SCPP genes, with the exception of SPP1, individually form one or two gene cluster(s) on each corresponding genome. Furthermore, we constructed phylogenetic trees using maximum likelihood method to estimate their evolution. The phylogenetic topology mostly supports two subclasses in some species, such as Cyprinus carpio, Sinocyclocheilus anshuiensis, S. grahamin, and S. rhinocerous, but not in the other examined fishes. By comparing the gene structures of recently reported candidate genes, SCPP1 and SCPP5, for determining scale phenotypes, we found that the hypothesis is suitable for Astyanax mexicanus, but denied by S. anshuiensis, even though they are both sparsely scaled for cave adaptation. Thus, we conclude that, although different fish species display similar scale phenotypes, the underlying genetic changes however might be diverse. In summary, this paper accelerates the recognition of the SCPP family in teleosts for potential scale evolution.

DENTAL ENAMEL FORMATION AND IMPLICATIONS FOR ORAL HEALTH AND DISEASE

Dental enamel is the hardest and most mineralized tissue in extinct and extant vertebrate species and provides maximum durability that allows teeth to function as weapons and/or tools as well as for food processing. Enamel development and mineralization is an intricate process tightly regulated by cells of the enamel organ called ameloblasts. These heavily polarized cells form a monolayer around the developing enamel tissue and move as a single forming front in specified directions as they lay down a proteinaceous matrix that serves as a template for crystal growth. Ameloblasts maintain intercellular connections creating a semi-permeable barrier that at one end (basal/proximal) receives nutrients and ions from blood vessels, and at the opposite end (secretory/apical/distal) forms extracellular crystals within specified pH conditions. In this unique environment, ameloblasts orchestrate crystal growth via multiple cellular activities including modulating the transport of minerals and ions, pH regulation, proteolysis, and endocytosis. In many vertebrates, the bulk of the enamel tissue volume is first formed and subsequently mineralized by these same cells as they retransform their morphology and function. Cell death by apoptosis and regression are the fates of many ameloblasts following enamel maturation, and what cells remain of the enamel organ are shed during tooth eruption, or are incorporated into the tooth's epithelial attachment to the oral gingiva. In this review, we examine key aspects of dental enamel formation, from its developmental genesis to the ever-increasing wealth of data on the mechanisms mediating ionic transport, as well as the clinical outcomes resulting from abnormal ameloblast function.

SCPP Genes and Their Relatives in Gar: Rapid Expansion of Mineralization Genes in Osteichthyans

Gar is an actinopterygian that has bone, dentin, enameloid, and ganoin (enamel) in teeth and/or scales. Mineralization of these tissues involves genes encoding various secretory calcium-binding phosphoproteins (SCPPs) in osteichthyans, but no SCPP genes have been identified in chondrichthyans to date. In the gar genome, we identified 38 SCPP genes, seven of which encode "acidic-residue-rich" proteins and 31 encode "Pro/Gln (P/Q) rich" proteins. These gar SCPP genes constitute the largest known repertoire, including many newly identified P/Q-rich genes expressed in teeth and/or scales. Among gar SCPP genes, six acidic and three P/Q-rich genes were identified as orthologs of sarcopterygian genes. The sarcopterygian orthologs of most of these acidic genes are involved in bone and/or dentin formation, and sarcopterygian orthologs of all three P/Q-rich genes participate in enamel formation. The finding of these genes in gar suggests that an elaborate SCPP gene-based genetic system for tissue mineralization was already present in stem osteichthyans. While SCPP genes have been thought to originate from ancient SPARCL1, SPARCL1L1 appears to be more closely related to these genes, because it established a structure similar to acidic SCPP genes probably in stem gnathostomes, perhaps at about the same time with the origin of tissue mineralization. Assuming enamel evolved in stem osteichthyans, all P/Q-rich SCPP genes likely arose within the osteichthyan lineage. Furthermore, the absence of acidic SCPP genes in chondrichthyans might be explained by the secondary loss of earliest acidic genes. It appears that many SCPP genes expanded rapidly in stem osteichthyans and in basal actinopterygians.

The genome of the largest bony fish, ocean sunfish (Mola mola), provides insights into its fast growth rate

BACKGROUND

The ocean sunfish (Mola mola), which can grow up to a length of 2.7 m and weigh 2.3 tons, is the world's largest bony fish. It has an extremely fast growth rate and its endoskeleton is mainly composed of cartilage. Another unique feature of the sunfish is its lack of a caudal fin, which is replaced by a broad and stiff lobe that results in the characteristic truncated appearance of the fish.

RESULTS

To gain insights into the genomic basis of these phenotypic traits, we sequenced the sunfish genome and performed a comparative analysis with other teleost genomes. Several sunfish genes involved in the growth hormone and insulin-like growth factor 1 (GH/IGF1) axis signalling pathway were found to be under positive selection or accelerated evolution, which might explain its fast growth rate and large body size. A number of genes associated with the extracellular matrix, some of which are involved in the regulation of bone and cartilage development, have also undergone positive selection or accelerated evolution. A comparison of the sunfish genome with that of the pufferfish (fugu), which has a caudal fin, revealed that the sunfish contains more homeobox (Hox) genes although both genomes contain seven Hox clusters. Thus, caudal fin loss in sunfish is not associated with the loss of a specific Hox gene.

CONCLUSIONS

Our analyses provide insights into the molecular basis of the fast growth rate and large size of the ocean sunfish. The high-quality genome assembly generated in this study should facilitate further studies of this 'natural mutant'.

The channel catfish genome sequence provides insights into the evolution of scale formation in teleosts

Catfish represent 12% of teleost or 6.3% of all vertebrate species, and are of enormous economic value. Here we report a high-quality reference genome sequence of channel catfish (Ictalurus punctatus), the major aquaculture species in the US. The reference genome sequence was validated by genetic mapping of 54,000 SNPs, and annotated with 26,661 predicted protein-coding genes. Through comparative analysis of genomes and transcriptomes of scaled and scaleless fish and scale regeneration experiments, we address the genomic basis for the most striking physical characteristic of catfish, the evolutionary loss of scales and provide evidence that lack of secretory calcium-binding phosphoproteins accounts for the evolutionary loss of scales in catfish. The channel catfish reference genome sequence, along with two additional genome sequences and transcriptomes of scaled catfishes, provide crucial resources for evolutionary and biological studies. This work also demonstrates the power of comparative subtraction of candidate genes for traits of structural significance.

The spotted gar genome illuminates vertebrate evolution and facilitates human-teleost comparisons

To connect human biology to fish biomedical models, we sequenced the genome of spotted gar (Lepisosteus oculatus), whose lineage diverged from teleosts before teleost genome duplication (TGD). The slowly evolving gar genome has conserved in content and size many entire chromosomes from bony vertebrate ancestors. Gar bridges teleosts to tetrapods by illuminating the evolution of immunity, mineralization and development (mediated, for example, by Hox, ParaHox and microRNA genes). Numerous conserved noncoding elements (CNEs; often cis regulatory) undetectable in direct human-teleost comparisons become apparent using gar: functional studies uncovered conserved roles for such cryptic CNEs, facilitating annotation of sequences identified in human genome-wide association studies. Transcriptomic analyses showed that the sums of expression domains and expression levels for duplicated teleost genes often approximate the patterns and levels of expression for gar genes, consistent with subfunctionalization. The gar genome provides a resource for understanding evolution after genome duplication, the origin of vertebrate genomes and the function of human regulatory sequences.

The role of the SIBLING, Bone Sialoprotein in skeletal biology - Contribution of mouse experimental genetics

Bone Sialoprotein (BSP) is a member of the "Small Integrin-Binding Ligand N-linked Glycoproteins" (SIBLING) extracellular matrix protein family of mineralized tissues. BSP has been less studied than other SIBLING proteins such as Osteopontin (OPN), which is coexpressed with it in several skeletal cell types. Here we review the contribution of genetically engineered mice (BSP gene knockout and overexpression) to the understanding of the role of BSP in the bone organ. The studies made so far highlight the role of BSP in skeletal mineralization, as well as its importance for proper osteoblast and osteoclast differentiation and activity, most prominently in primary/repair bone. The absence of BSP also affects the local environment of the bone tissue, in particular hematopoiesis and vascularization. Interestingly, lack of BSP induces an overexpression of OPN, and the cognate protein could be responsible for some aspects of the BSP gene knockout skeletal phenotype, while replacing BSP for some of its functions. Such interplay between the partly overlapping functions of SIBLING proteins, as well as the network of cross-regulations in which they are involved should now be the focus of further work.

New genomic and fossil data illuminate the origin of enamel

Enamel, the hardest vertebrate tissue, covers the teeth of almost all sarcopterygians (lobe-finned bony fishes and tetrapods) as well as the scales and dermal bones of many fossil lobe-fins. Enamel deposition requires an organic matrix containing the unique enamel matrix proteins (EMPs) amelogenin (AMEL), enamelin (ENAM) and ameloblastin (AMBN). Chondrichthyans (cartilaginous fishes) lack both enamel and EMP genes. Many fossil and a few living non-teleost actinopterygians (ray-finned bony fishes) such as the gar, Lepisosteus, have scales and dermal bones covered with a proposed enamel homologue called ganoine. However, no gene or transcript data for EMPs have been described from actinopterygians. Here we show that Psarolepis romeri, a bony fish from the the Early Devonian period, combines enamel-covered dermal odontodes on scales and skull bones with teeth of naked dentine, and that Lepisosteus oculatus (the spotted gar) has enam and ambn genes that are expressed in the skin, probably associated with ganoine formation. The genetic evidence strengthens the hypothesis that ganoine is homologous with enamel. The fossil evidence, further supported by the Silurian bony fish Andreolepis, which has enamel-covered scales but teeth and odontodes on its dermal bones made of naked dentine, indicates that this tissue originated on the dermal skeleton, probably on the scales. It subsequently underwent heterotopic expansion across two highly conserved patterning boundaries (scales/head-shoulder and dermal/oral) within the odontode skeleton.

Evolutionary analysis of selective constraints identifies ameloblastin (AMBN) as a potential candidate for amelogenesis imperfecta

Background

Ameloblastin (AMBN) is a phosphorylated, proline/glutamine-rich protein secreted during enamel formation. Previous studies have revealed that this enamel matrix protein was present early in vertebrate evolution and certainly plays important roles during enamel formation although its precise functions remain unclear. We performed evolutionary analyses of AMBN in order to (i) identify residues and motifs important for the protein function, (ii) predict mutations responsible for genetic diseases, and (iii) understand its molecular evolution in mammals.

Results

In silico searches retrieved 56 complete sequences in public databases that were aligned and analyzed computationally. We showed that AMBN is globally evolving under moderate purifying selection in mammals and contains a strong phylogenetic signal. In addition, our analyses revealed codons evolving under significant positive selection. Evidence for positive selection acting on AMBN was observed in catarrhine primates and the aye-aye. We also found that (i) an additional translation initiation site was recruited in the ancestral placental AMBN, (ii) a short exon was duplicated several times in various species including catarrhine primates, and (iii) several polyadenylation sites are present.

Conclusions

AMBN possesses many positions, which have been subjected to strong selective pressure for 200 million years. These positions correspond to several cleavage sites and hydroxylated, O-glycosylated, and phosphorylated residues. We predict that these conserved positions would be potentially responsible for enamel disorder if substituted. Some motifs that were previously identified as potentially important functionally were confirmed, and we found two, highly conserved, new motifs, the function of which should be tested in the near future. This study illustrates the power of evolutionary analyses for characterizing the functional constraints acting on proteins with yet uncharacterized structure.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0431-0) contains supplementary material, which is available to authorized users.

[Enamel: a unique self-assembling in mineral world]

Enamel is a unique tissue in vertebrates, acellular, formed on a labile scaffolding matrix and hypermineralized. The ameloblasts are epithelial cells in charge of amelogenesis. They secrete a number of matrix proteins degraded by enzymes during enamel mineralization. This ordered cellular and extracellular events imply that any genetic or environmental perturbation will produce indelible and recognizable defects. The specificity of defects will indicate the affected cellular process. Thus, depending on the specificity of alterations, the teratogenic event can be retrospectively established. Advances in the field allow to use enamel defects as diagnostic tools for molecular disorders. The multifunctionality of enamel peptides is presently identified from their chemical roles in mineralization to cell signaling, constituting a source of concrete innovations in regenerative medicine.

SCPP genes in the coelacanth: tissue mineralization genes shared by sarcopterygians

The coelacanth is the basal-most extant sarcopterygian that has teeth and tooth-like structures, comprising bone, dentin, and enamel or enameloid. Formation of these tissues involves many members of the secretory calcium-binding protein (SCPP) family. In tetrapods, acidic-residue-rich SCPPs are used in mineralization of bone and dentin, whereas Pro/Gln-rich SCPPs participate in enamel formation. Teleosts also employ many SCPPs for tissue mineralization. Nevertheless, the repertoire of SCPPs is largely different in teleosts and tetrapods; hence, filling this gap would be critical to elucidate early evolution of mineralized tissues in osteichthyans. In the present study, we searched for SCPP genes in the coelacanth genome and identified 11, of which two have clear orthologs in both tetrapods and teleosts, seven only in tetrapods, and two in neither of them. Given the divergence times of these vertebrate lineages, our discovery of this many SCPP genes shared between the coelacanth and tetrapods, but not with teleosts, suggests a complicated evolutionary scheme of SCPP genes in early osteichthyans. Our investigation also revealed both conserved and derived characteristics of SCPPs in the coelacanth and other vertebrates. Notably, acidic SCPPs independently evolved various acidic repeats in different lineages, while maintaining high acidity, presumably important for interactions with calcium. Furthermore, the three Pro/Gln-rich SCPP genes, required for mineralizing enamel matrix and confirmed only in tetrapods, were all identified in the coelacanth, strongly suggesting that enamel is equivalent in the coelacanth and tetrapods. This finding corroborates the previous proposition that true enamel evolved much earlier than the origin of tetrapods.

Characterisation of secretory calcium-binding phosphoprotein-proline-glutamine-rich 1: a novel basal lamina component expressed at cell-tooth interfaces

Functional genomic screening of the rat enamel organ (EO) has led to the identification of a number of secreted proteins expressed during the maturation stage of amelogenesis, including amelotin (AMTN) and odontogenic ameloblast-associated (ODAM). In this study, we characterise the gene, protein and pattern of expression of a related protein called secretory calcium-binding phosphoprotein-proline-glutamine-rich 1 (SCPPPQ1). The Scpppq1 gene resides within the secretory calcium-binding phosphoprotein (Scpp) cluster. SCPPPQ1 is a highly conserved, 75-residue, secreted protein rich in proline, leucine, glutamine and phenylalanine. In silico data mining has revealed no correlation to any known sequences. Northern blotting of various rat tissues suggests that the expression of Scpppq1 is restricted to tooth and associated tissues. Immunohistochemical analyses show that the protein is expressed during the late maturation stage of amelogenesis and in the junctional epithelium where it localises to an atypical basal lamina at the cell-tooth interface. This discrete localisation suggests that SCPPPQ1, together with AMTN and ODAM, participates in structuring the basal lamina and in mediating attachment of epithelia cells to mineralised tooth surfaces.