Molecular pathology of vertebral deformities in hyperthermic Atlantic salmon (Salmo salar)

Industrial fish farming makes use of intensive production regimes in an effort to decrease production time and costs. Elevated water temperatures are commonly applied, often without explicit control of factors like nutrition, water quality, densities and vaccination. The intensive rearing systems are unfortunately correlated with deformities affecting both skeletal and soft tissues [1, 2]. In teleosts, hyperthermia can induce vertebral deformities both during the embryonic development and after the vertebral column has been established [3–5]

The teleost vertebral body is built using a minimal bone mass to reduce negative buoyancy [6]. In salmon, the vertebral body comprises four mineralized or ossified layers. Formation of the different layers involves the balanced and highly regulated formation of bone and cartilaginous structures through patterns of mineralization and matrix deposition [7]. The specialized architecture makes it vulnerable to alterations in its tissue composition. Intramembranous ossification occurs by coordinated processes of production, maturation and mineralization of osteoid matrix [8]. Initially osteoblasts produce a thickening osteoid seam by collagen deposition without mineralization. This is followed by an increase in the mineralization rate and the final stage where collagen synthesis decreases and mineralization continues until the osteoid seam is fully mineralized. As part of the process, mineralization time lag appears to be required for allowing modifications of the osteoid so that it is able to support mineralization [9]. Indeed, fast growing Atlantic salmon has been shown to exhibit low vertebral mineral content and mechanical strength, together with an increased risk of developing vertebral deformities [10, 11].

Skeletal growth depends upon the dynamic equilibrium between cartilage production and bone apposition rate [12]. Ontogeny and growth of the vertebral column is under control of regulatory mechanisms involving transcription factors, signaling molecules and extracellular matrix proteins. The pathways of chondrocyte and osteoblast differentiation are interconnected during vertebral formation and must be coordinated. In particular, regulatory proteins, like the transcription factors Sox9, Runx2, Osterix, Twist and Mef2c have distinct functions both in the establishment of the vertebral bodies and later in the differentiation and maturation of specific skeletal cell types (review [13]). Similarly, signaling molecules like bone morphogenetic proteins (Bmp2 and Bmp4), and hedgehog proteins (Ihh and Shh) plays different roles both during cell differentiation and skeletal tissue ontogeny [14–16]. Osteoblasts and chondrocytes secrete the collagen fibers and ground substances of bone and cartilage. These cells are also responsible for the mineralization of the matrix through secretion of specialized molecules, such as Alkaline phosphatase (ALP), Osteocalcin and Osteonectin that binds inorganic minerals [17, 18]. A widely accepted view is that the spatial restriction of ECM mineralization to bone is explained by osteoblast-specific gene products that initiate the formation of hydroxyapatite crystals (Ca₁₀[PO₄]⁶[OH]²) [19]. The requirement for specifically expressed genes in osteoblasts (e.g. col1, osteocalcin and osteonectin) and chondrocytes (e.g. col2 and col10) to initiate the formation of matrix or control the growth of hydroxyapatite crystals is supported by numerous studies [18, 20, 21]. Furthermore, Matrix metalloproteinases (MMPs) and Tartrate-resistant acid phosphatase (TRAP) are involved in degradation of ECM and in the bone remodeling process performed by the osteoclasts [22].

In this work, 20 skeletal genes were used to study the effect of long term hyperthermic exposure on vertebral development and growth in Atlantic salmon. Fish exposed to high temperature (high intensive regime) had a significant higher incidence of deformities than fish from the same origin reared under a conservative temperature regime (low intensive regime). The study was aimed at exposing differences in risk level between the groups, rather than elaborating the pathologies of deformed vertebrae, hence, the study concentrated on phenotypically normal fish from both temperatures. Significant changes in gene transcription were found between phenotypically normal vertebrae of both groups, including down-regulation of genes encoding proteins important for mineralization. Further, in situ hybridization (ISH) and histological staining revealed phenotypical and functional changes in the arch centra. Our results are of basic interest for understanding bone metabolism and deformities, as well as a tool for assessing fish welfare in practical farming.

