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

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.

Common features of cartilage maturation are not conserved in an amphibian model

BACKGROUND

Mouse, chick, and zebrafish undergo a highly conserved program of cartilage maturation during endochondral ossification (bone formation via a cartilage template). Standard histological and molecular features of cartilage maturation are chondrocyte hypertrophy, downregulation of the chondrogenic markers Sox9 and Col2a1, and upregulation of Col10a1. We tested whether cartilage maturation is conserved in an amphibian, the western clawed frog Xenopus tropicalis, using in situ hybridization for standard markers and a novel laser-capture microdissection RNAseq data set. We also functionally tested whether thyroid hormone drives cartilage maturation in X tropicalis, as it does in other vertebrates.

RESULTS

The developing frog humerus mostly followed the standard progression of cartilage maturation. Chondrocytes gradually became hypertrophic as col2a1 and sox9 were eventually down-regulated, but col10a1 was not up-regulated. However, the expression levels of several genes associated with the early formation of cartilage, such as acan, sox5, and col9a2, remained highly expressed even as humeral chondrocytes matured. Greater deviances were observed in head cartilages, including the ceratohyal, which underwent hypertrophy within hours of becoming cartilaginous, maintained relatively high levels of col2a1 and sox9, and lacked col10a1 expression. Interestingly, treating frog larvae with thyroid hormone antagonists did not specifically reduce head cartilage hypertrophy, resulting rather in a global developmental delay.

CONCLUSION

These data reveal that basic cartilage maturation features in the head, and to a lesser extent in the limb, are not conserved in X tropicalis. Future work revealing how frogs deviate from the standard cartilage maturation program might shed light on both evolutionary and health studies.

Transcriptomic landscape of Atlantic salmon (Salmo salar L.) skin

In this study, we present the first spatial transcriptomic atlas of Atlantic salmon skin using the Visium Spatial Gene Expression protocol. We utilized frozen skin tissue from 4 distinct sites, namely the operculum, pectoral and caudal fins, and scaly skin at the flank of the fish close to the lateral line, obtained from 2 Atlantic salmon (150 g). High-quality frozen tissue sections were obtained by embedding tissue in optimal cutting temperature media prior to freezing and sectioning. Further, we generated libraries and spatial transcriptomic maps, achieving a minimum of 80 million reads per sample with mapping efficiencies ranging from 79.3 to 89.4%. Our analysis revealed the detection of over 80,000 transcripts and nearly 30,000 genes in each sample. Among the tissue types observed in the skin, the epithelial tissues exhibited the highest number of transcripts (unique molecular identifier counts), followed by muscle tissue, loose and fibrous connective tissue, and bone. Notably, the widest nodes in the transcriptome network were shared among the epithelial clusters, while dermal tissues showed less consistency, which is likely attributable to the presence of multiple cell types at different body locations. Additionally, we identified collagen type 1 as the most prominent gene family in the skin, while keratins were found to be abundant in the epithelial tissue. Furthermore, we successfully identified gene markers specific to epithelial tissue, bone, and mesenchyme. To validate their expression patterns, we conducted a meta-analysis of the microarray database, which confirmed high expression levels of these markers in mucosal organs, skin, gills, and the olfactory rosette.

Exposure of primary osteoblasts to combined magnetic and electric fields induced spatiotemporal endochondral ossification characteristic gene- and protein expression profiles

Purpose

Molecular processes in primary osteoblasts were analyzed in response to magnetic and electric field exposure to examine its potential impact on bone healing.

Methods

Primary osteoblasts were exposed to a combination of a magnetic field and an additional electric field (EFMF) (20 Hz, 700 mV, 5 mT, continuous sinusoids) in vitro. mRNA- and protein-expressions were assessed during a time interval of 21 days and compared with expression data obtained from control osteoblasts.

Results

We observed an autonomous osteoblast differentiation process in vitro under the chosen cultivation conditions. The initial proliferative phase was characterized by a constitutively high mRNA expression of extracellular matrix proteins. Concurrent EFMF exposure resulted in significanly increased cell proliferation (fold change: 1.25) and reduced mRNA-expressions of matrix components (0.5–0.75). The following reorganization of the extracellular matrix is prerequisite for matrix mineralization and is characterised by increased Ca²⁺ deposition (1.44). On molecular level EFMF exposure led to a significant decreased thrombospondin 1 (THBS1) mRNA- (0.81) and protein- (0.54) expression, which in turn reduced the TGFß1-dependent mRNA- (0.68) and protein- (0.5) expression of transforming growth factor beta induced (ßIG-H3) significantly, an inhibitor of endochondral ossification. Consequently, EFMF exposure stimulated the expression of genes characteristic for endochondral ossification, such as collagen type 10, A1 (1.50), osteopontin (1.50) and acellular communication network factor 3 (NOV) (1.45).

