- Research article
- Open Access
Distinct patterns of notochord mineralization in zebrafish coincide with the localization of Osteocalcin isoform 1 during early vertebral centra formation
BMC Developmental Biologyvolume 12, Article number: 28 (2012)
In chondrichthyans, basal osteichthyans and tetrapods, vertebral bodies have cartilaginous anlagen that subsequently mineralize (chondrichthyans) or ossify (osteichthyans). Chondrocytes that form the vertebral centra derive from somites. In teleost fish, vertebral centrum formation starts in the absence of cartilage, through direct mineralization of the notochord sheath. In a second step, the notochord is surrounded by somite-derived intramembranous bone. In several small teleost species, including zebrafish (Danio rerio), even haemal and neural arches form directly as intramembranous bone and only modified caudalmost arches remain cartilaginous. This study compares initial patterns of mineralization in different regions of the vertebral column in zebrafish. We ask if the absence or presence of cartilaginous arches influences the pattern of notochord sheath mineralization.
To reveal which cells are involved in mineralization of the notochord sheath we identify proliferating cells, we trace mineralization on the histological level and we analyze cell ultrastructure by TEM. Moreover, we localize proteins and genes that are typically expressed by skeletogenic cells such as Collagen type II, Alkaline phosphatase (ALP) and Osteocalcin (Oc). Mineralization of abdominal and caudal vertebrae starts with a complete ring within the notochord sheath and prior to the formation of the bony arches. In contrast, notochord mineralization of caudal fin centra starts with a broad ventral mineral deposition, associated with the bases of the modified cartilaginous arches. Similar, arch-related, patterns of mineralization occur in teleosts that maintain cartilaginous arches throughout the spine.
Throughout the entire vertebral column, we were able to co-localize ALP-positive signal with chordacentrum mineralization sites, as well as Collagen II and Oc protein accumulation in the mineralizing notochord sheath. In the caudal fin region, ALP and Oc signals were clearly produced both by the notochord epithelium and cells outside the notochord, the cartilaginous arches. Based on immunostaining, real time PCR and oc2:gfp transgenic fish, we identify Oc in the mineralizing notochord sheath as osteocalcin isoform 1 (Oc1).
If notochord mineralization occurs prior to arch formation, mineralization of the notochord sheath is ring-shaped. If notochord mineralization occurs after cartilaginous arch formation, mineralization of the notochord sheath starts at the insertion point of the arches, with a basiventral origin. The presence of ALP and Oc1, not only in cells outside the notochord, but also in the notochord epithelium, suggests an active role of the notochord in the mineralization process. The same may apply to Col II-positive chondrocytes of the caudalmost haemal arches that show ALP activity and Oc1 accumulation, since these chondrocytes do not mineralize their own cartilage matrix. Even without cartilaginous preformed vertebral centra, the cartilaginous arches may have an inductive role in vertebral centrum formation, possibly contributing to the distinct mineralization patterns of zebrafish vertebral column and caudal fin vertebral fusion.
The structural units of the vertebral column are the vertebrae, composed of neural and haemal arches and the vertebral body itself, the centrum. Vertebral bodies are joined by intervertebral tissue primarily derived from the notochord [1–3].
The vertebral column results from a strictly controlled segmentation process that occurs in all vertebrate species and is associated to two main structures, the notochord and the somites [4, 5]. The somites are epithelialized spheres of mesoderm that develop on either side of the neural tube, give rise to dermis, skeletal musculature (dermomyotome) and vertebrae (sclerotome) [6, 7]. While the somite contribution to vertebra formation has been extensively studied, particularly in birds [6, 8], the role of the notochord has received less attention. In teleosts, the majority of extant vertebrates, the notochord is composed of a core of large, vacuolated chordocytes, and an epithelial layer of chordoblasts that secrete the notochord sheath . The notochord sheath is a stratified structure, composed of a thin external membrane, with high elastin content, covering a thicker collagenous layer [10, 11].
In chondrichthyans and most osteichthyans, including tetrapods, but not in teleosts, vertebral bodies have a cartilaginous anlage that subsequently either mineralizes or is replaced by bone. In teleosts such as zebrafish, vertebral centra are formed in the absence of cartilage [12, 13]. Indeed, teleost vertebral centra form through the mineralization of the notochord sheath (chordacentrum), which is then surrounded by somite-derived intramembranous bone (autocentrum) [2, 12, 14–16].
In teleosts, such as Atlantic salmon (Salmo salar) or zebrafish (Danio rerio), the notochord plays an important role in early life stages, as its mechanical function is only replaced by the vertebral column in the postembryonic life [2, 17, 18]. In Atlantic salmon, the initial mineralization of the chordacentrum has been described to be associated with cells of the notochord epithelium (chordoblasts), while bone formation, by sclerotome-derived cells (autocentrum), is a second step . In contrast, Inohaya and co-workers suggest that in medaka (Oryzias latipes), only sclerotome-derived cells are involved in chordacentrum and autocentrum mineralization, with no role of chordoblasts . Yet, a recent study  shows that, also in medaka, with conditional ablation of osterix-positive osteoblasts, notochord sheath mineralization is maintained. Therefore, the main cellular and molecular determinants involved in early notochord mineralization are still under debate. Furthermore, the role of extracellular matrix proteins such as osteocalcin, generally expressed by mature and resting osteoblasts  and by hypertrophic chondrocytes [17, 22–24], in that process, remains unclear. In several teleosts, including zebrafish, two Osteocalcin genes (Oc1 and Oc2) have been identified .
In amniotes the vertebral column is divided into five main regions whereas the vertebral column of teleosts is often only subdivided into two main regions, abdominal and caudal (e.g., [26–28]). However, also in zebrafish, several regions can be recognized within the vertebral column. In particular, regions that contain the most anterior and the most posterior vertebrae are highly specialized [15, 29]. Regional differences are not only apparent at the morphological level but also regarding the tendency of vertebrae to fuse. While zebrafish vertebral bodies usually display no pathological fusion , caudal fin vertebrae undergo several fusions as part of regular development [15, 30–32]. Yet, other teleosts, such as Atlantic salmon and other farmed species, are known to suffer frequent pathological vertebral fusions [33, 34]. Whether regional differences in vertebrae morphology and mineralization relate to the susceptibility to fuse remains an open question.
This study aims to characterize mineralization patterns in different regions of the vertebral column (abdominal, caudal, caudal fin region) using various methods to reveal mineral deposition. Subsequently, these patterns are compared with the histogenesis of the arches, revealed through Collagen type II immunostaining, and to the timing of centrum formation. We characterize the proliferation of notochord cells and we also localize proteins related to mineralization, such as Alkaline phosphatase and Osteocalcin. We here provide the first evidence for the early presence of Osteocalcin 1 in mineralizing chordacentra. Finally, we discuss a possible association between timing of centra formation, the mineralization pattern and occurrence of vertebral fusion.
Centrum mineralization follows arch formation in the caudal fin region
The fully developed zebrafish vertebral column displays a mode of 33 vertebral bodies (Figure 1) including 4 centra in the Weberian apparatus, 10 abdominal vertebrae, 15 caudal vertebrae, and 3 to 4 caudal fin vertebrae (depending on the presence of a fourth preural) .
The first centra to mineralize (Figure 2a) are 3 and 4 (4.0 mm TL) in an anterior-to-posterior direction, followed by the mineralization of centrum 5, and later centra 1 and 2, in specimens of 4.4 mm TL. The mineralization of the remaining abdominal and caudal centra follows an anterior-to-posterior direction (Figure 2b). The mineralization of the caudal fin centra starts with the U1 from the compound centrum [PU1++U1], in specimens of 5.5 mm TL (Figure 2b). PU2 is the last centrum to be formed, in specimens of 6.7 mm TL, preceded by PU3 formation (Figure 2c). By then, the complex cartilaginous structure that forms the Weberian apparatus is already present in the anterior part of the vertebral column.
All Weberian, abdominal and caudal centra mineralize before development of the arches (Figure 3a). In the Weberian region, the haemal and neural arches develop as cartilaginous anlagen. These are Collagen type II positive, except for the associated spines that have no cartilaginous precursor. In the abdominal and caudal regions, the arches develop as intramembranous bones, with no cartilaginous anlage.
In contrast, in the modified caudal fin vertebrae, centrum mineralization occurs after arch formation. The (modified) arches associated to these centra have a cartilaginous anlage, also positive for Collagen type II. Staining of these arches was particularly strong at the bases (attachment areas to the notochord sheath) and distal parts, where no perichondral mineralization is present throughout the individual’s life span. Collagen type II immunostaining also showed protein accumulation in the notochord sheath, although less evident than in the cartilage (Figure 3a’).
In addition to the different timing of centrum mineralization versus arch formation, the present results also show a clear pattern concerning mineral expansion within the centra (Figure 3b, b’, c, c’).
With the exception of C1 and C2, vertebral bodies in the Weberian, abdominal and caudal regions start to mineralize in the form of a ring-shaped mineralized structure that expands in both anterior and posterior directions (Figure 3b and b’). These centra will be referred to as “ring centra”. In none of these vertebral bodies, have arches ever been observed prior to centrum mineralization. Nevertheless, the origin of mineralization matches the position of the myosepta, i.e. the site where neural and haemal arches will develop at a later stage. From this point, centrum mineralization expands one quarter anteriorly and three quarters posteriorly, showing clear incremental mineralization fronts (Figure 3b’).
Different from this pattern, in the caudal fin (PU4/3 to U2) region, mineralization starts as a broad ventral mineral deposition (basiventral origin) (Figure 3c and c’), which expands dorsally. Here, mineralization starts and expands from the attachment sites of the already developed modified cartilaginous haemal arches to the notochord. The developing centra show no apparent growth fronts as seen in the “ring centra” (Figure 3c’).
Inner and outer cell distribution during chorda- and autocentrum formation
Irrespective of the shape of the mineralized centrum, (“ring” and caudal fin centra) (Figure 4a, b), the origin of mineralization was shown both with von Kossa staining (Figure 4c, c’, d, d’) and TEM (Figure 4e-f), to occur always within the notochord sheath, thus establishing the chordacentrum. Once the chordacentrum is fully formed, perichordal bone is deposited mostly at the anterior and posterior edges, in this way establishing the vertebral endplates of the centra, thus forming the autocentrum (Figure 5a-d). Autocentrum formation occurs irrespective of the initial mineralization pattern within the notochord sheath (“ring” vs. basiventral origin). Yet, in the caudal fin centra, decreased bone formation is observed in the areas where cartilaginous arches attach to the notochord (Figure 5c-d).
During chordacentrum formation, notochord inner cells or chordoblasts do not accumulate in the prospective intervertebral space (Figure 4a, b), unlike what has been observed in other teleost species. However, once the centrum is mineralized, a higher number of chordoblasts can be seen in the non-mineralized intervertebral space (Figure 5a, c, e). By then, the chordoblasts are oblong, larger in size and are oriented in a dorsal-ventral direction (Figure 5e), as opposed to the squamous morphology displayed prior to and during chordacentrum development. In addition, chordoblasts proliferate in the intervertebral space, as shown by PCNA staining. Here, notochord sheath thickening occurs (Figure 5f). This pattern of chordoblast distribution is observed throughout the notochord, regardless of the type of centrum and its mineralization pattern (Figure 5a-d).
During chordacentrum formation, contrary to the inner notochord cells, the assumed sclerotome-derived notochord outer cells  in the “ring” vertebrae are visible as few scattered nuclei located around the notochord (Figure 4e). These cells are squamous and without apparent morphological features of active osteoblasts. In the caudal fin region, the cartilaginous haemal arches of the PUs (parhypural and hypurals 1, 2 and 4) directly attach to the notochord sheath with no space for other cells in-between the notochord sheath and the chondrocytes of the arches (Figure 4f).
During autocentrum formation, osteoblasts can be clearly identified in association with bone deposited around the “ring” chordacentra (Figure 5a). Major areas of the caudal fin chordacentra are covered by the cartilaginous arches attached to the notochord sheath surface (Figure 5c). Yet some osteoblasts can be identified wherever the arches are absent, lining the edges of the prospective vertebral endplates.
Osteocalcin in chordacentrum formation
In order to clarify possible factors in chordacentrum mineralization, both in “ring” and caudal fin centrum, we have localized Alkaline phosphatase (ALP) and Osteocalcin (Figure 6).
In “ring” centra, ALP-positive signal is co-localized with the site of chordacentrum mineralization prior to any bone deposition outside the notochord sheath. The flattened shape of both outer and inner notochord cells makes it difficult to clearly distinguish the origin of the ALP signal, due to its membrane-bound characteristic (Figure 6a). In caudal fin centra, the arch chondrocytes, particularly those adjacent to the mineralizing notochord sheath, co-localize with a strong ALP signal (Figure 6b).
For Osteocalcin immunostaining, we used an antibody raised against Argyrosomus regius Oc1, which was previously validated . It was considered to detect a single protein in zebrafish, since a western blot analysis gave a single band of approximately 5kDa. Nevertheless we cannot guarantee that it distinguishes zebrafish Oc1 from Oc2. However, since Oc2 appears only after day 7 post-fertilization in zebrafish as shown by qPCR analysis (Figure 6l, m), immunostaining in younger fish, such as during formation of chordacentra C1 to C5 (Figure 2a) could only represent sites of Oc1 protein accumulation.
In “ring” (Figure 6c) and caudal fin centra (Figure 6d), Oc1 protein accumulation is observed in the mineralizing matrix of the notochord sheath, even at very early stages of chordacentrum formation. Oc1 accumulation co-localizes with chordacentrum mineralization, reflecting the exact same pattern as observed with alizarin red staining (Figure 3b, b’, c, c’). Moreover, in later stages of centrum formation, we can detect progression fronts in “ring” centra (Figure 6e) while in caudal fin centra we observe homogenous protein accumulation (Figure 6f). In the caudal fin region, Oc1 is also detected in the vertebral arches and caudal fin rays (Figure 6d, f), which are examples of perichondral and intramembranous mineralization, respectively.
Analysis of semi-thin sections of the immunostained samples confirms accumulation of Osteocalcin in the notochord sheath and not around the notochord, during early (Figure 6g, g’, h, h’) and late (Figure 6i, j) stages of chordacentrum formation. In caudal fin centra (Figure 6h, j), the arch chondrocytes adjacent to the mineralizing notochord sheath also co-localize with Oc1 protein (Figure 6h’), as observed for ALP. No perichondral ossification occurred in the arch at sites positive for ALP and Oc1, suggesting that these signals are not associated to arch ossification.
Observation of the transgenic fish Tg(oc2:gfp; osx:mcherry) showed that no oc2 expression is present during early stages of centrum formation. In fact, oc2 is completely absent during chordacentrum formation, appearing first in intramembranous formation of abdominal and caudal arches, a process previously described as mediated by mature osteoblasts (Figure 6k, k’, k”). Analysis of relative oc1 and oc2 expression levels during zebrafish development shows that, while oc1 maternal transcripts are detected in early developmental stages, when no mature osteoblasts can yet be detected (Figure 6l), oc2 is only expressed after 7dpf (Figure 6m). Chordacentrum formation of the first vertebrae starts at 4.0 mm TL, as previously indicated, which, despite the natural variation in growth, corresponds to fish younger than 7 days.
Zebrafish and medaka vertebral column development and morphology have been extensively described (e.g.[1, 12, 19, 29, 31, 37]). Nevertheless, many questions remain unanswered. In particular, how do distinct mineralization patterns relate to regional differences in arch histogenesis and to susceptibility for vertebral fusion? Why do vertebrae in the caudal fin region fuse as part of the developmental process, in contrast to anterior regions where vertebrae remain individualized [15, 30]?
Here, we describe two distinct mineralization patterns of vertebral centra associated with specific vertebral regions. We characterize these patterns based on (1) how mineral is being deposited throughout the centrum, (2) timing of development (prior to or after arch formation) and (3) arch histogenesis (intramembranous or cartilaginous origin). While describing this process we localized accumulation of oc1 protein in mineralizing notochord sheath, suggesting it could play a function in early mineralization events during chordacentrum formation.
Notochord segmentation, arch contribution to centra and association with fusion
The first mineralization of the vertebral bodies, independent of the vertebral type (Weberian, abdominal, caudal and caudal fin), occurs within the notochord sheath, establishing the chordacentrum. These are the first indications of segmented units in the zebrafish notochord. Other authors have already described that early centrum mineralization in teleosts occurs within the notochord sheath [2, 12, 13, 19, 38, 39]. In addition, the notochord may actively contribute to segmentation and centrum formation [2, 12, 20].
Thus, while the chordacentrum is always established through mineralization of the notochord sheath, the pattern of mineralization differs between caudal fin and remaining vertebral regions: (1) a caudal fin centrum mineralizes after arch formation, while a “ring” centrum is always formed prior to haemal or neural arches, (2) the modified arches from caudal fin vertebrae are composed of cartilage while “ring” centra outside the Weberian apparatus are associated to arches formed by intramembranous ossification; and (3) caudal fin centra display a basiventral origin of mineralization, while all other vertebrae start with a “ring”-shaped mineralization.
Arch histogenesis appears to be a clear mark that distinguishes caudal fin vertebrae from the remaining regions. A similar, regional distinct pattern of mineralization has been observed in medaka, with cartilaginous arches in the caudal fin region and intramembranous arches in the anterior part of the vertebral column [19, 40]. Although notochord mineralization of anterior centra in medaka starts dorsal, subsequently chordacentrum mineralization will proceed as “ring” centra . Yet, “ring” notochord-based centrum formation, is not an unique feature of zebrafish and medaka, occurring also in other taxa . Considering caudal fin centra of zebrafish and medaka, both display similar patterns of mineralization , linked to the presence of modified cartilaginous arches. Many basal ray-finned fishes have cartilaginous arches directly attached to an unconstricted notochord  suggesting that cartilaginous arches are an ancestral feature. Therefore, the caudal fin endoskeleton of zebrafish, and also medaka, appear to display the evolutionarily conserved condition, with cartilaginous arch formation.
Mineralization in the caudal fin centra shows a clear association with attachment sites of the cartilaginous haemal arches to the notochord (basiventral mineralization pattern). The association between basiventral mineralization and presence of cartilaginous haemal arches has also been described in other teleost species, such as the Goldeye (Hiodon alosoides) , Atlantic salmon  and White seabream (Diplodus sargus) . In these species too, centrum mineralization originates exactly where the cartilaginous haemal arches are attached to the notochord. Mineralization then progresses from a basiventral origin through bilateral wedge-shaped areas that meet dorsally.
Moreover, we have identified mineralization related proteins in the chondrocytes of the cartilaginous arches, in the immediate vicinity of the area of the incipient mineralization of the notochord. This supports the possible role of the cartilaginous arches in early centrum establishment, a hypothesis that will be discussed in more detail in the next section.
The cartilaginous arches appear to be strictly associated not only with chordacentrum formation but also with the start of autocentrum formation, given that the arches are directly attached to the notochord, affecting bone deposition around the notochord.
The type of arches and associated pattern of centrum mineralization appears to be correlated not only to the formation of the vertebral bodies but also to their ability to fuse or stay as individual units (Figure 7). In the medaka bis (biaxial symmetries) mutant, centra mineralize at sites where ectopic hypurals attach to the notochord, leading to the development of fused centra . Thus, it is tempting to suggest from the zebrafish and medaka data, that the cartilaginous arches have a role in centrum mineralization. Perhaps the predisposition to develop vertebral fusions is related to a specific mineralization pattern with basiventral origin (cartilaginous arch associated).
Distribution of notochord inner and outer cells during centrum formation
In Atlantic salmon, notochord sheath mineralization is preceded by an apparent segmentation of the chordoblast layer . The same holds for medaka (PEW personal observations). In contrast, differential chordoblast distribution in zebrafish is only clearly observed after chordacentrum mineralization, with few scattered squamous cells in the mineralized area and large, dorso-ventrally oriented, highly proliferative cells in the intervertebral region.
Regarding notochord outer cells, it appears that the segmental pattern of centra in the caudal fin region is associated with the anatomical pattern of arches. In contrast, no association can be found in “ring” centra between outer cell distribution and early chordacentrum formation. The latter is in agreement with the data from zebrafish fused somite mutants (fss), where paraxial mesoderm lacks a proper segmentation . This mutant shows fusion of vertebral arches, while vertebral centra are individualized, suggesting that centrum segmentation does not rely on sclerotome patterning only. The resegmentation process in zebrafish has been suggested to follow non-lineage restricted compartments, where one somite contributes to more than one vertebral body, independently of the antero-posterior somitic domains (leaky-resegmentation) .
Although mineralization is not preceded by a clear pattern of inner or outer cell distribution (exceptions are arch-chondrocytes in the caudal fin), these cells may nevertheless contribute to centrum mineralization. In Atlantic salmon, initial centrum mineralization has been related to chordoblasts, followed by sclerotome-derived bone formation [2, 10]. Different from Atlantic salmon, in medaka sclerotome-derived osteoblasts, with no input of chordoblasts, have been described to form vertebral centra [19, 35]. Yet, recent data show that conditional ablation of osteoblasts in medaka maintains notochord mineralization and even leads to vertebral fusion . This suggested that notochord cells have the ability to induce mineralization of the notochord sheath. This is consistent with studies in zebrafish [12, 45] that show that notochord cell ablation prevents early centrum mineralization.
Osteocalcin and chordacentrum formation
Osteocalcin is the most abundant noncollagenous protein in bone of vertebrates, including all teleosts analyzed so far . Although Osteocalcin was already described in the vertebral bodies of teleosts including zebrafish [47–50], no previous reports demonstrated Oc in early notochordal mineralization events. Thus, the presence of Oc in stages of chordacentrum mineralization has to be considered in a perspective different from that of regular bone formation.
Within teleosts, Osteocalcin is represented by two distinct isoforms (Oc1 and Oc2) [25, 51], encoded by two different genes, which were proposed to originate from a duplication event . Analysis of relative gene expression levels by qPCR for oc1 and oc2 showed that oc2 expression is first detected after 7dpf, while oc1 is maternally transcribed, suggesting a different role of the two isoforms. In fact, oc1 is detected prior to any osteoblast differentiation or mineralization event. This leads to the hypothesis that oc1 must be involved in processes prior to bone formation. Detection of Oc1 by immunostaining shows that this protein accumulates exactly where mineralization of the chordacentrum occurs. In the “ring” centra, the origin of the secreted protein is difficult to assess, due to the squamous morphology of the cells of the notochord epithelium and the inconspicuous presence of cells outside the notochord. In the caudal fin region, ALP positive chondrocytes, directly adjoining the mineralizing notochord sheath, are positive for Oc immunostaining, as well as notochord inner cells. In the Tg(oc2:gfp; osx:mcherry) fish, oc2 expression was located at the bone of arches of the “ring” centra, that form via intramembranous bone formation. This time point corresponds to the phase of autocentrum mineralization, as previously described in medaka . Therefore, our results suggest that Oc1, but not Oc2, is connected to chordacentrum mineralization. Likewise, in larval and juvenile Atlantic cod , oc1 and oc2 are expressed in different cells associated with mineralizing structures. In Atlantic salmon, some mineralized areas with reduced expression of osteocalcin (oc1) were considered to be associated with expression of a second isoform that the authors did not identify . Laizé and co-workers  described Oc2 propeptide as serine-rich, potentially phosphorylated and containing acid residues. These characteristics allow this protein to bind numerous calcium ions. If Osteocalcin has an active role in the process of mineralization or if it is a bone-related hormone that passively binds to calcium [54, 55], is a current debate and must be further explored.
In conclusion, we have analyzed distinct mineralization patterns in different regions of the zebrafish vertebral column and related them to the timing and nature of vertebral arch formation. We have examined the distribution of cells within and outside the notochord, and critically assessed their potential contribution to chordacentrum and autocentrum formation. If notochord mineralization occurs prior to arch formation, the mineralization of the notochord sheath is ring-shaped. If notochord mineralization occurs after cartilaginous arch formation, mineralization of the notochord sheath starts at the insertion point of the arches, with a basiventral origin. Finally, we have detected Oc1 in the notochord sheath during chordacentrum formation, whereas oc2 was only expressed later in association with bone formation. These results emphasize the need to better characterize the function of both osteocalcin genes in zebrafish and throughout the vertebrate clade.
Ethics statement on animal experiments
Animal handling and experiments were legally accredited by the Portuguese Direcção Geral de Veterinária (DGV) and all the experimental procedures involving animals followed the EU (Directive 86/609/CEE) and National (Directives 1005/92 from October 23, 466/95 from May 17 and 1 1131/97 from November 7) legislation for animal experimentation and welfare.
Sampling and wholemount skeletal analysis
To establish a timeline of events for early vertebrae development, an ontogenic series of Danio rerio was studied. Zebrafish eggs were obtained from natural spawning of wild-type fish and larvae were maintained and reared at standard conditions . Fish were collected at size intervals of 0.1 mm, between 4.0 and 8.0 mm of total length (TL). Similar-sized specimens were always selected among age-matched fish in order to avoid additional interspecific variation.
Larvae were euthanized with a lethal dose of MS222 (Sigma-Aldrich, St. Louis, MO) and fixed for 24 h in 4% buffered paraformaldehyde at 4ºC.
Early mineralization stages of the vertebral bodies were observed following 10 minutes staining with 0.01% alizarin red S in 70% ethanol solution. Specimens were observed under a Zeiss Axio Imager microscope, through fluorescence detection .
Light and transmission electron microscopy (TEM)
Zebrafish specimens of 5.0 and 11.0 mm TL with no pre-staining were fixed in a mixture of 1.5% glutaraldehyde and 1.5% paraformaldehyde in 0.1 M cacodylate buffer and processed for embedding in Epon, according to procedures previously described . Parasagittal 1 μm semi-thin sections of the sites of interest were stained with toluidine blue for 1–2 min (0.2% toluidine blue, 2% Na2CO3), rinsed with water, air-dried and mounted with DPX (Fluka, Buchs, Switzerland). Ultrathin sections were contrasted with uranyl acetate and lead citrate and observed with a Jeol JEM 1010 (Jeol, Tokyo, Japan), operating at 60 kV. Images were digitized using a DITABIS system (Pforzheim, Germany).
Identification of sites of mineral deposition and Alkaline phosphatase activity (ALP)
To determine sites of mineral deposition and of Alkaline phosphatase (ALP) activity during centrum formation, specimens ranging from 5.0 to 6.5 mm TL were dehydrated and embedded in glycol methacrylate, as previously described , and sectioned at 5 μm.
For detection of minerals, a silver nitrate coupling method, according to von Kossa , was used. Sections were counterstained with neutral red. Slides were examined for the appearance of optimal signal, rinsed in distilled water and mounted with DPX (Sigma-Aldrich, St. Louis, MO).
Alkaline phosphatase activity was detected using a previously described method  with some modification. Briefly, slides were incubated in staining buffer (100 mM Tris/HCl, pH 9.5; 50mM MgCl2; 100 mM NaCl) for 10 min. ALP was demonstrated by adding 4-nitro blue tetrazolium chloride (NBT; 0.38 mg ml-1; Roche) and 5-bromo-4-chloro-3-indolyl phosphate (BCIP; 0.175 mg ml-1; Roche) to the buffer solution, for a maximum of 30 min at 37ºC, in the dark. Negative control slides were incubated for 10 min at 85ºC in staining buffer prior to addition of NBT and BCIP. Slides were examined for the appearance of optimal signal-noise ratio, counterstained with Nuclear Fast Red Solution (Sigma-Aldrich, St. Louis, MO) and mounted in DPX (Sigma-Aldrich, St. Louis, MO).
Immunolocalization of PCNA, Osteocalcin and Collagen type II
For detection of cell proliferation through Proliferating cell nuclear antigen (PCNA) immunolocalization, 6.2 mm TL fish were paraffin-embedded according to routine procedures and serially sectioned at 7 μm. A mouse monoclonal primary antibody produced against human antigen (Sigma-Aldrich, St. Louis, MO, dilution 1:1000) was used  together with a goat polyclonal anti-mouse immunoglobulin biotin-coupled secondary antibody (Dako, dilution 1:500).
Immunohistochemical detection of Osteocalcin and Collagen type II was performed on fish of 4 to 6.5 mm TL as wholemount staining, according to a procedure previously described , using an anti-meagre (Argyrosomus regius) Osteocalcin rabbit polyclonal primary antibody, previously validated for zebrafish . For Collagen type II an anti-chicken mouse antibody (II-II6B3, Hybridoma Bank) previously validated for zebrafish , was used. Control fish were treated with secondary antibody alone.
Osteocalcin stained specimens were subsequently dehydrated, embedded in epon, serially sectioned at 4 μm and mounted with DPX (Sigma-Aldrich, St. Louis, MO), to achieve a more detailed analysis. Wholemount fish stained for Collagen type II were observed using a Leica MZ Apo stereomicroscope. Sections were observed using a Zeiss Axio Imager Microscope and photographed using an Axiocam MRC videocamera.
The transgenic zebrafish line Tg(oc2:gfp; osx:mcherry), previously described by , was provided by Stefan Schulte-Merker from Hubrecht Institute-KNAW and University Medical Centre (Uppsalalaan 8, 3584 CT Utrecht, The Netherlands).
RNA extraction and oc1/oc2quantitative expression analysis
Total RNA was extracted from a pool of up to twenty specimens following the Chomczynski and Sacchi method . Samples were collected at stages 8-32 or 128 cells, at shield stage, at 12, 24, 48, 72 and 96 hours post-fertilization (hpf) and at 5, 7, 15, 30, 45 and 60 days post-fertilization (dpf). To prevent genomic DNA contamination, following RNA extraction, samples were DNase treated with RQ1 DNase (Promega, Madison, WI) according to the manufacturer’s protocol, followed by a phenol-chlorophorm purification step. RNA quantity and quality was determined by spectrophotometry (NanoDrop ND-1000, Thermo Scientific) and electrophoresis.
cDNA was prepared from 1μg of total RNA from each sample using Moloney-murine leukemia virus (M-MLV) reverse transcriptase (Invitrogen), RNase Out (Invitrogen), and an oligo(dT) adapter (5'-ACGCGTCGACCTCGAGATCGATG(T)13-3') in a total volume of 20μl for 1h at 37ºC. For quantitative real-time PCR a StepOnePlus 96 (Applied Biosystems) was used together with 1× SsoFastTM EvaGreen Supermix (BioRad, Richmond, USA), 0.2 μM of forward (Oc1: 5′-GAAGCGAACATGAAGAGTCTGACAGTCC-3′; Oc2: 5′-CCAACTCCGCATCAGACTCCGCATCA-3′) and reverse (Oc1- 5′-TTTATAGGCGGCGATGATTCC-3′; Oc2: 5′-AGCAACACTCCGCTTCAGCAGCACAT-3′) primers and 200 ng of reverse-transcribed RNA. PCR conditions were 15 min at 95°C followed by 40–50 amplification cycles (each cycle is 30 s at 95°C, 15 s at 68°C). Expression values were normalized using the Elongation factor 1α (Ef1α) housekeeping gene . Gel electrophoresis and melt curve analysis were used to confirm specific product formation.
All data were statistically analyzed using a non-parametric Kruskal-Wallis test, followed by Wilcoxon multiple comparison with the GraphPad Prism software.
The subdivision of the vertebral column in different regions is based on the nomenclatures of Arratia et al. , Bird and Mabee  and Nybelin  (Figure 1). From anterior to posterior, vertebrae were identified as (1) Weberian vertebrae, (2) abdominal vertebrae (rib-bearing with open haemal arches and without haemal spines), (3) caudal vertebrae (with haemal arches closed), and (4) caudal fin vertebrae (with modified haemal and neural arches and respective spines). The nomenclature of the Weberian apparatus follows . The caudal fin vertebrae nomenclature follows Patterson , adapted by Arratia and Schultze  and Bensimon-Brito et al. . The caudal fin vertebral region contains preural and ural vertebral bodies. Preural vertebral bodies (abbreviated as PU) have haemal and neural arches that support the caudal fin. Ural vertebral bodies (briefly urals, abbreviated as U) support the caudal fin with modified haemal arches that do not enclose the caudal artery. A plus sign is added in superscript when a vertebral centrum is derived from a fusion event (compound centrum), as described previously .
Morin-Kensicki EM, Melancon E, Eisen JS: Segmental relationship between somites and vertebral column in zebrafish. Development. 2002, 129: 3851-3386.
Grotmol S, Kryvi H, Nordvik K, Totland GK: Notochord segmentation may lay down the pathway for the development of the vertebral bodies in the Atlantic salmon. Anat Embryol. 2003, 207: 263-272. 10.1007/s00429-003-0349-y.
Turnpenny PD, Alman B, Cornier AS, Giampietro PF, Offiah A, Tassy O, Pourquié O, Kusumi K, Dunwoodie S: Abnormal vertebral segmentation and the notch signaling pathway in man. Dev Dyn. 2007, 236 (6): 1456-1474. 10.1002/dvdy.21182.
Hall BK: Evolutionary consequences of skeletal differentiation. Am Zool. 1975, 15: 329-350.
Kaplan KM, Spivak JM, Bendo JA: Embryology of the spine and associated congenital abnormalities. Spine J. 2005, 5 (5): 564-576. 10.1016/j.spinee.2004.10.044.
Christ B, Huang R, Wilting J: The development of the avian vertebral column. Anat Embryol. 2000, 202: 179-194. 10.1007/s004290000114.
Kahane N, Cinnamon Y, Kalcheim C: The origin and fate of pioneer myotomal cells in the avian embryo. Mech Dev. 1998, 74: 59-73. 10.1016/S0925-4773(98)00066-5.
Christ B, Huang R, Scaal M: Formation and differentiation of the avian sclerotome. Anat Embryol. 2004, 208: 333-350.
Stemple DL: Structure and function of the notochord: an essential organ for chordate development. Development. 2005, 132: 2503-2512. 10.1242/dev.01812.
Grotmol S, Nordvik K, Kryvi H, Totland GK: A segmental pattern of alkaline phosphatase activity within the notochord coincides with the initial formation of the vertebral bodies. J Anat. 2005, 206: 427-436. 10.1111/j.1469-7580.2005.00408.x.
Schmitz RJ: Comparative ultrastructure of the cellular components of the unconstricted notochord in the sturgeon and the lungfish. J Morphol. 1998, 236: 75-104. 10.1002/(SICI)1097-4687(199805)236:2<75::AID-JMOR1>3.0.CO;2-N.
Fleming A, Keynes R, Tannahill D: A central role for the notochord in vertebral patterning. Development. 2004, 131: 873-880. 10.1242/dev.00952.
Huxley TH: Observations on the development of some parts of the skeleton of fishes. Q J Microsc Sci. 1859, 7: 33-46.
François Y: Structure et développement de la vertèbre de Salmo et des téléostéens. Archives de Zoologie Experimentale et Generale. 1966, 107: 287-328.
Bensimon-Brito A, Cancela ML, Huysseune A, Witten PE: Vestiges, rudiments and fusion events: the zebrafish caudal fin endoskeleton in an evo-devo perspective. Evol Dev. 2012, 14 (1): 116-127.
Arratia G, Schultze H-P, Casciotta J: Vertebral column and associated elements in Dipnoans and comparison with other fishes: Development and homology. J Morphol. 2001, 250: 101-172. 10.1002/jmor.1062.
Gavaia PJ, Simes DC, Ortiz-Delgado JB, Viegas CSB, Pinto JP, Kelsh RN, Sarasquete MC, Cancela ML: Osteocalcin and matrix Gla protein in zebrafish (Danio rerio) and Senegal sole (Solea senegalensis): Comparative gene and protein expression during larval development through adulthood. Gene Expr Patterns. 2006, 6 (6): 637-652. 10.1016/j.modgep.2005.11.010.
Lauder GV: On the relationship of the myotome to the axial skeleton in vertebrate evolution. Paleobiology. 1980, 6 (1): 51-56.
Inohaya K, Takano Y, Kudo A: The Teleost intervertebral region acts as a growth center of the centrum: In vivo visualization of osteoblasts and their progenitors in transgenic fish. Dev Dyn. 2007, 236: 3031-3046. 10.1002/dvdy.21329.
Willems B, Büttner A, Huysseune A, Renn J, Witten PE, Winkler C: Conditional ablation of osteoblasts in medaka. Dev Biol. 2012, 364: 128-137. 10.1016/j.ydbio.2012.01.023.
Huysseune A, Takle H, Soenens M, Taerwe K, Witten PE: Unique and shared gene expression patterns in Atlantic salmon (Salmo salar) tooth development. Developmental Genes and Evolution. 2008, 218: 427-437. 10.1007/s00427-008-0237-9.
Pullig O, Weseloh G, Ronneberger D-L, Käkönen S-M, Swoboda B: Chondrocyte differentiation in Human osteoarthritis: Expression of osteocalcin in normal and osteoarthritic cartilage and bone. Calcif Tissue Int. 2000, 67: 230-240. 10.1007/s002230001108.
Aizawa T, Roach HI, Kokubun S, Tanaka Y: Changes in the expression of Fas, osteonectin and osteocalcin with age in the rabbit growth plate. The Journal of Bone and Joint Surgery (Br). 1998, 80-B (5): 880-887.
Neugebauer BM, Monroe MA, Broess M, Gerstenfeld LC, Hauschka PV: Characterization of structural sequences in the chicken osteocalcin gene: Expression of osteocalcin by maturing osteoblasts and by hypertrophic chondrocytes in vitro. J Bone Miner Res. 1995, 10 (1): 157-163.
Laizé V, Viegas CSB, Price PA, Cancela ML: Identification of an osteocalcin isoform in fish with a large acidic prodomain. J Biol Chem. 2006, 281 (22): 15037-15043. 10.1074/jbc.M600373200.
Fjelldal PG, Hansen T, Breck O, Ornsrud R, Lock EJ, Waagbø R, Wargelius A, Witten PE: Vertebral deformities in farmed Atlantic salmon (Salmo salar L.) – etiology and pathology. J Appl Ichthyol. 2012, 28 (3): 433-440. 10.1111/j.1439-0426.2012.01980.x.
Ward AB, Brainerd EL: Evolution of axial patterning in elongate fishes. Biol J Linn Soc. 2007, 90: 97-116. 10.1111/j.1095-8312.2007.00714.x.
Wargelius A, Fjelldal PG, Hansen T: Heat shock during early somitogenesis induces caudal vertebral column defects in Atlantic salmon (Salmo salar). Dev Genes Evol. 2005, 215 (7): 350-357. 10.1007/s00427-005-0482-0.
Bird NC, Hernandez LP: Building an evolutionary innovation: differential growth in the modified vertebral elements of the zebrafish Weberian apparatus. Zoology. 2009, 112 (2): 97-112. 10.1016/j.zool.2008.05.003.
Ferreri F, Nicolais C, Boglione C, Bertolini B: Skeletal characterization of wild and reared zebrafish: anomalies and meristic characters. J Fish Biol. 2000, 56 (5): 1115-1128. 10.1111/j.1095-8649.2000.tb02127.x.
Bird NC, Mabee PM: Developmental morphology of the axial skeleton of the zebrafish, Danio rerio (Ostariophysi: Cyprinidae). Dev Dyn. 2003, 228 (3): 337-357. 10.1002/dvdy.10387.
Bensimon-Brito A, Cancela ML, Huysseune A, Witten PE: The zebrafish (Danio rerio) caudal complex: a model to study vertebral body fusion. J Appl Ichthyol. 2010, 26 (2): 235-238. 10.1111/j.1439-0426.2010.01412.x.
Witten PE, Obach A, Huysseune A, Baeverfjord G: Vertebrae fusion in Atlantic salmon (Salmo salar): development, aggravation and pathways of containment. Aquaculture. 2006, 258 (1–4): 164-172.
Fjelldal PG, Nordgarden U, Wargelius A, Taranger GL, Waagbø R, Olsen RE: Effects of vegetable feed ingredients on bone health in Atlantic salmon. J Appl Ichthyol. 2010, 26 (2): 327-333. 10.1111/j.1439-0426.2010.01430.x.
Ekanayake S, Hall BK: The development of acellularity of the vertebral bone of the Japanese Medaka, Oryzias latipes (Teleostei; Cyprinidontidae). J Morphol. 1987, 193: 253-261. 10.1002/jmor.1051930304.
Simes DC, Williamson MK, Schaff BJ, Gavaia PJ, Ingleton PM, Price PA, Cancela ML: Characterization of osteocalcin (BGP) and matrix gla protein (MGP) fish specific antibodies: validation for immunodetection studies in lower vertebrates. Calcif Tissue Int. 2004, 74: 170-180. 10.1007/s00223-003-0079-4.
Du SJ, Frenkel V, Kindschi G, Zohar Y: Visualizing normal and defective bone development in Zebrafish embryos using the fluorescent chromophore calcein. Dev Biol. 2001, 238: 239-246. 10.1006/dbio.2001.0390.
Nordvik K, Kryvi H, Totland GK, Grotmol S: The salmon vertebral body develops through mineralization of two preformed tissues that are encompassed by two layers of bone. J Anat. 2005, 206: 103-114. 10.1111/j.1469-7580.2005.00372.x.
Kölliker A: Weitere Beobachtungen über die Wirbel der Selachier, insbesondere über die Wirbel der Lamnoidei, nebst allgemeinen Bemerkungen über die Bildung der Wirbel der Plagiostomen. 1863, Frankfurt: Verlag HL. Bröuner
Hibiya K, Katsumoto T, Kondo T, Kitabayashi I, Kudo A: Brpf1, a subunit of the MOZ histone acetyl transferase complex, maintains expression of anterior and posterior Hox genes for proper patterning of craniofacial and caudal skeletons. Dev Biol. 2009, 329: 176-190. 10.1016/j.ydbio.2009.02.021.
Laerm J: The origin and homology of the chondrostean vertebral centrum. Can J Zool. 1979, 57 (3): 475-485. 10.1139/z79-058.
Schultze H-P, Arratia G: Reevaluation of the caudal skeleton of some actinopterygian fishes: II. Hiodon, Elops, and Albula. J Morphol. 1988, 195 (3): 257-303. 10.1002/jmor.1051950304.
Koumoundouros G, Sfakianakis DG, Maingot E, Divanach P, Kentouri M: Osteological development of the vertebral column and of the fins in Diplodus sargus (Teleostei: Perciformes: Sparidae). Mar Biol. 2001, 139: 853-862. 10.1007/s002270100645.
van Eeden FJM, Granato M, Schach U, Brand M, Furutani-Seiki M, Haffter P, Hammerschmidt M, Heisenberg C-P, Jiang Y-J, Kane DA, et al: Mutations affecting somite formation and patterning in the zebrafish, Danio rerio. Development. 1996, 123: 153-164.
Fleming A, Keynes RJ, Tannahill D: The role of the notochord in vertebral column formation. J Anat. 2001, 199 (1–2): 177-180.
Nishimoto SK, Waite JH, Nishimoto M, Kriwacki RW: Structure, activity, and distribution of fish osteocalcin. J Biol Chem. 2003, 278 (14): 11843-11848. 10.1074/jbc.M211449200.
Simes DC, Williamson MK, Ortiz-Delgado JB, Viegas CSB, Price PA, Cancela ML: Purification of Matrix Gla Protein from a marine teleost fish, Argyrosomus regius: Calcified cartilage and not bone as the primary site of MGP accumulation in fish. J Bone Miner Res. 2003, 18 (2): 244-259. 10.1359/jbmr.2003.18.2.244.
Ytteborg E, Togersen J, Baeverfjord G, Takle H: Morphological and molecular characterization of developing vertebral fusions using a teleost model. BMC Physiol. 2010, 10 (13): 1-15.
Krossøy C, Ornsrud R, Wargelius A: Differential gene expression of bgp and mgp in trabecular and compact bone of Atlantic salmon (Salmo salar L.) verebrae. J Anat. 2009, 215: 663-672. 10.1111/j.1469-7580.2009.01153.x.
Pinto JP, Ohresser MCP, Cancela ML: Cloning of the bone Gla protein gene from the teleost fish Sparus aurata. Evidence for overall conservation in gene organization and bone-specific expression from fish to man. Gene Expr Patterns. 2001, 270: 77-91.
Desbois C, Hogue DA, Karsenty G: The mouse osteocalcin gene cluster contains three genes with two separate spatial and temporal patterns of expression. J Biol Chem. 1994, 269: 1183-1190.
Renn J, Winkler C: Characterization of collagen type 10a1 and osteocalcin in early and mature osteoblasts during skeleton formation in medaka. J Appl Ichthyol. 2010, 26: 196-201. 10.1111/j.1439-0426.2010.01404.x.
Lie KK, Moren M: Retinoic acid induces two osteocalcin isoforms and inhibits markers of osteoclast activity in Atlantic cod (Gadus morhua) ex vivo cultured craniofacial tissues. Comparative Biochemistry and Physiology, Part A. 2012, 161: 174-184. 10.1016/j.cbpa.2011.10.023.
Confavreux CB, Levine RL, Karsenty G: A paradigm of integrative physiology, the crosstalk between bone and energy metabolisms. Mol Cell Endocrinol. 2009, 310: 21-29. 10.1016/j.mce.2009.04.004.
Ferron M, Hinoi E, Karsenty G, Ducy P: Osteocalcin differentially regulates cell and adipocyte gene expression and affects the development of metabolic diseases in wild-type mice. Proceedings of the National Academy of Sciences USA. 2009, 105 (13): 5266-5270.
Westerfield M: The zebrafish book. A guide for the laboratory use of zebrafish Danio rerio. 2000, Eugene: University of Oregon Press, 4
Connolly MH, Yelick PC: High-throughput methods for visualizing the teleost skeleton: capturing autofluorescence of alizarin red. J Appl Ichthyol. 2010, 26: 274-277. 10.1111/j.1439-0426.2010.01419.x.
Huysseune A, Sire JY: Development of cartilage and bone tissues of the anterior part of the mandible in cichlid fish: a light and TEM study. Anat Rec. 1992, 233: 357-375. 10.1002/ar.1092330304.
Witten PE, Hansen A, Hall BK: Features of mono and multinucleated bone resorbing cells of the zebrafish Danio rerio and their contribution to skeletal development, remodeling and growth. J Morphol. 2001, 250: 197-207. 10.1002/jmor.1065.
Schenk RK, Olah AJ, Herrmann W: Preparation of calcified tissues for light microscopy. Methods of Calcified Tissue Preparation. Edited by: Dickson GR. 1984, Amsterdam: Elsevier, 1-56.
Miyake T, Cameron AM, Hall BK: Stagespecific expression patterns of alkaline phosphatase during development of the first arch skeleton in inbred C57BL/6 mouse embryos. J Anat. 1997, 190: 239-260. 10.1046/j.1469-7580.1997.19020239.x.
Ortego LS, Hawkins WE, Walker WW, Krol RM, Benson WH: Detection of Proliferating Cell Nuclear antigen in tissues of three small fish species. Biotech Histochem. 1994, 69 (6): 317-323. 10.3109/10520299409106312.
Verstraeten B, Sanders E, Huysseune A: Whole Mount Immunohistochemistry and In Situ Hybridization of Larval and Adult Zebrafish Dental Tissues. In: Methods Mol Biol. vol. 2012, 887: 179-191. 10.1007/978-1-61779-860-3_16.
Clément A, Wiweger M, Hardt S, Rusch MA, Selleck SB, Chien C-B, Roehl HH, et al: Regulation of Zebrafish Skeletogenesis by ext2/dackel and papst1/pinscher. PLoS Genet. 2008, 4: e1000136-10.1371/journal.pgen.1000136.
Knopf F, Hammond C, Chekuru A, Kurth T, Hans S, Weber CW, Mahatma G, Fisher S, Brand M, Schulte-Merker S, et al: Bone regenerates via dedifferentiation of osteoblasts in the zebrafish fin. Dev Cell. 2011, 20: 713-724. 10.1016/j.devcel.2011.04.014.
Chomczynski P, Sacchi N: Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987, 162: 156-159.
Tang R, Dodd A, Lai D, Mcnabb WC, Love DR: Validation of Zebrafish (Danio rerio) reference genes for quantitative real-time RT-PCR normalization. Acta Biochim Biophys Sin. 2007, 39 (5): 384-390. 10.1111/j.1745-7270.2007.00283.x.
Nybelin O: Zur morphologie und terminologie des schwanzskelettes der Actinopterygier. Arkiv fur Zoologi. 1963, 15: 485-516.
Patterson C: The caudal skeleton in Lower Liassic pholidophorid fishes. Bulletin of the British Museum (Natural History) Geology. 1968, 16 (5): 201-239.
Arratia G, Schultze H-P: Reevaluation of the caudal skeleton of certain actinopterygian fishes. III. Salmonidae. Homologization of caudal skeletal structures. J Morphol. 1992, 214: 187-249. 10.1002/jmor.1052140209.
A. Bensimon-Brito acknowledges a FCT PhD fellowship SFRH/BD/40573/2007 and is presently the recipient of an ERC Reseal (2007-StG-208631) fellowship, J. Cardeira is the recipient of a CCMAR fellowship CCMAR/BI/0029/2009 and A. Huysseune and P.E. Witten acknowledge a grant from FWO 3G.0040.08. LC acknowledges CCMAR funding. The authors thank Stefan Schulte-Merker for providing the transgenic zebrafish line Tg(oc2:gfp; osx:mcherry). For assistance in sample collection we thank G. Dionísio and R. Peres dos Santos. For technical assistance we thank Barbara Verstraeten, Myriam Claeys, Tommy D’heuvaert and Jasper Dewit.
The funding agencies had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. The authors have declared no competing interests.
Conceived and designed the experiments: ABB, MLC, AH, PEW. Performed the experiments: ABB and JC. Analyzed the data: ABB, JC, MLC, AH, PEW. Contributed reagents/materials/analysis tools: MLC, AH. PEW. Contributed to the writing of the paper: ABB, MLC, AH, PEW. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.