In the present study we analyzed and compared Atlantic salmon vertebrae from high and low temperature intensity regimes. Rate of development and growth was influenced by temperature regime as observed through SGR and time of sampling. The development from fertilization to first feeding lasted 5 months in the low intensive regime at 6°C, compared to 3 months in the high intensive regime at 10°C. Juveniles of the high intensive group also grew more rapidly after start-feeding than the low intensive group, where the former reached 2 g in 6 weeks after first feeding, 15 g in 3 months and 60 g in 7 months after first feeding, at a rearing temperature of 16°C. In comparison, the low intensive group at rearing temperature of 10°C reached similar sizes in 11 weeks, 5 months and 10 months, respectively. Accordingly, after start-feeding fish from the high intensive temperature regime displayed a higher SGR than the low temperature fish, 2.82 and 1.96 respectively.

Radiography, morphology and mineral analyses

On radiography analysis, the incidence of fish with skeletal abnormalities at 2 g size was 4.0 ± 2.8% and 10.0 ± 1.7% in the low and high intensive groups, respectively (n.s.; not significant). At 15 g size, the difference was more pronounced, 3.4 ± 2.0% and 17.9 ± 1.3% (p < 0.001). At the final sampling at 60 g size, 8 ± 1.4% of the fish in the low intensive group displayed some degree of skeletal pathology compared to 28.1 ± 2.3% in the high intensive group (p < 0.0001), results are shown in figure 1.

Morphometric analyses of vertebral shape demonstrated that fish classified as having a normal phenotype in both groups had more or less regularly shaped vertebrae, but that there was a difference in length-height proportion of vertebrae between fish from the two temperature regimes. Measurements on X-ray images showed that vertebral bodies from the high intensive groups were significantly shorter in craniocaudal direction compared to those from the low intensive groups. The ratios for the high and low intensive group were at 2 g 0.68 ± 0.02 and 0.76 ± 0.02, at 15 g 0.78 ± 0.03 and 0.89 ± 0.06 and at 60 g 0.86 ± 0.01 and 0.94 ± 0.01, respectively (p < 0.001). Examples of vertebral columns with normal phenotype from the high and low intensive group at 15 g are shown in figure 2.

Due to the built-in image contrast enhancement procedures of the semi-digital X-ray system, evaluation of skeletal mineralization as judged by radio density in images was impaired. Nevertheless, a lower contrast in skeletal structures was observed in the high intensity fish, in particular at the 15 g sampling, indicative of a lower mineralization rate at this stage.

Whole body mineral content at the end of the experiment (60 g size) showed low values for Ca, P and Zn content for both temperature regimes (Table 1), with no significant differences between treatments. There was a small, but significant lower level of whole body Fe and Na in the high intensive group. All Fe and Na values were lower than reference values [23], but in correspondence with Ca, P and Zn values, they were within a range which is commonly seen in commercially reared salmon.

	Low intensive group	High intensive group
P	3828 ± 162	3886 ± 285
Ca	3655 ± 341	3700 ± 677
Mg	302.3 ± 8	299 ± 7
Na	1303 ± 70	1084 ± 87*
Fe	9.7 ± 1.2	8.1 ± 0.2*
Mn	0.86 ± 0.05	0.8 ± 0.09
Zn	31.8 ± 2.6	33.5 ± 1.7
Cu	0.8 ± 0.12	0.7 ± 0.09

Mineral content (ppm) in whole body of Atlantic salmon at 60 g from low and high intensive temperature group. Total mineral content was measured in whole fish. Levels were low, but within the range of what are common in commercially reared salmon. Significant numbers (P = 0.05) indicated by *, n = 40 (means ± st.dev).

Quantitative vertebral mRNA expression

The skeletal genes were divided into three groups according to function; ECM constituents, transcription factors, and signaling molecules (Figure 3A-C).