Conclusions

In vitro exposure of osteoblasts to EFMF supports cell differentiation and induces gene- and protein-expression patterns characteristic for endochondral ossification during bone fracture healing in vivo.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40634-022-00477-9.

Probiotics Enhance Bone Growth and Rescue BMP Inhibition: New Transgenic Zebrafish Lines to Study Bone Health

Zebrafish larvae, especially gene-specific mutants and transgenic lines, are increasingly used to study vertebrate skeletal development and human pathologies such as osteoporosis, osteopetrosis and osteoarthritis. Probiotics have been recognized in recent years as a prophylactic treatment for various bone health issues in humans. Here, we present two new zebrafish transgenic lines containing the coding sequences for fluorescent proteins inserted into the endogenous genes for sp7 and col10a1a with larvae displaying fluorescence in developing osteoblasts and the bone extracellular matrix (mineralized or non-mineralized), respectively. Furthermore, we use these transgenic lines to show that exposure to two different probiotics, Bacillus subtilis and Lactococcus lactis, leads to an increase in osteoblast formation and bone matrix growth and mineralization. Gene expression analysis revealed the effect of the probiotics, particularly Bacillus subtilis, in modulating several skeletal development genes, such as runx2, sp7, spp1 and col10a1a, further supporting their ability to improve bone health. Bacillus subtilis was the more potent probiotic able to significantly reverse the inhibition of bone matrix formation when larvae were exposed to a BMP inhibitor (LDN212854).

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.

Mineralized Cartilage and Bone-Like Tissues in Chondrichthyans Offer Potential Insights Into the Evolution and Development of Mineralized Tissues in the Vertebrate Endoskeleton

The impregnation of biominerals into the extracellular matrix of living organisms, a process termed biomineralization, gives rise to diverse mineralized (or calcified) tissues in vertebrates. Preservation of mineralized tissues in the fossil record has provided insights into the evolutionary history of vertebrates and their skeletons. However, current understanding of the vertebrate skeleton and of the processes underlying its formation is biased towards biomedical models such as the tetrapods mouse and chick. Chondrichthyans (sharks, skates, rays, and chimaeras) and osteichthyans are the only vertebrate groups with extant (living) representatives that have a mineralized skeleton, but the basal phylogenetic position of chondrichthyans could potentially offer unique insights into skeletal evolution. For example, bone is a vertebrate novelty, but the internal supporting skeleton (endoskeleton) of extant chondrichthyans is commonly described as lacking bone. The molecular and developmental basis for this assertion is yet to be tested. Subperichondral tissues in the endoskeleton of some chondrichthyans display mineralization patterns and histological and molecular features of bone, thereby challenging the notion that extant chondrichthyans lack endoskeletal bone. Additionally, the chondrichthyan endoskeleton demonstrates some unique features and others that are potentially homologous with other vertebrates, including a polygonal mineralization pattern, a trabecular mineralization pattern, and an unconstricted perichordal sheath. Because of the basal phylogenetic position of chondrichthyans among all other extant vertebrates with a mineralized skeleton, developmental and molecular studies of chondrichthyans are critical to flesh out the evolution of vertebrate skeletal tissues, but only a handful of such studies have been carried out to date. This review discusses morphological and molecular features of chondrichthyan endoskeletal tissues and cell types, ultimately emphasizing how comparative embryology and transcriptomics can reveal homology of mineralized skeletal tissues (and their cell types) between chondrichthyans and other vertebrates.

Enzymatic Approach in Calcium Phosphate Biomineralization: A Contribution to Reconcile the Physicochemical with the Physiological View