ECM constituents included genes involved in bone matrix production and mineralization and 7 out of 9 of these genes were found to be down-regulated in high intensive group at 2 and 15 g (Figure 3A). Transcription of col1a1, osteocalcin, decorin, osteonectin, mmp9 and mmp13 were reduced in the high intensive group compared to the low intensive group. Col2a1 transcription was also down-regulated at both developmental stages, however the values were insignificant. Osteocalcin was severely down-regulated in 2 g high intensive group. Converse transcription profiles could be observed for col10a1 and alp between 2 g and 15 g fish; col10a1 was down-regulated at 2 g and up-regulated at 15 g whereas alp was up-regulated at 2 g and down-regulated at 15 g.

Temporal changes in transcription factor mRNA expression were found between high and low temperature group, and all genes except sox9 showed opposite expression at 2 and 15 g (Figure 3B). In the high intensive group, sox9 was down-regulated at 2 g (n.s.) and 15 g, but more pronounced in the latter. Investigation of the two osteoblast markers runx2 and osterix, revealed opposite mRNA expression levels at 2 and 15 g. Runx2 was up-regulated at 2 g (n.s), but down-regulated at 15 g. On the contrary, osterix was down-regulated (n.s.) at 2 g, but up-regulated at 15 g. Mef2c and twist was also down-regulated at 2 g, while up-regulated at 15 g (twist n.s.).

Signaling molecules included bmp2, bmp4, shh and ihh. Expression analysis of mRNA for signaling molecules showed statistically significant differences in expression levels between the temperature regimes and all transcripts were found more abundant in the 15 g group when compared to 2 g vertebrae. Bmp2 was the only up-regulated signaling molecule at 2 g, while all signaling genes were up-regulated at 15 g (Figure 3C).

To further examine changes in chondrocyte recruitment and structure between the temperature regimes, we included platelet derived growth factor receptor b (pdgfrb) and vimentin, because of their importance in proliferation and the cytoskeleton, respectively [24, 25]. Both transcripts were significantly down-regulated in 2 g, while significantly up-regulated at 15 g (Figure 3A, B).

In summary, we found that out of the 20 genes we analyzed, 8 were down-regulated in both temperature groups (col1a1, col2a1, osteocalcin, sox9, decorin, osteonectin, mmp9 and mmp13), 9 genes were up-regulated in the 15 g high intensive group, but down-regulated at 2 g (bmp4, col10a1, osterix, ihh, shh, mef2c, twist, vimentin and pdgfrb). And finally, alp and runx2 were up-regulated at 2 g but down-regulated at 15 g.

Vertebral tissue morphology and spatial mRNA expression

In areas where osteoblasts secrete the osteoid matrix, a generally stronger ISH signals was apparent in the low intensive group for all probes. The osteogenic marker gene col1a showed distinct staining to osteoblasts at the growth zone of the endbones of the vertebral bodies from fish of both temperature regimes (Figure 4A-C). Moreover, col1a signal was identified in the bone lining osteoblast cells situated at the lateral surfaces of the trabeculae and along the rims of the vertebral bodies. Investigation of osteocalcin mRNA revealed an expression pattern similar to col1a, with staining of cells in the osteogenous areas and in bone lining osteoblasts and apical surfaces of the trabeculae (Figure 4D-G). Specifically high osteocalcin signal was detected in the proliferative osteoblast growth zones on the endbones of the vertebral bodies. Osteonectin mRNA was detected in the osteogenic growth zone of the endbones and lining the exterior part of the vertebral body (Figure 5A-D). The chondrocytic marker col2a, hybridized heavily to chordoblasts in the notochord (Figure 5E-G), whereas col10a was detected in a continuous layer of cells along the rims of the vertebral body (Figure 5J-L).

Alizarin red S and toluidine blue stained chondrocytes in the arch centra and revealed distinct morphological differences between vertebrae from the two temperature groups. The low intensive group was defined by distinct sub-groups of chondrocytes in the different maturational stages i.e. resting, proliferating and hypertrophic. In contrast, the equivalent chondrocytes were more distorted in the high intensive group (Figure 6A, B). ISH analysis of col2a, col10a and osteonectin enabled classification of the different chondrocytes into distinct sub-populations of maturational development. Col2a hybridized to resting and pre-hypertrophic chondrocytes in two distinct bands of both low and high intensive group, but the mRNA expression was more evenly distributed in all cells of the latter group (Figure 5H, I). There were also generally less proliferating chondrocytes that tended to be less compact in this group. In proliferating chondrocytes we detected strong col2a mRNA expression in the high intensive group, but no expression in the low intensive group. Analysis of col10a showed restriction to the pre-hypertrophic and hypertrophic chondrocytes located in the deep cartilage zone (Figure 5M, N). Osteonectin was also expressed in chondrocytes and the signal increased towards the hypertrophic chondrocytes (Figure 5A-C). The pre-hypertrophic chondrocyte zone was found to be expanded in the high intensive fish and both col10a1 and osteonectin showed an expanded expression domain corresponding to an increased hypertrophic zone. No signal was detected in any of the samples hybridized with sense probes (data not shown).

In normal spinal columns from the low intensive group, positive TRAP staining was detected at the ossifying boarders of the hypertrophic chondrocytes in the arch centra. No positive staining was detected in samples from the high intensive group (Figure 7A, B).

The presented study aims at describing the molecular pathology underlying the development of vertebral deformities in Atlantic salmon reared at a high temperature regime that promotes fast growth during the early life stages. Within the period investigated, vertebral bodies form and develop and the skeletal tissue mineralizes. Rearing at high temperatures resulted in higher frequencies of vertebral deformities, as expected. The vertebral pathology observed in this study was most likely induced both during the embryonic development and after start-feeding, since the incidence of deformities continued to increase throughout the experiment after the first radiographic examination at 2 g. Similar temperature regimes before and after start-feeding have independently been shown to induce vertebral defects in juvenile salmon [1, 26]. However, whereas high temperatures during embryonic development is commonly related to somitic segmentation failure, deformities later in development may possibly be linked to fast growth induced by elevated temperatures and the impact this might have on the natural maturation and ontogeny of the vertebral bodies [3, 10, 11]. This causative relation has been shown for fast growing underyearling smolt that has a higher incidence of vertebral deformities than slower growing yearling smolt [27]. Further, morphometric analyses showed that elevated water temperature and faster growth is manifested by a difference in length-height proportion of vertebrae between fish from the two temperature regimes. Similar decrease in length-height proportion was described for the rapid-growing underyearling smolt [27]. Radiographic observations indicated a lower level of mineralization of osteoid tissues in the high temperature fish. However, we could not find any pronounced altered mineral content between the two temperature regimes. The observed values were low compared to reference values [23], but in a range commonly observed in commercially reared salmon. Apparently, whole body mineral analysis seems insufficient to assess problems related to the development of spinal deformities.

To determine whether the difference in likelihood of developing vertebral deformities between the two groups could be traced back to an altered gene transcription, we examined the expression of selected skeletal mRNAs in phenotypical normal salmon fry at 2 and 15 g. Histological examination of 15 g fish was included to improve interpretation of the transcriptional data. The selected genes showed conservation and similar spatial expression with those examined in other vertebrates, supporting that most of the factors and pathways that control skeletal formation are highly conserved in vertebrates.