Biomineralization is the process by which organisms produce hard inorganic matter from soft tissues with outstanding control of mineral deposition in time and space. For this purpose, organisms deploy a sophisticated "toolkit" that has resulted in significant evolutionary innovations, for which calcium phosphate (CaP) is the biomineral selected for the skeleton of vertebrates. While CaP mineral formation in aqueous media can be investigated by studying thermodynamics and kinetics of phase transitions in supersaturated solutions, biogenic mineralization requires coping with the inherent complexity of biological systems. This mainly includes compartmentalization and homeostatic processes used by organisms to regulate key physiological factors, including temperature, pH and ion concentration. A detailed analysis of the literature shows the emergence of two main views describing the mechanism of CaP biomineralization. The first one, more dedicated to the study of in vivo systems and supported by researchers in physiology, often involves matrix vesicles (MVs). The second one, more investigated by the physicochemistry community, involves collagen intrafibrillar mineralization particularly through in vitro acellular models. Herein, we show that there is an obvious need in the biological systems to control both where and when the mineral forms through an in-depth survey of the mechanism of CaP mineralization. This necessity could gather both communities of physiologists and physicochemists under a common interest for an enzymatic approach to better describe CaP biomineralization. Both homogeneous and heterogeneous enzymatic catalyses are conceivable for these systems, and a few preliminary promising results on CaP mineralization for both types of enzymatic catalysis are reported in this work. Through them, we aim to describe the relevance of our point of view and the likely findings that could be obtained when adding an enzymatic approach to the already rich and creative research field dealing with CaP mineralization. This complementary approach could lead to a better understanding of the biomineralization mechanism and inspire the biomimetic design of new materials.

Expression Profiling and Functional Analysis of Candidate Col10a1 Regulators Identified by the TRAP Program

Hypertrophic chondrocytes and their specific marker, the type X collagen gene (Col10a1), are critical components of endochondral bone formation during skeletal development. We previously found that Runx2 is an indispensable mouse Col10a1 gene regulator and identified many other transcription factors (TFs) that potentially interact with the 150-bp Col10a1 cis-enhancer. However, the roles of these candidate TFs in Col10a1 expression and chondrocyte hypertrophy have not been elucidated. Here, we focus on 32 candidate TFs recently identified by analyzing the 150-bp Col10a1 enhancer using the transcription factor affinity prediction (TRAP) program. We found that 12 TFs (Hoxa3, Lsx, Evx2, Dlx5, S8, Pax2, Egr2, Mef2a, Barhl2, GKlf, Sox17, and Crx) were significantly upregulated and four TFs (Lhx4, Tbx5, Mef2c, and Hb9) were significantly downregulated in hypertrophic MCT cells, which show upregulation of Col10a1 expression. Most of the differential expression pattern of these TFs conformed with the results obtained from ATDC5 cell model and primary mouse chondrocytes. Notably, Tbx5 was downregulated upon Col10a1 upregulation, overexpression of Tbx5 decreased Col10a1 expression, and knock-down of Tbx5 increased Col10a1 expression in hypertrophic chondrocytes, suggesting that Tbx5 is a negative regulator of Col10a1. We further generated a stable Tbx5-overexpressing ATDC5 cell line and ColX-Tbx5 transgenic mice driven by Col10a1-specific enhancers and promoters. Tbx5 overexpression decreased Col10a1 expression in ATDC5 cells cultured as early as day 7 and in limb tissue on post-natal day 1. Slightly weaker alkaline phosphatase staining was also observed in cell culture on day 7 and in limb digits on embryonic day 17.5, indicating mildly delayed ossification. Further characterization of these candidate Col10a1 transcriptional regulators could help identify novel therapeutic targets for skeletal diseases associated with abnormal chondrocyte hypertrophy.

Diversity and Evolution of Mineralized Skeletal Tissues in Chondrichthyans

The diversity of skeletal tissues in extant vertebrates includes mineralized and unmineralized structures made of bone, cartilage, or tissues of intermediate nature. This variability, together with the diverse nature of skeletal tissues in fossil species question the origin of skeletonization in early vertebrates. In particular, the study of skeletal tissues in cartilaginous fishes is currently mostly restrained to tessellated cartilage, a derived form of mineralized cartilage that evolved at the origin of this group. In this work, we describe the architectural and histological diversity of neural arch mineralization in cartilaginous fishes. The observed variations in the architecture include tessellated cartilage, with or without more massive sites of mineralization, and continuously mineralized neural arches devoid of tesserae. The histology of these various architectures always includes globular mineralization that takes place in the cartilaginous matrix. In many instances, the mineralized structures also include a fibrous component that seems to emerge from the perichondrium and they may display intermediate features, ranging from partly cartilaginous to mostly fibrous matrix, similar to fibrocartilage. Among these perichondrial mineralized tissues is also found, in few species, a lamellar arrangement of the mineralized extracellular matrix. The evolution of the mineralized tissues in cartilaginous fishes is discussed in light of current knowledge of their phylogenetic relationships.

Evolutionary developmental genetics of teeth and odontodes in jawed vertebrates: a perspective from the study of elasmobranchs