The lower transcription of ECM genes such as col1a1, osteocalcin, osteonectin and decorin suggests a defect in the late maturation of osteoblasts [17, 21, 28]. The correlation to impaired mineralization is supported by the shorter vertebral bodies in the high intensive groups throughout the study, as well as the impaired mineralization indicated by low contrast observed on X-ray. Col1a1 is the primary ECM component secreted by osteoblasts in the trabecular bone and growth plate and defects in the synthesis of col1 or type 1 procollagen have been found in several heritable disorders of connective tissue. Likewise, defects in the assembly of Col1 fibrils have been reported to cause abnormally thin and branched structures [29]. Decreased diameter and crosslink density of the collagen fibers have been suggested to reduce thermal stability of collagen and thereby the tissues ability to support load during elevated temperatures [30]. In chum salmon, Oncorhynchus keta, the denaturation temperature of collagen type 1 from skin has been reported to be about 19°C [31]. The collagen fibres are further organized and stabilized by a range of non-collagenous proteins, which functions by linking other proteins and minerals to the ECM scaffold. Decorin, which belongs to the small leucine rich repeat proteoglycan (PG) group (SLRPs) is involved in determining the mature collagen fibril structural phenotype and tissue function by facilitating protein-protein interaction with a range of other matrix components (mainly collagen fibres) and with the mineral phase during the formation of calcified tissues [32]. As a result, decorin has been shown to increase tensile strength of the collagen-decorin fiber [33]. Further, osteonectin is a phosphorylated glycoprotein that binds to collagen fibrils, calcium, and hydroxyapatite, linking the bone mineral and collagen phases and perhaps initiating active mineralization in normal skeletal tissue [17, 34]. Osteonectin-null mice display decreased trabecular bone volume and have bone of lesser stiffness than control mice [35]. Osteocalcin mRNA expression also serves as a useful molecular marker of mineralization because it is associated with the maturation of bone cells and mineralization [18, 36]. Alp is another marker gene for bone cell maturation and mineralization. Inhibition of alp activation, by for example heat or by gene knockout, inhibits calcification and causes mineralization defects in cultured bone cells and mice [37]. In addition, mutations in the alp gene lead to hypophosphatasia, in which bone matrix formation occurs, but mineralization is inhibited [38]. Our results showed that alp was down-regulated in the high intensive 15 g group, but up-regulated in 2 g fish. This may indicate that alp is a limiting factor for mineralization after long term exposure to the high temperature regime. Altogether, the simultaneous down-regulation of genes encoding structural proteins taking part in the bone matrix and mineralization strongly supports an assumption that disturbances of these processes constitute an important part of the mechanisms of development of vertebral deformities.

As for the ECM genes involved in osteoblast development and mineralization, high intensive temperature treatment had a significant effect on the transcription of transcription factors and signaling molecules involved in these processes. Intriguingly, Runx2 and Osterix, known as master regulators of osteoblast differentiation [39, 40], exhibited opposite mRNA expression levels at 2 and 15 g. Runx2-null mice have osteoblast differentiation arrested [41], while osterix-null mice embryos have a significant reduction of col1 expression and do not express the late osteoblast specific marker osteocalcin [39]. In addition, we analyzed the bHLH transcription factor twist. This gene works as a negative regulator of osteoblastogenesis by inhibiting expression of genes downstream of runx2[42]. At 2 g when osterix and twist was down-regulated while runx2 was up-regulated, osteocalcin was heavily down-regulated as was col1a1. The mRNA expression pattern was inverted at 15 g. Then osterix and twist was up-regulated and runx2 down-regulated, while osteocalcin and col1a1 were weakly down-regulated. Linking these results to the pathways involved in osteoblast development, the required simultaneous activation of osterix and runx2 did not appear at 2 g or at 15 g. However, Osterix function downstream of Runx2 during osteoblast differentiation, but may be regulated by Bmp2 in a Runx2-independent pathway [43]. Bmp2 can induce ectopic bone and cartilage formation in adult vertebrates [44]. Spinella-Jaegle et al [16] found that cooperation between Bmp2 and Shh was necessary to promote a strong induction of the osteoblast marker alp in human mesenchymal cell lines. At both 2 and 15 g, bmp2 was highly up-regulated in the high intensive group, possibly as a response to the low ECM mRNA expression and under-mineralized tissue. In addition, osterix and shh was up-regulated at 15 g, as was bmp4. Bmp4 treatment has been shown to stimulate new bone formation and is also expressed in osteoblasts prior to formation of mineralized bone nodules [45, 46]. However, in comparison to Spinella-Jaegles in vitro findings, we did not detect an increase in alp mRNA expression. Further, we detected a weaker signal of osteocalcin and osteonectin in osteoblasts from the ISH of the high intensive group at 15 g. Hence, despite the possible attempt of bmp2 to restore bone formation and mineralization, there was still lower transcription of ECM components in the high intensive group at 15 g. Summarized, our results may indicate that osteoblast proliferation and mineralization were restrained in the fast growing group.

The percentage of deformities significantly increased in the high intensive group from 2 g till 15 g, while the percentage was stable in the low intensive group. Hence, this period seems to involve important steps for the developmental fate of deformities. Between these two size stages we observed a change in expression pattern, from a downregulated to an upregulated transcription, of 9 genes, where 8 of them are involved in chondrogenesis. This suggested that chondrocytes go through changes in this period that could be important for the development of the observed pathologies.

In vertebrates as mouse and human [47, 48], the growth zones of long bones consists of well defined layers of progenitor, proliferative and hypertrophic chondrocytes [49]. These chondrocytes differ in their morphology, proliferation abilities and secretion of ECM components. For example, transcription of col2a1 is characteristic for the proliferative state whereas col10a1 is restricted to the hypertrophic state [47, 50]. ISH of these genes revealed that 15 g Atlantic salmon raised at the low intensive regime also had distinct sub-populations of progenitor, proliferative and hypertrophic chondrocytes at the growth zone of the neural and haemal arches. On the contrary, more distorted layers were found in Atlantic salmon raised at the high intensive regime. Moreover, an increased zone of hypertrophic chondrocytes was found in the proximity of the mineralized bone matrix in the high intensive group. Once these hypertrophic chondrocytes are fully differentiated, matrix calcification would normally be initiated [12]. However, we could not identify any variance in mineralization at the ossifying borders of the hypertrophic chondrocytes when examined by histological Alizarin red S staining.

The increased zone of hypertrophic chondrocytes in the high intensive group and the up-regulated transcription of hypertrophic marker genes suggest an arrest prior to the final maturation of chondrocytes. Thus, these chondrocytes seems unable to initiate mineralization. The chondrocyte hypertrophy marker col10a1 and its activator mef2c[51] were both up-regulated at 15 g in the high intensive group. Moreover, ihh, a repressor of terminal hypertrophic differentiation [52], was found to be highly up-regulated, whereas sox9, which is involved in early chondrocyte differentiation, and its downstream structural protein col2a[53], were down-regulated. The severely down-regulation of runx2 at 15 g is of interest, since runx2-null mice embryos have a narrow zone of proliferating chondrocytes and a wide zone of hypertrophic chondrocytes [54]. In addition, bmp4, which was up-regulated at 15 g, has been shown to accelerate the hypertrophic maturation process [55]. Interestingly, we also found an up-regulated expression of pdgfrb mRNA at 15 g. Kieswetter and collaborators [25] have reported that chondrocytes respond to PDGF by enhancing proliferation and cartilage matrix production while maintaining the cells in a less mature phenotype; corroborating our findings that the chondrocytes are some how arrested in the late hypertrophic stage at 15 g with a reduced possibility of completing the endochondral ossification process with calcified bone as end product. Similar findings have also been shown in rat ulnae, where loading was associated with an increased hypertrophic zone in the growth plate [56], but mineralization rate was suppressed [57]. Another interesting comparative pathological condition to our findings in salmon is tibial dyschondroplasia (TD), a metabolic disease of young poultry that affects the growth of bone and cartilage. The lesion is morphologically characterized by an accumulation of chondrocytes that appear to be unable to differentiate past a pre-hypertrophic stage [58]. TD often occurs in broilers and other poultry that have been bred for fast growth rates. The tibial cartilage does not mature enough to ossify, which leaves the growth plate prone to fracture, infection, and deformed bone development.