Most extant vertebrates display a high variety of tooth and tooth-like organs (odontodes) that vary in shape, position over the body and nature of composing tissues. The development of these structures is known to involve similar genetic cascades and teeth and odontodes are believed to share a common evolutionary history. Gene expression patterns have previously been compared between mammalian and teleost tooth development but we highlight how the comparative framework was not always properly defined to deal with different tooth types or tooth developmental stages. Larger-scale comparative analyses also included cartilaginous fishes: sharks display oral teeth and dermal scales for which the gene expression during development started to be investigated in the small-spotted catshark Scyliorhinus canicula during the past decade. We report several descriptive approaches to analyse the embryonic tooth and caudal scale gene expressions in S. canicula. We compare these expressions wih the ones reported in mouse molars and teleost oral and pharyngeal teeth and highlight contributions and biases that arise from these interspecific comparisons. We finally discuss the evolutionary processes that can explain the observed intra and interspecific similarities and divergences in the genetic cascades involved in tooth and odontode development in jawed vertebrates.

The multiscale architecture of tessellated cartilage and its relation to function

When describing the architecture and ultrastructure of animal skeletons, introductory biology, anatomy and histology textbooks typically focus on the few bone and cartilage types prevalent in humans. In reality, cartilage and bone are far more diverse in the animal kingdom, particularly within fishes (Chondrichthyes and Actinopterygii), where cartilage and bone types are characterized by features that are anomalous or even pathological in human skeletons. This review discusses the curious and complex architectures of shark and ray tessellated cartilage, highlighting similarities and differences with their mammalian skeletal tissue counterparts. By synthesizing older anatomical literature with recent high-resolution structural and materials characterization work, this review frames emerging pictures of form-function relationships in this tissue and of the evolution and true diversity of cartilage and bone.

Evolution of Matrix Gla and Bone Gla Protein Genes in Jawed Vertebrates

Matrix Gla protein (Mgp) and bone Gla protein (Bgp) are vitamin-K dependent proteins that bind calcium in their γ-carboxylated versions in mammals. They are recognized as positive (Bgp) or negative (Mgp and Bgp) regulators of biomineralization in a number of tissues, including skeletal tissues of bony vertebrates. The Mgp/Bgp gene family is poorly known in cartilaginous fishes, which precludes the understanding of the evolution of the biomineralization toolkit at the emergence of jawed vertebrates. Here we took advantage of recently released genomic and transcriptomic data in cartilaginous fishes and described the genomic loci and gene expression patterns of the Mgp/Bgp gene family. We identified three genes, Mgp1, Mgp2, and Bgp, in cartilaginous fishes instead of the single previously reported Mgp gene. We describe their genomic loci, resulting in a dynamic evolutionary scenario for this gene family including several events of local (tandem) duplications, but also of translocation events, along jawed vertebrate evolution. We describe the expression patterns of Mgp1, Mgp2, and Bgp in embryonic stages covering organogenesis in the small-spotted catshark Scyliorhinus canicula and present a comparative analysis with Mgp/Bgp family members previously described in bony vertebrates, highlighting ancestral features such as early embryonic, soft tissues, and neuronal expressions, but also derived features of cartilaginous fishes such as expression in fin supporting fibers. Our results support an ancestral function of Mgp in skeletal mineralization and a later derived function of Bgp in skeletal development that may be related to the divergence of bony vertebrates.

Laminin-Modified Dental Pulp Extracellular Matrix for Dental Pulp Regeneration

Native dental pulp extracellular matrix (DPEM) has proven to be an effective biomaterial for dental pulp regeneration. However, as a significant extracellular matrix glycoprotein, partial laminins were lost during the decellularization process, which were essential for odontoblast differentiation. Thereby, this study investigated the feasibility of LN supplementation to improve the surface of DPEM for odontoblast layer regeneration. The influences of laminin on cell adhesion and odontogenic differentiation were evaluated in vitro. Then, we fabricated laminin-modified DPEM based on the physical coating strategy and observed the location and persistency of laminin coating by immunofluorescent staining. Finally, laminin-modified DPEM combined with treated dentin matrix (TDM) was transplanted in orthotopic jaw bone of beagles (n = 3) to assess the effect of LNs on dental pulp tissue regeneration. The in vitro results showed that laminins could improve the adhesion of dental pulp stem cells (DPSCs) and promoted DPSCs toward odontogenic differentiation. Continuous odontoblastic layer-like structure was observed in laminin-modified DPEM group, expressing the markers for odontoblastogenesis, dentine matrix protein-1 (DMP-1) and dentin sialophosphoprotein (DSPP). Overall, these studies demonstrate that the supplementation of laminins to DPEM contributes to the odontogenic differentiation of cells and to the formation of odontoblast layer in dental pulp regeneration.

Endoskeletal mineralization in chimaera and a comparative guide to tessellated cartilage in chondrichthyan fishes (sharks, rays and chimaera)

An accepted uniting character of modern cartilaginous fishes (sharks, rays, chimaera) is the presence of a mineralized, skeletal crust, tiled by numerous minute plates called tesserae. Tesserae have, however, never been demonstrated in modern chimaera and it is debated whether the skeleton mineralizes at all. We show for the first time that tessellated cartilage was not lost in chimaera, as has been previously postulated, and is in many ways similar to that of sharks and rays. Tesserae in Chimaera monstrosa are less regular in shape and size in comparison to the general scheme of polygonal tesserae in sharks and rays, yet share several features with them. For example, Chimaera tesserae, like those of elasmobranchs, possess both intertesseral joints (unmineralized regions, where fibrous tissue links adjacent tesserae) and recurring patterns of local mineral density variation (e.g. Liesegang lines, hypermineralized ‘spokes’), reflecting periodic accretion of mineral at tesseral edges as tesserae grow. Chimaera monstrosa's tesserae, however, appear to lack the internal cell networks that characterize tesserae in elasmobranchs, indicating fundamental differences among chondrichthyan groups in how calcification is controlled. By compiling and comparing recent ultrastructure data on tesserae, we also provide a synthesized, up-to-date and comparative glossary on tessellated cartilage, as well as a perspective on the current state of research into the topic, offering benchmark context for future research into modern and extinct vertebrate skeletal tissues.

Adult chondrogenesis and spontaneous cartilage repair in the skate, Leucoraja erinacea

Mammalian articular cartilage is an avascular tissue with poor capacity for spontaneous repair. Here, we show that embryonic development of cartilage in the skate (Leucoraja erinacea) mirrors that of mammals, with developing chondrocytes co-expressing genes encoding the transcription factors Sox5, Sox6 and Sox9. However, in skate, transcriptional features of developing cartilage persist into adulthood, both in peripheral chondrocytes and in cells of the fibrous perichondrium that ensheaths the skeleton. Using pulse-chase label retention experiments and multiplexed in situ hybridization, we identify a population of cycling Sox5/6/9+ perichondral progenitor cells that generate new cartilage during adult growth, and we show that persistence of chondrogenesis in adult skates correlates with ability to spontaneously repair cartilage injuries. Skates therefore offer a unique model for adult chondrogenesis and cartilage repair and may serve as inspiration for novel cell-based therapies for skeletal pathologies, such as osteoarthritis.

Co-aligned chondrocytes: Zonal morphological variation and structured arrangement of cell lacunae in tessellated cartilage

In most vertebrates the embryonic cartilaginous skeleton is replaced by bone during development. During this process, cartilage cells (chondrocytes) mineralize the extracellular matrix and undergo apoptosis, giving way to bone cells (osteocytes). In contrast, sharks and rays (elasmobranchs) have cartilaginous skeletons throughout life, where only the surface mineralizes, forming a layer of tiles (tesserae). Elasmobranch chondrocytes, unlike those of other vertebrates, survive cartilage mineralization and are maintained alive in spaces (lacunae) within tesserae. However, the functions of the chondrocytes in the mineralized tissue remain unknown. Applying a custom analysis workflow to high-resolution synchrotron microCT scans of tesserae, we characterize the morphologies and arrangements of stingray chondrocyte lacunae, using lacunar morphology as a proxy for chondrocyte morphology. We show that the cell density is comparable in unmineralized and mineralized tissue and that cells maintain similar volume even when they have been incorporated into tesserae. Our findings support previous hypotheses that elasmobranch chondrocytes, unlike those of other taxa, do not proliferate, hypertrophy or undergo apoptosis during mineralization. Tessera lacunae show zonal variation in their shapes, being flatter further from and more spherical closer to the unmineralized cartilage matrix, and larger in the center of tesserae. The lacunae show pronounced organization into parallel layers and strong orientation toward neighboring tesserae. Tesserae also exhibit local variation in lacunar density, with the density considerably higher near pores passing through the tesseral layer, suggesting pores and cells interact, and that pores may contain a nutrient source. We propose that the different lacunar types reflect the stages of the tesserae formation process, while also representing local variation in tissue architecture and cell function. Lacunae are linked by small passages (canaliculi) in the matrix to form elongated series at the tesseral periphery and tight clusters in the center of tesserae, creating a rich connectivity among cells. The network arrangement and the shape variation of chondrocytes in tesserae indicate that cells may interact within and between tesserae and manage mineralization differently from chondrocytes in other vertebrates, perhaps performing analogous roles to osteocytes in bone.