The observed shorter phenotype of vertebral bodies from the high intensive group might have been a consequence of higher mechanical load in fast growing fish coincidental with a lower transcription of supportive ECM components. Together with the up-regulation of hypertrophic genes in high intensive fish at 15 g, we also found increased transcription of vimentin. Vimentin filaments have been shown to regulate the swelling pressure of chondrocytes [59] and strengthen resistance to mechanical stress [60]. Hence, the increased activation of vimentin and the increased proportion of hypertrophic chondrocytes in the high intensive temperature group at 15 g may reflect an adaptation to the fast growth by prioritizing maturation of chondrocytes that are more resistant to mechanical stress. At 2 g, however, the reduced level of vimentin mRNAs might possibly be linked to the mal-adaptive down-regulation of chondrocytic genes in high intensive group. Indeed, disruption of vimentin filaments has been shown to result in loss of cell contact with the surrounding matrix which may alter the signaling dynamics of the cell and in effect shut down transcriptional events [24].

Mineralizing hypertrophic chondrocytes acquire and express most of the phenotypic characteristics of osteoblasts, including high Alp activity and expression of osteonectin and osteocalcin[61]. These phenotypic traits shared with osteoblasts may be needed to bring about the final phase of endochondral ossification and replace mineralized cartilage with bone [62]. They may also permit mineralized cartilage to act as bone-like structural tissue and allow for a transition from cartilage to bone. In contrast to the down-regulated transcription of osteonectin and osteocalcin, as determined by real time qPCR, we observed an increased transcription pattern of these genes in the arch centra in the high intensive group by ISH. We also observed a tendency of lower transcription of the same genes in osteoblasts of the high intensive group. However, establishment of a calcifiable matrix requires degradation of some matrix molecules. Endochondral bone formation includes the participation of MMPs, which degrade cartilage matrix and allow vascular invasion [63]. At least two proteases are involved in this process; MMP13 which regulates remodeling of the hypertrophic cartilage matrix and MMP9 which has a role in vascularisation of the growth plate [64, 65]. When analyzing these MMPs in salmon vertebral columns, a significant down-regulation of both mmp9 and mmp13 in the high intensive group at 2 g were observed. At 15 g, mmp13 mRNA expression decreased even more, while mmp9 was significantly up-regulated. Indeed, MMP13 is known as the dominant collagenase in cartilage and its absence cause delay in endochondral ossification [65]. Further supporting the hypothesis that endochondral ossification was in some way delayed in the spinal columns from the high intensive group, runx2 deficiency has been shown to inhibit mmp expression [66] and lead to mild disturbances of chondrocyte differentiation, as discussed above. In addition, TRAP activity, essential for completing endochondral ossification [22], was absent in the erosive front of cartilage in neural and heamal arches of spinal columns from the high temperature group.

The presented results contribute to the understanding of the mechanisms involved in development of temperature-induced vertebral pathology by describing changes in vertebral tissue not yet manifesting pathological deviations. Our results strongly indicate that temperature induced fast growth is severely affecting gene transcription in osteoblasts and chondrocytes, leading to a change in the tissue structure and composition. The data presented here indicate that both production of bone and cartilage were disrupted when promoting fast growth using elevated temperature. It is not unlikely that this disequilibrium is involved in the higher rate of deformities observed in the high intensive group. Importantly, management control of deformities and health in general demands precise tools and knowledge to depict any problem as early as possible in the production line. The defined markers of bone and cartilage cell differentiation and matrix formation can be used to investigate how the progression of skeletogenesis is modulated by a variety of factors. Even though differences in the two experimental groups were undetectable externally, rearing at increased temperatures induced consistent transcriptional changes in several genes that correlated with the higher risk of developing deformities later in ontogeny. Hence, this article reveals the potential use of gene transcription profiling as a prognostic approach in aquaculture.